Skip to main content

Expression analysis of genes associated with human osteosarcoma tumors shows correlation of RUNX2 overexpression with poor response to chemotherapy

Abstract

Background

Human osteosarcoma is the most common pediatric bone tumor. There is limited understanding of the molecular mechanisms underlying osteosarcoma oncogenesis, and a lack of good diagnostic as well as prognostic clinical markers for this disease. Recent discoveries have highlighted a potential role of a number of genes including: RECQL4, DOCK5, SPP1, RUNX2, RB1, CDKN1A, P53, IBSP, LSAMP, MYC, TNFRSF1B, BMP2, HISTH2BE, FOS, CCNB1, and CDC5L.

Methods

Our objective was to assess relative expression levels of these 16 genes as potential biomarkers of osteosarcoma oncogenesis and chemotherapy response in human tumors. We performed quantitative expression analysis in a panel of 22 human osteosarcoma tumors with differential response to chemotherapy, and 5 normal human osteoblasts.

Results

RECQL4, SPP1, RUNX2, and IBSP were significantly overexpressed, and DOCK5, CDKN1A, RB1, P53, and LSAMP showed significant loss of expression relative to normal osteoblasts. In addition to being overexpressed in osteosarcoma tumor samples relative to normal osteoblasts, RUNX2 was the only gene of the 16 to show significant overexpression in tumors that had a poor response to chemotherapy relative to good responders.

Conclusion

These data underscore the loss of tumor suppressive pathways and activation of specific oncogenic mechanisms associated with osteosarcoma oncogenesis, while drawing attention to the role of RUNX2 expression as a potential biomarker of chemotherapy failure in osteosarcoma.

Peer Review reports

Background

Osteosarcoma is the most common pediatric tumor of the bone. Clinically, osteosarcoma has a bimodal distribution, with the majority of patients developing the disease during the period of active bone growth in early adolescence. The treatment generally involves surgery, often involving a loss of limb, and adjuvant chemotherapy. The best prognostic marker for osteosarcoma is the response to chemotherapy, where good response to chemotherapy is associated with an overall more favorable patient outcome and survival [1, 2].

At the molecular level osteosarcoma is characterized by a high level of genomic instability, highly heterogenous karyotypes both intra- and inter-tumor, and gross changes in gene expression [3–9]. Human osteosarcoma tumors often have osteoblast-like features but may vary within a broad range of epithelial mesenchymal lineages reflective of their poorly differentiated phenotype [10–12]. Therefore, assessment of the molecular changes in osteosarcoma tumors relative to normal osteoblasts can provide important insights concerning gene expression changes associated with both osteosarcoma oncogenesis and with molecular alterations governing differential clinical response to treatment.

The genetic change most commonly associated with osteosarcoma is the loss of the TP53 tumor suppressor gene through either genetic mutation or loss of gene expression [13–15]. Patients with Li-Fraumeni syndrome, which results from loss of TP53, have a strong predisposition to developing osteosarcoma [16]. Another tumor suppressor gene whose loss of expression is linked to osteosarcoma is RB1 [17]. In our recent studies, we have also shown that these genes play a central role in osteosarcoma-related gene expression networks both in human osteosarcoma cell lines [18] and tumor tissues [19]. We used a unique bioinformatic integrative whole-genome approach to map the genetic and epigenetic changes in osteosarcoma tumors and to identify gene networks related to osteosarcoma oncogenesis. Changes that showed the most significant associations with osteosarcoma gene networks included: overexpression and the most significant copy number gain of the chromosome 6p21.1 RUNX2 locus, loss of expression and genomic loss of the DOCK5 and TNFRSF10A loci at chromosome 8p21.1-21.3, and hypomethylation, copy number gain, and overexpression of the HISTH2BE gene at chromosome 1q21. Other genes that showed deregulated expression and significant contribution to osteosarcoma gene networks included overexpressed SPP1, IBSP, BMP2, and c-MYC, and uderexpressed CDKN1A, LSAMP, and CCNB1. Another gene that is thought to play a role in osteosarcoma and has been shown to be overexpressed is the FOS proto-oncogene [20, 21]. Similarly, CDC5L has been recently proposed as the putative oncogene at the 6p21 locus in osteosarcoma [22]. Finally, a DNA repair gene, RECQL4, has been shown to be overexpressed, and its level of overexpression correlates with overall genomic instability in osteosarcoma [23].

Microarray analysis and reverse-transcriptase polymerase chain reaction (RT-PCR) are useful for the molecular classification of tumors and for deriving biological mechanisms that underpin differential prognosis for patients with various types of cancer, including osteosarcoma [24–27]. The use of gene-expression profiling in clinical practice is however limited by the large number of genes that need to be analyzed and by the lack of reproducibility of various array platforms and interpretative methods [28]. Quantitative RT-PCR methods can be readily applied to RNA derived from formalin-fixed, paraffin-embedded (FFPE) pathological specimens, are reproducible and may be highly applicable in clinical practice [29], particularly for a rare tumor such as osteosarcoma in which access to frozen tissue is often limited. RT-PCR can only typically be used to analyze a small number of genes. Therefore it is important to select gene subsets for detailed analyses in which multiple lines of evidence implicate clinical utility. In previous studies, we [18, 19], and others [22, 30, 31] have performed microarray analyses of osteosarcoma cell lines and tissue samples and identified a series of genes with strong potential as biomarkers with clinical utility. Thus, the objective of the current study was to examine expression profiles of RECQL4, DOCK5, SPP1, RUNX2, RB1, CDKN1A, TP53, IBSP, LSAMP, MYC, TNFRSF1B, BMP2, HISTH2BE, FOS, CCNB1, and CDC5L genes in a cohort of osteosarcoma tumors and normal human osteoblasts. As a result, we show that RECQL4, SPP1, RUNX2, and IBSP are significantly overexpressed, and DOCK5, CDKN1A, RB1, P53, and LSAMP show significant loss of expression relative to normal osteoblasts. We also show that RUNX2 was the only gene with significant overexpression in tumors with an unfavorable response to chemotherapy relative to favorable responders.

Methods

Tissue samples

The collection of frozen tissue specimens (n = 15), archival formalin-fixed, paraffin-embedded osteosarcoma sections (n = 7), and clinicopathological data was obtained and handled in accordance with the Hospital for Sick Children Research Ethics guidelines (Toronto, Canada). This was a retrospective study of chemotherapy-naive biopsy samples collected sequentially between 1996 and 2005, and all specimens presented a tumor content higher than 90%. All patients were subjected during treatment to standard regimens for osteosarcoma, comprising cisplatin, doxorubicin, and methotrexate. The patient tumor specimens were revised at the time of study by the pathologist (P.T.). The details of the cases are presented in Table 1. The Huvos grading system was used to rate the level of tumor necrosis following preoperative chemotherapy: Grade I, little or no effect of chemotherapy noted; Grade II, partial response to chemotherapy, with between 50% and 90% necrosis; Grade III, greater than 90% necrosis; and Grade IV, no viable tumor cells are apparent [32]. The good responders are patients with necrosis ≥ 90% [33]. Normal human osteoblasts that were isolated from surgical bone specimens from five healthy individuals were obtained from Promocell (Heidelberg, Germany).

Table 1 Descriptive and histopathological features of the tumor cohort and normal osteoblasts.

RNA isolation

Total RNA from snap-frozen tissue (5 normal human osteoblasts and 15 tumors) was extracted and purified using the TRIzol Reagent method according to the manufactures protocol (Invitrogen, Carlsbad, CA, USA). FFPE tissues (7 tumors) were deparaffinized with xylene, washed with ETOH, and digested with a proteinase K buffer [34]. Total RNA was extracted and purified as above with TRIzol reagent. The RNA quality was good for all samples as assessed by BioAnalyzer RNA 600 Nano Kit (Agilent Technologies, Palo Alto, CA). Additionally, the overall trends for expression in tumors relative to normal osteoblasts were similar in frozen and FFPE samples corroborating the quality of the extracted RNA (Additional file 1).

Quantification of mRNA Expression

Quantitative real-time PCR (qRT-PCR) was used to quantify mRNA expression levels of 16 genes (RUNX2, DOCK5, TNFRS1B, HISTH2BE, P21, SSP1, P53, IBSP, CCNB1, BMP2, LSAMP, RB1, FOS, MYC, RECQL4, and CDC5L). Briefly, 1-2 ug of total RNA was converted to cDNA using GeneAmp Gold RNA PCR Core Kit (Applied Biosystems, Foster City, CA, USA), as per manufacturers recommendations. Primers were designed to specifically amplify templates of approximately 90-130 nucleotides overlapping exon boundaries of 3' terminal exons using the Primer-Blast software http://www.ncbi.nlm.nih.gov/tools/primer-blast/. The primers were subsequently tested by both PCR and qRT-PCR for specificity and single band amplification. The sequences of the PCR primer pairs used for each gene are shown in Additional file 2. The qRT-PCR assays for a particular gene were undertaken at the same time for all samples under identical conditions, in duplicate. The cycling conditions were as follows: 95°C for 2 min, 40 cycles of 95°C for 15 sec and 60°C for 45 sec, with a final extension 72°C for 5 min.

The mRNA expression levels were determined using Platinum SYBR Green qPCR Supermix-UDG with Rox (Invitrogen, Carlsbad, CA, USA), and the Applied Biosystems Prism 7900 Sequence Detection System (PE Applied Biosystems, Inc., Foster City, CA). The relative expression level of the genes of interest was computed relative to the endogenous control, phosphoglycerate kinase (PGK), to normalize for variances in the quality of RNA and the amount of input cDNA. Additionally, we validated our experimental conditions by analyzing the expression of two genes identified in our previous microarray study [19]; A2M, the highest overexpressing gene of the set and our positive control for the current study; and SLC14A, the gene with the lowest expression of the set and our negative control for the current study (Additional file 3). The mRNA expression levels for each sample were determined as fold-change values relative to the mean baseline expression levels for five human osteoblasts (HOBs), using the delta delta Ct method of analysis [35].

Statistical analysis

Results of the delta-delta Ct analysis were log10 transformed and imported to Partek Genomic Suite software (Additional file 1). The tumor samples were grouped by Huvos grade into those with favorable response to chemotherapy (Grades I and II) and unfavorable response (Grades III and IV); and were compared either as a group to normal human osteoblasts (i.e. tumor vs. normal), or to each other (i.e. unfavorable vs. favorable). Differences in p-values between groups were obtained using the non-parametric rank-sum Mann-Whitney test, and fold change differences between groups were obtained using the 1-way ANOVA tool (Additional file 4) using the Partek Genomic Suite software. This study was designed as a confirmatory analysis based on specific genes that were previously shown to have significant expression changes in steosarcoma, and thus multiple test correction was not applicable [36]. This analysis is designed to assess the correlation of specific gene expression as an individual parameter against the ostosarcoma phenotype, and is not meant to assess these genes as a group, gene network, or a multiple gene signature.

Results

The pathology of 22 osteogenic sarcomas and follow-up biopsies were analysed to determine response to chemotherapy and percent of necrotic tissue. This allowed tumors to be grouped as good responders (Huvos grades III and IV with favorable responses to chemotherapy) and poor responders (Huvos grades I and II with unfavorable responses to chemotherapy) (Table 1). Eleven tumors were identified to have 95% or more necrosis as a result of chemotherapy and were labelled as favorable responders, and the remaining 11 tumors were characterized as unfavorable responders. The baseline control consisted of a panel of five normal human osteoblast samples. The majority of tumors displayed osteoblastic histology, and most tumors and all osteoblast controls were of femoral origin.

In order to quantitatively assess the expression of the target genes (RECQL4, DOCK5, SPP1, RUNX2, RB1, CDKN1A, TP53, IBSP, LSAMP, MYC, TNFRSF1B, BMP2, HISTH2BE, FOS, CCNB1, and CDC5L) we performed qRT-PCR on the tumor cohort and human osteoblast samples (Additional file 1). Statistical analysis of these data revealed significant changes in a number of genes (Table 2). Tumors displayed significant overexpression of RECQL4, SPP1, RUNX2, and IBSP genes and loss of expression of DOCK5, CDKN1A, RB1, TP53, and LSAMP (p < 0.05) (Figure 1). The highest level of overexpression was measured in SPP1 with 113-fold overexpression, while the largest reduction of expression of 36-fold was evident in the DOCK5 gene. Comparison of tumors with unfavorable response to ch emotherapy to favorable responders revealed RUNX2 as the only significant gene (p = 0.03). On average unfavorable responders to chemotherapy showed 3.3-fold increase in the RUNX2 gene expression relative to favorable responders. Furthermore, RUNX2 expression showed a trend towards overexpression going from normal osteoblasts to favorable responders to chemotherapy and then to unfavorable responders to chemotherapy (Figure 1). The tumor sample #256 that exhibited the worst response to chemotherapy, also showed highest levels (113-fold) of RUNX2 overexpression (Additional file 1).

Table 2 Statistical evaluation of the tumor-specific and chemotherapy response-related gene expression signatures.
Figure 1
figure 1

Expression analysis of osteosarcoma-related genes. qRT-PCR levels of gene expression of 16 osteosarcoma-related genes in 5 normal human osteoblasts and 22 human osteosarcoma samples are shown on y-axis. The samples are grouped based on the response to chemotherapy status on y-axis. Corresponding box and whiskers plots representing the mean, 25th and 75th percentile (boxes), and 10th and 90th percentile (whiskers) are also shown. From left to right, and top to bottom of the panel, the plots are placed in the order of Mann-Whitney p-value significance (tumor vs. normal). HOB: normal human osteoblasts; Good: Favorable response to chemotherapy (Huvos grades III and IV); Poor: Unfavorable response to chemotherapy (Huvos grades I and II).

The remaining genes showed no significant changes in expression in tumors relative to normal osteoblasts. MYC, BMP2, and FOS show an associative trend when overexpressed, and TNFRSF10A shows a trend in underexpression in the tumor cohort. Notably, CDC5L showed essentially no change in expression levels relative to normal osteoblasts.

Discussion

Transformation of normal cells and initiation of tumorigenesis involves a combination of genetic and epigenetic changes [37]. Progressive acquisition of such changes ultimately results in destabilization of the genome, deregulation of gene expression pathways and activation of oncogenic gene expression networks. Identification of key genes, or network "nodes", will provide a more comprehensive understanding of tumorigenic processes, and provide more effective diagnostic, prognostic and therapeutic markers. In our recent integrative epi/genomic studies of osteosarcoma cell lines [18] and tumor genomes [19] we identified a number of such genes. A survey of the current literature allowed us to augment this list to a total of sixteen genes whose expression levels were assessed in an expanded tumor cohort. By comparing the gene expression levels to a panel of normal human osteoblasts, it allowed us to identify changes that are likely to be involved in osteosarcoma tumorigenesis. As a result, we identified significant disruptions of gene expression in nine of these genes, including loss of expression of DOCK5, CDKN1A, RB1, TP53, and LSAMP genes, and overexpression of RECQL4, SPP1, RUNX2, and IBSP genes.

Our study demonstrates a significant deregulation of proteins in osteosarcoma that are important effectors in the cell cycle and in differentiation. Significantly, we detected loss of TP53 expression, which may also play a role in loss of expression of CDKN1A, which encodes the cyclin-dependent kinase inhibitor 1A (p21) and whose expression is activated by p53 [38]. Loss of DOCK5 may play a similar role as it was shown recently that DOCK5 expression is essential for bone differentiation, from precursor osteoclasts [39]. Interestingly, in our recent study we showed that DOCK5 is located in the most significant region of copy number loss in osteosarcoma 8p21.2-p21.3 [19] along with the TNFRSF10A gene, for which we see an overall trend of loss of expression in our current study. TNFRSF10A is a receptor activated by tumor necrosis factor-related apoptosis inducing ligand TNFSF10 (also known as TRAIL), and is involved in the transduction of cell death signal and induction of cell apoptosis, which is mostly independent of p53 signalling [40].

A significant region of copy number loss at 3q13.31 has been identified in our previous study [19], and has also been observed by another group in 56% of osteosarcomas. This region of copy number loss was shown to correlate with loss of expression and hypermethylation of the LSAMP gene, and the authors proposed LSAMP as a novel tumor suppressor in osteosarcoma [31]. Our data agree with these findings and show significant loss of expression of LSAMP in the majority of our osteosarcoma samples. LSAMP codes for a neuronal surface glycoprotein found in cortical and subcortical regions of the limbic system, but it is currently unclear how this gene may be related to osteosarcoma tumorigenesis.

RUNX2 was one of the genes overexpressed in our set of tissue samples, and the only gene whose overexpression was significantly related to poor response to chemotherapy in osteosarcomas. RUNX2 is a member of the Runx transcription factor family consisting of RUNX1, RUNX2, and RUNX3 which function in the development of a number of tissues [41]. Of the three proteins, both RUNX1 and RUNX2 have been associated with oncogenesis. RUNX1 and RUNX2 upregulate LGALS3 (galectin-3) [42], a protein which suppresses anoikis and drug-induced apoptosis [43] and whose expression is correlated with metastasis in osteosarcoma [44] and progression in glioma [45]. Similarly, in lymphoma, overexpression of RUNX2 and MYC results in the "collaboration" of the two corresponding proteins to attenuate apoptosis and promote proliferation [46]. In developing osteoblasts, expression of RUNX2 normally decreases during maturation [47], and overexpression of the gene leads to a higher rate of bone turnover [48]. In bone metastatic breast cancer, RUNX2 promotes cancer cell survival and growth by activating expression of IHH and interacts with the TGFβ/BMP signal transduction pathway to parathyroid hormone-related protein (PTHrP) [49]. Thus the finding of elevated expression of RUNX2 in osteosarcomas with an unfavorable chemotherapy response is consistent with its oncogenic potential noted in other studies

In addition to RUNX2, three other genes, SPP1, RECQL4, and IBSP showed significant overexpression. SPP1 (osteopontin), like RUNX2, is a member of the BMP-signalling protein family. It shows the highest overexpression (113-fold) in our analysis, and has previously been shown to be significantly overexpressed in osteosarcoma tumors [50] and cell lines [51]. Antisense knockdown of SPP1 RNA in osteosarcoma cells results in inhibition of in vivo tumorigenesis in mice. These findings suggest that overexpression of SPP1 plays a role in osteosarcoma tumorigenesis, in particular, in cells lacking expression of cell cycle regulators and differentiation-related genes, as discussed earlier. RECQL4 is a gene whose protein product is involved in repair of DNA double stranded breaks and deregulation of its expression was recently shown to be strongly correlated with genomic instability in osteosarcoma [23]. Our data are consistent with this in further reaffirming the association of RECQL4 overexpression with osteosarcoma tumorigenesis. The final gene which showed significant overexpression in tumors relative to normal cells, IBSP, also known as bone sialoprotein, is a marker of terminal differentiation of bone. In normal osteoblasts RUNX2 and HDAC3 have been shown to suppress IBSP, and upon terminal differentiation loss, of RUNX2 expression derepresses IBSP and allows for terminal differentiation [52]. In tumors, however, although IBSP shows significant overexpression relative to osteoblasts, concurrent overexpression of RUNX2 indicates possible disruption of the terminal differentiation process.

By comparing our tumor panel, the majority of which (17/22) are the osteoblastic histological subtype, we identified gene expression changes associated with osteosarcoma oncogenesis. However, a possible limitation of this study is that some of the genetic associations described may not be present in osteosarcoma lineages that arise from more primitive cells of origin, including mesenchymal precursors [53–55]. In addition to the associations of gene expression with oncogenesis, we showed that the RUNX2 gene displays significant increase in expression in tumors with poor response to chemotherapy relative to the good responders. Our results may be reflective of either increased levels of gene expression in individual tumor cells during disease progression, or alternatively an increased proportion of cell lineages with RUNX2 expression in these genetically highly heterogenous cells. Both scenarios would allow selective advantage to the evolving cell lineages during tumorigenesis. Detailed immunohistochemical and imaging experiments will be required to further delineate these possibilities. It also remains to be determined if this correlation is evident at the protein level, and if so, RUNX2 may be a good histological marker for chemotherapy response, which is currently the best predictor of overall outcome in patients with osteosarcoma.

Conclusion

These data underscore the loss of tumor suppressive pathways the deregulation of cell cycle control proteins, and the activation of specific oncogenic mechanisms associated with osteosarcoma oncogenesis. Our results also draw attention to the role of RUNX2 expression as a potential biomarker of chemotherapy failure in osteosarcoma.

References

  1. Sakamoto A, Iwamoto Y: Current status and perspectives regarding the treatment of osteo-sarcoma: chemotherapy. Rev Recent Clin Trials. 2008, 3 (3): 228-231. 10.2174/157488708785700267.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Gorlick R: Osteosarcoma: clinical practice and the expanding role of biology. Journal of musculoskeletal & neuronal interactions. 2002, 2 (6): 549-551.

    CAS  Google Scholar 

  3. Batanian JR, Cavalli LR, Aldosari NM, Ma E, Sotelo-Avila C, Ramos MB, Rone JD, Thorpe CM, Haddad BR: Evaluation of paediatric osteosarcomas by classic cytogenetic and CGH analyses. Mol Pathol. 2002, 55 (6): 389-393. 10.1136/mp.55.6.389.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Lim G, Karaskova J, Beheshti B, Vukovic B, Bayani J, Selvarajah S, Watson SK, Lam WL, Zielenska M, Squire JA: An integrated mBAND and submegabase resolution tiling set (SMRT) CGH array analysis of focal amplification, microdeletions, and ladder structures consistent with breakage-fusion-bridge cycle events in osteosarcoma. Genes Chromosomes Cancer. 2005, 42 (4): 392-403. 10.1002/gcc.20157.

    Article  CAS  PubMed  Google Scholar 

  5. Lim G, Karaskova J, Vukovic B, Bayani J, Beheshti B, Bernardini M, Squire JA, Zielenska M: Combined spectral karyotyping, multicolor banding, and microarray comparative genomic hybridization analysis provides a detailed characterization of complex structural chromosomal rearrangements associated with gene amplification in the osteosarcoma cell line MG-63. Cancer Genet Cytogenet. 2004, 153 (2): 158-164. 10.1016/j.cancergencyto.2004.01.016.

    Article  CAS  PubMed  Google Scholar 

  6. Sandberg AA, Bridge JA: Updates on the cytogenetics and molecular genetics of bone and soft tissue tumors: osteosarcoma and related tumors. Cancer Genet Cytogenet. 2003, 145 (1): 1-30. 10.1016/S0165-4608(03)00105-5.

    Article  CAS  PubMed  Google Scholar 

  7. Selvarajah S, Yoshimoto M, Maire G, Paderova J, Bayani J, Squire JA, Zielenska M: Identification of cryptic microaberrations in osteosarcoma by high-definition oligonucleotide array comparative genomic hybridization. Cancer Genet Cytogenet. 2007, 179 (1): 52-61. 10.1016/j.cancergencyto.2007.08.003.

    Article  CAS  PubMed  Google Scholar 

  8. Squire JA, Pei J, Marrano P, Beheshti B, Bayani J, Lim G, Moldovan L, Zielenska M: High-resolution mapping of amplifications and deletions in pediatric osteosarcoma by use of CGH analysis of cDNA microarrays. Genes Chromosomes Cancer. 2003, 38 (3): 215-225. 10.1002/gcc.10273.

    Article  CAS  PubMed  Google Scholar 

  9. Zielenska M, Marrano P, Thorner P, Pei J, Beheshti B, Ho M, Bayani J, Liu Y, Sun BC, Squire JA, et al: High-resolution cDNA microarray CGH mapping of genomic imbalances in osteosarcoma using formalin-fixed paraffin-embedded tissue. Cytogenet Genome Res. 2004, 107 (1-2): 77-82. 10.1159/000079574.

    Article  CAS  PubMed  Google Scholar 

  10. Tang N, Song WX, Luo J, Haydon RC, He TC: Osteosarcoma development and stem cell differentiation. Clinical orthopaedics and related research. 2008, 466 (9): 2114-2130. 10.1007/s11999-008-0335-z.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Thomas D, Kansara M: Epigenetic modifications in osteogenic differentiation and transformation. J Cell Biochem. 2006, 98 (4): 757-769. 10.1002/jcb.20850.

    Article  CAS  PubMed  Google Scholar 

  12. Walkley CR, Qudsi R, Sankaran VG, Perry JA, Gostissa M, Roth SI, Rodda SJ, Snay E, Dunning P, Fahey FH, et al: Conditional mouse osteosarcoma, dependent on p53 loss and potentiated by loss of Rb, mimics the human disease. Genes Dev. 2008, 22 (12): 1662-1676. 10.1101/gad.1656808.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Gokgoz N, Wunder JS, Mousses S, Eskandarian S, Bell RS, Andrulis IL: Comparison of p53 mutations in patients with localized osteosarcoma and metastatic osteosarcoma. Cancer. 2001, 92 (8): 2181-2189. 10.1002/1097-0142(20011015)92:8<2181::AID-CNCR1561>3.0.CO;2-3.

    Article  CAS  PubMed  Google Scholar 

  14. Miller CW, Aslo A, Tsay C, Slamon D, Ishizaki K, Toguchida J, Yamamuro T, Lampkin B, Koeffler HP: Frequency and structure of p53 rearrangements in human osteosarcoma. Cancer research. 1990, 50 (24): 7950-7954.

    CAS  PubMed  Google Scholar 

  15. Radig K, Schneider-Stock R, Oda Y, Neumann W, Mittler U, Roessner A: Mutation spectrum of p53 gene in highly malignant human osteosarcomas. Gen Diagn Pathol. 1996, 142 (1): 25-32.

    CAS  PubMed  Google Scholar 

  16. Fuchs B, Pritchard DJ: Etiology of osteosarcoma. Clinical orthopaedics and related research. 2002, 397: 40-52. 10.1097/00003086-200204000-00007.

    Article  PubMed  Google Scholar 

  17. Deshpande A, Hinds PW: The retinoblastoma protein in osteoblast differentiation and osteosarcoma. Curr Mol Med. 2006, 6 (7): 809-817.

    CAS  PubMed  Google Scholar 

  18. Sadikovic B, Yoshimoto M, Al-Romaih K, Maire G, Zielenska M, Squire JA: In vitro analysis of integrated global high-resolution DNA methylation profiling with genomic imbalance and gene expression in osteosarcoma. PLoS ONE. 2008, 3 (7): e2834-10.1371/journal.pone.0002834.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Sadikovic B, Yoshimoto M, Chilton-MacNeill S, Thorner P, Squire JA, Zielenska M: Identification of interactive networks of gene expression associated with osteosarcoma oncogenesis by integrated molecular profiling. Hum Mol Genet. 2009, 18 (11): 1962-1975. 10.1093/hmg/ddp117.

    Article  CAS  PubMed  Google Scholar 

  20. Gamberi G, Benassi MS, Bohling T, Ragazzini P, Molendini L, Sollazzo MR, Pompetti F, Merli M, Magagnoli G, Balladelli A, et al: C-myc and c-fos in human osteosarcoma: prognostic value of mRNA and protein expression. Oncology. 1998, 55 (6): 556-563. 10.1159/000011912.

    Article  CAS  PubMed  Google Scholar 

  21. Wu JX, Carpenter PM, Gresens C, Keh R, Niman H, Morris JW, Mercola D: The proto-oncogene c-fos is over-expressed in the majority of human osteosarcomas. Oncogene. 1990, 5 (7): 989-1000.

    CAS  PubMed  Google Scholar 

  22. Lu XY, Lu Y, Zhao YJ, Jaeweon K, Kang J, Xiao-Nan L, Ge G, Meyer R, Perlaky L, Hicks J, et al: Cell cycle regulator gene CDC5L, a potential target for 6p12-p21 amplicon in osteosarcoma. Mol Cancer Res. 2008, 6 (6): 937-946. 10.1158/1541-7786.MCR-07-2115.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Maire G, Yoshimoto M, Chilton-MacNeill S, Thorner PS, Zielenska M, Squire JA: Recurrent RECQL4 imbalance and increased gene expression levels are associated with structural chromosomal instability in sporadic osteosarcoma. Neoplasia. 2009, 11 (3): 260-268. 263p following 268

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Davicioni E, Wai DH, Anderson MJ: Diagnostic and prognostic sarcoma signatures. Molecular Diagnosis and Therapy. 2008, 12 (6): 359-374.

    Article  CAS  PubMed  Google Scholar 

  25. Cervigne NK, Reis PP, Machado J, Sadikovic B, Bradley G, Galloni NN, Pintilie M, Jurisica I, Gilbert R, Gullane P, et al: Identification of a microRNA signature associated with progression of leukoplakia to oral carcinoma. Human Molecular Genetics. 2009

    Google Scholar 

  26. Parsons DW, Jones S, Zhang X, Lin JC, Leary RJ, Angenendt P, Mankoo P, Carter H, Siu IM, Gallia GL, et al: An integrated genomic analysis of human glioblastoma multiforme. Science. 2008, 321 (5897): 1807-1812. 10.1126/science.1164382.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. McLendon R, Friedman A, Bigner D, Van Meir EG, Brat DJ, Mastrogianakis M, Olson JJ, Mikkelsen T, Lehman N, Aldape K, et al: Comprehensive genomic characterization defines human glioblastoma genes and core pathways. Nature. 2008

    Google Scholar 

  28. Rogers S, Cambrosio A: Making a new technology work: the standardization and regulation of microarrays. Yale Journal of Biology and Medicine. 2007, 80 (4): 165-178.

    CAS  PubMed  Google Scholar 

  29. Farragher SM, Tanney A, Kennedy RD, Paul Harkin D: RNA expression analysis from formalin fixed paraffin embedded tissues. Histochemistry and Cell Biology. 2008, 130 (3): 435-445. 10.1007/s00418-008-0479-7.

    Article  CAS  PubMed  Google Scholar 

  30. Li G, Zhang W, Zeng H, Chen L, Wang W, Liu J, Zhang Z, Cai Z: An integrative multi-platform analysis for discovering biomarkers of osteosarcoma. BMC Cancer. 2009, 9: 150-10.1186/1471-2407-9-150.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Kresse SH, Ohnstad HO, Paulsen EB, Bjerkehagen B, Szuhai K, Serra M, Schaefer KL, Myklebost O, Meza-Zepeda LA: LSAMP, a novel candidate tumor suppressor gene in human osteosarcomas, identified by array comparative genomic hybridization. Genes, Chromosomes and Cancer. 2009, 48 (8): 679-693. 10.1002/gcc.20675.

    Article  CAS  PubMed  Google Scholar 

  32. Rosen G, Caparros B, Huvos AG, Kosloff C, Nirenberg A, Cacavio A, Marcove RC, Lane JM, Mehta B, Urban C: Preoperative chemotherapy for osteogenic sarcoma: selection of postoperative adjuvant chemotherapy based on the response of the primary tumor to preoperative chemotherapy. Cancer. 1982, 49 (6): 1221-1230. 10.1002/1097-0142(19820315)49:6<1221::AID-CNCR2820490625>3.0.CO;2-E.

    Article  CAS  PubMed  Google Scholar 

  33. Picci P, Sangiorgi L, Rougraff BT, Neff JR, Casadei R, Campanacci M: Relationship of chemotherapy-induced necrosis and surgical margins to local recurrence in osteosarcoma. J Clin Oncol. 1994, 12 (12): 2699-2705.

    CAS  PubMed  Google Scholar 

  34. Fritsch MK, Bridge JA, Schuster AE, Perlman EJ, Argani P: Performance characteristics of a reverse transcriptase-polymerase chain reaction assay for the detection of tumor-specific fusion transcripts from archival tissue. Pediatr Dev Pathol. 2003, 6 (1): 43-53. 10.1007/s10024-002-0013-1.

    Article  CAS  PubMed  Google Scholar 

  35. Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001, 25 (4): 402-408. 10.1006/meth.2001.1262.

    Article  CAS  PubMed  Google Scholar 

  36. Bender R, Lange S: Adjusting for multiple testing--when and how?. Journal of clinical epidemiology. 2001, 54 (4): 343-349. 10.1016/S0895-4356(00)00314-0.

    Article  CAS  PubMed  Google Scholar 

  37. Sadikovic B, Al-Romaih K, Squire JA, Zielenska M: Cause and consequences of genetic and epigenetic alterations in human cancer. Curr Genomics. 2008, 9 (6): 394-408. 10.2174/138920208785699580.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. el-Deiry WS, Harper JW, O'Connor PM, Velculescu VE, Canman CE, Jackman J, Pietenpol JA, Burrell M, Hill DE, Wang Y, et al: WAF1/CIP1 is induced in p53-mediated G1 arrest and apoptosis. Cancer research. 1994, 54 (5): 1169-1174.

    CAS  PubMed  Google Scholar 

  39. Brazier H, Stephens S, Ory S, Fort P, Morrison N, Blangy A: Expression profile of RhoGTPases and RhoGEFs during RANKL-stimulated osteoclastogenesis: identification of essential genes in osteoclasts. J Bone Miner Res. 2006, 21 (9): 1387-1398. 10.1359/jbmr.060613.

    Article  CAS  PubMed  Google Scholar 

  40. Carlo-Stella C, Lavazza C, Locatelli A, Vigano L, Gianni AM, Gianni L: Targeting TRAIL agonistic receptors for cancer therapy. Clin Cancer Res. 2007, 13 (8): 2313-2317. 10.1158/1078-0432.CCR-06-2774.

    Article  CAS  PubMed  Google Scholar 

  41. Lian JB, Javed A, Zaidi SK, Lengner C, Montecino M, van Wijnen AJ, Stein JL, Stein GS: Regulatory controls for osteoblast growth and differentiation: role of Runx/Cbfa/AML factors. Critical Reviews in Eukaryotic Gene Expression. 2004, 14 (1-2): 1-41. 10.1615/CritRevEukaryotGeneExpr.v14.i12.10.

    Article  CAS  PubMed  Google Scholar 

  42. Zhang HY, Jin L, Stilling GA, Ruebel KH, Coonse K, Tanizaki Y, Raz A, Lloyd RV: RUNX1 and RUNX2 upregulate Galectin-3 expression in human pituitary tumors. Endocrine. 2009, 35 (1): 101-111. 10.1007/s12020-008-9129-z.

    Article  CAS  PubMed  Google Scholar 

  43. Nakahara S, Oka N, Raz A: On the role of galectin-3 in cancer apoptosis. Apoptosis. 2005, 10 (2): 267-275. 10.1007/s10495-005-0801-y.

    Article  CAS  PubMed  Google Scholar 

  44. Khanna C, Khan J, Nguyen P, Prehn J, Caylor J, Yeung C, Trepel J, Meltzer P, Helman L: Metastasis-associated differences in gene expression in a murine model of osteosarcoma. Cancer research. 2001, 61 (9): 3750-3759.

    CAS  PubMed  Google Scholar 

  45. Vladimirova V, Waha A, Luckerath K, Pesheva P, Probstmeier R: Runx2 is expressed in human glioma cells and mediates the expression of galectin-3. Journal of Neuroscience Research. 2008, 86 (11): 2450-2461. 10.1002/jnr.21686.

    Article  CAS  PubMed  Google Scholar 

  46. Blyth K, Vaillant F, Hanlon L, Mackay N, Bell M, Jenkins A, Neil JC, Cameron ER: Runx2 and MYC collaborate in lymphoma development by suppressing apoptotic and growth arrest pathways in vivo. Cancer research. 2006, 66 (4): 2195-2201. 10.1158/0008-5472.CAN-05-3558.

    Article  CAS  PubMed  Google Scholar 

  47. Pregizer S, Baniwal SK, Yan X, Borok Z, Frenkel B: Progressive recruitment of Runx2 to genomic targets despite decreasing expression during osteoblast differentiation. Journal of Cellular Biochemistry. 2008, 105 (4): 965-970. 10.1002/jcb.21900.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Maruyama Z, Yoshida CA, Furuichi T, Amizuka N, Ito M, Fukuyama R, Miyazaki T, Kitaura H, Nakamura K, Fujita T, et al: Runx2 determines bone maturity and turnover rate in postnatal bone development and is involved in bone loss in estrogen deficiency. Developmental Dynamics. 2007, 236 (7): 1876-1890. 10.1002/dvdy.21187.

    Article  CAS  PubMed  Google Scholar 

  49. Pratap J, Wixted JJ, Gaur T, Zaidi SK, Dobson J, Gokul KD, Hussain S, van Wijnen AJ, Stein JL, Stein GS, et al: Runx2 transcriptional activation of Indian Hedgehog and a downstream bone metastatic pathway in breast cancer cells. Cancer research. 2008, 68 (19): 7795-7802. 10.1158/0008-5472.CAN-08-1078.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Dalla-Torre CA, Yoshimoto M, Lee CH, Joshua AM, de Toledo SR, Petrilli AS, Andrade JA, Chilton-MacNeill S, Zielenska M, Squire JA: Effects of THBS3, SPARC and SPP1 expression on biological behavior and survival in patients with osteosarcoma. BMC Cancer. 2006, 6: 237-10.1186/1471-2407-6-237.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Luo X, Chen J, Song WX, Tang N, Luo J, Deng ZL, Sharff KA, He G, Bi Y, He BC, et al: Osteogenic BMPs promote tumor growth of human osteosarcomas that harbor differentiation defects. Lab Invest. 2008, 88 (12): 1264-1277. 10.1038/labinvest.2008.98.

    Article  CAS  PubMed  Google Scholar 

  52. Lamour V, Detry C, Sanchez C, Henrotin Y, Castronovo V, Bellahcene A: Runx2- and histone deacetylase 3-mediated repression is relieved in differentiating human osteoblast cells to allow high bone sialoprotein expression. J Biol Chem. 2007, 282 (50): 36240-36249. 10.1074/jbc.M705833200.

    Article  CAS  PubMed  Google Scholar 

  53. Cleton-Jansen AM, Anninga JK, Briaire-de Bruijn IH, Romeo S, Oosting J, Egeler RM, Gelderblom H, Taminiau AH, Hogendoorn PC: Profiling of high-grade central osteosarcoma and its putative progenitor cells identifies tumourigenic pathways. British journal of cancer. 2009, 101 (12): 2064-10.1038/sj.bjc.6605482.

    Article  PubMed  PubMed Central  Google Scholar 

  54. Berger M, Muraro M, Fagioli F, Ferrari S: Osteosarcoma derived from donor stem cells carrying the Norrie's disease gene. N Engl J Med. 2008, 359 (23): 2502-2504. 10.1056/NEJMc0807172.

    Article  CAS  PubMed  Google Scholar 

  55. Tolar J, Nauta AJ, Osborn MJ, Panoskaltsis Mortari A, McElmurry RT, Bell S, Xia L, Zhou N, Riddle M, Schroeder TM, et al: Sarcoma derived from cultured mesenchymal stem cells. Stem Cells. 2007, 25 (2): 371-379. 10.1634/stemcells.2005-0620.

    Article  CAS  PubMed  Google Scholar 

Pre-publication history

Download references

Acknowledgements

This work was supported by the Canadian Cancer Society/National Cancer Institute of Canada [grant #016215]. B.S. is a Fellow of the National Cancer Institute of Canada (Terry Fox Foundation) and Restracomp Fellow at the Hospital for Sick Children, Toronto, Canada.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Jeremy Squire.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

BS: Conceived and designed experiments, performed experiments, analysed data, drafted the manuscript. PT: Performed experiments. SCM: Performed experiments. JM: Revised the manuscript. NKC: Performed experiments, analysed data. JS: Conceived and designed experiments, revised the manuscript. MZ: Conceived and designed experiments, revised the manuscript. All authors have read and approved the final manuscript.

Electronic supplementary material

Additional file 1: Primer sequences. Primer sequences for Real-Time PCR (XLS 28 KB)

12885_2009_2001_MOESM2_ESM.XLS

Additional file 2: Gene expression levels in osteosarcoma tumors and normal osteoblasts. Values of gene expression are represented as delta-delta Ct (log10). (XLS 21 KB)

12885_2009_2001_MOESM3_ESM.TIFF

Additional file 3: Controls for qRT-PCR analyses. Respectively, A2M and SLC14A were the highest- and lowest- expressed genes in a previous microarray study performed by our group [19], and they were used to validate our conditions for the qRT-PCR experiments of the current study. (TIFF 1 MB)

12885_2009_2001_MOESM4_ESM.XLS

Additional file 4: Detailed results of the statistical analyses of gene expression. Tables represent values calculated using Mann-Whitney and 1-way ANOVA statistical analyses respectively (XLS 30 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Sadikovic, B., Thorner, P., Chilton-MacNeill, S. et al. Expression analysis of genes associated with human osteosarcoma tumors shows correlation of RUNX2 overexpression with poor response to chemotherapy. BMC Cancer 10, 202 (2010). https://doi.org/10.1186/1471-2407-10-202

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2407-10-202

Keywords