miR-486-5p expression is regulated by DNA methylation in osteosarcoma.

BACKGROUND: Osteosarcoma is the most common primary malignant tumour of bone occurring in children and young adolescents and is characterised by complex genetic and epigenetic changes. The miRNA miR-486-5p has been shown to be downregulated in osteosarcoma and in cancer in general. RESULTS: To investigate if the mir-486 locus is epigenetically regulated, we integrated DNA methylation and miR-486-5p expression data using cohorts of osteosarcoma cell lines and patient samples. A CpG island in the promoter of the ANK1 host gene of mir-486 was shown to be highly methylated in osteosarcoma cell lines as determined by methylation-specific PCR and direct bisulfite sequencing. High methylation levels were seen for osteosarcoma patient samples, xenografts and cell lines based on quantitative methylation-specific PCR. 5-Aza-2'-deoxycytidine treatment of osteosarcoma cell lines caused induction of miR-486-5p and ANK1, indicating common epigenetic regulation in osteosarcoma cell lines. When overexpressed, miR-486-5p affected cell morphology. CONCLUSIONS: miR-486-5p represents a highly cancer relevant, epigenetically regulated miRNA in osteosarcoma, and this knowledge contributes to the understanding of osteosarcoma biology.

Use of liquid biopsies to monitor disease progression in a sarcoma patient: a case report.

BACKGROUND: Many patients experience local recurrence or metastases after receiving potentially curative treatment, and early detection of these events is important for disease control. Recent technological advances make it possible to use blood plasma containing circulating cell-free tumour DNA (ctDNA) as a liquid biopsy. In this case report we show how serial liquid biopsies can be used to monitor the disease course and detect disease recurrence in a sarcoma patient. CASE PRESENTATION: A 55-year-old male presented with a rapidly growing, painful palpable mass in the left groin region, and a biopsy revealed a high-grade malignant spindle cell sarcoma. No metastases were detected on radiologic imaging scans. Using targeted resequencing with a custom 900 cancer gene panel, eight somatic mutations among them KRAS and NF1, were identified in the primary tumour. Targeted resequencing of plasma cell-free DNA (ctDNA) collected before and after surgery and at disease progression confirmed the presence of six of eight mutations at all three time points. The ctDNA level, estimated from the somatic allele frequencies of these six mutations, was high in plasma taken at the time of surgery, at levels similar to the primary tumour. Detection of low levels of ctDNA three days after surgery indicated persistent microscopic disease. Repeated radiologic imaging six weeks postoperatively showed widespread metastatic disease in the lungs, skeleton and the pelvic region. At this time point there was a dramatic increase in the ctDNA level, reflecting the disease progression of the patient. The patient had an unusually aggressive cancer, and succumbed to the disease 13 weeks after surgery. CONCLUSIONS: This case report demonstrated that targeted resequencing of ctDNA from longitudinal collected plasma can be used to monitor disease progression in a soft tissue sarcoma patient, including manifestation of metastatic disease. The ctDNA represented the genomic profile of the tumour, supporting clinical use of liquid biopsies to identify tumour-specific mutations as well as recurrent disease.

Global gene expression profiling of human osteosarcomas reveals metastasis-associated chemokine pattern.

Global gene expression analysis was performed on a panel of 23 osteosarcoma samples of primary and metastatic origin using the Applied Biosystems Gene Expression Array System. When comparing the primary tumours with the metastases, we found a significantly increased expression of genes involved in immunological processes, for example coding for cytokines and chemokines, in the metastatic samples. In addition, a comparison of the gene expression in primary samples from patients with or without metastases demonstrated that patients who later developed metastases had high expression of the chemokine (C-X-C motif) receptor 4 (CXCR4), similar to the metastatic samples, suggesting that these signal molecules play an important role in promoting metastasis. Increased knowledge of mechanisms and interactions between specified molecular signalling pathways in osteosarcomas could lead to a more rational strategy for development of targeted therapy.

DNA copy number changes in high-grade malignant peripheral nerve sheath tumors by array CGH.

BACKGROUND: Malignant peripheral nerve sheath tumors (MPNSTs) are rare and highly aggressive soft tissue tumors showing complex chromosomal aberrations. In order to identify recurrent chromosomal regions of gain and loss, and thereby novel gene targets of potential importance for MPNST development and/or progression, we have analyzed DNA copy number changes in seven high-grade MPNSTs using microarray-based comparative genomic hybridization (array CGH). RESULTS: Considerable more gains than losses were observed, and the most frequent minimal recurrent regions of gain included 1q24.1-q24.2, 1q24.3-q25.1, 8p23.1-p12, 9q34.11-q34.13 and 17q23.2-q25.3, all gained in five of seven samples. The 17q23.2-q25.3 region was gained in all five patients with poor outcome and not in the two patients with disease-free survival. cDNA microarray analysis and quantitative real-time reverse transcription PCR were used to investigate expression of genes located within these regions. The gene lysyl oxidase-like 2 (LOXL2) was identified as a candidate target for the 8p23.1-p12 gain. Within 17q, the genes topoisomerase II-alpha (TOP2A), ets variant gene 4 (E1A enhancer binding protein, E1AF) (ETV4) and baculoviral IAP repeat-containing 5 (survivin) (BIRC5) showed increased expression in all samples compared to two benign tumors. Increased expression of these genes has previously been associated with poor survival in other malignancies, and for TOP2A, in MPNSTs as well. In addition, we have analyzed the expression of five micro RNAs located within the 17q23.2-q25.3 region, but none of them showed high expression levels compared to the benign tumors. CONCLUSION: Our study shows the potential of using DNA copy number changes obtained by array CGH to predict the prognosis of MPNST patients. Although no clear correlations between the expression level and patient outcome were observed, the genes TOP2A, ETV4 and BIRC5 are interesting candidate targets for the 17q gain associated with poor survival.

Array comparative genomic hybridization reveals distinct DNA copy number differences between gastrointestinal stromal tumors and leiomyosarcomas.

Leiomyosarcomas are spindle cell tumors showing smooth muscle differentiation. Until recently, most gastrointestinal stromal tumors (GIST) were also classified as smooth muscle tumors, but now GISTs are recognized as a separate entity, defined as spindle cell and/or epithelioid tumors localized in the gastrointestinal tract. Using microarray-based comparative genomic hybridization (array CGH), we have created a detailed map of DNA copy number changes for 7 GISTs and 12 leiomyosarcomas. Considerable gains and losses of chromosomal segments were observed in both tumor types. The most frequent aberration observed in GISTs was loss of chromosomes 14 and 22, with minimal recurrent regions in 14q11.2-q32.33 (71% of the tumors) and 22q12.2-q13.31 (100%). In leiomyosarcomas, frequent loss of chromosome 10 and 13q was observed, with minimal recurrent regions in 10q21.3 (75%) and 13q14.2-q14.3 (75%). Recurrent high-level amplification of 17p13.1-p11.2 was detected in leiomyosarcomas. Expression profiling using cDNA microarrays revealed four candidate genes in this region with high expression (AURKB, SREBF1, MFAP4, and FLJ10847). Altered expression of AURKB and SREBF1 has been observed previously in other malignancies. Hierarchical clustering of all samples separated GISTs and leiomyosarcomas into two distinct clusters. Statistical analysis identified six chromosomal regions, 1p36.11-p13.1, 9q21.11-9q34.3, 14q11.2-q23.2, 14q31.3-q32.33, 15q24.3-q26.3, and 22q11.21-q13.31, which were significantly different in copy number between GISTs and leiomyosarcomas. Our results show the potential of using array comparative genomic hybridization to classify histologically similar tumors such as GISTs and leiomyosarcomas.

Evaluation of commercial DNA and RNA extraction methods for high-throughput sequencing of FFPE samples.

Nucleic acid material of adequate quality is crucial for successful high-throughput sequencing (HTS) analysis. DNA and RNA isolated from archival FFPE material are frequently degraded and not readily amplifiable due to chemical damage introduced during fixation. To identify optimal nucleic acid extraction kits, DNA and RNA quantity, quality and performance in HTS applications were evaluated. DNA and RNA were isolated from five sarcoma archival FFPE blocks, using eight extraction protocols from seven kits from three different commercial vendors. For DNA extraction, the truXTRAC FFPE DNA kit from Covaris gave higher yields and better amplifiable DNA, but all protocols gave comparable HTS library yields using Agilent SureSelect XT and performed well in downstream variant calling. For RNA extraction, all protocols gave comparable yields and amplifiable RNA. However, for fusion gene detection using the Archer FusionPlex Sarcoma Assay, the truXTRAC FFPE RNA kit from Covaris and Agencourt FormaPure kit from Beckman Coulter showed the highest percentage of unique read-pairs, providing higher complexity of HTS data and more frequent detection of recurrent fusion genes. truXTRAC simultaneous DNA and RNA extraction gave similar outputs as individual protocols. These findings show that although successful HTS libraries could be generated in most cases, the different protocols gave variable quantity and quality for FFPE nucleic acid extraction. Selecting the optimal procedure is highly valuable and may generate results in borderline quality specimens.

Noninvasive Detection of ctDNA Reveals Intratumor Heterogeneity and Is Associated with Tumor Burden in Gastrointestinal Stromal Tumor.

Molecular analysis of circulating tumor DNA (ctDNA) has a large potential for clinical application by capturing tumor-specific aberrations through noninvasive sampling. In gastrointestinal stromal tumor (GIST), analysis of KIT and PDGFRA mutations is important for therapeutic decisions, but the invasiveness of traditional biopsies limits the possibilities for repeated sampling. Using targeted next-generation sequencing, we have analyzed circulating cell-free DNA from 50 GIST patients. Tumor-specific mutations were detected in 16 of 44 plasma samples (36%) from treatment-naïve patients and in three of six (50%) patients treated with tyrosine kinase inhibitors. A significant association between detection of ctDNA and the modified National Institutes of Health risk classification was found. All patients with metastatic disease had detectable ctDNA, and tumor burden was the most important detection determinant. Median tumor size was 13.4 cm for patients with detectable mutation in plasma compared with 4.4 cm in patients without detectable mutation (P = 0.006). ctDNA analysis of a patient with disease progression on imatinib revealed that multiple resistance mutations were synchronously present, and detailed analysis of tumor tissue showed that these were spatially distributed in the primary tumor. Plasma samples taken throughout the course of treatment demonstrated that clonal evolution can be monitored over time. In conclusion, we have shown that detection of GIST-specific mutations in plasma is particularly feasible for patients with high tumor burden. In such cases, we have demonstrated that mutational analysis by use of liquid biopsies can capture the molecular heterogeneity of the whole tumor, and may guide treatment decisions during progression. Mol Cancer Ther; 17(11); 2473-80. ©2018 AACR.

Adipocyte differentiation of human bone marrow-derived stromal cells is modulated by microRNA-155, microRNA-221, and microRNA-222.

Human mesenchymal stromal cells (hMSCs) are capable of limited self-renewal and multilineage differentiation in vitro. Several studies have demonstrated that microRNAs (miRNAs, miRs), post-transcriptional modifiers of mRNA stability and protein translation, play crucial roles in the regulation of these complex processes. To gain knowledge regarding the role of miRNAs in human adipocyte differentiation, we examined the miRNA expression profile of the immortalized human bone marrow-derived stromal cell line hMSC-Tert20. Such a model system has the advantage of a reproducible and consistent phenotype while maintaining important properties of the primary donor cells, including the potential to differentiate to adipocytes, osteoblasts, and chondrocytes. We identified 12 miRNAs that were differentially expressed during adipogenesis, of which several have been previously shown to play important roles in adipocyte biology. Among these, the expression of miRNA-155, miRNA-221, and miRNA-222 decreased during the adipogenic program of both immortalized and primary hMSCs, suggesting that they act as negative regulators of differentiation. Interestingly, ectopic expression of the miRNAs significantly inhibited adipogenesis and repressed induction of the master regulators PPARγ and CCAAT/enhancer-binding protein alpha. Our study provides the first experimental evidence that miRNA-155, miRNA-221, and miRNA-222 have an important function in human adipocyte differentiation, and that their downregulation is necessary to relieve the repression of genes crucial for this process.

Modulation of the osteosarcoma expression phenotype by microRNAs.

BACKGROUND: Osteosarcomas are the most common primary malignant tumors of bone and show multiple and complex genomic aberrations. miRNAs are non-coding RNAs capable of regulating gene expression at the post transcriptional level, and miRNAs and their target genes may represent novel therapeutic targets or biomarkers for osteosarcoma. In order to investigate the involvement of miRNAs in osteosarcoma development, global microarray analyses of a panel of 19 human osteosarcoma cell lines was performed. PRINCIPAL FINDINGS: We identified 177 miRNAs that were differentially expressed in osteosarcoma cell lines relative to normal bone. Among these, miR-126/miR-126*, miR-142-3p, miR-150, miR-223, miR-486-5p and members of the miR-1/miR-133a, miR-144/miR-451, miR-195/miR-497 and miR-206/miR-133b clusters were found to be downregulated in osteosarcoma cell lines. All miRNAs in the paralogous clusters miR-17-92, miR-106b-25 and miR-106a-92 were overexpressed. Furthermore, the upregulated miRNAs included miR-9/miR-9*, miR-21*, miR-31/miR-31*, miR-196a/miR-196b, miR-374a and members of the miR-29 and miR-130/301 families. The most interesting inversely correlated miRNA/mRNA pairs in osteosarcoma cell lines included miR-9/TGFBR2 and miR-29/p85α regulatory subunit of PI3K. PTEN mRNA correlated inversely with miR-92a and members of the miR-17 and miR-130/301 families. Expression profiles of selected miRNAs were confirmed in clinical samples. A set of miRNAs, miR-1, miR-18a, miR-18b, miR-19b, miR-31, miR-126, miR-142-3p, miR-133b, miR-144, miR-195, miR-223, miR-451 and miR-497 was identified with an intermediate expression level in osteosarcoma clinical samples compared to osteoblasts and bone, which may reflect the differentiation level of osteosarcoma relative to the undifferentiated osteoblast and fully differentiated normal bone. SIGNIFICANCE: This study provides an integrated analysis of miRNA and mRNA in osteosarcoma, and gives new insight into the complex genetic mechanisms of osteosarcoma development and progression.

Integrative analysis reveals relationships of genetic and epigenetic alterations in osteosarcoma.

BACKGROUND: Osteosarcomas are the most common non-haematological primary malignant tumours of bone, and all conventional osteosarcomas are high-grade tumours showing complex genomic aberrations. We have integrated genome-wide genetic and epigenetic profiles from the EuroBoNeT panel of 19 human osteosarcoma cell lines based on microarray technologies. PRINCIPAL FINDINGS: The cell lines showed complex patterns of DNA copy number changes, where genomic copy number gains were significantly associated with gene-rich regions and losses with gene-poor regions. By integrating the datasets, 350 genes were identified as having two types of aberrations (gain/over-expression, hypo-methylation/over-expression, loss/under-expression or hyper-methylation/under-expression) using a recurrence threshold of 6/19 (>30%) cell lines. The genes showed in general alterations in either DNA copy number or DNA methylation, both within individual samples and across the sample panel. These 350 genes are involved in embryonic skeletal system development and morphogenesis, as well as remodelling of extracellular matrix. The aberrations of three selected genes, CXCL5, DLX5 and RUNX2, were validated in five cell lines and five tumour samples using PCR techniques. Several genes were hyper-methylated and under-expressed compared to normal osteoblasts, and expression could be reactivated by demethylation using 5-Aza-2'-deoxycytidine treatment for four genes tested; AKAP12, CXCL5, EFEMP1 and IL11RA. Globally, there was as expected a significant positive association between gain and over-expression, loss and under-expression as well as hyper-methylation and under-expression, but gain was also associated with hyper-methylation and under-expression, suggesting that hyper-methylation may oppose the effects of increased copy number for detrimental genes. CONCLUSIONS: Integrative analysis of genome-wide genetic and epigenetic alterations identified dependencies and relationships between DNA copy number, DNA methylation and mRNA expression in osteosarcomas, contributing to better understanding of osteosarcoma biology.

mRNA expression profiles of primary high-grade central osteosarcoma are preserved in cell lines and xenografts.

BACKGROUND: Conventional high-grade osteosarcoma is a primary malignant bone tumor, which is most prevalent in adolescence. Survival rates of osteosarcoma patients have not improved significantly in the last 25 years. Aiming to increase this survival rate, a variety of model systems are used to study osteosarcomagenesis and to test new therapeutic agents. Such model systems are typically generated from an osteosarcoma primary tumor, but undergo many changes due to culturing or interactions with a different host species, which may result in differences in gene expression between primary tumor cells, and tumor cells from the model system. We aimed to investigate whether gene expression profiles of osteosarcoma cell lines and xenografts are still comparable to those of the primary tumor. METHODS: We performed genome-wide mRNA expression profiling on osteosarcoma biopsies (n = 76), cell lines (n = 13), and xenografts (n = 18). Osteosarcoma can be subdivided into several histological subtypes, of which osteoblastic, chondroblastic, and fibroblastic osteosarcoma are the most frequent ones. Using nearest shrunken centroids classification, we generated an expression signature that can predict the histological subtype of osteosarcoma biopsies. RESULTS: The expression signature, which consisted of 24 probes encoding for 22 genes, predicted the histological subtype of osteosarcoma biopsies with a misclassification error of 15%. Histological subtypes of the two osteosarcoma model systems, i.e. osteosarcoma cell lines and xenografts, were predicted with similar misclassification error rates (15% and 11%, respectively). CONCLUSIONS: Based on the preservation of mRNA expression profiles that are characteristic for the histological subtype we propose that these model systems are representative for the primary tumor from which they are derived.