Epigenetic changes in osteosarcoma

Retraction

Bao G, Pan W, Huang J, Zhou T. K-RasG12V/T35S -ERK1/2 pathway regulates H2BS14ph through Mst1 to facilitate the advancement of breast cancer cells. BioFactors. 2023;49:202. https://doi.org/10.1002/biof.1589 This article, published online on 28 November 2019 in Wiley Online Library, has been retracted by agreement between the International Union of Biochemistry and Molecular Biology, the Editor in Chief (Dr. Angelo Azzi), and Wiley Periodicals LLC. The retraction has been agreed following an investigation based on allegations raised by a third party. Evidence for image manipulation was found in figures 1, 4, 5, and 6. As a result, the conclusions of this article are considered to be invalid.

Histone Deacetylase Inhibitor Trichostatin A Suppresses Cell Proliferation and Induces Apoptosis by Regulating the PI3K/AKT Signalling Pathway in Gastric Cancer Cells

BACKGROUND

Gastric cancer, a common malignant tumour worldwide, has a relatively poor prognosis and is a serious threat to human health. Histone Deacetylase Inhibitors (HDACi) are anticancer agents that are known to affect the cell growth of different cancer types. Trichostatin A (TSA) selectively inhibits the class I and II mammalian Histone Deacetylase (HDAC) family enzymes and regulates many cell processes. Still, the underlying mechanisms of HDACs are not fully understood in gastric cancer.

OBJECTIVE

This study aims to investigate the antitumor effect and the mechanism of growth modulation of gastric cancer cells by TSA.

METHODS

The cell proliferation of gastric cancer cells was measured by MTT and BrdU immunofluorescence assays. Soft agar assay was used to detect the colony formation ability of gastric cancer cells. Flow cytometry was used to examine cell cycle and apoptosis. Western blot was employed to detect protein expression of target factors.

RESULTS

TSA inhibits the proliferation of MKN-45 and SGC-7901 cells and leads to significant repression of colony number and size. Flow cytometry assays show TSA induces cell cycle arrest at G1 phase and apoptosis, and TSA effects the expression of related factors in the mitochondrial apoptotic signalling and cell cycle-related regulatory pathways. Furthermore, TSA increased histone H3K27 acetylation and downregulated the expression of PI3K and p-AKT.

CONCLUSION

Downregulating PI3K/AKT pathway activation is involved in TSA-mediated proliferation inhibition of gastric cancer.

The histone acetyltransferase HBO1 functions as a novel oncogenic gene in osteosarcoma

HBO1 (KAT7 or MYST2) is a histone acetyltransferase that acetylates H3 and H4 histones.

Methods: HBO1 expression was tested in human OS tissues and cells. Genetic strategies, including shRNA, CRISPR/Cas9 and overexpression constructs, were applied to exogenously alter HBO1 expression in OS cells. The HBO1 inhibitor WM-3835 was utilized to block HBO1 activation.

Results:HBO1 mRNA and protein expression is significantly elevated in OS tissues and cells. In established (MG63/U2OS lines) and primary human OS cells, shRNA-mediated HBO1 silencing and CRISPR/Cas9-induced HBO1 knockout were able to potently inhibit cell viability, growth, proliferation, as well as cell migration and invasion. Significant increase of apoptosis was detected in HBO1-silenced/knockout OS cells. Conversely, ectopic HBO1 overexpression promoted OS cell proliferation and migration. We identified ZNF384 (zinc finger protein 384) as a potential transcription factor of HBO1. Increased binding between ZNF384 and HBO1 promoter was detected in OS cell and tissues, whereas ZNF384 silencing via shRNA downregulated HBO1 and produced significant anti-OS cell activity. In vivo, intratumoral injection of HBO1 shRNA lentivirus silenced HBO1 and inhibited OS xenograft growth in mice. Furthermore, growth of HBO1-knockout OS xenografts was significantly slower than the control xenografts. WM-3835, a novel and high-specific small molecule HBO1 inhibitor, was able to potently suppressed OS cell proliferation and migration, and led to apoptosis activation. Furthermore, intraperitoneal injection of a single dose of WM-3835 potently inhibited OS xenograft growth in SCID mice.

Conclusion: HBO1 overexpression promotes OS cell growth in vitro and in vivo.

Insights in Osteosarcoma by Proton Nuclear Magnetic Resonance Serum Metabonomics

Pediatric osteosarcoma outcomes have improved over the last decades; however, patients who do not achieve a full resection of the tumor, even after aggressive chemotherapy, have the worst prognosis. At a genetic level, osteosarcoma presents many alterations, but there is scarce information on alterations at metabolomic levels. Therefore, an untargeted nuclear magnetic resonance metabonomic approach was used to reveal blood serum alterations, when samples were taken from 21 patients with osteosarcoma aged from 12–20 (18, 86%) to 43 (3, 14%) years before any anticancer therapy were collected. The results showed that metabolites differed greatly between osteosarcoma and healthy control serum samples, especially in lipids, aromatic amino acids (phenylalanine and tyrosine), and histidine concentrations. Besides, most of the loading plots point to protons of the fatty acyls (-CH₃ and -CH₂-) from very-low- and low-density lipoproteins and cholesterol, as crucial metabolites for discrimination of the patients with osteosarcoma from the healthy samples. The relevance of blood lipids in osteosarcoma was highlighted when analyzed together with the somatic mutations disclosed in tumor samples from the same cohort of patients, where six genes linked to the cholesterol metabolism were found being altered too. The high consistency of the discrimination between osteosarcoma and healthy control blood serum suggests that nuclear magnetic resonance could be successfully applied for osteosarcoma diagnostic and prognostic purposes, which could ameliorate the clinical efficacy of therapy.

Multifaceted Functions and Novel Insight Into the Regulatory Role of RNA N6-Methyladenosine Modification in Musculoskeletal Disorders

RNA modifications have emerged as key regulators of transcript expression in diverse physiological and pathological processes. As one of the most prevalent types of RNA modifications, N⁶-methyladenosine (m⁶A) has become the highlight in modulation of various diseases through interfering RNA splicing, translation, nuclear export, and decay. In many cases, the detailed functions of m⁶A in cellular processes and diseases remain unclear. Notably, recent studies have determined the relationship between m⁶A modification and musculoskeletal disorders containing osteosarcoma, osteoarthritis, rheumatoid arthritis, osteoporosis, etc. Herein, this review comprehensively summarizes the recent advances of m⁶A modification in pathogenesis and progression of musculoskeletal diseases. Specifically, the underlying molecular mechanisms, detection technologies, regulatory functions, clinical implications, and future perspectives of m⁶A in musculoskeletal disorders are discussed, with the aim to provide a novel insight into their association.

Integrated Genome-Wide Methylation and Expression Analyses Reveal Key Regulators in Osteosarcoma

Osteosarcoma (OS) is one of the most common types of primary bone tumors in early adolescence with unsatisfied prognosis. Aberrant DNA methylation had been demonstrated to be related to tumorigenesis and progression of multiple cancers and could serve as the potential biomarkers for the prognosis of human cancers. In conclusion, this study identified 18 downregulated hypomethylation genes and 52 upregulated hypomethylation genes in OS by integrating the analysis the GSE97529 and GSE42572 datasets. Bioinformatics analysis revealed that OS-specific methylated genes were involved in regulating multiple biological processes, including chemical synaptic transmission, transcription, response to drug, and regulating immune response. KEGG pathway analysis showed that OS-specific methylated genes were associated with the regulation of Hippo, cAMP calcium, MAPK, and Wnt signaling pathways. By analyzing R2 datasets, this study showed that the dysregulation of these OS-specific methylated genes was associated with the metastasis-free survival time in patients with OS, including CBLN4, ANKMY1, BZW1, KRTCAP3, GZMB, KRTDAP, LY9, PFKFB2, PTPN22, and CLDN7. This study provided a better understanding of the molecular mechanisms underlying the progression and OS and novel biomarkers for the prognosis of OS.

Identification of aberrantly methylated differentially expressed genes targeted by differentially expressed miRNA in osteosarcoma

Background

Osteosarcoma (OS) is the most common primary bone tumors diagnosed in children and adolescents. Recent studies have shown a prognostic role of DNA methylation in various cancers, including OS. The aim of this study was to identify the aberrantly methylated genes that are prognostically relevant in OS.

Methods

The differentially expressed mRNAs, miRNAs and methylated genes (DEGs, DEMs and DMGs respectively) were screened from various GEO databases, and the potential target genes of the DEMs were predicted by the RNA22 program. The protein-protein interaction (PPI) networks were constructed using the STRING database and visualized by Cytoscape software. The functional enrichment and survival analyses of the screened genes was performed using the R software.

Results

Forty-seven downregulated hypermethylated genes and three upregulated hypomethylated genes were identified that were enriched in cell activation, migration and proliferation functions, and were involved in cancer-related pathways like JAK-STAT and PI3K-AKT. Eight downregulated hypermethylated tumor suppressor genes (TSGs) were identified among the screened genes based on the TSGene database. These hub genes are likely involved in OS genesis, progression and metastasis, and are potential prognostic biomarkers and therapeutic targets.

Conclusions

TSGs including PYCARD, STAT5A, CXCL12 and CXCL14 were aberrantly methylated in OS, and are potential prognostic biomarkers and therapeutic targets. Our findings provide new insights into the role of methylation in OS progression.

K-RasG12V/Y40C-PI3K/AKT pathway regulates H1.4S35ph through PKA to promote the occurrence and development of osteosarcoma cancer

Background: Osteosarcoma is prevalent in children and adolescents. H1.4 modification is involved in various types of cancers. Ras pathway is often activated in human cancers. Herein, we explored the effects of Ras pathway through H1.4^S35ph. Methods: Osteosarcoma cancer cell line MG-63 was transfected with Ras gene with G12V and Y40C site mutation. The phosphorylation of H1.4^S35 and AKT was detected by Western blot. Cell viability, cell colonies and migration were analyzed by MTT assay, soft-agar colony formation assay and Transwell assay, respectively. The expression of Ras pathway downstream factors and PKA was detected by qRT-PCR. The relationship between Ras and downstream factors was detected by ChIP. The cell cycle progression was measured by flow cytometry. Results: Transfection with Ras^G12V/Y40C decreased H1.4^S35ph expression while switched on p-AKT^Ser473. Ras^G12V/Y40C increased cell viability, colony numbers and migration while H1.4^S35E (H1.4^S35ph overexpression) led to the opposite results. The regulation of Ras^G12V/Y40C and H1.4^S35E on Ras downstream factors was contrary to each other. Results demonstrated a positive relationship between PKA with H1.4^S35ph with Ras^G12V/Y40C down-regulated both. However, PKA and MDM2 revealed negative regulation with Ras^G12V/Y40C transfection up-regulated MDM2. Conclusion: Ras^G12V/Y40C-PI3K/AKT signal pathway decreased H1.4^S35ph through down-regulation of PKA while up-regulation of MDM2 in MG-63 cells. Highlights H1.4^S35ph is regulated by K-Ras^G12V/Y40-PI3K/AKT in MG-63 cells; Overexpression of H1.4^S35ph regulates MG-63 cell growth; H1.4^S35ph regulates Ras downstream factors; K-Ras^G12V/Y40C-PI3K/AKT activity induces PKA degradation to down-regulate H1.4^S35ph; K-Ras^G12V/Y40C-PI3K/AKT activity involves in PKA degradation via MDM2.

CBX2 is a functional target of miRNA let-7a and acts as a tumor promoter in osteosarcoma

Osteosarcoma is the most common type of primary malignant tumor of skeletal with poor prognosis in children and adolescents. Accumulating evidence indicates that CBX2 is overexpressed in multiple human neoplasm and play a critical role in tumorigenesis and progression. However, its functional role and upstream regulation mechanism in osteosarcoma remain unknown. In the present study, tissue microarray (TMA) analysis was performed to determine the association between CBX2 expression and clinical prognosis of osteosarcoma patients by immunohistochemistry. We also investigated the functional role of CBX2 using small interfering RNA (siRNA) in vitro and in vivo. Additionally, we confirmed the direct binding between CBX2 and let‐7a via qPCR, western blot and luciferase reporter assay. We found that CBX2 is dramatically upregulated in osteosarcoma tissues and high CBX2 expression was correlated with metastasis, recurrence, and chemotherapy response, as well as unfavorable prognosis in patients with osteosarcoma. Similar results were observed in a sarcoma cohort from The Cancer Genome Atlas (TCGA) dataset. Further experiments revealed that CBX2 knockdown significantly impeded osteosarcoma cell proliferation and invasion ability in vitro, and suppressed the tumor growth in tumor xenografts model. Mechanistically, we confirmed that CBX2 is a functional target of miRNA let‐7a. Overexpression of let‐7a inhibits osteosarcoma cell proliferation, which was reversed by CBX2 overexpression. Taken together, our study demonstrates that let‐7a/CBX2 plays a crucial role in osteosarcoma progression. CBX2 could serve as a promising prognostic biomarker and potential therapeutic target for osteosarcoma patients.

Insulin/IGF-1R, SIRT1, and FOXOs Pathways-An Intriguing Interaction Platform for Bone and Osteosarcoma

Aging is a substantial risk factor for the development of osteoarthritis (OA) and, probably, an essential substrate for the development of neoplastic disease of the bone, such as osteosarcoma, which is the most common malignant mesenchymal primary bone tumor. Genetic studies have established that the insulin/insulin-like growth factor 1 (IGF-1)/phosphatidylinositol-3 kinase (PI3K)/AKT (Protein Kinase B) signal transduction pathway is involved across species, including nematodes, fruit flies, and mammals. SIRT1, a phylogenetically-conserved family of deacetylases, seems to play pleiotropic effects in epithelial malignancies of the liver and interact with the IGF-1/PI3K/AKT signal transduction pathway. Some of the most critical processes in degenerative conditions may indeed include the insulin/IGF1R and SIRT1 signaling pathways as well as some specific transcription factors. The Forkhead box O (FOXO) transcription factors (FOXOs) control diverse cellular functions, such as metabolism, longevity, and cell death. FOXOs play a critical role in the IGF-1/PI3K/AKT signal transduction pathway. FOXOs can indeed be modulated to reduce age-related diseases. FOXOs have advantageous inhibitory effects on fibroblast and myofibroblast activation, which are accompanied by a subsequent excessive production of extracellular matrix. FOXOs can block or decrease the fibrosis levels in numerous organs. Previously, we observed a correlation between nuclear FOXO3 and high caspase-8 expression, which induces cellular apoptosis in response to harmful external stimuli. In this perspective, we emphasize the current advances and interactions involving the insulin/IGF1R, SIRT1, and FOXOs pathways in the bone and osteosarcoma for a better understanding of the mechanisms potentially underpinning tissue degeneration and tumorigenesis.

Screening of candidate key genes associated with human osteosarcoma using bioinformatics analysis

The aim of the present study was to identify the key genes associated with osteosarcoma (OS) using a bioinformatics approach. Microarray data (GSE36004) was downloaded from the Gene Expression Omnibus database, including 19 OS cell lines and 6 normal controls. Differentially expressed genes (DEGs) in the OS cell lines were identified using the Limma package, and differentially methylated regions were screened with methyAnalysis in R. Copy number analysis was performed and genes with copy number gains/losses were further screened using DNAcopy and cghMCR packages. Functional enrichment analyses were performed using the Database for Annotation, Visualization and Integrated Discovery online tool, and protein-protein interactions were identified based on information obtained from the Search Tool for the Retrieval of Interacting Genes database. A total of 47 downregulated genes were screened in hyper-methylated regions, including the fragment crystallizable (Fc) region of immunoglobulin E, high affinity I, receptor for; γ polypeptide (FCER1G), leptin (LEP) and feline Gardner-Rasheed sarcoma viral oncogene homolog (FGR). In addition, a total of 17 upregulated genes, including the TPase family, AAA domain containing 2 (ATAD2) and cyclin-dependent kinase 4 (CDK4), exhibited copy number gains, while 5 downregulated genes, including Rho GTPase activating protein 9 (ARHGAP9) and major histocompatibility complex, class II, DO α (HLA-DOA), exhibited copy number losses. These results indicate that hyper-methylation of FCER1G, LEP, and FGR may serve a crucial function in the development of OS. In addition, copy number alterations of these DEGs, including ATAD2, CDK4, ARHGAP9 and HLA-DOA, may also contribute to OS progression. These DEGs may be candidate targets for the diagnosis and treatment of this disease.

Histone Deacetylase Inhibitor Trichostatin a Promotes the Apoptosis of Osteosarcoma Cells through p53 Signaling Pathway Activation

Purpose: The purpose of this study was to investigate the profile of histone deacetylase (HDAC) activity and expression in osteosarcoma cells and tissues from osteosarcoma patients and to examine the mechanism by which a histone deacetylase (HDAC) inhibitor, Trichostatin A (TSA), promotes the apoptosis of osteosarcoma cells. Methods: HDAC activity and histone acetyltransferase (HAT) activity were determined in nuclear extracts of MG63 cells, hFOB 1.19 cells and tissues from 6 patients with primary osteosarcoma. The protein expression of Class I HDACs (1, 2, 3 and 8) and the activation of the p53 signaling pathway were examined by Western blot. Cell growth and apoptosis were determined by 3-(4, 5-dimethyl-2-thiazolyl)-2H-tetrazolium bromide (MTT) assay and flow cytometry, respectively. Results: Nuclear HDAC activity and class I HDAC expression were significantly higher in MG63 cells than in hFOB 1.19 cells, and a similar trend was observed in the human osteosarcoma tissues compared with the paired adjacent non-cancerous tissues. TSA significantly inhibited the growth of MG63 cells and promoted apoptosis in a dose-dependent manner through p53 signaling pathway activation. Conclusion: Class I HDACs play a central role in the pathogenesis of osteosarcoma, and HDAC inhibitors may thus have promise as new therapeutic agents against osteosarcoma.

Activin A regulates proliferation, invasion and migration in osteosarcoma cells

Activin A is a member of the TGF‑β superfamily. Previous studies have demonstrated that activin A exhibited pluripotent effects in several tumours. However, the roles of activin A signaling in osteosarcoma pathogenesis have not been previously investigated. Therefore, the present study aimed to investigate the effects of activin A on osteosarcoma cell proliferation, invasion and migration. Firstly, the expression of activin A in osteosarcoma cell lines (MG63, SaOS‑2 and U2OS) and a human osteoblastic cell line (hFOB1.19) was detected using reverse transcription quantitative polymerase chain reaction and western blotting. Activin A was upregulated in osteosarcoma cell lines compared with hFOB1.19 cells. To investigate the effects of activin A on osteosarcoma cell proliferation, invasion and migration, MG63 cells were generated in which activin A was either overexpressed or depleted. MTT staining, propidium iodide staining and a Transwell assay were used to analyze the cell cycle, proliferation, invasion and migration of MG63 cells, respectively. The results of the present study revealed that the abilities of proliferation, invasion and migration were suppressed in MG63 cells in which activin A was depleted, while they were enhanced in activin A-overexpressing cells. In conclusion, the results of the present study suggested that activin A may facilitate proliferation, invasion and migration of osteosarcoma cells, and it may therefore be a potential target for the treatment of osteosarcoma.

Separate and combined effects of DNMT and HDAC inhibitors in treating human multi-drug resistant osteosarcoma HosDXR150 cell line

Understanding the molecular mechanisms underlying multi-drug resistance (MDR) is one of the major challenges in current cancer research. A phenomenon which is common to both intrinsic and acquired resistance, is the aberrant alteration of gene expression in drug-resistant cancers. Although such dysregulation depends on many possible causes, an epigenetic characterization is considered a main driver. Recent studies have suggested a direct role for epigenetic inactivation of genes in determining tumor chemo-sensitivity. We investigated the effects of the inhibition of DNA methyltransferase (DNMT) and hystone deacethylase (HDAC), considered to reverse the epigenetic aberrations and lead to the re-expression of de novo methylated genes in MDR osteosarcoma (OS) cells. Based on our analysis of the HosDXR150 cell line, we found that in order to reduce cell proliferation, co-treatment of MDR OS cells with DNMT (5-Aza-dC, DAC) and HDAC (Trichostatin A, TSA) inhibitors is more effective than relying on each treatment alone. In re-expressing epigenetically silenced genes induced by treatments, a very specific regulation takes place which suggests that methylation and de-acetylation have occurred either separately or simultaneously to determine MDR OS phenotype. In particular, functional relationships have been reported after measuring differential gene expression, indicating that MDR OS cells acquired growth and survival advantage by simultaneous epigenetic inactivation of both multiple p53-independent apoptotic signals and osteoblast differentiation pathways. Furthermore, co-treatment results more efficient in inducing the re-expression of some main pathways according to the computed enrichment, thus emphasizing its potential towards representing an effective therapeutic option for MDR OS.

Epigenetic changes in carcinogenesis of gallbladder

Gallbladder cancer (GBC) is a lethal and a common malignancy affecting mostly females. There are restricted high incidence pockets across the world and in northern India highest incidence of GBC is reported from the Gangetic belt. The etiology of this disease remains largely unknown though several risk factors have been stated. The genetic aberrations in GBC involving mutations in tumor suppressor genes and oncogenes have been reported in literature. However, there is scarcity of data regarding epigenetic changes that may also be involved in gallbladder carcinogenesis. This review attempts to summarize our current understanding of the epigenetic changes in GBC.

MicroRNA 181a improves proliferation and invasion, suppresses apoptosis of osteosarcoma cell

MicroRNA 181a (miR-181a) was found dysregulated in a variety of human cancers and significantly associated with clinical outcome of cancer patients. However, the direct role of miR-181a has not yet been characterized in osteosarcoma progression. This study was aimed at investigating the effects of miR-181a on osteosarcoma cell biological behavior. First, the expression of miR-181a in osteosarcoma cell lines (MG63, HOS, SaOS-2, and U2OS) and a human osteoblastic cell line (hFOB1.19) was detected by qRT-PCR. Results showed that miR-181a was overexpressed in osteosarcoma cell lines compared to human osteoblastic cell line (hFOB1.19). To investigate the effects of miR-181a on proliferation, apoptosis, and invasion of osteosarcoma cells, we generated human osteosarcoma MG63 cells in which miR-181a was either overexpressed or depleted. The MG63 cell viability, cycle, apoptosis, and invasive ability were analyzed by 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide staining, propidium iodide (PI) staining, Annexin V-FITC/PI double staining, and Transwell invasion experiment, respectively. The results showed that MG63 cell viability, proliferation, and invasive abilities were suppressed, and the apoptosis was enhanced in the group with underexpression of miR-181a. The viability, proliferation, and invasive abilities were improved, and the apoptosis was inhibited in the group with overexpression of miR-181a. The results from Western blotting indicated that miR-181a might be associated with the up-regulation of bcl-2 and matrix metalloproteinase 9 and the down-regulation of tissue inhibitor of metalloproteinases-3 and p21 in MG63 cells. Taken together, our results suggested that miR-181a might facilitate proliferation and invasion and suppress apoptosis of osteosarcoma cells, which might be a potential target for the treatment of osteosarcoma.

Reduced mitochondrial DNA copy number in Chinese patients with osteosarcoma

A plethora of somatic mutations and germline variations in mitochondrial DNA (mtDNA) have been increasingly reported in numerous cancer entities including osteosarcoma. However, it remains largely unclear whether mtDNA copy number changes occur during the multistep process of osteosarcoma carcinogenesis. For this purpose, we determined quantitative mtDNA levels in 31 primary osteosarcoma specimens and 5 normal bone tissue samples using a real-time polymerase chain reaction assay. Our data showed that the average mtDNA amount was significantly reduced in osteosarcoma tissues compared with normal bone controls. The copy number of mtDNA was statistically associated with tumor metastasis. There was an approximately 2-fold decrease of mtDNA quantity in tumors with metastasis than that in low-grade tumors without metastasis. Furthermore, change in mtDNA content was linked with somatic mutations in the D-loop regulatory region. Tumors carrying somatic D-loop mutations, at the polycytidine stretch between nucleotide positions 303 and 309 or close to the replication origin sites of the heavy strand, had significantly lowered mtDNA levels in comparison with those without mutations. Taken together, these results provide evidence for the first time that reduced mtDNA content may be critically implicated in the development and/or progression of osteosarcoma. Somatic D-loop mutation is likely one key factor among others leading to altered mtDNA amount in osteosarcoma.

Somatic mitochondrial DNA mutations in Chinese patients with osteosarcoma

Somatic mutations in mitochondrial DNA (mtDNA) have been long proposed to drive the pathogenesis and progression of human malignancies. Previous investigations have revealed a high frequency of somatic mutations in the D-loop control region of mtDNA in osteosarcoma. However, little is known with regard to whether or not somatic mutations also occur in the coding regions of mtDNA in osteosarcoma. To test this possibility, in the present study we screened somatic mutations over the full-length mitochondrial genome of 31 osteosarcoma tumour tissue samples, and corresponding peripheral blood samples from the same cohort of patients. We detected a sum of 11 somatic mutations in the mtDNA coding regions in our series. Nine of them were missense or frameshift mutations that have the potential to hamper mitochondrial respiratory function. In combination with our earlier observations on the D-loop fragment, 71.0% (22/31) of patients with osteosarcoma carried at least one somatic mtDNA mutation, and a total of 40 somatic mutations were identified. Amongst them, 29 (72.5%) were located in the D-loop region, two (5%) were in the sequences of the tRNA genes, two (5%) were in the mitochondrial ATP synthase subunit 6 gene and seven (17.5%) occurred in genes encoding components of the mitochondrial respiratory complexes. In addition, somatic mtDNA mutation was not closely associated with the clinicopathological characteristics of osteosarcoma. Together, these findings suggest that somatic mutations are highly prevalent events in both coding and non-coding regions of mtDNA in osteosarcoma. Some missense and frameshift mutations are putatively harmful to proper mitochondrial activity and might play vital roles in osteosarcoma carcinogenesis.

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.