Original Research ARTICLE
Effect of m6A RNA Methylation Regulators on Malignant Progression and Prognosis in Renal Clear Cell Carcinoma
- 1Department of Urology, The First Affiliated Hospital of Chongqing Medical University, Chongqing, China
- 2Department of Oral Maxillofacial-Head and Neck Oncology, Shanghai Ninth People's Hospital, College of Stomatology, Shanghai Jiao Tong University School of Medicine, Shanghai, China
- 3Department of Head and Neck Cancer Center, Chongqing University Cancer Hospital & Chongqing Cancer Institute & Chongqing Cancer Hospital, Chongqing, China
Objectives: This study aims to explore the roles of 13 m6A RNA methylation regulators in clear cell renal cell carcinoma (ccRCC), and identify a risk signature and prognostic values of m6A RNA methylation regulators in ccRCC.
Materials and Methods: RNA sequence data of ccRCC was obtained from The Cancer Genome Atlas (TCGA) database. Differentially expressed of 13 m6A RNA methylation regulators in ccRCC stratified by different clinicopathological characteristics were unveiled using “limma” package in R version 3.6.0. Cox regression and LASSO analyses were conducted to identify the powerful independent prognostic factors in ccRCC associated with overall survival (OS). Protein-protein interaction (PPI) network and correlation analyses of the 13 m6A RNA methylation regulators were performed using “STRING” and R package, respectively. Principal component analysis (PCA) was also done using R. In addition, gene ontology (GO), GSEA and Kyoto Encyclopedia of Genes and Genomes pathways were used to functionally annotate the differentially expressed genes in different subgroups.
Results: Most of the 13 m6A RNA methylation regulators are differentially expressed in ccRCC tissue samples stratified by different clinicopathological characteristics in 537 patients. Next, a risk signature for predicting prognosis of ccRCC patients, was established based on two powerful independent prognostic m6A RNA methylation regulators (METTL14 and METTL3). Then, two subgroups (cluster1 and 2) were identified by consensus clustering to the two powerful independent factors and the cluster1 had a poorer prognosis than cluster2. Furthermore, the genes in cluster1 were significantly enriched in cancer-related pathways, biological process, and hallmarks, including “cell adhesion molecules (CAMs),” “leukocyte migration,” “Wnt/β-catenin signaling,” and so on.
Conclusion: M6A RNA methylation regulators play important roles in the initiation and progression of ccRCC and provide a novel sight to understand m6A RNA modification in ccRCC.
Renal cell carcinoma (RCC) is one of the most common adult genitourinary cancer, with ~73,820 newly diagnosed cases and nearly 14,770 of RCC-related mortality in the United States (1). Clear cell renal cell carcinoma (ccRCC) is the most common histological subtype of RCC, account for 70% (2). Although the diagnostic techniques of renal cell carcinoma are improved in the past decades, ~30% ccRCC patients have already developed metastases or local progression at the diagnosis (3). As ccRCC resistance to radiotherapy and chemotherapy, surgery is still an effective method for the treatment of ccRCC at present (4, 5). However, around one third of patients with localized RCC experience local recurrence and metastasis after surgical treatment (6). Complicated biological processes and unintelligible molecular mechanism lead to unfavorable prognosis of ccRCC patients. Hence, there is an urgent need to explore the molecular mechanism of ccRCC and identify novel targets for treatment and intervention for ccRCC patients.
Currently, the functions of RNAs in various cellular processer attract more and more attentions and accumulating studies are emerged in this field during the past decade. Up till now, many chemical modifications have been identified in distinct types of RNAs, including tRNA, miRNA, mRNA, long non-coding RNA, and others (7–12). These RNA modifications are reported to be several forms, such as N7-methyladenosine, 5-methylcytosine, N6-methyladenosine (m6A), and 2′-O-methylaion (7, 11, 13). M6A modification has been proved to be the most widespread, abundant and conserved form of mRNA methylation in eukaryotes (11, 14, 15). Generally speaking, m6A is enriched in 3′untranslated terminal region (UTR), near long internal exons and stop codons (16), thereby leading to alterations of RNA transcription, translation, metabolism, and processing (7, 17–20). As we known that RNA modification is mediated by a methyltransferase complex associated with three homologous factors named “readers” (YTHDC1, YTHDC2, YTHDF1, YTHDF2, and HNRNPC), “writers” (METTL3, METTL14, WTAP, KIAA1429, RBM15, and ZC3H13) and “erasers” (FTO and ALKBH5) (8, 21–26). M6A RNA methylation exhibits a reversible and dynamic biological process through the regulation of writers, readers and erasers. The identification of the m6A RNA methylation regulators bring us a new insight for the role of m6A modification in the regulation of gene expression (12, 27, 28).
Accumulating evidences have suggested the fact that the dysregulated expression of m6A RNA methylation regulators are involved in the initiation and development of human cancers. It was reported that METTL3, induced by hepatitis B X-interacting protein, could facilitate the proliferation of breast cancer through inhibiting the expression of let-7g (29). High expression of YTHDF1 was found to be associated with advanced stages and unfavorable prognosis in hepatocellular carcinoma (30). The key role of m6A RNA methylation regulators are also noted in different cancers, such as lung cancer (31), bladder cancer (32), and nasopharyngeal carcinoma (33). Recently, a research focused on the mutations and copy number variants of 10 m6A regulatory genes in ccRCC and found copy number variants of these regulatory genes closely correlated with pathologic stage and prognosis of patients with ccRCC (34). However, there is still a lack of comprehensive analyses of m6A RNA methylation modification in ccRCC with various clinicopathological features.
In the present study, we assessed the correlation of the expression levels of 13 m6A RNA regulators and clinicopathological characteristic in ccRCC based on the data obtained from The Cancer Genome Atlas (TCGA) (n = 537). We discovered that m6A RNA methylation regulators served as important roles in the initiation and progression of ccRCC, and according to two powerful independent prognostic m6A methylation regulators, a risk signature was established to classify the prognosis of ccRCC.
Materials and Methods
RNA sequence data from 537 patients with ccRCC were retrieved from The Cancer Genome Atlas (TCGA) database (https://cancergenome.nih.gov/) in 2019 and the values of these RNA sequence data were normalized by expectation-maximization. Patient clinical information was obtained using the Data Transfer Tool (provided by GDC Apps) (TCGA sample IDs and RNA-Seq information were presented in Supplementary Table 1). Clinicopathological information for the 537 ccRCC patients was summarized in Table 1. This study met the publication guidelines stated by TCGA (https://cancergenome.nih.Gov/publications/publicationguidelines). All data used in the study was obtained from TCGA, and hence ethics approval and informed consent were not required.
Selection of m6A RNA Methylation Regulators
According to study published by Yang Y et al. (8), 13 m6A RNA methylation regulators, including ALKBH5, FTO, YTHDC1, YTHDC2, YTHDF1, YTHDF2, HNRNPC, METTL3, METTL14, WTAP, KIAA1429, RBM15, and ZC3H13, were used for our analysis. Then, the correlation between the expression of these m6A RNA methylation regulators and different clinicopathological characteristics were evaluated based on the data from TCGA.
In order to explore the roles of the m6A RNA methylation regulators in ccRCC, our study was designed and analyzed according to the flow chart (Figure 1). The relationship of expressions of these 13 regulators and distinct clinicopathological characteristics in ccRCC was analyzed using the “limma” package (http://www.bioconductor.org/packages/release/bioc/html/limma.html) with a cut-off criteria of p < 0.05. Then, the differentially expressed m6A methylation regulators between tumor tissues and normal tissues were verified by two gene expression profiles (GSE14994 and GSE15641), which were downloaded from Gene Expression Omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/geo/). The GSE14994 datasets included 59 ccRCC samples and 11 normal samples, while GSE5641 datasets contained 32 ccRCC samples and 23 normal samples. The differentially expressed m6A methylation regulators between tumor tissues and normal tissues was also analyzed using the “limma” package (http://www.bioconductor.org/packages/release/bioc/html/limma.html) with a cut-off criteria of p < 0.05. Next, a PPI network of the 13 m6A RNA methylation regulators was constructed by using the Search Tool for the Retrieval of Interacting Genes (STRING, http://string.embl.de/). The combined score higher than 0.70 was regarded statistical significance. The correlation analysis was also performed by R package. Then, the prognostic m6A RNA methylation regulators were identified using univariate Cox regression analysis. Based on the results of univariate analysis, seven regulators highly correlated with overall survival (OS) (p < 0.05) were selected for the LASSO Cox regression analysis (35). Then, 2 m6A regulators were identified as the powerful independent prognostic factors by LASSO analysis. In addition, the prognostic values of these two regulators were verified by the Kaplan Meier plotter (www.kmplot.com), an online tool based on Gene Expression Omnibus database (GEO), European Genome-phenome Archive (EGA), and TCGA (36). The hazard ratio (HR) with 95% confidence intervals and log-rank P-value were calculated. Log-rank p > 0.05 was considered statistical significance. HR > 1, gene expression was negatively associated with OS., while HR < 1, gene expression was positively associated with OS. At last, two regulators' coefficients were identified according to the best penalty parameter λ. The risk score (RS) was estimated using the following formula:
where n represents the number of modules RNAs; Coef (i) d is the coefficient; X(i) denotes the z-score-transformed relative expression level for each gene identified by LASSO analysis. When the RS for a given sample was less than the mean RS of all samples, the latter was considered the low risk sample, otherwise, it was considered a high-risk sample. Survival curves in the high-risk and low-risk groups were estimated using the Kaplan-Meier method. Additionally, the receiver operating characteristic (ROC) curves and area under the ROC curves (AUC values) were used to access sensitivity and specificity.
Furthermore, two different subgroups (cluster1 and 2) were identified using “Consensus ClusterPlus” package (https://www.bioconductor.org/packages/release/bioc/html/ConsensusClusterPlus.html) with resample rate of 80%, 50 iterations and Pearson correlation. Principal component analysis (PCA) was performed by the R package for R version 3.6.0 to assess the gene expression patterns in the two ccRCC subgroups. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEEG) pathway analyses were done with the R packages in R version 3.6.0 to functionally annotate the differentially expressed genes in different subgroups. Gene Set Enrichment Analysis (GSEA) was also conducted to study the functions associated with different subgroups of ccRCC. Hallmark gene set “h.all.v6.0.symbol.gmt” was applied in GSEA analysis.
The t-test was performed to investigate the distribution of risk score in patients grouped by grade or classification. Univariate and multivariate Cox regression analyses were applied to identified the prognostic factors and different clinicopathological characteristics. Survival curves was plotted by using the “survival” package in R. Long-rank test was used to assess statistical significance. All statistical results with p < 0.05 were regard to be statistically significant.
Expression of m6A RNA Methylation Regulators Is Associated With Clinicopathological Characteristics in ccRCC
To better understand the important roles of m6A RNA methylation regulators in oncogenesis and progression, we explored expression levels of m6A RNA methylation regulators in different tissue samples, including tumor status (normal and tumor), WHO grade [low grade (including Grade I and II) and high grade (including Grade III and IV)], as well as pathological stage (early stage including Stage I and II, later stage including Stage III and IV) tissue samples. The results are shown as heatmaps in Figures 2A–C, suggesting that the expressions of most m6A methylation regulators are significantly related to tumor status. Ten m6A RNA methylation regulators were significantly abnormally expressed in ccRCC tissues samples. Compared to normal tissue samples, ZC3H13, METTL14, and YTHDF2 were down regulated, while FTO, ALKBH5, WTAP, METTL3, YTHDC2, KIAA1429, and RBM15 were up regulated in ccRCC tissue samples (Figure 2D). These differentially expressed regulators were verified by GSE15641 and GSE15641 datasets. The results of validation suggested that FTO, WTAP, RBM15, and ZC3H13 were significantly abnormally expressed in ccRCC tissues (Figures 2G,H), which was consistent with our previous result. However, ZC3H13 was down regulated in ccRCC samples in GSE14994 datasets while the opposite result was observed in ccRCC samples in both GSE15641 datasets and TCGA. This may be due to the different study populations and different scales.
Figure 2. Expression of m6A RNA methylation regulators in ccRCC with different clinicopathological characteristics. (A–C) The heatmaps of 13 m6A RNA methylation regulators in different tissue samples (A for tumor status, B for WHO grade, and C for pathological stage). (D–F) The expression of 13 m6A RNA methylation regulators in different tissue samples (D for tumor status, E for WHO grade, and F for pathological stage). (G,H) The expression of m6A RNA methylation regulators in tumor tissues and normal tissues (G for GSE15641 and H for GSE14994). *p < 0.05, **p < 0.01, and ***p < 0.001; the red fusiformis represents tumor tissue and the blue fusiformis represents normal tissue.
Next, we investigated the relationship between expression of m6A RNA methylation regulators, WHO grade (Figure 2E) and pathological stage (Figure 2F) in ccRCC, respectively. The results demonstrated that low expression levels of FTO, METTLE14, YTHDC1, ZC3H13, and KIAA1429 were significantly associated with both high grade and later pathological stage. Besides that, HNRNPC, METTL3, YTHDF2, and RBM15 were also observed down regulation in high grade tissue samples.
Interaction and Correlation of 13 m6A RNA Methylation Regulators in ccRCC
The interrelationships between the 13 m6A RNA methylation genes were retrieved from STRING database to construct PPI network (Figure 3A) and their correlations were also analyzed using “corrplot” package in R software (Figure 3B). As a result, we found that there were close and complicated interrelationships between the six writers (Figure 3A). The expressions of the six writers were significantly associated with each other except METTL3 and ZC3H13, METTL3 and KIAA1429, ZC3H13 and WTAP in ccRCC (Figure 3B). Few interactions were observed between the five readers in PPI network (Figure 3A). However, the significant correlations between HNRNPC, YTHDC2, YTHDC1, YTHDF1, and YTHDF2 were observed in ccRCC (Figure 3B). There were evidences supporting the interaction between FTO and ALKBH5 in PPI network (Figure 3A), and the expressions of FTO and ALKBH5 were positively associated with each other in ccRCC (Figure 3B). In addition, the expressions of FTO, KIAA1429, ZC3H13, METTL14, and YTHDC1 were highly related to each other (Figure 3B), which agreed with that the expressions of these genes mentioned above being negatively associated with the increasing malignancy of ccRCC.
Figure 3. Interaction among m6A RNA methylation regulators. (A) The PPI network of the 13 m6A methylation regulators constructed using STRING. (B) Spearman correlation analysis of the 13 m6A methylation regulators.
Prognostic Value of m6A RNA Methylation Regulators and a Risk Signature Established Based on Two Identified m6A RNA Methylation Regulators
In order to investigate the prognostic value of these 13 m6A RNA methylation regulators in ccRCC, univariate Cox regression analysis was performed based on the expression levels of the regulators from TCGA (Figure 4A). As a result, we found that seven out of the 13 regulators were significantly associated with overall survival (OS) (p < 0.05). Among the seven regulators, only METTL3 was considered as risky gene with HR> 1, while FTO, METTL14, YTHDC1, YTHDC2, ZC3H13, and KIAA1429 were considered as protective genes with HR <1.
Figure 4. Risk signature with 2 m6A RNA methylation regulators. (A–C) The process of constructing the signature based on 2 m6A RNA methylation regulators. The Hazard ratio (HR), 95% confidence interval (CI) estimated by univariate Cox regression (A) and the coefficients estimated by multivariate Cox regression via LASSO are presented (B,C). (D,E) The overall survival curves of the two regulators estimated by the Kaplan Meier plotter. Log-rank p < 0.05 was considered statistical significance. HR > 1, gene expression was negatively associated with OS., while HR < 1, gene expression was positively associated with OS. (F) Kaplan-Meier overall survival (OS) curves for patients in high- and low-risk group divided according to the risk score. (G) ROC analysis and AUC value of the ROC curve suggesting the sensitivity and specificity for risk signature.
To identified the most powerful prognostic m6A RNA methylation regulators, the last absolute shrinkage and selection operator (LASSO) Cox regression analysis to the seven prognosis-related genes was conducted (Figure 4B) and the coefficient of each independent prognostic gene was shown in Figure 4C. The LASSO results showed that two regulators (METTL14 and METTL3) were the powerful prognostic factors. Then, to verify the two regulators, their prognostic values were verified by the Kaplan Meier plotter. The results showed that ccRCC patients with high METTL14 expression had favorable prognosis (Figure 4D), while ccRCC patients with high METTL3 expression had bad prognosis (Figure 4E), which was in accordance with the LASSSO results and this partly strength the reliability of our findings.
Based on the powerful prognostic factors (METTL14 and METTL3), a risk signature was constructed. Then, the risk score was estimated based on the coefficients obtained from the LASSO analysis. In order to test the prognostic role of the two-gene risk signature, ccRCC patients from TCGA (n = 525) were assigned into high-risk and low-risk groups according to the median risk score (Table 1). The result of survival analysis showed that the high-risk group had significantly shorter survival time compared to low-risk group (Figure 4F). The 5-years OS was 44.3% in high-risk group and 75.2% in low-risk group. Time-dependent ROC curve was used to assess the sensitivity and specificity of the prediction and the result showed that AUC values was 0.704 (Figure 4G), suggesting well-prediction performances.
Prognostic Risk Score Indicated Strong Associations With Clinical Characteristics in ccRCC
The expression levels of the two identified m6A RNA methylation regulators in high-risk and low-risk group were presented in the heatmap (Figure 5A). The results showed that there were significant differences between the high-risk and low-risk groups in term of grade (p < 0.001), pathological stage (p < 0.001), stage T (p < 0.001), and stage M (p < 0.01). The relationship between the risk score and each clinicopathological characteristics was also explored in the present study. The results showed that the risk score has significant differences between patients divided by WHO grade and pathological stage (Figures 5B,C). Furthermore, the ROC curve showed that the risk score (AUC = 0.708) (Figure 5D) was better than grade (AUC = 0.655) in predicting 5-years survival rate (Figure 5E) and was similar to pathological stage (AUC = 0.709) in prediction for 5-years survival rate (Figure 5F).
Figure 5. Relation between the risk score and clinicopathological characteristics, as well as prognostic value of the risk signature in patients stratified by WHO grade and pathological stage. (A) The heatmap shows the expression of the 2 m6A RNA methylation regulators in high-risk and low-risk ccRCC. The distribution of clinicopathological characteristics was compared between the high-risk and low-risk groups. *p < 0.05, **p < 0.01 and ***p < 0.001. (B,C) The distribution of risk score in patients stratified by WHO grade and pathological stage. *p < 0.05, **p < 0.01, and ***p < 0.001. (D–F) ROC curves showed the predictive efficiency of the risk score (D), WHO grade (E), and pathological stage (F). (G) Univariate Cox regression analysis of the associated between clinicopathological factors (including risk score) and overall survival of patients. (H) Multivariate Cox regression analysis of the associated between clinicopathological factors (including risk score) and overall survival of patients. (I) Kaplan-Meier overall survival curve for patients with low grade and high grade. (J) Kaplan-Meier overall survival curve for patients with early grade and later pathological stage.
To test whether the risk signature was an independent prognostic factor, univariate and multivariate Cox regression analyses were performed. As a result, the age at diagnosis, grade, pathological stage, stage TNM and risk score were associated with OS in univariate analysis (Figure 5G) and only risk score and age at diagnosis were still significantly related to OS (p < 0.05) in multivariate Cox regression analysis (Figure 5H). The prognostic value of the risk signature for different grades and pathological stage was also analyzed using Kaplan-Meier curve. The results showed that patients with low risk scores had a significantly longer survival time than those with high scores in high grade (Figure 5I) and later stage (Figure 5J). These results demonstrated that the risk score retrieved from m6A RNA methylation regulators could be served as an independent prognostic factor in ccRCC.
Consensus Clustering of Two Independent Prognostic m6A RNA Methylation Regulators Identified Two Clusters of ccRCC With Different Clinical Outcomes
According to the expression similarity of the 2 m6A RNA methylation regulators (METTL14 and METTL3) identified above as the powerful independent prognostic factors, k = 4 could be the optimal choice with clustering increasing from k = 2–9 (Figures 6A,B). However, we observed that only when k = 2, the interference between subgroups was minimal (Supplementary Image 1). Hence, k = 2 was used for consensus clustering analysis and two subgroups named cluster1 and cluster2 were identified. We found that patients in cluster1 had a significantly shorter overall survival (OS) than those in cluster2 (Figure 6C).
Figure 6. Overall survival of ccRCC patients in cluster 1 and 2 subgroups. (A) The tracking plot for k = 2 to k = 9. (B) Relative change in area under CDF curve for k = 2–9. (C) Kaplan-Meier overall survival (OS) curve for ccRCC patients in cluster 1 and 2. (D) Principal component analysis of the total RNA expression profile. ccRCC in cluster1 are marked in red. (E,F) Functional annotation of the genes with different expression in the cluster 1 using KEGG pathway (E) and GO terms (F).
Moreover, we performed principal component analysis (PCA) for comparison of the transcriptional profile between cluster1 and cluster2. The result suggested that there was a significant distinction between the two subgroups (Figure 6D). Then, the genes significantly upregulated (fold change> 2 and p < 0.05) or downregulated (fold change < −2 and p < 0.05) were selected for gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEEG) pathway analysis. The results of KEEG analysis showed that differentially expressed genes were mainly enriched in malignancy-associated pathways, including “cytokine-cytokine receptor interaction,” “cell adhesion molecules (CAMs),” “systemic adhesion lupus erythematosus,” “chemokine signaling pathway,” “viral protein interaction with cytokine and cytokine receptor,” “complement and coagulation cascades,” “th1 and th2 cell differentiation,” and so on (Figure 6E). GO analysis results were associated with cancer-related biological processes, including “leukocyte migration,” “regulation of leukocyte activation,” “regulation of immune effector process,” and so on (Figure 6F). Next, the gene set enrichment analysis (GSEA) was also performed and the result suggested that the malignant hallmarks of cancer, containing allograft rejection (NES = 1.94, normalized p = 0.008, FDR q = 0.070), interferon gamma response (NES = 1.81, normalized p = 0.038, FDR q = 0.134), IL6/STAT3 signaling (NES = 1.77, normalized p = 0.024, FDR q = 0.125), inflammatory response (NES = 1.74, normalized p = 0.032, FDR q = 0.117), Wnt/β-catenin signaling (NES = 1.68, normalized p = 0.038, FDR q = 0.113), and P53 pathway (NES = 1.57, normalized p = 0.041, FDR q = 0.177) had significant correlations with the cluster1 (Figure 7). These results above showed the two clusters identified based on the two powerful independent prognostic factors were closely associated with the malignancy of ccRCC.
Figure 7. GSEA showed that genes with higher expression in cluster1 were enriched for hallmarks of malignant tumors: (A) allograft rejection, (B) interferon gamma response, (C) IL6/STAT3 signaling, (D) inflammatory response, (E) Wnt/β-catenin signaling, and (F) P53 pathway.
CcRCC is one of the most popular common form of renal cancer worldwide. There are no effective therapeutic strategies for ccRCC patients with advanced stage or metastasis and the rate of 5-years disease-free survival in patients with metastasis is only 12%. M6A RNA methylation was demonstrated to be closely associated with tumorigenesis and overall survival (OS) in patients with renal cell carcinoma (37). In the present study, the expressions of m6A RNA methylation regulators are also highly correlated to the tumor progression and prognosis in ccRCC. Firstly, 10 of 13 m6A RNA methylation regulators were identified as playing important roles in the development of ccRCC, including ZC3H13, METTL14 and YTHDF2, FTO, ALKBH5, WTAP, METTL3, YTHDC2, KIAA1429, and RBM15. We also found that the expression of FTO, METTLE14, YTHDC1, ZC3H13, and KIAA1429 was negatively related to grade and stage classification. Next, a prognostic risk signature of two identified m6A RNA methylation regulators was retrieved based on the data from TCGA. Then, according to the prognostic risk signature, patients with ccRCC were assigned into high-risk and low-risk groups. Then, two ccRCC subgroups (cluster1 and cluster1) were also identified by consensus clustering according to the expressions of the 2 m6A RNA methylation regulators (METTL14 and METTL3), which were selected for construction of the risk signature. We found that the cluster1 and cluster2 not only affected the prognosis of ccRCC patients and clinicopathological characteristics but were also associated with the cancer-related pathways, pivotal biological processes and hallmarks in ccRCC.
Recently, more and more studies have suggested that the abnormal expression of the m6A RNA methylation regulators is involved in human cancer. The writer WTAP was reported to be a oncogene in different cancers, such as AML (38), glioblastoma (39), and pancreatic ductal adenocarcinoma (40). The eraser ALKBH5 was observed to be up regulated in glioblastoma stem-like cell (GSCs), leading to the initiation and development of glioblastoma (41). Interestingly, the reader YTHDF2 was suggested to facilitate the migration of prostate cancer in vitro (42) but inhibit invasion and migration in pancreatic cancer (43). These researches showed that alteration of the m6A RNA methylation regulators was observed in cancers and the same regulator could function as diverse roles in the initiation and progression of various cancers.
As shown in the section of results, the expression of the 13 m6A RNA methylation regulators in ccRCC was observed to associated with different clinicopathological characteristics. In term of m6A RNA methylation writers, the expression of KIAA1429, ZC3H13, and METTL14 were significantly decreased in ccRCC patients with high grade and later stage. Surprisingly, KIAA1429 was observed to be up regulated in ccRCC tissue samples compared with normal tissue samples. As an m6A RNA methylation eraser, the expression of FTO was found to be similar to that of KIAA1429, up regulated in ccRCC patients while down regulating in ccRCC patients with high grade and later stage. This may be because KIAA1429 and FTO could exert various functions at different stage of ccRCC tumorigenesis and development. KIAA1429 is also named as vir-like m6A methyltransferase associated or VIRMA and was reported to promote the proliferation, migration and invasion of HCC by regulating mRNA methylation levels of ID2 in cell lines (44). Another research also showed that KIAA1429 was significantly up regulated in seminomas, while down regulated in non-neminomatous tumors (45). FTO was reported to be significantly associated with bad prognosis in ccRCC patients and its prognosis value increased as the increase of stages (46). Low expression of ZC3H13 was associated with the progression of colorectal cancer by inactivating Ras-ERK signaling pathway (47). A recent study reported that METTL14 could mediate P2RX6 mRNA and protein level, promoting renal cancer cells migration and invasion via ATP-induced Ca2+ influx modulating ERK1/2 phosphorylation and MMP9 signal pathway in vitro and in vivo assays (48). METTL14 was also suggested to be an inhibitor of tumor metastasis and considered as a favorable factor in HCC via mediating m6A-dependent miRNA processing (49). However, another research demonstrated that METTL14 functioned as an oncogene to promote tumorigenesis by mRNA m6A modification in leukemia (50). FTO was demonstrated to promote the proliferation, transformation and survival of acute myeloid leukemia cells in vivo and in vitro (51). As one of the two powerful prognostic m6A RNA methylation regulators, METTL3 mRNA and protein was reported to be downregulated in RCC tissues by quantitative real-time PCR (qRT-PCR) and western blot, which was in accordance with our results. Then, the researchers also found that low level of METTL3 was associated with larger tumor size and higher histological grade in vivo study, and could promote RCC cell proliferation, migration and invasion function and induced G0/G1 arrest in vitro study. In addition, significant changes in epithelial-to-mesenchymal transition (EMT) and PI3K-Akt-mTOR pathways were noted in their research. This indicated that METTL3 could exert its function via EMT and PI3K-Akt-mTOR pathways, which is worth further investigating. Moreover, the researchers found that upregulated expression of METTL3 always associated with favorable prognosis in ccRCC patients (37). Zhou et al. (34) also suggested that low level of METTL3 was correlated with some cancer-related biological processes, including adipogenesis, mTOR pathway and reactive oxygen species. In summary, dysregulation of m6A RNA methylation regulators were highly correlated with tumor progression and prognosis in ccRCC. These findings could also provide novel methods for the treatment and prevention of ccRCC.
Next, GO and KEEG analysis associated with m6A RNA methylation regulators were also conducted in this study. Several biological processes and pathways correlated with the tumorigenesis and progression of ccRCC were identified including “cytokine-cytokine receptor interaction,” “cell adhesion molecules (CAMs),” “complement and coagulation cascades,” “leukocyte migration,” “regulation of leukocyte activation,” “regulation of immune effector process,” and so on. It was known that the expression or function of cytokines, including cellular self-renewal, senescence, migration, and apoptosis was often altered in tumor tissues compared with healthy tissues (52–54). CAMs was found to promote the metastasis of head and neck squamous cell carcinoma via detaching the tumor cells from primary tumor and decreasing the intercellular adhesion (55). More and more evidences suggested that complement component and complement activated product could facilitate the growth and angiogenesis of tumor, as well as immunosuppression (56–58). A recent research (59) reported that the abnormal downregulated genes were possibly associated with the gallbladder cancer progression through the complement and coagulation cascades. Leukocyte migration was suggested to be activated in the progression of cancer (60). A study was showed that immune system was closely associated with the progression of hepatocellular carcinoma (61). Here, we revealed that the m6A RNA methylation regulators were involved in many biological processes and signaling pathways, suggesting their important roles in initiation and development of ccRCC.
In this study, we found that the prognostic signature obtained using 2 m6A RNA methylation regulators (METTL14 and METTL3) had significant value in ccRCC. The ROC for risk signature identified above exhibited a satisfied prediction performance in ccRCC. Additionally, the prognostic value of signature for low grade and high grade, early stage and later stage ccRCC was also noted. High-risk Patients even with the same high grade or later stage had a significantly shorter survival time than low-risk patients. Risk signature in the present study may be benefit for physicians to more precisely estimated individualized survival prediction.
However, several limitations in this study should be acknowledged. First, the number of normal samples (72) is much less than that of tumor samples (539). It may affect the reliability of our results. Second, the present study is purely computational, future experimental and clinical data are need to validate our results. Third, some important clinical parameters of ccRCC patients, such as treatment strategy, vascular invasion and surgical margin, are not available from TCGA. Finally, patients in our study are mainly Americans and this may lead to the risk of potential selection bias.
In conclusion, our study revealed that the expression of the 13 m6A RNA methylation regulators are closely correlated with the malignant clinicopathological characteristics of ccRCC and are also highly related to the upregulated expression of genes enriched in the biological processes and pathways that facilitate the malignant development of ccRCC. Our findings may be regarded as important evidence for the role of m6A RNA methylation in ccRCC. Future experimental and clinical studies are necessary to produce a solid confirmation of our results.
Data Availability Statement
Publicly available datasets were analyzed in this study. This data can be found here: The Cancer Genome Atlas (TCGA) database.
This study met the publication guidelines stated by TCGA (https://cancergenome.nih.Gov/publications/publicationguidelines). All data used in the study was obtained from TCGA, and hence ethics approval and informed consent were not required.
Informed consent was obtained from all individual participants included in the study. And the present research meets the publication guidelines provided by TCGA (https://cancergenome.nih.Gov/publications/publicationguidelines).
JW, WH, and XG contributed conception and design of the study. JW and CZ organized the database. JW and CZ performed the statistical analysis. JW wrote the first draft of the manuscript. CZ wrote sections of the manuscript, contributed to manuscript revision, and proofreading. All authors contributed to manuscript revision, read, and approved the submitted version.
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.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2020.00003/full#supplementary-material
Supplementary Image 1. Consensus clustering matrix for k = 2 (A), k = 3 (B), and k = 4 (C).
Supplementary Table 1. TCGA sample IDs and RNA-Seq information.
2. Moch H, Cubilla AL, Humphrey PA, Reuter VE, Ulbright TM. The 2016 WHO classification of tumours of the urinary system and male genital organs-part A: renal, penile, and testicular tumours. Eur Urol. (2016) 70:93–105. doi: 10.1016/j.eururo.2016.02.029
4. Cavaliere C, D'Aniello C, Pepa CD, Pisconti S, Berretta M, Facchini G. Current and emerging treatments for metastatic renal cell carcinoma. Curr Cancer Drug Targets. (2018) 18:468–79. doi: 10.2174/1568009617666170209094030
5. Miao D, Margolis CA, Gao W, Voss MH, Li W, Martini DJ, et al. Genomic correlates of response to immune checkpoint therapies in clear cell renal cell carcinoma. Science. (2018) 359:801–6. doi: 10.1126/science.aan5951
6. Brierley JD, Gospodarowicz MK, Wittekind C editors. TNM Classification of Malignant Tumours, 8th edn. Wiley-Blackwell (2016). p. 187–210. Available online at: https://www.wiley.com/enes/TNM+Classification+of+Malignant+Tumours%2C+8th+Edition-p-9781119263579
10. Yang D, Qiao J, Wang G, Lan Y, Li G, Guo X, et al. N6-Methyladenosine modification of lincRNA 1281 is critically required for mESC differentiation potential. Nucleic Acids Res. (2018) 46:3906–20. doi: 10.1093/nar/gky130
15. Wang S, Sun C, Li J, Zhang E, Ma Z, Xu W, et al. Roles of RNA methylation by means of N(6)-methyladenosine (m(6)A) in human cancers. Cancer Lett. (2017) 408:112–20. doi: 10.1016/j.canlet.2017.08.030
19. Slobodin B, Han R, Calderone V, Vrielink J, Loayza-Puch F, Elkon R, et al. Transcription impacts the efficiency of mRNA translation via co-transcriptional N6-adenosine methylation. Cell. (2017) 169:326–37 e12. doi: 10.1016/j.cell.2017.03.031
20. Ke S, Pandya-Jones A, Saito Y, Fak JJ, Vagbo CB, Geula S, et al. m(6)A mRNA modifications are deposited in nascent pre-mRNA and are not required for splicing but do specify cytoplasmic turnover. Genes Dev. (2017) 31:990–1006. doi: 10.1101/gad.301036.117
21. Scholler E, Weichmann F, Treiber T, Ringle S, Treiber N, Flatley A, et al. Interactions, localization, and phosphorylation of the m(6)A generating METTL3-METTL14-WTAP complex. RNA. (2018) 24:499–512. doi: 10.1261/rna.064063.117
22. Tang C, Klukovich R, Peng H, Wang Z, Yu T, Zhang Y, et al. ALKBH5-dependent m6A demethylation controls splicing and stability of long 3'-UTR mRNAs in male germ cells. Proc Natl Acad Sci USA. (2018) 115:E325–33. doi: 10.1073/pnas.1717794115
23. Wojtas MN, Pandey RR, Mendel M, Homolka D, Sachidanandam R, Pillai RS. Regulation of m(6)A transcripts by the 3'–>5' RNA helicase YTHDC2 is essential for a successful meiotic program in the mammalian germline. Mol Cell. (2017) 68:374–87 e12. doi: 10.1016/j.molcel.2017.09.021
25. Wu R, Yao Y, Jiang Q, Cai M, Liu Q, Wang Y, et al. Epigallocatechin gallate targets FTO and inhibits adipogenesis in an mRNA m(6)A-YTHDF2-dependent manner. Int J Obes. (2018) 42:1378–88. doi: 10.1038/s41366-018-0082-5
29. Cai X, Wang X, Cao C, Gao Y, Zhang S, Yang Z, et al. HBXIP-elevated methyltransferase METTL3 promotes the progression of breast cancer via inhibiting tumor suppressor let-7g. Cancer Lett. (2018) 415:11–9. doi: 10.1016/j.canlet.2017.11.018
30. 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
31. Liu J, Eckert MA, Harada BT, Liu SM, Lu Z, Yu K, et al. m(6)A mRNA methylation regulates AKT activity to promote the proliferation and tumorigenicity of endometrial cancer. Nat Cell Biol. (2018) 20:1074–83. doi: 10.1038/s41556-018-0174-4
32. Cheng M, Sheng L, Gao Q, Xiong Q, Zhang H, Wu M, et al. The m(6)A methyltransferase METTL3 promotes bladder cancer progression via AFF4/NF-kappaB/MYC signaling network. Oncogene. (2019) 38:3667–80. doi: 10.1038/s41388-019-0683-z
34. Zhou J, Wang J, Hong B, Ma K, Xie H, Li L, et al. Gene signatures and prognostic values of m6A regulators in clear cell renal cell carcinoma - a retrospective study using TCGA database. Aging. (2019) 11:1633–47. doi: 10.18632/aging.101856
35. Ramsay IS, Ma S, Fisher M, Loewy RL, Ragland JD, Niendam T, et al. Model selection and prediction of outcomes in recent onset schizophrenia patients who undergo cognitive training. Schizoph Res Cogn. (2018) 11:1– 5. doi: 10.1016/j.scog.2017.10.001
36. Nagy A, Lanczky A, Menyhart O, Gyorffy B. Validation of miRNA prognostic power in hepatocellular carcinoma using expression data of independent datasets. Scient Rep. (2018) 8:9227. doi: 10.1038/s41598-018-29514-3
37. Li X, Tang J, Huang W, Wang F, Li P, Qin C, et al. The M6A methyltransferase METTL3: acting as a tumor suppressor in renal cell carcinoma. Oncotarget. (2017) 8:96103–16. doi: 10.18632/oncotarget.21726
40. Li BQ, Huang S, Shao QQ, Sun J, Zhou L, You L, et al. WT1-associated protein is a novel prognostic factor in pancreatic ductal adenocarcinoma. Oncol Lett. (2017) 13:2531–8. doi: 10.3892/ol.2017.5784
41. Zhang S, Zhao BS, Zhou A, Lin K, Zheng S, Lu Z, et al. m(6)A 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
42. Li J, Meng S, Xu M, Wang S, He L, Xu X, et al. Downregulation of N(6)-methyladenosine binding YTHDF2 protein mediated by miR-493-3p suppresses prostate cancer by elevating N(6)-methyladenosine levels. Oncotarget. (2018) 9:3752–64. doi: 10.18632/oncotarget.23365
43. Chen J, Sun Y, Xu X, Wang D, He J, Zhou H, et al. YTH domain family 2 orchestrates epithelial-mesenchymal transition/proliferation dichotomy in pancreatic cancer cells. Cell Cycle. (2017) 16:2259–71. doi: 10.1080/15384101.2017.1380125
44. Cheng X, Li M, Rao X, Zhang W, Li X, Wang L, et al. KIAA1429 regulates the migration and invasion of hepatocellular carcinoma by altering m6A modification of ID2 mRNOncoTargets A, and therapy. (2019) 12:3421–8. doi: 10.2147/OTT.S180954
45. Lobo J, Costa AL, Cantante M, Guimaraes R, Lopes P, Antunes L, et al. m(6)A RNA modification and its writer/reader VIRMA/YTHDF3 in testicular germ cell tumors: a role in seminoma phenotype maintenance. J Transl Med. (2019) 17:79. doi: 10.1186/s12967-019-1837-z
46. Wen L, Yu Y, Lv H, He Y, Yang B. FTO mRNA expression in the lower quartile is associated with bad prognosis in clear cell renal cell carcinoma based on TCGA data mining. Ann Diagn Pathol. (2019) 38:1–5. doi: 10.1016/j.anndiagpath.2018.10.009
47. Zhu D, Zhou J, Zhao J, Jiang G, Zhang X, Zhang Y, et al. ZC3H13 suppresses colorectal cancer proliferation and invasion via inactivating Ras-ERK signaling. J Cell Physiol. (2019) 234:8899–907. doi: 10.1002/jcp.27551
48. Gong D, Zhang J, Chen Y, Xu Y, Ma J, Hu G, et al. The m(6)A-suppressed P2RX6 activation promotes renal cancer cells migration and invasion through ATP-induced Ca(2+) influx modulating ERK1/2 phosphorylation and MMP9 signaling pathway. J Exp Clin Cancer Res. (2019) 38:233. doi: 10.1186/s13046-019-1223-y
49. Ma JZ, Yang F, Zhou CC, Liu F, Yuan JH, Wang F, et al. METTL14 suppresses the metastatic potential of hepatocellular carcinoma by modulating N(6) -methyladenosine-dependent primary MicroRNA processing. Hepatology. (2017) 65:529–43. doi: 10.1002/hep.28885
50. 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 m(6)A modification. Cell Stem Cell. (2018) 22:191–205 e9. doi: 10.1016/j.stem.2017.11.016
51. Li Z, Weng H, Su R, Weng X, Zuo Z, Li C, et al. FTO plays an oncogenic role in acute myeloid leukemia as a N(6)-methyladenosine RNA demethylase. Cancer Cell. (2017) 31:127–41. doi: 10.1016/j.ccell.2016.11.017
52. Martins-Filho A, Jammal MP, Micheli DC, Tavares-Murta BM, Etchebehere RM, Murta EFC, et al. Role of intracystic cytokines and nitric oxide in ovarian neoplasms. Scand J Immunol. (2017) 86:462–70. doi: 10.1111/sji.12617
55. Valentiner U, Knips J, Pries R, Clauditz T, Munscher A, Sauter G, et al. Selectin binding sites are involved in cell adhesive properties of head and neck squamous cell carcinoma. Cancers. (2019) 11:1672. doi: 10.3390/cancers11111672
59. Kong L, Wu Q, Zhao L, Ye J, Li N, Yang H. Identification of messenger and long noncoding RNAs associated with gallbladder cancer via gene expression profile analysis. J Cell Biochem. (2019) 120:19377–87. doi: 10.1002/jcb.28953
Keywords: m6A methylation, methyltransferase, demethylases, epigenetics, prognostic signature
Citation: Wang J, Zhang C, He W and Gou X (2020) Effect of m6A RNA Methylation Regulators on Malignant Progression and Prognosis in Renal Clear Cell Carcinoma. Front. Oncol. 10:3. doi: 10.3389/fonc.2020.00003
Received: 14 July 2019; Accepted: 06 January 2020;
Published: 24 January 2020.
Edited by:Xiaojie Tan, Second Military Medical University, China
Copyright © 2020 Wang, Zhang, He and Gou. 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: Chengyao Zhang, firstname.lastname@example.org