Kinome and mRNA expression profiling of high-grade osteosarcoma cell lines implies Akt signaling as possible target for therapy
© Kuijjer et al.; licensee BioMed Central Ltd. 2014
Received: 24 April 2013
Accepted: 14 January 2014
Published: 21 January 2014
High-grade osteosarcoma is a primary malignant bone tumor mostly occurring in adolescents and young adults, with a second peak at middle age. Overall survival is approximately 60%, and has not significantly increased since the introduction of neoadjuvant chemotherapy in the 1970s. The genomic profile of high-grade osteosarcoma is complex and heterogeneous. Integration of different types of genome-wide data may be advantageous in extracting relevant information from the large number of aberrations detected in this tumor.
We analyzed genome-wide gene expression data of osteosarcoma cell lines and integrated these data with a kinome screen. Data were analyzed in statistical language R, using LIMMA for detection of differential expression/phosphorylation. We subsequently used Ingenuity Pathways Analysis to determine deregulated pathways in both data types.
Gene set enrichment indicated that pathways important in genomic stability are highly deregulated in these tumors, with many genes showing upregulation, which could be used as a prognostic marker, and with kinases phosphorylating peptides in these pathways. Akt and AMPK signaling were identified as active and inactive, respectively. As these pathways have an opposite role on mTORC1 signaling, we set out to inhibit Akt kinases with the allosteric Akt inhibitor MK-2206. This resulted in inhibition of proliferation of osteosarcoma cell lines U-2 OS and HOS, but not of 143B, which harbors a KRAS oncogenic transformation.
We identified both overexpression and hyperphosphorylation in pathways playing a role in genomic stability. Kinome profiling identified active Akt signaling, which could inhibit proliferation in 2/3 osteosarcoma cell lines. Inhibition of PI3K/Akt/mTORC1 signaling may be effective in osteosarcoma, but further studies are required to determine whether this pathway is active in a substantial subgroup of this heterogeneous tumor.
KeywordsOsteosarcoma Tumor cell lines Kinome profiling Gene expression profiling Genomic instability Bone tumor
High-grade osteosarcoma is the most prevalent primary malignant bone tumor. Most frequently, the long bones of adolescents and young adults are affected, with a yearly incidence of approximately 5 cases per million per year . Patients are generally treated with high doses of neoadjuvant chemotherapy to prevent the outgrowth of micrometastases. In 15-25% of all patients, however, metastatic disease is clinically detectable at diagnosis and despite the intensive treatment, 45% of all patients develop distant metastases, the leading cause of death of osteosarcoma patients [2, 3]. The introduction of neoadjuvant chemotherapy in the 1970s has increased survival from 10-20% to approximately 60%. However, survival has reached a plateau, and new treatments are urgently needed [4–6]. Osteosarcoma is an extremely genomically unstable tumor, with karyotypes harboring numerous numerical and structural changes [7, 8]. In addition, osteosarcoma genotypes show a considerable degree of heterogeneity, both intra- and intertumoral. Both the complex genotype and its heterogeneity render it difficult to determine which genomic alterations are important in osteosarcomagenesis, as not all alterations may lead to a difference in mRNA, protein levels, or enzyme activity in the tumor tissue. Integration of different data types is therefore of particular relevance for studying a heterogeneous tumor with a complex genomic profile such as osteosarcoma. Genomic and expression data of osteosarcoma tumor samples have been integrated by different groups, and many of the reported recurrent osteosarcoma driver genes play a role in cell cycle regulation and maintenance of genomic stability [9, 10]. Yet, even though recurrent driver genes may provide knowledge on what pathways are affected that help tumor cells survive, such driver genes may not always be accessible as targets for treatment. This especially holds for pathways involved in genetic stability, since the damage is already done.
Oncogenic kinases are often active in tumor cells, and a number of kinases can be pharmacologically inhibited. Therapies targeting oncogenic kinases have provided promising results in inhibiting proliferation of cancer cells, and some kinases have been targeted in preclinical and clinical studies in childhood sarcomas (as reviewed in Wachtel et al. ), e.g. IGF1R and mTOR [12, 13]. An unbiased approach to identify active kinases in cancer is to perform kinome-wide screens. Such screens have previously been effectively used in other types of sarcoma and have led to the detection of specific targets for treatment [14, 15]. As combining the analysis of different data types using systems biology approaches can give a more complete impression of the state of a tumor cell, we set out to integrate genome-wide gene expression data of osteosarcoma cell lines with kinome profiling data. Osteosarcoma cell lines are widely available and have been shown to be representative for the tumor of origin, both on a genome-wide as on a functional level, and are therefore a good model to study osteosarcoma preclinically [9, 16].
We previously have performed genome-wide expression analysis on a panel of 19 osteosarcoma cell lines . In the present study, we compared these expression profiles with the different putative progenitor cells of osteosarcoma – mesenchymal stem cells (MSCs) and osteoblasts – in order to define the common denominator pathways that are deregulated in osteosarcoma. We then integrated expression data with a serine/threonine (Ser/Thr) kinome screen, to determine whether pathways with enrichment of differentially expressed genes show enrichment in of hyperphosphorylation as well. In order to detect overactive kinases in osteosarcoma, which may be potential targets for treatment, we identified the most significant pathways by a single-way analysis of the kinome profiling data.
Osteosarcoma cell lines were previously characterized and described . Human bone-marrow-derived MSCs were obtained from two osteosarcoma patients, and were characterized and handled as described . For kinome profiling of osteosarcoma versus MSCs, cells were cultured in Dulbecco’s Modified Eagle Medium (DMEM; Invitrogen, Carlsbad, CA, USA), supplemented with 10% fetal bovine serum (Greiner Bio-one, Frickenhausen, Germany), in order to eliminate differences in kinase activity caused by culture conditions. For inhibition experiments and kinome profiling of inhibition experiments, osteosarcoma cell lines 143B, U-2 OS, and HOS were maintained in RPMI 1640 supplemented with 10% fetal calf serum (both from Invitrogen, Carlsbad, CA). The human pre-B acute lymphoblastic leukemia cell line NALM-6 cell line was kindly provided by Mw. N. Duinkerken (Department of Hematology, Leiden University Medical Center, the Netherlands), and was maintained in Iscove’s Modified Dulbecco’s Medium (IMDM) supplemented with GlutaMAX-1 (Life Technologies, Carlsbad, CA) and 10% fetal bovine serum (Greiner Bio-one, Frickenhausen, Germany). All cells were regularly tested for mycoplasm and were genotyped before and after experiments using the Powerplex 1.2 system (Promega, Leiden, the Netherlands), as described previously , and using CellID STR profiling (Promega, Leiden, the Netherlands). Latest genotyping results are added in Additional file 1). Cell lines corresponded to the entries in the ATCC (http://www.atcc.org) and DSMZ (http://www.dsmz.de) databases.
Kinome profiling was performed on osteosarcoma cell lines 143B and U-2 OS and on two MSCs – MSC001 and MSC006. Cells at 80% confluence were washed twice with Phosphate buffered Saline and lysed with M-PER Mammalian Extraction Buffer, supplemented with Halt Phosphatase Inhibitor Cocktail and EDTA free Halt Protease Inhibitor Cocktail (Pierce Biotechnology, Rockford, IL), according to the manufacture’s protocol. Cells were incubated on ice for at least 30 minutes before collecting the lysates and centrifuging these for 15 minutes at 4°C at > 10,000× g. Protein concentration was measured using a detergent-compatible Protein Assay (Bio-Rad Laboratories, Hercules, CA) according to the manufacturer’s protocol. Samples were snap-frozen and stored at −70°C.
MK-2206 was dissolved in DMSO at a concentration of 10 mM and stored at −20°C. For 143B, U-2 OS, and HOS, 2,000, 4,000, and 2,000 cells/well respectively, were plated in a 96-wells plate. NALM-6, a human pre-B acute lymphoblastic leukemia (ALL) cell line, was included as a positive control, as ALL cell lines have been shown to be highly sensitive to MK-2206 . This cell line grows in suspension and was plated at 50,000 cells/well. After 24 hrs, MK-2206 was added in triplicate in different concentrations – 0 nM, 0.5 nM, 1 nM, 5 nM, 10 nM, 50 nM, 100 nM, 500 nM, 1 μM, 5 μM, and 10 μM. For 143B and HOS, the effect of concentrations of 2, 3, 4, and 5 nM was assessed as well. Cells were grown in the presence of inhibitor for 120 hours. Cell proliferation was determined by incubating the cells with reagent WST-1 (Roche, Basel, Switzerland) for 2 hrs and subsequently measured using a Wallac 1420 VICTOR2 (Perkin Elmer, Waltham, MA). Data were analyzed in Graphpad Prism 5.01 (http://www.graphpad.com). Relative IC50s were calculated using results from the different concentrations up to the highest dose where toxicity was not yet present. The results shown are representative results from at least three independent experiments.
Genome-wide gene expression profiling
We analyzed our previously published data of osteosarcoma cell lines (n = 19), MSCs (n = 12), and osteoblasts (n = 3) (GEO superseries, accession number GSE42352) . Microarray data processing and quality control were performed in the statistical language R version 2.15  as described previously .
Kinome profiling was performed on 1 μg of cell lysate on the serine/threonine (Ser/Thr) Kinase PamChip® peptide microarrays (PamGene, ‘s-Hertogenbosch, the Netherlands) according to the manufacturer’s protocol, essentially as described in Hilhorst et al. . This peptide microarray comprises 142 peptide sequences derived from human phosphorylation sites. Peptide phosphorylation is detected in time with a mixture of fluorescently labeled anti-phosphoserine/threonine antibodies. We used at least three technical replicates for each MSC line, and four technical replicates for the osteosarcoma cell lines. Images were taken every 5 minutes, over the course of 60 minutes. Signal quantification on phosphorylated peptides was performed in BioNavigator software (PamGene International, ‘s Hertogenbosch, the Netherlands). Subsequently, data were normalized in R  using the vsn package . Median signals at 60 minutes of incubation with the cell lysates were analyzed in Bioconductor  package array QualityMetrics  to identify poor quality samples, which were removed from further analysis. Technical replicates of good quality were averaged. To determine whether these data were reproducible, we analyzed data from different cycles (0, 10, 20, 30, 40, 50, and 60 minutes incubation with cell lysates).
In the second kinome profiling experiment we compared lysates of untreated cells with lysates of cells treated with MK-2206. Different treatment durations and concentrations were used – no treatment, treatment for 5, 30, 180, and 960 minutes with 1 μM MK-2206, and treatment for 180 minutes with 10 μM of the drug. Kinome profiling was performed as described above, with the difference that we used 1–5 technical replicates per condition. Of this experiment, we analyzed signals at 30 minutes of incubation with the lysates.
Statistical analyses of microarray data
We performed LIMMA analysis  in order to determine differential mRNA expression between osteosarcoma cell lines (n = 19) and control cell lines – MSCs (n = 12) and osteoblasts (n = 3) and to determine differential phosphorylation of peptides on the PamChip® microarray between osteosarcoma cell lines (n = 2) and MSCs (n = 2). We used a Benjamini and Hochberg False Discovery Rate (FDR) of 0.05 as cut-off for significance. Kinome profiling signals obtained for the different treatment conditions were analyzed in a paired approach, in which signals from untreated cells were subtracted from the signals from treated cells. For both kinome profiling experiments, we used a cut-off of 0.1 for the absolute log fold change (logFC). Heatmaps were generated using the function heatmap.2 of R package gplots.
In order to reveal pathways which were significantly affected on mRNA levels in osteosarcoma cell lines, we intersected the toptables obtained by LIMMA analysis of osteosarcoma cell lines versus MSCs and of osteosarcoma cell lines versus osteoblasts. Gene symbols for all probes were imported into the software Ingenuity Pathways Analysis (IPA, Ingenuity Systems, http://www.ingenuity.com), together with FDR adjusted P-values (adjP) and average logFCs. Only the gene symbols of probes that were both significantly upregulated or both significantly downregulated in osteosarcoma cell lines as compared with MSCs and with OBs (adjP < 0.05) were selected to be considered as significantly differentially expressed in the IPA analysis. For differential phosphorylation, we imported the results from the LIMMA analysis on kinome profiling data, with a cut-off of 0.05 for adjusted P-value and a cut-off of 0.1 for logFC. The significance of the association between the data set and the canonical pathways was measured as described previously . Pathways with adjP < 0.05 were considered to be significantly affected. In addition, transcription factor analyses were performed on gene expression data in IPA in order to predict activated or inhibited transcription factors based on expression of target genes, returning p-values (with a cut-off of 0.05 for significance) and regulation z-scores.
Genome-wide gene expression profiling of high-grade osteosarcoma cell lines
Gene expression is altered in pathways regulating genomic stability
Kinome profiling of osteosarcoma cell lines
Altered phosphorylation in genomic stability pathways
PI3K/Akt and AMPK signaling in osteosarcoma
MK-2206 inhibits proliferation of U-2 OS and HOS, but not of 143B
Different phosphorylation patterns upon treatment with MK-2206
Osteosarcoma is a highly genomically unstable tumor. The identification of specific molecular targets that drive oncogenesis and that might be targets for therapy may thereby be hampered. Genome-wide gene expression profiling of high-grade osteosarcoma cell lines, in fact, showed an enrichment of differential expression in pathways important in genomic stability (Figure 2), with a role in cell cycle and checkpoint regulation (e.g. p53 signaling, G1/S and G2/M checkpoint regulation), DNA damage response (e.g. ATM signaling, role of BRCA1 in DNA damage response), and purine/pyrimidine metabolism. Most significantly differentially expressed genes in these pathways were upregulated, for example DNA-PK, BRCA1, and CDC25A. Some downregulated genes were detected as well, such as CDKN1A, which has an inhibitory role on cell cycle progression, and genes downstream of TP53 (e.g. THBS1 and SERPINE1, encoding TSP1 and PAI-1, respectively). Expression levels of genes in these pathways in osteosarcoma pre-treatment biopsies correlated with survival, as was previously reported on the same dataset  by using the CIN25 signature .
IPA transcription factor analysis showed that MYC was the most significantly activated (z-score of 6.294), and TP53 the most significantly inactivated (z-score of −7.660) transcription factor. Other highly predicted activated transcription factors were e.g. E2F1/2/3 (Additional file 6). These different transcription factors are known to be affected in osteosarcoma [7, 9, 30]. The role of these transcription factors in cell cycle progression further confirms the importance of these pathways in osteosarcoma. Important to note is that we took a different approach to determine significantly altered pathways from in our previous publications [9, 31]. We only used overlapping genes with same pattern of expression (both significant up- or downregulation) in osteosarcoma cell lines versus both control sets. This approach ensured us that all genes detected in the enriched pathways are significantly up- or downregulated in both comparisons, while our previous analyses described pathways which are significantly altered, but for which the gene list per pathway accounting for the significant effect may be different. We specifically took this more conservative approach for our current study, because we wanted to directly compare the expression levels and kinase activities of the specific players in each pathway. We also hypothesized that, using a method testing the overall aberration of a pathway, it would be more difficult to pick up specific players to inhibit pharmacologically. The pathways we detected with this analysis – pathways playing a role in cell cycling and genomic instability – were, as expected, also significantly affected in the less conservative globaltest analysis (which tests groups of genes instead of single genes) reported in our recent BMC Cancer publication  (data not shown).
Given the extreme genomic instability which is notorious in osteosarcoma and has led to the formulation of a novel genetic mechanism, chromothripsis , it is not surprising that the most prominent pathways are associated with this signature. Unfortunately pharmacological targeting of genomic instability is a challenge. Kinome-wide screens have previously led to the detection of specific targets for treatment in other sarcoma types [14, 15], and as such a screen can complement us with extra information on aberrations in the pathways we detected with gene expression analyses, we performed kinome profiling of osteosarcoma cell lysates. Since the pathways that were shown to be significantly affected on mRNA expression mostly contained Ser/Thr kinases, we selected a Ser/Thr peptide microarray – the Ser/Thr PamChip®. Pathway analysis on kinome profiling data showed that 50% of the pathways that were significant on gene expression data were also significantly enriched in differential phosphorylation signals (Figure 4). All significant peptides were higher phosphorylated in osteosarcoma cell lines, except for a peptide present in CREB1. Since most of these peptides showed higher phosphorylation, we expect these pathways to be highly active, demonstrating higher cell cycling of the tumor cells, and deregulated responses to DNA damage.
We next determined the most significantly affected pathways in the kinome data from the entire IPA canonical pathways database, and detected deregulation of the PI3K/Akt and AMPK signaling pathways. Molecules downstream of Akt kinases showed higher phosphorylation (Figure 5), while downstream of AMPK, lower levels of phosphorylation were detected (Additional file 9). Akt and AMPK act antagonistically to regulate mTOR signaling through inhibitory and activating phosphorylation of TSC2, respectively . The Akt pathway is one of the most commonly affected pathways in cancer, with active PI3K/Akt signaling leading to excessive cell growth and proliferation [34, 35]. Inhibition of this pathway by targeting mTOR with agents such as rapamycin is effective in some cancer types . In a recent phase II trial in bone and soft tissue sarcomas, inhibition of mTOR with ridaforolimus resulted in better progression-free survival . Inhibiting mTOR can, however, also activate a strong negative feedback loop from S6K1 to enhance Akt signaling [34, 36]. It may, therefore, be more effective to inhibit Akt itself. Inhibition of Akt was recently tested in a panel of xenografts of different pediatric cancers, and was most effective in osteosarcoma, with significant differences in event-free survival in 6/6 xenografts . In addition, AMPK activators suppress growth of cell lines of various tumor types .
We treated osteosarcoma cell lines with the allosteric Akt inhibitor MK-2206 (Selleck Chemicals LLC, Houston, TX). Inhibition of proliferation was dose-dependent in U-2 OS (IC50 2.5 μM), but not in 143B (Figure 6). Important to note is that active Akt signaling can be detected by kinome profiling in this cell line, but this does not necessarily imply that this pathway can also be fully inhibited, for example in the case that downstream actors in the same pathway cause a survival benefit for the cell line. As 143B is derived from the HOS cell line with a KRAS oncogenic transformation, we determined inhibitory effects of MK-2206 on HOS as well. HOS responded to MK-2206 in a similar manner as U-2 OS (IC50 2.6 μM). This suggests that constitutive Ras/Raf/ERK signaling causes insensitivity to inhibition of the Akt pathway to MK-2206. Kinome profiling of cells treated with MK-2206 resulted in different phosphorylation patterns in 143B and U-2 OS of peptides of molecules in the PI3K/Akt pathway (Figure 7). Differences between these cell lines were found in BAD Ser-99, of which phosphorylation was inhibited after treatment with MK-2206 in the responsive cell line U-2 OS, but stimulated in 143B, and in BAD Ser-118, where an opposite pattern was detected. BAD Ser-99 is the major site of Akt phosphorylation, while Ser-118 is the major site of PKA phosphorylation . Opposite patterns were also detected for TP53 Thr-18 and CDKN1A Thr-145/Ser-146, of which CDKN1A Thr-145 can also be directly phosphorylated by Akt. These results suggest that activity of other kinases may be affected by inhibition of Akt using MK-2206, or by MK-2206 itself. This depends on the cellular context, as we otherwise would not have expected to detect any differences in a paired analysis for the different conditions in each cell type.
An important finding of our studies is that the PI3K/Akt and AMPK signaling pathways were detected with kinome profiling, while mRNA expression profiling did not result in the identification of these pathways. This suggests that in osteosarcoma, these pathways are regulated by phosphorylation rather than by transcriptional activity. Copy number and mRNA expression levels of Akt family members and their upstream players did not provide us with a possible mechanism for elevated Akt activity, although PTEN showed lower, but not significantly lower, gene expression levels in both cell lines as compared with the two MSC controls (data not shown). Gene expression and protein synthesis imply a long time commitment of a cell to potential activation of its synthesized proteins. Phosphorylation, on the other hand, provides a very rapid way to mobilize extra catalytic power for a short time, and allows fine-tuning of the activation of a pathway to the needs of a cell. This difference in time scale emphasizes the importance of applying different platforms for the analysis of a complex tumor as high-grade osteosarcoma.
In summary, this study shows that genomic stability pathways are deregulated on both mRNA and kinome levels, with most significantly affected genes being upregulated and/or phosphorylated. Akt was detected as most probably overactive in osteosarcoma, as downstream peptides were hyperphosphorylated as compared with MSCs. Akt inhibitor MK-2206 could inhibit 2/3 osteosarcoma cell lines. Based on these results, we conclude that attenuating the PI3K/Akt/mTOR pathway may be effective in a subset of osteosarcomas.
Description of additional files
The following additional files are available with the online version of this paper. Additional file 1 (.xls) includes the latest genotyping results of cell lines 143B and U2OS. Additional file 2 (.pdf) is a figure depicting unsupervised clustering of gene expression data. Additional file 3 (.pdf) is a figure showing differentially expressed genes in osteosarcoma cell lines versus control cell cultures. Additional file 4 (.pdf) depicts unsupervised clustering of all genes present in the significantly affected pathways determined by IPA analysis. Additional file 5 (.pdf) depicts Kaplan-Meier analysis of the different clusters detected in Additional file 4. Additional file 6 (.xls) is a table including results from the transcription factor activity prediction analysis in IPA. Additional file 7 (.pdf) is a Venn diagram showing significantly differentially phosphorylated peptides over time. Additional file 8 (.pdf) shows unsupervised clustering of technical replicates used in the kinome profiling experiment. Additional file 9 (.pdf) illustrates significant differential phosphorylation in the AMPK signaling pathway. Additional file 10 (.pdf) depicts distances between kinome profiling data of treated and untreated osteosarcoma cells using unsupervised clustering.
Acute lymphoblastic leukemia
False discovery rate
FDR adjusted P-value
Half maximal inhibitory concentration
Ingenuity pathways analysis
Log fold change
Mesenchymal stem cell
The authors would like to thank N. Duinkerken for the NALM-6 cell line. This study was funded by EuroBoNeT (LSHC-CT-2006-018814), the Dutch Cancer Society (KWF, 2008–4060), the Netherlands Organization for Health Research and Development (92003–399).
- Raymond AK, Ayala AG, Knuutila S: Conventional osteosarcoma. World Health Classification of Tumours. Pathology and Genetics of Tumours of Soft Tissue and Bone. Edited by: Fletcher CDM, Unni KK, Mertens F. 2002, Lyon: IARC Press, 264-270.Google Scholar
- Bacci G, Longhi A, Versari M, Mercuri M, Briccoli A, Picci P: Prognostic factors for osteosarcoma of the extremity treated with neoadjuvant chemotherapy - 15-year experience in 789 patients treated at a single institution. Cancer. 2006, 106: 1154-1161. 10.1002/cncr.21724.View ArticlePubMedGoogle Scholar
- Buddingh EP, Anninga JK, Versteegh MI, Taminiau AH, Egeler RM, van Rijswijk CS, Hogendoorn PCW, Lankester AC, Gelderblom H: Prognostic factors in pulmonary metastasized high-grade osteosarcoma. Pediatr Blood Cancer. 2010, 54: 216-221.PubMedGoogle Scholar
- Allison DC, Carney SC, Ahlmann ER, Hendifar A, Chawla S, Fedenko A, Angeles C, Menendez LR: A meta-analysis of osteosarcoma outcomes in the modern medical era. Sarcoma. 2012, 2012: 704872.View ArticlePubMedPubMed CentralGoogle Scholar
- Anninga JK, Gelderblom H, Fiocco M, Kroep JR, Taminiau AHM, Hogendoorn PCW, Egeler RM: Chemotherapeutic adjuvant treatment for osteosarcoma: where do we stand?. Eur J Cancer. 2011, 47: 2431-2445. 10.1016/j.ejca.2011.05.030.View ArticlePubMedGoogle Scholar
- Hattinger CM, Pasello M, Ferrari S, Picci P, Serra M: Emerging drugs for high-grade osteosarcoma. Expert Opin Emerg Drugs. 2010, 15: 615-634. 10.1517/14728214.2010.505603.View ArticlePubMedGoogle Scholar
- Cleton-Jansen AM, Buerger H, Hogendoom PCW: Central high-grade osteosarcoma of bone: diagnostic and genetic considerations. Curr Diagn Pathol. 2005, 11: 390-399. 10.1016/j.cdip.2005.08.005.View ArticleGoogle Scholar
- Szuhai K, Cleton-Jansen AM, Hogendoorn PC, Bovee JV: Molecular pathology and its diagnostic use in bone tumors. Cancer Genet. 2012, 205: 193-204. 10.1016/j.cancergen.2012.04.001.View ArticlePubMedGoogle Scholar
- Kuijjer ML, Rydbeck H, Kresse SH, Buddingh EP, Lid AB, Roelofs H, Burger H, Myklebost O, Hogendoorn PCW, Meza-Zepeda LA, Cleton-Jansen AM: Identification of osteosarcoma driver genes by integrative analysis of copy number and gene expression data. Genes Chromosomes Cancer. 2012, 51: 696-706. 10.1002/gcc.21956.View ArticlePubMedGoogle Scholar
- Lockwood WW, Stack D, Morris T, Grehan D, O’Keane C, Stewart GL, Cumiskey J, Lam WL, Squire JA, Thomas DM, O’Sullivan MJ: Cyclin E1 is amplified and overexpressed in osteosarcoma. J Mol Diagn. 2011, 13: 289-296. 10.1016/j.jmoldx.2010.11.020.View ArticlePubMedPubMed CentralGoogle Scholar
- Wachtel M, Schafer BW: Targets for cancer therapy in childhood sarcomas. Cancer Treat Rev. 2010, 36: 318-327. 10.1016/j.ctrv.2010.02.007.View ArticlePubMedGoogle Scholar
- Kolb EA, Kamara D, Zhang W, Lin J, Hingorani P, Baker L, Houghton P, Gorlick R: R1507, a fully human monoclonal antibody targeting IGF-1R, is effective alone and in combination with rapamycin in inhibiting growth of osteosarcoma xenografts. Pediatr Blood Cancer. 2010, 55: 67-75.View ArticlePubMedGoogle Scholar
- Chawla SP, Staddon AP, Baker LH, Schuetze SM, Tolcher AW, D’Amato GZ, Blay JY, Mita MM, Sankhala KK, Berk L, Rivera VM, Clackson T, Loewy JW, Haluska FG, Demetri GD: Phase II study of the mammalian target of rapamycin inhibitor ridaforolimus in patients with advanced bone and soft tissue sarcomas. J Clin Oncol. 2012, 30: 78-84. 10.1200/JCO.2011.35.6329.View ArticlePubMedGoogle Scholar
- Schrage YM, Briaire-de Bruijn IH, de Miranda NF, van OJ, Taminiau AHM, van Wezel T, Hogendoorn PCW, Bovee JVMG: Kinome profiling of chondrosarcoma reveals SRC-pathway activity and dasatinib as option for treatment. Cancer Res. 2009, 69: 6216-6222. 10.1158/0008-5472.CAN-08-4801.View ArticlePubMedGoogle Scholar
- Willems SM, Schrage YM, Bruijn IH, Szuhai K, Hogendoorn PCW, Bovee JVMG: Kinome profiling of myxoid liposarcoma reveals NF-kappaB-pathway kinase activity and casein kinase II inhibition as a potential treatment option. Mol Cancer. 2010, 9: 257.View ArticlePubMedPubMed CentralGoogle Scholar
- Mohseny AB, Machado I, Cai Y, Schaefer KL, Serra M, Hogendoorn PCW, Llombart-Bosch A, Cleton-Jansen AM: Functional characterization of osteosarcoma cell lines provides representative models to study the human disease. Lab Invest. 2011, 91: 1195-1205. 10.1038/labinvest.2011.72.View ArticlePubMedGoogle Scholar
- Ottaviano L, Schaefer KL, Gajewski M, Huckenbeck W, Baldus S, Rogel U, Mackintosh C, de Alava E, Myklebost O, Kresse SH, Meza-Zepeda LA, Serra M, Cleton-Jansen AM, Hogendoorn PCW, Buerger H, Aigner T, Gabbert HE, Poremba C: Molecular characterization of commonly used cell lines for bone tumor research: a trans-European EuroBoNet effort. Genes Chromosomes Cancer. 2010, 49: 40-51. 10.1002/gcc.20717.View ArticlePubMedGoogle Scholar
- Cleton-Jansen AM, Anninga JK, Briaire-de Bruijn I, Romeo S, Oosting J, Egeler RM, Gelderblom H, Taminiau AHM, Hogendoorn PCW: Profiling of high-grade central osteosarcoma and its putative progenitor cells identifies tumourigenic pathways. Br J Cancer. 2009, 101: 2064-10.1038/sj.bjc.6605482.View ArticlePubMedPubMed CentralGoogle Scholar
- Gorlick R, Maris JM, Houghton PJ, Lock R, Carol H, Kurmasheva RT, Kolb EA, Keir ST, Reynolds CP, Kang MH, Billups CA, Smith MA: Testing of the Akt/PKB inhibitor MK-2206 by the pediatric preclinical testing program. Pediatr Blood Cancer. 2011, 59: 518-524.View ArticlePubMedPubMed CentralGoogle Scholar
- R Development Core Team: R: A language and environment for statistical computing, reference index version 2.15.0. 2011, Vienna, Austria: R Foundation for Statistical ComputingGoogle Scholar
- Buddingh EP, Kuijjer ML, Duim RA, Burger H, Agelopoulos K, Myklebost O, Serra M, Mertens F, Hogendoorn PCW, Lankester AC, Cleton-Jansen AM: Tumor-infiltrating macrophages are associated with metastasis suppression in high-grade osteosarcoma: a rationale for treatment with macrophage-activating agents. Clin Cancer Res. 2011, 17: 2110-2119. 10.1158/1078-0432.CCR-10-2047.View ArticlePubMedGoogle Scholar
- Hilhorst R, Houkes L, Mommersteeg M, Musch J, van den Berg A, Ruijtenbeek R: Peptide microarrays for profiling of serine/threonine kinase activity of recombinant kinases and lysates of cells and tissue samples. Gene Regulation: Methods and Protocols, Methods in Molecular Biology. Edited by: Bina M. 2013, New York: Springer Science+Business Media, LLCGoogle Scholar
- Smyth GK: Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004, 3 (3).
- Huber W, Von HA, Sultmann H, Poustka A, Vingron M: Variance stabilization applied to microarray data calibration and to the quantification of differential expression. Bioinformatics. 2002, 18 (1): S96-S104.View ArticlePubMedGoogle Scholar
- Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, Hornik K, Hothorn T, Huber W, Iacus S, Irizarry R, Leisch F, Li C, Maechler M, Rossini AJ, Sawitzki G, Smith C, Smyth G, Tierney L, Yang JY, Zhang J: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004, 5: R80-10.1186/gb-2004-5-10-r80.View ArticlePubMedPubMed CentralGoogle Scholar
- Kauffmann A, Gentleman R, Huber W: arrayQualityMetrics–a bioconductor package for quality assessment of microarray data. Bioinformatics. 2009, 25: 415-416. 10.1093/bioinformatics/btn647.View ArticlePubMedGoogle Scholar
- Mohseny AB, Cai Y, Kuijjer M, Xiao W, van den Akker B, de Andrea CE, Jacobs R, ten Dijke P, Hogendoorn PCW, Cleton-Jansen AM: The activities of Smad and Gli mediated signalling pathways in high-grade conventional osteosarcoma. Eur J Cancer. 2012, 48: 3429-3438. 10.1016/j.ejca.2012.06.018.View ArticlePubMedGoogle Scholar
- Rhim JS, Cho HY, Huebner RJ: Non-producer human cells induced by murine sarcoma virus. Int J Cancer. 1975, 15: 23-29. 10.1002/ijc.2910150104.View ArticlePubMedGoogle Scholar
- Carter SL, Eklund AC, Kohane IS, Harris LN, Szallasi Z: A signature of chromosomal instability inferred from gene expression profiles predicts clinical outcome in multiple human cancers. Nat Genet. 2006, 38: 1043-1048. 10.1038/ng1861.View ArticlePubMedGoogle Scholar
- Mohseny AB, Tieken C, van der Velden PA, Szuhai K, de AC, Hogendoorn PCW, Cleton-Jansen AM: Small deletions but not methylation underlie CDKN2A/p16 loss of expression in conventional osteosarcoma. Genes Chromosomes Cancer. 2010, 49: 1095-1103. 10.1002/gcc.20817.View ArticlePubMedGoogle Scholar
- Kuijjer ML, Peterse EF, van den Akker BE, Briaire-de Bruijn IH, Serra M, Meza-Zepeda LA, Myklebost O, Hassan AB, Hogendoorn PC, Cleton-Jansen AM: IR/IGF1R signaling as potential target for treatment of high-grade osteosarcoma. BMC Cancer. 2013, 13: 245-10.1186/1471-2407-13-245.View ArticlePubMedPubMed CentralGoogle Scholar
- Stephens PJ, Greenman CD, Fu B, Yang F, Bignell GR, Mudie LJ, Pleasance ED, Lau KW, Beare D, Stebbings LA, McLaren S, Lin ML, McBride DJ, Varela I, Nik-Zainal S, Leroy C, Jia M, Menzies A, Butler AP, Teague JW, Quail MA, Burton J, Swerdlow H, Carter NP, Morsberger LA, Iacobuzio-Donahue C, Follows GA, Green AR, Flanagan AM, Stratton MR, Futreal PA, Campbell PJ: Massive genomic rearrangement acquired in a single catastrophic event during cancer development. Cell. 2011, 144: 27-40. 10.1016/j.cell.2010.11.055.View ArticlePubMedPubMed CentralGoogle Scholar
- Memmott RM, Dennis PA: Akt-dependent and -independent mechanisms of mTOR regulation in cancer. Cell Signal. 2009, 21: 656-664. 10.1016/j.cellsig.2009.01.004.View ArticlePubMedPubMed CentralGoogle Scholar
- Engelman JA, Luo J, Cantley LC: The evolution of phosphatidylinositol 3-kinases as regulators of growth and metabolism. Nat Rev Genet. 2006, 7: 606-619. 10.1038/nrg1879.View ArticlePubMedGoogle Scholar
- Manning BD, Cantley LC: AKT/PKB signaling: navigating downstream. Cell. 2007, 129: 1261-1274. 10.1016/j.cell.2007.06.009.View ArticlePubMedPubMed CentralGoogle Scholar
- Guertin DA, Sabatini DM: Defining the role of mTOR in cancer. Cancer Cell. 2007, 12: 9-22. 10.1016/j.ccr.2007.05.008.View ArticlePubMedGoogle Scholar
- Vakana E, Altman JK, Platanias LC: Targeting AMPK in the treatment of malignancies. J Cell Biochem. 2012, 113: 404-409. 10.1002/jcb.23369.View ArticlePubMedGoogle Scholar
- Hornbeck PV, Chabra I, Kornhauser JM, Skrzypek E, Zhang B: PhosphoSite: a bioinformatics resource dedicated to physiological protein phosphorylation. Proteomics. 2004, 4: 1551-1561. 10.1002/pmic.200300772.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1755-8794/7/4/prepub
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.