Analysis of METTL3 and METTL14 in hepatocellular carcinoma

Xiangxiang Liu; Jian Qin; Tianyi Gao; Chenmeng Li; Xiaoxiang Chen; Kaixuan Zeng; Mu Xu; Bangshun He; Bei Pan; Xueni Xu; Yuqin Pan; Huiling Sun; Tao Xu; Shukui Wang

doi:10.18632/aging.103959

Research Paper Volume 12, Issue 21 pp 21638—21659

Analysis of METTL3 and METTL14 in hepatocellular carcinoma

Xiangxiang Liu^1, , Jian Qin^2, , Tianyi Gao^2, , Chenmeng Li^1, , Xiaoxiang Chen^2, , Kaixuan Zeng^1, , Mu Xu^2, , Bangshun He^2, , Bei Pan^2, , Xueni Xu^1, , Yuqin Pan^2, , Huiling Sun^2, , Tao Xu^2, , Shukui Wang^1,3, ,

¹ School of Medicine, Southeast University, Nanjing 210096, Jiangsu, China
² General Clinical Research Center, Nanjing First Hospital, Nanjing Medical University, Nanjing 210006, Jiangsu, China
³ Jiangsu Collaborative Innovation Center on Cancer Personalized Medicine, Nanjing Medical University, Nanjing 211100, Jiangsu, China

Received: April 10, 2020 Accepted: August 1, 2020 Published: November 6, 2020

https://doi.org/10.18632/aging.103959
How to Cite

Copyright: © 2020 Liu et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Abstract

N6-methyladenosine (m6A) RNA methylation is the most prevalent modification of messenger RNAs (mRNAs) and catalyzed by a multicomponent methyltransferase complex (MTC), among which methyltransferase-like 3 (METTL3) and METTL14 are two core molecules. However, METTL3 and METTL14 play opposite regulatory roles in hepatocellular carcinoma (HCC). Based on The Cancer Genome Atlas (TCGA) database and Gene Expression Omnibus (GEO) database, we conducted a multi-omics analysis of METTL3 and METTL14 in HCC, including RNA-sequencing, m6ARIP-sequencing, and ribosome-sequencing profiles. We found that the expression and prognostic value of METTL3 and METTL14 are opposite in HCC. Besides, after METTL3 and METTL14 knockdown, most of the dysregulated mRNAs, signaling pathways and biological processes are distinct in HCC, which partly explains the contrary regulatory role of METTL3 and METTL14. Intriguingly, these mRNAs whose stability or translation efficiency are influenced by METTL3 or METTL14 in an m6A dependent manner, jointly regulate multiple signaling pathways and biological processes, which supports the cooperative role of METTL3 and METTL14 in catalyzing m6A modification. In conclusion, our study further clarified the contradictory role of METTL3 and METTL14 in HCC.

Introduction

N6-methyladenosine (m6A) RNA methylation, the most prevalent modification of messenger RNAs (mRNAs), accounts for almost half of the total methylated ribonucleotides [1]. In general, m6A modification is present in the transcripts of over 7,000 genes in mammalian cells, and it prefers to occur at the consensus RRACH motif (R = G or A; H = A, C, or U). Transcriptome-wide m6A site mapping reveals more details on its localization which preferentially enriches at coding sequence (CDS), around stop codon, and 3’untranslated region (3’UTR) in the transcriptomes [2, 3]. The formation of m6A is catalyzed by a multicomponent methyltransferase complex (MTC), among which methyltransferase-like 3 (METTL3), METTL14, and Wilms’ tumor 1-associating protein (WTAP), Vir Like M6A Methyltransferase Associated (KIAA1429), RNA Binding Motif Protein 15 (RBM15), and Zinc Finger CCCH-Type Containing 13 (ZC3H13) have been detected [4–8]. The m6A is a reversible modification that can also be removed by RNA demethylases, including fat mass and obesity-associated protein (FTO) and alkylated DNA repair protein alkB homolog 5 (ALKBH5) [9, 10]. In addition, m6A modified mRNAs can be bound by multiple specific RNA binding proteins, of which the known ones are YTH Domain-Containing Protein 1 (YTHDC1), YTH Domain Family, Member 1/2/3 (YTHDF1/2/3), Insulin-Like Growth Factor 2 MRNA Binding Protein 1/2/3 (IGF2BP1/2/3) [11–15].

METTL3 and METTL14, two core components of MTC, colocalize in nuclear speckles, and catalyze the covalent transfer of a methyl group to adenine in a heterodimer form [16]. The unusual m6A modification caused by differentially expressed METTL13 or METTL14 plays a critical role in the malignant progression of various cancers, such as bladder cancer, gastric cancer, and hepatocellular carcinoma (HCC) [17–19]. Interestingly, emerging evidence indicated an opposite regulatory role of METTL3 and METTL14 in several cancers, such as glioblastoma [20, 21], HCC [19, 22], and colorectal cancer (CRC) [23, 24]. METTL3 was demonstrated to be upregulated in CRC and to facilitate CRC progression by maintaining the SRY-Box 2 (SOX2) mRNA stability [23]. However, our recent study proved that METTL14 is significantly downregulated in CRC and suppresses cell growth and metastasis via regulating primary miR-375 processing [25]. Consistently, METTL3 serves as an oncogene, but METTL14 is a tumor suppressor in HCC [19, 22]. Study in vitro showed that both METTL3 and METTL14 have a methyltransferase domain to methylate RNA and the m6A methyltransferase activity is much higher in METTL3/METTL14 complex than that in either subunit alone [4]. However, the crystal structure and biochemical evidence suggested that METTL3, rather than METTL14, is the unique catalytic subunit, and METTL14 functions in structural stabilization and RNA substrate recognition [26]. Therefore, it is urgently needed to further clarify the characterization of METTL3 and METTL14 in cancers.

In this study, based on The Cancer Genome Atlas (TCGA) database and Gene Expression Omnibus (GEO) database, we conducted a multi-omics analysis of METTL3 and METTL14 in HCC. We validated the reverse expression and prognostic value of METTL3 and METTL14 in HCC and discovered that most of the mRNAs and associated signaling pathways and biological processes regulated by METTL3 and METTL14 are different, which may be partly responsible for their contrary functions in HCC. However, these mRNAs whose stability or translation efficiency (TE) are affected by METTL3 or METTL14 in an m6A dependent manner, jointly regulate multiple signaling pathways and biological processes, such as TGF-beta signaling pathway, protein ubiquitination, and cell cycle.

Results

The opposite expression and prognostic value of METTL3 and METTL14 in HCC

To verify the expression of METTL3 and METTL14 in HCC, we analyzed the TCGA database and two GEO datasets. TCGA database and GSE14520 [27] analysis showed an increased expression of METTL3 in HCC tissues compared to normal liver tissues (NTs) (Figure 1A, 1B). In contrast, The TCGA database and GSE54236 [28] analysis exhibited a downregulated expression of METTL14 in HCC tissues (Figure 1C and 1D). IHC analysis validated the overexpressed METTL3 and downregulated METTL14 expression in HCC (Figure 1E). Besides, we discovered that HCC patients with higher METTL3 expression have shorter overall survival (OS) time, relapse-free survival (RFS) time, progression-free survival (PFS) time, and disease-specific survival (DSS) time compared to those with low METTL3 expression (Figure 1F). Conversely, HCC patients with low METTL14 expression undergo poorer OS rate, RFS rate, PFS rate, and DSS rate compared to those with high METTL14 expression (Figure 1G). These results demonstrated an opposite expression and prognostic value of METTL3 and METTL14 in HCC.

Figure 1. The opposite expression and prognostic value of METTL3 and METTL14 in HCC. (A, B) The expression of METTL3 in HCC tissues and NTs based on the TCGA database and GSE14520 analysis. (C, D) The expression of METTL14 in HCC tissues and NTs based on the TCGA database and GSE14520 analysis. (E) IHC analysis of METTL3 and METTL14 in HCC tissues and adjacent NTs. (F) The associations between METTL3 expression and OS, RFS, PFS, and DSS of HCC patients. (G) The associations between METTL14 expression and OS, RFS, PFS, and DSS of HCC patients. * p<0.05, *** p<0.001.

These differentially expressed genes regulated by METTL3 and METTL14 knockdown participate in different signaling pathways and biological processes

To explore the regulatory role of METTL3 and METTL14 on mRNA expression, we analyzed GSE90642 and GSE37001 datasets [2, 15] The RNA sequencing data of GSE90642 showed that 329 and 530 mRNAs were downregulated and overexpressed in HepG2 cells after METTL14 knockdown, respectively (Figure 2A). In addition, METTL14 knockdown resulted in 705 mRNAs with downregulated m6A modification. Interestingly, 270 mRNAs with upregulated m6A modification were also observed (Figure 2A). GSE37001 analysis revealed that 1147 downregulated mRNAs and 812 overexpressed mRNAs after METTL3 knockdown in HepG2 cells (Figure 2A). Intriguingly, when overlapping these differentially expressed genes (DEGs), we found that only a small number of DEGs (n=101) were co-regulated by METTL3 and METTL14, of which 51% even showed an opposite expression (Figure 2C). Thereafter, we conducted a KEGG analysis of these DEGs. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that the DEGs regulated by METTL14 (M14DEGs) were significantly enriched in ten signaling pathways, such as MAPK, PI3K-Akt, and TGF-beta signaling pathways (Figure 2D). Meanwhile, these DEGs regulated by METTL3 (M3DEGs) were significantly enriched in twenty signaling pathways, such as p53, TNF, and TGF-beta signaling pathways (Figure 2E). As expected, only three signaling pathways were jointly regulated by METTL3 and METTL14, including MAPK, TNF, and TGF-beta signaling pathways. Subsequently, we carried out Gene Ontology (GO) analysis of M14DEGs and M3DEGs, respectively (Additional file 1: Supplementary Table 1, 2). In line with KEGG analysis results, only a small part of biological processes (n=22) was jointly regulated by METTL3 and METTL14, such as cell migration (GO:0016477) (Figure 2F, 2G). Taken together, we concluded that most M3DEGs and M14DEGs as well as associated signaling pathways and biological processes are different in HCC.

The DEGs regulated by METTL3 and METTL14 knockdown participate in different signaling pathways and biological processes. (A, B) The DEGs regulated by METTL3 and METTL14 knockdown. (C) Integrated analysis of M3DEGs and M14DEGs. (D, E) KEGG pathway analysis of M3DEGs and M14DEGs. (F) The number of biological processes regulated by M3DEGs and M14DEGs. (G) The common biological processes regulated by M3DEGs and M14DEGs.

Figure 2. The DEGs regulated by METTL3 and METTL14 knockdown participate in different signaling pathways and biological processes. (A, B) The DEGs regulated by METTL3 and METTL14 knockdown. (C) Integrated analysis of M3DEGs and M14DEGs. (D, E) KEGG pathway analysis of M3DEGs and M14DEGs. (F) The number of biological processes regulated by M3DEGs and M14DEGs. (G) The common biological processes regulated by M3DEGs and M14DEGs.

The signaling pathways and biological processes of m6A modified M3DEGs

To investigate whether those M3DEGs were m6A modified genes, we first analyzed m6ARIP-sequencing data in HepG2 cells by using GSE37003 [2] which revealed a total of 12424 m6A peaks in 7180 transcripts. Chromosome location analysis showed that these m6A peaks were enriched in all chromosomes (Figure 3A). Besides, most transcripts (82%) had only one or two m6A peaks (Figure 3B). Consistent with previous reports, most m6A peaks (85%) located in CDS, around stop codon, and 3' UTR of mRNAs (Figure 3C). Unexpectedly, when we overlapped m6A modified transcripts with M3DEGs, we found that only 37% of M3DEGs had m6A peaks (Figure 3D). KEGG pathway analysis showed that these m6A modified M3DEGs were enriched in 18 signaling pathways, such as p53 signaling pathway, signaling pathways regulating pluripotency of stem cells, and TGF-beta signaling pathway (Figure 3E). After that, we conducted a protein-protein interaction (PPI) network of these m6A modified M3DEGs which detected four functional molecular clusters (Figure 3F). KEGG pathway analysis of these clusters exhibited three common signaling pathways, including cell cycle and p53 signaling pathway of cluster 2, Glycosaminoglycan biosynthesis-chondroitin sulfate/dermatan sulfate of cluster 3. In addition, Cluster 1 and cluster 4 revealed two unique pathways, including Ubiquitin mediated proteolysis and the MAPK signaling pathway (Figure 3F). Consistent with KEGG pathway analysis, GO enrichment analysis showed that cluster 1 was involved in protein polyubiquitination and ubiquitination, that cluster 2 was involved in the cell cycle, that cluster 3 was involved in dermatan sulfate biosynthetic process, and that cluster 4 was involved in MAPK activity (Figure 3G). Moreover, all top 30 hub genes among m6A modified M3DEGs belonged to cluster 1 and cluster 2 (Figure 3H). Meanwhile, the expression of most hub genes was positively correlated with METTL3 expression in HCC (87%) and significantly associated with the OS of HCC patients (73%) (Additional file 2: Supplementary Table 3, 4). Based on the above data, we believed that METTL3 directly regulates a small number of genes’ expression in an m6A dependent manner and that these m6A modified M3DEGs are mainly involved in protein ubiquitination and cell cycle.

Figure 3. The signaling pathways and biological processes of m6A modified M3DEGs. (A) The chromosome location of transcripts having m6A peaks in HepG2 cell. (B) The number of genes with various m6A peaks. (C) The distribution of m6A peaks in transcripts. (D) Integrated analysis of M3DEGs and mRNAs with m6A peaks. (E) KEGG pathway analysis of m6A modified M3DEGs. (F, G) KEGG pathway and GO enrichment analysis of functional molecular clusters among the PPI network of m6A modified M3DEGs. (H) The top 30 hub genes among m6A modified M3DEGs.

The signaling pathways and biological processes of m6A modified M14DEGs

To analyze the m6A modified M14DEGs, we overlapped m6A modified genes with M14DEGs. We found that 54% of M14DEGs were m6A modified mRNAs (Figure 4A). KEGG pathway analysis showed that these m6A modified M14DEGs were enriched in only 4 signaling pathways, including the MAPK signaling pathway, Hippo signaling pathway, Endocytosis, and TGF-beta signaling pathway (Figure 4B). PPI network of these m6A modified M14DEGs revealed three functional molecular clusters (Figure 4C). KEGG pathway analysis of cluster 3 showed a common signaling pathway, Endocytosis. Interestingly, cluster 1 and cluster 3 revealed three unique pathways, including Ribosome biogenesis in eukaryotes, RNA transport, and Focal adhesion (Figure 4C). GO enrichment analysis showed that cluster 1 was associated with protein polyubiquitination and ubiquitination, and that cluster 2 was associated with regulation of epithelial cell proliferation, and that cluster 3 was associated with cell proliferation, migration, and adhesion (Figure 4D). Furthermore, all top 30 hub genes among m6A modified M14DEGs belonged to cluster 1 and cluster 2 (Figure 4E). Moreover, consistently, the expression of most hub genes was positively correlated with METTL14 expression in HCC (73%) and significantly associated with the OS of HCC patients (50%) (Additional file 3: Supplementary Table 5, 6). Overall, our results indicated that METTL14, liking METTL3, directly regulates part of genes’ expression in an m6A dependent manner and that these m6A modified M14DEGs mainly participate in protein ubiquitination and cell proliferation.

Figure 4. The signaling pathways and biological processes of m6A modified M14DEGs. (A) Integrated analysis of M14DEGs and mRNAs with m6A peaks. (B) KEGG pathway analysis of m6A modified M14DEGs. (C, D) KEGG pathway and GO enrichment analysis of functional molecular clusters among the PPI network of m6A modified M14DEGs. (E) The top 30 hub genes among m6A modified M14DEGs.

These genes with changed TE regulated by METTL3 and METTL14 knockdown participate in distinct signaling pathways and biological processes

To further investigate the regulatory role of METTL3 and METTL14 on mRNA translation, we analyzed GSE63591 [12] and GSE121952 [29] datasets. Ribosome sequencing data of GSE63591 showed that the TE of 690 mRNAs was significantly downregulated while the TE of 1330 mRNAs was significantly upregulated after METTL3 knockdown (Figure 5A). GSE121952 analysis exhibited that after METTL14 knockdown, the TE of 844 and 1613 mRNAs was significantly downregulated and upregulated, respectively (Figure 5B). As expected, when we overlapped these TEGs, we found that only a small part (n=163) was regulated by METTL3 and METTL14 collectively (Figure 5C). As shown in Figure 5D, KEGG pathway analysis showed that these TEGs regulated by METTL3 (M3TEGs) were significantly enriched in 23 signaling pathways, such as Lysosome, cell cycle, and p53 signaling pathway, and that these TEGs regulated by METTL14 (M14TEGs) were significantly enriched in 21 distinct signaling pathways, such as cAMP signaling pathway, Spliceosome, and Protein digestion and absorption (Figure 5E). Subsequently, GO enrichment analysis of these TEGs was analyzed (Additional file 4: Supplementary Table 7, 8), which demonstrated that only a small part of biological processes (n=5) was regulated by METTL3 and METTL14 corporately, including regulation of transcription (GO:0000122, GO:0045944, GO:0045892, and GO:0045893) and cell proliferation (GO:0008284) (Figure 5F, 5G). Together, we proved that most of M3TEGs and M14TEGs are different and involved with distinct signaling pathways and biological processes in HCC.

The TEGs regulated by METTL3 and METTL14 knockdown participate in distinct signaling pathways and biological processes. (A, B) The TEGs regulated by METTL3 and METTL14 knockdown. (C) Integrated analysis of M3TEGs and M14TEGs. (D, E) KEGG pathway analysis of M3TEGs and M14TEGs. (F) The number of biological processes regulated by M3TEGs and M14TEGs. (G) The common biological processes regulated by M3TEGs and M14TEGs.

Figure 5. The TEGs regulated by METTL3 and METTL14 knockdown participate in distinct signaling pathways and biological processes. (A, B) The TEGs regulated by METTL3 and METTL14 knockdown. (C) Integrated analysis of M3TEGs and M14TEGs. (D, E) KEGG pathway analysis of M3TEGs and M14TEGs. (F) The number of biological processes regulated by M3TEGs and M14TEGs. (G) The common biological processes regulated by M3TEGs and M14TEGs.

The signaling pathways and biological processes of m6A modified M3TEGs

To analyze these M3TEGs with m6A modification, we overlapped m6A modified transcripts with M3TEGs, which showed that 44% of M3TEGs were directly regulated by m6A modification (Figure 6A). KEGG pathway analysis showed that these m6A modified M3TEGs were significantly associated with 15 signaling pathways, including Lysosome, Cell cycle, TGF-beta signaling pathway, Chronic myeloid leukemia, Proteoglycans in cancer, Focal adhesion, Pathways in cancer, Transcription misregulation in cancer, Signaling pathways regulating pluripotency of stem cells, Progesterone-mediated oocyte maturation, Bladder cancer, p53 signaling pathway, Melanoma, Notch signaling pathway, and Colorectal cancer (Figure 6B). PPI network of these m6A modified M3TEGs revealed three functional molecular clusters (Figure 6C). KEGG pathway analysis of cluster 2 and cluster 3 exhibited five common signaling pathways, including Cell cycle, Lysosome, Bladder cancer, Pathways in cancer, and Focal adhesion. Besides, cluster 1 and cluster 3 revealed several unique pathways, including Ubiquitin mediated proteolysis, PI3K-Akt signaling pathway, Non-small cell lung cancer, and Spliceosome (Figure 6C). GO analysis showed that cluster 1 regulated protein polyubiquitination and ubiquitination, and that cluster 2 regulated cell cycle, and that cluster 3 regulated cell division and adhesion (Figure 6D). Furthermore, all top 30 hub genes among m6A modified M3TEGs were enriched in cluster 1 and cluster 2 (Figure 6E). In addition, the expression of most hub genes was positively correlated with METTL3 expression in HCC (97%) and significantly associated with the OS of HCC patients (67%) (Additional file 5: Supplementary Table 9, 10). The data described above suggested that METTL3 directly regulates part of mRNAs TE through catalyzing m6A modification, and that like m6A modified M3DEGs, these m6A modified M3TEGs are also mainly involved with protein ubiquitination and cell cycle.

Figure 6. The signaling pathways and biological processes of m6A modified M3TEGs. (A) Integrated analysis of M3TEGs and mRNAs with m6A peaks. (B) KEGG pathway analysis of m6A modified M3TEGs. (C, D) KEGG pathway and GO enrichment analysis of functional molecular clusters among PPI network of m6A modified M3TEGs. (E) The top 30 hub genes among m6A modified M3TEGs.

The signaling pathways and biological processes of m6A modified M14TEGs

To explore the m6A modified M14TEGs, we overlapped m6A modified genes with M14TEGs, which showed that only 20% of M14TEGs were m6A modified (Figure 7A). KEGG pathway analysis exhibited that these m6A modified M14TEGs were significantly enriched in 7 signaling pathways, including pathways in cancer, RNA transport, Spliceosome, Ribosome biogenesis in eukaryotes, Shigellosis, Basal cell carcinoma, and Protein digestion and absorption (Figure 7B). PPI network of these m6A modified M14TEGs revealed four functional molecular clusters (Figure 7C). KEGG pathway analysis of cluster 1 and cluster 3 showed two common signaling pathways, including Ribosome biogenesis in eukaryotes and Spliceosome. Additionally, cluster 2 and cluster 4 added two unique pathways, including Cell cycle and Ubiquitin mediated proteolysis (Figure 7C). GO enrichment analysis showed that cluster 1 participated in rRNA processing and ribosomal large subunit biogenesis, and that cluster 2 participated in cell division and microtubule-based movement, and that cluster 3 participated in mRNA splice and processing, and that cluster 4 participated in ubiquitin-dependent protein catabolic process, protein polyubiquitination and ubiquitination (Figure 7D). Moreover, all top 30 hub genes among m6A modified M14TEGs were enriched in cluster 1 and cluster 3 (Figure 7E). Meanwhile, the expression of most hub genes was positively correlated with METTL3 expression in HCC (90%) and associated with the OS of HCC patients (60%) (Additional file 6: Supplementary Table 11, 12). All in all, our data indicated that METTL14 directly regulates a small number of mRNAs TE through catalyzing m6A modification and that m6A modified M14TEGs not only are associated with ribosome biogenesis and mRNA splice but also regulate cell cycle and protein ubiquitination in collaboration with M3TEGs.

Figure 7. The signaling pathways and biological processes of m6A modified M14TEGs. (A) Integrated analysis of M14TEGs and mRNAs with m6A peaks. (B) KEGG pathway analysis of m6A modified M14TEGs. (C, D) KEGG pathway and GO enrichment analysis of functional molecular clusters among the PPI network of m6A modified M14TEGs. (E) The top 30 hub genes among m6A modified M14TEGs.

Discussion

Despite chemical modifications of DNA and histones, epigenetic regulation also contains hundreds of distinct post-transcriptional modifications in cellular RNAs among which m6A is the most abundant one occurring in eukaryotic mRNAs [30, 31]. Through its direct binding with m6A binding proteins, m6A influences almost every step in the process of an mRNA molecule, from splicing, stability, and export to translation [13, 14, 32, 33]. As core molecules of MTC, METTL3 and METTL14 often perform an opposite effect on tumor process in same cancer [19, 22]. The crystal structure analysis revealed that only METTL3, but not METTL14, carries out the catalytic role and that METTL14 mainly enhances the activity of METTL3 in part through stabilizing structure [26]. In addition, METTL3 was reported to promote the translation of important oncogenes in human cancers independent of its catalytic activity and m6A binding proteins [34]. Moreover, a cancer-associated mutation in METTL14 not only decreased methyltransferase activity but also decreased the substrate specificity of the MTC, such that both the consensus sequence GGACU and the non-consensus sequence GGAUU were methylated at similar efficiencies [4]. Hence, we hypothesized that apart from the synergistic effect of METTL3 and METTL14 on catalyzing m6A modification, each molecule has its unique functions.

In this study, we found that although METTL3 or METTL14 knockdown resulted in the abnormal expression of hundreds of mRNAs in HepG2 cells, only a small part of DEGs was affected by METTL3 and METTL14 together, and even half of the expression of DEGs had the opposite change. In agreement, the analysis of the KEGG pathway and GO enrichment showed that METTL3 and METTL14 jointly regulated few signal pathways and biological processes. The above data suggested that in HCC, most of the targets and biological functions of METTL3 and METTL14 are different. Further analysis revealed that only 37% of M3DEGs and 54% of M14DEGs have m6A peaks. KEGG pathway analysis showed that these m6A modified DEGs only participated in half of the signaling pathways regulated by METTL3 and METTL14, in which the TGF-beta signaling pathway was the same one. The TGF-beta signaling pathway is an important oncogenic pathway in cancers. Our previous study demonstrated that TGF-beta promotes the metastasis of triple-negative breast cancer [35]. It is worth noting that both hub M3DEGs and M14DEGs modified by m6A were involved in protein ubiquitination. In fact, protein ubiquitination is closely related to the TGF-beta signaling pathway. For instance, the four and a half LIM-only protein 2 (FHL2) activates TGF-beta signaling by regulating the ubiquitination of the E3 ligase [36]. Our results implied that M3DEGs and M14DEGs partake in several signaling pathways and biological processes collectively in an m6A dependent manner.

Moreover, we further analyzed the effect of METTL3 and METTL14 on mRNA translation. In accordance with transcriptome analysis, METTL3 or METTL14 downregulation resulted in hundreds of TEGs. However, only a few TEGs were jointly controlled by METTL3 and METTL14. Surprisingly, KEGG pathway analysis revealed that M3TEGs and M14TEGs did not share any signaling pathway. Consistently, GO enrichment analysis showed that even fewer biological processes were co-regulated by METTL3 and METTL14 through influencing mRNA TE. These data indicated that M3TEGs and M14TEGs participate in distinct signaling pathways and biological processes. Furthermore, we found that 44% of M3TEGs and 20% of M14TEGs were m6A modified. Interestingly, most signaling pathways of M3TEGs were the same as that of M3DEGs, such as the p53 signaling pathway and cell cycle. In addition, like M3DEGs, M3TEGs also had significant enrichment in protein ubiquitination and cell cycle. An uncontrolled cell cycle is a critical stimulus to tumor progression. For example, the integration of genomic and transcriptional features in pancreatic cancer reveals an increased cell cycle progression in metastasis [37]. Indeed, protein ubiquitination is also significantly associated with cell cycle transition [38, 39]. Hence, we hypothesized that METTL3 may regulate specific mRNAs related to protein ubiquitination and cell cycle at the transcription and translation levels dependent on m6A modification, to promote the progress of HCC. Unlike M3TEGs, M14TEGs were mainly involved in mRNA splicing, and part of M14TEGs and M14DEGs co-regulated ribosome biogenesis in eukaryotes in HCC. Intriguingly, regulation of protein ubiquitination was also an important biological process of M14TEGs. Based on our results, we speculated that m6A modification regulated by METTL3 and METTL14 might be closely related to anomalous protein ubiquitination in HCC. Besides, when analyzing the hub genes in M3DEGs, M14DEGs, M3TEGs, and M14TEGs, we found that expression of most hub genes was positively correlated with METTL3 expression or METTL14 expression and served as unfavorable predictors of HCC patients’ survival.

M6A modification is a comprehensive and context-dependent biological process [40]. Transcripts with m6A modification occur in different destinies if they bind to various m6A binding proteins. For instance, IGF2BPs stabilize, but YTHDF2 degrades m6A modified mRNAs [13, 32]. Besides, miRNAs can regulate m6A abundance by modulating METTL3 binding to mRNAs [41]. Intriguingly, METTL3 and METTL14 also modulate N6-methyladenosine-dependent primary miRNA processing [19, 42]. Moreover, miRNAs can also regulate mRNAs expression through direct binding to the 3'UTRs of targets by inducing mRNAs degradation and/or translational repression [43]. Therefore, it is limited to study the role of METTL3 and METTL14 in cancer only from the perspective of catalyzing m6A modification.

There are still some deficiencies in our study. First, to identified m6A modified mRNAs, we just overlapped dysregulated mRNAs with those having m6A peaks. Whether these dysregulated mRNAs are m6A modified ones should be validated through m6ARIP-PCR. Second, the cell line used in the GSE63591 dataset is HeLa cells. Although m6A is a conserved modification [4], similar investigations in HCC cells ought to be done to verify our results. Third, our analysis only described the role of METTL3 and METTL14 in HCC, whether this hypothesis applies to other types of cancer, such as colorectal cancer and glioma, needs further investigations. Forth, the expression of hub genes, especially these in M3TEGs and M14TEGs, should be validated through IHC analysis.

In conclusion, our study further clarified the characteristics of METTL3 and METTL14 in HCC which is beneficial to further research of METTL3 and METTL14.

Materials and Methods

Study cohort

The METTL3 and METTL14 expression data in NTs and HCC tissues were retrieved from the TCGA database (https://cancergenome.nih.gov/) and the GEO database (https://www.ncbi.nlm.nih.gov/geo/) (GSE14520 and GSE54236). TCGA database of HCC contains 50 matched NTs and HCC tissues. The GSE14520 dataset includes 214 NTs and 225 HCC tissues, while GSE54236 consists of 80 NTs and 81 HCC tissues. The RNA sequencing data of METTL3 knockdown and METTL14 knockdown in HepG2 cells were enrolled from GSE37001 and GSE90642, both METTL3 knockdown and METTL14 knockdown resulted from small interfering RNA (siRNA) transfection. The m6A peaks enriched in transcripts in HepG2 cells were analyzed and deposited in GSE37003. In addition, ribosome sequencing profiles about the effect of METTL3 knockdown and METTL14 knockdown on mRNA TE was analyzed by using GSE63591 and GSE121952. siRNA and short hairpin RNA (shRNA) was used to decrease METTL3 and METTL14 expression, respectively.

Immunohistochemistry

Immunohistochemistry analysis for METTL3 and METTL14 was performed on paraffin sections using a primary antibody against METTL3 (1:2000, Proteintech, 15073-1-AP), METTL14 (1:1500, Proteintech, 26158-1-AP), and a horseradish peroxidase conjugated IgG (1:500; Invitrogen). Three high power fields (400×magnification) were randomly selected from 30 paired HCC tissues and ANTs. For histological scoring, the degree of positivity was initially classified according to scoring both the proportion of positive staining tumor cells and the staining intensities. Scores representing the proportion of positively stained tumor cells were graded as: 0 (<10%); 1 (11%-25%); 2 (26%-50%); 3 (51%-75%) and 4 (>75%). The intensity of staining was determined as: 0 (no staining); 1 (weak staining = light yellow); 2 (moderate staining = yellow brown); and 3 (strong staining = brown). The staining index (SI) was calculated as the product of staining intensity × percentage of positive tumor cells, resulting in scores of 0, 1, 2, 3, 4, 6, 8, 9, and 12. The reactivity degree was assessed by at least two pathologists independently. Informed consent was obtained from patients and the study was approved by the ethics committee of Nanjing First Hospital.

Bioinformatics analysis

The differentially expressed METTL3 and METTL14 between HCC tissues and NTs were analyzed by using R software with the “limma” package (http://www.bioconductor.org/packages/release/bioc/html/limma.html). In TCGA database, GSE14520, and GSE54236, genes with |logFC| > 1 and p-value < 0.05 were regarded as the significantly dysregulated ones. The correlations between the two genes in HCC were identified by using the TCGA database. The associations between gene expression and prognosis of HCC patients were explored in the Kaplan-Meier plotter database (http://kmplot.com/analysis/) which is an online tool based on GEO, European Genome-phenome Archive (EGA), and TCGA databases with the median expression value as the cut-off value. The hazard ratio (HR) and log-rank p-value were calculated. Log-rank p < 0.05 was considered to be statistically significant. HR > 1 means that gene expression was negatively correlated with prognosis, while HR < 1 shows a positive correlation. Due to the different analysis platforms, we identified that in GSE90642, DEGs were those with |logFC| > 0.585 and p-value < 0.05, while in GSE37001, DEGs were those with |logFC| > 1 and p-value < 0.05. In GSE63591, TEGs were those with |logFC| > 2 and p-value < 0.05. However, in GSE121952, TEGs were those with |logFC| > 1 and p-value < 0.05. The DAVID website (https://david.ncifcrf.gov) was used to carry out GO enrichment and KEGG pathway analysis of putative targets of METTL3 and METTL14. PPI networks were constructed by using the Search Tool for the Retrieval of Interacting Genes (STRING) online tool (http://string.embl.de/) and further optimized in Cytoscape. The combined score higher than 0.40 was regarded as statistical significance. Moreover, the functional molecular complexes in the PPI network were identified automatedly by using Molecular Complex Detection (MCODE) app in Cytoscape, and the MOCODE score higher than 5.0 was regarded as statistical significance. In addition, the top 30 hub genes among PPI network were selected by the Cytohubba app with Maximal Clique Centrality (MCC) method.

Statistical analysis

Data were expressed as mean ± SD (standard deviation) and performed by using GraphPad Prism 8 (GraphPad, USA) software. The difference between groups was tested by two-tail Student’s paired or unpaired t-test. p < 0.05 was considered to be statistically significant.

Supplementary Materials

Supplementary Table 1

Supplementary Table 2

Supplementary Tables 3, 4, 5 and 6

Supplementary Table 7

Supplementary Table 8

Supplementary Tables 9, 10, 11 and 12

Abbreviations

m6A: N6-methyladenosine; MTC: methyltransferase complex; METTL3: methyltransferase-like 3; METTL14: methyltransferase-like 14; HCC: hepatocellular carcinoma; TCGA: The Cancer Genome Atlas; GEO: Gene Expression Omnibus; DEGs: differentially expressed genes; mRNAs: messenger RNAs; 3’UTR: 3’untranslated region; WTAP: Wilms’ tumor 1-associating protein; KIAA1429: Vir Like M6A Methyltransferase Associated; RBM15: RNA Binding Motif Protein 15; ZC3H13: Zinc Finger CCCH-Type Containing 13; FTO: fat mass and obesity-associated protein; ALKBH5: alkylated DNA repair protein alkB homolog 5; YTHDC1: YTH Domain-Containing Protein 1; YTHDF1/2/3: YTH Domain Family, Member 1/2/3; IGF2BP1/2/3: Insulin Like Growth Factor 2 MRNA Binding Protein 1/2/3; CRC: colorectal cancer; SOX2: SRY-Box 2; siRNA: small interfere RNA; TE: translational efficiency; shRNA: short hairpin RNA; OS: overall survival; RFS: relapse-free survival; PFS: progression-free survival; DSS: disease-specific survival; EGA: European Genome-phenome Archive; HR: hazard ratio; TEG: genes with significantly changed TE; GO: Gene ontology; KEGG: Kyoto Encyclopedia of Genes and Genomes; STRING: Search Tool for the Retrieval of Interacting Genes; MCODE: Molecular Complex Detection; MCC: Maximal Clique Centrality; M14DEGs: DEGs regulated by METTL14; M3DEGs: DEGs regulated by METTL3; M3TEGs: TEGs regulated by METTL3; M14TEGs: TEGs regulated by METTL14; FHL2: four and a half LIM-only protein 2.

Author Contributions

LXX, WSK, QJ, and GTY discussed and designed this study; LCM, ZKX, PB, XM, and XXN analyzed TCGA and GEO databases; LXX and CXX wrote the manuscript. HBS, PYQ, XT, and SHL revised the manuscript. All authors read and approved the final manuscript.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Funding

This project was supported by grants from The National Nature Science Foundation of China (No. 81972806), Jiangsu Provincial Key Research and Development Plan (BE2019614), Key Project of Science and Technology Development of Nanjing Medicine (ZDX16001) to SKW; The National Nature Science Foundation of China (No. 81802093) to HLS; Innovation team of Jiangsu provincial health-strengthening engineering by science and education (CXTDB2017008); Jiangsu Youth Medical Talents Training Project to BSH(QNRC2016066) and YQP (QNRC2016074); Jiangsu 333 High-level Talents Cultivating Project (Gastric cancer, no. BRA201702).

References

1. Desrosiers R, Friderici K, Rottman F. Identification of methylated nucleosides in messenger RNA from novikoff hepatoma cells. Proc Natl Acad Sci USA. 1974; 71:3971–75. https://doi.org/10.1073/pnas.71.10.3971 [PubMed]
2. Dominissini D, Moshitch-Moshkovitz S, Schwartz S, Salmon-Divon M, Ungar L, Osenberg S, Cesarkas K, Jacob-Hirsch J, Amariglio N, Kupiec M, Sorek R, Rechavi G. Topology of the human and mouse m6A RNA methylomes revealed by m6A-seq. Nature. 2012; 485:201–06. https://doi.org/10.1038/nature11112 [PubMed]
3. Bodi Z, Button JD, Grierson D, Fray RG. Yeast targets for mRNA methylation. Nucleic Acids Res. 2010; 38:5327–35. https://doi.org/10.1093/nar/gkq266 [PubMed]
4. Liu J, Yue Y, Han D, Wang X, Fu Y, Zhang L, Jia G, Yu M, Lu Z, Deng X, Dai Q, Chen W, He C. A METTL3-METTL14 complex mediates mammalian nuclear RNA N6-adenosine methylation. Nat Chem Biol. 2014; 10:93–95. https://doi.org/10.1038/nchembio.1432 [PubMed]
5. Ping XL, Sun BF, Wang L, Xiao W, Yang X, Wang WJ, Adhikari S, Shi Y, Lv Y, Chen YS, Zhao X, Li A, Yang Y, et al. Mammalian WTAP is a regulatory subunit of the RNA N6-methyladenosine methyltransferase. Cell Res. 2014; 24:177–89. https://doi.org/10.1038/cr.2014.3 [PubMed]
6. Lan T, Li H, Zhang D, Xu L, Liu H, Hao X, Yan X, Liao H, Chen X, Xie K, Li J, Liao M, Huang J, et al. KIAA1429 contributes to liver cancer progression through N6-methyladenosine-dependent post-transcriptional modification of GATA3. Mol Cancer. 2019; 18:186. https://doi.org/10.1186/s12943-019-1106-z [PubMed]
7. Patil DP, Chen CK, Pickering BF, Chow A, Jackson C, Guttman M, Jaffrey SR. m⁶A RNA methylation promotes XIST-mediated transcriptional repression. Nature. 2016; 537:369–73. https://doi.org/10.1038/nature19342 [PubMed]
8. Knuckles P, Lence T, Haussmann IU, Jacob D, Kreim N, Carl SH, Masiello I, Hares T, Villaseñor R, Hess D, Andrade-Navarro MA, Biggiogera M, Helm M, et al. Zc3h13/flacc is required for adenosine methylation by bridging the mRNA-binding factor Rbm15/spenito to the m⁶A machinery component wtap/Fl(2)d. Genes Dev. 2018; 32:415–29. https://doi.org/10.1101/gad.309146.117 [PubMed]
9. Mathiyalagan P, Adamiak M, Mayourian J, Sassi Y, Liang Y, Agarwal N, Jha D, Zhang S, Kohlbrenner E, Chepurko E, Chen J, Trivieri MG, Singh R, et al. FTO-dependent N⁶-methyladenosine regulates cardiac function during remodeling and repair. Circulation. 2019; 139:518–32. https://doi.org/10.1161/CIRCULATIONAHA.118.033794 [PubMed]
10. Li XC, Jin F, Wang BY, Yin XJ, Hong W, Tian FJ. The m6A demethylase ALKBH5 controls trophoblast invasion at the maternal-fetal interface by regulating the stability of CYR61 mRNA. Theranostics. 2019; 9:3853–65. https://doi.org/10.7150/thno.31868 [PubMed]
11. Xiao W, Adhikari S, Dahal U, Chen YS, Hao YJ, Sun BF, Sun HY, Li A, Ping XL, Lai WY, Wang X, Ma HL, Huang CM, et al. Nuclear m(6)A reader YTHDC1 regulates mRNA splicing. Mol Cell. 2016; 61:507–19. https://doi.org/10.1016/j.molcel.2016.01.012 [PubMed]
12. Wang X, Zhao BS, Roundtree IA, Lu Z, Han D, Ma H, Weng X, Chen K, Shi H, He C. N(6)-methyladenosine modulates messenger RNA translation efficiency. Cell. 2015; 161:1388–99. https://doi.org/10.1016/j.cell.2015.05.014 [PubMed]
13. Du H, Zhao Y, He J, Zhang Y, Xi H, Liu M, Ma J, Wu L. YTHDF2 destabilizes m(6)A-containing RNA through direct recruitment of the CCR4-NOT deadenylase complex. Nat Commun. 2016; 7:12626. https://doi.org/10.1038/ncomms12626 [PubMed]
14. Shi H, Wang X, Lu Z, Zhao BS, Ma H, Hsu PJ, Liu C, He C. YTHDF3 facilitates translation and decay of N⁶-methyladenosine-modified RNA. Cell Res. 2017; 27:315–28. https://doi.org/10.1038/cr.2017.15 [PubMed]
15. Huang H, Weng H, Sun W, Qin X, Shi H, Wu H, Zhao BS, Mesquita A, Liu C, Yuan CL, Hu YC, Hüttelmaier S, Skibbe JR, et al. Recognition of RNA N⁶-methyladenosine by IGF2BP proteins enhances mRNA stability and translation. Nat Cell Biol. 2018; 20:285–95. https://doi.org/10.1038/s41556-018-0045-z [PubMed]
16. Wang P, Doxtader KA, Nam Y. Structural basis for cooperative function of Mettl3 and Mettl14 methyltransferases. Mol Cell. 2016; 63:306–17. https://doi.org/10.1016/j.molcel.2016.05.041 [PubMed]
17. Han J, Wang JZ, Yang X, Yu H, Zhou R, Lu HC, Yuan WB, Lu JC, Zhou ZJ, Lu Q, Wei JF, Yang H. METTL3 promote tumor proliferation of bladder cancer by accelerating pri-miR221/222 maturation in m6A-dependent manner. Mol Cancer. 2019; 18:110. https://doi.org/10.1186/s12943-019-1036-9 [PubMed]
18. Wang Q, Chen C, Ding Q, Zhao Y, Wang Z, Chen J, Jiang Z, Zhang Y, Xu G, Zhang J, Zhou J, Sun B, Zou X, Wang S. METTL3-mediated m⁶A modification of HDGF mRNA promotes gastric cancer progression and has prognostic significance. Gut. 2020; 69:1193–205. https://doi.org/10.1136/gutjnl-2019-319639 [PubMed]
19. Ma JZ, Yang F, Zhou CC, Liu F, Yuan JH, Wang F, Wang TT, Xu QG, Zhou WP, Sun SH. METTL14 suppresses the metastatic potential of hepatocellular carcinoma by modulating N⁶-methyladenosine-dependent primary MicroRNA processing. Hepatology. 2017; 65:529–43. https://doi.org/10.1002/hep.28885 [PubMed]
20. Cui Q, Shi H, Ye P, Li L, Qu Q, Sun G, Sun G, Lu Z, Huang Y, Yang CG, Riggs AD, He C, Shi Y. m⁶A RNA methylation regulates the self-renewal and tumorigenesis of glioblastoma stem cells. Cell Rep. 2017; 18:2622–34. https://doi.org/10.1016/j.celrep.2017.02.059 [PubMed]
21. Visvanathan A, Patil V, Arora A, Hegde AS, Arivazhagan A, Santosh V, Somasundaram K. Essential role of METTL3-mediated m⁶A modification in glioma stem-like cells maintenance and radioresistance. Oncogene. 2018; 37:522–33. https://doi.org/10.1038/onc.2017.351 [PubMed]
22. Chen M, Wei L, Law CT, Tsang FH, Shen J, Cheng CL, Tsang LH, Ho DW, Chiu DK, Lee JM, Wong CC, Ng IO, Wong CM. RNA N6-methyladenosine methyltransferase-like 3 promotes liver cancer progression through YTHDF2-dependent posttranscriptional silencing of SOCS2. Hepatology. 2018; 67:2254–70. https://doi.org/10.1002/hep.29683 [PubMed]
23. Li T, Hu PS, Zuo Z, Lin JF, Li X, Wu QN, Chen ZH, Zeng ZL, Wang F, Zheng J, Chen D, Li B, Kang TB, et al. METTL3 facilitates tumor progression via an m⁶A-IGF2BP2-dependent mechanism in colorectal carcinoma. Mol Cancer. 2019; 18:112. https://doi.org/10.1186/s12943-019-1038-7 [PubMed]
24. Yang X, Zhang S, He C, Xue P, Zhang L, He Z, Zang L, Feng B, Sun J, Zheng M. METTL14 suppresses proliferation and metastasis of colorectal cancer by down-regulating oncogenic long non-coding RNA XIST. Mol Cancer. 2020; 19:46. https://doi.org/10.1186/s12943-020-1146-4 [PubMed]
25. Chen X, Xu M, Xu X, Zeng K, Liu X, Sun L, Pan B, He B, Pan Y, Sun H, Xia X, Wang S. METTL14 suppresses CRC progression via regulating N6-methyladenosine-dependent primary miR-375 processing. Mol Ther. 2020; 28:599–612. https://doi.org/10.1016/j.ymthe.2019.11.016 [PubMed]
26. Zhou KI, Pan T. Structures of the m(6)A methyltransferase complex: two subunits with distinct but coordinated roles. Mol Cell. 2016; 63:183–85. https://doi.org/10.1016/j.molcel.2016.07.005 [PubMed]
27. Roessler S, Jia HL, Budhu A, Forgues M, Ye QH, Lee JS, Thorgeirsson SS, Sun Z, Tang ZY, Qin LX, Wang XW. A unique metastasis gene signature enables prediction of tumor relapse in early-stage hepatocellular carcinoma patients. Cancer Res. 2010; 70:10202–12. https://doi.org/10.1158/0008-5472.CAN-10-2607 [PubMed]
28. Villa E, Critelli R, Lei B, Marzocchi G, Cammà C, Giannelli G, Pontisso P, Cabibbo G, Enea M, Colopi S, Caporali C, Pollicino T, Milosa F, et al. Neoangiogenesis-related genes are hallmarks of fast-growing hepatocellular carcinomas and worst survival. Results from a prospective study. Gut. 2016; 65:861–69. https://doi.org/10.1136/gutjnl-2014-308483 [PubMed]
29. Huang H, Weng H, Zhou K, Wu T, Zhao BS, Sun M, Chen Z, Deng X, Xiao G, Auer F, Klemm L, Wu H, Zuo Z, et al. Histone H3 trimethylation at lysine 36 guides m⁶A RNA modification co-transcriptionally. Nature. 2019; 567:414–19. https://doi.org/10.1038/s41586-019-1016-7 [PubMed]
30. Wu R, Jiang D, Wang Y, Wang X. N (6)-methyladenosine (m(6)A) methylation in mRNA with a dynamic and reversible epigenetic modification. Mol Biotechnol. 2016; 58:450–59. https://doi.org/10.1007/s12033-016-9947-9 [PubMed]
31. Fu Y, Dominissini D, Rechavi G, He C. Gene expression regulation mediated through reversible m⁶A RNA methylation. Nat Rev Genet. 2014; 15:293–306. https://doi.org/10.1038/nrg3724 [PubMed]
32. Müller S, Glaß M, Singh AK, Haase J, Bley N, Fuchs T, Lederer M, Dahl A, Huang H, Chen J, Posern G, Hüttelmaier S. IGF2BP1 promotes SRF-dependent transcription in cancer in a m6A- and miRNA-dependent manner. Nucleic Acids Res. 2019; 47:375–90. https://doi.org/10.1093/nar/gky1012 [PubMed]
33. Xu C, Wang X, Liu K, Roundtree IA, Tempel W, Li Y, Lu Z, He C, Min J. Structural basis for selective binding of m6A RNA by the YTHDC1 YTH domain. Nat Chem Biol. 2014; 10:927–29. https://doi.org/10.1038/nchembio.1654 [PubMed]
34. Lin S, Choe J, Du P, Triboulet R, Gregory RI. The m(6)A methyltransferase METTL3 promotes translation in human cancer cells. Mol Cell. 2016; 62:335–45. https://doi.org/10.1016/j.molcel.2016.03.021 [PubMed]
35. Zeng K, He B, Yang BB, Xu T, Chen X, Xu M, Liu X, Sun H, Pan Y, Wang S. The pro-metastasis effect of circANKS1B in breast cancer. Mol Cancer. 2018; 17:160. https://doi.org/10.1186/s12943-018-0914-x [PubMed]
36. Xia T, Lévy L, Levillayer F, Jia B, Li G, Neuveut C, Buendia MA, Lan K, Wei Y. The four and a half LIM-only protein 2 (FHL2) activates transforming growth factor β (TGF-β) signaling by regulating ubiquitination of the E3 ligase arkadia. J Biol Chem. 2013; 288:1785–94. https://doi.org/10.1074/jbc.M112.439760 [PubMed]
37. Connor AA, Denroche RE, Jang GH, Lemire M, Zhang A, Chan-Seng-Yue M, Wilson G, Grant RC, Merico D, Lungu I, Bartlett JM, Chadwick D, Liang SB, et al. Integration of genomic and transcriptional features in pancreatic cancer reveals increased cell cycle progression in metastases. Cancer Cell. 2019; 35:267–82.e7. https://doi.org/10.1016/j.ccell.2018.12.010 [PubMed]
38. Silmon de Monerri NC, Yakubu RR, Chen AL, Bradley PJ, Nieves E, Weiss LM, Kim K. The ubiquitin proteome of toxoplasma gondii reveals roles for protein ubiquitination in cell-cycle transitions. Cell Host Microbe. 2015; 18:621–33. https://doi.org/10.1016/j.chom.2015.10.014 [PubMed]
39. Sisinni L, Maddalena F, Condelli V, Pannone G, Simeon V, Li Bergolis V, Lopes E, Piscazzi A, Matassa DS, Mazzoccoli C, Nozza F, Lettini G, Amoroso MR, et al. TRAP1 controls cell cycle G2-M transition through the regulation of CDK1 and MAD2 expression/ubiquitination. J Pathol. 2017; 243:123–34. https://doi.org/10.1002/path.4936 [PubMed]
40. Zhou C, Molinie B, Daneshvar K, Pondick JV, Wang J, Van Wittenberghe N, Xing Y, Giallourakis CC, Mullen AC. Genome-wide maps of m6A circRNAs identify widespread and cell-type-specific methylation patterns that are distinct from mRNAs. Cell Rep. 2017; 20:2262–76. https://doi.org/10.1016/j.celrep.2017.08.027 [PubMed]
41. Chen T, Hao YJ, Zhang Y, Li MM, Wang M, Han W, Wu Y, Lv Y, Hao J, Wang L, Li A, Yang Y, Jin KX, et al. m(6)A RNA methylation is regulated by microRNAs and promotes reprogramming to pluripotency. Cell Stem Cell. 2015; 16:289–301. https://doi.org/10.1016/j.stem.2015.01.016 [PubMed]
42. Peng W, Li J, Chen R, Gu Q, Yang P, Qian W, Ji D, Wang Q, Zhang Z, Tang J, Sun Y. Upregulated METTL3 promotes metastasis of colorectal cancer via miR-1246/SPRED2/MAPK signaling pathway. J Exp Clin Cancer Res. 2019; 38:393. https://doi.org/10.1186/s13046-019-1408-4 [PubMed]
43. Fabian MR, Sonenberg N, Filipowicz W. Regulation of mRNA translation and stability by microRNAs. Annu Rev Biochem. 2010; 79:351–79. https://doi.org/10.1146/annurev-biochem-060308-103103 [PubMed]

Research Paper Volume 12, Issue 21 pp 21638—21659