Clinical and Molecular Analysis of Pathologic Fracture-associated Osteosarcoma: MicroRNA profile Is Different and Correlates with Prognosis.

BACKGROUND: MicroRNAs are small, noncoding RNAs that regulate the expression of posttranslational genes. The presence of some specific microRNAs has been associated with increased risk of both local recurrence and metastasis and worse survival in patients with osteosarcoma. Pathologic fractures in osteosarcoma are considered to be more the manifestation of a neoplasm with a more aggressive biological behavior than the cause itself of worse prognosis. However, this has not been proved at the biological or molecular level. Currently, there has not been a microRNA profiling study of patients who have osteosarcoma with and without pathologic fractures that has described differences in terms of microRNA profiling between these two groups and their correlation with biologic behavior. QUESTIONS/PURPOSES: (1) In patients with osteosarcoma of the extremities, how do the microRNA profiles of those with and without pathologic fractures compare? (2) What relationship do microRNAs have with local recurrence, risk of metastasis, disease-specific survival, and overall survival in osteosarcoma patients with pathologic fractures? METHODS: Between 1994 and 2013, 217 patients were diagnosed and treated at our institution for osteosarcoma of the extremities. Patients were excluded if (1) they underwent oncologic resection of the osteosarcoma at an outside institution (two patients) or (2) they were diagnosed with an extraskeletal osteosarcoma (29 patients) or (3) they had less than 1 year of clinical follow-up and no oncologic outcome (local recurrence, metastasis, or death) (four patients). A total of 182 patients were eligible. Of those, 143 were high-grade osteosarcomas. After evaluation of tumor samples before chemotherapy treatment, a total of 80 consecutive samples were selected for sequencing. Demographic and clinical comparison between the sequenced and non-sequenced patients did not demonstrate any differences, confirming that both groups were comparable. Diagnostic samples from the extremities of 80 patients with high-grade extremity osteosarcomas who had not yet received chemotherapy underwent microRNA sequencing for an ongoing large-scale osteosarcoma genome profiling project at our institution. Six samples were removed after a second look by a musculoskeletal pathologist who verified cellularity and quality of samples to be sequenced, leaving a total of 74 patients. Of these, two samples were removed as they were confirmed to be pelvic tumors in a second check after sequencing. The final study sample was 72 patients (11 patients with pathologic fractures and 61 without). Sequencing data were correlated with fractures and local recurrence, risk of metastasis, disease-specific survival, and overall survival through Kaplan-Meier analyses. RESULTS: Several microRNAs were expressed differently between the two groups. Among the markers with the highest differential expression (edgeR and DESeq algorithms), Hsa-mIR 656-3p, hsa-miR 493-5p, and hsa-miR 381-3p were upregulated in patients with pathologic fractures, whereas hsa-miR 363, hsa-miR 885-5p, and has-miR 20b-5p were downregulated. The highest differential expression fracture and nonfracture-associated microRNA markers also distinguished groups of patients with different metastasis risk, a well as different disease-specific and overall survival. Furthermore, the profile of pathologic fractures demonstrated a higher differential expression for microRNA markers that were previously associated with a higher risk of metastasis and lower survival rates in patients with osteosarcoma. CONCLUSIONS: In patients who have osteosarcoma, the microRNA profiles of those with pathologic fractures are different than of patients without pathologic fractures. The highest differential expression mircroRNA molecules in patients with pathologic fractures predict also higher risk of metastatic disease as well as worse disease-specific survival and overall survival. Furthermore, we found higher differential expression of microRNAs in the pathologic fracture group previously associated with poor prognosis. The higher risk of metastasis and poorer overall survival in patients with pathologic fractures is inherent to tumor aggressive biologic behavior. It is plausible that the fracture itself is not the direct cause of worse prognosis but another manifestation of tumor biologic aggressiveness. Identification of these molecules through liquid biopsies may help to determine which patients may benefit from surgery before fractures occur. The same technology can be applied to identify patterns of response to conventional chemotherapy, assisting in more specific and accurate systemic therapy. LEVEL OF EVIDENCE LEVEL: III, prognostic study.

Expression and Clinical Implication of Autophagy-Associated Protein p62 in Osteosarcoma.

PURPOSE: Recent studies highlight the role of autophagy in cancer tumorigenesis, recurrence, metastasis, and chemoresistance. p62 is an adapter protein that is crucial for the autophagy pathway. In this study, we will describe the expression of p62 and its correlation with clinic prognosis in osteosarcoma. METHODS: Western blot was used to test the expression of p62 in osteosarcoma cell lines (U2OS, KHOS, MG63, Saos-2, U2OSR2, KHOSR2, and 143B). A tissue microarray (TMA) was analyzed by immunohistochemistry to determine the expression levels of p62 in osteosarcoma patients and evaluate any correlation between p62 and clinical characteristics in osteosarcoma patients. RESULTS: p62 was expressed differently in all cell lines. The TMA also showed differential expression in osteosarcoma tissues. Seventy-five of 79 (94.9%) patient tissues exhibited p62 immunostaining, ranging from no staining (4 of 97, 5.1%) to 1+ staining (40 of 79, 50.6%), 2+ staining (17 of 79, 21.5%), and 3+ staining (18 of 79, 22.8%). The low staining (1+) was classified as the p62 weak group (50.6%), the medium staining (2+) and intense staining (3+) were classified as the p62 strong group (44.3%). Analyzing the clinical data of the osteosarcoma TMA, we found that the 5-year survival rate of patients with weak p62 expression was significantly lower than that of the patients with strong p62 expression (p = 0.0165). Furthermore, the decreased p62 expression may be associated with higher metastatic and chemoresistant rates in osteosarcoma patients. CONCLUSION: Our results suggest that p62 may be an effective predictor of prognosis and a potential target for therapy in osteosarcoma.

Expression and role of autophagy-associated p62 (SQSTM1) in multidrug resistant ovarian cancer.

OBJECTIVE: Multidrug resistance is the major cause of treatment failure in ovarian cancer. p62 (SQSTM1) is a multifunctional protein involved in multiple cellular processes including proliferation, drug sensitivity and autophagy-associated cancer cell growth. However, the role of p62 in drug resistance remains controversial. METHODS: In this study, we examined p62 expression by immunohistochemistry in a unique ovarian cancer tissue microarray (TMA), which was constructed with paired primary, metastatic, and recurrent tumor tissues. The expression levels of p62 and autophagy related proteins were evaluated in two panels of human cancer cell lines by western blot. Cell viabilities were determined by MTT assay after exposure ovarian cancer cells to different concentrations of paclitaxel alone or in combination with autophagy inhibitors. RESULTS: Both the metastatic and recurrent tumor tissues expressed less p62 than the patient-matched primary tumor. A significant inverse correlation has been found between p62 expression and both the disease-free survival and overall survival. Additionally, multidrug resistant cancer cell lines expressed lower levels of p62 as compared with their parental drug sensitive cell lines. Importantly, inhibition of autophagy enhanced paclitaxel sensitivity in drug resistant ovarian cancer cells. Furthermore, the wound healing assay exhibited that the inhibition of autophagy significantly decreased resistant ovarian cancer cell migration in vitro. CONCLUSION: Our findings highlight the potential of p62 as a new prognostic marker for ovarian cancer patients and p62's associated autophagy pathway may be a promising therapeutic target to prevent metastasis, recurrence and to reverse drug resistance in ovarian cancer.

Potentials of Long Noncoding RNAs (LncRNAs) in Sarcoma: From Biomarkers to Therapeutic Targets.

Sarcoma includes some of the most heterogeneous tumors, which make the diagnosis, prognosis and treatment of these rare yet diverse neoplasms especially challenging. Long noncoding RNAs (lncRNAs) are important regulators of cancer initiation and progression, which implies their potential as neoteric prognostic and diagnostic markers in cancer, including sarcoma. A relationship between lncRNAs and sarcoma pathogenesis and progression is emerging. Recent studies demonstrate that lncRNAs influence sarcoma cell proliferation, metastasis, and drug resistance. Additionally, lncRNA expression profiles are predictive of sarcoma prognosis. In this review, we summarize contemporary advances in the research of lncRNA biogenesis and functions in sarcoma. We also highlight the potential for lncRNAs to become innovative diagnostic and prognostic biomarkers as well as therapeutic targets in sarcoma.

HERVs establish a distinct molecular subtype in stage II/III colorectal cancer with poor outcome.

Colorectal cancer (CRC) is one of the most lethal malignancies. The extreme heterogeneity in survival rate is driving the need for new prognostic biomarkers. Human endogenous retroviruses (hERVs) have been suggested to influence tumor progression, oncogenesis and elicit an immune response. We examined multiple next-generation sequencing (NGS)-derived biomarkers in 114 CRC patients with paired whole-exome and whole-transcriptome sequencing (WES and WTS, respectively). First, we demonstrate that the median expression of hERVs can serve as a potential biomarker for prognosis, relapse, and resistance to chemotherapy in stage II and III CRC. We show that hERV expression and CD8+ tumor-infiltrating T-lymphocytes (TILs) synergistically stratify overall and relapse-free survival (OS and RFS): the median OS of the CD8-/hERV+ subgroup was 29.8 months compared with 37.5 months for other subgroups (HR = 4.4, log-rank P < 0.001). Combing NGS-based biomarkers (hERV/CD8 status) with clinicopathological factors provided a better prediction of patient survival compared to clinicopathological factors alone. Moreover, we explored the association between genomic and transcriptomic features of tumors with high hERV expression and establish this subtype as distinct from previously described consensus molecular subtypes of CRC. Overall, our results underscore a previously unknown role for hERVs in leading to a more aggressive subtype of CRC.

MicroRNA-mRNA networks define translatable molecular outcome phenotypes in osteosarcoma.

There is a lack of well validated prognostic biomarkers in osteosarcoma, a rare, recalcitrant disease for which treatment standards have not changed in over 20 years. We performed microRNA sequencing in 74 frozen osteosarcoma biopsy samples, constituting the largest single center translationally analyzed osteosarcoma cohort to date, and we separately analyzed a multi-omic dataset from a large NCI supported national cooperative group cohort. We validated the prognostic value of candidate microRNA signatures and contextualized them in relevant transcriptomic and epigenomic networks. Our results reveal the existence of molecularly defined phenotypes associated with outcome independent of clinicopathologic features. Through machine learning based integrative pharmacogenomic analysis, the microRNA biomarkers identify novel therapeutics for stratified application in osteosarcoma. The previously unrecognized osteosarcoma subtypes with distinct clinical courses and response to therapy could be translatable for discerning patients appropriate for more intensified, less intensified, or alternate therapeutic regimens.

Author Correction: Intragenic antagonistic roles of protein and circRNA in tumorigenesis.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Intragenic antagonistic roles of protein and circRNA in tumorigenesis.

circRNAs arise from back splicing events during mRNA processing, and when deregulated can play an active role in cancer. Here we characterize a new circRNA (circPOK) encoded by the Zbtb7a gene (also kown as POKEMON, LRF) in the context of mesenchymal tumor progression. circPOK functions as a non-coding proto-oncogenic RNA independently and antithetically to its linear transcript counterpart, which acts as a tumor suppressor by encoding the Pokemon transcription factor. We find that circPOK regulates pro-proliferative and pro-angiogenic factors by co-activation of the ILF2/3 complex. Importantly, the expression of Pokemon protein and circRNA is aberrantly uncoupled in cancer through differential post-transcriptional regulation. Thus, we identify a novel type of genetic unit, the iRegulon, that yields biochemically distinct RNA products, circular and linear, with diverse and antithetical functions. Our findings further expand the cellular repertoire towards the control of normal biological outputs, while aberrant expression of such components may underlie disease pathogenesis including cancer.

Advances in chromosomal translocations and fusion genes in sarcomas and potential therapeutic applications.

Chromosomal translocations and fusion genes are very common in human cancer especially in subtypes of sarcomas, such as rhabdomyosarcoma, Ewing's sarcoma, synovial sarcoma and liposarcoma. The discovery of novel chromosomal translocations and fusion genes in different tumors are due to the advancement of next-generation sequencing (NGS) technologies such as whole genome sequencing. Recently, many novel chromosomal translocations and gene fusions have been identified in different types of sarcoma through NGS approaches. In addition to previously known sarcoma fusion genes, these novel specific fusion genes and associated molecular events represent important targets for novel therapeutic approaches in the treatment of sarcomas. This review focuses on recent advances in chromosomal translocations and fusion genes in sarcomas and their potential therapeutic applications in the treatment of sarcomas.

Inhibition of cyclin-dependent kinase 4 as a potential therapeutic strategy for treatment of synovial sarcoma.

Synovial sarcoma is a highly aggressive but rare form of soft tissue malignancy that primarily affects the extremities of the arms or legs, for which current chemotherapeutic agents have not been proven to be very effective. The cyclin-dependent kinase 4/6-retinoblastoma protein (CDK4/6-Rb) pathway of cell cycle control is known to be aberrant in a large proportion of cancers. Recently, CDK4 inhibitors have successfully been used pre-clinically for the treatment of many human cancers, and in 2015, following the success of clinical trials, the FDA approved the first selective CDK4/6 inhibitor, palbociclib, for the treatment of endocrine therapy resistant breast cancers. However, the expression and therapeutic potential of targeting CDK4 in synovial sarcoma remains unclear. In the present study, we report that CDK4 is highly expressed in human synovial sarcoma, and high CDK4 expressions are associated with poor prognosis in sarcomas patients and the clinical stage and the TNM grade in synovial sarcoma patients. Knockdown of CDK4 with specific small interference RNAs inhibits cell proliferation and enhances apoptotic effects in synovial sarcoma cells. CDK4 inhibitor palbociclib suppresses synovial sarcoma cell proliferation and growth in a dose and time-dependent manner. Palbociclib also inhibits the CDK4/6-Rb signaling pathway and promotes cell apoptosis without changing CDK4/6 protein levels, suggesting that palbociclib only represses the hyper-activation, not the expression of CDK4/6. Flow cytometry analysis reveals that palbociclib induces G1 cell-cycle arrest and apoptotic effects by targeting the CDK4/6-Rb pathway in synovial sarcoma cells. Furthermore, wound healing assays demonstrate that inhibition of the CDK4/6-Rb pathway by palbociclib significantly decreases synovial sarcoma cell migration in vitro. Our study highlights the importance of the CDK4/6-Rb pathway in human synovial sarcoma pathogenesis, and the role of the current selective CDK4/6 inhibitor, palbociclib, as a potential promising targeted therapeutic agent in the treatment of human synovial sarcoma.

Aberration of p73 Promoter Methylation in Chondrosarcoma.

BACKGROUND: p73 is a tumor-suppressor gene with significant homology to p53. Abnormal promoter methylation of p73 is present in different types of cancer. However, the promoter methylation status of p73 in chondrosarcoma (CS) is unknown. MATERIALS AND METHODS: p73 promoter methylation status was evaluated by quantitative polymerase chain reaction (PCR), p73 protein expression by western blot, and the relationship between p73 methylation and clinical data was analyzed. RESULTS: In 42 tumor tissues with CS, we found that three cases (7%) maintained methylation levels between 51% and 75%, and 39 cases (93%) had levels between 76% and 100%. p73 methylation level was significantly (p<0.05) positively associated with histological grade. Loss of p73 protein expression was correlated with high methylation of the p73 promoter; p73 expression was restored after exposure to a demethylating drug. CONCLUSION: p73 is epigenetically silenced in CS due to promoter methylation, which suggests the utility of p73 methylation as a biomarker.

Expression and Therapeutic Potential of SOX9 in Chordoma.

Purpose: Conventional chemotherapeutic agents are ineffective in the treatment of chordoma. We investigated the functional roles and therapeutic relevance of the sex-determining region Y (SRY)-box 9 (SOX9) in chordoma.Experimental Design: SOX9 expression was examined by immunohistochemistry (IHC) using 50 chordoma tissue samples. SOX9 expression in chordoma cell lines was examined by Western blot and immunofluorescent assays. We used synthetic human SOX9 siRNA to inhibit the expression of SOX9. Cell proliferation ability and cytotoxicity of inhibiting SOX9 were assessed by 3-(4, 5-dimethylthiazolyl-2)-2, 5-diphenyltetrazolium bromide (MTT) and clonogenic assays. The effect of SOX9 knockdown on chordoma cell motility was evaluated by a wound-healing assay and a Transwell invasion chamber assay. Knockdown of SOX9 induced apoptosis, cell-cycle arrest, as well as decreased expression of cancer stem cell markers were determined by Western blot and flow cytometric assays. The effect of the combination of SOX9 siRNA and the chemotherapeutic drug doxorubicin/cisplatin on chordoma cells was assessed by an MTT assay.Results: Tissue microarray and IHC analysis showed that SOX9 is broadly expressed in chordomas and that higher expression levels of SOX9 correlated with a poor prognosis. RNA interference (RNAi)-mediated knockdown of SOX9 inhibited chordoma cell growth, decreased cell motility, and induced apoptosis as well as cell-cycle arrest. Moreover, the combination of SOX9 inhibition and chemotherapeutic drugs had an enhanced anti-cancer effect on chordoma cells.Conclusions: Our results demonstrate that SOX9 plays a crucial role in chordoma. Targeting SOX9 provides a new rationale for treatment of chordoma. Clin Cancer Res; 23(17); 5176-86. ©2017 AACR.

An imprinted non-coding genomic cluster at 14q32 defines clinically relevant molecular subtypes in osteosarcoma across multiple independent datasets.

BACKGROUND: A microRNA (miRNA) collection on the imprinted 14q32 MEG3 region has been associated with outcome in osteosarcoma. We assessed the clinical utility of this miRNA set and their association with methylation status. METHODS: We integrated coding and non-coding RNA data from three independent annotated clinical osteosarcoma cohorts (n = 65, n = 27, and n = 25) and miRNA and methylation data from one in vitro (19 cell lines) and one clinical (NCI Therapeutically Applicable Research to Generate Effective Treatments (TARGET) osteosarcoma dataset, n = 80) dataset. We used time-dependent receiver operating characteristic (tdROC) analysis to evaluate the clinical value of candidate miRNA profiles and machine learning approaches to compare the coding and non-coding transcriptional programs of high- and low-risk osteosarcoma tumors and high- versus low-aggressiveness cell lines. In the cell line and TARGET datasets, we also studied the methylation patterns of the MEG3 imprinting control region on 14q32 and their association with miRNA expression and tumor aggressiveness. RESULTS: In the tdROC analysis, miRNA sets on 14q32 showed strong discriminatory power for recurrence and survival in the three clinical datasets. High- or low-risk tumor classification was robust to using different microRNA sets or classification methods. Machine learning approaches showed that genome-wide miRNA profiles and miRNA regulatory networks were quite different between the two outcome groups and mRNA profiles categorized the samples in a manner concordant with the miRNAs, suggesting potential molecular subtypes. Further, miRNA expression patterns were reproducible in comparing high-aggressiveness versus low-aggressiveness cell lines. Methylation patterns in the MEG3 differentially methylated region (DMR) also distinguished high-aggressiveness from low-aggressiveness cell lines and were associated with expression of several 14q32 miRNAs in both the cell lines and the large TARGET clinical dataset. Within the limits of available CpG array coverage, we observed a potential methylation-sensitive regulation of the non-coding RNA cluster by CTCF, a known enhancer-blocking factor. CONCLUSIONS: Loss of imprinting/methylation changes in the 14q32 non-coding region defines reproducible previously unrecognized osteosarcoma subtypes with distinct transcriptional programs and biologic and clinical behavior. Future studies will define the precise relationship between 14q32 imprinting, non-coding RNA expression, genomic enhancer binding, and tumor aggressiveness, with possible therapeutic implications for both early- and advanced-stage patients.

Getting to the edge: protein dynamical networks as a new frontier in plant-microbe interactions.

A systems perspective on diverse phenotypes, mechanisms of infection, and responses to environmental stresses can lead to considerable advances in agriculture and medicine. A significant promise of systems biology within plants is the development of disease-resistant crop varieties, which would maximize yield output for food, clothing, building materials, and biofuel production. A systems or "-omics" perspective frames the next frontier in the search for enhanced knowledge of plant network biology. The functional understanding of network structure and dynamics is vital to expanding our knowledge of how the intercellular communication processes are executed. This review article will systematically discuss various levels of organization of systems biology beginning with the building blocks termed "-omes" and ending with complex transcriptional and protein-protein interaction networks. We will also highlight the prevailing computational modeling approaches of biological regulatory network dynamics. The latest developments in the "-omics" approach will be reviewed and discussed to underline and highlight novel technologies and research directions in plant network biology.