Genes regulated in metastatic osteosarcoma: evaluation by microarray analysis in four human and two mouse cell line systems

ARHGAP44 expression is associated with the metastasis of osteosarcoma and is a promising prognostic biomarker

Osteosarcoma (OS) is the most domain primary malignant bone tumor. Treatment resistances and metastases result in a decreasing 5-year overall survival rate of OS. However, Rho GTPase-activating protein 44 (ARHGAP44) has not been well studied in OS. The OS patient data were obtained from Therapeutically Applicable Research to Generate Effective Treatments and Gene Expression Omnibus databases. We utilized Survival and Survminer package for survival analysis based on Kaplan-Meier method. The association between ARHGAP44 expression with the prognosis of OS was determined by Wilcoxon rank-sum test and multivariate Cox regression analysis. The real-time polymerase chain reaction and western blotting were conducted to validate the results. Gene set enrichment analysis was done to find significant Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. The OS sample ARHGAP44 expression level was significantly higher than that in normal samples, which was validated in cell lines. High ARHGAP44 expression was associated with metastasis of OS. The OS patients with high ARHGAP44 expression had worse prognosis compared with low ARHGAP44 expression OS patients. In total, 10 KEGG pathways significantly activated in high ARHGAP44 expression OS patients, such as Hedgehog signaling pathway, Steroid biosynthesis, and so on. In summary, high ARHGAP44 expression was closely correlated with the metastasis and poor prognosis of OS. ARHGAP44 was a potential prognostic biomarker for OS.

Single-cell transcriptomics reveals the regulative roles of cancer associated fibroblasts in tumor immune microenvironment of recurrent osteosarcoma

Rationale: Osteosarcoma (OS) is the most common primary bone tumor with a poor prognosis, but the detailed mechanism is still unclear. A comprehensive investigation of tumor microenvironment (TME) of OS might help find effective anti-tumor strategies. Single-cell transcriptomics is a powerful new tool to explore TME. Therefore, this study is designed to investigate the TME and gene expression pattern of primary and recurrent OS at the single-cell level.

Methods: The single-cell RNA sequencing and bioinformatic analysis were conducted to investigate the cellular constitution of primary, recurrent, and lung metastatic OS lesions according to the datasets of GSE152048 and GSE162454. TIMER database was used to investigate the role of LOX in the prognosis of sarcoma. The functions of related cells and markers were further confirmed by in vitro and in vivo experiments.

Results: Cancer associated fibroblasts (CAFs) were found with a higher infiltrating level in recurrent OS, and were enriched in the epithelial-mesenchymal transition (EMT) pathway. CAFs showed remarkably increased expression of LOX, which might lead to EMT and poor prognosis of OS. Mechanically, LOX regulated the function of CAFs and macrophage polarization to remodel the tumor immune microenvironment. Moreover, LOX inhibitor could inhibit migration and promote apoptosis of OS both in vitro and in vivo.

Conclusions: This study revealed the heterogeneity of recurrent OS and highlighted an innovative mechanism that CAFs regulate EMT of OS via LOX. Targeting LOX of CAFs showed promising efficacy in remodeling TME and treating recurrent OS.

Transcription factor c-Myb: novel prognostic factor in osteosarcoma

The transcription factor c-Myb is an oncoprotein promoting cell proliferation and survival when aberrantly activated/expressed, thus contributing to malignant transformation. Overexpression of c-Myb has been found in leukemias, breast, colon and adenoid cystic carcinoma. Recent studies revealed its expression also in osteosarcoma cell lines and suggested its functional importance during bone development. However, the relevance of c-Myb in control of osteosarcoma progression remains unknown. A retrospective clinical study was carried out to assess a relationship between c-Myb expression in archival osteosarcoma tissues and prognosis in a cohort of high-grade osteosarcoma patients. In addition, MYB was depleted in metastatic osteosarcoma cell lines SAOS-2 LM5 and 143B and their growth, chemosensitivity, migration and metastatic activity were determined. Immunohistochemical analysis revealed that high c-Myb expression was significantly associated with poor overall survival in the cohort and metastatic progression in young patients. Increased level of c-Myb was detected in metastatic osteosarcoma cell lines and its depletion suppressed their growth, colony-forming capacity, migration and chemoresistance in vitro in a cell line-dependent manner. MYB knock-out resulted in reduced metastatic activity of both SAOS-2 LM5 and 143B cell lines in immunodeficient mice. Transcriptomic analysis revealed the c-Myb-driven functional programs enriched for genes involved in the regulation of cell growth, stress response, cell adhesion and cell differentiation/morphogenesis. Wnt signaling pathway was identified as c-Myb target in osteosarcoma cells. Taken together, we identified c-Myb as a negative prognostic factor in osteosarcoma and showed its involvement in the regulation of osteosarcoma cell growth, chemosensitivity, migration and metastatic activity.

LOX upregulates FAK phosphorylation to promote metastasis in osteosarcoma

Osteosarcoma is a malignant bone tumor that commonly occurs in the pediatric population. Despite the use of chemotherapy and surgery, metastasis remains to be the leading cause of death in patients with osteosarcoma. We have previously reported that cytoplasmic mislocalization of p27 is associated with a poor outcome in osteosarcoma. In this study, we further show that lysyl oxidase (LOX) expression was associated with p27 mislocalization. LOX is an enigmatic molecule that acts as a tumor suppressor or a metastatic promoter; however, its role in osteosarcoma is still unclear. Hence, we performed both in vitro and in vivo analyses to dissect the role of LOX in osteosarcoma. The result of our survival analysis indicated that LOX expression significantly correlated with a poor outcome in osteosarcoma with or without controlling for the initial metastasis status (P < 0.05). Functionally, we found that higher LOX expression promoted osteosarcoma cell proliferation, migration, and invasiveness in vitro and produced a higher number of mice with pulmonary metastases in an orthotopic xenograft mouse model. Mechanistically, phospho-FAK was increased in osteosarcoma cells with high LOX expression. Our results further showed that FAK inhibition significantly reduced tumor cell proliferation and migration in vitro as well as LOX-mediated metastasis in mice. Together, our findings suggest that there is a novel link between p27 mislocalization and LOX expression. LOX plays a pivotal role in osteosarcoma metastasis by upregulating FAK phosphorylation. FAK inhibition may constitute a promising therapeutic strategy to reduce the development of metastasis in osteosarcoma with LOX overexpression.

Osteosarcoma cells in early and late stages as cancer in vitro progression model for assessing the responsiveness of cells to silver nanoparticles

Understanding of biology of osteosarcoma malignant progression is indispensable for enhancement of conventional chemotherapy by the use of silver nanoparticles (AgNPs). We presented an in vitro model of cancer progression closely resembling processes occurring in vivo in terms of protein profile. A comparison of cytotoxic and genotoxic potential of AgNPs in Saos-2 cells in early stages of cancerous progression (early passages) with the cells in advanced stages (late passages) demonstrated significantly reduced responsiveness of the late passage cells to nanoparticles toxicity. It was also confirmed by proteome analysis as we identified considerably higher number of differentially expressed proteins in Saos-2 cells in early passages compared to the late passage cells. Our studies showed that the ability of AgNPs as potential drug carriers to deliver a medication and/or to evoke toxic effects might be significantly diminished in advanced stages of cancer progression.

Identification of potential key genes and functional role of CENPF in osteosarcoma using bioinformatics and experimental analysis

Osteosarcoma, which arises from bone tissue, is considered to be one of the most common types of cancer in children and teenagers. As the etiology of osteosarcoma has not been fully elucidated, the overall prognosis for patients is generally poor. In recent years, the development of bioinformatical technology has allowed researchers to identify numerous molecular biological characteristics associated with the prognosis of osteosarcoma using online databases. In the present study, Gene Expression Omnibus (GEO) database was used and three microarray datasets were obtained. The GEO2R web tool was utilized and differentially expressed genes (DEGs) in osteosarcoma tissue were identified. Venn analysis was performed to determine the intersection of the DEG profiles. DEGs were analyzed by Gene Ontology function and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis. Protein-protein interactions (PPIs) between these DEGs were analyzed using the Search Tool for the Retrieval of Interacting Genes database, and the PPI network was then visualized using Cytoscape software. The top ten genes were identified based on measurement of degree, density of maximum neighborhood component, maximal clique centrality and mononuclear cell counts in the PPI network, and five overlapping genes [origin recognition complex subunit 6 (ORC6), IGF-binding protein 5 (IGFBP5), minichromosome maintenance 10 replication initiation factor (MCM10), MET proto-oncogene, receptor tyrosine kinase (MET) and centromere protein F (CENPF)] were identified. Additionally, three module networks were analyzed by Molecular Complex Detection (MCODE), and six key genes [ORC6, MCM10, DEP domain containing 1 (DEPDC1), CENPF, TIMELESS interacting protein (TIPIN) and shugoshin 1 (SGOL1)] were screened. Combined with the results from Cytoscape and MCODE, eight hub genes (ORC6, MCM10, DEPDC1, CENPF, TIPIN, SGOL1, MET and IGFBP5) were obtained. Furthermore, Kaplan-Meier plotter survival analysis was used to evaluate the prognostic value of these eight hub genes in patients with osteosarcoma. Oncomine and GEPIA databases were applied to further confirm the expression levels of hub genes in tissue. Finally, the functional roles of the core gene CENPF were investigated using Cell Counting Kit-8, wound healing and Transwell assays, which indicated that CENPF knockdown inhibited the proliferation, migration and invasion of osteosarcoma cells. These results provided potential prognostic markers, as well as a basis for further investigation of the mechanism underlying osteosarcoma.

Identification of abnormally high expression of POGZ as a new biomarker associated with a poor prognosis in osteosarcoma

Osteosarcoma (OS) is the most prevalent malignant bone tumor in children and young adults. There is an urgent need for a novel biomarker related to the prognosis of OS. We performed a meta-analysis incorporating six independent datasets and performed a survival analysis with one independent dataset GSE21257 in the GEO database for gene screening. The results revealed that one potential biomarker related to OS survival, POGZ was the most significantly upregulated gene. We also verified that the POGZ was overexpressed in clinical samples. The survival analysis revealed that POGZ is associated with a poor prognosis in OS. Moreover, flow cytometry analysis of isolated OS cells demonstrated that OS cells were arrested in the G1 phase after POGZ knockdown. The RNA-seq results indicated that POGZ was co-expressed with CCNE1 and CCNB1. Pathway analysis showed that genes associated with high expression levels of POGZ were related to the cell cycle pathway. A cell model was constructed to detect the effects of POGZ. After POGZ knockdown, OS cell proliferation, invasion and migration were all decreased. Therefore, POGZ is an important gene for evaluating the prognosis of OS patients and is a potential therapeutic target.

RNA sequencing and functional studies of patient-derived cells reveal that neurexin-1 and regulators of this pathway are associated with poor outcomes in Ewing sarcoma

Purpose

The development of biomarkers and molecularly targeted therapies for patients with Ewing sarcoma (ES) in order to minimise morbidity and improve outcome is urgently needed. Here, we set out to isolate and characterise patient-derived ES primary cell cultures and daughter cancer stem-like cells (CSCs) to identify biomarkers of high-risk disease and candidate therapeutic targets.

Methods

Thirty-two patient-derived primary cultures were established from treatment-naïve tumours and primary ES-CSCs isolated from these cultures using functional methods. By RNA-sequencing we analysed the transcriptome of ES patient-derived cells (n = 24) and ES-CSCs (n = 11) to identify the most abundant and differentially expressed genes (DEGs). Expression of the top DEG(s) in ES-CSCs compared to ES cells was validated at both RNA and protein levels. The functional and prognostic potential of the most significant gene (neurexin-1) was investigated using knock-down studies and immunohistochemistry of two independent tumour cohorts.

Results

ES-CSCs were isolated from all primary cell cultures, consistent with the premise that ES is a CSC driven cancer. Transcriptional profiling confirmed that these cells were of mesenchymal origin, revealed novel cell surface targets for therapy that regulate cell-extracellular matrix interactions and identified candidate drivers of progression and relapse. High expression of neurexin-1 and low levels of regulators of its activity, APBA1 and NLGN4X, were associated with poor event-free and overall survival rates. Knock-down of neurexin-1 decreased viable cell numbers and spheroid formation.

Conclusions

Genes that regulate extracellular interactions, including neurexin-1, are candidate therapeutic targets in ES. High levels of neurexin-1 at diagnosis are associated with poor outcome and identify patients with localised disease that will relapse. These patients could benefit from more intensive or novel treatment modalities. The prognostic significance of neurexin-1 should be validated independently.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13402-021-00619-8.

Targeting Mechanotransduction in Osteosarcoma: A Comparative Oncology Perspective

Mechanotransduction is the process in which cells can convert extracellular mechanical stimuli into biochemical changes within a cell. While this a normal process for physiological development and function in many organ systems, tumour cells can exploit this process to promote tumour progression. Here we summarise the current state of knowledge of mechanotransduction in osteosarcoma (OSA), the most common primary bone tumour, referencing both human and canine models and other similar mesenchymal malignancies (e.g., Ewing sarcoma). Specifically, we discuss the mechanical properties of OSA cells, the pathways that these cells utilise to respond to external mechanical cues, and mechanotransduction-targeting strategies tested in OSA so far. We point out gaps in the literature and propose avenues to address them. Understanding how the physical microenvironment influences cell signalling and behaviour will lead to the improved design of strategies to target the mechanical vulnerabilities of OSA cells.

Molecular Biology of Osteosarcoma

Osteosarcoma (OS) is the most frequent primary bone cancer in children and adolescents and the third most frequent in adults. Many inherited germline mutations are responsible for syndromes that predispose to osteosarcomas including Li Fraumeni syndrome, retinoblastoma syndrome, Werner syndrome, Bloom syndrome or Diamond-Blackfan anemia. TP53 is the most frequently altered gene in osteosarcoma. Among other genes mutated in more than 10% of OS cases, c-Myc plays a role in OS development and promotes cell invasion by activating MEK-ERK pathways. Several genomic studies showed frequent alterations in the RB gene in pediatric OS patients. Osteosarcoma driver mutations have been reported in NOTCH1, FOS, NF2, WIF1, BRCA2, APC, PTCH1 and PRKAR1A genes. Some miRNAs such as miR-21, -34a, -143, -148a, -195a, -199a-3p and -382 regulate the pathogenic activity of MAPK and PI3K/Akt-signaling pathways in osteosarcoma. CD133+ osteosarcoma cells have been shown to exhibit stem-like gene expression and can be tumor-initiating cells and play a role in metastasis and development of drug resistance. Although currently osteosarcoma treatment is based on adriamycin chemoregimens and surgery, there are several potential targeted therapies in development. First of all, activity and safety of cabozantinib in osteosarcoma were studied, as well as sorafenib and pazopanib. Finally, novel bifunctional molecules, of potential imaging and osteosarcoma targeting applications may be used in the future.

A risk signature-based on metastasis-associated genes to predict survival of patients with osteosarcoma

Osteosarcoma (OS) is the most common primary solid malignant bone tumor, and its metastasis is a prominent cause of high mortality in patients. In this study, a prognosis risk signature was constructed based on metastasis-associated genes. Four microarrays datasets with clinical information were downloaded from Gene Expression Omnibus, and 256 metastasis-associated genes were identified by limma package. Further, a protein-protein interaction network was constructed, and survival analysis was performed using data from the Therapeutically Applicable Research to Generate Effective Treatments data matrix, identifying 19 genes correlated with prognosis. Six genes were selected by the least absolute shrinkage and selection operator regression for multivariate cox analysis. Finally, a three-gene (MYC, CPE, and LY86) risk signature was constructed, and datasets GSE21257 and GSE16091 were used to validate the prediction efficiency of the signature. The survival times of low- and high-risk groups were significantly different in the training set and validation set. Additionally, gene set enrichment analysis revealed that the genes in the signature may affect the cell cycle, gap junctions, and interleukin-6 production. Therefore, the three-gene survival risk signature could potentially predict the prognosis of patients with OS. Further, proteins encoded by CPE and LY86 may provide novel insights into the prediction of OS prognosis and therapeutic targets.

Metastatic immune infiltrates correlate with those of the primary tumour in canine osteosarcoma

Our lack of understanding of the immune microenvironment in canine osteosarcoma (cOSA) has limited the identification of potential immunotherapeutic targets. In particular, our ability to utilize readily available tissue from a dog's primary tumour to predict the type and extent of immune response in their pulmonary metastatic lesions is unknown. We, therefore, collected 21 matched pairs of primary tumours and pulmonary metastatic lesions from dogs with OSA and performed immunohistochemistry to quantify T-lymphocyte (CD3), FOXP3+ cell, B-lymphocyte (Pax-5), and CD204+ macrophage infiltration. We found that T-lymphocytes and FOXP3+ infiltrates in primary tumours positively correlated with that of metastatic lesions (ρ = 0.512, P = 0.038 and ρ = 0.698, P = 0.007, respectively), while a strong trend existed for CD204+ infiltrates (ρ = 0.404, P = 0.087). We also observed T- and B-lymphocytes, and CD204+ macrophages to be significantly higher in a dog's pulmonary metastasis compared to their primary tumour (P = 0.018, P = 0.018, P = 0.016, respectively), while FOXP3+ cells were only significantly higher in metastases when all primary tumour and metastasis lesions were compared without pairing (P = 0.036). Together, these findings suggest that the metastatic immune microenvironment may be influenced by that of the primary cOSA, and that primary tumour immune biomarkers could potentially be applied to predict immunotherapeutic responses in gross metastatic disease. We, therefore, provide a rationale for the treatment of cOSA pulmonary metastases with immunotherapeutics that enhance the anti-tumour activity of these immune cells, particularly in dogs with moderate to high immune cell infiltration in their primary tumours.

The relationship between metastatic potential and in vitro mechanical properties of osteosarcoma cells

Osteosarcoma is the most frequent primary tumor of bone and is characterized by its high tendency to metastasize in lungs. Although treatment in cases of early diagnosis results in a 5-yr survival rate of nearly 60%, the prognosis for patients with secondary lesions at diagnosis is poor, and their 5-yr survival rate remains below 30%. In the present work, we have used a number of analytical methods to investigate the impact of increased metastatic potential on the biophysical properties and force generation of osteosarcoma cells. With that aim, we used two paired osteosarcoma cell lines, with each one comprising a parental line with low metastatic potential and its experimentally selected, highly metastatic form. Mechanical characterization was performed by means of atomic force microscopy, tensile biaxial deformation, and real-time deformability, and cell traction was measured using two-dimensional and micropost-based traction force microscopy. Our results reveal that the low metastatic osteosarcoma cells display larger spreading sizes and generate higher forces than the experimentally selected, highly malignant variants. In turn, the outcome of cell stiffness measurements strongly depends on the method used and the state of the probed cell, indicating that only a set of phenotyping methods provides the full picture of cell mechanics.

A risk score model for the prediction of osteosarcoma metastasis

Osteosarcoma is the most common primary solid malignancy of the bone, and its high mortality usually correlates with early metastasis. In this study, we developed a risk score model to help predict metastasis at the time of diagnosis. We downloaded and mined four expression profile datasets associated with osteosarcoma metastasis from the Gene Expression Omnibus. After data normalization, we performed LASSO logistic regression analysis together with 10-fold cross validation using the GSE21257 dataset. A combination of eight genes (RAB1,CLEC3B,FCGBP,RNASE3,MDL1,ALOX5AP,VMO1 and ALPK3) were identified as being associated with osteosarcoma metastasis. These genes were put into a gene risk score model, and the prediction efficiency of the model was then validated using three independent datasets (GSE33383, GSE66673, and GSE49003) by plotting receiver operating characteristic curves. The expression levels of the eight genes in all datasets were shown as heatmaps, and gene ontology gene annotation and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis were performed. These eight genes play a role in cancer-related biological processes, such as apoptosis and biosynthetic processes. Our results may aid in elucidating the possible mechanisms of osteosarcoma metastasis, and may help to facilitate the individual management of patients with osteosarcoma after treatment.

Identifying osteosarcoma metastasis associated genes by weighted gene co-expression network analysis (WGCNA)

Osteosarcoma (OS), the most common malignant bone tumor, accounts for the heavy healthy threat in the period of children and adolescents. OS occurrence usually correlates with early metastasis and high death rate. This study aimed to better understand the mechanism of OS metastasis.Based on Gene Expression Omnibus (GEO) database, we downloaded 4 expression profile data sets associated with OS metastasis, and selected differential expressed genes. Weighted gene co-expression network analysis (WGCNA) approach allowed us to investigate the most OS metastasis-correlated module. Gene Ontology functional and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were used to give annotation of selected OS metastasis-associated genes.We select 897 differential expressed genes from OS metastasis and OS non-metastasis groups. Based on these selected genes, WGCNA further explored 142 genes included in the most OS metastasis-correlated module. Gene Ontology functional and KEGG pathway enrichment analyses showed that significantly OS metastasis-associated genes were involved in pathway correlated with insulin-like growth factor binding.Our research figured out several potential molecules participating in metastasis process and factors acting as biomarker. With this study, we could better explore the mechanism of OS metastasis and further discover more therapy targets.

Bisphosphonate-Functionalized Hydroxyapatite Nanoparticles for the Delivery of the Bromodomain Inhibitor JQ1 in the Treatment of Osteosarcoma

Osteosarcoma (OS) is one of the most common neoplasia among children, and its survival statistics have been stagnating since the combinatorial anticancer therapy triad was first introduced. Here, we report on the assessment of the effect of hydroxyapatite (HAp) nanoparticles loaded with medronate, the simplest bisphosphonate, as a bone-targeting agent and JQ1, a small-molecule bromodomain inhibitor, as a chemotherapeutic in different 2D and 3D K7M2 OS in vitro models. Both additives decreased the crystallinity of HAp, but the effect was more intense for medronate because of its higher affinity for HAp. As the result of PO₄^3--NH⁺ binding, JQ1 shielded the surface phosphates of HAp and pushed its surface charge to more positive values, exhibiting the opposite effect from calcium-blocking medronate. In contrast to the faster and more exponential release of JQ1 from monetite, its release from HAp nanoparticles followed a zero-order kinetics, but 98% of the payload was released after 48 h. The apoptotic effect of HAp nanoparticles loaded with JQ1, with medronate and with both JQ1 and medronate, was selective in 2D culture: pronounced against the OS cells and nonexistent against the healthy fibroblasts. While OS cell invasion was significantly inhibited by all of the JQ1-containing HAp formulations, that is, with and without medronate, all of the combinations of the targeting compound, medronate, and the chemotherapeutic, JQ1, delivered using HAp, but not HAp alone, inhibited OS cell migration from the tumor spheroids. JQ1 delivered using HAp had an effect on tumor migration, invasion, and apoptosis even at extremely low, subnanomolar concentrations, at which no effect of JQ1 per se was observed, meaning that this form of delivery could help achieve a multifold increase of this drug's efficacy. More than 80% of OS cells internalized JQ1-loaded HAp nanoparticles after 24 h of coincubation, suggesting that this augmentation of the activity of JQ1 may be due to the intracellular delivery and sustained release of the drug enabled by HAp. In addition to the reduction of the OS cell viability, the reduction of the migration and invasion radii was observed in OS tumor spheroids challenged with even JQ1-free medronate-functionalized HAp nanoparticles, demonstrating a definite anticancer activity of medronate alone when combined with HAp. The effect of medronate-functionalized JQ1-loaded HAp nanoparticles was most noticeable against OS cells differentiated into an osteoblastic lineage, in which case they surpassed in effect pure JQ1 and medronate-free compositions. The activity of JQ1 was mediated through increased Ezrin expression and decreased RUNX2 expression and was MYC and FOSL1 independent, but these patterns of gene expression changed in cells challenged with the nanoparticulate form of delivery, having been accompanied by the upregulation of RUNX2 and downregulation of Ezrin in OS cells treated with medronate-functionalized JQ1-loaded HAp nanoparticles.

The critical role of the ZNF217 oncogene in promoting breast cancer metastasis to the bone

Bone metastasis affects >70% of patients with advanced breast cancer. However, the molecular mechanisms underlying this process remain unclear. On the basis of analysis of clinical datasets, and in vitro and in vivo experiments, we report that the ZNF217 oncogene is a crucial mediator and indicator of bone metastasis. Patients with high ZNF217 mRNA expression levels in primary breast tumours had a higher risk of developing bone metastases. MDA-MB-231 breast cancer cells stably transfected with ZNF217 (MDA-MB-231-ZNF217) showed the dysregulated expression of a set of genes with bone-homing and metastasis characteristics, which overlapped with two previously described 'osteolytic bone metastasis' gene signatures, while also highlighting the bone morphogenetic protein (BMP) pathway. The latter was activated in MDA-MB-231-ZNF217 cells, and its silencing by inhibitors (Noggin and LDN-193189) was sufficient to rescue ZNF217-dependent cell migration, invasion or chemotaxis towards the bone environment. Finally, by using non-invasive multimodal in vivo imaging, we found that ZNF217 increases the metastatic growth rate in the bone and accelerates the development of severe osteolytic lesions. Altogether, the findings of this study highlight ZNF217 as an indicator of the emergence of breast cancer bone metastasis; future therapies targeting ZNF217 and/or the BMP signalling pathway may be beneficial by preventing the development of bone metastases. Copyright © 2017 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Characterization of the metastatic phenotype of a panel of established osteosarcoma cells

Osteosarcoma (OS) is the most common bone tumor in pediatric patients. Metastasis is a major cause of mortality and morbidity. The rarity of this disease coupled with the challenges of drug development for metastatic cancers have slowed the delivery of improvements in long-term outcomes for these patients. In this study, we collected 18 OS cell lines, confirmed their expression of bone markers and complex karyotypes, and characterized their in vivo tumorgenicity and metastatic potential. Since prior reports included conflicting descriptions of the metastatic and in vivo phenotypes of these models, there was a need for a comparative assessment of metastatic phenotypes using identical procedures in the hands of a single investigative group. We expect that this single characterization will accelerate the study of this metastatic cancer. Using these models we evaluated the expression of six previously reported metastasis-related OS genes. Ezrin was the only gene consistently differentially expressed in all the pairs of high/low metatstatic OS cells. We then used a subtractive gene expression approach of the high and low human metastatic cells to identify novel genes that may be involved in OS metastasis. PHLDA1 (pleckstrin homology-like domain, family A) was identified as one of the genes more highly expressed in the high metastatic compared to low metastatic cells. Knocking down PHLDA1 with siRNA or shRNA resulted in down regulation of the activities of MAPKs (ERK1/2), c-Jun N-terminal kinases (JNK), and p38 mitogen-activated protein kinases (MAPKs). Reducing the expression of PHLDA1 also delayed OS metastasis progression in mouse xenograft models.

Bioinformatics analysis of differentially expressed pathways related to the metastatic characteristics of osteosarcoma

In this study, gene expression data of osteosarcoma (OSA) were analyzed to identify metastasis-related biological pathways. Four gene expression data sets (GSE21257, GSE9508, GSE49003 and GSE66673) were downloaded from Gene Expression Omnibus (GEO). An analysis of differentially expressed genes (DEGs) was performed using the Significance Analysis of Microarray (SAM) method. Gene expression levels were converted into scores of pathways by the Functional Analysis of Individual Microarray Expression (FAIME) algorithm and the differentially expressed pathways (DEPs) were then disclosed by a t-test. The distinguishing and prediction ability of the DEPs for metastatic and non-metastatic OSA was further confirmed using the principal component analysis (PCA) method and 3 gene expression data sets (GSE9508, GSE49003 and GSE66673) based on the support vector machines (SVM) model. A total of 616 downregulated and 681 upregulated genes were identified in the data set, GSE21257. The DEGs could not be used to distinguish metastatic OSA from non-metastatic OSA, as shown by PCA. Thus, an analysis of DEPs was further performed, resulting in 14 DEPs, such as NRAS signaling, Toll-like receptor (TLR) signaling, matrix metalloproteinase (MMP) regulation of cytokines and tumor necrosis factor receptor-associated factor (TRAF)-mediated interferon regulatory factor 7 (IRF7) activation. Cluster analysis indicated that these pathways could be used to distinguish between metastatic OSA from non-metastatic OSA. The prediction accuracy was 91, 66.7 and 87.5% for the data sets, GSE9508, GSE49003 and GSE66673, respectively. The results of PCA further validated that the DEPs could be used to distinguish metastatic OSA from non-metastatic OSA. On the whole, several DEPs were identified in metastatic OSA compared with non-metastatic OSA. Further studies on these pathways and relevant genes may help to enhance our understanding of the molecular mechanisms underlying metastasis and may thus aid in the development of novel therapies.

MicroRNAs in the pathobiology of sarcomas

Sarcomas are a rare and heterogeneous group of tumors. The last decade has witnessed extensive efforts to understand the pathobiology of many aggressive sarcoma types. In parallel, we have also begun to unravel the complex gene regulation processes mediated by microRNAs (miRNAs) in sarcomas and other cancers, discovering that microRNAs have critical roles in the majority of both oncogenic and tumor suppressor signaling networks. Expression profiles and a greater understanding of the biologic roles of microRNAs and other noncoding RNAs have considerably expanded our current knowledge and provided key pathobiological insights into many sarcomas, and helped identify novel therapeutic targets. The limited number of sarcoma patients in each sarcoma type and their heterogeneity pose distinct challenges in translating this knowledge into the clinic. It will be critical to prioritize these novel targets and choose those that have a broad applicability. A small group of microRNAs have conserved roles across many types of sarcomas and other cancers. Therapies that target these key microRNA-gene signaling and regulatory networks, in combination with standard of care treatment, may be the pivotal component in significantly improving treatment outcomes in patients with sarcoma or other cancers.

Genomic instability of osteosarcoma cell lines in culture: impact on the prediction of metastasis relevant genes

Background

Osteosarcoma is a rare but highly malignant cancer of the bone. As a consequence, the number of established cell lines used for experimental in vitro and in vivo osteosarcoma research is limited and the value of these cell lines relies on their stability during culture. Here we investigated the stability in gene expression by microarray analysis and array genomic hybridization of three low metastatic cell lines and derivatives thereof with increased metastatic potential using cells of different passages.

Principal Findings

The osteosarcoma cell lines showed altered gene expression during in vitro culture, and it was more pronounced in two metastatic cell lines compared to the respective parental cells. Chromosomal instability contributed in part to the altered gene expression in SAOS and LM5 cells with low and high metastatic potential. To identify metastasis-relevant genes in a background of passage-dependent altered gene expression, genes involved in "Pathways in cancer" that were consistently regulated under all passage comparisons were evaluated. Genes belonging to "Hedgehog signaling pathway" and "Wnt signaling pathway" were significantly up-regulated, and IHH, WNT10B and TCF7 were found up-regulated in all three metastatic compared to the parental cell lines.

Conclusions

Considerable instability during culture in terms of gene expression and chromosomal aberrations was observed in osteosarcoma cell lines. The use of cells from different passages and a search for genes consistently regulated in early and late passages allows the analysis of metastasis-relevant genes despite the observed instability in gene expression in osteosarcoma cell lines during culture.

Tumourigenic canine osteosarcoma cell lines associated with frizzled-6 up-regulation and enhanced side population cell frequency

An increased serum alkaline phosphatase concentration is known to be associated with a negative prognosis in canine and human osteosarcoma. To expand upon previous studies regarding the biological relevance of increased serum alkaline phosphatase as a negative prognostic factor, xenogeneic heterotopic transplants were performed using six canine primary osteosarcoma cell lines generated from patients with differing serum alkaline phosphatase concentrations (three normal and three increased). Three of the six cell lines were capable of generating tumours and tumour formation was independent of the serum alkaline phosphatase status of the cell line. Microarray analysis identified 379 genes as being differentially expressed between the tumourigenic and non-tumourigenic cell lines. Frizzled-6 was upregulated to the greatest extent (7.78-fold) in tumourigenic cell lines compared with non-tumourigenic cell lines. Frizzled-6, a co-receptor for Wnt ligands has been associated with enhanced tumour-initiating cells and poor prognosis for other tumours. The increased expression of frizzled-6 was confirmed by quantitative reverse transcription polymerase chain reaction (QPCR) and Western blot analysis. Additionally, the tumourigenic cell lines also had an increase in the percentage of side population cells compared with non-tumourigenic cell lines (5.89% versus 1.58%, respectively). There were no differences in tumourigenicity, frizzled-6 or percentage of side population cells noted between osteosarcoma cell lines generated from patients of differing serum alkaline phosphatase concentration. However, to our knowledge this is the first study to identified frizzled-6 as a possible marker of osteosarcoma cell populations with enhanced tumourigenicity and side population cells. Future work will focus on defining the role of frizzled-6 in osteosarcoma tumourigenesis and tumour-initiating cells.

NKD2, a negative regulator of Wnt signaling, suppresses tumor growth and metastasis in osteosarcoma

Osteosarcoma (OS) is the most frequent pediatric malignant bone tumor that has a high propensity for metastases. Through osteoblast-specific alteration of p53 status, we developed a genetically engineered mouse model of localized and metastatic OS to gain an understanding into the molecular pathogenesis of OS. Microarray analysis of both localized tumors and metastatic tumors identified the downregulation of the naked cuticle homolog 2 (NKD2) gene, a negative regulator of Wnt signaling. Overexpression of NKD2 in metastatic human and mouse OS cells significantly decreases cell proliferation, migration and invasion ability in vitro and drastically diminishes OS tumor growth and metastasis in vivo, whereas downregulation enhances migratory and invasive potential. Evaluation of NKD2-overexpressing tumors revealed upregulation of tumor-suppressor genes and downregulation of molecules involved in blood vessel formation and cell migration. Furthermore, assessment of primary human OS revealed downregulation of NKD2 in metastatic and recurrent OS. Finally, we provide biological evidence that use of small-molecule inhibitors targeting the Wnt pathway can have therapeutic efficacy in decreasing metastatic properties in OS. Our studies provide compelling evidence that downregulation of NKD2 expression and alterations in associated regulated pathways have a significant role in driving OS tumor growth and metastasis.

Reduced Latency in the Metastatic Niche Contributes to the More Aggressive Phenotype of LM8 Compared to Dunn Osteosarcoma Cells

Metastasis is the major cause of death of osteosarcoma patients and its diagnosis remains difficult. In preclinical studies, however, forced expression of reporter genes in osteosarcoma cells has remarkably improved the detection of micrometastases and, consequently, the quality of the studies. We recently showed that Dunn cells equipped with a lacZ reporter gene disseminated from subcutaneous primary tumors as frequently as their highly metastatic subline LM8, but only LM8 cells grew to macrometastases. In the present time-course study, tail-vein-injected Dunn and LM8 cells settled within 24 h at the same frequency in the lung, liver, and kidney of mice. Furthermore, Dunn cells also grew to macrometastases, but, compared to LM8, with a delay of two weeks in lung and one week in liver and kidney tissue, consistent with prolonged survival of the mice. Dunn- and LM8-cell-derived ovary and spine metastases occurred less frequently. In vitro, Dunn cells showed less invasiveness and stronger contact inhibition and intercellular adhesion than LM8 cells and several cancer- and dormancy-related genes were differentially expressed. In conclusion, Dunn cells, compared to LM8, have a similar capability but a longer latency to form macrometastases and provide an interesting new experimental system to study tumor cell dormancy.