Comparing the epigenetic landscape in myonuclei purified with a PCM1 antibody from a fast/glycolytic and a slow/oxidative muscle.

Muscle cells have different phenotypes adapted to different usage, and can be grossly divided into fast/glycolytic and slow/oxidative types. While most muscles contain a mixture of such fiber types, we aimed at providing a genome-wide analysis of the epigenetic landscape by ChIP-Seq in two muscle extremes, the fast/glycolytic extensor digitorum longus (EDL) and slow/oxidative soleus muscles. Muscle is a heterogeneous tissue where up to 60% of the nuclei can be of a different origin. Since cellular homogeneity is critical in epigenome-wide association studies we developed a new method for purifying skeletal muscle nuclei from whole tissue, based on the nuclear envelope protein Pericentriolar material 1 (PCM1) being a specific marker for myonuclei. Using antibody labelling and a magnetic-assisted sorting approach, we were able to sort out myonuclei with 95% purity in muscles from mice, rats and humans. The sorting eliminated influence from the other cell types in the tissue and improved the myo-specific signal. A genome-wide comparison of the epigenetic landscape in EDL and soleus reflected the differences in the functional properties of the two muscles, and revealed distinct regulatory programs involving distal enhancers, including a glycolytic super-enhancer in the EDL. The two muscles were also regulated by different sets of transcription factors; e.g. in soleus, binding sites for MEF2C, NFATC2 and PPARA were enriched, while in EDL MYOD1 and SIX1 binding sites were found to be overrepresented. In addition, more novel transcription factors for muscle regulation such as members of the MAF family, ZFX and ZBTB14 were identified.

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.

Strategies to inhibit FGFR4 V550L-driven rhabdomyosarcoma.

BACKGROUND: Rhabdomyosarcoma (RMS) is a paediatric cancer driven either by fusion proteins (e.g., PAX3-FOXO1) or by mutations in key signalling molecules (e.g., RAS or FGFR4). Despite the latter providing opportunities for precision medicine approaches in RMS, there are currently no such treatments implemented in the clinic. METHODS: We evaluated biologic properties and targeting strategies for the FGFR4 V550L activating mutation in RMS559 cells, which have a high allelic fraction of this mutation and are oncogenically dependent on FGFR4 signalling. Signalling and trafficking of FGFR4 V550L were characterised by confocal microscopy and proteomics. Drug effects were determined by live-cell imaging, MTS assay, and in a mouse model. RESULTS: Among recently developed FGFR4-specific inhibitors, FGF401 inhibited FGFR4 V550L-dependent signalling and cell proliferation at low nanomolar concentrations. Two other FGFR4 inhibitors, BLU9931 and H3B6527, lacked potent activity against FGFR4 V550L. Alternate targeting strategies were identified by RMS559 phosphoproteomic analyses, demonstrating that RAS/MAPK and PI3K/AKT are essential druggable pathways downstream of FGFR4 V550L. Furthermore, we found that FGFR4 V550L is HSP90-dependent, and HSP90 inhibitors efficiently impeded RMS559 proliferation. In a RMS559 mouse xenograft model, the pan-FGFR inhibitor, LY2874455, did not efficiently inhibit growth, whereas FGF401 potently abrogated growth. CONCLUSIONS: Our results pave the way for precision medicine approaches against FGFR4 V550L-driven RMS.

Validation of a novel risk score to predict early and late recurrence in solitary fibrous tumour.

BACKGROUND: Current risk models in solitary fibrous tumour (SFT) were developed using cohorts with short follow-up and cannot reliably identify low-risk patients. We recently developed a novel risk model (G-score) to account for both early and late recurrences. Here, we aimed to validate the G-score in a large international cohort with long-term follow-up. METHODS: Data were collected from nine sarcoma referral centres worldwide. Recurrence-free interval (RFi) was the primary endpoint. RESULTS: The cohort comprised 318 patients with localised extrameningeal SFTs. Disease recurrence occurred in 96 patients (33%). The estimated 5-year RFi rate was 72%, and the 10-year RFi rate was 52%. G-score precisely predicted recurrence risk with estimated 10-year RFi rate of 84% in low risk, 54% in intermediate risk and 36% in high risk (p < 0.001; C-index 0.691). The mDemicco (p < 0.001; C-index 0.749) and SalasOS (p < 0.001; C-index 0.674) models also predicted RFi but identified low-risk patients less accurate with 10-year RFi rates of 72% and 70%, respectively. CONCLUSIONS: G-score is a highly significant predictor of early and late recurrence in SFT and is superior to other models to predict patients at low risk of relapse. A less intensive follow-up schedule could be considered for patients at low recurrence risk according to G-score.

Pembrolizumab in advanced osteosarcoma: results of a single-arm, open-label, phase 2 trial.

AIM: To evaluate the activity and safety of the PD-1 antibody pembrolizumab in adult patients with advanced osteosarcoma. MATERIAL AND METHODS: The study was a single-arm, open-label, phase 2 trial in patients with unresectable, relapsed osteosarcoma. The primary endpoint was clinical benefit rate (CBR) at 18 weeks of treatment, defined as complete response, partial response, or stable disease using RECIST v1.1. The trial had a Simon´s two-stage design, and ≥ 3 of 12 patients with clinical benefit in stage 1 were required to proceed to stage 2. The trial is registered with ClinicalTrials.gov, number NCT03013127. NanoString analysis was performed to explore tumor gene expression signatures and pathways. RESULTS: Twelve patients were enrolled and received study treatment. No patients had clinical benefit at 18 weeks of treatment, and patient enrollment was stopped after completion of stage 1. Estimated median progression-free survival was 1.7 months (95% CI 1.2-2.2). At time of data cut-off, 11 patients were deceased due to osteosarcoma. Median overall survival was 6.6 months (95% CI 3.8-9.3). No treatment-related deaths or drug-related grade 3 or 4 adverse events were observed. PD-L1 expression was positive in one of 11 evaluable tumor samples, and the positive sample was from a patient with a mixed treatment response. CONCLUSION: In this phase 2 study in advanced osteosarcoma, pembrolizumab was well-tolerated but did not show clinically significant antitumor activity. Future trials with immunomodulatory agents in osteosarcoma should explore combination strategies in patients selected based on molecular profiles associated with response.

Integrating Anatomical, Molecular and Clinical Risk Factors in Gastrointestinal Stromal Tumor of the Stomach.

BACKGROUND: Adjuvant imatinib for 3 years is recommended to patients with high-risk gastrointestinal stromal tumor (GIST). Risk stratification is inaccurate, and risk assessments are further complicated by the increased use of neoadjuvant treatment. Anatomical criteria for prognostication have not been investigated. METHODS: Clinical, molecular, and anatomical variables were retrospectively studied in a population-based cohort of 295 patients with gastric GIST resected between 2000 and 2018. Gastric subsite was divided into the upper, middle, and lower thirds. Growth pattern was classified as luminal, exophytic, or transmural based on imaging and surgical reports. RESULTS: Of 113 tumors in the upper third of the stomach, 103 (91.2%) were KIT mutated, 7 (6.2%) were PDGFRA mutated, and 104 (92.0%) harbored genotypes sensitive to imatinib. Transmural tumors were strongly associated with a high mitotic index. Five-year recurrence-free survival (RFS) was 71% for patients with transmural tumors versus 96% with luminal or exophytic tumors (hazard ratio [HR] 8.45, 95% confidence interval [CI] 3.69-19.36; p < 0.001), and, in high-risk patients, 5-year RFS was 46% for patients with transmural tumors versus 83% with luminal or exophytic tumors (HR 4.47, 95% CI 1.71-11.66; p = 0.001). Among 134 patients with tumors > 5 cm, there were 29 recurrences. Only five patients with exophytic or luminal tumors had recurrent disease, of whom four had tumor rupture. Five-year RFS for patients with exophytic/luminal tumors >5 cm without rupture was 98%. CONCLUSIONS: In the upper third, over 90% of tumors were sensitive to imatinib. Patients with exophytic or luminal tumors without rupture, irrespective of size, had an excellent prognosis and may not benefit from adjuvant therapy.

Sumoylation of Rap1 mediates the recruitment of TFIID to promote transcription of ribosomal protein genes.

Transcription factors are abundant Sumo targets, yet the global distribution of Sumo along the chromatin and its physiological relevance in transcription are poorly understood. Using Saccharomyces cerevisiae, we determined the genome-wide localization of Sumo along the chromatin. We discovered that Sumo-enriched genes are almost exclusively involved in translation, such as tRNA genes and ribosomal protein genes (RPGs). Genome-wide expression analysis showed that Sumo positively regulates their transcription. We also discovered that the Sumo consensus motif at RPG promoters is identical to the DNA binding motif of the transcription factor Rap1. We demonstrate that Rap1 is a molecular target of Sumo and that sumoylation of Rap1 is important for cell viability. Furthermore, Rap1 sumoylation promotes recruitment of the basal transcription machinery, and sumoylation of Rap1 cooperates with the target of rapamycin kinase complex 1 (TORC1) pathway to promote RPG transcription. Strikingly, our data reveal that sumoylation of Rap1 functions in a homeostatic feedback loop that sustains RPG transcription during translational stress. Taken together, Sumo regulates the cellular translational capacity by promoting transcription of tRNA genes and RPGs.

Preclinical Evaluation of Vemurafenib as Therapy for BRAFV600E Mutated Sarcomas.

The BRAFV600E mutation, which in melanoma is targetable with vemurafenib, is also found in sarcomas and we here evaluate the therapeutic potential in sarcoma cell lines. METHODS: Four sarcoma cell lines harboring the BRAFV600E mutation, representing liposarcomas (SA-4 and SW872), Ewing sarcoma (A673) and atypical synovial sarcoma (SW982), were treated with vemurafenib and the effects on cell growth, apoptosis, cell cycle progression and cell signaling were determined. RESULTS: Vemurafenib induced a strong cytostatic effect in SA-4 cells, mainly due to cell cycle arrest, whereas only moderate levels of apoptosis were observed. However, a high dose was required compared to BRAFV600E mutated melanoma cells, and removal of vemurafenib demonstrated that the continuous presence of drug was required for sustained growth inhibition. A limited growth inhibition was observed in the other three cell lines. Protein analyses demonstrated reduced phosphorylation of ERK during treatment with vemurafenib in all the four sarcoma cell lines confirming that the MAPK pathway is active in these cell lines, and that the pathway can be inhibited by vemurafenib, but also that these cells can proliferate despite this. CONCLUSIONS: These findings indicate that vemurafenib alone would not be an efficient therapy against BRAFV600E mutated sarcomas. However, further investigations of combination with other drugs are warranted.

Multi-omics of 34 colorectal cancer cell lines - a resource for biomedical studies.

BACKGROUND: Colorectal cancer (CRC) cell lines are widely used pre-clinical model systems. Comprehensive insights into their molecular characteristics may improve model selection for biomedical studies. METHODS: We have performed DNA, RNA and protein profiling of 34 cell lines, including (i) targeted deep sequencing (n = 612 genes) to detect single nucleotide variants and insertions/deletions; (ii) high resolution DNA copy number profiling; (iii) gene expression profiling at exon resolution; (iv) small RNA expression profiling by deep sequencing; and (v) protein expression analysis (n = 297 proteins) by reverse phase protein microarrays. RESULTS: The cell lines were stratified according to the key molecular subtypes of CRC and data were integrated at two or more levels by computational analyses. We confirm that the frequencies and patterns of DNA aberrations are associated with genomic instability phenotypes and that the cell lines recapitulate the genomic profiles of primary carcinomas. Intrinsic expression subgroups are distinct from genomic subtypes, but consistent at the gene-, microRNA- and protein-level and dominated by two distinct clusters; colon-like cell lines characterized by expression of gastro-intestinal differentiation markers and undifferentiated cell lines showing upregulation of epithelial-mesenchymal transition and TGFß signatures. This sample split was concordant with the gene expression-based consensus molecular subtypes of primary tumors. Approximately » of the genes had consistent regulation at the DNA copy number and gene expression level, while expression of gene-protein pairs in general was strongly correlated. Consistent high-level DNA copy number amplification and outlier gene- and protein- expression was found for several oncogenes in individual cell lines, including MYC and ERBB2. CONCLUSIONS: This study expands the view of CRC cell lines as accurate molecular models of primary carcinomas, and we present integrated multi-level molecular data of 34 widely used cell lines in easily accessible formats, providing a resource for preclinical studies in CRC.

Monitoring multiple myeloma by quantification of recurrent mutations in serum.

Circulating tumor DNA is a promising biomarker to monitor tumor load and genome alterations. We explored the presence of circulating tumor DNA in multiple myeloma patients and its relation to disease activity during long-term follow-up. We used digital droplet polymerase chain reaction analysis to monitor recurrent mutations, mainly in mitogen activated protein kinase pathway genes NRAS, KRAS and BRAF Mutations were identified by next-generation sequencing or polymerase chain reaction analysis of bone marrow plasma cells, and their presence analyzed in 251 archived serum samples obtained from 20 patients during a period of up to 7 years. In 17 of 18 patients, mutations identified in bone marrow during active disease were also found in a time-matched serum sample. The concentration of mutated alleles in serum correlated with the fraction in bone marrow plasma cells (r=0.507, n=34, P<0.002). There was a striking covariation between circulating mutation levels and M protein in ten out of 11 patients with sequential samples. When relapse evaluation by circulating tumor DNA and M protein could be directly compared, the circulating tumor DNA showed relapse earlier in two patients (3 and 9 months), later in one patient (4 months) and in three patients there was no difference. In three patients with transformation to aggressive disease, the concentrations of mutations in serum increased up to 400 times, an increase that was not seen for the M protein. In conclusion, circulating tumor DNA in myeloma is a multi-faceted biomarker reflecting mutated cells, total tumor mass and transformation to a more aggressive disease. Its properties are both similar and complementary to M protein.

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.

VOLIN and KJON-Two novel hyperdiploid myeloma cell lines.

Multiple myeloma can be divided into two distinct genetic subgroups: hyperdiploid (HRD) or nonhyperdiploid (NHRD) myeloma. Myeloma cell lines are important tools to study myeloma cell biology and are commonly used for preclinical screening and testing of new drugs. With few exceptions human myeloma cell lines are derived from NHRD patients, even though about half of the patients have HRD myeloma. Thus, there is a need for cell lines of HRD origin to enable more representative preclinical studies. Here, we present two novel myeloma cell lines, VOLIN and KJON. Both of them were derived from patients with HRD disease and shared the same genotype as their corresponding primary tumors. The cell lines' chromosomal content, genetic aberrations, gene expression, immunophenotype as well as some of their growth characteristics are described. Neither of the cell lines was found to harbor immunoglobulin heavy chain translocations. The VOLIN cell line was established from a bone marrow aspirate and KJON from peripheral blood. We propose that these unique cell lines may be used as tools to increase our understanding of myeloma cell biology. © 2016 Wiley Periodicals, Inc.

The regulatory landscape of osteogenic differentiation.

Differentiation of osteoblasts from mesenchymal stem cells (MSCs) is an integral part of bone development and homeostasis, and may when improperly regulated cause disease such as bone cancer or osteoporosis. Using unbiased high-throughput methods we here characterize the landscape of global changes in gene expression, histone modifications, and DNA methylation upon differentiation of human MSCs to the osteogenic lineage. Furthermore, we provide a first genome-wide characterization of DNA binding sites of the bone master regulatory transcription factor Runt-related transcription factor 2 (RUNX2) in human osteoblasts, revealing target genes associated with regulation of proliferation, migration, apoptosis, and with a significant overlap with p53 regulated genes. These findings expand on emerging evidence of a role for RUNX2 in cancer, including bone metastases, and the p53 regulatory network. We further demonstrate that RUNX2 binds to distant regulatory elements, promoters, and with high frequency to gene 3' ends. Finally, we identify TEAD2 and GTF2I as novel regulators of osteogenesis.

Cdc28 kinase activity regulates the basal transcription machinery at a subset of genes.

The cyclin-dependent kinase Cdc28 is the master regulator of the cell cycle in Saccharomyces cerevisiae. Cdc28 initiates the cell cycle by activating cell-cycle-specific transcription factors that switch on a transcriptional program during late G1 phase. Cdc28 also has a cell-cycle-independent, direct function in regulating basal transcription, which does not require its catalytic activity. However, the exact role of Cdc28 in basal transcription remains poorly understood, and a function for its kinase activity has not been fully explored. Here we show that the catalytic activity of Cdc28 is important for basal transcription. Using a chemical-genetic screen for mutants that specifically require the kinase activity of Cdc28 for viability, we identified a plethora of basal transcription factors. In particular, CDC28 interacts genetically with genes encoding kinases that phosphorylate the C-terminal domain of RNA polymerase II, such as KIN28. ChIP followed by high-throughput sequencing (ChIP-seq) revealed that Cdc28 localizes to at least 200 genes, primarily with functions in cellular homeostasis, such as the plasma membrane proton pump PMA1. Transcription of PMA1 peaks early in the cell cycle, even though the promoter sequences of PMA1 (as well as the other Cdc28-enriched ORFs) lack cell-cycle elements, and PMA1 does not recruit Swi4/6-dependent cell-cycle box-binding factor/MluI cell-cycle box binding factor complexes. Finally, we found that recruitment of Cdc28 and Kin28 to PMA1 is mutually dependent and that the activity of both kinases is required for full phosphorylation of C-terminal domain-Ser5, for efficient transcription, and for mRNA capping. Our results reveal a mechanism of cell-cycle-dependent regulation of basal transcription.

Performance comparison of four exome capture systems for deep sequencing.

BACKGROUND: Recent developments in deep (next-generation) sequencing technologies are significantly impacting medical research. The global analysis of protein coding regions in genomes of interest by whole exome sequencing is a widely used application. Many technologies for exome capture are commercially available; here we compare the performance of four of them: NimbleGen's SeqCap EZ v3.0, Agilent's SureSelect v4.0, Illumina's TruSeq Exome, and Illumina's Nextera Exome, all applied to the same human tumor DNA sample. RESULTS: Each capture technology was evaluated for its coverage of different exome databases, target coverage efficiency, GC bias, sensitivity in single nucleotide variant detection, sensitivity in small indel detection, and technical reproducibility. In general, all technologies performed well; however, our data demonstrated small, but consistent differences between the four capture technologies. Illumina technologies cover more bases in coding and untranslated regions. Furthermore, whereas most of the technologies provide reduced coverage in regions with low or high GC content, the Nextera technology tends to bias towards target regions with high GC content. CONCLUSIONS: We show key differences in performance between the four technologies. Our data should help researchers who are planning exome sequencing to select appropriate exome capture technology for their particular application.

Reexpression of LSAMP inhibits tumor growth in a preclinical osteosarcoma model.

BACKGROUND: Osteosarcomas are the most common primary malignant tumors of bone, showing complex chromosomal rearrangements with multiple gains and losses. A frequent deletion within the chromosomal region 3q13.31 has been identified by us and others, and is mainly reported to be present in osteosarcomas. The purpose of the study was to further characterize the frequency and the extent of the deletion in an extended panel of osteosarcoma samples, and the expression level of the affected genes within the region. We have identified LSAMP as the target gene for the deletion, and have studied the functional implications of LSAMP-reexpression. METHODS: LSAMP copy number, expression level and protein level were investigated by quantitative PCR and western blotting in an osteosarcoma panel. The expression of LSAMP was restored in an osteosarcoma cell line, and differences in proliferation rate, tumor formation, gene expression, migration rate, differentiation capabilities, cell cycle distribution and apoptosis were investigated by metabolic dyes, tumor formation in vivo, gene expression profiling, time-lapse photography, differentiation techniques and flow cytometry, respectively. RESULTS: We found reduced copy number of LSAMP in 45/76 osteosarcoma samples, reduced expression level in 25/42 samples and protein expression in 9/42 samples. By restoring the expression of LSAMP in a cell line with a homozygous deletion of the gene, the proliferation rate in vitro was significantly reduced and tumor growth in vivo was significantly delayed. In response to reexpression of LSAMP, mRNA expression profiling revealed consistent upregulation of the genes hairy and enhancer of split 1 (HES1), cancer/testis antigen 2 (CTAG2) and kruppel-like factor 10 (KLF10). CONCLUSIONS: The high frequency and the specificity of the deletion indicate that it is important for the development of osteosarcomas. The deletion targets the tumor suppressor LSAMP, and based on the functional evidence, the tumor suppressor function of LSAMP is most likely exerted by reducing the proliferation rate of the tumor cells, possibly by indirectly upregulating one or more of the genes HES1, CTAG2 or KLF10. To our knowledge, this study describes novel functions of LSAMP, a first step to understanding the functional role of this specific deletion in osteosarcomas.

Evolutionary mode and timing of dissemination of high-grade serous carcinomas.

Dissemination within the peritoneal cavity is a main determinant of poor patient outcomes from high-grade serous carcinomas (HGSCs). The dissemination process is poorly understood from a cancer evolutionary perspective. We reconstructed the evolutionary trajectories across a median of 5 tumor sites and regions from each of 23 patients based on deep whole-exome sequencing. Polyclonal cancer origin was detected in 1 patient. Ovarian tumors had more complex subclonal architectures than other intraperitoneal tumors in each patient, which indicated that tumors developed earlier in the ovaries. Three common modes of dissemination were identified, including monoclonal or polyclonal dissemination of monophyletic (linear) or polyphyletic (branched) subclones. Mutation profiles of initial or disseminated clones varied greatly among cancers, but recurrent mutations were found in 7 cancer-critical genes, including TP53, BRCA1, BRCA2, and DNMT3A, and in the PI3K/AKT1 pathway. Disseminated clones developed late in the evolutionary trajectory models of most cancers, in particular in cancers with DNA damage repair deficiency. Polyclonal dissemination was predicted to occur predominantly as a single and rapid wave, but chemotherapy exposure was associated with higher genomic diversity of disseminated clones. In conclusion, we described three common evolutionary dissemination modes across HGSCs and proposed factors associated with dissemination diversity.

Systematic analysis of RNA-binding proteins identifies targetable therapeutic vulnerabilities in osteosarcoma.

Osteosarcoma is the most common primary malignant bone tumor with a strong tendency to metastasize, limiting the prognosis of affected patients. Genomic, epigenomic and transcriptomic analyses have demonstrated the exquisite molecular complexity of this tumor, but have not sufficiently defined the underlying mechanisms or identified promising therapeutic targets. To systematically explore RNA-protein interactions relevant to OS, we define the RNA interactomes together with the full proteome and the transcriptome of cells from five malignant bone tumors (four osteosarcomata and one malignant giant cell tumor of the bone) and from normal mesenchymal stem cells and osteoblasts. These analyses uncover both systematic changes of the RNA-binding activities of defined RNA-binding proteins common to all osteosarcomata and individual alterations that are observed in only a subset of tumors. Functional analyses reveal a particular vulnerability of these tumors to translation inhibition and a positive feedback loop involving the RBP IGF2BP3 and the transcription factor Myc which affects cellular translation and OS cell viability. Our results thus provide insight into potentially clinically relevant RNA-binding protein-dependent mechanisms of osteosarcoma.

Genome-wide map of quantified epigenetic changes during in vitro chondrogenic differentiation of primary human mesenchymal stem cells.

BACKGROUND: For safe clinical application of engineered cartilage made from mesenchymal stem cells (MSCs), molecular mechanisms for chondrogenic differentiation must be known in detail. Changes in gene expression and extracellular matrix synthesis have been extensively studied, but the epigenomic modifications underlying these changes have not been described. To this end we performed whole-genome chromatin immunoprecipitation and deep sequencing to quantify six histone modifications, reduced representation bisulphite sequencing to quantify DNA methylation and mRNA microarrays to quantify gene expression before and after 7 days of chondrogenic differentiation of MSCs in an alginate scaffold. To add to the clinical relevance of our observations, the study is based on primary bone marrow-derived MSCs from four donors, allowing us to investigate inter-individual variations. RESULTS: We see two levels of relationship between epigenetic marking and gene expression. First, a large number of genes ontogenetically linked to MSC properties and the musculoskeletal system are epigenetically prepatterned by moderate changes in H3K4me3 and H3K9ac near transcription start sites. Most of these genes remain transcriptionally unaltered. Second, transcriptionally upregulated genes, more closely associated with chondrogenesis, are marked by H3K36me3 in gene bodies, highly increased H3K4me3 and H3K9ac on promoters and 5' end of genes, and increased H3K27ac and H3K4me1 marking in at least one enhancer region per upregulated gene. Within the 7-day time frame, changes in promoter DNA methylation do not correlate significantly with changes in gene expression. Inter-donor variability analysis shows high level of similarity between the donors for this data set. CONCLUSIONS: Histone modifications, rather than DNA methylation, provide the primary epigenetic control of early differentiation of MSCs towards the chondrogenic lineage.