Profiling the tyrosine phosphoproteome of different mouse mammary tumour models reveals distinct, model-specific signalling networks and conserved oncogenic pathways
Breast Cancer Research volume 16, Article number: 437 (2014)
Although aberrant tyrosine kinase signalling characterises particular breast cancer subtypes, a global analysis of tyrosine phosphorylation in mouse models of breast cancer has not been undertaken to date. This may identify conserved oncogenic pathways and potential therapeutic targets.
We applied an immunoaffinity/mass spectrometry workflow to three mouse models: murine stem cell virus-Neu, expressing truncated Neu, the rat orthologue of human epidermal growth factor receptor 2, Her2 (HER2); mouse mammary tumour virus-polyoma virus middle T antigen (PyMT); and the p53−/− transplant model (p53). Pathways and protein–protein interaction networks were identified by bioinformatics analysis. Molecular mechanisms underpinning differences in tyrosine phosphorylation were characterised by Western blot analysis and array comparative genomic hybridisation. The functional role of mesenchymal–epithelial transition factor (Met) in a subset of p53-null tumours was interrogated using a selective tyrosine kinase inhibitor (TKI), small interfering RNA (siRNA)–mediated knockdown and cell proliferation assays.
The three models could be distinguished on the basis of tyrosine phosphorylation signatures and signalling networks. HER2 tumours exhibited a protein–protein interaction network centred on avian erythroblastic leukaemia viral oncogene homologue 2 (Erbb2), epidermal growth factor receptor and platelet-derived growth factor receptor α, and they displayed enhanced tyrosine phosphorylation of ERBB receptor feedback inhibitor 1. In contrast, the PyMT network displayed significant enrichment for components of the phosphatidylinositol 3-kinase signalling pathway, whereas p53 tumours exhibited increased tyrosine phosphorylation of Met and components or regulators of the cytoskeleton and shared signalling network characteristics with basal and claudin-low breast cancer cells. A subset of p53 tumours displayed markedly elevated cellular tyrosine phosphorylation and Met expression, as well as Met gene amplification. Treatment of cultured p53-null cells exhibiting Met amplification with a selective Met TKI abrogated aberrant tyrosine phosphorylation and blocked cell proliferation. The effects on proliferation were recapitulated when Met was knocked down using siRNA. Additional subtypes of p53 tumours exhibited increased tyrosine phosphorylation of other oncogenes, including Peak1/SgK269 and Prex2.
This study provides network-level insights into signalling in the breast cancer models utilised and demonstrates that comparative phosphoproteomics can identify conserved oncogenic signalling pathways. The Met-amplified, p53-null tumours provide a new preclinical model for a subset of triple-negative breast cancers.
It is well established that tyrosine kinase signalling pathways play major roles in breast cancer development and progression. Evidence for this includes (1) the amplification of the receptor tyrosine kinase (RTK) human epidermal growth factor receptor 2 (HER2)/avian erythroblastic leukaemia viral oncogene homologue 2 (erbB2) in approximately 25% of breast cancers ; (2) the presence of activating mutations in PIK3CA, encoding the p110α catalytic subunit of phosphatidylinositol 3-kinase (PI3K), in approximately 30% of cases ; and dysregulation of other signalling intermediates such as the non-RTK Src , the docking protein Gab2  and the serine/threonine kinase Akt3 . In addition, an array of data from cell line and mouse models of breast cancer confirm the oncogenic or tumour suppressor roles of particular tyrosine-phosphorylated signalling proteins -. These findings have led to the development of small-molecule– and antibody-based targeted therapies that have either entered the clinic or are currently undergoing clinical trials, including the use of trastuzumab to treat HER2-positive breast cancer .
Gene expression profiling has revealed that breast cancer can be subclassified into luminal A, luminal B, HER2, basal and claudin-low subtypes, which differ in terms of patient prognosis and response to therapy . Basal and claudin-low breast cancers present a major clinical challenge, as their frequent triple-negative phenotype (lacking HER2, oestrogen and progesterone receptors) confers intrinsic resistance to HER2-targeted and endocrine therapy . The HER2 subtype is characterised by amplification of the corresponding gene, and other breast cancer subtypes also exhibit characteristic perturbations in tyrosine kinase signalling. For example, increased expression of the RTKs Met and epidermal growth factor receptor (Egfr) are associated with basal and triple-negative breast cancers -, and, at least in cell line models, the basal and claudin-low subtypes are characterised by a prominent Src family kinase (SFK)–governed network ,. Further characterisation of the signalling networks associated with different breast cancer subtypes may identify novel therapeutic targets and new applications for existing therapies, as well as biomarkers that aid patient stratification for therapy.
Genetically engineered mouse (GEM) models have greatly contributed to our understanding of the mechanisms and functions of particular oncogenes and tumour suppressors implicated in breast cancer, as well as biological processes involved in tumour progression . These include GEM models in which normal or activated forms of neu (the rat homologue of HER2/erbB2) are expressed in the mammary gland , as well as the polyoma virus middle T antigen (PyMT) model . PyMT mimics an activated RTK, localising to the plasma membrane, as well as in intracellular membranes, where it is phosphorylated on specific tyrosine residues by particular SFKs. Phosphorylation of PyMT creates binding sites for Shc (at Y250), the p85 subunit of PI3K (at Y315) and phospholipase Cγ1 (PLCγ1, at Y322), and key roles for the Y250 and Y315/Y322 sites in mammary tumourigenesis have been demonstrated . Interestingly, transcript profiling has demonstrated that mouse mammary tumour virus (MMTV)-Neu and MMTV-PyMT model tumours are relatively homogeneous and exhibit gene expression similarities to human luminal-type cancers . In contrast, the p53-null transplant model of mammary tumourigenesis is characterised by tumours exhibiting histological and molecular heterogeneity, as well as genetic instability ,. Indeed, transcript profiling has demonstrated that p53-null tumours can be classified into basal-like, luminal and claudin-low subtypes characterised by distinct genomic copy number changes .
An important concept emanating from research using GEM models is that comparative genomic and transcriptomic strategies, wherein particular mouse tumours are compared to human breast cancer subtypes, can be utilised identify conserved mechanisms essential for disease development and progression ,,,. In this study, we undertook global tyrosine phosphorylation profiling of three GEM models of breast cancer: a HER2 model featuring expression of an activated form of the receptor lacking the extracellular domain , as well as the MMTV-PyMT and p53-null transplant models. This enabled characterisation of the tyrosine phosphorylation-based signalling networks characteristic of each tumour type, revealed similarities and differences between these tumour models, and identified oncogenic pathways conserved in the human disease.
Neu-pMIL was used to generate the HER2 breast cancer model. This expresses a truncated form of Neu (approximately 647 to 1,260 amino acids, the rodent orthologue of Her2) with elevated activity due to truncation of the extracellular domain . Truncated Neu was subcloned from Neu-pLJ by digestion using Sal1 followed by subcloning into pENTR2B (Invitrogen, Mulgrave, Victoria, Australia). The gene was subsequently cloned into pMIL, a derivative of pMig  that expresses a luciferase marker, to generate Neu-pMIL, in which expression of Neu is driven by the murine stem cell virus (MSCV) promoter. The Gateway LR recombination reaction was used as per the manufacturer's instructions (Invitrogen).
Generation of tumours
MSCV-Neu (HER2) tumours
Primary mammary epithelial cells from FVB/N mice were cultured and retrovirally transduced with truncated Neu encoded by Neu-pMIL as described previously ,,. Cells were transplanted into the cleared mammary fat pad of naïve FVB/N recipients within 2 days of retroviral transduction.
Tp53-null (p53) mice from The Jackson Laboratory (stock no. 002899; Bar Harbor, ME, USA) were on an FVB/N background. Mammary epithelium from 8-week-old mice was transplanted to the cleared mammary fat pad of naïve BL/6 (clone 5101) or FVB/N (all other clones) recipients as described previously ,.
Mice were aged for up to 12 months to permit spontaneous tumour formation . Tumour fragments were passaged into the cleared mammary fat pad of naïve mice as previously described  (Additional file 1: Figure S1).
PyMT tumours were derived from MMTV-PyMT-transgenic mice on an FVB/N background as previously described . All animal work was approved by the animal ethics committee of Garvan Institute of Medical Research, St Vincent's Hospital.
Phosphotyrosine peptide enrichment
Resected mouse mammary tumours were lysed in 8 M urea, 20 mM 2-[4-(2-hydroxyethyl)piperazin-1-yl]ethanesulfonic acid (HEPES), 2.5 mM sodium pyrophosphate, 1 mM β-glycerol phosphate, 1 mM sodium orthovanadate, 1 mM ethylenediaminetetraacetic acid and 1 mM Tris(2-carboxyethyl)phosphine, pH 8.0. Lysates were cleared by sonication and centrifugation prior to protease digestion, phosphotyrosine (pY) immunoprecipitation (IP) and nano-liquid chromatography tandem mass spectrometry (nano-LC-MS/MS). Lysates were quantitated, and 20 mg of each sample were diluted to a final concentration of 1 M urea with 20 mM HEPES and digested overnight with trypsin at room temperature. Heavy proline (+6 Da)– and alanine (+4 Da)–labelled synthetic standard pY peptides of MK14, elongation factor Tu and EGFR were spiked into each sample at 5 pmol, 500 fmol and 10 pmol, respectively, to enable normalisation of label-free quantitative values. Peptides were then desalted and concentrated using C18 Sep-Pak columns (Waters, Milford, MA, USA), and lyophilised. pY peptides were subject to IP using the PhosphoScan procedure (Cell Signaling Technology, Beverly, MA, USA) with pY100 antibodies as previously described . pY-enriched peptides were dried using a vacuum centrifuge and stored at −80°C prior to nano-LC-MS/MS analysis. All samples were analysed by nano-LC-MS/MS using two technical replicates.
Nano–liquid chromatography tandem mass spectrometry
pY peptides were resuspended in 15 μl of MS buffer containing 1% formic acid, 2% acetonitrile (ACN) and 0.05% heptafluorobutyric acid and subjected to nano-LC-MS/MS. Peptides were separated by nano-LC using an UltiMate 3000 high-performance liquid chromatography and autosampler system (Dionex, Sunnyvale, CA, USA). Samples were concentrated and desalted onto a C18 micro-precolumn (500 μm × 2 mm; Michrom Bioresources, Auburn, CA, USA) with H2O:ACN (98:2, 0.05% trifluoroacetic acid) at 15 μl/min for 4 minutes. The micro-precolumn was then switched online with a nano-C18 column (75 μm × about 10 cm, 5 μm, 200 Å Magic; Michrom) and the reverse phase nano-eluent was subjected to positive nanoflow electrospray analysis in information-dependent acquisition (IDA) mode (250 nl/min for 30 minutes). Mass spectra were acquired on an Orbitrap Velos mass spectrometer (Thermo Electron, Waltham, MA, USA). In IDA mode, survey scans in the mass-to-charge ratio (m/z) range 350 to 1,750 were acquired with lock mass enabled. Up to the 15 most abundant ions (>5,000 counts) with charge states greater than +2 were sequentially isolated and further subjected to MS/MS fragmentation within the linear ion trap using collisionally induced dissociation. MS/MS spectra were accumulated with an activation time of 30 milliseconds at a target value of 30,000 ions and a resolution of 60,000. m/z ratios selected for MS/MS were dynamically excluded for 30 seconds.
The nano-LC-MS/MS .raw files were processed with MaxQuant software (version 220.127.116.11), which uses the Andromeda algorithm for data processing, database searching and protein identification . Extracted peak lists were searched against the UniProtKB/Swiss-Prot Mus musculus database (Version 2010_10) containing 35,052 entries (including common contaminants) and a proportionally sized decoy database for false discovery rate (FDR) generation. The following search parameters were selected: fixed cysteine carbamidomethylation modification, variable methionine oxidation modification, variable protein N-acetylation, and variable phosphorylation of serine, threonine and tyrosine. A minimum peptide length of six amino acids and up to two missed cleavages were allowed. The initial first search mass tolerance was 20 ppm for precursor ions and 0.5 Da for fragment ions. Matches between runs were enabled with default settings and label-free quantitation enabled. The FDR was limited to 1% for both protein and peptide identification. For identification of the standard peptides, additional variable modifications of heavy proline (+6 Da) and alanine (+4 Da) were enabled.
Phosphotyrosine peptide quantitation
Raw pY peptide spectral intensities were extracted from the 'Evidence' output files generated in MaxQuant. These intensities were then normalised against pY peptide intensities for the heavy-labelled spiked-in peptide standards. These normalised intensities were then log10-transformed prior to all bioinformatics analyses, and their direct comparison enabled the relative quantitation of pY sites between the p53, PyMT and HER2 GEM models of breast cancer.
Raw MaxQuant output files were subjected to the following filtering criteria prior to their inclusion in further bioinformatics analysis: (1) contaminants and reversed matches were removed, ensuring that all identifications were reported with an FDR of <1%; and (2) pY sites were filtered for a minimum localisation probability confidence of 0.75.
The k-nearest neighbour  method was applied to impute missing values if a pY site was present in more than 50% of all samples; otherwise, the minimal value across all samples or the overall mean per sample for that pY site was used to replace missing values.
Clustering and heat map generation
Heat maps were generated using the R package gplots with hierarchical clustering applied to generate dendrograms for pY sites and samples. Euclidean distance with complete linkage was used, and the values were standardised at rows. The Bioconductor package limma was used to identify differences between the cancer models based on their pY-site expression levels. Empirical Bayes moderated P-values for contrasting and comparing expression between all samples were corrected for multiple testing by the Strimmer method as implemented in the R package fdrtool . For all tests, Q-values less than 0.05 were considered significant. pY sites were thought of as class-specific (Her2-specific, p53-specific, PyMT-specific) on the basis of their contrast significance and the exhibition of increased spectral intensity for a particular sample.
Random forest classifier
The random forest classifier  was built using a two-step approach with the R package randomForest 2.15.0. In the first step, a classifier was built using all 381 consistently identified pY sites (that is, ≥75% of samples) as variables. In the second step, 17 pY sites with the importance score > 1.2 for each tumour model, which were calculated from the first classifier using a permutation-based approach, were used to build a new classifier. The out-of-bag estimate of error rates is 4.76% for both classifiers.
The protein–protein interaction data were extracted from Protein Interaction Network Analysis (PINA) ,; substrate–kinase relationship information was obtained from PhosphoSitePlus ; and mouse-human orthologues extracted from MSOAR . The networks were generated using a Cytoscape PINA4MS plugin (HC Lee, M Pinese, L Bourbon, I Rooman, RJ Daly, AV Biankin, J Wu, manuscript under preparation. URL: http://apps.cytoscape.org/apps/pina4ms).
The KOBAS web server , was applied for pathway enrichment analysis. The hypergeometric test was used to calculate the statistical significance of each pathway, followed by multiple test correction using the Benjamini–Hochberg method .
Antibodies and immunoblotting
Antibodies against the following proteins were used, with immunoblotting carried out as previously described : IRS1 (Upstate Biotechnology, Lake Placid, NY, USA); IRS1 pY612 (Invitrogen), Stat3, Crkl, Jak-1, P85α, CDK16, Cav1, ErbB3, phospho-ErbB3 (Tyr1328; Tyr1325 in mouse), ErbB2, Met, Met pY1234/1235 and pY (P-Tyr-100) (Cell Signaling Technology); phospho-p85 (Tyr467) (GeneTex, Irvine, CA, USA); platelet-derived growth factor receptor α (PDGFRα) and β-actin (Santa Cruz Biotechnology, Santa Cruz, CA, USA); ERBB receptor feedback inhibitor 1 (Errfi1) (Sigma-Aldrich, St Louis, MO, USA).
Mouse tumour samples and liver tissue samples were harvested, snap-frozen and processed using the QIAamp DNA Mini Kit (QIAGEN, Valencia, CA, USA) according to the manufacturer's instructions. Briefly, up to 25 mg of frozen tissue was ground and digested with proteinase K overnight at 56°C. DNA was purified on a QIAamp spin column, eluted in 200 μl of water and assessed for quality using a NanoDrop spectrophotometer (Thermo Scientific, Waltham, MA, USA). Samples were cleaned up by sodium acetate/ethanol precipitation to ensure 260/280 and 260/230 absorbance ratios >1.8. DNA integrity was checked using gel electrophoresis.
Array comparative genomic hybridisation
Array comparative genomic hybridisation (aCGH) was performed by The Ramaciotti Centre for Gene Function Analysis (Randwick, NSW, Australia) on differentially labelled tumour DNA (Cy5) and sex-matched reference liver DNA (Cy3) using the Agilent SurePrint G3 mouse CGH 4 × 180 K microarray platform (Agilent Technologies, Santa Clara, CA, USA). The aCGH data were analysed using circular binary segmentation  to translate intensity measurements into regions of equal copy numbers. The raw array data were median-normalised; duplicate probes were averaged; and the resulting data were smoothed to remove outliers prior to segmentation. The gain and loss status for each region was assigned by determining plateaus in an ordered plot of all segments by their mean log2 ratios. By this method, we denoted deletions to be segments with a mean value less than −0.5, amplifications to have mean value between 0.5 and 1.5 and highly amplified regions to have a mean value greater than 1.5. Using this information, we merged segmental values across the genome to create a common set of copy number levels for each individual tumour. Analysis was carried out using the Bioconductor DNAcopy package  in R Project R2.15.2 software .
Met inhibitor experiments
Tumours were harvested and cellular preparations made using mechanical disruption and collagenase digestion as described previously . Cells were cultured in vitro for two passages in mouse mammary epithelial cell media consisting of Dulbecco's modified Eagle's medium/Ham's F12 nutrient mixture supplemented with foetal calf serum, mouse epidermal growth factor, hydrocortisone, HEPES, insulin, gentamicin and penicillin/streptomycin , then trypsinised and seeded into 96-well plates at a density of 6,000 cells/well. The following day, cells were treated with the Met inhibitor PHA-665752 (Selleck Chemicals, Houston, TX, USA) at concentrations of 0.25 μM, 0.5 μM and 1 μM or vehicle treated with 0.1% dimethyl sulphoxide (DMSO; Sigma-Aldrich). Cell proliferation was assessed and recorded using the IncuCyte ZOOM system (Essen BioScience, Ann Arbor, MI, USA). Briefly, cells were cultured for 5 days while images of cells inside each well (three images per well across four replicate wells) were captured every 2 hours. Proliferation was measured at each time point as a percentage of confluence, as determined using the IncuCyte software on the basis of automated analysis of captured images. Data are representative of two replicate experiments. For immunoblotting studies, cells at 70% to 90% confluence were treated with either vehicle control (DMSO) or Met TKI (1 μM final concentration) for 1 hour and then lysed in 1% (vol/vol) Triton X-100 lysis buffer as previously described .
Met knockdown and cell proliferation assay
Met knockdown was performed using either individual (ON-TARGETplus Met siRNA 1 and 2) or SMARTpool siRNA targeting mouse Met (Dharmacon, Lafayette, CO, USA). The ON-TARGETplus nontargeting pool was used as a negative control. Briefly, cells were seeded into 96-well plates at a density of 4,500 cells/well 16 hours prior to transfection with 20 nM of siRNA using Lipofectamine 2000 reagent (Invitrogen). Cell proliferation was assessed by MTS assay ([3-(4,5-dimethyl-2-yl)-5-(3-carboxymethoxyphenyl)-2-(4-sulfophenyl)-2H-tetrazolium, inner salt; Promega, Madison, WI, USA), according to the manufacturer's protocol, 72 hours posttransfection.
Tyrosine phosphorylation profiling of mouse mammary tumour models
Because the mammary gland develops postnatally, murine mammary epithelium can be transferred into the surgically cleared fat pads of juvenile isogenic naïve recipients. The use of GEM as donors, or viral transduction of cells after extraction, allows the rapid generation of chimeric models in which implanted transgenic or gene-knockout mammary epithelium develops in an otherwise wild-type host ,. Exploiting this strategy, we performed phosphoproteomic profiling using a combined immunoaffinity/LC-MS/MS approach across multiple MSCV-Neu (HER2), Tp53-null (p53) and PyMT tumours, as illustrated in Figure 1A. The lineage of each tumour from its original clone is outlined in Additional file 1: Figure S1A. In total, 856 unique pY sites were identified and quantitated from the tumour sample population with a FDR <1%. Following high-confidence localisation site filtering at ≥0.75, 763 unique pY sites from 506 unique phosphoproteins remained, with approximately 63% being common across all tumour types (Figure 1B and Additional file 2: Table S1).
These pY sites mapped to 506 unique proteins. These proteins, which were grouped based on Gene Ontology functional categories, are illustrated in Additional file 3: Figure S2A. Additional file 3: Figure S2B shows an expanded section of the kinases identified in the three tumour types. The largest proportion of kinases comprised serine/threonine kinases, followed by non-RTKs and RTKs. The least represented group comprised the dual-specificity kinases. The total pY site intensities were determined for the non-RTKs, which are shown in Additional file 3: Figure S2C. The most intense signal from the pY sites of non-RTKs was derived from SFKs (Yes, Fyn, Lyn and Hck), followed by focal adhesion kinase 1. When this approach was applied to RTKs, the most intense signal was derived from Met, followed by EphA2 (Additional file 3: Figure S2D).
To compare tyrosine phosphorylation patterns between the different tumour types and identify potential similarities, we undertook unsupervised hierarchical clustering. This revealed that two p53 tumour samples formed a distinct subgroup characterised by markedly elevated tyrosine phosphorylation of many pY sites (Additional file 4: Figure S3). The remaining tumour samples were distributed over a second major subgroup and exhibited a complex pattern of relationships. Of note, five PyMT tumours were closely related, and the majority of the p53 and HER2 tumours clustered together (Additional file 4: Figure S3). Also, although some tumours of shared origin, such as p53 tumours 1203 and 1204, clustered together, others, such as HER2 tumours HER 1197 and 3051, were well separated (Additional file 1: Figure S1B), indicating that divergence can occur following serial transplantation.
To generate a list of pY sites characteristic of each tumour type, we first identified pY sites detected in at least 75% of the corresponding samples. Therefore, we selected pY sites that were identified in at least (1) six of eight p53 tumour samples, (2) six of eight PyMT tumour samples or (3) four of five Her2 tumour samples. Following this selection-filtering process, 381 pY sites remained (Figure 1C and Additional file 5: Table S2). These sites mapped to 256 proteins, with 177 of them being found in all three tumour types. Upon unsupervised hierarchical clustering of this refined data set, all of the PyMT tumours were clustered together (Additional file 6: Figure S4). However, p53 tumours 1203 and 1204 still formed a distinct subgroup, and the HER2 tumours were divided into two groups, one related to PyMT tumours and the other to p53 tumours (Additional file 6: Figure S4).
To further refine the list of sites upregulated in a particular tumour type versus the other two, an analysis of variance (ANOVA)–based approach was implemented. Thirty-two pY sites were upregulated in HER2 tumours, eighty-six in p53 tumours and eighty in PyMT tumours (Additional file 5: Table S2). Of note, certain proteins exhibited site-selective variation in phosphorylation between different tumour types, examples being Tln1 and Cbl. Elevated levels of Erbb3 pY1325 and p85 pY467 in PyMT tumours were confirmed by Western blot analysis with appropriate phosphospecific antibodies (Additional file 7: Figure S5). Immunoblotting for selected targets revealed that, for certain pY sites characteristic of a given subtype, elevated phosphorylation was accompanied by increased expression of the given protein. This was exemplified by Irs1 and Pik3r1 (p85) in PyMT tumours (Figure 2A and B). However, for others, the increased phosphorylation must reflect a change in relative phosphorylation, as for Erbb3 in PyMT tumours, Cdk16 in p53 tumours and Erffi1 in HER2 tumours (Figure 2C to E). References for the sites characterised in this figure are provided in PhosphoSite .
We also developed and applied a random forest classifier to identify 17 highly variable pY sites that can be used to distinguish between the tumour types (Figure 3). Amongst these, elevated phosphorylation of Pik3r1, Pik3r3, a G protein–coupled receptor (Gprc5c) and an uncharacterised protein (2310030G06Rik) was specific to PyMT tumours, and p53 tumours were characterised by increased phosphorylation of Vim, Ptrf, Vcl, Fgr, Met and a protein related to heat shock protein 70 kDa. Elevated tyrosine phosphorylation of Errfi1 characterised HER2 tumours.
Pathway and protein–protein interaction analysis
Lists of pY sites enriched in each tumour type were then subjected to bioinformatics analysis to identify associated pathways or biological processes (Table 1). Strikingly, marked differences were observed between p53, HER2 and PyMT tumours. The former exhibited an enrichment of processes linked to cytoskeletal reorganisation, with the top five significant terms being leucocyte transendothelial migration, bacterial invasion of epithelial cells, regulation of actin cytoskeleton, focal adhesion and muscle contraction. The only term significantly enriched for HER2 tumours was the ErbB signalling pathway, whereas the PyMT tumours were characterised by processes involving PI3K signalling, including type 2 diabetes mellitus, the FcεRI signalling pathway, acute myeloid leukaemia, Fcγ receptor–mediated phagocytosis and the neurotrophin signalling pathway.
Significantly upregulated pY sites (determined by ANOVA) characteristic of each tumour type were analysed using PINA , (Figures 4, 5 and 6). Where available, consensus lists were searched against mouse data sets containing previously reported interaction data; however, in cases where interactions were not reported in mouse data sets, human orthologues were used to model potential interactions. In addition, kinase substrate relationships were extracted from PhosphoSitePlus . A key characteristic of p53 tumours was the presence of a major protein–protein interaction hub composed of several proteins involved in cytoskeletal organisation (Actn1, Actn4, Sept7, Sept 9 and Plec), as well as a further hub centred on Itgb1 (Figure 4). The latter involved numerous proteins localised to focal adhesions, including Bcar1, Tln1, Pxn and Vcl. These findings are in agreement with the corresponding pathway analysis (Table 1). Furthermore, there are notable similarities between this network and that present in human basal and claudin-low breast cancer cell lines , including increased tyrosine phosphorylation of Met, Cav1, Bcar1 and Tln1. This is consistent with previous reports that subsets of p53 tumours exhibit gene expression signatures characteristic of these breast cancer subgroups ,. In contrast, HER2 tumours exhibited a protein–protein interaction network that reflected the relationship between HER2 and its associating substrates (for example, Cbl, Ptpn11, Plcγ1), but it is noteworthy that Egfr and Pdgfrα, as well as the feedback inhibitor Errfi1, are also components of this network (Figure 5). Consistent with the pathway analysis, PyMT tumours exhibit a major interaction hub that centres on the PI3K regulatory subunit Pik3r1 (Figure 6). Interestingly, in addition to PyMT itself, this hub involves two RTKs (Kit and Erbb3) and two docking proteins (Irs1 and Gab1) that are well-characterised binding partners for Pik3r1. Although Src is required for efficient induction of mammary tumours by PyMT , increased phosphorylation of Y416 on the Src activation loop (a peptide sequence conserved amongst SFKs) or Src-selective peptides was not detected in PyMT tumours. However, three additional SFKs—Yes, Lyn and Hck—exhibited increased phosphorylation on unique peptide sequences and resided within the network (Figure 6). The former is known to interact with PyMT  but is not essential for PyMT-induced mammary tumourigenesis .
Subclassification of p53 tumours
In a previous study on the Tp53-null mouse mammary tumour model, researchers demonstrated that the generated tumours exhibit considerable heterogeneity and can be classified by transcript profiling into a variety of molecular subtypes, including one luminal, one claudin-low and two basal-like subtypes . To interrogate this potential heterogeneity at the level of tyrosine phosphorylation patterns, we undertook unsupervised hierarchical clustering of the p53 tumours based on all phosphorylation sites detected in these tumours (Additional file 8: Figure S6). Tumours 1203 and 1204 clustered together and were characterised by markedly elevated (greater than tenfold) tyrosine phosphorylation of approximately 100 proteins compared to the remainder of p53 tumours (Additional file 9: Table S3). Tumour 14845 also exhibited relatively high tyrosine phosphorylation of many proteins, but not to the same extent as tumours 1203 and 1204, and tumours 3049, 1201 and 1202 displayed elevated phosphorylation of more distinct protein subsets. Tumours 1205 and 1206 were characterised by relatively low levels of tyrosine phosphorylation. Hypothesizing that the extremely high levels of tyrosine phosphorylation in tumours 1203 and 1204 might reflect upregulation of a particular tyrosine kinase in these two tumours, we noted that several of the upregulated pY sites were from Met, specifically Y1228, Y1232, Y1233, Y1347 and Y1354 (Figure 7). In addition, Y6 and Y14 on Cav1 also exhibited increased phosphorylation, which was of interest, given the close proximity of the Met and Cav1 genes on mouse chromosome 6 and coamplification of these two genes in particular tumour types -. Immunoblotting for Met and Cav1 revealed markedly increased total protein expression in the p53 samples from tumours 1203 and 1204 (Figure 8A and B). In addition, aCGH was performed on the p53 tumours. This revealed that, though p53 tumours 1201, 1205, 1206 and 14845 exhibited small (two- to threefold) increases in copy number for regions of chromosome 6 harbouring Met and Cav1, tumours 1203 and 1204 displayed 20- and 25-fold amplifications, respectively, of a chromosomal segment spanning these two genes (Figure 8C and D and Additional file 10: Table S4).
Having demonstrated that tumours 1203 and 1204 are characterised by Met gene amplification, we compared the tyrosine phosphorylation profile of these two tumours versus the other p53-null tumours in more detail. Consistent with enhanced Met signalling, tyrosine phosphorylation of Hgs, which regulates Met endocytic trafficking , was markedly elevated (Additional file 9: Table S3). These tumours were also characterised by increased tyrosine phosphorylation of Ptrf/cavin, which is localised with Cav1 to caveolae, and of the RTK Axl and the docking protein Dock1 (Additional file 9: Table S3). Interestingly, pathway analysis revealed that tumours 1203/1204 demonstrated significant enrichment for the term glycolysis/gluconeogenesis, reflecting increased tyrosine phosphorylation of several metabolic enzymes, including Pgam1, Pkm2 and Pgk1 (Additional file 11: Table S5).
To gain insight into the functional role of Met in the p53 tumours, we established cell lines from tumours 1201 and 1206 (low-level Met gene amplification) and tumour 1204 (high-level Met amplification). Pretreatment of these cell lines with the Met TKI PHA-665752 decreased Met phosphorylation on Y1234/Y1235 and markedly reduced the elevated cellular tyrosine phosphorylation levels in tumour 1204 cells (Figure 9A). Importantly, though administration of the Met TKI at a concentration (1 μM) that ablated Met tyrosine phosphorylation had no effect on proliferation in tumour 1201 cells (Figure 9B), this treatment markedly attenuated proliferation of tumour 1204 cells (Figure 9C). As an alternative and complementary approach, we knocked down Met expression using selective siRNAs (Additional file 12: Figure S7). Met knockdown had only a modest effect in the 1201 and 1206 tumour cell lines, but markedly attenuated proliferation in the 1204 tumour cell line with high-level Met amplification (P < 0.005 for comparisons of Met knockdowns in the 1201 or 1206 tumour cell line versus the 1204 tumour cell line). These data indicate that the extremely high levels of tyrosine phosphorylation detected in tumour 1204 (Figure 7) are dependent on Met activity, and that cells from these tumours are addicted to Met activation for their proliferation in vitro.
Interrogation of the tyrosine phosphorylation patterns characteristic of the other p53-null subgroups revealed that phosphorylation of the atypical kinase Peak1/SgK269  was elevated in tumour 14845 and to a lesser extent the subgroup containing 3049, 1202 and 1201, whereas phosphorylation of Prex2 characterised tumours 14845 and 3049, where it was accompanied by enhanced phosphorylation of Pik3r2 (Figure 7). Taken in association with our demonstration of Met amplification in tumours 1203 and 1204, these findings highlight the power of tyrosine phosphorylation profiling to identify candidate oncogenic drivers .
Specific genetically modified mouse models of breast cancer have made major contributions to our understanding of the role of particular oncogenes in the development of this malignancy . In addition, characterisation of tumours arising in these models by approaches such as transcript profiling and aCGH have identified molecular aberrations conserved between mouse and human mammary tumours that are likely to represent key 'driver' events ,,. However, to date, no comprehensive characterisation of the tyrosine kinase signalling networks present in such models has yet been undertaken. In this study, we used MS-based global tyrosine phosphorylation profiling to characterise three commonly used mouse models of breast cancer, which revealed distinct differences between these models and features conserved with subtypes of the human disease.
For an erbB2 model, we used a truncated form of Neu , that models a form of erbB2 lacking the extracellular domain present in approximately 30% of HER2-positive breast cancer patients, for whom it is associated with a worse prognosis . This form exhibits potent oncogenic activity . These tumours were characterised by increased phosphorylation of several signalling molecules known to associate with erbB2, including Ptpn11, Plcγ1, Cbl and CrkL, as well as its heterodimerisation partner Egfr. Given the known heterodimerisation of full-length erbB2 and erbB3, our finding that phosphorylation of erbB3 Y1325 was not elevated in the majority of HER2 tumours (expressing truncated neu), despite increased erbB3 expression (Figure 2), is surprising and warrants further investigation in cell culture models and clinical specimens. Additional novel findings include enhanced tyrosine phosphorylation of Errfi1 and Pdgfrα. The former represents a negative feedback regulator of erbB signalling that inhibits EGFR and erbB2 kinase activity by binding directly to the receptor kinase domain and blocking the formation of the activating dimer interface . However, although increased expression of Errfi1 is associated with high erbB receptor signalling , expression of Errfi1 was similar between HER2 and p53 tumours, indicating that it is the relative tyrosine phosphorylation of Errfi1 on Y393/Y394 that is altered. Importantly, phosphorylation of Errfi1 on Y394 reduces its ability to inhibit Egfr , indicating that this modification may act to attenuate negative regulation of HER2 signalling by Errfi1 in these tumours. Consequently, it will be of great interest to determine how expression of Errfi1, and also its tyrosine phosphorylation, relate to HER2 expression and outcome in breast cancer patients.
The finding of increased Pdgfrα phosphorylation is of interest because autocrine Pdgf action is associated with TGF-β-induced epithelial–mesenchymal transition (EMT) in MMTV-Neu mammary tumours and because increased Pdgfrα and Pdgfrβ expression occur in late-stage, invasive human breast cancers . Researchers in a previous study demonstrated that a truncated form of HER2, similar to that used in this study, specifically induced the expression of a set of genes involved in cancer cell spread and associated with poor patient prognosis . Thus, the induction of Pdgfrα signalling in our model may reflect the oncogenic potency of the truncated form of Neu utilised and the occurrence of EMT.
A previous study on the PyMT model revealed a dependence on PI3K activation for efficient tumourigenesis . Although that finding was revealed by mutation of the p85 binding sites on PyMT , our pathway and protein–protein interaction analyses on the tyrosine-phosphorylated proteins selectively enriched in these tumours also support a key role for PI3 kinase signalling in this tumour type. These results highlight the presence of other p85 binding proteins in these tumours in addition to PyMT itself, including erbB3, Gab1 and Irs1. All of the latter three proteins have been characterised as positive regulators of PI3K signalling . In addition, enhanced tyrosine phosphorylation of p85a on several sites was detected in these tumours. Tyrosine phosphorylation of p85a on Y688 has been reported to relieve its inhibitory activity on PI3K ; however, we did not detect phosphorylation on this site in our present study, and the functional role of the p85a phosphorylation events enhanced in PyMT tumours requires further characterisation.
An important question is the identity of the tyrosine kinases driving PyMT tyrosine phosphorylation. It is well established that PyMT associates with, and is phosphorylated by, the SFKs Src and Yes, and also that expression of Src is essential for efficient PyMT-induced tumourigenesis ,. However, in a recent article, researchers reported that PyMT associates with the RTK erbB2 and the catalytically impaired erbB3 and that both erbB2 activity and erbB3 expression are essential for PyMT transforming activity . These data, together with those published by the Muller group , indicate that erbB2 may act upstream of Src in this model. However, we could detect only erbB3-derived, and not erbB2-derived, tyrosine-phosphorylated peptides in the majority of PyMT tumours. A potential explanation is that one or more RTKs may substitute for erbB2 in this model. Because Egfr phosphorylation was detected and phosphorylation of Kit was enhanced in this tumour type, these RTKs represent potential candidates.
The p53 tumour model exhibits heterogeneity at the histological and gene expression levels, and, on the basis of the latter parameter, can be classified into different molecular subgroups, including basal-like, claudin-low and luminal subtypes . Although we profiled only eight p53 tumours, our data are consistent with the ability of this model to generate tumours with a basal or claudin-low phenotype. First, there was notable overlap between the identified signalling networks in these tumours and basal and claudin-low breast cancer cells . Second, tyrosine-phosphorylated proteins characteristic of p53 tumours exhibited an enrichment for pathways associated with cytoskeletal organisation, in concordance with the observed mutational spectrum reported for triple-negative breast cancers . Third, two of the p53 tumours exhibited markedly enhanced tyrosine phosphorylation of many cellular proteins, including Met, and this was associated with high-level amplification of a segment of chromosome 6 that contains the Met and Cav1 genes. This is consistent with high expression and amplification of the MET gene in human basal breast cancers ,, and amplification of the Met/Cav1 locus in basal-like mammary tumours in mice with mammary-specific deletion of Lfng . In addition, researchers in a recent study demonstrated that transgenic Met overexpression in the mouse mammary gland cooperates with conditional p53 loss to drive the formation of claudin-low mammary tumours, and they showed that claudin-low tumours developing in mice without the MMTV-Met transgene and only conditional mammary gland-specific deletion of p53 exhibited amplification of the endogenous Met gene . These data are in concordance with our data demonstrating Met amplification in a subset of p53-null mammary tumours and highlight the potential utility of transplantable tumours 1203 and 1204 and the MMTV-Metmt;Trp53fl/+;Cre model generated by the Park group  as preclinical models in which to test Met-directed therapies for triple-negative breast cancer. Indeed, our data indicate that p53-null cells with high-level Met amplification are dependent on Met activation for their continued proliferation.
Our finding that tumours 1203 and 1204 demonstrated increased tyrosine phosphorylation of several enzymes involved in glycolysis and gluconeogenesis indicates that this tumour subset is likely to exhibit altered metabolic properties. Indeed, tyrosine phosphorylation of both Pkm2 Y105 and Pgam1 Y26 regulates their activity and, by modulating aerobic glycolysis and anabolic biosynthesis, promotes tumourigenesis ,. In light of these findings, it will be important to determine whether MET-amplified human breast cancers exhibit corresponding metabolic changes and whether the altered metabolism confers an 'Achilles heel' amenable to therapeutic intervention.
The identification of two p53 tumours exhibiting Met gene amplification highlights how mouse models of mammary tumourigenesis can be used to identify breast cancer oncogenic 'drivers'. Of note, our phosphoproteomic profiling of p53 tumours identified markedly enhanced tyrosine phosphorylation of Peak1/SgK269 and Prex2 in other p53-null tumour subsets. The former protein is an atypical kinase recently demonstrated to exhibit characteristics of a breast cancer oncogene, and site-selective Peak1/SgK269 tyrosine phosphorylation promotes biological activity . Prex2 is a negative regulator of PTEN (phosphatase and tensin homologue) and is mutated in melanoma , although the role of Prex2 tyrosine phosphorylation is yet to be determined. Consequently, these two proteins represent strong candidates for further functional characterisation and evaluation as therapeutic targets and biomarkers.
Our finding that three commonly utilised mouse models of breast cancer are characterised by distinct tyrosine phosphorylation–based signalling networks provides important insights into mechanisms of mammary tumour development and progression. The identification of Met as an oncogenic driver in a subset of p53-null tumours highlights how comparative phosphoproteomics can identify conserved oncogenic signalling pathways and identifies a new preclinical model for a subset of triple-negative breast cancers.
Array comparative genomic hybridisation
Analysis of variance
False discovery rate
Genetically engineered mouse
G protein–coupled receptor family C, group 5, member C
- k-NN k:
Mouse mammary tumour virus
Tandem mass spectrometry
Nano-liquid chromatography tandem mass spectrometry
Protein interaction network analysis
Polyoma virus middle T antigen
Receptor tyrosine kinase
Src family kinase
Yarden Y, Pines G: The ERBB network: at last, cancer therapy meets systems biology. Nat Rev Cancer. 2012, 12: 553-563. 10.1038/nrc3309.
Kalinsky K, Jacks LM, Heguy A, Patil S, Drobnjak M, Bhanot UK, Hedvat CV, Traina TA, Solit D, Gerald W, Moynahan ME: PIK3CA mutation associates with improved outcome in breast cancer. Clin Cancer Res. 2009, 15: 5049-5059. 10.1158/1078-0432.CCR-09-0632.
Creedon H, Brunton VG: Src kinase inhibitors: promising cancer therapeutics?. Crit Rev Oncog. 2012, 17: 145-159. 10.1615/CritRevOncog.v17.i2.20.
Celebi JT, Adams SJ, Aydin IT: GAB2—a scaffolding protein in cancer. Mol Cancer Res. 2012, 10: 1265-1270. 10.1158/1541-7786.MCR-12-0352.
Banerji S, Cibulskis K, Rangel-Escareno C, Brown KK, Carter SL, Frederick AM, Lawrence MS, Sivachenko AY, Sougnez C, Zou L, Cortes ML, Fernandez-Lopez JC, Peng S, Ardlie KG, Auclair D, Bautista-Piña V, Duke F, Francis J, Jung J, Maffuz-Aziz A, Onofrio RC, Parkin M, Pho NH, Quintanar-Jurado V, Ramos AH, Rebollar-Vega R, Rodriguez-Cuevas S, Romero-Cordoba SL, Schumacher SE, Stransky N, et al: Sequence analysis of mutations and translocations across breast cancer subtypes. Nature. 2012, 486: 405-409. 10.1038/nature11154.
Bentires-Alj M, Gil SG, Chan R, Wang ZC, Wang Y, Imanaka N, Harris LN, Richardson A, Neel BG, Gu H: A role for the scaffolding adapter GAB2 in breast cancer. Nat Med. 2006, 12: 114-121. 10.1038/nm1341.
Brough R, Frankum JR, Sims D, Mackay A, Mendes-Pereira AM, Bajrami I, Costa-Cabral S, Rafiq R, Ahmad AS, Cerone MA, Natrajan R, Sharpe R, Shiu KK, Wetterskog D, Dedes KJ, Lambros MB, Rawjee T, Linardopoulos S, Reis-Filho JS, Turner NC, Lord CJ, Ashworth A: Functional viability profiles of breast cancer. Cancer Discov. 2011, 1: 260-273. 10.1158/2159-8290.CD-11-0107.
Ursini-Siegel J, Schade B, Cardiff RD, Muller WJ: Insights from transgenic mouse models of ERBB2-induced breast cancer. Nat Rev Cancer. 2007, 7: 389-397. 10.1038/nrc2127.
Yuan W, Stawiski E, Janakiraman V, Chan E, Durinck S, Edgar KA, Kljavin NM, Rivers CS, Gnad F, Roose-Girma M, Haverty PM, Fedorowicz G, Heldens S, Soriano RH, Zhang Z, Wallin JJ, Johnson L, Merchant M, Modrusan Z, Stern HM, Seshagiri S: Conditional activation of Pik3caH1047R in a knock-in mouse model promotes mammary tumorigenesis and emergence of mutations. Oncogene. 2012, 32: 318-326. 10.1038/onc.2012.53.
Hudis CA: Trastuzumab—mechanism of action and use in clinical practice. N Engl J Med. 2007, 357: 39-51. 10.1056/NEJMra043186.
Brenton JD, Carey LA, Ahmed AA, Caldas C: Molecular classification and molecular forecasting of breast cancer: ready for clinical application?. J Clin Oncol. 2005, 23: 7350-7360. 10.1200/JCO.2005.03.3845.
Rakha EA, Reis-Filho JS, Ellis IO: Basal-like breast cancer: a critical review. J Clin Oncol. 2008, 26: 2568-2581. 10.1200/JCO.2007.13.1748.
Graveel CR, DeGroot JD, Su Y, Koeman J, Dykema K, Leung S, Snider J, Davies SR, Swiatek PJ, Cottingham S, Watson MA, Ellis MJ, Sigler RE, Furge KA, Vande Woude GF: Met induces diverse mammary carcinomas in mice and is associated with human basal breast cancer. Proc Natl Acad Sci U S A. 2009, 106: 12909-12914. 10.1073/pnas.0810403106.
Nielsen TO, Hsu FD, Jensen K, Cheang M, Karaca G, Hu Z, Hernandez-Boussard T, Livasy C, Cowan D, Dressler L, Akslen LA, Ragaz J, Gown AM, Gilks CB, van de Rijn M, Perou CM: Immunohistochemical and clinical characterization of the basal-like subtype of invasive breast carcinoma. Clin Cancer Res. 2004, 10: 5367-5374. 10.1158/1078-0432.CCR-04-0220.
Ponzo MG, Lesurf R, Petkiewicz S, O'Malley FP, Pinnaduwage D, Andrulis IL, Bull SB, Chughtai N, Zuo D, Souleimanova M, Germain D, Omeroglu A, Cardiff RD, Hallett M, Park M: Met induces mammary tumors with diverse histologies and is associated with poor outcome and human basal breast cancer. Proc Natl Acad Sci U S A. 2009, 106: 12903-12908. 10.1073/pnas.0810402106.
Hochgräfe F, Zhang L, O'Toole SA, Browne BC, Pinese M, Porta Cubas A, Lehrbach GM, Croucher DR, Rickwood D, Boulghourjian A, Shearer R, Nair R, Swarbrick A, Faratian D, Mullen P, Harrison DJ, Biankin AV, Sutherland RL, Raftery MJ, Daly RJ: Tyrosine phosphorylation profiling reveals the signaling network characteristics of basal breast cancer cells. Cancer Res. 2010, 70: 9391-9401. 10.1158/0008-5472.CAN-10-0911.
Grigoriadis A, Mackay A, Noel E, Wu PJ, Natrajan R, Frankum J, Reis-Filho JS, Tutt A: Molecular characterisation of cell line models for triple-negative breast cancers. BMC Genomics. 2012, 13: 619-10.1186/1471-2164-13-619.
Vargo-Gogola T, Rosen JM: Modelling breast cancer: one size does not fit all. Nat Rev Cancer. 2007, 7: 659-672. 10.1038/nrc2193.
Fluck MM, Schaffhausen BS: Lessons in signaling and tumorigenesis from polyomavirus middle T antigen. Microbiol Mol Biol Rev. 2009, 73: 542-563. 10.1128/MMBR.00009-09.
Webster MA, Hutchinson JN, Rauh MJ, Muthuswamy SK, Anton M, Tortorice CG, Cardiff RD, Graham FL, Hassell JA, Muller WJ: Requirement for both Shc and phosphatidylinositol 3’ kinase signaling pathways in polyomavirus middle T-mediated mammary tumorigenesis. Mol Cell Biol. 1998, 18: 2344-2359.
Herschkowitz JI, Simin K, Weigman VJ, Mikaelian I, Usary J, Hu Z, Rasmussen KE, Jones LP, Assefnia S, Chandrasekharan S, Backlund MG, Yin Y, Khramtsov AI, Bastein R, Quackenbush J, Glazer RI, Brown PH, Green JE, Kopelovich L, Furth PA, Palazzo JP, Olopade OI, Bernard PS, Churchill GA, Van Dyke T, Perou CM: Identification of conserved gene expression features between murine mammary carcinoma models and human breast tumors. Genome Biol. 2007, 8: R76-10.1186/gb-2007-8-5-r76.
Herschkowitz JI, Zhao W, Zhang M, Usary J, Murrow G, Edwards D, Knezevic J, Greene SB, Darr D, Troester MA, Hilsenbeck SG, Medina D, Perou CM, Rosen JM: Comparative oncogenomics identifies breast tumors enriched in functional tumor-initiating cells. Proc Natl Acad Sci U S A. 2012, 109: 2778-2783. 10.1073/pnas.1018862108.
Jerry DJ, Kittrell FS, Kuperwasser C, Laucirica R, Dickinson ES, Bonilla PJ, Butel JS, Medina D: A mammary-specific model demonstrates the role of the p53 tumor suppressor gene in tumor development. Oncogene. 2000, 19: 1052-1058. 10.1038/sj.onc.1203270.
Bennett CN, Green JE: Unlocking the power of cross-species genomic analyses: identification of evolutionarily conserved breast cancer networks and validation of preclinical models. Breast Cancer Res. 2008, 10: 213-10.1186/bcr2125.
Bennett CN, Green JE: Genomic analyses as a guide to target identification and preclinical testing of mouse models of breast cancer. Toxicol Pathol. 2010, 38: 88-95. 10.1177/0192623309357074.
Nair R, Roden DL, Teo WS, McFarland A, Junankar S, Ye S, Nguyen A, Yang J, Nikolic I, Hui M, Morey A, Shah J, Pfefferle AD, Usary J, Selinger C, Baker LA, Armstrong N, Cowley MJ, Naylor MJ, Ormandy CJ, Lakhani SR, Herschkowitz JI, Perou CM, Kaplan W, O'Toole SA, Swarbrick A: c-Myc and Her2 cooperate to drive a stem-like phenotype with poor prognosis in breast cancer. Oncogene. 2013, 33: 3992-4002. 10.1038/onc.2013.368.
Bargmann CI, Weinberg RA: Increased tyrosine kinase activity associated with the protein encoded by the activated neu oncogene. Proc Natl Acad Sci U S A. 1988, 85: 5394-5398. 10.1073/pnas.85.15.5394.
Van Parijs L, Refaeli Y, Lord JD, Nelson BH, Abbas AK, Baltimore D: Uncoupling IL-2 signals that regulate T cell proliferation, survival, and Fas-mediated activation-induced cell death. Immunity. 1999, 11: 281-288. 10.1016/S1074-7613(00)80103-X.
Swarbrick A, Roy E, Allen T, Bishop JM: Id1 cooperates with oncogenic Ras to induce metastatic mammary carcinoma by subversion of the cellular senescence response. Proc Natl Acad Sci U S A. 2008, 105: 5402-5407. 10.1073/pnas.0801505105.
Welm AL, Kim S, Welm BE, Bishop JM: MET and MYC cooperate in mammary tumorigenesis. Proc Natl Acad Sci U S A. 2005, 102: 4324-4329. 10.1073/pnas.0500470102.
Zhang M, Behbod F, Atkinson RL, Landis MD, Kittrell F, Edwards D, Medina D, Tsimelzon A, Hilsenbeck S, Green JE, Michalowska AM, Rosen JM: Identification of tumor-initiating cells in a p53-null mouse model of breast cancer. Cancer Res. 2008, 68: 4674-4682. 10.1158/0008-5472.CAN-07-6353.
Guy CT, Cardiff RD, Muller WJ: Induction of mammary tumors by expression of polyomavirus middle T oncogene: a transgenic mouse model for metastatic disease. Mol Cell Biol. 1992, 12: 954-961.
Cox J, Neuhauser N, Michalski A, Scheltema RA, Olsen JV, Mann M: Andromeda: a peptide search engine integrated into the MaxQuant environment. J Proteome Res. 2011, 10: 1794-1805. 10.1021/pr101065j.
Troyanskaya O, Cantor M, Sherlock G, Brown P, Hastie T, Tibshirani R, Botstein D, Altman RB: Missing value estimation methods for DNA microarrays. Bioinformatics. 2001, 17: 520-525. 10.1093/bioinformatics/17.6.520.
Strimmer K: A unified approach to false discovery rate estimation. BMC Bioinformatics. 2008, 9: 303-10.1186/1471-2105-9-303.
Breiman L: Random forests. Mach Learn. 2001, 45: 5-32. 10.1023/A:1010933404324.
Cowley MJ, Pinese M, Kassahn KS, Waddell N, Pearson JV, Grimmond SM, Biankin AV, Hautaniemi S, Wu J: PINA v2.0: mining interactome modules. Nucleic Acids Res. 2012, 40: D862-D865. 10.1093/nar/gkr967.
Wu J, Vallenius T, Ovaska K, Westermarck J, Mäkelä TP, Hautaniemi S: Integrated network analysis platform for protein-protein interactions. Nat Methods. 2009, 6: 75-77. 10.1038/nmeth.1282.
Hornbeck PV, Kornhauser JM, Tkachev S, Zhang B, Skrzypek E, Murray B, Latham V, Sullivan M: PhosphoSitePlus: a comprehensive resource for investigating the structure and function of experimentally determined post-translational modifications in man and mouse. Nucleic Acids Res. 2012, 40: D261-D270. 10.1093/nar/gkr1122.
Shi G, Zhang L, Jiang T: MSOAR 2.0: incorporating tandem duplications into ortholog assignment based on genome rearrangement. BMC Bioinformatics. 2010, 11: 10-10.1186/1471-2105-11-10.
Wu J, Mao X, Cai T, Luo J, Wei L: KOBAS server: a web-based platform for automated annotation and pathway identification. Nucleic Acids Res. 2006, 34: W720-W724. 10.1093/nar/gkl167.
Xie C, Mao X, Huang J, Ding Y, Wu J, Dong S, Kong L, Gao G, Li CY, Wei L: KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases. Nucleic Acids Res. 2011, 39: W316-W322. 10.1093/nar/gkr483.
Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Series B Methodol. 1995, 57: 289-300.
Janes PW, Daly RJ, deFazio A, Sutherland RL: Activation of the Ras signalling pathway in human breast cancer cells overexpressing erbB-2. Oncogene. 1994, 9: 3601-3608.
Olshen AB, Venkatraman ES, Lucito R, Wigler M: Circular binary segmentation for the analysis of array-based DNA copy number data. Biostatistics. 2004, 5: 557-572. 10.1093/biostatistics/kxh008.
Bioconductor Open Source Software for Bioinformatics: DNAcopy: DNA copy number data analysis.  (accessed 15 September 2014)., [http://www.bioconductor.org/packages/release/bioc/html/DNAcopy.html]
R Project for Statistical Computing: R-2.15.2.  (accessed 15 September 2014)., [http://cran.r-project.org/bin/windows/base/old/2.15.2/]
Brummer T, Schramek D, Hayes VM, Bennett HL, Caldon CE, Musgrove EA, Daly RJ: Increased proliferation and altered growth factor dependence of human mammary epithelial cells overexpressing the Gab2 docking protein. J Biol Chem. 2006, 281: 626-637. 10.1074/jbc.M509567200.
PhosphoSitePlus.  (accessed 15 September 2014)., [http://www.phosphosite.org/]
Guy CT, Muthuswamy SK, Cardiff RD, Soriano P, Muller WJ: Activation of the c-Src tyrosine kinase is required for the induction of mammary tumors in transgenic mice. Genes Dev. 1994, 8: 23-32. 10.1101/gad.8.1.23.
Hellman A, Zlotorynski E, Scherer SW, Cheung J, Vincent JB, Smith DI, Trakhtenbrot L, Kerem B: A role for common fragile site induction in amplification of human oncogenes. Cancer Cell. 2002, 1: 89-97. 10.1016/S1535-6108(02)00017-X.
Xu K, Usary J, Kousis PC, Prat A, Wang DY, Adams JR, Wang W, Loch AJ, Deng T, Zhao W, Cardiff RD, Yoon K, Gaiano N, Ling V, Beyene J, Zacksenhaus E, Gridley T, Leong WL, Guidos CJ, Perou CM, Egan SE: Lunatic fringe deficiency cooperates with the Met/Caveolin gene amplicon to induce basal-like breast cancer. Cancer Cell. 2012, 21: 626-641. 10.1016/j.ccr.2012.03.041.
Smolen GA, Muir B, Mohapatra G, Barmettler A, Kim WJ, Rivera MN, Haserlat SM, Okimoto RA, Kwak E, Dahiya S, Garber JE, Bell DW, Sgroi DC, Chin L, Deng CX, Haber DA: Frequent met oncogene amplification in a Brca1/Trp53 mouse model of mammary tumorigenesis. Cancer Res. 2006, 66: 3452-3455. 10.1158/0008-5472.CAN-05-4181.
Abella JV, Peschard P, Naujokas MA, Lin T, Saucier C, Urbé S, Park M: Met/hepatocyte growth factor receptor ubiquitination suppresses transformation and is required for Hrs phosphorylation. Mol Cell Biol. 2005, 25: 9632-9645. 10.1128/MCB.25.21.9632-9645.2005.
Croucher DR, Hochgräfe F, Zhang L, Liu L, Lyons RJ, Rickwood D, Tactacan CM, Browne BC, Ali N, Chan H, Shearer R, Gallego-Ortega D, Saunders DN, Swarbrick A, Daly RJ: Involvement of Lyn and the atypical kinase SgK269/PEAK1 in a basal breast cancer signaling pathway. Cancer Res. 2013, 73: 1969-1980. 10.1158/0008-5472.CAN-12-1472.
Rikova K, Guo A, Zeng Q, Possemato A, Yu J, Haack H, Nardone J, Lee K, Reeves C, Li Y, Hu Y, Tan Z, Stokes M, Sullivan L, Mitchell J, Wetzel R, MacNeill J, Ren JM, Yuan J, Bakalarski CE, Villen J, Kornhauser JM, Smith B, Li D, Zhou X, Gygi SP, Gu TL, Polakiewicz RD, Rush J, Comb MJ: Global survey of phosphotyrosine signaling identifies oncogenic kinases in lung cancer. Cell. 2007, 131: 1190-1203. 10.1016/j.cell.2007.11.025.
Hodgson JG, Malek T, Bornstein S, Hariono S, Ginzinger DG, Muller WJ, Gray JW: Copy number aberrations in mouse breast tumors reveal loci and genes important in tumorigenic receptor tyrosine kinase signaling. Cancer Res. 2005, 65: 9695-9704. 10.1158/0008-5472.CAN-05-0755.
Sáez R, Molina MA, Ramsey EE, Rojo F, Keenan EJ, Albanell J, Lluch A, García-Conde J, Baselga J, Clinton GM: p95HER-2 predicts worse outcome in patients with HER-2-positive breast cancer. Clin Cancer Res. 2006, 12: 424-431. 10.1158/1078-0432.CCR-05-1807.
Pedersen K, Angelini PD, Laos S, Bach-Faig A, Cunningham MP, Ferrer-Ramón C, Luque-García A, García-Castillo J, Parra-Palau JL, Scaltriti M, Ramón y Cajal S, Baselga J, Arribas J: A naturally occurring HER2 carboxy-terminal fragment promotes mammary tumor growth and metastasis. Mol Cell Biol. 2009, 29: 3319-3331. 10.1128/MCB.01803-08.
Zhang X, Pickin KA, Bose R, Jura N, Cole PA, Kuriyan J: Inhibition of the EGF receptor by binding of MIG6 to an activating kinase domain interface. Nature. 2007, 450: 741-744. 10.1038/nature05998.
Nagashima T, Ushikoshi-Nakayama R, Suenaga A, Ide K, Yumoto N, Naruo Y, Takahashi K, Saeki Y, Taiji M, Tanaka H, Tsai SF, Hatakeyama M: Mutation of epidermal growth factor receptor is associated with MIG6 expression. FEBS J. 2009, 276: 5239-5251. 10.1111/j.1742-4658.2009.07220.x.
Wang Z, Raines LL, Hooy RM, Roberson H, Leahy DJ, Cole PA: Tyrosine phosphorylation of Mig6 reduces its inhibition of the epidermal growth factor receptor. ACS Chem Biol. 2013, 8: 2372-2376. 10.1021/cb4005707.
Jechlinger M, Sommer A, Moriggl R, Seither P, Kraut N, Capodiecci P, Donovan M, Cordon-Cardo C, Beug H, Grünert S: Autocrine PDGFR signaling promotes mammary cancer metastasis. J Clin Invest. 2006, 116: 1561-1570. 10.1172/JCI24652.
Lemmon MA, Schlessinger J: Cell signaling by receptor tyrosine kinases. Cell. 2010, 141: 1117-1134. 10.1016/j.cell.2010.06.011.
Cuevas BD, Lu Y, Mao M, Zhang J, LaPushin R, Siminovitch K, Mills GB: Tyrosine phosphorylation of p85 relieves its inhibitory activity on phosphatidylinositol 3-kinase. J Biol Chem. 2001, 276: 27455-27461. 10.1074/jbc.M100556200.
Cook RS, Garrett JT, Sánchez V, Stanford JC, Young C, Chakrabarty A, Rinehart C, Zhang Y, Wu Y, Greenberger L, Horak ID, Arteaga CL: ErbB3 ablation impairs PI3K/Akt-dependent mammary tumorigenesis. Cancer Res. 2011, 71: 3941-3951. 10.1158/0008-5472.CAN-10-3775.
Shah SP, Roth A, Goya R, Oloumi A, Ha G, Zhao Y, Turashvili G, Ding J, Tse K, Haffari G, Bashashati A, Prentice LM, Khattra J, Burleigh A, Yap D, Bernard V, McPherson A, Shumansky K, Crisan A, Giuliany R, Heravi-Moussavi A, Rosner J, Lai D, Birol I, Varhol R, Tam A, Dhalla N, Zeng T, Ma K, Chan SK, et al: The clonal and mutational evolution spectrum of primary triple-negative breast cancers. Nature. 2012, 486: 395-399.
Network CGA: Comprehensive molecular portraits of human breast tumours. Nature. 2012, 490: 61-70. 10.1038/nature11453.
Knight JF, Lesurf R, Zhao H, Pinnaduwage D, Davis RR, Saleh SM, Zuo D, Naujokas MA, Chughtai N, Herschkowitz JI, Prat A, Mulligan AM, Muller WJ, Cardiff RD, Gregg JP, Andrulis IL, Hallett MT, Park M: Met synergizes with p53 loss to induce mammary tumors that possess features of claudin-low breast cancer. Proc Natl Acad Sci U S A. 2013, 110: E1301-E1310. 10.1073/pnas.1210353110.
Hitosugi T, Kang S, Vander Heiden MG, Chung TW, Elf S, Lythgoe K, Dong S, Lonial S, Wang X, Chen GZ, Xie J, Gu TL, Polakiewicz RD, Roesel JL, Boggon TJ, Khuri FR, Gilliland DG, Cantley LC, Kaufman J, Chen J: Tyrosine phosphorylation inhibits PKM2 to promote the Warburg effect and tumor growth. Sci Signal. 2009, 2: ra73-
Hitosugi T, Zhou L, Fan J, Elf S, Zhang L, Xie J, Wang Y, Gu TL, Alečković M, LeRoy G, Kang Y, Kang HB, Seo JH, Shan C, Jin P, Gong W, Lonial S, Arellano ML, Khoury HJ, Chen GZ, Shin DM, Khuri FR, Boggon TJ, Kang S, He C, Chen J: Tyr26 phosphorylation of PGAM1 provides a metabolic advantage to tumours by stabilizing the active conformation. Nat Commun. 2013, 4: 1790-10.1038/ncomms2759.
Berger MF, Hodis E, Heffernan TP, Deribe YL, Lawrence MS, Protopopov A, Ivanova E, Watson IR, Nickerson E, Ghosh P, Zhang H, Zeid R, Ren X, Cibulskis K, Sivachenko AY, Wagle N, Sucker A, Sougnez C, Onofrio R, Ambrogio L, Auclair D, Fennell T, Carter SL, Drier Y, Stojanov P, Singer MA, Voet D, Jing R, Saksena G, Barretina J, et al: Melanoma genome sequencing reveals frequent PREX2 mutations. Nature. 2012, 485: 502-506.
This work was supported by the National Health and Medical Research Council of Australia, the National Breast Cancer Foundation and the Cancer Council New South Wales. AS is a Career Development Fellow of the National Health and Medical Research Council of Australia.
The authors declare that they have no competing interests.
RJD conceived and designed the study. NAA, FH, HC, RN, SY, LZ and RJL undertook the experiments. JW, MP, HCL and NA performed bioinformatics analysis on the generated data sets. AS, CJO and SJC assisted with experimental design and data interpretation. NAA, AS and RJD wrote the paper. All authors read and approved the paper before submission.
Electronic supplementary material
Additional file 1: Figure S1.: (A) The lineage of the different tumour samples characterised. (B) Clustering relationships in comparison to tumour of origin for the xenografts. The coloured lines indicate which tumour transplants share a common origin. The coloured ovals indicate instances where only one transplant tumour from a particular clonal line was characterised. (PDF 123 KB)
Additional file 2: Table S1.: 763 pY sites identified from the pY profiling of p53, Her2 and PyMT mammary mouse tumours. (XLSX 166 KB)
Additional file 3: Figure S2.: (A) Gene Ontology annotation of the phosphoproteins identified from the three mouse tumour models. (B) Proportion of different kinase types identified (based on number and subclassification). (C) Contribution of different non-receptor tyrosine kinases, based on pY peptide spectral intensity. (D) Contribution of different receptor tyrosine kinases based on pY peptide spectral intensity. (PDF 1 MB)
Additional file 4: Figure S3.: Unsupervised hierarchical clustering using all 763 identified pY sites. (Imputation: k-NN + column mean). (PDF 245 KB)
Additional file 6: Figure S4.: Unsupervised hierarchical clustering using the 381 pY sites detected in at least 75% of one tumour type. (Imputation: k-NN+ row min). (PDF 175 KB)
Additional file 7: Figure S5.: (A) and (B) Immunoblot quantitation of phosphorylation on specific phosphosites compared with relative tyrosine phosphorylation on the same sites determined by MS. Data are presented as in Figure 2, and colour key and histogram apply to (A) and (B). The heat maps indicate the intensity of phosphorylation on the following specific sites: ErbB3 (pY1325) and p85 (pY467). The lane corresponding to tumour 1186 has been removed due to sample degradation. Because phosphorylation of Y1325 was not detected in HER2 tumour 3051 by MS, the positive signal obtained by Western blotting may reflect cross-reactivity of the commercial antibody with a different phosphorylation site on erbB3 or another erbB receptor. (PDF 611 KB)
Additional file 8: Figure S6.: Unsupervised hierarchical clustering using the 707 pY sites identified in p53 tumours. (Imputation: k-NN + column mean). (PDF 149 KB)
Additional file 9: Table S3.: Proteins showing >10-fold increased phosphorylation in 1203/4 versus the other p53 tumours. (XLSX 103 KB)
Additional file 11: Table S5.: Pathway analysis of proteins showing >10-fold increased phosphorylation in 1203/4 versus other p53 tumours. (XLSX 23 KB)
Additional file 12: Figure S7.: (A) MET knockdown using siRNA in cell lines derived from p53-null mouse tumours. Western blots were incubated with the indicated antibodies. # in 1201 indicates nonspecific band. (B) Effect of MET knockdown on cell proliferation by MTS assay measuring relative absorbance. At least three independent sets of replicate experiments were performed with the same trend observed. Data from one representative experiment are shown. *P < 0.05 and **P < 0.005 by t-test (siRNA treated vs control). 'Control' refers to ON-TARGETplus nontargeting siRNA, and 'Pool' refers to SMARTpool consisting of four individual siRNAs targeting mouse MET. (PDF 947 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Ali, N.A., Wu, J., Hochgräfe, F. et al. Profiling the tyrosine phosphoproteome of different mouse mammary tumour models reveals distinct, model-specific signalling networks and conserved oncogenic pathways. Breast Cancer Res 16, 437 (2014). https://doi.org/10.1186/s13058-014-0437-3