Identification of osteosarcoma driver genes by integrative analysis of copy number and gene expression data

Aneuploidy as a driver of human cancer

Aneuploidy, an abnormal chromosome composition, is a major contributor to cancer development and progression and an important determinant of cancer therapeutic responses and clinical outcomes. Despite being recognized as a hallmark of human cancer, the exact role of aneuploidy as a 'driver' of cancer is still largely unknown. Identifying the specific genetic elements that underlie the recurrence of common aneuploidies remains a major challenge of cancer genetics. In this Review, we discuss recurrent aneuploidies and their function as drivers of tumor development. We then delve into the context-dependent identification and functional characterization of the driver genes underlying driver aneuploidies and examine emerging strategies to uncover these driver genes using cancer genomics data and cancer models. Lastly, we explore opportunities for targeting driver aneuploidies in cancer by leveraging the functional consequences of these common genetic alterations.

Genomic Landscape of Osteosarcoma of Bone in an Older-Aged Patient Population and Analysis of Possible Etiologies Based on Molecular Signature

Background: Osteosarcoma (OS), the most common primary malignant bone tumor, occurs mostly in the pediatric and adolescent (P/A) population where it has been subject to intense study whereas OS arising in the older-aged adult population has undergone less scrutiny. Materials and Methods: In this study, we assess the molecular aberrations detected in eight older adult patients (>59 years of age) with OS of bone by whole-exome sequencing (WES) on formalin-fixed, paraffin-embedded tissue and quantified the contributions of endogenous and exogenous mutational processes to tumor mutational burden and to tumorigenesis through computational analysis. Results: We identified 86 clinically significant somatic mutations. TP53 mutations occurred in OSs of three patients and one patient harbored a pathogenic germline mutation of TP53. Loss-of-heterozygosity of DNA-damage repair genes occurred in all six tumors evaluated. Computational analysis of single nucleotide variants within each tumor detected eight distinct mutagenic processes of which age-associated mutational processes, thiopurine chemotherapy, and defective homologous DNA recombination repair contributed the most to both tumor mutation burden and tumor pathogenesis. Conclusion: The genomic landscape of our older OS patients deciphered by WES is extremely diverse with only 15% of mutated somatic genes uncovered in our study previously described in P/A-enriched OS studies. Endogenous age-related mutagenic processes, defective DNA homologous recombination repair, and exogenous effects of chemotherapy are mainly responsible for pathogenic mutations in OS occurring in our cohort.

Targeting WNT5B and WNT10B in osteosarcoma

WNT signaling regulates osteosarcoma proliferation. However, there is controversy in the field of osteosarcoma as to whether WNT signaling is pro- or anti-tumorigenic. WNT-targeting therapeutics, both activators and inhibitors, are compared. WNT5B, a β-catenin-independent ligand, and WNT10B, a β-catenin-dependent WNT ligand, are each expressed in osteosarcomas, but they are not expressed in the same tumors. Furthermore, WNT10B and WNT5B regulate different histological subtypes of osteosarcomas. Using WNT signaling modulators as therapeutics may depend on the WNT ligand and/or the activated signaling pathway.

WNT5B drives osteosarcoma stemness, chemoresistance and metastasis

BACKGROUND

Treatment for osteosarcoma, a paediatric bone cancer with no therapeutic advances in over three decades, is limited by a lack of targeted therapies. Osteosarcoma frequently metastasises to the lungs, and only 20% of patients survive 5 years after the diagnosis of metastatic disease. We found that WNT5B is the most abundant WNT expressed in osteosarcoma tumours and its expression correlates with metastasis, histologic subtype and reduced survival.

METHODS

Using tumor-spheroids to model cancer stem-like cells, we performed qPCR, immunoblotting, and immunofluorescence to monitor changes in gene and protein expression. Additionally, we measured sphere size, migration and forming efficiency to monitor phenotypic changes. Therefore, we characterised WNT5B's relevance to cancer stem-like cells, metastasis, and chemoresistance and evaluated its potential as a therapeutic target.

RESULTS

In osteosarcoma cell lines and patient-derived spheres, WNT5B is enriched in stem cells and induces the expression of the stemness gene SOX2. WNT5B promotes sphere size, sphere-forming efficiency, and cell proliferation, migration, and chemoresistance to methotrexate (but not cisplatin or doxorubicin) in spheres formed from conventional cell lines and patient-derived xenografts. In vivo, WNT5B increased osteosarcoma lung and liver metastasis and inhibited the glycosaminoglycan hyaluronic acid via upregulation of hyaluronidase 1 (HYAL1), leading to changes in the tumour microenvironment. Further, we identified that WNT5B mRNA and protein correlate with the receptor ROR1 in primary tumours. Targeting WNT5B through inhibition of WNT/ROR1 signalling with an antibody to ROR1 reduced stemness properties, including chemoresistance, sphere size and SOX2 expression.

CONCLUSIONS

Together, these data define WNT5B's role in driving osteosarcoma cancer stem cell expansion and methotrexate resistance and provide evidence that the WNT5B pathway is a promising candidate for treating osteosarcoma patients.

KEY POINTS

WNT5B expression is high in osteosarcoma stem cells leading to increased stem cell proliferation and migration through SOX2. WNT5B expression in stem cells increases rates of osteosarcoma metastasis to the lungs and liver in vivo. The hyaluronic acid degradation enzyme HYAL1 is regulated by WNT5B in osteosarcoma contributing to metastasis. Inhibition of WNT5B with a ROR1 antibody decreases osteosarcoma stemness.

Integrative transcriptome analysis identifies MYBL2 as a poor prognosis marker for osteosarcoma and a pan-cancer marker of immune infiltration

MYBL2 (MYB proto-oncogene like 2) is an emerging prognostic marker for malignant tumors, and its potential role in osteosarcoma and its relationship with immune infiltration in pan-cancer is yet to be elucidated. We constructed a transcription factor activity profile of osteosarcoma using the single-cell regulatory network inference algorithm based on single-cell RNA sequencing data obtained from the Gene Expression Omnibus. Subsequently, we calculated the extent of MYBL2 activation in malignant proliferative osteoblasts. We also explored the association between MYBL2 and chemotherapy resistance in osteosarcoma. Furthermore, we systematically correlated MYBL2 with immunological signatures in the tumor microenvironment in pan-cancer, including immune cell infiltration, immune checkpoints, and tumor immunotherapy prognosis. Finally, we developed and validated a risk score (MRGS), derived an osteosarcoma risk score nomogram based on MRGS, and tested its ability to predict prognosis. MYBL2 and gene enrichment analyses in osteosarcoma and pan-cancer revealed that MYBL2 was positively correlated with cell proliferation and tumor immune pathways. MYBL2 expression positively correlated with SLC19A1 in pan-cancer and osteosarcoma cell lines. Pan-cancer immune infiltration analysis revealed that MYBL2 was correlated with myeloid-derived suppressor cells, Th2 cell infiltration, CD276, RELT gene expression, and tumor mutation burden. In summary, MYBL2 regulates proliferation, progression, and immune infiltration in osteosarcoma and pan-cancer. Therefore, we found that MYBL2 could be used as a potential marker for predicting the osteosarcoma prognosis. Patients with osteosarcoma and high MYBL2 expression are theoretically more sensitive to methotrexate. An osteosarcoma prognostic nomogram can provide new ideas in the search for osteosarcoma prognostic markers.

Osteosarcoma cells exhibit functional interactions with stromal cells, fostering a lung microenvironment conducive to the establishment of metastatic tumor cells

BACKGROUND

Osteosarcoma (OS) stands out as the most common bone tumor, with approximately 20% of the patients receiving a diagnosis of metastatic OS at their initial assessment. A significant challenge lies in the frequent existence of undetected metastases during the initial diagnosis. Mesenchymal stem cells (MSCs) possess unique abilities that facilitate tumor growth, and their interaction with OS cells is crucial for metastatic spread.

METHODS AND RESULTS

We demonstrated that, in vitro, MSCs exhibited a heightened migration response toward the secretome of non-metastatic OS cells. When challenged to a secretome derived from lungs preloaded with OS cells, MSCs exhibited greater migration toward lungs colonized with metastatic OS cells. Moreover, in vivo, MSCs displayed preferential migratory and homing behavior toward lungs colonized by metastatic OS cells. Metastatic OS cells, in turn, demonstrated an increased migratory response to the MSCs' secretome. This behavior was associated with heightened cathepsin D (CTSD) expression and the release of active metalloproteinase 2 (MMP2) by metastatic OS cells.

CONCLUSIONS

Our assessment focused on two complementary tumor capabilities crucial to metastatic spread, emphasizing the significance of inherent cell features. The findings underscore the pivotal role of signaling integration within the niche, with a complex interplay of migratory responses among established OS cells in the lungs, prometastatic OS cells in the primary tumor, and circulating MSCs. Pulmonary metastases continue to be a significant factor contributing to OS mortality. Understanding these mechanisms and identifying differentially expressed genes is essential for pinpointing markers and targets to manage metastatic spread and improve outcomes for patients with OS.

SKP2 Knockout in Rb1/p53-Deficient Mouse Models of Osteosarcoma Induces Immune Infiltration and Drives a Transcriptional Program with a Favorable Prognosis

Osteosarcoma is an aggressive bone malignancy with a poor prognosis. One putative proto-oncogene in osteosarcoma is SKP2, encoding a substrate recognition factor of the SCF E3 ubiquitin ligase. We previously demonstrated that Skp2 knockout in murine osteosarcoma improved survival and delayed tumorigenesis. Here, we performed RNA sequencing (RNA-seq) on tumors from a transgenic osteosarcoma mouse model with conditional Trp53 and Rb1 knockouts in the osteoblast lineage ("DKO": Osx1-Cre;Rb1lox/lox;p53lox/lox) and a triple-knockout model with additional Skp2 germline knockout ("TKO": Osx1-Cre;Rb1lox/lox;p53lox/lox;Skp2-/-), followed by qPCR and immunohistochemistry validation. To investigate the clinical implications of our results, we analyzed a human osteosarcoma patient cohort ("NCI-TARGET OS") with RNA-seq and clinical data. We found large differences in gene expression after SKP2 knockout. Surprisingly, we observed increased expression of genes related to immune microenvironment infiltration in TKO tumors, especially the signature genes for macrophages and to a lesser extent, T cells, B cells, and vascular cells. We also uncovered a set of relevant transcription factors that may mediate these changes. In osteosarcoma patient cohorts, high expression of genes upregulated in TKO was correlated with favorable overall survival, which was largely explained by the macrophage gene signatures. This relationship was further supported by our finding that SKP2 expression was negatively correlated with macrophage infiltration in the NCI-TARGET osteosarcoma and the TCGA Sarcoma cohorts. Overall, our findings indicate that SKP2 may mediate immune exclusion from the osteosarcoma tumor microenvironment, suggesting that SKP2 modulation in osteosarcoma may induce antitumor immune activation.

Osteosarcoma Multi-Omics Landscape and Subtypes

Osteosarcoma (OS) is the most common primary bone malignancy that exhibits remarkable histologic diversity and genetic heterogeneity. The complex nature of osteosarcoma has confounded precise molecular categorization, prognosis, and prediction for this disease. In this study, we performed a comprehensive multiplatform analysis on 86 osteosarcoma tumors, including somatic copy-number alteration, gene expression and methylation, and identified three molecularly distinct and clinically relevant subtypes of osteosarcoma. The subgrouping criteria was validated on another cohort of osteosarcoma tumors. Previously unappreciated osteosarcoma-type-specific changes in specific genes' copy number, expression and methylation were revealed based on the subgrouping. The subgrouping and novel gene signatures provide insights into refining osteosarcoma therapy and relationships to other types of cancer.

Overview of CD70 as a Potential Therapeutic Target for Osteosarcoma

Osteosarcoma is a primary malignant bone tumor. Effective chemotherapy regimens for refractory disease are scarce, accounting for no improvement in survival. Immune-based cell therapies have emerged as novel alternatives. However, advancements with these therapies have been seen mostly when immune cells are armed to target specific tumor Ags. Recent studies identified cluster of differentiation 70 (CD70) as a promising target to osteosarcoma particularly because CD70 is highly expressed in osteosarcoma lung metastases (Pahl et al. 2015. Cancer Cell Int. 15: 31), and its overexpression by tumors has been correlated with immune evasion and tumor proliferation (Yang et al. 2007. Blood 110: 2537-2544). However, the limited knowledge of the overall CD70 expression within normal tissues and the potential for off-target effect pose several challenges (Flieswasser et al. 2022. J. Exp. Clin. Cancer Res. 41: 12). Nonetheless, CD70-based clinical trials are currently ongoing and are preliminarily showing promising results for patients with osteosarcoma. The present review sheds light on the recent literature on CD70 as it relates to osteosarcoma and highlights the benefits and challenges of targeting this pathway.

Large Scale Comparative Deconvolution Analysis of the Canine and Human Osteosarcoma Tumor Microenvironment Uncovers Conserved Clinically Relevant Subtypes

Osteosarcoma is a relatively rare but aggressive cancer of the bones with a shortage of effective biomarkers. Although less common in humans, Osteosarcomas are fairly common in adult pet dogs and have been shown to share many similarities with their human analogs. In this work, we analyze bulk transcriptomic data of 213 primary and 100 metastatic Osteosarcoma samples from 210 pet dogs enrolled in nation-wide clinical trials to uncover three Tumor Microenvironment (TME)-based subtypes: Immune Enriched (IE), Immune Enriched Dense Extra-Cellular Matrix-like (IE-ECM) and Immune Desert (ID) with distinct cell type compositions, oncogenic pathway activity and chromosomal instability. Furthermore, leveraging bulk transcriptomic data of canine primary tumors and their matched metastases from different sites, we characterize how the Osteosarcoma TME evolves from primary to metastatic disease in a standard of care clinical setting and assess its overall impact on clinical outcomes of canines. Most importantly, we find that TME-based subtypes of canine Osteosarcomas are conserved in humans and predictive of progression free survival outcomes of human patients, independently of known prognostic biomarkers such as presence of metastatic disease at diagnosis and percent necrosis following chemotherapy. In summary, these results demonstrate the power of using canines to model the human Osteosarcoma TME and discover novel biomarkers for clinical translation.

MYC regulates CSF1 expression via microRNA 17/20a to modulate tumor-associated macrophages in osteosarcoma

Osteosarcoma (OS) is the most common primary bone tumor of childhood. Approximately 20%-30% of OSs carry amplification of chromosome 8q24, which harbors the oncogene c-MYC and correlates with a poor prognosis. To understand the mechanisms that underlie the ability of MYC to alter both the tumor and its surrounding tumor microenvironment (TME), we generated and molecularly characterized an osteoblast-specific Cre-Lox-Stop-Lox-c-MycT58A p53fl/+ knockin genetically engineered mouse model (GEMM). Phenotypically, the Myc-knockin GEMM had rapid tumor development with a high incidence of metastasis. MYC-dependent gene signatures in our murine model demonstrated significant homology to the human hyperactivated MYC OS. We established that hyperactivation of MYC led to an immune-depleted TME in OS demonstrated by the reduced number of leukocytes, particularly macrophages. MYC hyperactivation led to the downregulation of macrophage colony-stimulating factor 1, through increased microRNA 17/20a expression, causing a reduction of macrophage population in the TME of OS. Furthermore, we developed cell lines from the GEMM tumors, including a degradation tag-MYC model system, which validated our MYC-dependent findings both in vitro and in vivo. Our studies utilized innovative and clinically relevant models to identify a potentially novel molecular mechanism through which MYC regulates the profile and function of the OS immune landscape.

The TLR3 L412F polymorphism prevents TLR3-mediated tumor cell death induction in pediatric sarcomas

Toll-like receptor 3 (TLR3) is a pattern recognition receptor mainly known for its role in innate immune response to infection. Indeed, binding of double-stranded RNA (dsRNA) to TLR3 triggers a pro-inflammatory cascade leading to cytokine release and immune cell activation. Its anti-tumoral potential has emerged progressively, associated with a direct impact on tumor cell death induction and with an indirect action on immune system reactivation. Accordingly, TLR3 agonists are currently being tested in clinical trials for several adult cancers. Meanwhile, TLR3 variants have been linked to auto-immune disorders, and as risk factors of viral infection and cancers. However, aside from neuroblastoma, TLR3 role in childhood cancers has not been evaluated. Here, by integrating public transcriptomic data of pediatric tumors, we unveil that high TLR3 expression is largely associated with a better prognosis in childhood sarcomas. Using osteosarcomas and rhabdomyosarcomas as models, we show that TLR3 efficiently drives tumor cell death in vitro and induces tumor regression in vivo. Interestingly, this anti-tumoral effect was lost in cells expressing the homozygous TLR3 L412F polymorphism, which is enriched in a rhabdomyosarcomas cohort. Thus, our results demonstrate the therapeutic potential associated with the targeting of TLR3 in pediatric sarcomas, but also the need to stratify patients eligible for this clinical approach with respect to the TLR3 variants expressed.

SKP2 knockout in Rb1/p53 deficient mouse models of osteosarcoma induces immune infiltration and drives a transcriptional program with a favorable prognosis

Purpose

Osteosarcoma (OS) is an aggressive bone malignancy with a poor prognosis. One putative proto-oncogene in OS is SKP2, encoding a substrate recognition factor of the SCF E3 ubiquitin ligase. We previously demonstrated that SKP2 knockout in murine OS improved survival and delayed tumorigenesis. Here we aim to define the SKP2 drives transcriptional program and its clinical implication in OS.

Experimental Design

We performed RNA-sequencing (RNA-seq) on tumors from a transgenic OS mouse model with conditional Trp53 and Rb1 knockouts in the osteoblast lineage ("DKO": Osx1-Cre;Rb1lox/lox;p53lox/lox) and a triple-knockout model with additional Skp2 germline knockout ("TKO": Osx1-Cre;Rb1lox/lox;p53lox/lox;SKP2-/-). We validated our RNA-seq findings using qPCR and immunohistochemistry. To investigate the clinical implications of our results, we analyzed a human OS patient cohort ("NCI-TARGET OS") with RNA-seq and clinical data.

Results

We found large differences in gene expression after SKP2 knockout. Strikingly, we observed increased expression of genes related to immune microenvironment infiltration in TKO tumors. We observed significant increases in signature genes for macrophages and to a lesser extent, T cells, B cells and vascular cells. We also uncovered a set of relevant transcription factors that may mediate the changes. In OS patient cohorts, high expression of genes upregulated in TKO was correlated with favorable overall survival, which was largely explained by the macrophage gene signatures. This relationship was further supported by our finding that SKP2 expression was negatively correlated with macrophage infiltration in the NCI-TARGET OS and the TCGA Sarcoma cohort.

Conclusion

Our findings indicate that SKP2 may mediate immune exclusion from the OS tumor microenvironment, suggesting that SKP2 modulation in OS may induce anti-tumor immune activation.

DNA methylation patterns suggest the involvement of DNMT3B and TET1 in osteosarcoma development

DNA methylation may be involved in the development of osteosarcomas. Osteosarcomas commonly arise during the bone growth and remodeling in puberty, making it plausible to infer the involvement of epigenetic alterations in their development. As a highly studied epigenetic mechanism, we investigated DNA methylation and related genetic variants in 28 primary osteosarcomas aiming to identify deregulated driver alterations. Methylation and genomic data were obtained using the Illumina HM450K beadchips and the TruSight One sequencing panel, respectively. Aberrant DNA methylation was spread throughout the osteosarcomas genomes. We identified 3146 differentially methylated CpGs comparing osteosarcomas and bone tissue samples, with high methylation heterogeneity, global hypomethylation and focal hypermethylation at CpG islands. Differentially methylated regions (DMR) were detected in 585 loci (319 hypomethylated and 266 hypermethylated), mapped to the promoter regions of 350 genes. These DMR genes were enriched for biological processes related to skeletal system morphogenesis, proliferation, inflammatory response, and signal transduction. Both methylation and expression data were validated in independent groups of cases. Six tumor suppressor genes harbored deletions or promoter hypermethylation (DLEC1, GJB2, HIC1, MIR149, PAX6, and WNT5A), and four oncogenes presented gains or hypomethylation (ASPSCR1, NOTCH4, PRDM16, and RUNX3). Our analysis also revealed hypomethylation at 6p22, a region that contains several histone genes. Copy-number changes in DNMT3B (gain) and TET1 (loss), as well as overexpression of DNMT3B in osteosarcomas provide a possible explanation for the observed phenotype of CpG island hypermethylation. While the detected open-sea hypomethylation likely contributes to the well-known osteosarcoma genomic instability, enriched CpG island hypermethylation suggests an underlying mechanism possibly driven by overexpression of DNMT3B likely resulting in silencing of tumor suppressors and DNA repair genes.

Nuclear receptor modulators inhibit osteosarcoma cell proliferation and tumour growth by regulating the mTOR signaling pathway

Osteosarcoma is the most common primary malignant bone tumour in children and adolescents. Chemoresistance leads to poor responses to conventional therapy in patients with osteosarcoma. The discovery of novel effective therapeutic targets and drugs is still the main focus of osteosarcoma research. Nuclear receptors (NRs) have shown substantial promise as novel therapeutic targets for various cancers. In the present study, we performed a drug screen using 29 chemicals that specifically target 17 NRs in several different human osteosarcoma and osteoblast cell lines. The retinoic acid receptor beta (RARb) antagonist LE135, peroxisome proliferator activated receptor gamma (PPARg) antagonist T0070907, liver X receptor (LXR) agonist T0901317 and Rev-Erba agonist SR9011 significantly inhibited the proliferation of malignant osteosarcoma cells (U2OS, HOS-MNNG and Saos-2 cells) but did not inhibit the growth of normal osteoblasts. The effects of these NR modulators on osteosarcoma cells occurred in a dose-dependent manner and were not observed in NR-knockout osteosarcoma cells. These NR modulators also significantly inhibited osteosarcoma growth in vivo and enhanced the antitumour effect of doxorubicin (DOX). Transcriptomic and immunoblotting results showed that these NR modulators may inhibit the growth of osteosarcoma cells by regulating the PI3K/AKT/mTOR and ERK/mTOR pathways. DDIT4, which blocks mTOR activation, was identified as one of the common downstream target genes of these NRs. DDIT4 knockout significantly attenuated the inhibitory effects of these NR modulators on osteosarcoma cell growth. Together, our results revealed that modulators of RARb, PPARg, LXRs and Rev-Erba inhibit osteosarcoma growth both in vitro and in vivo through the mTOR signaling pathway, suggesting that treatment with these NR modulators is a novel potential therapeutic strategy.

Pan-Cancer analysis of TERT and Validation in Osteosarcoma Cell Lines

PURPOSE

To investigate the possible functions of TERT in pan-cancer and OS.

METHODS

First, differential TERT gene expression analysis was conducted using multi-omics data integrative analyses, including differential expression, prognosis, the correlation between infiltrating inflammatory immune cells, and mutation in pan-cancer. Furthermore, differential TERT gene expression analysis was conducted using mRNA expression profiles related to OS based on the GEO Datasets. Various differentially expressed genes were chosen based on a fitness threshold for further investigations. Finally, the function of the TERT gene was assessed in OS cells, including cellular proliferation, migration, and metastasis.

RESULTS

Pan-cancer research demonstrated that variable expression of TERT was not only associated with numerous types of human cancer but was also intimately linked to DNA methylation. Bioinformatic investigation revealed a link between the differential expression of TERT with immune cell infiltration in the tumor microenvironment (TME). In vitro studies indicated that inhibition of TERT decreased OS cell proliferation, motility, and metastasis.

CONCLUSION

TERT may serve as a useful genomic biomarker for the diagnosis and prediction of pan-cancer and as a prospective therapeutic target for the treatment of OS.

RNA-seq and Single-Cell Transcriptome Analyses of TRAIL Receptors Gene Expression in Human Osteosarcoma Cells and Tissues

Osteosarcoma (OS) is the most common primary cancer in the skeletal system, characterized by a high incidence of lung metastasis, local recurrence and death. Systemic treatment of this aggressive cancer has not improved significantly since the introduction of chemotherapy regimens, underscoring a critical need for new treatment strategies. TRAIL receptors have long been proposed to be therapeutic targets for cancer treatment, but their role in osteosarcoma remains unclear. In this study, we investigated the expression profile of four TRAIL receptors in human OS cells using total RNA-seq and single-cell RNA-seq (scRNA-seq). The results revealed that TNFRSF10B and TNFRSF10D but not TNFRSF10A and TNFRSF10C are differentially expressed in human OS cells as compared to normal cells. At the single cell level by scRNA-seq analyses, TNFRSF10B, TNFRSF10D, TNFRSF10A and TNFRSF10C are most abundantly expressed in endothelial cells of OS tissues among nine distinct cell clusters. Notably, in osteoblastic OS cells, TNFRSF10B is most abundantly expressed, followed by TNFRSF10D, TNFRSF10A and TNFRSF10C. Similarly, in an OS cell line U2-OS using RNA-seq, TNFRSF10B is most abundantly expressed, followed by TNFRSF10D, TNFRSF10A and TNFRSF10C. According to the TARGET online database, poor patient outcomes were associated with low expression of TNFRSF10C. These results could provide a new perspective to design novel therapeutic targets of TRAIL receptors for the diagnosis, prognosis and treatment of OS and other cancers.

CCT6A knockdown suppresses osteosarcoma cell growth and Akt pathway activation in vitro

We assessed the role of the protein-coding gene chaperonin-containing TCP1 subunit 6A (CCT6A) in osteosarcoma, as this is currently unknown. Using data from the R2 online genomic analysis and visualization application, we found that CCT6A messenger ribonucleic acid (RNA) expression is increased in osteosarcoma tissue and cells. Transfection of CCT6A small interfering RNA into cultured osteosarcoma cells revealed that CCT6A knockdown attenuates cell growth, cell viability, cell survival, and induced apoptosis and cell cycle progression at the G0/G1 phases. Moreover, CCT6A knockdown downregulated phospho-protein kinase B (p-Akt), cyclinD1 and B-cell lymphoma-2, whereas upregulated Bcl-2-associated X-protein expression. Thus, CCT6A knockdown inhibits cell proliferation, induces cell apoptosis, and suppresses the Akt pathway.

SLC7A8 coding for LAT2 is associated with early disease progression in osteosarcoma and transports doxorubicin

Background: Despite (neo) adjuvant chemotherapy with cisplatin, doxorubicin and methotrexate, some patients with primary osteosarcoma progress during first-line systemic treatment and have a poor prognosis. In this study, we investigated whether patients with early disease progression (EDP), are characterized by a distinctive pharmacogenetic profile.

Methods and Findings: Germline DNA from 287 Dutch high-grade osteosarcoma patients was genotyped using the DMET Plus array (containing 1,936 genetic markers in 231 drug metabolism and transporter genes). Associations between genetic variants and EDP were assessed using logistic regression models and associated variants (p <0.05) were validated in independent cohorts of 146 (Spain and United Kingdom) and 28 patients (Australia). In the association analyses, EDP was significantly associated with an SLC7A8 locus and was independently validated (meta-analysis validation cohorts: OR 0.19 [0.06–0.55], p = 0.002). The functional relevance of the top hits was explored by immunohistochemistry staining and an in vitro transport models. SLC7A8 encodes for the L-type amino acid transporter 2 (LAT2). Transport assays in HEK293 cells overexpressing LAT2 showed that doxorubicin, but not cisplatin and methotrexate, is a substrate for LAT2 (p < 0.0001). Finally, SLC7A8 mRNA expression analysis and LAT2 immunohistochemistry of osteosarcoma tissue showed that the lack of LAT2 expression is a prognostic factor of poor prognosis and reduced overall survival in patients without metastases (p = 0.0099 and p = 0.14, resp.).

Conclusion: This study identified a novel locus in SLC7A8 to be associated with EDP in osteosarcoma. Functional studies indicate LAT2-mediates uptake of doxorubicin, which could give new opportunities to personalize treatment of osteosarcoma patients.

UHRF1 overexpression promotes osteosarcoma metastasis through altered exosome production and AMPK/SEMA3E suppression

Loss-of-function mutations at the retinoblastoma (RB1) gene are associated with increased mortality, metastasis, and poor therapeutic outcome in several cancers, including osteosarcoma. However, the mechanism(s) through which RB1 loss worsens clinical outcome remains understudied. Ubiquitin-like with PHD and Ring Finger domains 1 (UHRF1) has been identified as a critical downstream effector of the RB/E2F signaling pathway that is overexpressed in various cancers. Here, we determined the role and regulatory mechanisms of UHRF1 in rendering osteosarcoma cells more aggressive. Higher UHRF1 expression correlated with malignancy in osteosarcoma cell lines, clinical samples, and genetically engineered mouse models. Gain- and loss-of-function assays revealed that UHRF1 has cell-intrinsic and extrinsic functions promoting cell proliferation, migration, invasion, angiogenesis, and metastasis. UHRF1 overexpression induced angiogenesis by suppressing AMPK activation and Semaphorin 3E (SEMA3E) expression. Further, UHRF1-mediated migration and metastasis resulted, at least in part, through altered expression of extracellular vesicles and their cargo, including urokinase-type plasminogen activator (uPA). Novel osteosarcoma genetically engineered mouse models confirmed that knocking out Uhrf1 considerably decreased metastasis and reversed the poorer survival associated with Rb1 loss. This presents a new mechanistic insight into RB1 loss-associated poor prognosis and novel oncogenic roles of UHRF1 in the regulation of angiogenesis and exosome secretion, both critical for osteosarcoma metastasis. This provides substantial support for targeting UHRF1 or its downstream effectors as novel therapeutic options to improve current treatment for osteosarcoma.

Integrated DNA Copy Number and Expression Profiling Identifies IGF1R as a Prognostic Biomarker in Pediatric Osteosarcoma

Osteosarcoma is a primary malignant bone tumor arising from bone-forming mesenchymal cells in children and adolescents. Despite efforts to understand the biology of the disease and identify novel therapeutics, the survival of osteosarcoma patients remains dismal. We have concurrently profiled the copy number and gene expression of 226 osteosarcoma samples as part of the Strategic Partnering to Evaluate Cancer Signatures (SPECS) initiative. Our results demonstrate the heterogeneous landscape of osteosarcoma in younger populations by showing the presence of genome-wide copy number abnormalities occurring both recurrently among samples and in a high frequency. Insulin growth factor receptor 1 (IGF1R) is a receptor tyrosine kinase which binds IGF1 and IGF2 to activate downstream pathways involved in cell apoptosis and proliferation. We identify prevalent amplification of IGF1R corresponding with increased gene expression in patients with poor survival outcomes. Our results substantiate previously tenuously associated copy number abnormalities identified in smaller datasets (13q34+, 20p13+, 4q35-, 20q13.33-), and indicate the significance of high fibroblast growth factor receptor 2 (FGFR2) expression in distinguishing patients with poor prognosis. FGFR2 is involved in cellular proliferation processes such as division, growth and angiogenesis. In summary, our findings demonstrate the prognostic significance of several genes associated with osteosarcoma pathogenesis.

New dawn for cancer cell death: Emerging role of lipid metabolism

BACKGROUND

Resistance to cell death, a protective mechanism for removing damaged cells, is a "Hallmark of Cancer" that is essential for cancer progression. Increasing attention to cancer lipid metabolism has revealed a number of pathways that induce cancer cell death.

SCOPE OF REVIEW

We summarize emerging concepts regarding lipid metabolic reprogramming in cancer that is mainly involved in lipid uptake and trafficking, de novo synthesis and esterification, fatty acid synthesis and oxidation, lipogenesis, and lipolysis. During carcinogenesis and progression, continuous metabolic adaptations are co-opted by cancer cells, to maximize their fitness to the ever-changing environmental. Lipid metabolism and the epigenetic modifying enzymes interact in a bidirectional manner which involves regulating cancer cell death. Moreover, lipids in the tumor microenvironment play unique roles beyond metabolic requirements that promote cancer progression. Finally, we posit potential therapeutic strategies targeting lipid metabolism to improve treatment efficacy and survival of cancer patient.

MAJOR CONCLUSIONS

The profound comprehension of past findings, current trends, and future research directions on resistance to cancer cell death will facilitate the development of novel therapeutic strategies targeting the lipid metabolism.

Protein Expression of AEBP1, MCM4, and FABP4 Differentiate Osteogenic, Adipogenic, and Mesenchymal Stromal Stem Cells

Mesenchymal stem cells (MSCs) gain an increasing focus in the field of regenerative medicine due to their differentiation abilities into chondrocytes, adipocytes, and osteoblastic cells. However, it is apparent that the transformation processes are extremely complex and cause cellular heterogeneity. The study aimed to characterize differences between MSCs and cells after adipogenic (AD) or osteoblastic (OB) differentiation at the proteome level. Comparative proteomic profiling was performed using tandem mass spectrometry in data-independent acquisition mode. Proteins were quantified by deep neural networks in library-free mode and correlated to the Molecular Signature Database (MSigDB) hallmark gene set collections for functional annotation. We analyzed 4108 proteins across all samples, which revealed a distinct clustering between MSCs and cell differentiation states. Protein expression profiling identified activation of the Peroxisome proliferator-activated receptors (PPARs) signaling pathway after AD. In addition, two distinct protein marker panels could be defined for osteoblastic and adipocytic cell lineages. Hereby, overexpression of AEBP1 and MCM4 for OB as well as of FABP4 for AD was detected as the most promising molecular markers. Combination of deep neural network and machine-learning algorithms with data-independent mass spectrometry distinguish MSCs and cell lineages after adipogenic or osteoblastic differentiation. We identified specific proteins as the molecular basis for bone formation, which could be used for regenerative medicine in the future.

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.

Hsp90 regulates the tumorigenic function of tyrosine protein kinase in osteosarcoma

Despite recent advances in diagnosis and treatment, osteosarcoma remains as the most common bone cancer in children and is associated with poor prognosis. Growing evidence has supported dysregulation of threonine and tyrosine protein kinase (TTK) expression as a hallmark of multiple cancers, however, its function in osteosarcoma remains to be elucidated. In the present study, we found that TTK was frequently overexpressed in osteosarcoma and associated with increased tumour growth and progression. Moreover, using both in vitro and in vivo assays, we provided evidence that TTK level was regulated by a molecular chaperone, heat shock protein 90 (Hsp90). Hsp90 directly interacted with TTK and prevents proteasome-dependent TTK degradation, leading to the accumulation of TTK in osteosarcoma cells. Elevated TTK promoted cancer cell proliferation and survival by activating cell-cycle progression and inhibiting apoptosis. Consistently, depletion of TTK by Hsp90 inhibition induced cell-cycle arrest, generated aneuploidy and eventually resulted in apoptotic cancer cell death. Together, our study revealed an important Hsp90-TTK regulatory axis in osteosarcoma cells to promote cancer cell growth and survival. These findings expand our knowledge on osteosarcoma pathogenesis and offer novel therapeutic options for clinical practice.

A Zebrafish Model of Metastatic Colonization Pinpoints Cellular Mechanisms of Circulating Tumor Cell Extravasation

Metastasis is a multistep process in which cells must detach, migrate/invade local structures, intravasate, circulate, extravasate, and colonize. A full understanding of the complexity of this process has been limited by the lack of ability to study these steps in isolation with detailed molecular analyses. Leveraging a comparative oncology approach, we injected canine osteosarcoma cells into the circulation of transgenic zebrafish with fluorescent blood vessels in a biologically dynamic metastasis extravasation model. Circulating tumor cell clusters that successfully extravasated the vasculature as multicellular units were isolated under intravital imaging (n = 6). These extravasation-positive tumor cell clusters sublines were then molecularly profiled by RNA-Seq. Using a systems-level analysis, we pinpointed the downregulation of KRAS signaling, immune pathways, and extracellular matrix (ECM) organization as enriched in extravasated cells (p < 0.05). Within the extracellular matrix remodeling pathway, we identified versican (VCAN) as consistently upregulated and central to the ECM gene regulatory network (p < 0.05). Versican expression is prognostic for a poorer metastasis-free and overall survival in patients with osteosarcoma. Together, our results provide a novel experimental framework to study discrete steps in the metastatic process. Using this system, we identify the versican/ECM network dysregulation as a potential contributor to osteosarcoma circulating tumor cell metastasis.

Identification of PSMD14 as a potential novel prognosis biomarker and therapeutic target for osteosarcoma

BACKGROUND

Osteosarcoma is the most common primary bone tumor. The survival rate of osteosarcoma patients has not significantly increased in the past decades. Uncovering the mechanisms of malignancy, progression, and metastasis will shed light on the development of new therapeutic targets and treatment for osteosarcoma.

AIM

The aim of this study is to identify potential osteosarcoma biomarker and/or therapeutic targets by using integrated bioinformatics analysis.

METHODS AND RESULTS

We utilized existing gene expression datasets to identify differential expressed genes (DEGs) that could serve as osteosarcoma biomarkers or even as therapeutic targets. We found 48 DEGs were overlapped in three datasets. Among these 48 DEGs, PSMD14 was on the top of the up-regulated gene list. We further found that higher PSMD14 expression was correlated with higher risk group (younger age group, ≤20.83 years of age), metastasis within 5 years and higher grade of tumor. Higher PSMD14 expression in osteosarcoma had positive correlation with higher infiltration of CD8+ T cells, neutrophils and myeloid dendritic cells. Kaplan-Myer survival data further revealed that higher expression of PSMD14 predicted significantly worse prognosis (p = .013). Gene set enrichment analysis was further performed for the DEGs related to PSMD14 in osteosarcoma. We found that lower PSMD14 expression group had more immune responses such as interferon γ, α responses, inflammation response etc. However, the higher PSMD14 expression group had more cell proliferation-related biological processes, such as G2M checkpoints and Myc targets. Through establishing protein-protein interaction networks using PSMD14 related DEGs, we identified 10 hub genes that were all ribosomal proteins. These hub genes may play roles in osteosarcoma tumorigenesis, progression and/or metastasis.

CONCLUSION

We identified PSMD14 gene as a possible osteosarcoma biomarker, and/or a possible therapeutic target.

The Immune Landscape of Osteosarcoma: Implications for Prognosis and Treatment Response

Osteosarcoma (OS) is a high-grade malignant stromal tumor composed of mesenchymal cells producing osteoid and immature bone, with a peak of incidence in the second decade of life. Hence, although relatively rare, the social impact of this neoplasm is particularly relevant. Differently from carcinomas, molecular genetics and the role of the tumor microenvironment in the development and progression of OS are mainly unknown. Indeed, while the tumor microenvironment has been widely studied in other solid tumor types and its contribution to tumor progression has been definitely established, tumor-stroma interaction in OS has been quite neglected for years. Only recently have new insights been gained, also thanks to the availability of new technologies and bioinformatics tools. A better understanding of the cross-talk between the bone microenvironment, including immune and stromal cells, and OS will be key not only for a deeper knowledge of osteosarcoma pathophysiology, but also for the development of novel therapeutic strategies. In this review, we summarize the current knowledge about the tumor microenvironment in OS, mainly focusing on immune cells, discussing their role and implication for disease prognosis and treatment response.

PLEK2 promotes osteosarcoma tumorigenesis and metastasis by activating the PI3K/AKT signaling pathway

Increasing evidence suggest that pleckstrin-2 (PLEK2) acts as an oncogene in several malignancies. The present study aimed to investigate the effects of PLEK2 on osteosarcoma (OS) tumorigenesis and metastasis. PLEK2 expression in OS was analyzed via bioinformatics, reverse transcription-quantitative PCR, western blot and immunohistochemistry analyses. The Cell Counting Kit-8 (CCK-8), colony formation and EdU assays were performed to assess the role of PLEK2 in OS cell proliferation. The pro-metastatic effects of PLEK2 were assessed via the Transwell and wound healing assays. In addition, the PLEK2 downstream pathway was analyzed via bioinformatics analysis and verified via western blot analysis. The results demonstrated that PLEK2 expression was upregulated in both OS cell lines and specimens. The results of the CCK-8, colony formation and EdU assays demonstrated that PLEK2 promoted OS cell proliferation in vitro. The in vivo experiments further demonstrated that PLEK2 knockdown significantly suppressed OS growth. In addition, the Transwell and wound healing assays indicated that PLEK2 promoted OS invasiveness in vitro, which was induced by the activation of the epithelial-to-mesenchymal transition process. Bioinformatics analysis revealed that PLEK2 can activate the phosphoinositide 3-kinase (PI3K)/protein kinase B (AKT)/mechanistic target of rapamycin (mTOR) pathway, which was verified via western blot analysis. Taken together, the results of the present study suggest that PLEK2 may play a tumor-promoting role in OS via the PI3K/AKT signaling pathway.

The Clinical Implications of Tumor Mutational Burden in Osteosarcoma

Background

Osteosarcoma (OTS) is aggressive bone malignancy without well-recognized prognosis biomarker. Tumor mutational burden (TMB) has been proved as effective biomarker in predicting clinical outcomes in several cancer types. However, its prognostic value in OTS remains unknown. In this study, we aim to evaluate the implication of TMB in OTS patients.

Methods

To depict the landscape of somatic mutations in OTS, we performed Whole-Exome Sequencing (WES) on 31 OTS tissue samples and corresponding White Blood Cells (WBCs) as matched control. TMB was calculated as the total number of somatic alterations in coding regions normalized to the per sequenced genomic megabase (~30.4Mb in WES). The prognostic values of TMB were evaluated by Kaplan-Meier methods and Cox regression models.

Results

The median age was 16.0 years at diagnosis, and 54.8% of patients were male. The most common genetic alterations were mainly involved in cell cycle and DNA damage response and repair, including H3F3A, TP53, MYC, and CDKN2A/B. The median progression-free survival (PFS) was 775.5 days in TMB-High (defined as third quartile of TMB value, <2.565) versus 351 days in TMB-Low (<2.565). All patients with TMB-High are PFS-Long (>400 days), while 36.4% of all patients with TMB-Low were PFS-Long (P=0.003). TMB were significantly greater in PFS-Long than in PFS-Short (<400 days) (P=0.002). Moreover, the median overall survival (OS) was 1,307 days in TMB-High versus 672.5 days in TMB-Low. Furthermore, TMB-High group had significantly improved PFS (P=0.04) and OS (P=0.03).

Conclusions

TMB-High can be used as prognostic marker for OTS. Our findings demonstrate that TMB may be helpful in combination with traditionally clinicopathologic risk factors to optimize risk stratification and guide treatment decisions.

BMPR2 expression level is correlated with low immune infiltration and predicts metastasis and poor survival in osteosarcoma

Osteosarcoma is the most common malignant bone tumor in adolescents and young adults, and identifying biomarkers for prognosis and therapy is necessary. Bone morphogenetic protein receptor 2 (BMPR2) is involved in various cellular functions, including cell adhesion, proliferation and invasion, inflammation, apoptosis and metastatic spread. However, the correlation between BMPR2 expression levels and prognosis and tumor-infiltrating immune cells in osteosarcoma is not well understood. In the present study, the expression level of BMPR2 was investigated using the Oncomine and R2 databases. The association between the expression level of BMPR2 and the clinical prognosis of patients with cancer was analyzed using the R2 database. The relationship between the expression level of BMPR2 and immune cell infiltration in the stroma of osteosarcoma was assessed using the Tumor Immune Estimation Resource (TIMER) and CIBERSORT. The correlations between BMPR2 expression level and infiltrated immune cell gene marker sets in osteosarcoma were validated in the TIMER and R2 databases. Analysis of a cohort of patients with osteosarcoma revealed that BMPR2 expression was significantly higher in osteosarcoma compared with in normal tissue and was correlated with poor prognosis. M0 macrophages, M2 macrophages, resting mast, γ δ T and CD8+ T cells were the top five immune cells with the highest degrees of infiltration in osteosarcoma. In addition, BMPR2 expression level showed a significant negative correlation with the gene markers of CD8+ T cells, monocytes and M2 macrophages. Low levels of infiltrating CD8+ T cells, monocytes or M2 macrophages in osteosarcoma was significantly associated with poor survival. These data suggested that CD8+ T cells, monocytes and M2 macrophages play significant roles in the establishment of the immune microenvironment of osteosarcoma. High BMPR2 expression was associated with poor prognosis and low infiltration of CD8+ T cells, monocytes and M2 macrophages in osteosarcoma. Hence, BMPR2 can be considered a biomarker of the immune infiltration, metastasis and prognosis of osteosarcoma.

Systems Biology Approach Identifies Prognostic Signatures of Poor Overall Survival and Guides the Prioritization of Novel BET-CHK1 Combination Therapy for Osteosarcoma

Osteosarcoma (OS) patients exhibit poor overall survival, partly due to copy number variations (CNVs) resulting in dysregulated gene expression and therapeutic resistance. To identify actionable prognostic signatures of poor overall survival, we employed a systems biology approach using public databases to integrate CNVs, gene expression, and survival outcomes in pediatric, adolescent, and young adult OS patients. Chromosome 8 was a hotspot for poor prognostic signatures. The MYC-RAD21 copy number gain (8q24) correlated with increased gene expression and poor overall survival in 90% of the patients (n = 85). MYC and RAD21 play a role in replication-stress, which is a therapeutically actionable network. We prioritized replication-stress regulators, bromodomain and extra-terminal proteins (BETs), and CHK1, in order to test the hypothesis that the inhibition of BET + CHK1 in MYC-RAD21+ pediatric OS models would be efficacious and safe. We demonstrate that MYC-RAD21+ pediatric OS cell lines were sensitive to the inhibition of BET (BETi) and CHK1 (CHK1i) at clinically achievable concentrations. While the potentiation of CHK1i-mediated effects by BETi was BET-BRD4-dependent, MYC expression was BET-BRD4-independent. In MYC-RAD21+ pediatric OS xenografts, BETi + CHK1i significantly decreased tumor growth, increased survival, and was well tolerated. Therefore, targeting replication stress is a promising strategy to pursue as a therapeutic option for this devastating disease.

Identification of potential crucial genes and key pathways in osteosarcoma

Background

The aim of this study is to identify the potential pathogenic and metastasis-related differentially expressed genes (DEGs) in osteosarcoma through bioinformatic analysis of Gene Expression Omnibus (GEO) database.

Results

Gene expression profiles of GSE14359, GSE16088, and GSE33383, in total 112 osteosarcoma tissue samples and 7 osteoblasts, were analyzed. Seventy-four normal-primary DEGs (NPDEGs) and 764 primary-metastatic DEGs (PMDEGs) were screened. VAMP8, A2M, HLA-DRA, SPARCL1, HLA-DQA1, APOC1 and AQP1 were identified continuously upregulating during the oncogenesis and metastasis of osteosarcoma. The enriched functions and pathways of NPDEGs include procession and presentation of antigens, activation of MHC class II receptors and phagocytosis. The enriched functions and pathways of PMDEGs include mitotic nuclear division, cell adhesion molecules (CAMs) and focal adhesion. With protein-protein interaction (PPI) network analyzed by Molecular Complex Detection (MCODE) plug-in of Cytoscape software, one hub NPDEG (HLA-DRA) and 7 hub PMDEGs (CDK1, CDK20, CCNB1, MTIF2, MRPS7, VEGFA and EGF) were eventually selected, and the most significant pathways in NPDEGs module and PMDEGs module were enriched in the procession and presentation of exogenous peptide antigen via MHC class II and the nuclear division, respectively.

Conclusions

By integrated bioinformatic analysis, numerous DEGs related to osteosarcoma were screened, and the hub DEGs identified in this study are possibly part of the potential biomarkers for osteosarcoma. However, further experimental studies are still necessary to elucidate the biological function and mechanism of these genes.

Experimental Study of Somatic Variants of Osteosarcoma by Whole-Exome Sequencing

Background

This study aimed to investigate the role of gene mutation site distribution, biological function, pathway enrichment, and gene association analysis in the occurrence, development, and migration of osteosarcoma.

Material/Methods

Somatic mutation screening was performed using the whole-exome sequencing of osteosarcoma samples, and the distribution of mutations was demonstrated by Circos diagrams. Metascape was used to analyze the GO and KEGG signal pathway enrichment of the genes harboring protein coding alterations, and GeneMANIA was used to analyze the interaction of mutated genes.

Results

The results showed that the protein coding alterations were found throughout the whole genome in 3 osteosarcoma samples. A large number of identical or related biological processes and pathways were found in osteosarcoma samples. The GeneMANIA analysis of the 10 mutations shared by 3 samples showed that the target gene minichromosome maintenance complex component 4 (MCM4) and 3 lateral genes were most functional, and were all related to DNA replication. The analysis of GO and KEGG signal pathway enrichment showed that the mutated genes were involved mainly in tumor-related metabolic pathways. Three mutated genes were involved in the cell process, and 2 mutated genes were involved in the metabolic process. Known driver gene mutations were also observed in the samples.

Conclusions

The gene analysis confirmed that patients with osteosarcoma had a wide range of common gene mutations related to each other, which are involved in tumor-related metabolic pathways. These findings provide a basis for further gene-targeted therapy and pathway research.

PRMT5-TRIM21 interaction regulates the senescence of osteosarcoma cells by targeting the TXNIP/p21 axis

Osteosarcoma (OS) is the most common bone malignancy in adolescents and has poor clinical outcomes. Protein arginine methyltransferase 5 (PRMT5) has recently been shown to be aberrantly expressed in various cancers, yet its role in OS remains elusive. Here, we found that PRMT5 was overexpressed in OS and its overexpression predicted poor clinical outcomes. PRMT5 knockdown significantly triggered pronounced senescence in OS cells, as evidenced by the increase in senescence-associated β-galactosidase (SA-β-gal)-stained cells, induction of p21 expression, and upregulation of senescence-associated secretory phenotype (SASP) gene expression. In addition, we found that PRMT5 plays a key role in regulating DNA damaging agents-induced OS cell senescence, possibly, via affecting the repair of DNA damage. Furthermore, we found that TXNIP acts as a key factor mediating PRMT5 depletion-induced DNA damage and cellular senescence. Mechanistically, TRIM21, which interacts with PRMT5, was essential for the regulation of TXNIP/p21 expression. In summary, we propose a model in which PRMT5, by interaction with TRIM21, plays a key role in regulating the TXNIP/p21 axis during senescence in OS cells. The present findings suggest that PRMT5 overexpression in OS cells might confer resistance to chemotherapy and that targeting the PRMT5/TRIM21/TXNIP signaling may enhance the therapeutic efficacy in OS.

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.

Integrative genomic analysis of matched primary and metastatic pediatric osteosarcoma

Despite being the most common childhood bone tumor, the genomic characterization of osteosarcoma remains incomplete. In particular, very few osteosarcoma metastases have been sequenced to date, critical to better understand mechanisms of progression and evolution in this tumor. We performed an integrated whole genome and exome sequencing analysis of paired primary and metastatic pediatric osteosarcoma specimens to identify recurrent genomic alterations. Sequencing of 13 osteosarcoma patients including 13 primary, 10 metastatic, and 3 locally recurring tumors revealed a highly heterogeneous mutational landscape, including cases of hypermutation and microsatellite instability positivity, but with virtually no recurrent alterations except for mutations involving the tumor suppressor genes RB1 and TP53. At the germline level, we detected alterations in multiple cancer related genes in the majority of the cohort, including those potentially disrupting DNA damage response pathways. Metastases retained only a minimal number of short variants from their corresponding primary tumors, while copy number alterations showed higher conservation. One recurrently amplified gene, KDR, was highly expressed in advanced cases and associated with poor prognosis. © 2019 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

A meta-analysis of multiple matched aCGH/expression cancer datasets reveals regulatory relationships and pathway enrichment of potential oncogenes

The copy numbers of genes in cancer samples are often highly disrupted and form a natural amplification/deletion experiment encompassing multiple genes. Matched array comparative genomics and transcriptomics datasets from such samples can be used to predict inter-chromosomal gene regulatory relationships. Previously we published the database METAMATCHED, comprising the results from such an analysis of a large number of publically available cancer datasets. Here we investigate genes in the database which are unusual in that their copy number exhibits consistent heterogeneous disruption in a high proportion of the cancer datasets. We assess the potential relevance of these genes to the pathology of the cancer samples, in light of their predicted regulatory relationships and enriched biological pathways. A network-based method was used to identify enriched pathways from the genes’ inferred targets. The analysis predicts both known and new regulator-target interactions and pathway memberships. We examine examples in detail, in particular the gene POGZ, which is disrupted in many of the cancer datasets and has an unusually large number of predicted targets, from which the network analysis predicts membership of cancer related pathways. The results suggest close involvement in known cancer pathways of genes exhibiting consistent heterogeneous copy number disruption. Further experimental work would clarify their relevance to tumor biology. The results of the analysis presented in the database METAMATCHED, and included here as an R archive file, constitute a large number of predicted regulatory relationships and pathway memberships which we anticipate will be useful in informing such experiments.

Targeted transcriptional profiling of the tumor microenvironment reveals lymphocyte exclusion and vascular dysfunction in metastatic osteosarcoma

Osteosarcoma (OS) is the most common bone tumor in pediatric and adolescent/young adult patients yet little is known about the microenvironment that supports this aggressive disease. We have used targeted gene expression profiling and immunohistochemistry to characterize the microenvironment of metastatic and non-metastatic OS specimens from pediatric patients exhibiting poor histologic response to chemotherapy. Our results indicate that metastatic specimens exhibit lymphocyte exclusion as T cells are confined to the periphery of the pulmonary lesions. Furthermore, our data provides evidence of vascular dysfunction in metastatic OS indicated by increased expression of VEGFA, an increased ANGPT2:ANGPT1 gene expression ratio, and decreased expression of SELE, the gene encoding the adhesion molecule E-selectin. Moreover, correlation analyses show an inverse relationship between lymphocyte abundance and markers of vascular dysfunction exclusively in the metastatic specimens. Together, our data shows that the non-metastatic OS specimens demonstrate increased expression of various immunotherapeutic targets in comparison metastatic specimens and identifies vascular dysfunction and lymphocyte exclusion as important processes for therapeutic intervention in metastatic disease.

Regulation of TEAD Transcription Factors in Cancer Biology

Transcriptional enhanced associate domain (TEAD) transcription factors play important roles during development, cell proliferation, regeneration, and tissue homeostasis. TEAD integrates with and coordinates various signal transduction pathways including Hippo, Wnt, transforming growth factor beta (TGFβ), and epidermal growth factor receptor (EGFR) pathways. TEAD deregulation affects well-established cancer genes such as KRAS, BRAF, LKB1, NF2, and MYC, and its transcriptional output plays an important role in tumor progression, metastasis, cancer metabolism, immunity, and drug resistance. To date, TEADs have been recognized to be key transcription factors of the Hippo pathway. Therefore, most studies are focused on the Hippo kinases and YAP/TAZ, whereas the Hippo-dependent and Hippo-independent regulators and regulations governing TEAD only emerged recently. Deregulation of the TEAD transcriptional output plays important roles in tumor progression and serves as a prognostic biomarker due to high correlation with clinicopathological parameters in human malignancies. In addition, discovering the molecular mechanisms of TEAD, such as post-translational modifications and nucleocytoplasmic shuttling, represents an important means of modulating TEAD transcriptional activity. Collectively, this review highlights the role of TEAD in multistep-tumorigenesis by interacting with upstream oncogenic signaling pathways and controlling downstream target genes, which provides unprecedented insight and rationale into developing TEAD-targeted anticancer therapeutics.

Co-Expression Network Analysis Identified Gene Signatures in Osteosarcoma as a Predictive Tool for Lung Metastasis and Survival

Osteosarcoma (OS) is the most common primary bone tumor, whose poor prognosis is mainly due to lung metastasis. The aim of this study is to build a practical and valid diagnostic test that can predict the risk of OS metastasis and progression. We performed weighted gene co-expression network analysis (WGCNA) on GSE21257 from the Gene Expression Omnibus (GEO) database, which contains microarray data of biopsies from OS patients. In these modules, the highest association was found between the blue module and metastasis stage (r = -0.52) by Pearson's correlation analysis. Based on Least Absolute Shrinkage and Selection Operator (LASSO) Cox regression, we derived eight clinically significant genes and constructed an eight-gene signature for metastasis status. It showed great efficacy to distinguish metastasis from non-metastasis (AUC = 0.886) and the results were validated in The Cancer Genome Atlas (TCGA) database. Functional enrichment analysis of hub genes showed that their biological processes focused on immune-related pathways, suggesting the important roles of immune cells, immune pathways and the tumor microenvironment in metastasis development. In conclusion, we discovered an efficient gene signature with great efficacy to distinguish metastasis status, which may help improve early diagnosis and treatment, enhancing the clinical outcomes of OS patients. Besides we created an effective protocol to seek for several hub genes in high-throughput data by combining WGCNA and LASSO Cox regression.

Identification of CSPG4 as a promising target for translational combinatorial approaches in osteosarcoma

Background:

Osteosarcoma (OSA) is a highly metastatic pediatric bone tumor. Adjuvant chemotherapy and surgical resection represent standard treatments; however, the prognosis is still poor. Effective strategies are urgently needed. Chondroitin sulfate proteoglycan (CSPG)4 is a transmembrane proteoglycan with a low expression in normal tissues but high expression in several solid tumors, where it plays a central tumorigenic role. Therefore, it represents a promising therapeutic target. The high homology between human and canine CSPG4 and the recognized translational power of canine tumors as preclinical models for human malignancies prompted us to evaluate CSPG4 expression and the consequences of its immune-targeting for both human and canine OSA treatment.

Methods:

We analyzed CSPG4 overexpression in human and canine OSA samples and its significance for the survival of OSA patients. We exploited functional in vitro experiments to assess the antitumor potential of CSPG4 immune-targeting.

Results:

CSPG4 is overexpressed in OSA and has possible clinical implications as suggested by an evident correlation between CSPG4 overexpression and a shorter survival for both OSA-affected humans and dogs. The potential of CSPG4 immune-targeting for OSA treatment came from the ability of anti-CSPG4 monoclonal antibodies and sera, derived from human-CSPG4-DNA vaccinated canine patients, to significantly inhibit human and canine CSPG4-positive OSA cell proliferation, migration, and osteospheres generation. Moreover, CSPG4 immune-targeting has been shown to potentiate the effect of doxorubicin.

Conclusions:

Overall, these results provide the rationale to investigate the CSPG4 immune-targeting as a promising weapon for the treatment of CSPG4-positive OSA canine patients, to be successfully translated to a human setting.

E-Cadherin Represses Anchorage-Independent Growth in Sarcomas through Both Signaling and Mechanical Mechanisms

CDH1 (also known as E-cadherin), an epithelial-specific cell-cell adhesion molecule, plays multiple roles in maintaining adherens junctions, regulating migration and invasion, and mediating intracellular signaling. Downregulation of E-cadherin is a hallmark of epithelial-to-mesenchymal transition (EMT) and correlates with poor prognosis in multiple carcinomas. Conversely, upregulation of E-cadherin is prognostic for improved survival in sarcomas. Yet, despite the prognostic benefit of E-cadherin expression in sarcoma, the mechanistic significance of E-cadherin in sarcomas remains poorly understood. Here, by combining mathematical models with wet-bench experiments, we identify the core regulatory networks mediated by E-cadherin in sarcomas, and decipher their functional consequences. Unlike carcinomas, E-cadherin overexpression in sarcomas does not induce a mesenchymal-to-epithelial transition (MET). However, E-cadherin acts to reduce both anchorage-independent growth and spheroid formation of sarcoma cells. Ectopic E-cadherin expression acts to downregulate phosphorylated CREB1 (p-CREB) and the transcription factor, TBX2, to inhibit anchorage-independent growth. RNAi-mediated knockdown of TBX2 phenocopies the effect of E-cadherin on CREB levels and restores sensitivity to anchorage-independent growth in sarcoma cells. Beyond its signaling role, E-cadherin expression in sarcoma cells can also strengthen cell-cell adhesion and restricts spheroid growth through mechanical action. Together, our results demonstrate that E-cadherin inhibits sarcoma aggressiveness by preventing anchorage-independent growth. IMPLICATIONS: We highlight how E-cadherin can restrict aggressive behavior in sarcomas through both biochemical signaling and biomechanical effects.

Integration of genomic copy number variations and chemotherapy-response biomarkers in pediatric sarcoma

Background

While most pediatric sarcomas respond to front-line therapy, some bone sarcomas do not show radiographic response like soft-tissue sarcomas (rhabdomyosarccomas) but do show 90% necrosis. Though, new therapies are urgently needed to improve survival and quality of life in pediatric patients with sarcomas. Complex chromosomal aberrations such as amplifications and deletions of DNA sequences are frequently observed in pediatric sarcomas. Evaluation of copy number variations (CNVs) associated with pediatric sarcoma patients at the time of diagnosis or following therapy offers an opportunity to assess dysregulated molecular targets and signaling pathways that may drive sarcoma development, progression, or relapse. The objective of this study was to utilize publicly available data sets to identify potential predictive biomarkers of chemotherapeutic response in pediatric Osteosarcoma (OS), Rhabdomyosarcoma (RMS) and Ewing’s Sarcoma Family of Tumors (ESFTs) based on CNVs following chemotherapy (OS n = 117, RMS n = 64, ESFTs n = 25 tumor biopsies).

Methods

There were 206 CNV profiles derived from pediatric sarcoma biopsies collected from the public databases TARGET and NCBI-Gene Expression Omnibus (GEO). Through our comparative genomic analyses of OS, RMS, and ESFTs and 22,255 healthy individuals called from the Database of Genomic Variants (DGV), we identified CNVs (amplifications and deletions) pattern of genomic instability in these pediatric sarcomas. By integrating CNVs of Cancer Cell Line Encyclopedia (CCLE) identified in the pool of genes with drug-response data from sarcoma cell lines (n = 27) from Cancer Therapeutics Response Portal (CTRP) Version 2, potential predictive biomarkers of therapeutic response were identified.

Results

Genes associated with survival and/recurrence of these sarcomas with statistical significance were found on long arm of chromosome 8 and smaller aberrations were also identified at chromosomes 1q, 12q and x in OS, RMS, and ESFTs. A pool of 63 genes that harbored amplifications and/or deletions were frequently associated with recurrence across OS, RMS, and ESFTs. Correlation analysis of CNVs from CCLE with drug-response data of CTRP in 27 sarcoma cell lines, 33 CNVs out of 63 genes correlated with either sensitivity or resistance to 17 chemotherapies from which actionable CNV signatures such as IGF1R, MYC, MAPK1, ATF1, and MDM2 were identified. These CNV signatures could potentially be used to delineate patient populations that will respond versus those that will not respond to a particular chemotherapy.

Conclusions

The large-scale analyses of CNV-drug screening provides a platform to evaluate genetic alterations across aggressive pediatric sarcomas. Additionally, this study provides novel insights into the potential utilization of CNVs as not only prognostic but also as predictive biomarkers of therapeutic response. Information obtained in this study may help guide and prioritize patient-specific therapeutic options in pediatric bone and soft-tissue sarcomas.

Electronic supplementary material

The online version of this article (10.1186/s12920-018-0456-5) contains supplementary material, which is available to authorized users.

The microRNAs miR-449a and miR-424 suppress osteosarcoma by targeting cyclin A2 expression

MicroRNAs of the miR-16 and miR-34 families have been reported to inhibit cell cycle progression, and their loss has been linked to oncogenic transformation. Utilizing a high-throughput, genome-wide screen for miRNAs and mRNAs that are differentially regulated in osteosarcoma (OS) cell lines, we report that miR-449a and miR-424, belonging to the miR-34 and miR-16 families, respectively, target the major S/G₂ phase cyclin, cyclin A2 (CCNA2), in a bipartite manner. We found that the 3'-UTR of CCNA2 is recognized by miR-449a, whereas the CCNA2 coding region is targeted by miR-424. Of note, we observed loss of both miR-449a and miR-424 in OS, resulting in derepression of CCNA2 and appearance of aggressive cancer phenotypes. Ectopic expression of miR-449a and miR-424 significantly decreased cyclin A2 levels and inhibited proliferation rate, migratory potential, and colony-forming ability of OS cells. To further probe the roles of miR-449a and miR-424 in OS, we developed an OS mouse model by intraosseous injection of U2OS cells into the tibia bone of NOD-scid mice, which indicated that miR-449a and miR-424 co-expression suppresses tumor growth. On the basis of this discovery, we analyzed the gene expression of human OS biopsy samples, revealing that miR-449a and miR-424 are both down-regulated, whereas cyclin A2 is significantly up-regulated in these OS samples. In summary, the findings in our study highlight that cyclin A2 repression by miRNAs of the miR-16 and miR-34 families is lost in aggressive OS.

The cancer-related transcription factor RUNX2 modulates expression and secretion of the matricellular protein osteopontin in osteosarcoma cells to promote adhesion to endothelial pulmonary cells and lung metastasis

Osteosarcomas are bone tumors that frequently metastasize to the lung. Aberrant expression of the transcription factor, runt-related transcription factor 2 (RUNX2), is a key pathological feature in osteosarcoma and associated with loss of p53 and miR-34 expression. Elevated RUNX2 may transcriptionally activate genes mediating tumor progression and metastasis, including the RUNX2 target gene osteopontin (OPN/SPP1). This gene encodes a secreted matricellular protein produced by osteoblasts to regulate bone matrix remodeling and tissue calcification. Here we investigated whether and how the RUNX2/OPN axis regulates lung metastasis of osteosarcoma. Importantly, RUNX2 depletion attenuates lung metastasis of osteosarcoma cells in vivo. Using next-generation RNA-sequencing, protein-based assays, as well as the loss- and gain-of-function approaches in selected osteosarcoma cell lines, we show that osteopontin messenger RNA levels closely correlate with RUNX2 expression and that RUNX2 controls the levels of secreted osteopontin. Elevated osteopontin levels promote heterotypic cell-cell adhesion of osteosarcoma cells to human pulmonary microvascular endothelial cells, but not in the presence of neutralizing antibodies. Collectively, these findings indicate that the RUNX2/OPN axis regulates the ability of osteosarcoma cells to attach to pulmonary endothelial cells as a key step in metastasis of osteosarcoma cells to the lung.

HACE1 is a potential tumor suppressor in osteosarcoma

Osteosarcoma is a malignant bone sarcoma characterized by extensive genomic disruption and a propensity for metastatic spread. Osteoid production suggests a close relationship with normal osteoblasts, and the latter are the presumptive cell of origin of this disease. The HACE1 gene, localized to human chromosome 6q21, encodes the HACE1 HECT E3 ligase, a tumor suppressor in diverse tumors that acts in part by targeting the activated form of RAC1 GTPase for proteasomal degradation. Disruption or loss of 6q21 is relatively common in osteosarcomas, and Hace1-/-/Tp53+/- mice frequently develop osteosarcomas, in contrast to Tp53+/- mice, which do not. This suggests an unexplored link between HACE1 loss and osteosarcoma. Here we compared HACE1 expression in normal osteoblasts and osteosarcoma cell lines in vitro by western blotting and quantitative RT-PCR, and in human osteosarcoma specimens by immunohistochemistry. Both HACE1 transcript and protein levels were reduced in osteosarcoma compared to osteoblasts in vitro. Reduced HACE1 expression in osteosarcoma tumors was observed in 76% of cases and associated with high-grade lesions. Further, clonally derived pairs of high and low metastatic osteosarcoma cell lines showed significant downregulation in the high compared to corresponding low metastatic cells. Ectopic expression of HACE1 markedly inhibited anchorage-independent growth and cell motility of HACE1 osteosarcoma cell lines, and was associated with reduced RAC1 activation and decreased reactive oxygen species (ROS). Finally, HACE1 overexpression blocked osteosarcoma xenograft growth and dramatically reduced pulmonary metastases. These findings point to a potential tumor suppressor function for HACE1 in osteosarcoma.

Genome-Informed Targeted Therapy for Osteosarcoma

Osteosarcoma is a highly aggressive cancer for which treatment has remained essentially unchanged for more than 30 years. Osteosarcoma is characterized by widespread and recurrent somatic copy-number alterations (SCNA) and structural rearrangements. In contrast, few recurrent point mutations in protein-coding genes have been identified, suggesting that genes within SCNAs are key oncogenic drivers in this disease. SCNAs and structural rearrangements are highly heterogeneous across osteosarcoma cases, suggesting the need for a genome-informed approach to targeted therapy. To identify patient-specific candidate drivers, we used a simple heuristic based on degree and rank order of copy-number amplification (identified by whole-genome sequencing) and changes in gene expression as identified by RNA sequencing. Using patient-derived tumor xenografts, we demonstrate that targeting of patient-specific SCNAs leads to significant decrease in tumor burden, providing a road map for genome-informed treatment of osteosarcoma. SIGNIFICANCE: Osteosarcoma is treated with a chemotherapy regimen established 30 years ago. Although osteosarcoma is genomically complex, we hypothesized that tumor-specific dependencies could be identified within SCNAs. Using patient-derived tumor xenografts, we found a high degree of response for "genome-matched" therapies, demonstrating the utility of a targeted genome-informed approach.This article is highlighted in the In This Issue feature, p. 1.

Role of Activator Protein-1 Complex on the Phenotype of Human Osteosarcomas Generated from Mesenchymal Stem Cells

Osteosarcoma (OS) is a highly aggressive bone tumor that usually arises intramedullary at the extremities of long bones. Due to the fact that the peak of incidence is in the growth spurt of adolescence, the specific anatomical location, and the heterogeneity of cells, it is believed that osteosarcomagenesis is a process associated with bone development. Different studies in murine models showed that the tumor-initiating cell in OS could be an uncommitted mesenchymal stem cell (MSC) developing in a specific bone microenvironment. However, only a few studies have reported transgene-induced human MSCs transformation and mostly obtained undifferentiated sarcomas. In our study, we demonstrate that activator protein 1 family members induce osteosarcomagenesis in immortalized hMSC. c-JUN or c-JUN/c-FOS overexpression act as tumorigenic factors generating OS with fibroblastic or pleomorphic osteoblastic phenotypes, respectively. Stem Cells 2018;36:1487-1500.

Cross-species genomics identifies DLG2 as a tumor suppressor in osteosarcoma

Leveraging the conserved cancer genomes across mammals has the potential to transform driver gene discovery in orphan cancers. Here, we combine cross-species genomics with validation across human–dog–mouse systems to uncover a new bone tumor suppressor gene. Comparative genomics of spontaneous human and dog osteosarcomas (OS) expose Disks Large Homolog 2 (DLG2) as a tumor suppressor candidate. DLG2 copy number loss occurs in 42% of human and 56% of canine OS. Functional validation through pertinent human and canine OS DLG2-deficient cell lines identifies a regulatory role of DLG2 in cell division, migration and tumorigenesis. Moreover, osteoblast-specific deletion of Dlg2 in a clinically relevant genetically engineered mouse model leads to acceleration of OS development, establishing DLG2 as a critical determinant of OS. This widely applicable cross-species approach serves as a platform to expedite the search of cancer drivers in rare human malignancies, offering new targets for cancer therapy.