Jianhao Li1,2†
Benchen Rao1,2†
Jie Yang3†Liwen Liu1,2Maoxin Huang4
Xin Liu1,2Guangying Cui1,2Chao Li5Qicai Han6Hao Yang6Xichun Cui2*
Ranran Sun1,2*- 1Precision Medicine Center, The First Affiliated Hospital of Zhengzhou University, Zhengzhou, China
- 2Key Laboratory of Clinical Medicine, The First Affiliated Hospital of Zhengzhou University, Zhengzhou, China
- 3Department of Orthopedics, Zhengzhou Central Affiliated Hospital to Zhengzhou University, Zhengzhou, China
- 4Dermatology Department, The First Affiliated Hospital of Zhengzhou University, Zhengzhou, China
- 5Department of Orthopedics, The First Affiliated Hospital of Zhengzhou University, Zhengzhou, China
- 6Department of Bone and Soft Tissue, The Affiliated Cancer Hospital of Zhengzhou University, Zhengzhou, China
Background: Osteosarcoma (OS) is the most common primary bone tumor. The disease has a poor prognosis due to the delay in the diagnosis and the development of metastasis. N6-Methyladenosine (m6A)-related regulators play an essential role in various tumors. In this study, a comprehensive analysis was conducted to elucidate the relationship between the expression profiles of m6A-related molecules and the clinical outcome of OS patients.
Materials and Methods: Public genome datasets and a tissue microarray (TMA) cohort were used to analyze the mRNA and protein expression levels of m6A regulators. Next, immunofluorescence (IF) analysis was used to determine the subcellular localization of m6A-related molecules. Kaplan–Meier and Cox regression analyses were performed to confirm the prognostic value of m6A-related molecules in OS. A comprehensive bioinformatic analysis was conducted to identify the potential molecular mechanisms mediated by m6A modification in OS.
Results: We found that m6A-related regulator expression was dysregulated in OS tissues, especially in metastatic tumor tissues. Low expression of METTL3, METTL14, and YTHDF2 and high expression of KIAA1429 and HNRNPA2B1 were significantly associated with poor prognosis in the TMA cohort. Simultaneously, the genome meta-cohort analysis revealed that low expression of FTO and METTL14 and high expression of METTL3, HNRNPA2B1, and YTHDF3 were associated with poor prognosis in OS. Cox regression analysis showed that HNRNPA2B1 might be an independent risk factor for OS. Bioinformatic analysis indicated that m6A regulators might be involved in OS progression through humoral immune response and cell cycle pathways.
Conclusion: M6A-related regulators are frequently dysregulated and correlate with metastasis and prognosis in OS. M6A-related regulators may serve as novel therapeutic targets and prognostic biomarkers for OS.
Introduction
Human osteosarcoma (OS) is one of the most common aggressive bone cancers, and it has higher incidence and mortality rates in teenagers; it frequently affects distal femur, proximal tibia, and humerus (1, 2). Despite significant progress in therapeutic strategies against OS over the past few decades, its prognosis remains poor due to the delay in the diagnosis and the development of metastasis (3). Thus, it is critical to clarify novel biomarkers and to ensure the effective treatment of OS patients.
N6-Methyl adenosine (m6A) is a posttranscriptional modification of RNA and is the most prevalent internal chemical modification of mRNAs (4). M6A network components have been well characterized into three subtypes: “writers,” “readers,” and “erasers” (5). The m6A modification is facilitated by writers and removed by erasers; moreover, it can recruit specific reader proteins (6). M6A-related regulators are involved in various physiological and pathological processes through the regulation of RNA stability, mRNA splicing and translation, and microRNA processing (7–9). Additionally, m6A is involved in the initiation and progression of cancers, including liver cancer, breast cancer, glioma, cervical cancer, colorectal cancer, and hepatoblastoma (10–15). However, the clinical value and potential mechanism of m6A-related regulators in OS are still unclear.
In this study, we evaluated the expression status and prognostic value of m6A-related proteins in OS based on public genome database analysis and tissue microarray (TMA) analysis. We found that dysregulated expression of several m6A-related proteins was frequently and closely associated with tumor metastasis and clinical outcomes in OS. Moreover, consensus clustering for m6A-related regulators was performed to identify the clusters with a better prognosis. Mechanistically, Gene Set Enrichment Analysis (GSEA) and Gene Set Variation Analysis (GSVA) indicated that the dysregulated m6A-related regulators might be associated with cell cycle and humoral immune response pathways. In conclusion, we revealed that m6A-related regulators play a critical role in the development of OS and proposed that m6A-related regulators could be potential therapeutic targets.
Materials and Methods
Datasets
The RNA-seq transcriptome data and clinical information from four OS datasets (GSE12865, GSE21257, GSE42352, GSE39055, TARGET-OS) were collected from the GEO (http://www.ncbi.nlm.nih.gov/geo/) and TARGET (https://ocg.cancer.gov/programs/target) databases, based on the inclusion and exclusion criteria of the GEO and TARGET databases. Inclusion criteria: datasets involving human osteosarcoma and expression profiling by array. Exclusion criteria: datasets with a sample size smaller than 10 (2, 11, 16, 17). The characteristics of the datasets are presented in Table S1. Subsequently, three gene expression matrix files (GSE21257, GSE39055, TARGET-OS) with survival follow-up data were merged into one meta-cohort. Next, the “sva” package of the R software was used to remove the batch effect (Figure S1). A total of 175 patients with survival follow-up information were divided into higher and lower expression groups using the best cutoff value of the log-rank test. Kaplan–Meier survival curves were used to analyze survival rates. Univariate and multivariate analyses were conducted based on the Cox model. The raw data were analyzed by BRB-array tools as previously described (18).
Tissue Samples
Studies involving human samples were approved (approval number: 2011-KY-047) by the Ethics Committee of the First Affiliated Hospital of Zhengzhou University (ZZU cohort), Zhengzhou Central Affiliated Hospital of Zhengzhou University, and the Affiliated Cancer Hospital of Zhengzhou University. All legal guardians of the children signed an informed consent. OS samples from the three hospitals were combined to create the TMA cohort containing 120 OS tissues and 65 surrounding non-tumorous tissues.
Immunohistochemistry (IHC)
Immunohistochemistry (IHC) of m6A-related genes was performed as described previously (19, 20). Briefly, 5-μm-thick TMA sections were deparaffinized and then treated with hydrogen peroxide to quench endogenous peroxidase activity. Subsequently, the sections were incubated overnight with anti-human m6A-related protein antibodies (1:200, Abcam, USA) at 4°C. Next, the immunoreactive cells were detected by SignalStain® DAB (CST, USA), then counterstained with Haematoxylin QS (Vector Laboratories). Two experienced pathologists, who were blinded to the clinicopathological data, evaluated the immunostaining samples separately. A semi-quantitative scoring system was established based on different staining intensities, and the proportion of positive cells was scored as follows: 0, none; 1+, <25%; 2+, 25–50%; 3+, 50–75%; and 4+, 75–100%. The staining intensity was scored as follows: 0, none; 1+, weak; 2+, medium; and 3+, strong. The total score was calculated by multiplying two subscores, and the samples with scores of 0–6 were regarded as low expression, whereas the other ratings were regarded as high expression during statistical analysis. Antibody information is listed in Table S2.
Cell Lines and Culture
The human OS cell lines U2-OS and KHOS-240S were purchased from ATCC (Manassas, USA). Cells were maintained at 37°C in a humidified atmosphere of 5% CO2 in Dulbecco's modified Eagle medium (DMEM) supplemented with 10% fetal bovine serum (FBS) (Gibco, New York, NY, USA) and 100 U/mL penicillin/streptomycin (Corning, New York, NY, USA). The STR reports for the cell lines U2-OS and KHOS-240S are presented in Supplemental Materials 1–2.
Immunofluorescence Assay
OS cells cultured in 24-well plates were fixed in 4% paraformaldehyde and permeabilized in 1% Triton X-100 phosphate-buffered saline (PBS). After blocking with 1% bovine serum albumin (BSA), the cells were incubated with primary antibodies at 4°C overnight and then with appropriate corresponding secondary antibodies (Jackson ImmunoResearch Inc., USA) at room temperature for 30 min. Finally, the nuclei were counterstained with DAPI (Beyotime, China), and images were obtained with a Zeiss Axio microscope (Zeiss, Oberkochen, Germany). Detailed information on the antibodies used in this study is listed in Table S2.
Statistical Analysis
All statistical analyses were performed using the R statistical package (R version 3.5.3) unless otherwise stated. The differences between the two independent groups were analyzed using Student's t test (unpaired, two-tailed) or Permutation test when there were fewer than three samples in either group (21). Kaplan–Meier overall survival analysis was performed with a log-rank test. Univariate Cox regression analysis was used to indicate the relationship between the different variables and survival. The correlation was evaluated using the two-tailed Pearson test. We clustered OS patients into different clusters with “ConsensusClusterPlus” (22). Additionally, 150 clinically actionable genes were obtained from a recent publication (23). Subsequently, the protein–protein interactions among m6A regulators and 150 clinically actionable genes were identified based on the STRING (https://string-db.org/) interaction database (24). Cytoscape software was used to visualize the interactions. GSVA was performed using the Bioconductor R package “GSVA” (25). GSEA was conducted using clusterProfiler, an R/Bioconductor package (26). In all cases, P < 0.05 was considered statistically significant expression is frequently dysregulated.
Results
M6A-Related Gene Expression Is Frequently Dysregulated in Osteosarcoma
To determine the significant biological function of m6A-related regulators in tumorigenesis and development, the expression pattern of m6A-related genes was analyzed through the GEO database (GSE42352) (Figure 1A). The mRNA expression levels of YTHDF2, YTHDF1, HNRNPC, FTO, METTL3, RBM15, HNRNPA2B1, and YTHDC1 were higher in OS cells than in normal cell lines. In addition, differential expression analysis of GSE12865 through a permutation test further confirmed that the expression of RBM15 was significantly upregulated in the tumor tissues (Figure 1B). Subsequently, the relationship between tumor metastasis and the expression level of m6A-related regulators was examined in GSE21257 and GSE42352 (Figures 1C,D). The overexpression of RBM15B, METTL14, and HNRNPA2B1 is significantly related to tumor metastasis. In conclusion, these results suggested that dysregulated m6A-related regulators were associated with tumorigenesis and metastasis in OS.
Figure 1. m6A-related regulators expression status in osteosarcoma. (A) The mRNA expression level of m6A-related regulators in normal and OS cell lines of GSE42352. (B) The mRNA expression level of m6A-related molecules in tumor and non-tumor tissues of GSE12865. (C,D) The correlation between the mRNA expression of m6A-related regulators and tumor metastasis.
Subcellular Location of m6A-Associated Proteins in Osteosarcoma Cell Lines
The diverse subcellular location of proteins may reflect on the different functions of m6A-related regulators in OS cells. Therefore, IF was performed to determine the subcellular location of m6A-related proteins in U2-OS and KHOS-240S cell lines. We found that all the m6A “writers” had intense nuclear and weak cytoplasmic staining in the two OS cell lines (Figure 2A). Moreover, m6A “readers” HNRNPC and HNRNPA2B1 were detected only in the nucleus, whereas YTHDF1 and YTHDF2 had weak nuclear and intense cytoplasmic staining. The fluorescence signal of YTHDF3 was intense in the nuclear and cytoplasm while the staining intensity of YTHDC1 was weak in the cytoplasm and nucleus (Figure 2B). As m6A methylation erasers, FTO and ALKBH5 were moderately expressed in the cytoplasm and nucleus (Figure 2C). The details of the subcellular localization of the m6A-related proteins in OS cells were presented in Table S3.
Figure 2. Subcellular location of m6A-related molecules in two osteosarcoma cell lines. (A) Immunofluorescence of m6A “writers” in U2-OS and KHOS-240S cell lines. (B) Immunofluorescence of m6A “readers” in U2-OS and KHOS-240S cell lines. (C) Immunofluorescence of m6A “erasers” in U2-OS and KHOS-240S cell lines.
Dysregulated m6A-Related Regulators Are Associated With Poor Prognosis in Osteosarcoma
To investigate the association between m6A-related protein expression and the clinical outcome of OS patients, a TMA cohort containing 120 OS tissues and 65 surrounding non-tumorous tissues was employed (Figures 3A, 4A). Differential expression analysis indicated that the protein expression levels of KIAA1429, RBM15, HNRNPC, HNRNPA2B1, YTHDF1, YTHDF2, and YTHDC1 were higher in tumor samples than in non-tumor samples (Figures 3B, 4B). Moreover, osteosarcoma patients were divided into two groups (high and low expression groups) based on each m6A-related protein IHC staining intensity and survival analysis was performed. The results revealed that high expression of KIAA1429 and HNRNPA2B1 was significantly associated with poor overall survival rates, while low expression of METTL3, METTL14, and YTHDF2 was related to poor prognosis in OS (Figures 3C, 4C).
Figure 3. The relationship between the expression level of m6A-related “writers” and poor prognosis in osteosarcoma. (A) Representative IHC high and low staining of m6A “writers” in OS tissues. (B) Comparison of m6A-related “writers” relative expression between OS tissues and nontumor tissues in TMA-combined cohort. (C) Kaplan–Meier analysis of the correlation between m6A “writers” expression and overall survival of OS patients in the TMA-combined cohort.
Figure 4. The relationship between the expression level of m6A related “reader” and “erasers” and poor prognosis in osteosarcoma. (A) Representative IHC high and low staining of m6A “reader” and “erasers” in OS tissues. (B) Comparison of m6A “reader” and “erasers” relative expression between OS tissues and non-tumor tissues in combined-TMA cohort. (C) Kaplan–Meier analysis of the correlation between m6A “reader” and “erasers” expression and overall survival of OS patients in the TMA-combined cohort.
Furthermore, an OS genome meta-cohort containing 175 OS tissues was constructed based on three independent datasets (GSE21257, GSE39055, and TARGET-OS, Figure S1). Overall survival analysis was performed, and the results indicated that the patients with high expression of METTL3, HNRNPA2B1, YTHDF3, and FTO and low expression of METTL14 had a significant shorter OS rate (Figures 5A–C). Furthermore, univariate analysis based on the TMA cohort suggested that the high expression of KIAA1429 and HNRNPA2B1 and the low expression of METTL3 and METTL14 were potential independent prognostic risk factors for OS patients. Moreover, the prognostic value of HNRNPA2B1 was validated in the genome meta-cohort, and KIAA1429, METTL3, and METTL14 did not exhibit significant prognostic potential (Figures 6A,B). In conclusion, dysregulated m6A-related regulators may predict poor prognosis in OS.
Figure 5. The correlation between m6A-related regulator expression and overall survival of osteosarcoma patients in the meta-cohort. (A) Kaplan–Meier analysis showing the correlation between m6A “writers” expression levels and the OS patients' survival rates. (B) Kaplan–Meier analysis showing the correlation between m6A “readers” expression levels and the OS patients' survival rates. (C) Kaplan–Meier analysis showing the correlation between m6A “erasers” expression levels and the OS patients' survival rates.
Figure 6. High expression level of HNRNPA2B1 is correlated with poor prognosis of osteosarcoma. (A) Univariate Cox regression analysis of OS patients' overall survival in meta-cohort. (B) Univariate Cox regression analysis of OS patient's overall survival in TMA combined OS cohort.
A Cluster Associated With a Favorable Prognosis Was Identified Through Consensus Clustering of m6A-Related Genes
To further explore the prognostic value of m6A-related gene expression patterns and to understand their underlying mechanisms in OS, we analyzed the relationships among 15 m6A-related regulators. Then, we identified a cluster associated with favorable outcomes and explored its possible mechanisms based on the genome meta-cohort.
Considering the similarity of the biological functions of m6A-related regulators, we analyzed the correlation (Figure 7A) and interaction (Figure 7B) among the 15 m6A-related genes. The results indicated that there were higher correlations among HNRNPA2B1, YTHDF2, RBM15, and ALKBH5. However, there was an inverse correlation between METTL16 and HNRNPA2B1, YTHDF2, and RBM15 expression in OS. Moreover, the protein–protein interaction networks showed that these m6A-related proteins interacted with each other frequently (Figure 7B). Additionally, the interaction between m6A-related proteins and 150 clinically actionable genes was examined to understand the clinical significance of m6A-related proteins. The results showed that m6A-related proteins frequently interacted with some of these genes (Figure 7B). Conclusively, the biological functions of 15 m6A-related genes were closely related in OS.
Figure 7. Interaction among m6A RNA–related regulators, differential survival status, and functional annotation of OS in three clusters. (A) Pearson correlation analysis of the 15 m6A regulators. (B) Protein–protein interactions among m6A regulators and clinically actionable genes obtained from the STRING database. (C) Kaplan–Meier overall survival curves for 175 osteosarcoma patients. (D) Heatmap and survival status of the three clusters defined by the m6A related regulators' consensus expression. (E) Three-cluster scatter plot of different survival time and status. (F,G) GSEA revealed that genes with higher expression in the cluster 1 subgroup were enriched for humoral immune response. (H,I) GSEA revealed the relationship between three clusters in osteosarcoma and cell cycle pathway. (J,K) GSEA revealed the relationship between three clusters in osteosarcoma and transition of the mitotic cell cycle pathway. (L,M) GSVA among three clusters using hallmark gene sets.
Furthermore, three subgroups were clustered by consensus clustering based on the expression of m6A-related genes in the genome meta-cohort. Patients were divided into three clusters with K = 3. Survival analysis showed that cluster 1 had a significantly longer overall survival than cluster 2 and cluster 3 (Figures 7C–E). Subsequently, GSEA revealed that the humoral immune response was significantly enriched in cluster 1 which exhibited good prognosis (Figures 7F,G), while the cell cycle and G2/M transition of mitotic cell cycle pathways were significantly enriched in cluster 2/3 which exhibited poor prognosis (Figures 7H–K). Similarly, the results of GSVA indicated that the dysregulated m6A-related genes might promote mitotic spindle and G2M checkpoint pathways, thus resulting in poor prognosis in OS (Figures 7L,M). Overall, these findings indicated that m6A modification may be involved in OS progression by regulating the humoral immune response and cell cycle–related pathways in OS.
Discussion
Many studies have demonstrated that dysregulated m6A modification is associated with human diseases such as obesity, type 2 diabetes mellitus, and infertility (27–29). In addition, the dysregulated expression of m6A-related genes is closely related to multiple tumors, including OS. To date, there have been only two reports addressing the functional role of m6A in osteosarcoma. Wang et al. demonstrated that METTL3 and ALKBH5 were closely associated with doxorubicin resistance in OS (30). Miao et al. reported that upregulated METTL3 promoted osteosarcoma progression by regulating LEF1 (31). However, the expression status and functional role of m6A-related genes in tumorigenesis and progression of OS have not yet been systematically analyzed.
In this study, we comprehensively analyzed the mRNA and protein expression levels of m6A-related regulators between OS and normal tissues based on two relatively large-scale OS cohorts (Figure 8). Both cohorts confirmed that the expression levels of RBM15, HNRNPC, HNRNPA2B1, YTHDF1, and YTHDC1 were increased in OS tissues. Consistent with our results, previous studies have demonstrated that RBM15 is upregulated in chronic myelogenous leukemia (32). HNRNPC and HNRNPA2B1 are upregulated in breast cancer (33, 34). YTHDF1 is upregulated in ovarian cancer, lung cancer, hepatocellular carcinoma, and colorectal cancer (35–38). YTHDC1 expression is increased in Kaposi's Sarcoma (39). In general, these findings suggest that m6A-related regulators are frequently dysregulated in OS. These dysregulated m6A-related regulators may be involved in the tumorigenesis and progression of OS.
Figure 8. The mechanism of m6A modification in osteosarcoma cells. m6A of RNA methylation is dynamically regulated by binding proteins “writers” (WTAP, KIAA1429, RBM15, RBM15B, METTL3, METTL14, and METTL16), “readers” (YTHDF1, YTHDF2, YTHFDF3, YTHDC1, and HNRNPs) and “erasers” (FTO and ALKBH5). “writers” add m6A modification on RNA, “readers” are proteins that bind to m6A modifications and exert various functions including pre-mRNA splicing, pri-miRNA processing, nuclear export, RNA translation modulation and mRNA decay. “erasers” could serve as demethylases.
Numerous studies have reported that dysregulated m6A-related molecule expression is closely associated with the outcome of tumor patients (23). Nevertheless, the impact of m6A methylation on OS prognosis is still unclear. For the first time, we systematically analyzed the prognostic value of m6A-related regulators in OS at both on the mRNA and protein levels. We found that high expression of HNRNPA2B1 and low expression of METTL14 were significantly associated with poor overall survival in the genome meta-cohort. Similarly, upregulated HNRNPA2B1 and downregulated METTL4 were positively related to low survival rates in the TMA cohort. Consistent with our results, HNRNPA2B1 has been reported to act as an oncogene in breast cancer, pancreatic cancer, lung cancer, hepatocellular carcinoma, and cervical cancer (34, 40–43), while METTL14 functions as a suppressor gene in glioma and breast cancer (44, 45). In contrast, the potential role of METTL14 as an oncogene has also been observed in acute myelocytic leukemia and hepatocellular carcinoma (44, 46). In summary, the expression pattern of m6A-related regulators may be considered a potential prognostic biomarker in OS.
M6A-related regulators can play a crucial role in cancer progression by targeting downstream molecules in an m6A-dependent manner. For instance, METTL3 promotes tumor progression by regulating LEF1 and SOCS2 (13, 31); ALKBH5 promotes tumor cell proliferation by increasing FOXM1 (47); YTHDF1 silences the drug-resistance gene AKR1C1 (36); and KIAAL429 targets and downregulates GATA3, which further contributes to liver cancer progression (48). To explore the underlying mechanism involved in m6A modification–mediated OS progression, we performed a systematic bioinformatic analysis based on an OS genome meta-cohort. Frequent cross talk among m6A-related regulators was observed in OS. Subsequently, we identified three clusters by consensus clustering based on the expression of m6A-related genes and found that cluster 1 had a significantly better prognosis than the others. Bioinformatic analysis revealed that cluster 1 was closely associated with the activated humoral immune response and suppressed cell cycle–related pathways. Aberrant regulation of the immune response and cell cycle–related pathways is considered a primary causative node in cancer progression (49, 50). Interestingly, recent literature has suggested that m6A modification can regulate the cell cycle and humoral immune response pathways (51, 52). It is therefore not surprising that the potential regulatory network between m6A modification and immune response or cell cycle–related pathways, which need further investigation.
However, our present study has some shortcomings. First, we only analyzed the relationship between dysregulated m6A-related regulators and clinical features in OS but did not verify these findings through in vivo and in vitro experiments. Second, we identified that cluster 1 was associated with a better prognosis. However, these findings need further validation in more independent cohorts. Third, we only explored the underlying mechanism based on bioinformatic prediction. M6A modification–mediated aberrant activation or suppression pathways in OS remains a subject for further study.
Conclusion
In summary, we have identified for the first time that m6A-related regulators are frequently dysregulated and closely correlated with poor survival in OS. M6A modification–mediated aberrant activation of cell-cycle related pathways and suppression of immune response may be responsible for the crucial role of m6A in OS progression. These findings may present a promising diagnostic biomarker and a potential target therapeutic strategy for OS patients.
Data Availability Statement
Publicly available datasets were analyzed in this study. This data can be found here: GEO http://www.ncbi.nlm.nih.gov/geo/, TARGET /https://ocg.cancer.gov/programs/target.
Ethics Statement
The studies involving human participants were reviewed and approved by the Ethics Committee of the First Affiliated Hospital of Zhengzhou University, Zhengzhou Central Affiliated Hospital of Zhengzhou University, and the Affiliated Cancer Hospital of Zhengzhou University. Written informed consent to participate in this study was provided by the participants' legal guardian/next of kin.
Author Contributions
JL and MH performed all the experimental work. JY, CL, and HY collected osteosarcoma specimens. QH and XL participated in data analysis. RS, XC, and GC conceived and participated in the design of the study. The manuscript was written by JL and BR. All authors read and approved the final manuscript.
Funding
This study was supported by the National Natural Science Foundation of China (81702346 and 81702757), National S&T Major Project of China (2018ZX10301201-008), National Key Research and Development Program of China (2018YFC2000500), and China Postdoctoral Science Foundation (2017M610463 and 2018M632814).
Conflict of Interest
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.
Acknowledgments
We thank the many clinical doctors from the Pathology Department and Precision Medicine Center, First Affiliated Hospital of Zhengzhou University, who were involved in this study.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2020.00769/full#supplementary-material
References
1. Mirabello L, Troisi RJ, Savage SA. Osteosarcoma incidence and survival rates from 1973 to 2004: data from the surveillance, epidemiology, and end results program. Cancer. (2009) 115:1531–43. doi: 10.1002/cncr.24121
2. Buddingh EP, Kuijjer ML, Duims RA, Bürger H, Agelopoulos K, Myklebost O, et al. 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. doi: 10.1158/1078-0432.CCR-10-2047
3. Siegel RL, Miller KD, Jemal A. Cancer statistics, 2018 CA: a cancer journal for clinicians. ACS J. (2018).68:7–30. doi: 10.3322/caac.21442
4. Wu Y, Zhou C, Yuan Q. Role of DNA and RNA N6-adenine methylation in regulating stem cell fate. Curr Stem Cell Res Ther. (2018) 13:31–8. doi: 10.2174/1574888X12666170621125457
5. Yang Y, Hsu PJ, Chen YS, Yang YG. Dynamic transcriptomic mA decoration: writers, erasers, readers and functions in RNA metabolism. Cell Res. (2018) 28:616–24. doi: 10.1038/s41422-018-0040-8
6. Heyn H, Esteller M. An adenine code for DNA: a second life for N6-methyladenine. Cell. (2015) 161:710–3. doi: 10.1016/j.cell.2015.04.021
7. Zhao BS, Roundtree IA, He C. Post-transcriptional gene regulation by mRNA modifications. Nat Rev Mol Cell Biol. (2017) 18:31–42. doi: 10.1038/nrm.2016.132
8. Lence T, A khtar J, Bayer M, Schmid K, Spindler L, Ho CH, et al. mA modulates neuronal functions and sex determination in Drosophila. Nature. (2016) 540:242–7. doi: 10.1038/nature20568
9. Du Y, Hou G, Zhang H, Dou J, He J, Guo Y, et al. SUMOylation of the m6A-RNA methyltransferase METTL3 modulates its function. Nucleic Acids Res. (2018) 46:5195–208. doi: 10.1093/nar/gky156
10. Visvanathan A, Patil V, Arora A, Hegde AS, Arivazhagan A, Santosh V, et al. Essential role of METTL3-mediated mA modification in glioma stem-like cells maintenance and radioresistance. Oncogene. (2018) 37:522–33. doi: 10.1038/onc.2017.351
11. Sadikovic B, Yoshimoto M, Chilton-MacNeill S, Thorner P, Squire JA, Zielenska M. Identification of interactive networks of gene expression associated with osteosarcoma oncogenesis by integrated molecular profiling. Hum Mol Genet. (2009) 18:1962–75. doi: 10.1093/hmg/ddp117
12. Liu L, Liu X, Dong Z, Li J, Yu Y, Chen X, et al. N6-methyladenosine-related genomic targets are altered in breast cancer tissue and associated with poor survival. J Cancer. (2019) 10:5447–59. doi: 10.7150/jca.35053
13. Chen M, Wei L, Law CT, Tsang FH, Shen J, Cheng CL, et al. RNA N6-methyladenosine methyltransferase-like 3 promotes liver cancer progression through YTHDF2-dependent posttranscriptional silencing of SOCS2. Hepatology. (2018) 67:2254–70. doi: 10.1002/hep.29683
14. Cui X, Wang Z, Li J, Zhu J, Ren Z, Zhang D, et al. Cross talk between RNA N6-methyladenosine methyltransferase-like 3 and miR-186 regulates hepatoblastoma progression through Wnt/β-catenin signalling pathway. Cell Prolif. (2020) 53:e12768. doi: 10.1111/cpr.12768
15. Liu X, Liu L, Dong Z, Li J, Yu Y, Chen X, et al. Expression patterns and prognostic value of m6A-related genes in colorectal cancer. Am J Transl Res. (2019) 11:3972–91.
16. Kelly AD, Haibe-Kains B, Janeway KA, Hill KE, Howe E, Goldsmith J, et al. MicroRNA paraffin-based studies in osteosarcoma reveal reproducible independent prognostic profiles at 14q32. Genome Med. (2013) 5:2. doi: 10.1186/gm406
17. Kuijjer ML, van den Akker BE, Hilhorst R, Mommersteeg M, Buddingh EP, Serra M, et al. Kinome and mRNA expression profiling of high-grade osteosarcoma cell lines implies Akt signaling as possible target for therapy. BMC Med Genomics. (2014) 7:4. doi: 10.1186/1755-8794-7-4
18. Li J, Liu L, Liu X, Xu P, Hu Q, Yu Y. DDX11 the role of upregulated as a potential prognostic and diagnostic biomarker in lung adenocarcinoma. J Cancer. (2019) 10:4208–16. doi: 10.7150/jca.33457
19. He Y, Chen X, Yu Y, Li J, Hu Q, Xue C, et al. LDHA is a direct target of miR-30d-5p and contributes to aggressive progression of gallbladder carcinoma. Mol Carcinog. (2018) 57:772–83. doi: 10.1002/mc.22799
20. Chen J, Yu Y, Li H, Hu Q, Chen X, He Y, et al. Long non-coding RNA PVT1 promotes tumor progression by regulating the miR-143/HK2 axis in gallbladder cancer. Mol Cancer. (2019) 18:33. doi: 10.1186/s12943-019-0947-9
21. De Winter J. Using the student's t-test with extremely small sample sizes. Prac Assess Res Eval. (2013) 18:12. doi: 10.7275/e4r6-dj05
22. Wilkerson MD, Hayes DN. ConsensusClusterPlus: a class discovery tool with confidence assessments and item tracking. Bioinformatics. (2010) 26:1572–3. doi: 10.1093/bioinformatics/btq170
23. Li Y, Xiao J, Bai J, Tian Y, Qu Y, Chen X, et al. Molecular characterization and clinical relevance of mA regulators across 33 cancer types. Mol Cancer. (2019) 18:137. doi: 10.1186/s12943-019-1066-3
24. Szklarczyk D, Franceschini A, Wyder S, Forslund K, Heller D, Huerta-Cepas J, et al. STRING v10: protein-protein interaction networks, integrated over the tree of life. Nucleic Acids Res. (2015) 43:D447–52. doi: 10.1093/nar/gku1003
25. Hänzelmann S, Castelo R, Guinney J. GSVA: gene set variation analysis for microarray and RNA-seq data. BMC Bioinform. (2013) 14:7. doi: 10.1186/1471-2105-14-7
26. Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. (2012) 16:284–7. doi: 10.1089/omi.2011.0118
27. Deng X, Su R, Stanford S, Chen J. Critical enzymatic functions of FTO in obesity and cancer. Front Endocrinol. (2018) 9:396. doi: 10.3389/fendo.2018.00396
28. Ding C, Zou Q, Ding J, Ling M, Wang W, Li H, et al. Increased N6-methyladenosine causes infertility is associated with FTO expression. J Cell Physiol. (2018) 233:7055–66. doi: 10.1002/jcp.26507
29. Klungland A, Dahl JA. Dynamic RNA modifications in disease. Curr Opin Genet Dev. (2014) 26:47–52. doi: 10.1016/j.gde.2014.05.006
30. Wang Y, Zeng L, Liang C, Zan R, Ji W, Zhang Z, et al. Integrated analysis of transcriptome-wide mA methylome of osteosarcoma stem cells enriched by chemotherapy. Epigenomics. (2019) 11:1693–715. doi: 10.2217/epi-2019-0262
31. 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:719–25. doi: 10.1016/j.bbrc.2019.06.128
32. Yang Y, Wang S, Zhang Y, Zhu X. Biological effects of decreasing RBM15 on chronic myelogenous leukemia cells. Leuk Lymphoma. (2012) 53:2237–44. doi: 10.3109/10428194.2012.684350
33. Wu Y, Zhao W, Liu Y, Tan X, Li X, Zou Q, et al. Function of HNRNPC in breast cancer cells by controlling the dsRNA-induced interferon response. EMBO J. (2018) 37:e99017. doi: 10.15252/embj.201899017
34. Hu Y, Sun Z, Deng J, Hu B, Yan W, Wei H, et al. Splicing factor hnRNPA2B1 contributes to tumorigenic potential of breast cancer cells through STAT3 and ERK1/2 signaling pathway. Tumour Biol. (2017) 39:1010428317694318. doi: 10.1177/1010428317694318
35. Liu T, Wei Q, Jin J, Luo Q, Liu Y, Yang Y, et al. The m6A reader YTHDF1 promotes ovarian cancer progression via augmenting EIF3C translation. Nucleic Acids Res. (2020) 48:3816–31. doi: 10.1093/nar/gkaa048
36. Shi Y, Fan S, Wu M, Zuo Z, Li X, Jiang L, et al. YTHDF1 links hypoxia adaptation and non-small cell lung cancer progression. Nat Commun. (2019) 10:4892. doi: 10.1038/s41467-019-12801-6
37. Zhao X, Chen Y, Mao Q, Jiang X, Jiang W, Chen J, et al. Overexpression of YTHDF1 is associated with poor prognosis in patients with hepatocellular carcinoma. Cancer Biomark. (2018) 21:859–68. doi: 10.3233/CBM-170791
38. Bai Y, Yang C, Wu R, Huang L, Song S, Li W, et al. YTHDF1 regulates tumorigenicity and cancer stem cell-like activity in human colorectal carcinoma. Front Oncol. (2019) 9:332. doi: 10.3389/fonc.2019.00332
39. Ye F, Chen ER, Nilsen TW. Kaposi's sarcoma-associated herpesvirus utilizes and manipulates RNA N-adenosine methylation to promote lytic replication. J Virol. (2017) 91:e00466-17. doi: 10.1128/JVI.00466-17
40. Dai S, Zhang J, Huang S, Lou B, Fang B, Ye T, et al. HNRNPA2B1 regulates the epithelial-mesenchymal transition in pancreatic cancer cells through the ERK/snail signalling pathway. Cancer Cell Int. (2017) 17:12. doi: 10.1186/s12935-016-0368-4
41. Dowling P, Pollard D, Larkin A, Henry M, Meleady P, Gately K, et al. Abnormal levels of heterogeneous nuclear ribonucleoprotein A2B1 (hnRNPA2B1) in tumour tissue and blood samples from patients diagnosed with lung cancer. Mol Biosyst. (2015) 11:743–52. doi: 10.1039/C4MB00384E
42. Wang H, Liang L, Dong Q, Huan L, He J, Li B, et al. Long noncoding RNA miR503HG, a prognostic indicator, inhibits tumor metastasis by regulating the HNRNPA2B1/NF-κB pathway in hepatocellular carcinoma. Theranostics. (2018) 8:2814–29. doi: 10.7150/thno.23012
43. Shi X, Ran L, Liu Y, Zhong SH, Zhou PP, Liao MX, et al. Knockdown of hnRNP A2/B1 inhibits cell proliferation, invasion and cell cycle triggering apoptosis in cervical cancer via PI3K/AKT signaling pathway. Oncol Rep. (2018) 39:939–50. doi: 10.3892/or.2018.6195
44. Wang S, Sun C, Li J, Zhang E, Ma Z, Xu W, et al. Roles of RNA methylation by means of N-methyladenosine (mA) in human cancers. Cancer Lett. (2017) 408:112–20. doi: 10.1016/j.canlet.2017.08.030
45. Wu L, Wu D, Ning J, Liu W, Zhang D. Changes of N6-methyladenosine modulators promote breast cancer progression. BMC Cancer. (2019) 19:326. doi: 10.1186/s12885-019-5538-z
46. Weng H, Huang H, Wu H, Qin X, Zhao BS, Dong L, et al. METTL14 inhibits hematopoietic stem/progenitor differentiation and promotes leukemogenesis via mRNA mA modification. Cell Stem Cell. (2018) 22:191–205.e9. doi: 10.1016/j.stem.2017.11.016
47. Zhang S, Zhao BS, Zhou A, Lin K, Zheng S, Lu Z, et al. mA demethylase ALKBH5 maintains tumorigenicity of glioblastoma stem-like cells by sustaining FOXM1 expression and cell proliferation program. Cancer Cell. (2017) 31:591–606.e6. doi: 10.1016/j.ccell.2017.02.013
48. Lan T, Li H, Zhang D, Xu L, Liu H, Hao X, et al. KIAA1429 contributes to liver cancer progression through N6-methyladenosine-dependent post-transcriptional modification of GATA3. Mol Cancer. (2019) 18:186. doi: 10.1186/s12943-019-1106-z
49. Herbst RS, Soria JC, Kowanetz M, Fine GD, Hamid O, Gordon MS, et al. Predictive correlates of response to the anti-PD-L1 antibody MPDL3280A in cancer patients. Nature. (2014) 515:563–7. doi: 10.1038/nature14011
50. Otto T, Sicinski P. Cell cycle proteins as promising targets in cancer therapy. Nat Rev Cancer. (2017) 17:93–115. doi: 10.1038/nrc.2016.138
51. Kansara M, Teng MW, Smyth MJ, Thomas DM. Translational biology of osteosarcoma. Nat Rev Cancer. (2014) 14:722–35. doi: 10.1038/nrc3838
Keywords: N6-methyladenosine, osteosarcoma, tumor metastasis, prognosis, biomarkers
Citation: Li J, Rao B, Yang J, Liu L, Huang M, Liu X, Cui G, Li C, Han Q, Yang H, Cui X and Sun R (2020) Dysregulated m6A-Related Regulators Are Associated With Tumor Metastasis and Poor Prognosis in Osteosarcoma. Front. Oncol. 10:769. doi: 10.3389/fonc.2020.00769
Received: 27 December 2019; Accepted: 21 April 2020;
Published: 02 June 2020.
Edited by:
Rimas J. Orentas, Seattle Children's Research Institute, United StatesReviewed by:
Songmei Liu, Wuhan University, ChinaKarina Formoso, Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Argentina
Copyright © 2020 Li, Rao, Yang, Liu, Huang, Liu, Cui, Li, Han, Yang, Cui and Sun. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xichun Cui, xichun639@126.com; Ranran Sun, fccsunrr@zzu.edu.cn
†These authors have contributed equally to this work