Identification of METTL14 in Kidney Renal Clear Cell Carcinoma Using Bioinformatics Analysis

Wang, Qian; Zhang, Hao; Chen, Quanbing; Wan, Zhenghua; Gao, Xiaoyong; Qian, Wenhui

doi:https://doi.org/10.1155/2019/5648783

Disease Markers

On this page

Abstract Introduction Materials and Methods Results Discussion Conclusion Data Availability Conflicts of Interest Authors’ Contributions Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2019 | Article ID 5648783 | https://doi.org/10.1155/2019/5648783

Identification of METTL14 in Kidney Renal Clear Cell Carcinoma Using Bioinformatics Analysis

Qian Wang,¹Hao Zhang,²Quanbing Chen,³Zhenghua Wan,⁴Xiaoyong Gao,³and Wenhui Qian³

Academic Editor: Małgorzata Knaś

Received18 Aug 2019

Accepted04 Nov 2019

Published30 Dec 2019

Abstract

The kidney renal clear cell carcinoma (KIRC) with poor prognosis is the main histological subtype of the renal cell carcinoma, accounting for 80–90% of patients. Currently, the N6-methyladenosine (m6A) epitranscriptional modification draws much attention. The m6A RNA modification, the most plentiful internal modification of mRNAs and noncoding RNAs in the majority of eukaryotes, regulates mRNAs at different levels and is involved in disease occurrence and progression. The GTExPortal and TCGAportal were applied to investigate the METTL14 mRNA expression in different tissues and KIRC stages. The Human Protein Atlas was used to verify the location of METTL14 in KIRC tissues. The main microRNAs (miRNAs) related to KIRC were analyzed using OncoLnc and starBase, while corresponding circular RNAs (circRNAs) interacting with miRNAs were predicted via circBank; then, the METTL14-miRNA-circRNA interaction network was established. The level of methyltransferase-like 14 (METTL14) mRNA was significantly lower in KIRC tissues compared with normal kidney tissues, which was relative to clinical and pathological stages. circRNAs may regulate METTL14 mRNA as miRNAs sponge to affect the KIRC progression. METTL14 mRNA is likely to regulate PTEN mRNA expression via changing its m6A RNA modification level. METTL14 mRNA expression negatively correlated with the KIRC stages and positively correlated with KIRC patients’ overall survival, which has great potential to serve as a clinical biomarker in KIRC.

1. Introduction

The renal cell carcinoma (RCC) is the most common malignant tumor derived from the kidney. There are approximately 210,000 new cases diagnosed all over the world each year, accounting for 2–3% of all cancers. Currently, the kidney renal clear cell carcinoma (KIRC) is the main histological subtype of RCC and accounts for 80–90% of RCC patients. Yet the prognosis of KIRC is poor [1–3]. At present, the TNM stage and Fuhrman grade are usually used as predictors to assess the risk of patients with KIRC [4, 5]. However, sometimes the outcomes may still be variable for patients with similar clinical features or scores. Meanwhile, a lot of studies indicated that KIRC was heterogeneous in various aspects including clinicopathological, molecular, and cellular heterogeneity. Therefore, it is required and urgent to study the mechanism of KIRC and explore new useful molecular markers for diagnosis and prognosis.

There have been more than 100 kinds of RNA modifications identified in living organisms [6]. Some researches have widely reported certain types of RNA modifications in eukaryotic mRNA, containing N1-methyladenosine (m1A), N6-methyladenosine (m6A), and 5-methylcytosine (m5C), among which m6A was first discovered in the 1970s. Currently, the m6A epitranscriptional modification draws much attention. The m6A RNA modification is the most plentiful internal modification of mRNA and noncoding RNA in the majority of eukaryotes, as well as appreciably clusters around the stop codon and 3 untranslated region (3UTR) [7–9]. The physiological importance of m6A RNA modification has been confirmed by its pivotal roles in tissue development and differentiation [10–16]. m6A RNA modification regulates mRNA at different levels, such as structure, maturation, stability, splicing, export, translation, and decay [17]. Moreover, m6A RNA modification is also involved in cell fate decision, cell cycle regulation, cell differentiation, and circadian rhythm maintenance [18]. Numerous RNA-binding proteins (RBPs) are also affected by m6A RNA modification.

In this study, we explored the relationship between methyltransferase-like 14 (METTL14) and KIRC. We used UALCAN to get a heatmap of gene expression related to m6A RNA modification in KIRC and then found that METTL14 mRNA showed a differential expression in KIRC tissues compared with normal tissues. The formation of m6A RNA modification is catalyzed by the RNA methyltransferase complex including methyltransferase-like 3 (METTL3), METTL14, and Wilms tumor 1-associated protein (WTAP), which are the writers of the m6A marks [19]. GTExPortal, TCGAportal, and UALCAN were applied to investigate the expression of METTL14 mRNA in different tissues and KIRC stages.

Additionally, the Human Protein Atlas was used to verify the location of METTL14 in KIRC tissues. MicroRNAs (miRNAs) related to KIRC in PITA, miRmap, and miRanda and corresponding circular RNAs (circRNAs) interacting with miRNAs were analyzed via circBank; then, the METTL14-miRNA-circRNA interaction network was established. Further, both phosphatase and tensin homologue (PTEN) and eukaryotic translation initiation factor 3 subunit A (EIF3A) were found to be positively correlated with METTL14 in KIRC. PTEN is reported as a well-known tumor suppressor gene. EIF3A is a protein encoded by the EIF3A gene in human, which plays crucial roles in regulating translation of a subset of mRNAs, cell cycle progression, and cell proliferation [20]. In this study, PTEN and EIF3A both were found to have a positive correlation with METTL14 in KIRC.

2. Materials and Methods

2.1. Identification of Genes

The heatmap of gene expressions related to m6A in KIRC was explored in UALCAN (http://ualcan.path.uab.edu/index.html). The expressions of METTL14 mRNA in different tissues were further probed in GTExPortal and UALCAN. In addition, the Human Protein Atlas was used to investigate the location of METTL14 based on pathological slides of KIRC patients. The correlations between expressions of genes were analyzed in GEPIA 2 (http://gepia2.cancer-pku.cn/#index).

2.2. The Relationship between METTL14 and KIRC

Via TCGAportal (http://www.tcgaportal.org/), we investigated the expression level of METTL14 mRNA in KIRC stages, mutation types, and corresponding overall survival of KIRC patients. Then, MEXPRESS (http://mexpress.be/) was used to assess the pathological types of KIRC which were related to genders.

2.3. The Predictions of the Potential Mechanism

We found the m6A motif base sequence in which METTL14 binds to mRNA and corresponding location using RMBase v2.0 (http://rna.sysu.edu.cn/rmbase/). A total of 9 hub genes related to KIRC and the relationship between them were obtained from TCGAportal. This website was also used to investigate the correlation between METTL14, PTEN, and EIF3A mRNA expressions.

2.4. The Identification of Functional Analysis of METTL14

PITA, miRmap, and miRanda, three online tools, were used to explore miRNAs which could bind to METTL14 mRNA. Then, three final sets were merged to pick out common 14 miRNAs (e.g., miR-361, miR-20b, miR-17, miR-106b, and miR-93). OncoLnc (http://www.oncolnc.org/) and starBase (http://starbase.sysu.edu.cn/) were used to further analyze the 14 miRNAs, and finally, the 4 miRNAs were selected. circBank (http://www.circbank.cn/) was applied to find a total of 24 circRNAs interacting with the 4 miRNAs. Subsequently, we constructed the abovementioned miRNAs and circRNAs using Cytoscape to display their interaction network with METTL14 mRNA. Protein interactions among METTL14 and relative proteins were predicted on STRING (http://string-db.org/). GO function analysis for METTL14 mRNA was performed to investigate the molecular function and biological process of METTL14.

3. Results

3.1. METTL14 mRNA Was Significantly Downregulated in KIRC

The heatmap of gene expression related to m6A in KIRC showed that METTL14 mRNA has lower expression in KIRC tissues compared with adjacent normal tissues (Figure 1(a)). The METTL14 mRNA expression in different normal tissues, besides kidney tissues, was further probed to find that METTL14 mRNA has generally higher expression in normal tissues (Figure 1(b)). Additionally, expression of METTL14 mRNA was further confirmed to be significantly downregulated in KIRC tissues compared with adjacent normal kidney tissues (Figure 1(c)). Further, METTL14 was mainly located in the nucleus by immunohistochemistry (Figure 1(d)).

(a)

(b)

(c)

(d)

3.2. Patients with Lower METTL14 Expression Had a Worse Prognosis

The METTL14 mRNA expression decreased as the KIRC clinical stage progressed (Figure 2(a)). In addition, the METTL14 mRNA expression has a significantly negative correlation with KIRC pathological stages with male patients more prone to advanced pathological stages by MEXPRESS (Figure 2(b)). Importantly, KIRC patients with lower METTL14 mRNA expression had shorter overall survival (Figure 2(c)). In additionally, the in-frame mutation was the least among mutation types of METTL14 mRNA in KIRC, while the numbers of the other 4 types were close to the quantity of no mutation (Figure 2(d)).

(a)

(b)

(c)

(d)

3.3. The Potential Mechanism of METTL14 mRNA in KIRC

A total of 9 hub genes related to KIRC were listed, such as VHL, PMRM1, and SETD2. Among these genes, the mutation rate of PTEN mRNA increased as the METTL14 mRNA expression decreased (Figure 3(a)). Further, as the METTL14 mRNA expression increased, the PTEN mRNA expression in KIRC exhibited significant upregulation (Figure 3(b)). METTL14 could act as an m6A writer to affect the stability of some mRNAs and further regulate mRNA expressions [21]. The m6A motif base sequence in which METTL14 binds to mRNA was mainly GGACU, which was located in the junction between the coding sequence (CDS) and 3-untranslated region (3-UTR) of mRNAs (Figures 3(c) and 3(d)). EIF3A could act as an m6A reader to participate in m6A RNA modification progression [22]. EIF3A mRNA expression in KIRC tissues was lower than that in normal kidney tissues. A positive correlation between PTEN and METTL14 mRNA in KIRC was observed. We also noticed that as the METTL14 mRNA expression increased, the EIF3A mRNA expression in KIRC exhibited significant upregulation (Figure S1). In addition, patients with lower EIF3A expression had shorter overall survival (Figure 3(e)).

(a)

(b)

(c)

(d)

(e)

Figure 3

The potential mechanism of METTL14 in KIRC. (a) The relevance between METTL14 mRNA expression and the mutation rate of 9 hub genes in KIRC. (b) The expression of PTEN mRNA and METTL14 mRNA in KIRC had a significantly positive correlation. (c) The m6A motif base sequence in which METTL14 binds to mRNA was more likely to be GGACU (RMBase v2.0). (d) The m6A motif base sequence was located in the junction between CDS and 3UTR of mRNA with higher probability. The location of GGACU sequence in PTEN mRNA. (e) (A) The expression level of EIF3A mRNA had significantly differential expression in KIRC versus normal tissues. (B) The expression of EIF3A mRNA and PTEN mRNA in KIRC had a significantly positive correlation. (C, D) The TCGA survival analysis for EIF3A.

3.4. Functional and Pathway Enrichment Analysis

PITA, miRmap, and miRanda were first used to explore some miRNAs binding to METTL14 mRNA separately, and the three final sets were merged to get common 14 miRNAs (e.g., miR-361, miR-20b, miR-17, miR-106b, and miR-93) (Figure 4(a)). The further analysis of 14 miRNAs showed that the expression of miR-130a-3p, miR-130b-3p, miR-106b-5p, and miR-301a-3p had a significantly negative correlation with METTL14 mRNA, respectively, in KIRC (Figure 4(b)). The KIRC patients, who had higher miR-130a-3p, miR-130b-3p, miR-106b-5p, and miR-301a-3p expression level, had shorter overall survival (Figure 4(c)).

(a)

(b)

(c)

(d)

(e)

A total of 24 circRNAs interacting with the 4 miRNAs were exhibited. Then, an interaction network was constructed based on the abovementioned miRNAs and circRNAs (Figure 4(d)). Protein interaction among METTL14 and relative proteins indicated that METTL14 had potential to play a role by regulating the expression of proteins (Figure 4(e)).

3.5. Molecular Function and Biological Process

The results of GO function analysis for proteins interacting with METTL14 demonstrated that the top five molecular functions are DNA-directed RNA polymerase activity, nucleic acid binding, mRNA (2-O-methyladenosine-N6-)-methyltransferase activity, transferase activity, and RNA binding (Table 1), as well as the top five biological processes including 7-methylguanosine mRNA capping, positive regulation of viral transcription, transcription elongation from RNA polymerase II promoter, RNA splicing, and mRNA splicing via the spliceosome (Table 2).

4. Discussion

The renal cell carcinoma is a heterogeneous disease, among which the kidney renal clear cell carcinoma was one of the highest incidence subtypes with poor prognosis. Due to limited predictors assessing the risk, a part of KIRC patients might be diagnosed as having inaccurate grade, which would decrease the survival rate of patients with KIRC to some extent. Hence, it is crucial to identify new specific predictors for KIRC.

As a tumor suppressor, PTEN mutation happens frequently in a large number of cancers. The protein encoded by PTEN is phosphatidylinositol-3,4,5-trisphosphate 3-phosphatase, which contains a catalytic domain as well as a tensin-like domain similar to that of the dual-specificity protein tyrosine phosphatases. This protein preferentially dephosphorylates phosphoinositide substrates, which is different from most of the protein tyrosine phosphatases. It negatively regulates intracellular levels of phosphatidylinositol-3,4,5-trisphosphate and functions as a tumor suppressor by negatively regulating the AKT/PKB signaling pathway. METTL14 had a positive correlation with the tumor suppressor gene PTEN, so we speculated that METTL14 exerted a tumor suppressor effect in KIRC by stabilizing PTEN mRNA.

In the previous study, the formation of m6A RNA modification is a reversible process [23]. m6A “writers” with methyltransferase activity consisted of three individual proteins: METTL3, METTL14, and WTAP. FTO and alkB homolog 5 (ALKBH5), m6A demethylases, are called “erasers” [24, 25]. Another protein family is m6A “readers” which can recognize m6A RNA modification to modulate mRNA fate [26]. EIF3A is one of the m6A readers.

In our study, the mRNA expressions of PTEN and EIF3A both had a positive correlation with the expression of METTL14 mRNA in KIRC, suggesting that the interaction between EIF3A and METTL14 may produce synergistic effects and hence jointly regulate KIRC progression. We consider that proteins related to m6A RNA modification have regulatory roles in each other, such as readers in writers and erasers in readers. Additionally, we also doubt whether those proteins are regulated by other modifications like m1A and m5C. The potential mechanism needs further research.

There is another discovery worth talking. The four important functional miRNAs in this study have not been reported yet in KIRC. Given their meaningful characteristics in KIRC, it definitely deserves further exploring and researching.

5. Conclusion

Bioinformatics analysis demonstrated that downregulated METTL14 mRNA expression has a significantly negative correlation with the KIRC stages and a positive correlation with overall survival of KIRC patients. Furthermore, circRNAs may regulate METTL14 mRNA as miRNAs sponge to affect the progression of KIRC. METTL14 mRNA is likely to regulate PTEN mRNA expression via changing its m6A RNA modification level. The relationship between KIRC and m6A RNA modification is closely related. In conclusion, METTL14 has great potential to work as a clinical biomarker for KIRC in the future. Further study is needed to explore the value of METTL14 in the stages and prognosis of KIRC.

Data Availability

The data that support the findings of this study are available from the corresponding authors upon reasonable request.

Conflicts of Interest

The authors declare that they have no competing interests.

Authors’ Contributions

Qian Wang and Hao Zhang contributed equally to this work.

Supplementary Materials

Figure S1: the correlation between the expressions of EIF3A mRNA and METTL14 mRNA in KIRC. (Supplementary Materials)

References

B. Ljungberg, K. Bensalah, S. Canfield et al., “EAU guidelines on renal cell carcinoma: 2014 update,” European Urology, vol. 67, no. 5, pp. 913–924, 2015.
View at: Publisher Site | Google Scholar
A. M. Forsea, V. D. Marmol, E. D. Vries, E. E. Bailey, and A. C. Geller, “Melanoma incidence and mortality in Europe: new estimates, persistent disparities,” British Journal of Dermatology, vol. 167, no. 5, pp. 1124–1130, 2012.
View at: Publisher Site | Google Scholar
C. L. Chaffer and R. A. Weinberg, “A perspective on cancer cell metastasis,” Science, vol. 331, no. 6024, pp. 1559–1564, 2011.
View at: Publisher Site | Google Scholar
V. Ficarra, G. Novara, A. Galfano et al., “The ‘stage, size, grade and necrosis’ score is more accurate than the University of California Los Angeles Integrated Staging System for predicting cancer-specific survival in patients with clear cell renal cell carcinoma,” British Journal of Urology International, vol. 103, no. 2, pp. 165–170, 2010.
View at: Publisher Site | Google Scholar
R. Zigeuner and H. Gchromecki, “External validation of the Mayo Clinic stage, size, grade, and necrosis (SSIGN) score for clear-cell renal cell carcinoma in a single European centre applying routine pathology,” European Urology, vol. 57, no. 1, pp. 102–109, 2010.
View at: Publisher Site | Google Scholar
P. Boccaletto, M. A. Machnicka, E. Purta et al., “MODOMICS: a database of RNA modification pathways. 2017 update,” Nucleic Acids Research, vol. 46, no. D1, pp. D303–D307, 2018.
View at: Publisher Site | Google Scholar
D. Dominissini, S. Moshitch-Moshkovitz, S. Schwartz et al., “Topology of the human and mouse m6A RNA methylomes revealed by m6A-seq,” Nature, vol. 485, no. 7397, pp. 201–206, 2012.
View at: Publisher Site | Google Scholar
K. D. Meyer, Y. Saletore, P. Zumbo, O. Elemento, C. E. Mason, and S. R. Jaffrey, “Comprehensive analysis of mRNA methylation reveals enrichment in 3 UTRs and near stop codons,” Cell, vol. 149, no. 7, pp. 1635–1646, 2012.
View at: Publisher Site | Google Scholar
Z. Bodi, A. Bottley, N. Archer, S. T. May, and R. G. Fray, “Yeast m6A methylated mRNAs are enriched on translating ribosomes during meiosis, and under rapamycin treatment,” PLoS One, vol. 10, no. 7, article e0132090, 2015.
View at: Publisher Site | Google Scholar
Z. Lin, P. J. Hsu, X. Xing et al., “Mettl3-/Mettl14-mediated mRNA N (6)-methyladenosine modulates murine spermatogenesis,” Cell Research, vol. 27, no. 10, pp. 1216–1230, 2017.
View at: Publisher Site | Google Scholar
J. Lv, Y. Zhang, S. Gao et al., “Endothelial-specific m6A modulates mouse hematopoietic stem and progenitor cell development via Notch signaling,” Cell Research, vol. 28, no. 2, article cr2017143, pp. 249–252, 2017.
View at: Publisher Site | Google Scholar
L. P. Vu, B. F. Pickering, Y. Cheng et al., “The _N_ ⁶-methyladenosine (m⁶A)-forming enzyme METTL3 controls myeloid differentiation of normal hematopoietic and leukemia cells,” Nature Medicine, vol. 23, no. 11, pp. 1369–1376, 2017.
View at: Publisher Site | Google Scholar
J. Yu, M. Chen, H. Huang et al., “Dynamic m6A modification regulates local translation of mRNA in axons,” Nucleic Acids Research, vol. 46, no. 3, pp. 1412–1423, 2018.
View at: Publisher Site | Google Scholar
C. Zhang, Y. Chen, B. Sun et al., “m6A modulates haematopoietic stem and progenitor cell specification,” Nature, vol. 549, no. 7671, pp. 273–276, 2017.
View at: Publisher Site | Google Scholar
K. Xu, Y. Yang, G. H. Feng et al., “Mettl3-mediated m(6) A regulates spermatogonial differentiation and meiosis initiation,” Cell Research, vol. 27, no. 9, pp. 1100–1114, 2017.
View at: Publisher Site | Google Scholar
C. X. Wang, G. S. Cui, X. Liu et al., “METTL3-mediated m⁶A modification is required for cerebellar development,” PLOS Biology, vol. 16, no. 6, p. e2004880, 2018.
View at: Publisher Site | Google Scholar
Z. Liu and J. Zhang, “Human C-to-U coding RNA editing is largely nonadaptive,” Molecular Biology and Evolution, vol. 35, no. 4, pp. 963–969, 2018.
View at: Publisher Site | Google Scholar
Y. Wu, C. Zhou, and Q. Yuan, “Role of DNA and RNA N6-adenine methylation in regulating stem cell fate,” Current Stem Cell Research & Therapy, vol. 13, no. 1, 2017.
View at: Publisher Site | Google Scholar
Y. Wang, Y. Li, J. I. Toth, M. D. Petroski, Z. Zhang, and J. C. Zhao, “N6-methyladenosine modification destabilizes developmental regulators in embryonic stem cells,” Nature Cell Biology, vol. 16, no. 2, pp. 191–198, 2014.
View at: Publisher Site | Google Scholar
G. Ren, L. Yu, O. Yang, S. Liu, and D. Li, “The role of miR-497/EIF3A axis in TGFβ1-induced epithelial-mesenchymal transition and extracellular matrix in rat alveolar epithelial cells and pulmonary fibroblasts,” Journal of Cellular Biochemistry, vol. 118, no. 10, pp. 3401–3408, 2017.
View at: Publisher Site | Google Scholar
H. Weng, H. Huang, H. Wu et al., “METTL14 inhibits hematopoietic stem/progenitor differentiation and promotes leukemogenesis via mRNA m6A modification,” Cell Stem Cell, vol. 22, no. 2, pp. 191–205.e9, 2017.
View at: Publisher Site | Google Scholar
J. Tong, R. A. Flavell, and H. B. Li, “RNA m6A modification and its function in diseases,” Frontiers of Medicine, vol. 12, no. 4, pp. 481–489, 2018.
View at: Publisher Site | Google Scholar
G. Jia, Y. Fu, and C. He, “Reversible RNA adenosine methylation in biological regulation,” Trends in Genetics, vol. 29, no. 2, pp. 108–115, 2013.
View at: Publisher Site | Google Scholar
G. Jia, Y. Fu, X. Zhao et al., “N6-methyladenosine in nuclear RNA is a major substrate of the obesity-associated FTO,” Nature Chemical Biology, vol. 7, no. 12, pp. 885–887, 2011.
View at: Publisher Site | Google Scholar
G. Zheng, J. A. Dahl, Y. Niu et al., “ALKBH5 is a mammalian RNA demethylase that impacts RNA metabolism and mouse fertility,” Molecular Cell, vol. 49, no. 1, pp. 18–29, 2013.
View at: Publisher Site | Google Scholar
A. Li, Y. S. Chen, X. L. Ping et al., “Cytoplasmic m6A reader YTHDF3 promotes mRNA translation,” Cell Research, vol. 27, no. 3, pp. 444–447, 2017.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2019 Qian Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

1089

Downloads

1177

Citations