Imprinting defects at human 14q32 locus alters gene expression and is associated with the pathobiology of osteosarcoma

Biological Sample Collection to Advance Research and Treatment: A Fight Osteosarcoma Through European Research and Euro Ewing Consortium Statement

Osteosarcoma and Ewing sarcoma are bone tumors mostly diagnosed in children, adolescents, and young adults. Despite multimodal therapy, morbidity is high and survival rates remain low, especially in the metastatic disease setting. Trials investigating targeted therapies and immunotherapies have not been groundbreaking. Better understanding of biological subgroups, the role of the tumor immune microenvironment, factors that promote metastasis, and clinical biomarkers of prognosis and drug response are required to make progress. A prerequisite to achieve desired success is a thorough, systematic, and clinically linked biological analysis of patient samples, but disease rarity and tissue processing challenges such as logistics and infrastructure have contributed to a lack of relevant samples for clinical care and research. There is a need for a Europe-wide framework to be implemented for the adequate and minimal sampling, processing, storage, and analysis of patient samples. Two international panels of scientists, clinicians, and patient and parent advocates have formed the Fight Osteosarcoma Through European Research consortium and the Euro Ewing Consortium. The consortia shared their expertise and institutional practices to formulate new guidelines. We report new reference standards for adequate and minimally required sampling (time points, diagnostic samples, and liquid biopsy tubes), handling, and biobanking to enable advanced biological studies in bone sarcoma. We describe standards for analysis and annotation to drive collaboration and data harmonization with practical, legal, and ethical considerations. This position paper provides comprehensive guidelines that should become the new standards of care that will accelerate scientific progress, promote collaboration, and improve outcomes.

Bioinformatics analysis and validation of mesenchymal stem cells related gene MT1G in osteosarcoma

BACKGROUND

Osteosarcoma (OS) is a primary malignant bone tumor arising from mesenchymal cells. The standard clinical treatment for OS involves extensive tumor resection combined with neoadjuvant chemotherapy or radiotherapy. OS's invasiveness, lung metastasis, and drug resistance contribute to a low cure rate and poor prognosis with this treatment. Metallothionein 1G (MT1G), observed in various cancers, may serve as a potential therapeutic target for OS.

METHODS

OS samples in GSE33382 and TARGET datasets were selected as the test cohorts. As the external validation cohort, 13 OS tissues and 13 adjacent cancerous tissues from The Second Affiliated Hospital of Nanchang University were collected. Patients with OS were divided into high and low MT1G mRNA-expression groups; differentially expressed genes (DEGs) were identified as MT1G-related genes. The biological function of MT1G was annotated using Kyoto Encyclopedia of Genes and Genomes (KEGG), Gene Ontology (GO) and gene set enrichment analysis (GSEA). Gene expression correlation analysis and competing endogenous RNA (ceRNA) regulatory network construction were used to determine potential biological regulatory relationships of DEGs. Survival analysis assessed the prognostic value of MT1G.

RESULTS

MT1G expression increased in OS samples and presented higher in metastatic OS compared with non-metastatic OS. Functional analyses indicated that MT1G was mainly associated with spliceosome. A ceRNA network with DEGs was constructed. MT1G is an effective biomarker predicting survival and correlated with increased recurrence rates and poorer survival.

CONCLUSIONS

This research identified MT1G as a potential biomarker for OS prognosis, highlighting its potential as a therapy target.

Potential biomarkers for the early detection of bone metastases

The clinical manifestations of bone metastases are diversified while many sites remain asymptomatic at early stage. As the early diagnosis method is not perfect and the early symptoms of tumor bone metastasis are not typical, bone metastasis is not easy to be detected. Therefore, the search for bone metastasis-related markers is effective for timely detection of tumor bone metastases and the development of drugs to inhibit bone metastases. As a result, bone metastases can only be diagnosed when symptoms are found, increasing the risk of developing skeletal-related event (SREs), which significantly impairs the patient's quality of life. Therefore, the early diagnosis of bone metastases is of great importance for the treatment and prognosis of cancer patients. Changes of bone metabolism indexes appear earlier in bone metastases, but the traditional biochemical indexes of bone metabolism lack of specificity and could be interfered by many factors, which limits their application in the study of bone metastases. Some new biomarkers of bone metastases have good diagnostic value, such as proteins, ncRNAs, circulating tumor cells (CTCs). Therefore, this study mainly reviewed the initial diagnostic biomarkers of bone metastases which were expected to provide references for the early detection of bone metastases.

Deregulated expression of the 14q32 miRNA cluster in clear cell renal cancer cells

Clear cell renal cell carcinomas (ccRCC) are characterized by arm-wide chromosomal alterations. Loss at 14q is associated with disease aggressiveness in ccRCC, which responds poorly to chemotherapeutics. The 14q locus contains one of the largest miRNA clusters in the human genome; however, little is known about the contribution of these miRNAs to ccRCC pathogenesis. In this regard, we investigated the expression pattern of selected miRNAs at the 14q32 locus in TCGA kidney tumors and in ccRCC cell lines. We demonstrated that the miRNA cluster is downregulated in ccRCC (and cell lines) as well as in papillary kidney tumors relative to normal kidney tissues (and primary renal proximal tubule epithelial (RPTEC) cells). We demonstrated that agents modulating expression of DNMT1 (e.g., 5-Aza-deoxycytidine) could modulate 14q32 miRNA expression in ccRCC cell lines. Lysophosphatidic acid (LPA, a lysophospholipid mediator elevated in ccRCC) not only increased labile iron content but also modulated expression of a 14q32 miRNA. Through an overexpression approach targeting a subset of 14q32 miRNAs (specifically at subcluster A: miR-431-5p, miR-432-5p, miR-127-3p, and miR-433-3p) in 769-P cells, we uncovered changes in cellular viability and claudin-1, a tight junction marker. A global proteomic approach was implemented using these miRNA overexpressing cell lines which uncovered ATXN2 as a highly downregulated target. Collectively, these findings support a contribution of miRNAs at 14q32 in ccRCC pathogenesis.

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.

YBX1-interacting small RNAs and RUNX2 can be blocked in primary bone cancer using CADD522

Primary bone cancer (PBC) comprises several subtypes each underpinned by distinctive genetic drivers. This driver diversity produces novel morphological features and clinical behaviour that serendipitously makes PBC an excellent metastasis model. Here, we report that some transfer RNA-derived small RNAs termed tRNA fragments (tRFs) perform as a constitutive tumour suppressor mechanism by blunting a potential pro-metastatic protein-RNA interaction. This mechanism is reduced in PBC progression with a gradual loss of tRNAGly^TCC cleavage into 5' end tRF-Gly^TCC when comparing low-grade, intermediate-grade and high-grade patient tumours. We detected recurrent activation of miR-140 leading to upregulated RUNX2 expression in high-grade patient tumours. Both tRF-Gly^TCC and RUNX2 share a sequence motif in their 3' ends that matches the YBX1 recognition site known to stabilise pro-metastatic mRNAs. Investigating some aspects of this interaction network, gain- and loss-of-function experiments using small RNA mimics and antisense LNAs, respectively, showed that ectopic tRF-Gly^TCC reduced RUNX2 expression and dispersed 3D micromass architecture in vitro. iCLIP sequencing revealed YBX1 physical binding to the 3' UTR of RUNX2. The interaction between YBX1, tRF-Gly^TCC and RUNX2 led to the development of the RUNX2 inhibitor CADD522 as a PBC treatment. CADD522 assessment in vitro revealed significant effects on PBC cell behaviour. In xenograft mouse models, CADD522 as a single agent without surgery significantly reduced tumour volume, increased overall and metastasis-free survival and reduced cancer-induced bone disease. Our results provide insight into PBC molecular abnormalities that have led to the identification of new targets and a new therapeutic.

The emerging role of non-coding RNAs in the epigenetic regulation of pediatric cancers

Epigenetics is a process that involves the regulation of gene expression without altering the sequence of DNA. Numerous studies have documented that epigenetic mechanisms play a critical role in cell growth, differentiation, and cancer over the past decade. The well-known epigenetic modifications are either on DNA or at the histone proteins. Although several studies have focused on regulating gene expression by non-coding RNAs, the current understanding of their biological functions in various human diseases, particularly in cancers, is inadequate. Only about two percent of DNA is involved in coding the protein-coding genes, and leaving the rest 98 percent is non-coding and the scientific community regarded as junk or noise with no known purpose. Most non-coding RNAs are derived from such junk DNA and are known to be involved in various signaling pathways involving cancer initiation, progression, and the development of therapy resistance in many human cancer types. Recent studies have suggested that non-coding RNAs, especially microRNAs, piwi-interactingRNAs, and long non-coding RNAs, play a significant role in controlling epigenetic mechanism(s), indicating the potential effect of epigenetic modulation of non-coding RNAs on cancer progression. In this review article, we briefly presented epigenetic marks' characteristics, crosstalk between epigenetic modifications and microRNAs, piwi-interactingRNAs, and long non-coding RNAs to uncover the effect on the phenotype of pediatric cancers. Further, current knowledge on understanding the RNA epigenetics will help design novel therapeutics that target epigenetic regulatory networks to benefit cancer patients in the clinic.

Comparative analysis of genome-wide DNA methylation identifies patterns that associate with conserved transcriptional programs in osteosarcoma

Osteosarcoma is an aggressive tumor of the bone that primarily affects young adults and adolescents. Osteosarcoma is characterized by genomic chaos and heterogeneity. While inactivation of tumor protein p53 (TP53) is nearly universal other high frequency mutations or structural variations have not been identified. Despite this genomic heterogeneity, key conserved transcriptional programs associated with survival have been identified across human, canine and induced murine osteosarcoma. The epigenomic landscape, including DNA methylation, plays a key role in establishing transcriptional programs in all cell types. The role of epigenetic dysregulation has been studied in a variety of cancers but has yet to be explored at scale in osteosarcoma. Here we examined genome-wide DNA methylation patterns in 24 human and 44 canine osteosarcoma samples identifying groups of highly correlated DNA methylation marks in human and canine osteosarcoma samples. We also link specific DNA methylation patterns to key transcriptional programs in both human and canine osteosarcoma. Building on previous work, we built a DNA methylation-based measure for the presence and abundance of various immune cell types in osteosarcoma. Finally, we determined that the underlying state of the tumor, and not changes in cell composition, were the main driver of differences in DNA methylation across the human and canine samples. SIGNIFICANCE: Genome wide comparison of DNA methylation patterns in osteosarcoma across two species lays the ground work for the exploration of DNA methylation programs that help establish conserved transcriptional programs in the context of varied mutational landscapes.

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.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08346-6.

Essential Role of the 14q32 Encoded miRNAs in Endocrine Tumors

BACKGROUND

The 14q32 cluster is among the largest polycistronic miRNA clusters. miRNAs encoded here have been implicated in tumorigenesis of multiple organs including endocrine glands.

METHODS

Critical review of miRNA studies performed in endocrine tumors have been performed. The potential relevance of 14q32 miRNAs through investigating their targets, and integrating the knowledge provided by literature data and bioinformatics predictions have been indicated.

RESULTS

Pituitary adenoma, papillary thyroid cancer and a particular subset of pheochromocytoma and adrenocortical cancer are characterized by the downregulation of miRNAs encoded by the 14q32 cluster. Pancreas neuroendocrine tumors, most of the adrenocortical cancer and medullary thyroid cancer are particularly distinct, as 14q32 miRNAs were overexpressed. In pheochromocytoma and growth-hormone producing pituitary adenoma, however, both increased and decreased expression of 14q32 miRNAs cluster members were observed. In the background of this phenomenon methodological, technical and biological factors are hypothesized and discussed. The functions of 14q32 miRNAs were also revealed by bioinformatics and literature data mining.

CONCLUSIONS

14q32 miRNAs have a significant role in the tumorigenesis of endocrine organs. Regarding their stable expression in the circulation of healthy individuals, further investigation of 14q32 miRNAs could provide a potential for use as biomarkers (diagnostic or prognostic) in endocrine neoplasms.

The EMT transcription factor ZEB1 blocks osteoblastic differentiation in bone development and osteosarcoma

Osteosarcoma is an often-fatal mesenchyme-derived malignancy in children and young adults. Overexpression of EMT-transcription factors (EMT-TFs) has been associated with poor clinical outcome. Here, we demonstrated that the EMT-TF ZEB1 is able to block osteoblastic differentiation in normal bone development as well as in osteosarcoma cells. Consequently, overexpression of ZEB1 in osteosarcoma characterizes poorly differentiated, highly metastatic subgroups and its depletion induces differentiation of osteosarcoma cells. Overexpression of ZEB1 in osteosarcoma is frequently associated with silencing of the imprinted DLK-DIO3 locus, which encodes for microRNAs targeting ZEB1. Epigenetic reactivation of this locus in osteosarcoma cells reduces ZEB1 expression, induces differentiation, and sensitizes to standard treatment, thus indicating therapeutic options for ZEB1-driven osteosarcomas. © 2021 The Authors. The Journal of Pathology published by John Wiley & Sons, Ltd. on behalf of The Pathological Society of Great Britain and Ireland.

DNA Methylation Controls Metastasis-Suppressive 14q32-Encoded miRNAs

Expression of 14q32-encoded miRNAs is a favorable prognostic factor in patients with metastatic cancer. In this study, we used genomic inhibition of DNA methylation through disruption of DNA methyltransferases DNMT1 and DNMT3B and pharmacologic inhibition with 5-Aza-2'-deoxycytidine (5-Aza-dC, decitabine) to demonstrate that DNA methylation predominantly regulates expression of metastasis-suppressive miRNAs in the 14q32 cluster. DNA demethylation facilitated CCCTC-binding factor (CTCF) recruitment to the maternally expressed gene 3 differentially methylated region (MEG3-DMR), which acts as a cis-regulatory element for 14q32 miRNA expression. 5-Aza-dC activated demethylation of the MEG3-DMR and expression of 14q32 miRNAs, which suppressed adhesion, invasion, and migration (AIM) properties of metastatic tumor cells. Cancer cells with MEG3-DMR hypomethylation exhibited constitutive expression of 14q32 miRNAs and resistance to 5-Aza-dC-induced suppression of AIM. Expression of methylation-dependent 14q32 miRNAs suppressed metastatic colonization in preclinical models of lung and liver metastasis and correlated with improved clinical outcomes in patients with metastatic cancer. These findings implicate epigenetic modification via DNA methylation in the regulation of metastatic propensity through miRNA networks and identify a previously unrecognized action of decitabine on the activation of metastasis-suppressive miRNAs. SIGNIFICANCE: This study investigates epigenetic regulation of metastasis-suppressive miRNAs and the effect on metastasis.

Variants in the 14q32 miRNA cluster are associated with osteosarcoma risk in the Spanish population

Association studies in osteosarcoma risk found significant results in intergenic regions, suggesting that regions which do not codify for proteins could play an important role. The deregulation of microRNAs (miRNAs) has been already associated with osteosarcoma. Consequently, genetic variants affecting miRNA function could be associated with risk. This study aimed to evaluate the involvement of all genetic variants in pre-miRNAs described so far in relationship to the risk of osteosarcoma. We analyzed a total of 213 genetic variants in 206 pre-miRNAs in two cohorts of osteosarcoma patients (n = 100) and their corresponding controls (n = 256) from Spanish and Slovenian populations, using Goldengate Veracode technology (Illumina). Four polymorphisms in pre-miRNAs at 14q32 miRNA cluster were associated with osteosarcoma risk in the Spanish population (rs12894467, rs61992671, rs58834075 and rs12879262). Pathway enrichment analysis including target genes of these miRNAs pointed out the WNT signaling pathways overrepresented. Moreover, different single nucleotide polymorphism (SNP) effects between the two populations included were observed, suggesting the existence of population differences. In conclusion, 14q32 miRNA cluster seems to be a hotspot for osteosarcoma susceptibility in the Spanish population, but not in the Slovenian, which supports the idea of the existence of population differences in developing this disease.

DLK1-DIO3 imprinted locus deregulation in development, respiratory disease, and cancer

INTRODUCTION

The imprinted DLK1-DIO3 locus at 14q32.1-32.31 holds biological significance in fetal development, whereby imprinting errors are causal to developmental disorders. Emerging evidence has implicated this locus in other diseases including cancer, highlighting the biological parallels between fetal organ and tumour development. Areas covered: Controlled regulation of gene expression from the imprinted DLK1-DIO3 locus at 14q32.1-32.31 is crucial for proper fetal development. Deregulation of locus gene expression due to imprinting errors has been mechanistically linked to the developmental disorders Kagami-Ogata Syndrome and Temple Syndrome. In adult tissues, deregulation of locus genes has been associated with multiple malignancies although the causal genetic mechanisms remain largely uncharacterised. Here, we summarize the genetic mechanisms underlying the developmental disorders that arise as a result of improper locus imprinting and the resulting developmental phenotypes, emphasizing both the coding and noncoding components of the locus. We further highlight biological parallels common to both fetal development and disease, with a specific focus on lung development, respiratory disease, and lung cancer. Expert commentary: Many commonalities between respiratory and developmental defects have emerged with respect to the 14q32 locus, emphasizing the importance of studying the effects of imprinting on gene regulation patterns at this locus in both biological settings.

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.

Electronic supplementary material

The online version of this article (doi:10.1186/s13045-017-0465-4) contains supplementary material, which is available to authorized users.