- Open Access
RNA adenosine modifications related to prognosis and immune infiltration in osteosarcoma
Journal of Translational Medicine volume 20, Article number: 228 (2022)
RNA adenosine modifications, which are primarily mediated by “writer” enzymes (RMWs), play a key role in epigenetic regulation in various biological processes, including tumorigenesis. However, the expression and prognostic role of these genes in osteosarcoma (OS) remain unclear.
Univariate and multivariate Cox analyses were used to construct the RMW signature for OS using Target datasets. RMW expression in OS tissue was detected by qPCR analysis. Xcell and GSVA were used to determine the relationship between RMWs and immune infiltration. The DGIdb and CMap databases were used for drug prediction. In vivo and in vitro experiments showed that strophanthidin elicited antitumor activity against OS.
A 3-RMW (CSTF2, ADAR and WTAP) prognostic signature in OS was constructed using the Target dataset and verified using GEO datasets and 63 independent OS tissues via qPCR analysis. High-risk OS patients had poor overall survival, and the prognostic signature was an independent prognostic factor for OS. Functional studies showed that tumour-, metabolism-, cell cycle- and immune-related pathways were related to high risk. Next, we found that RMW-derived high-risk patients exhibited increased infiltration of M2 macrophages and cDCs. Furthermore, we predicted the potential drugs for OS using the DGIdb and CMap databases. In vivo and in vitro experiments showed that strophanthidin elicited antitumor activity against OS by repressing cell growth and inducing cell cycle arrest at the G1 phase.
The 3-RWM-based prognostic signature established in this study is a novel gene signature associated with immune infiltration, and strophanthidin was identified as a candidate therapy for OS by repressing OS cell growth and the cell cycle.
Human osteosarcoma (OS) is the most frequent aggressive bone cancer in children and adolescents . Despite improvements in multimodal therapies, the prognosis of OS remains poor (20–30%) due to the delay in diagnosis and the development of metastasis . Therefore, it is urgent to reveal novel biomarkers and ensure effective OS treatment.
Although RNA modification has been recognized for more than half a century, its cell biology remains largely unexplored . RNA adenosine modification is the most common type of RNA modification, including m1A and m6A modification. RNA editing includes adenosine-to-inosine (A-to-I) editing and alternative polyadenylation (APA) . RNA methylation accounts for approximately 60% of RNA modifications, whereas m1A and m6A modifications are common and abundant in RNA methylation modification by the methylation of the adenine base. These modifications regulate RNA stability and translation . m1A and m6A modifications play key roles in various cellular processes, thus leading to a variety of diseases, including cancer [5, 6]. Polyadenylation (APA) is a phenomenon in which one gene contains multiple polyadenylation (pA) sites to produce transcript isoforms at either the 3′-untranslated region (UTR) or coding regions . Another common posttranscriptional mechanism is RNA editing, which alters gene expression by regulating the nucleotides of transcripts . The most widespread type of RNA editing is A-to-I, which is catalysed by ADAR enzymes and alters the coding, folding, splicing, or transport of transcripts. Dysregulated A-to-I editing can lead to various diseases, including tumorigenesis . Studies have shown that many of these modifications could be related to a complex network that intersects each other . Recently, these RNA modifications have been reported as biomarkers and play key roles in the progression of tumour prognosis [11, 12]. However, their role in OS remains limited.
In this work, we identified RMW patterns and RNA modification-related DEG patterns that were associated with the prognosis of OS and related to immune infiltration. Next, we constructed and verified an RMW signature of OS, which was not only associated with prognosis but also with immune infiltration and the cell cycle. Finally, strophanthidin was predicted and verified as an effective drug for OS.
Materials and methods
The expression of OS patients with clinical information was obtained from the GSE21257 and Target datasets. Eighty-five OS patients and 53 OS patients with survival follow-up information were collected from the Target and GSE21257 datasets, respectively.
Consensus clustering analysis
To determine the RNA modification patterns in OS, 26 RMWs were selected for consensus clustering analysis using the Consensus Cluster Plus package, and k = 2 seemed to be the most appropriate selection.
The DEGs identified between the two RMW patterns were analysed using the “limma” package with adjusted p < 0.05. Survival-related DEGs were analysed using univariate Cox regression analysis and subsequently selected for a consensus clustering analysis using the Consensus Cluster Plus package, and k = 2 seemed to be the most appropriate selection.
The survival analysis of patients from the two patterns was analysed using the packages “survminer” and “survival”.
The relationship between patterns and clinical characteristics was analysed using the “RColorBrewer” package.
The RMW prognostic signature
The survival-related RMWs were analysed using univariate Cox regression analysis. The prognostic signature was analysed using lasso regression and multivariate Cox regression analysis. The risk score was calculated as follows:
Risk score = ∑iCoefficient RMWs *Expression RMWs.
The OS patients were divided into high- and low-risk groups based on the best cut-off value of the risk score from the TCGA dataset. The R “survival” package was used to assess the relationship between survival and the RMW signature. Univariate and multivariate Cox regressions were used to reveal the independent risk factors for OS using R software. ROC curves were employed to reveal the prognostic value using the R software “ROC” package. The nomogram and calibration curve were analysed using the R “rms” package.
RMW signature-related pathways
GSEA and GSVA were used to reveal the risk signature-related pathways using the Target dataset. GSEA software (Cambridge, MA, United States) was used to describe the pathways related to RMWs. For GSVA, the signalling pathway alterations between the high- and low-risk groups were analysed using the “GSVA” R package.
Immune infiltration in OS tissues
The immune infiltration of OS tissues from the Target database was evaluated using xCell (https://xcell.ucsf.edu/). Next, the immune-related signalling and immune cell infiltration in different RMW patterns, DEG patterns and high/low risk groups were analysed using the R “ggpubr” package.
Clinical samples and qPCR analysis
OS tissues were obtained from 63 OS patients (September 2018 to January 2020) from the Third Xiangya Hospital, Central South University. Our study was approved by the ethics committee of the Third Xiangya Hospital, Central South University.
Total RNA was collected from OS tissues and cells using TRIpure reagent (BioTek, VT, USA) and then reverse transcribed using a HiScript Q RT SuperMix kit (Vazyme, Nanjing, China). Then, qRT–PCR was performed to assess mRNA expression using SYBR Green Master Mix (CWBIO, Jiangsu, China) in an Applied Biosystems QuantStudio 3 Real-Time PCR System (Thermo Fisher Scientific, MA, USA) as previously described . The nucleotide sequences of primers for CSTF2, ADAR, WTAP and β-actin are listed in Table 1.
Human OS cell lines (MG63, HOS and U2OS), L02 and HUVEC cell lines were purchased from American Type Culture Collection (ATCC) (Manassas, VA, USA). All cells were cultured as previously described .
Cell viability assay
The MTT assay was used to detect cell viability. In short, the cells were cultured in 96-well plates overnight and then treated with strophanthidin (0.001 to 100 μM) (Santa Cruz, sc-215914A)/PBS for 24 or 48 h. Then, cell viability was detected using an MTT assay. The IC50 value was calculated using GraphPad Prism.
OS cells were incubated with strophanthidin for 48 h. Then, ~ 1000 cells were seeded in 6-well plates and cultured for 10 days for the colony forming assay.
Flow cytometric analysis
OS cells were incubated with strophanthidin for 48 h. Then, flow cytometry was used for the cell cycle assay as previously described .
Western blot analysis
OS cells were incubated with strophanthidin for 48 h. Then, the protein was collected using RIPA buffer with a protease inhibitor cocktail (Roche Applied Science, Indianapolis, USA). The proteins were separated and transferred to PVDF membranes (Millipore, Massachusetts, USA). Primary antibodies against GAPDH (ab8245, Abcam), CDK4 (ab108357, Abcam), CCND1 (ab16663, Abcam) and p53 (ab32389, Abcam) and secondary antibodies were used to detect target proteins.
Four- to six-week-old female BALB/c nude mice were maintained in SPF conditions. Nude mice were injected with HOS cells (1 × 106) into the left scapula subcutaneously. For the control group, the mice were administered 200 μl of saline po/day. After 10 days, nude mice were randomly assigned to 2 groups. For the treatment group, mice were treated with strophanthidin (0.5 mg/kg/day). The body weights and the tumour volume were measured every 4 days. After 24 days, all nude mice were sacrificed, and the tumour weights were detected. All animal experiments were approved by the ethics committee of Third Xiangya Hospital, Central South University.
All data from Target and GEO were analysed by R-3.6.1. The differential RMWs were evaluated using the limma and heatmap packages. Student’s t test was used for two groups of in vitro experiments, and chi-squared (χ2) tests were used for in vivo experiments. p < 0.05 was considered statistically significant.
RMW patterns associated with immune infiltration
As noted in previous studies, a total of 26 RNA modification “writers” (RMWs, including 7 m6A “writers”, 4 m1A “writers”, 3 A to I “writers”, and 12 APA “writers”) were collected to analyse the prognostic role in OS. Twenty-six RMWs were observed in OS tissues from the target dataset and GSE21257 dataset, as shown in Additional file 1: Fig. S1A and B. Next, 26 RMWs were used for consensus clustering analysis, and 2 clusters were identified (Fig. 1A). The survival analysis showed that patients in Cluster B had poor overall survival (Fig. 1B). The heatmap revealed RMW expression and clinical characteristics in OS (Fig. 1C). GSVA and xCell analysis showed that Cluster B showed higher immune-related signals and immune cell infiltration (Fig. 1D and E). The chi-square test results showed that RMW patterns were associated with survival state and metastatic characteristics (Fig. 1F).
DEGs patterns associated with immune infiltration
To further describe the functional role of the RMW patterns above, 316 DEGs were identified in Cluster B compared to Cluster A. GO and KEGG enrichment analyses showed that the DEGs were enriched in cancer-related signalling pathways and apoptosis-, senescence- and immune-related pathways (Additional file 1: Fig. S1C). Next, 16 of 316 DEGs were associated with overall survival using Cox analysis. The 16 DEGs were used for consensus clustering analysis, and 2 DEG clusters were identified (Fig. 2A). The survival analysis showed that patients in DEG Cluster 2 had poor overall survival (Fig. 2B). The heatmap revealed the expression of the 16 risk-related DEGs and clinical characteristics in OS (Fig. 2C). We found that most patients in Cluster A were in Cluster 1, and most patients in Cluster B were in Cluster 2. Consistent with the RMW pattern, Cluster 2 showed higher immune-related signals and immune cell infiltration (Fig. 2D and E), and the DEG patterns were associated with survival state and metastatic characteristics (Fig. 2F).
RMW-associated prognostic signature in OS
To reveal the prognosis of RMWs in OS, univariate and multivariate Cox analyses were used to construct the RMW signature in OS using Target datasets. Three RMWs were related to the overall survival of OS patients (Fig. 3A). Lasso regression analysis and multivariate Cox regression analysis were used to construct the risk model (Fig. 3B). The risk score was calculated as follows: risk score = (0.03583 × CSTF2) + (-0.011308 × ADAR) + (-0.01465 × WTAP) (Fig. 3C). The AUC value of the RMW risk model was 0.860, and the best cut-off value was used for the low/high-risk group (Fig. 3D). The survival results revealed the poor prognosis of OS patients with high expression of CSTF2 and low expression of ADAR and WTAP in both TCGA and GEO datasets (Additional file 1: Fig. S2).
Seventeen patients and 68 patients were included in the high- and low-risk groups, respectively (Fig. 3E). The PCA demonstrated the markedly different distribution of the OS patients with high/low risk (Fig. 3F). OS patients with high risk showed a poorer prognosis than those with low risk (Fig. 3G). The AUC values of the RMWs were 0.860, 0.839, and 0.790 for survival times of 1 year, 2 years and 3 years, respectively, in the Target dataset (Fig. 3H). The results indicate that the model was convincing, specific, and sensitive.
The relationship between RMWs signature and clinical characteristics
Next, we analysed the independent prognostic factor in OS. We found that Metastatic (P = 0.004), Primary_site (P = 0.022) and risk score (P = 0.003) were independent prognostic factors in the Target dataset (Fig. 4A and B). The ROC curve showed that the AUC values of Metastatic, RMWs, and Primary_site were 0.912, 0.860, and 0.544, respectively (Fig. 4C). We next analysed the relationship between the risk score and clinical characteristics. As shown in Fig. 4D, the risk score was higher in the metastatic group than in the nonmetastatic group. The relationships between risk RMWs and clinical characteristics in the target datasets are shown in Additional file 1: Fig. S3.
Next, we analysed survival in different clinical subgroups and low/high signatures. We found that a high-risk score was related to poor prognosis in age, metastasis and sex (Fig. 5A). Finally, the nomogram and calibration curves were applied to estimate the survival probabilities in OS (Fig. 5B), implying that the RMW signature provided an accurate prognosis of OS.
Validation of the RMW signature in the GEO dataset
We next validated the RMW signature in GSE21257 with 53 OS patients. The risk scores of OS patients in GSE21257 were calculated according to the formula results. Figure 6A shows 20 patients and 33 patients in the high- and low-risk groups, respectively. The PCA demonstrated the markedly different distribution of OS patients in the high/low-risk group (Fig. 6B). OS patients with high risk showed a poorer prognosis than those with low risk (Fig. 6C, p = 0.005117). The AUC values of the RMWs were 0.791, 0.794, and 0.702 for survival times of 1 year, 2 years and 3 years, respectively, in the GSE21257 dataset (Fig. 6D).
RMW expression in OS tissue detected by qPCR analysis
To verify our previous results, we collected 63 OS tissues to detect RMW expression using qPCR. Table 2 shows the clinical information of the OS patients. As shown in Fig. 6E, the survival analysis demonstrated that OS patients with high expression of CSTF2 and OS patients with low expression of ADAR and WTAP showed poor prognoses. Moreover, OS patients with high risk showed poor prognoses, and the AUC values at 1, 2, and 3 years were 0.869, 0.823 and 0.890, respectively. (Fig. 6F). CSTF2, ADAR and WTAP expression was also related to OS metastasis. The CSTF2 expression level was significantly higher and ADAR and WTAP expression levels were significantly lower in OS patients who died (Fig. 6G).
RMW signature-related immune infiltration in OS
To further reveal the potential molecular mechanisms caused by RMW-mediated risk, we analysed the DEGs and related signalling pathways in high- and low-risk patients. A total of 1699 DEGs (494 upregulated and 1205 downregulated) were detected between the high- and low-risk groups. GO analysis showed that DEGs participate in the regulation of the catabolic process, cytokine signalling in the immune system, and cell cycle process (Fig. 7A). Next, GSEA and GSVA were used to reveal the risk signature-related pathways. As shown in Fig. 7B, the immune-related signalling pathway and cell cycle-related signalling pathway were regulated by high/low risk. The GSVA results showed that cancer- and metabolism-related and immune-related signalling pathways were differentially expressed in the low/high-risk groups (Fig. 7C). xCell was used to analyse the immune infiltration in OS, and 34 immune cells and 3 immune-related scores are shown in Fig. 7D. The results showed that the infiltration of CD8 naïve T cells, M1 macrophages, M2 macrophages and cDCs was evidently increased in the high-risk group. Moreover, Pearson analysis revealed a correlation between the risk score and immune cell infiltration (Additional file 1: Fig. S4). However, Fig. 7E shows that the immune checkpoint blockade (ICB) molecules CD44, CD200, CD276 and KIR3DL1 were differentially expressed between the high/low groups.
Potential drugs for OS patients
To reveal the RMW-related hub genes, 12 prognosis-related genes were identified from 699 DEGs (Additional file 1: Table S1), and cancer-related transcription factors related to prognostic genes were identified, as shown in Additional file 1: Table S2 (R > 0.4 and p < 0.05). The 10 prognostic genes and 25 related transcription factor regulatory networks are shown in Fig. 8A. Then, the 35 hub genes were submitted to CMAP and DGldb for drug prediction. Fifty-five drugs were identified by the CMap database (Additional file 1: Table S3), and 1130 drugs were identified by the DGldb database. In total 10 drugs overlapped, including dyclonine, resveratrol, zaprinast, abamectin, thioguanosine, cefotaxime, mesoridazine, clioquinol, physostigmine and strophanthidin (Fig. 8B). The relationship between drugs and targets is described in Fig. 8C. Strophanthidin, a cardiac glycoside, was reported to be a promising anticancer agent by regulating the cell cycles of breast, lung and liver cancer cells . It is also an immune cell activator that activates ROR T receptors . Therefore, we selected strophanthidin for further functional analysis.
Strophanthidin is an effective anti-OS drug in vitro and in vivo
We next detected the effects of strophanthidin on cell proliferation in OS cells (MG63, HOS and U2OS). As shown in Fig. 9A, strophanthidin inhibited the proliferation ability of OS (MG63, HOS and U2OS) cells in a dose-dependent manner. The IC50 value of strophanthidin in HOS cells was much lower than those in MG63 and U2OS cells. However, 1 μM strophanthidin did not have significant effects on L02 cells or HUVECs (Fig. 9B). The colony assay also showed that 1 μM strophanthidin significantly reduced the number of colons in MG63, HOS and U2OS cells (Fig. 9C). We next analysed the role of strophanthidin in the OS cell cycle. As shown in Fig. 9D, strophanthidin (1 μM) treatment induced G1 arrest in OS cells. Furthermore, we detected the expression of cell cycle-related proteins in OS cells. To further reveal the inhibition of strophanthidin on OS growth, 0.5 mg/kg strophanthidin was used to treat OS mice. As shown in Fig. 9E-G, strophanthidin treatment reduced tumour growth in vivo.
OS is a commonly aggressive tumour in children and adolescents with poor survival rates. Recently, prognostic models for improving the prognosis of OS have attracted considerable attention. In this study, we constructed and verified the RMW patterns and prognostic signature in OS and shed light on strophanthidin as a novel therapeutic drug for OS patients.
RNA adenosine modification is the most common type of RNA modification and is involved in various pathophysiological processes, including tumorigenesis [17, 18]. m6A modification is an RNA modification with dynamic and reversible posttranscriptional characteristics that affects RNA processing, degradation and translation. Studies have shown that m6A modification levels are dysregulated in tumour tissues and are associated with the occurrence and progression of various tumours [19,20,21,22]. The “writers” play a key role in these RNA modifications, and some of the writers were reported as therapeutic targets of tumours. For example, METTL3, the m6A writer, is related to the prognosis and progression of acute myeloid leukaemia (AML), and METTL3 inhibition was reported as a therapeutic strategy for myeloid leukaemia . A-to-I RNA editing is another abundant RNA modification event affecting adenosines in mammals, playing a critical role in the pathogenesis of various tumours . Xu et al. revealed the role of A-to-I-edited miRNA cancer progression and highlighted the translational potential of edited miRNAs as a new class of cancer therapeutics . Moreover, other adenine-related RNA modifications (m1A methylation and APA) have recently also emerged as key players in cancer pathophysiology by regulating the expression of cancer-related genes . In OS, m6A modification of YAP is involved in the progression of OS . The relationship between other adenine-related RNA modifications and OS remains unknown. RNA modification plays an important role in various biological processes through interactions with various “writers” (RMWs). Recent studies have focused on RNA modification signatures to evaluate the prognosis of various tumours. Our study explored RMW patterns and RNA modification DEG patterns that were associated with the survival of OS patients. Next, the RMW prognostic signature was constructed and identified as an independent prognostic factor for OS. In conclusion, these results indicated the important role of RMWs in OS. Although previous studies revealed several prognostic models (including ferroptosis-related gene signatures and autophagy-related prognosis models) for OS [27,28,29,30], this is the first study of the RWM-related prognosis model and RMP patterns for OS. Moreover, the ROC of the RMW model was considerably increased compared with that of the ferroptosis-related gene signature (ROC = 8.06) and autophagy-related prognosis model (ROC = 8.38). These results indicated that the RMW signature is a more effective prognostic model for OS.
The immune microenvironment (TME) is related to the prognosis of various cancers, including OS [31,32,33]. The recruitment of stromal cells was reported to be associated with a poor patient prognosis. Infiltrating myeloid cells are common and important components of the TME, and they are reported to drive OS progression . Tumour-infiltrating myeloid cells include monocytes, dendritic cells (DCs), tumour-associated macrophages (TAMs), and neutrophils. In OS lesions, monocytes and macrophages are the most common myeloid cells, and DCs account for < 5% of myeloid cells . M2 TAMs are anti-inflammatory macrophages and are often associated with a worse prognosis . M1 TAMs are antitumor immune cells that express proinflammatory cytokines. Anti-PD-1 is an effective therapeutic strategy for OS due to the regulation of the infiltration of M1 and M2 macrophages in tumour tissue . Recently, DCs were identified as therapeutic targets of immunotherapy due to their powerful antigen-presenting features. Zhou et al. showed that the infiltration of CD1c+ DCs was significantly increased in metastatic OS . In this study, we used xCell to reveal the correlations between the risk signature and immune infiltration. We found that M1 macrophages, M2 macrophages and cDCs were significantly increased in OS patients at high risk. Recent studies shed new light on the regulation of RNA modification on the immune system . Furthermore, RNA adenosine modifications are associated with immunoregulation in tumour tissues . METTL3-mediated m6A modification promotes tumour growth and metastasis by macrophage reprogramming . METTL3-mediated m6A modification also affects resistance to chemotherapy by regulating M2-TAM infiltration . METTL3-mediated m6A methylation regulates dendritic cell activation . These results highlighted that RNA adenosine modifications might affect the prognosis of OS partly by regulating the immune microenvironment.
DGIdb and CMap, drug prediction databases, are widely used to screen drugs that can regulate certain target genes. In this study, we utilized these two databases to obtain the possible RMW-related drugs for OS. Ten drugs were identified, including dyclonine, resveratrol, zaprinast, abamectin, thioguanosine, cefotaxime, mesoridazine, clioquinol, physostigmine and strophanthidin. Dyclonine, an ALDH3A1 inhibitor, is widely used as a topical antipruritic agent . Recent studies reported that dyclonine enhances MG132-induced cell cytotoxicity in breast cancer . Resveratrol, an oestrogenic compound, was reported as a treatment for osteoporosis with anti-inflammatory and antioxidant properties . Recent studies also revealed the therapeutic potential of resveratrol on various tumours, including OS [45,46,47]. Zaprinast, a synthetic GPR35 agonist, was reported to rescue OVX-induced bone . It also exhibits antitumor activity in colon cancer, lung cancer, and melanoma [49,50,51]. Abamectin, an anthelmintic agent in animals, induces oxidative stress in cerebral and hepatic tissues . Thioguanosine has antitumor activities against hepatocellular carcinoma . Mesoridazine, an antipsychotic drug, has cardiac conduction side effects with fatal consequences in patients . Clioquinol has been reported as an antitumor drug for OS . Physostigmine has been reported to repress pancreatic cancer growth with anticholinergic toxicity . Strophanthidin, a cardiac glycoside, was reported to be a promising anticancer agent by regulating the cell cycles of breast, lung and liver cancer cells . It is also an immune cell activator that activates ROR T receptors . Therefore, we selected strophanthidin for further functional analysis. In this study, in vivo and in vitro experiments showed that strophanthidin elicited antitumor activity against OS by repressing cell proliferation and the cell cycle. As shown in Fig. 8, Smad3 was predicted to be the target of strophanthidin. Smad3 is a key mediator of TGF-β signalling, an important signalling pathway in tumorigenesis . Smad3 was also reported to be involved in the regulation of the cell cycle via the induction of CDK inhibitors [57, 58]. Based on the above research, we hypothesize that strophanthidin regulates the cell cycle of OS cells by targeting Smad3.
In summary, we analysed RNA modification “writers” (RMWs) in OS and revealed RMW patterns and signatures associated with the cell cycle and immune infiltration. Strophanthidin was identified and verified as a candidate therapy for OS by repressing the growth and cell cycle of OS cells.
Availability of data and materials
The data used to support the findings of this study are available either online or from the corresponding author upon request.
Matsuoka K, Bakiri L, Wolff LI, Linder M, Mikels-Vigdal A, Patiño-García A, Lecanda F, Hartmann C, Sibilia M, Wagner EF. Wnt signaling and Loxl2 promote aggressive osteosarcoma. Cell Res. 2020;30:885–901. https://doi.org/10.1038/s41422-020-0370-1.
Isakoff MS, Bielack SS, Meltzer P, Gorlick R. Osteosarcoma: current treatment and a collaborative pathway to success. J Clin Oncol. 2015;33:3029–35. https://doi.org/10.1200/jco.2014.59.4895.
Roundtree IA, Evans ME, Pan T, He C. Dynamic RNA modifications in gene expression regulation. Cell. 2017;169:1187–200. https://doi.org/10.1016/j.cell.2017.05.045.
Lee Y, Choe J, Park OH, Kim YK. Molecular mechanisms driving mRNA degradation by m(6)A modification. Trends Genet. 2020;36:177–88. https://doi.org/10.1016/j.tig.2019.12.007.
Wang T, Kong S, Tao M, Ju S. The potential role of RNA N6-methyladenosine in cancer progression. Mol Cancer. 2020;19:88. https://doi.org/10.1186/s12943-020-01204-7.
Chen XY, Zhang J, Zhu JS. The role of m(6)A RNA methylation in human cancer. Mol Cancer. 2019;18:103. https://doi.org/10.1186/s12943-019-1033-z.
Gruber AJ, Zavolan M. Alternative cleavage and polyadenylation in health and disease. Nat Rev Genet. 2019;20:599–614. https://doi.org/10.1038/s41576-019-0145-z.
Baysal BE, Sharma S, Hashemikhabir S, Janga SC. RNA editing in pathogenesis of cancer. Cancer Res. 2017;77:3733–9. https://doi.org/10.1158/0008-5472.can-17-0520.
Jain M, Jantsch MF, Licht K. The editor’s I on disease development. Trends Genet. 2019;35:903–13. https://doi.org/10.1016/j.tig.2019.09.004.
Tassinari V, Cesarini V, Tomaselli S, Ianniello Z, Silvestris DA, Ginistrelli LC, Martini M, De Angelis B, De Luca G, Vitiani LR, et al. ADAR1 is a new target of METTL3 and plays a pro-oncogenic role in glioblastoma by an editing-independent mechanism. Genome Biol. 2021;22:51. https://doi.org/10.1186/s13059-021-02271-9.
Li Y, Gu J, Xu F, Zhu Q, Chen Y, Ge D, Lu C. Molecular characterization, biological function, tumor microenvironment association and clinical significance of m6A regulators in lung adenocarcinoma. Brief Bioinform. 2020. https://doi.org/10.1093/bib/bbaa225.
Xu F, Huang X, Li Y, Chen Y, Lin L. m(6)A-related lncRNAs are potential biomarkers for predicting prognoses and immune responses in patients with LUAD. Mol Ther Nucleic Acids. 2021;24:780–91. https://doi.org/10.1016/j.omtn.2021.04.003.
Iliopoulos D, Bimpaki EI, Nesterova M, Stratakis CA. MicroRNA signature of primary pigmented nodular adrenocortical disease: clinical correlations and regulation of Wnt signaling. Cancer Res. 2009;69:3278–82. https://doi.org/10.1158/0008-5472.can-09-0155.
Pan Z, Cheng DD, Wei XJ, Li SJ, Guo H, Yang QC. Chitooligosaccharides inhibit tumor progression and induce autophagy through the activation of the p53/mTOR pathway in osteosarcoma. Carbohydr Polym. 2021;258: 117596. https://doi.org/10.1016/j.carbpol.2020.117596.
Reddy D, Ghosh P, Kumavath R. Strophanthidin attenuates MAPK, PI3K/AKT/mTOR, and Wnt/β-catenin signaling pathways in human cancers. Front Oncol. 2019;9:1469. https://doi.org/10.3389/fonc.2019.01469.
Karaś K, Sałkowska A, Walczak-Drzewiecka A, Ryba K, Dastych J, Bachorz RA, Ratajewski M. The cardenolides strophanthidin, digoxigenin and dihydroouabain act as activators of the human RORγ/RORγT receptors. Toxicol Lett. 2018;295:314–24. https://doi.org/10.1016/j.toxlet.2018.07.002.
McFadden MJ, Horner SM. N(6)-methyladenosine regulates host responses to viral infection. Trends Biochem Sci. 2021;46:366–77. https://doi.org/10.1016/j.tibs.2020.11.008.
Jiang X, Liu B, Nie Z, Duan L, Xiong Q, Jin Z, Yang C, Chen Y. The role of m6A modification in the biological functions and diseases. Signal Transduct Target Ther. 2021;6:74. https://doi.org/10.1038/s41392-020-00450-x.
Chen M, Wong CM. The emerging roles of N6-methyladenosine (m6A) deregulation in liver carcinogenesis. Mol Cancer. 2020;19:44. https://doi.org/10.1186/s12943-020-01172-y.
Wang M, Liu J, Zhao Y, He R, Xu X, Guo X, Li X, Xu S, Miao J, Guo J, et al. Upregulation of METTL14 mediates the elevation of PERP mRNA N(6) adenosine methylation promoting the growth and metastasis of pancreatic cancer. Mol Cancer. 2020;19:130. https://doi.org/10.1186/s12943-020-01249-8.
Chen X, Xu M, Xu X, Zeng K, Liu X, Pan B, Li C, Sun L, Qin J, Xu T, et al. METTL14-mediated N6-methyladenosine modification of SOX4 mRNA inhibits tumor metastasis in colorectal cancer. Mol Cancer. 2020;19:106. https://doi.org/10.1186/s12943-020-01220-7.
Yuan Y, Yan G, He M, Lei H, Li L, Wang Y, He X, Li G, Wang Q, Gao 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:60. https://doi.org/10.1038/s41419-020-03315-x.
Yankova E, Blackaby W, Albertella M, Rak J, De Braekeleer E, Tsagkogeorga G, Pilka ES, Aspris D, Leggate D, Hendrick AG, et al. Small molecule inhibition of METTL3 as a strategy against myeloid leukaemia. Nature. 2021. https://doi.org/10.1038/s41586-021-03536-w.
Han J, An O, Hong H, Chan THM, Song Y, Shen H, Tang SJ, Lin JS, Ng VHE, Tay DJT, et al. Suppression of adenosine-to-inosine (A-to-I) RNA editome by death associated protein 3 (DAP3) promotes cancer progression. Sci Adv. 2020;6:eaba5136. https://doi.org/10.1126/sciadv.aba5136.
Xu X, Wang Y, Mojumdar K, Zhou Z, Jeong KJ, Mangala LS, Yu S, Tsang YH, Rodriguez-Aguayo C, Lu Y, et al. A-to-I-edited miRNA-379-5p inhibits cancer cell proliferation through CD97-induced apoptosis. J Clin Invest. 2019;129:5343–56. https://doi.org/10.1172/jci123396.
Han X, Wang M, Zhao YL, Yang Y, Yang YG. RNA methylations in human cancers. Semin Cancer Biol. 2020. https://doi.org/10.1016/j.semcancer.2020.11.007.
Wu G, Zhang M. A novel risk score model based on eight genes and a nomogram for predicting overall survival of patients with osteosarcoma. BMC Cancer. 2020;20:456. https://doi.org/10.1186/s12885-020-06741-4.
Wan Y, Qu N, Yang Y, Ma J, Li Z, Zhang Z. Identification of a 3-gene signature based on differentially expressed invasion genes related to cancer molecular subtypes to predict the prognosis of osteosarcoma patients. Bioengineered. 2021;12:5916–31. https://doi.org/10.1080/21655979.2021.1971919.
Li J, Tang X, Du Y, Dong J, Zhao Z, Hu H, Song T, Guo J, Wang Y, Xu T, et al. Establishment of an autophagy-related clinical prognosis model for predicting the overall survival of osteosarcoma. Biomed Res Int. 2021;2021:5428425. https://doi.org/10.1155/2021/5428425.
Lei T, Qian H, Lei P, Hu Y. Ferroptosis-related gene signature associates with immunity and predicts prognosis accurately in patients with osteosarcoma. Cancer Sci. 2021;112:4785–98. https://doi.org/10.1111/cas.15131.
Wu CC, Beird HC, Andrew Livingston J, Advani S, Mitra A, Cao S, Reuben A, Ingram D, Wang WL, Ju Z, et al. Immuno-genomic landscape of osteosarcoma. Nat Commun. 2020;11:1008. https://doi.org/10.1038/s41467-020-14646-w.
Chen T, Zhao L. Patrolling monocytes inhibit osteosarcoma metastasis to the lung. Aging (Albany NY). 2020;12:23004–16. https://doi.org/10.18632/aging.104041.
Wang Y, Wang Y, Xu C, Liu Y, Huang Z. Identification of novel tumor-microenvironment-regulating factor that facilitates tumor immune infiltration in colon cancer. Mol Ther Nucleic Acids. 2020;22:236–50. https://doi.org/10.1016/j.omtn.2020.08.029.
Kansara M, Thomson K, Pang P, Dutour A, Mirabello L, Acher F, Pin JP, Demicco EG, Yan J, Teng MWL, et al. Infiltrating myeloid cells drive osteosarcoma progression via GRM4 regulation of IL23. Cancer Discov. 2019;9:1511–9. https://doi.org/10.1158/2159-8290.cd-19-0154.
Zhou Y, Yang D, Yang Q, Lv X, Huang W, Zhou Z, Wang Y, Zhang Z, Yuan T, Ding X, et al. Single-cell RNA landscape of intratumoral heterogeneity and immunosuppressive microenvironment in advanced osteosarcoma. Nat Commun. 2020;11:6322. https://doi.org/10.1038/s41467-020-20059-6.
Buddingh EP, Kuijjer ML, Duim RA, Bürger H, Agelopoulos K, Myklebost O, Serra M, Mertens F, Hogendoorn PC, Lankester AC, Cleton-Jansen AM. Tumor-infiltrating macrophages are associated with metastasis suppression in high-grade osteosarcoma: a rationale for treatment with macrophage activating agents. Clin Cancer Res. 2011;17:2110–9. https://doi.org/10.1158/1078-0432.ccr-10-2047.
Dhupkar P, Gordon N, Stewart J, Kleinerman ES. Anti-PD-1 therapy redirects macrophages from an M2 to an M1 phenotype inducing regression of OS lung metastases. Cancer Med. 2018;7:2654–64. https://doi.org/10.1002/cam4.1518.
Shulman Z, Stern-Ginossar N. The RNA modification N(6)-methyladenosine as a novel regulator of the immune system. Nat Immunol. 2020;21:501–12. https://doi.org/10.1038/s41590-020-0650-4.
Yin H, Zhang X, Yang P, Zhang X, Peng Y, Li D, Yu Y, Wu Y, Wang Y, Zhang J, et al. RNA m6A methylation orchestrates cancer growth and metastasis via macrophage reprogramming. Nat Commun. 2021;12:1394. https://doi.org/10.1038/s41467-021-21514-8.
Lan H, Liu Y, Liu J, Wang X, Guan Z, Du J, Jin K. Tumor-associated macrophages promote oxaliplatin resistance via METTL3-mediated m(6)A of TRAF5 and necroptosis in colorectal cancer. Mol Pharm. 2021;18:1026–37. https://doi.org/10.1021/acs.molpharmaceut.0c00961.
Wang H, Hu X, Huang M, Liu J, Gu Y, Ma L, Zhou Q, Cao X. Mettl3-mediated mRNA m(6)A methylation promotes dendritic cell activation. Nat Commun. 1898;2019:10. https://doi.org/10.1038/s41467-019-09903-6.
Morginson WJ, Rich CO, Eskelson YD, Kirkman LW, Utterback M, Burton AM, Coletti JM. Dyclonine hydrochloride: a new topical antipruritic agent. Postgrad Med. 1956;19:605–7. https://doi.org/10.1080/00325481.1956.11708352.
Ju D, Wang X, Xie Y. Dyclonine and alverine citrate enhance the cytotoxic effects of proteasome inhibitor MG132 on breast cancer cells. Int J Mol Med. 2009;23:205–9.
Abhishek Shah A, Shah A, Lewis S, Ghate V, Saklani R, Narayana Kalkura S, Baby C, Kumar Singh P, Nayak Y, Chourasia MK. Cyclodextrin based bone regenerative inclusion complex for resveratrol in postmenopausal osteoporosis. Eur J Pharm Biopharm. 2021. https://doi.org/10.1016/j.ejpb.2021.07.008.
Xu N, Wang L, Fu S, Jiang B. Resveratrol is cytotoxic and acts synergistically with NF-κB inhibition in osteosarcoma MG-63 cells. Arch Med Sci. 2021;17:166–76. https://doi.org/10.5114/aoms.2020.100777.
Li Y, Bäckesjö CM, Haldosén LA, Lindgren U. Resveratrol inhibits proliferation and promotes apoptosis of osteosarcoma cells. Eur J Pharmacol. 2009;609:13–8. https://doi.org/10.1016/j.ejphar.2009.03.004.
Alkhalaf M, Jaffal S. Potent antiproliferative effects of resveratrol on human osteosarcoma SJSA1 cells: novel cellular mechanisms involving the ERKs/p53 cascade. Free Radic Biol Med. 2006;41:318–25. https://doi.org/10.1016/j.freeradbiomed.2006.04.019.
Zhang Y, Shi T, He Y. GPR35 regulates osteogenesis via the Wnt/GSK3β/β-catenin signaling pathway. Biochem Biophys Res Commun. 2021;556:171–8. https://doi.org/10.1016/j.bbrc.2021.03.084.
Monz K, Maas-Kück K, Schumacher U, Schulz T, Hallmann R, Schnäker EM, Schneider SW, Prehm P. Inhibition of hyaluronan export attenuates cell migration and metastasis of human melanoma. J Cell Biochem. 2008;105:1260–6. https://doi.org/10.1002/jcb.21925.
Pitari GM, Di Guglielmo MD, Park J, Schulz S, Waldman SA. Guanylyl cyclase C agonists regulate progression through the cell cycle of human colon carcinoma cells. Proc Natl Acad Sci USA. 2001;98:7846–51. https://doi.org/10.1073/pnas.141124698.
Oguri T, Achiwa H, Sato S, Bessho Y, Takano Y, Miyazaki M, Muramatsu H, Maeda H, Niimi T, Ueda R. The determinants of sensitivity and acquired resistance to gemcitabine differ in non-small cell lung cancer: a role of ABCC5 in gemcitabine sensitivity. Mol Cancer Ther. 2006;5:1800–6. https://doi.org/10.1158/1535-7163.mct-06-0025.
Radi AM, Mohammed ET, Abushouk AI, Aleya L, Abdel-Daim MM. The effects of abamectin on oxidative stress and gene expression in rat liver and brain tissues: modulation by sesame oil and ascorbic acid. Sci Total Environ. 2020;701: 134882. https://doi.org/10.1016/j.scitotenv.2019.134882.
Chen MH, Yang WL, Lin KT, Liu CH, Liu YW, Huang KW, Chang PM, Lai JM, Hsu CN, Chao KM, et al. Gene expression-based chemical genomics identifies potential therapeutic drugs in hepatocellular carcinoma. PLoS ONE. 2011;6: e27186. https://doi.org/10.1371/journal.pone.0027186.
Matesic DF, Abifadel DN, Garcia EL, Jann MW. Effect of thioridazine on gap junction intercellular communication in connexin 43-expressing cells. Cell Biol Toxicol. 2006;22:257–68. https://doi.org/10.1007/s10565-006-0047-7.
Pfitzinger PL, Fangmann L, Wang K, Demir E, Gürlevik E, Fleischmann-Mundt B, Brooks J, D’Haese JG, Teller S, Hecker A, et al. Indirect cholinergic activation slows down pancreatic cancer growth and tumor-associated inflammation. J Exp Clin Cancer Res. 2020;39:289. https://doi.org/10.1186/s13046-020-01796-4.
Li W, Zhou J, Zhang Y, Zhang J, Li X, Yan Q, Han J, Hu F. Echinacoside exerts anti-tumor activity via the miR-503-3p/TGF-β1/Smad aixs in liver cancer. Cancer Cell Int. 2021;21:304. https://doi.org/10.1186/s12935-021-01890-3.
Lai W, Tang Y, Huang XR, Ming-Kuen Tang P, Xu A, Szalai AJ, Lou TQ, Lan HY. C-reactive protein promotes acute kidney injury via Smad3-dependent inhibition of CDK2/cyclin E. Kidney Int. 2016;90:610–26. https://doi.org/10.1016/j.kint.2016.06.010.
Hur J, Choi JI, Lee H, Nham P, Kim TW, Chae CW, Yun JY, Kang JA, Kang J, Lee SE, et al. CD82/KAI1 maintains the dormancy of long-term hematopoietic stem cells through interaction with DARC-expressing macrophages. Cell Stem Cell. 2016;18:508–21. https://doi.org/10.1016/j.stem.2016.01.013.
We appreciate the researchers and patients who participated in the Target and GEO datasets.
This work was supported by the National Natural Science Foundation of China (Grant Nos. 82172594, 81772866), the Wisdom Accumulation and Talent Cultivation Project of the Third Xiangya Hospital of Central South University (YX202001), the Key Research and Development Program of Hunan Province (Grant No. 2022SK2033), Clinical Research Center of Hunan Province (Grant No. 2020 SK4017) the National Key Research and Development Program of China (Grant No. 2018YFC1004800), and the National Undergraduate Training Program for Innovation and Entrepreneurship of the Third Xiangya Hospital of Central South University (2021105330209).
Ethics approval and consent to participate
Consent for publication
All authors provide consent for publication.
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Figure S1. A, Heatmap of RMWs in OS tissue from the Target dataset and GSE21257 dataset. B, GO and KEGG analyses of DEGs. Figure S2. Survival analysis of OS patients with differential expression of CSTF2, ADAR and WTAP in both TCGA and GEO datasets. Figure S3. The correlation between risk RMW expression and clinical characteristics. Figure S4. The relationship between the risk score and immune infiltration was analysed using Pearson’s correlation analysis. Table S1. Univariate Cox analysis revealed prognosis-related DEGs in the high- and low-risk groups. Table S2. The relationship between TF and prognosis-related DEGs. Table S3. Drug prediction using the CMap database.
About this article
Cite this article
Chen, S., Zeng, J., Huang, L. et al. RNA adenosine modifications related to prognosis and immune infiltration in osteosarcoma. J Transl Med 20, 228 (2022). https://doi.org/10.1186/s12967-022-03415-6
- RNA adenosine modifications