- Open Access
Identification of N6-methylandenosine related LncRNAs biomarkers associated with the overall survival of osteosarcoma
BMC Cancer volume 21, Article number: 1285 (2021)
Osteosarcoma (OS) is a differentiation disease caused by the genetic and epigenetic differentiation of mesenchymal stem cells into osteoblasts. OS is a common, highly malignant tumor in children and adolescents. Fifteen to 20 % of the patients find distant metastases at their first visit. The purpose of our study was to identify biomarkers for tracking the prognosis and treatment of OS to improve the survival rate of patients.
Materials and methods
In this study, which was based on Therapeutically Applicable Research to Generate Effective Treatments (TARGET), we searched for m6A related lncRNAs in OS. We constructed a network between lncRNA and m6A, and built an OS prognostic risk model.
We identified 14,581 lncRNAs by using the dataset from TARGET. We obtained 111 m6A-related lncRNAs through a Pearson correlation analysis. A network was built between lncRNA and m6A genes. Eight m6A-related lncRNAs associated with survival were identified through a univariate Cox analysis. A selection operator (LASSO) Cox regression was used to construct a prognostic risk model with six genes (RP11-286E11.1, LINC01426, AC010127.3, DLGAP1-AS2, RP4-657D16.3, AC002398.11) obtained through least absolute shrinkage.
We also discovered upregulated levels of DLGAP1-AS2 and m6A methylation in osteosarcoma tissues/cells compared with normal tissues/osteoblasts cells.
We constructed a risk score prognosis model of m6A-related lncRNAs (RP11-286E11.1, LINC01426, AC010127.3, DLGAP1-AS2, RP4-657D16.3, AC002398.11) using the dataset downloaded from TRAGET. We verified the value of the model by dividing all samples into test groups and training groups. However, the role of m6A-related lncRNAs in osteosarcoma needs to be further researched by cell and in vivo studies.
Osteosarcoma (OS), a disease that mainly affects children and adolescents, is a common malignant tumor that occurs in the mesenchymal tissue of the metaphysis of long bones [1,2,3]. Although extensive surgical resection combined with chemotherapy and radiotherapy has achieved certain good results, about 40–50% of patients experience lung metastasis [4, 5]. The five-year survival rate of the patients with lung metastases is only 28%. Thus, it is important to develop a novel therapeutic approach to effectively treatOS.
More than 60% of all RNA modifications are RNA methylation. N6-methyladenosine (m6A) is the most common form of modification of human messenger RNA (mRNA) and long non-coding RNA (lncRNA) [6, 7]. M6A modification mainly occurs on the adenine of the RRACH sequence, and its function is mainly related to writers, readers, and erasers. Writers (RBM15B, METTL3, METTL14, METTL16, RBM15, WTAP, VIRMA [KIA1499], and ZC3H13) are methyltransferases [8, 9]. The eraser is a demethylase that can reverse methylation and includes FTO and ALKBH5. The reader is signal transducers, with m6A binding protein that can bind to M6A, including FMR1, HNRNPC, HNRNPA2B1, IGFBP1, IGFBP2, IGFBP3, LRPPRC, RBMX, YTHDC1, YTHDC2, YTHDF1, YTHDF2, and YTHDF3. Accumulating evidence has shown that m6A modification plays an important role in tumorigenesis and tumor progression, such as hepatocellular carcinoma, lung cancer, and leukemia [10,11,12]. Wang et al. revealed that METTL3 promotes tumorigenesis by up-regulating m6A modification of APC . Qu et al. found that HBx and ALKBH5 promote hepatocellular carcinogenesis through a positive feedback loop . Cao et al. reported that IGF2BP2 can regulate insulin sensitivity and glucose metabolism to mediate cancer pathogenesis . Additionally, ALKBH5 has been shown to negatively regulate the osteogenic differentiation of mesenchymal stem cells through PRMT6 .
LncRNAs, a class of RNA transcripts, are greater than 200 nucleotides in length . Although they generally cannot encode proteins or peptides, lncRNAs play an important role in cell activation, proliferation, differentiation, apoptosis and metabolism by regulating gene expression and function . Recent research had shown that dysregulated lncRNA plays an important role in many diseases, such as tumor, cardiovascular diseases, and metabolic disorders [17, 18]. Dysregulated lncRNA is closely associated with tumor growth and metastasis observed in cancers such as breast, lung, liver, and colorectal cancers [19,20,21]. Some studies have suggested that lncRNAs may potentially serve as biomarkers and targets in diagnosing and treating tumors.
Bioinformatics is a new field of biological research, which involves the processing and analysis of biological data using mathematical, statistical and computational methods. Due to the large amount of data generated by new technologies such as genome sequencing and microarray chip technology, the traditional gene-by-gene method is not enough to meet the growth and demands of biological research. Therefore, bioinformatics is a valuable way to expand biological insight and promote the development of new therapeutic approaches. We searched for m6A -related lncRNAs, based on Therapeutically Applicable Research to Generate Effective Treatments (TARGET, https://ocg.cancer.gov/programs/target). Using bioinformatic and statistical analysis methods, we constructed a prognostic risk model to identify biomarkers related to OS prognosis and treatment to improve the survival rate of patients.
Data collection and data processing
We downloaded the OS gene expression profiles, and the corresponding clinical data from the TARGET database, which comprised a total of 84 OS samples. LncRNAs were identified using the Ensemble IDs of the genes. All lncRNAs were extracted base on the perl. Expression matrixes of 23 m6A genes were extracted from previous publications, including the expression data on writers (RBM15B, METTL3, METTL14, METTL16, RBM15, WTAP, VIRMA [KIA1499], and ZC3H13), readers (FMR1, HNRNPC, HNRNPA2B1, IGFBP1, IGFBP2, IGFBP3, LRPPRC, RBMX, YTHDC1, YTHDC2, YTHDF1, YTHDF2, and YTHDF3) and erasers (FTO and ALKBH5).
Identification of m6A-related lncRNAs
Using a Pearson correlation analysis, we searched for m6A-related lncRNAs in each gene (with the | Pearson R| > 0.3 and P < 0.01). The interaction betweenlncRNA and m6A was visualized through the igraph R package. We then used a univariate Cox regression analysis to mine the prognostic m6A-related lncRNAs (pFilter = 0.05).
Identification of cluster based on m6A-related lncRNAs and gene set enrichment analysis
The data matrix of m6A-related lncRNAs was inputted to ConsensusClusterPlus. Stability evidence for a given number of groups (k) and cluster assignments was outputted. We set the number of clusters, k, from 2 to 10, and found the ideal k from it. A heatmap was generated using the ‘pheatmap’ R packages. A Kaplan–Meier (KM) survival curve was drawn to assess prognoses between the clusters. Gene Set Enrichment Analysis has been widely applied to perform KEGG (Kyoto Encyclopedia of Genes and Genomes) and GO (Gene Ontology) analysis. The top ten pathways were shown.
Construction of prognosis risk score model
Utilizing a least absolute shrinkage and selection operator (LASSO) regression analysis, all m6A-related lncRNAs were analyzed to select the most optimal prognostic biomarkers and to construct the risk score model. A LASSO regression analysis is a penalty regression method that reduces overfitting by simultaneously performing shrinkage and model selection. Some coefficients can be compressed to 0 based by constructing a penalty function, resulting in a more refined model. Therefore, the superiority of subset shrinkage is retained. In addition, the LASSO model can perform biased estimation of multicollinearity data processing and realize the selection of variables in the estimation. The ideal method to solve the multicollinearity problem in a regression analysis is a LASSO regression analysis.
Validation of prognosis risk score model
In this study, we calculated the prognosis risk scores based on a formula as follows:
We divided all samples into low-risk groups and high-risk groups, using the median risk score as the cutoff value. To avoid random allocation bias, all OS patients were randomly grouped at a ratio of 0.5:0.5 into a training cohort and a test cohort. Using a KM survival analysis, and calculating the area under the ROC curve (AUC), and etc., we assessed the accuracy and efficiency of the risk score model. We performed the univariate and multivariate Cox survival analyses according to clinical factors such as sex, age, primary site of OS and metastasis to assess the predictive independence of the risk score model. In addition, KM survival analysis was applied to confirm the effectiveness of the risk score model for different clinical features.
RNA m6A quantification
The m6A RNA methylation assay kit (ab185912; Abcam, UK) was utilized to measure the m6A content. Briefly, extracted RNAs were coated on each well, followed by adding capture antibody and detection antibody. Finally, the m6A levels were detected by measuring the absorbance at 450 nm using a microplate reader.
RNA isolation and qRT-PCR
Total RNA was isolated using RNeasy Mini Kit (Qiagen, Germany), and reversely transcribed according to the standard of cDNA Synthesis Kit (Vazyme, China). Then, qRT-PCR was performed using SYBR Kit (Vazyme, China) to evaluate the mRNA levels of DLGAP1-AS2 and GAPDH according to the supplier’s instructions.
Flow chart of analysis
We designed a protocol to analyze m6A-related lncRNAs and constructed the prognosis risk score model (Fig. 1).
Screening m6A-related lncRNAs in OS patients
First, we identified 14,581 lncRNAs by using the OS dataset from TARGET. We then extracted the expression matrixes of 23 m6A-related lncRNA. A lncRNA was defined as an m6A-related lncRNA with an expression value was correlated with one or more of the 23 m6A-related genes (| Pearson R| > 0.3 and p < 0.01). We obtained 111 m6A-related lncRNAs through Pearson correlation analysis. The network was built between lncRNA and m6A (Fig. 2).
Identification of clusters according to m6A-related lncRNAs
Eight m6A-related lncRNAs were identified using univariate Cox analysis (RP11-799D4.4, RP11-286E11.1, LINC01426, AC010127.3, DLGAP1-AS2, RP4-657D16.3, RP1-178F15.5, and AC002398.11). The results of the 8 m6A-related lncRNAs are shown in Table 1. According to the consensus clustering analysis, the optimal number of clusters was k = 2 (Fig. 3A, B, C, D). Compared with the Cluster 2 (C2) subtype, most genes were highly expressed in the Cluster 1 (C1) subtype.
Prognostic analysis and biological functional analysis of the m6A-related lncRNAs
The KM survival curve was drawn to compare the survival of the two molecular subtypes. In the results, the survival rate of C2 was significantly better than that of C1 (P < 0.05; Fig. 4A). A heatmap of the relationships between the expression levels of these eight m6A-related lncRNAs and clinicopathological characteristics were shown in Fig. 4B. GSEA was performed to analyze the KEGG and GO analysis in C1 and C2(Fig. 4C,D).
Construction of the prognostic risk model
Using a LASSO-Cox regression analysis, we further selected appropriate m6A-related lncRNAs to maintain a high accuracy rate. We show the changing trajectory of each independent variable in Fig. 5A. The confidence interval of each lambda was shown in Fig. 5B. These figures showed that a total of six genes were selected as the target genes for subsequent analysis when the model reached the optimal value. The six m6A-related lncRNAs were obtained as the most suitable predictors for prognosis (RP11-286E11.1, LINC01426, AC010127.3, DLGAP1-AS2, RP4-657D16.3, AC002398.11). A risk scoring model to predict the prognosis of OS patients was constructed as follows: RP11-286E11.1 * 0.00097667835205502 + LINC01426 * (− 0.00245896608774788) + AC010127.3 * (0.00475697115852846) + DLGAP1-AS2 * (0.00138111340031212) + RP4-657D16.3 * (− 0.00419685380105118) + AC002398.11 * (− 0.000637455261090537).
Evaluation of the risk model
We randomly divided the 84 samples into two groups, 42 samples were used as the testing cohort samples, and the other 42 samples were used as the training cohort samples. The risk score distribution of the training cohort is shown in Fig. 6A, and the survival time of the training cohort is shown in Fig. 6B. The average area under the curve (AUC) of the training cohort reached 0.787 (Fig. 6C). The training cohort have 21 samples as high risk and 21 as low risk. There was a significant difference between high-risk samples and low-risk samples in the KM survival curves (P = 0.002, Fig. 6D). The heatmap showed the high-risk samples and low-risk samples had different m6A-related lncRNA expression levels (Fig. 6E). The risk score distribution of the test cohort is shown in Fig. 7A, and Fig. 7B shows the survival time of the test cohort. The AUC of test cohort was 0.816 (Fig. 7C). The test cohort had 17 high-risk samples and 25 low-risk samples.. There was a significant difference between high-risk samples and low-risk samples in the KM survival curves (P = 0.003, Fig. 7D). The heatmap shows that the high-risk samples and low-risk samples had different expression levels of m6A-related lncRNA (Fig. 7E).
Prognostic analysis of the risk model and clinical features
Using a univariate and multivariate Cox regression analysis, we analyzed clinical variables, such as age, gender, site, and metastasis (Fig. 8A, B, C, D). The P-value of the univariate and multivariate Cox survival analysis in the training group and the test groupwas < 0.05. The results suggested that the model was an independent risk factor of survival for OS patients. We have drawn the KM survival curves according to the clinical characteristics to compare the survival of the high-risk group and low-risk group. In the dataset, 38 patients were 14 years of age or younger, and 46 patients older than 14 years. In the group of age over 14, there are 18 patients as high risk and 28 patients as low risk (P < 0.001, Fig. 9A). In the group 14 years of age or younger, 20 patients were high risk and 18 patients as low risk (P = 0.010, Fig. 9B). Of the 84 OS patients, there were 37 women and 47 men. Seventeen women were classified into the high-risk group and 20 women were classified into the low-risk group (P < 0.001, Fig. 9C). Twenty-one men were classified into the high-risk group and 26 men were classified into the low-risk group (P = 0.028, Fig. 9D). Twenty patients had metastases, while the remaining 64 patients had no metastases. Among the patients with OS metastasis, 14 patients were in the high-risk group and 6 patients were in the low-risk group (P = 0.015, Fig. 9E). Among patients without metastasis, 24 patients were in the high-risk group and 40 patients s were in the low-risk group (P = 0.008, Fig. 9F). The primary OS site was the leg or foot in 77 patients, the arm or hand in 5 patients, and the pelvis in 2 patients., s. In the group of patients with the leg as the primary OS site, 35 patients s were high risk and 42 patientswere low risk (P < 0.001, Fig. 9G). Because the number of patients with the arm (Fig. 9H) or the pelvis (Fig. 9I) as the primary OS site was too small, the results were meaningless.
The m6A quantitative experiment showed that the m6A level of RNA in human osteosarcoma cell lines (MG63,143B) was higher than that in human osteoblast cell line(hFOB1.19) (P < 0.05) (Fig. 10). The qRT-PCR analysis results showed that DLGAP1-AS2 expression increased in tumor tissues (P = 0.07) (Fig. 11). The p value > 0.05 may be caused by insufficient tissue samples. Although the P value> 0.05, the expression level of DLGAP1-AS2 in each pair of tissues was higher in the tumor group than in the normal group.
Our study investigated the prognostic significance of m6A-related lncRNAs in 84 OS patients from the TARGET dataset. We identified a total of 111 m6A-related lncRNA in our study, and eight m6A-related lncRNAs were selected by conducting a univariate Cox analysis. We then performed a KEGG analysis and a GO analysis according to the cluster assignments.. Six m6A-related lncRNAs (RP11-286E11.1, LINC01426, AC010127.3, DLGAP1-AS2, RP4-657D16.3, and AC002398.11) were identified in a LASSO Cox analysis as the most suitable risk score model for prognosis. The application value of the model was tested with ROC curve analysis and KM survival curves. In addition, a univariate and multivariate Cox regression analysis also confirmed the validity of the model.
According to GO analysis, the m6A-related lncRNAs are enriched in histone H3-K27 methylation, pericentric heterochromatin, and phosphatase activator activity. Histone H3 mutations have been found to play a role in a variety of cancers, such as paediatric brain tumors . Penin et al. reported that pericentric heterochromatin of chromosome 9 was the primary target of HSF1 in both normal and tumor heat-shocked cells . Zheng et al. found that LINC00174 promotes the metastasis and growth of OS through upregulating slingshot protein phosphatase 2 expression . A KEGG analysis found that the m6A-related lncRNAs are enriched in the mTOR pathway. Mickymaray et al. found that rhaponticin can inhibit OS growth by inhibiting the PI3K-Akt-mTOR pathway . Liu et al. reported that miR-140 inhibits OS development by influencing ubiquitin-specific protease 22 and promoting p21 expression .
A number of studies have suggested that m6A modification may play a regulatory role in the development of cancer. M6A-related could serve as new therapeutic targets and prognostic biomarkers for OS . Miao et al. found that METTL3 regulates m6A levels of LEF1 and activates the Wnt/β-Catenin signaling pathway to promote OS progression . Zhou et al. demonstrated that METTL3 acts as an oncogene in OS by regulating ATAD2 . Liu et al. reported that METTL14 overexpression promotes OS cell apoptosis and inhibits tumor progression by activating caspase3 . Chen et al. reported that WTAP promotes OS growth and metastasis by inhibiting HMBOX1 expression . Yuan et al. found that ALKBH5 suppresses OS progression by inhibiting the pre-miR-181b-1/YAP signaling axis . However, the role of m6A-related lncRNAs in OS is still unclear. We found that our risk score model built with six m6A-related lncRNAs could accurately predict the prognosis of OS. Liu et al. found that LINC01426 facilitates the progression and stemness in lung adenocarcinoma . In addition, LINC01426 accelerates glioblastoma progression by regulating miR-345-3p/VAMP8 signaling axis . Liu et al. reported DLGAP1-AS2 is associated with poor prognosis in cholangiocarcinoma . In hepatocellular carcinoma cell, DLGAP1-AS2 knockdown inhibits cell metastasis by regulating miR-154-5p methylation . We also discovered upregulated levels of DLGAP1-AS2 and m6A methylation in osteosarcoma tissues/cells compared with normal tissues/osteoblasts cells. Among the six m6A-related lncRNAs used to construct the risk prognosis model, several genes have been reported to be associated with tumor prognosis such as respiratory system tumors and digestive system tumors [35, 36]. However, no research has reported the relationship between OS and these genes or how the lncRNAs interact with m6A-related genes. Therefore, the aim of our study is to identify lncRNAs associated with OS prognosis and provide new targets to prevent poor OS prognosis. A limitation is that we have not further researched the role of m6A-related lncRNAs in osteosarcoma through cell and in vivo studies. Using the dataset downloaded from TRAGET, we constructed a risk score prognosis model of m6A-related lncRNAs. By dividing all samples into test groups and training groups, we verified the value of the model. However, the role of m6A-related lncRNAs in osteosarcoma needs to be further researched by cell and in vivo studies.
Availability of data and materials
All data are fully available without restriction. In this study, publicly available datasets were analyzed. This data can be found here: Therapeutically Applicable Research to Generate Effective Treatments (TARGET, https://ocg.cancer.gov/programs/target).
Wen JF, Jiang YQ, Li C, Dai XK, Wu T, Yin WZ. LncRNA-SARCC sensitizes osteosarcoma to cisplatin through the miR-143-mediated glycolysis inhibition by targeting hexokinase 2. Cancer Biomark. 2020;28(2):231–46.
Sasaki R, Osaki M, Okada F. Microrna-based diagnosis and treatment of metastatic human osteosarcoma. Cancers. 2019;11(4):553.
Shimbo K, Miyaki S, Ishitobi H, Kato Y, Kubo T, Shimose S, et al. Exosome-formed synthetic microRNA-143 is transferred to osteosarcoma cells and inhibits their migration. Biochem Biophys Res Commun. 2014;445(2):381–7.
Hirahata M, Osaki M, Kanda Y, Sugimoto Y, Yoshioka Y, Kosaka N, et al. PAI-1, a target gene of miR-143, regulates invasion and metastasis by upregulating MMP-13 expression of human osteosarcoma. Cancer Med. 2016;5(5):892–902.
Yu X, Hu L, Li SY, Shen J, Wang DL, Xu RJ, et al. Long non-coding RNA taurine upregulated gene 1 promotes osteosarcoma cell metastasis by mediating HIF-1 alpha via miR-143-5p. Cell Death Dis. 2019;10(4):280.
Tu Z, Wu L, Wang P, Hu Q, Tao C, Li K, et al. N6-methylandenosine-related lncRNAs are potential biomarkers for predicting the overall survival of lower-grade glioma patients. Front Cell Dev Biol. 2020;8:642.
Li J, Wang F, Liu Y, Wang H, Ni B. N(6)-methyladenosine (m(6)A) in pancreatic cancer: regulatory mechanisms and future direction. Int J Biol Sci. 2021;17(9):2323–35.
Tang M, Lv Y. The role of N(6) -Methyladenosine modified circular RNA in pathophysiological processes. Int J Biol Sci. 2021;17(9):2262–77.
Lan Q, Liu PY, Bell JL, Wang JY, Hüttelmaier S, Zhang XD, et al. The emerging roles of RNA m(6)A methylation and demethylation as critical regulators of tumorigenesis, drug sensitivity, and resistance. Cancer Res. 2021;81(13):3431–40.
Zhang D, Ning J, Okon I, Zheng X, Satyanarayana G, Song P, et al. Suppression of m6A mRNA modification by DNA hypermethylated ALKBH5 aggravates the oncological behavior of KRAS mutation/LKB1 loss lung cancer. Cell Death Dis. 2021;12(6):518.
Cao J, Yan W, Ma X, Huang H, Yan H. Insulin-like growth factor 2 mRNA-binding protein 2 - a potential link between type 2 diabetes mellitus and cancer. J Clin Endocrinol Metab. 2021;106(10):2807–18.
Gong H, Liu L, Cui L, Ma H, Shen L. ALKBH5-mediated m6A-demethylation of USP1 regulated T-cell acute lymphoblastic leukemia cell glucocorticoid resistance by Aurora B. Mol Carcinog. 2021;60(9):644–57.
Wang W, Shao F, Yang X, Wang J, Zhu R, Yang Y, et al. METTL3 promotes tumour development by decreasing APC expression mediated by APC mRNA N(6)-methyladenosine-dependent YTHDF binding. Nat Commun. 2021;12(1):3803.
Qu S, Jin L, Huang H, Lin J, Gao W, Zeng Z. A positive-feedback loop between HBx and ALKBH5 promotes hepatocellular carcinogenesis. BMC Cancer. 2021;21(1):686.
Li Z, Wang P, Li J, Xie Z, Cen S, Li M, et al. The N(6)-methyladenosine demethylase ALKBH5 negatively regulates the osteogenic differentiation of mesenchymal stem cells through PRMT6. Cell Death Dis. 2021;12(6):578.
Chen W, Liu S, Wang F. Potential impact and mechanism of long non-coding RNAs on cancer and associated T cells. J Cancer. 2021;12(16):4873–82.
Guo R, Zou B, Liang Y, Bian J, Xu J, Zhou Q, et al. LncRNA RCAT1 promotes tumor progression and metastasis via miR-214-5p/E2F2 axis in renal cell carcinoma. Cell Death Dis. 2021;12(7):689.
Che D, Fang Z, Mai H, Xu Y, Fu L, Zhou H, et al. The lncRNA ANRIL gene rs2151280 GG genotype is associated with increased susceptibility to recurrent miscarriage in a southern Chinese population. J Inflamm Res. 2021;14:2865–72.
Liu X, Yin Z, Xu L, Liu H, Jiang L, Liu S, et al. Upregulation of LINC01426 promotes the progression and stemness in lung adenocarcinoma by enhancing the level of SHH protein to activate the hedgehog pathway. Cell Death Dis. 2021;12(2):173.
Hou T, Ye L, Wu S. Knockdown of LINC00504 inhibits the proliferation and invasion of breast cancer via the downregulation of miR-140-5p. OncoTargets Ther. 2021;14:3991–4003.
Chen K, Zhang Z, Yu A, Li J, Liu J, Zhang X. lncRNA DLGAP1-AS2 knockdown inhibits hepatocellular carcinoma cell migration and invasion by regulating miR-154-5p methylation. Biomed Res Int. 2020;2020:6575724.
Sievers P, Sill M, Schrimpf D, Stichel D, Reuss DE, Sturm D, et al. A subset of pediatric-type thalamic gliomas share a distinct DNA methylation profile, H3K27me3 loss and frequent alteration of EGFR. Neuro-oncology. 2021;23(1):34–43.
Penin J, Dufour S, Faure V, Fritah S, Seigneurin-Berny D, Col E, et al. Chromosome Y pericentric heterochromatin is a primary target of HSF1 in male cells. Chromosoma. 2021;130(1):53–60.
Zheng C, Li R, Zheng S, Fang H, Xu M, Zhong L. LINC00174 facilitates cell proliferation, cell migration and tumor growth of osteosarcoma via regulating the TGF-β/SMAD signaling pathway and upregulating SSH2 expression. Front Mol Biosci. 2021;8:697773.
Mickymaray S, Alfaiz FA, Paramasivam A, Veeraraghavan VP, Periadurai ND, Surapaneni KM, et al. Rhaponticin suppresses osteosarcoma through the inhibition of PI3K-Akt-mTOR pathway. Saudi J Biol Sci. 2021;28(7):3641–9.
Liu W, Wang D, Liu L, Wang L, Yan M. miR-140 inhibits osteosarcoma progression by impairing USP22-mediated LSD1 stabilization and promoting p21 expression. Mol Ther Nucleic Acids. 2021;24:436–48.
Li J, Rao B, Yang J, Liu L, Huang M, Liu X, et al. Dysregulated m6A-related regulators are associated with tumor metastasis and poor prognosis in osteosarcoma. Front Oncol. 2020;10:769.
Miao W, Chen J, Jia L, Ma J, Song D. The m6A methyltransferase METTL3 promotes osteosarcoma progression by regulating the m6A level of LEF1. Biochem Biophys Res Commun. 2019;516(3):719–25.
Zhou L, Yang C, Zhang N, Zhang X, Zhao T, Yu J. Silencing METTL3 inhibits the proliferation and invasion of osteosarcoma by regulating ATAD2. Biomed Pharmacother. 2020;125:109964.
Liu Z, Liu N, Huang Z, Wang W. METTL14 overexpression promotes osteosarcoma cell apoptosis and slows tumor progression via caspase 3 activation. Cancer Manag Res. 2020;12:12759–67.
Chen S, Li Y, Zhi S, Ding Z, Wang W, Peng Y, et al. WTAP promotes osteosarcoma tumorigenesis by repressing HMBOX1 expression in an m(6)A-dependent manner. Cell Death Dis. 2020;11(8):659.
Yuan Y, Yan G, He M, Lei H, Li L, Wang Y, et al. ALKBH5 suppresses tumor progression via an m(6)A-dependent epigenetic silencing of pre-miR-181b-1/YAP signaling axis in osteosarcoma. Cell Death Dis. 2021;12(1):60.
Cao J, Tang Z, Su Z. Long non-coding RNA LINC01426 facilitates glioblastoma progression via sponging miR-345-3p and upregulation of VAMP8. Cancer Cell Int. 2020;20:327.
Liu Z, Pan L, Yan X, Duan X. The long noncoding RNA DLGAP1-AS2 facilitates cholangiocarcinoma progression via miR-505 and GALNT10. FEBS Open Bio. 2021;11(2):413–22.
Miao W, Li N, Gu B, Yi G, Su Z, Cheng H. LncRNA DLGAP1-AS2 modulates glioma development by up-regulating YAP1 expression. J Biochem. 2020;167(4):411–8.
Dai J, Wang B, Zhao Y, Zuo X, Cui H, Chen X, et al. Long noncoding RNA LINC01426 sequesters microRNA-519d-5p to promote non-small cell lung cancer progression by increasing ETS1 expression. Cancer Manag Res. 2020;12:12697–708.
We are sincerely acknowledging the contributions from the TARGET.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Zhang, P., Xu, K., Wang, J. et al. Identification of N6-methylandenosine related LncRNAs biomarkers associated with the overall survival of osteosarcoma. BMC Cancer 21, 1285 (2021). https://doi.org/10.1186/s12885-021-09011-z