Integrative analysis of m3C associated genes reveals METTL2A as a potential oncogene in breast Cancer
Journal of Translational Medicine volume 20, Article number: 476 (2022)
RNA methylation modifications, especially m6A mRNA modification, are known to be extensively involved in tumor development. However, the relationship between N3-methylcytidine (m3C) related genes and tumorigenesis has rarely been studied. In this research, we found that m3C-related genes were expressed at different levels and affected patients’ prognosis across multiple cancer types from The Cancer Genome Atlas and multi-omics levels. Importantly, methyltransferase-like proteins 2A (METTL2A) had a high amplification frequency (~ 7%) in patients with breast invasive carcinoma (BRCA), and its overexpression was an independent predictor of poor overall survival. Enrichment analysis of associated genes revealed that METTL2A may activate DNA synthesis and cell proliferation pathways in BRCA cells. Through drug sensitivity analysis, Trifluridine, PD407824, and Taselisib were shown to be effective drugs for METTL2A-positive BRCA patients. Overall, our research conducts a holistic view of the expression level and prognostic signature of m3C-related genes with multiple malignancies. Importantly, METTL2A has been intensely explored as a potential oncogene in BRCA, to aid the development of potential drug agents for precision therapy in breast cancer patients.
With decades of epitranscriptomal studies, it is increasingly clear that the covalent modifications in RNAs are a crucial additional aspect of gene regulation. It is convinced that 163 post-transcriptional modifications of RNA contribute strongly to the diversity of functions fulfilled by RNA molecules, especially during the process of tumor progression . Methylations are among the most common RNA modifications, which include N6-methyladenosine (m6A) , N1-methyladenosine (m1A) , N5-methylcytidine (m5C) , and N3-methylcytidine (m3C) . Various methylation modifications have a wide range of effects on the RNA secondary structure folding, stability, and function, which influence the physiological processes and are closely related to numerous human diseases .
The m6A modification, one of the most well-known mRNA modifications, is extensively involved in all sorts of tumor development . In contrast, the relationship between m3C modification and tumorigenesis has very rarely been studied. The 3-methylcytidine modification at position 32 (m3C32) is discovered in fission yeast with two enzymes, Trm140 and Trm141, catalyze tRNAThr and tRNASer m3C32 modification, respectively . The N6-threonylcarbamoyladenosine (t6A) or N6-isopentenyladenosine (i6A) modification at position 37 (t6A37 or i6A37) of yeast tRNAThr or tRNASer function as a key determinant for the m3C process . While the m3C modification for tRNAthr/tRNAArg and tRNAser are catalyzed by Mettl2 and Mettl6 in mice, respectively. In addition, the deletion of Mettl8 has little influence on the abundance of tRNAs which supports the hypothesis that Mettl8 is an mRNA m3C methyltransferase rather than a tRNA . In human cells, N3-methylcytidine of RNA can be catalyzed by methyltransferase-like proteins (METTLs) that contain 4 members of METTL2A, METTL2B, METTL6, and METTL8 . METTL2A and METTL2B are two homologs of Mettl2 with only six amino acids differ. With G35 and t6A at position 37 (t6A37) on the anticodon loop, METTL2A catalyzes the m3C32 in human tRNAThr, while there is little m3C32 modification activity presented by METTL2B. METTL6 interacts with seryl-tRNA synthetase (SerRS, encoded by SARS1) to mediate the biogenesis of m3C32 modification in human tRNASer . In contrast to the preceding three cytoplasmic m3C32-modifying enzymes, it had been established that METTL8 functioned in mitochondria [10, 11].
There were several studies figuring out the vital important relationship between METTL6 and the prognostic of malignant tumor patients. It was convinced that knockdown of METTL6 could significantly decrease sensitivity of lung cancer cells to cisplatin . A bioinformatics analysis revealed that METTL6 tended to amplify in luminal breast tumors with highly proliferative capacity [13, 14]. Recently, it is demonstrated that the deletion of METTL6 showed a significant inhibition on the proliferation of liver cancer cell lines and the high expression of METTL6 in malignant patients tended to indicate a poor prognostic [15, 16]. The m3C32 modification of mitochondria tRNASer and tRNAThr were catalyzed by the mitochondrial protein METTL8 and the high expression of METTL8 contributed an enhanced respiratory chain activity in pancreatic cancer . The alternative mRNA splicing research reported that one of the isoforms of METTL8, METLL8-iso1, had been investigated to be transported to mitochondrial with the assistance of its N-terminal mitochondrial targeting sequence . Aside from the classical m3C32 modification of mitochondrial tRNAs, upregulated of METTL8 by Yin Yang 1 (YY1) transcription factor mediated the m3C modification on the mRNA of AT-rich interactive domain-containing protein 1A and attenuated its translation, which induced the migration of breast cancer cell lines . In addition, the SUMOylated METTL8 induces the R-loop formation and promotes tumorigenesis in colorectal cancer [18, 19]. Nevertheless, the functions and mechanisms of m3C enzymes in cancer progression remain poorly understood.
Herein, we conducted the landscape analysis of four m3C enzymes' expressions and their relationship with prognosis across multiple cancer types. We also analyzed the DNA copy numbers and mutations of the four m3C associated genes using multi-omic data from The Cancer Genome Atlas (TCGA). Furthermore, we determined the correlation of METTL2A expression with clinical parameters of patients with breast cancer and analyzed related signaling pathways and potential therapeutic agents of METTL2A in BRCA patients. Finally, we detected METTL2A expression in tumor tissues of BRCA using immunohistochemistry. The data and analyses suggest that m3C associated genes, especially METTL2A, play a key function and are tightly linked to the development of various cancers.
Analysis of m3C associated genes in pan-cancer
The gene expression data and survival profiles of 33 cancers from the University of North Carolina TCGA genome characterization center was downloaded from the UCSC Xena (https://xenabrowser.net/) in log2(x + 1) transformed FPKM normalized count form [20, 21]. Because some cancer types are lack corresponding normal tissue in TCGA, we carried out the differential expression analysis between tumor tissue and normal tissue in 24 cancer types. Gene id was transformed to gene symbol using the annotation file downloaded from the NCBI website (https://www.ncbi.nlm.nih.gov/) . The genes with a P < 0.05 were assumed to be significantly variational genes in tumor tissue or among various subtypes.
The survival analysis of the four m3C associated genes in all TCGA tumor patients was carried out using the GEPIA2 database (http://gepia2.cancer-pku.cn) . Then the cox proportional hazard ratio and Kaplan–Meier curve were obtained from this website. The cBioPortal for Cancer Genomics (https://www.cbioportal.org/) provided a convenient tool to analyze multidimensional cancer genomics data [24, 25]. A total of 10967 TCGA patients across 32 studies (lack of READ) were included in the mutation and copy number variation analysis with the keyword "TCGA pan-cancer atlas". The OncoPrint module of the cBioPortal database displayed a concise and compact graphical summary of genomic alterations in the four m3C associated genes across 32 cancer types. The package perturbation clustering for data integration and disease subtyping (PINSPlus) was utilized to analysis the classification of cancer into various subtypes [26, 27]. The number of clusters was set to be 3 to keep the consistency across various tumors. The clustering algorithm was k-means and the perturbation method was noise for default. The subtypes with different survival profiles were presented on the K-M plot and Pheatmap plot.
Gene expression profiling analysis of m3C associated genes in BRCA
Then the breast invasive carcinoma cohort was enrolled to explore the genomic alterations and their relation to the expression of the four m3C associated genes. The expression data and clinical data of 719 BRCA positive patient cohort were also downloaded from the UCSC Xena for predictive analysis [20, 21]. The gene expression data and clinical information of the GSE3744 (including 7 normal samples and 40 BRCA tumor samples based on GPL570) [28, 29], GSE1456 (including 159 BRCA patients based on GPL96 and GPL97) [30, 31], GSE3494 (including 251 BRCA patients based on GPL96 and GPL97), GSE4922 (including 289 BRCA patients based on GPL96 and GPL97) , and GSE6532 (including 327 BRCA patients based on GPL96, GPL97 and 87 BRCA patients based on GPL570) , datasets were downloaded from the Gene Expression Omnibus (GEO) (https://www.ncbi.nlm.nih.gov/geo/) . The expression data were all normalized using the robust multi-array average method. Using the annotation file of GPL96, GPL97, and GPL570, gene ids were transformed to gene symbols. Finally, five validation BRCA cohorts were established. Samples with missing clinical features were excluded when analyzing. In addition, the UALCAN tool was utilized to conduct METTL2A protein expression analysis with data from Clinical Proteomic Tumor Analysis Consortium (CPTAC) dataset [35, 36].
Immunohistochemistry and scoring
Samples from four patients with histological diagnosis of BRCA at the First Hospital of Changsha, China, during 2020, were retrospectively collected. The TNM stage were categorized according to the NCCN Guidelines Version 1.2022 Breast Cancer. The use of patient information and tissues was sanctioned by the Ethics Committee of the First Hospital of Changsha. The detail information for the breast cancer tissues was included in the raw data for Fig. 6.
Immunohistochemistry (IHC) staining was performed with an immunohistochemical kit DAB chromogenic agent (Servicebio, G1211). Briefly, formalin-fixed, paraffin-embedded tumors were sectioned, and slides were deparaffinized using xylenes and rehydrated in a graded alcohol system. Then the tissue slides were sequentially treated with citric acid (PH6.0) antigen retrieval buffer for antigen retrieval and 3% H2O2/MeOH for blocking endogenous peroxidase activity. After being blocked with 3% bovine serum albumin (BSA), the sections were incubated with anti-METTL2A antibody and treated with secondary antibody (HRP labeled) from the corresponding species of primary antibody. This step was followed by DAB chromogenic reaction and nucleus counterstaining with hematoxylin stain solution. Finally, the stained sections were dehydrated, sealed and visualized. The staining intensities of protein were analyzed with ImageJ (https://imagej.net/Fiji).
Gene ontology enrichment analysis (GO) and gene set enrichment analysis (GSEA)
Based on the METTL2A correlated genes, GO enrichment analysis was conducted through the package “clusterProfiler” in R.4.0.0 . The top 30 of enriched pathways with adj P < 0.05 were assumed to be enriched GO pathways and presented on the histogram. Additionally, GSEA analysis was carried out in the TCGA BRCA cohort to investigate the tumor hallmarks that were more common in the METTL2A high expression subgroup compared to the METTL2A low expression subgroup (divided according to the median expression of METTL2A). The annotated gene set was selected (hallmark gene sets) as the reference gene set. Gene set permutations were set to 1000 times to identify significantly different pathways. The GSEA pathways with nominal P < 0.05 and FDR < 0.05 were identified as a significantly enriched hallmark pathway of GSEA.
Drug sensitivity analysis
The drug sensitivity data of cell lines were acquired from the PRISM Repurposing dataset (https://www.theprismlab.org/) and the Genomics of Drug Sensitivity in Cancer (GDSC) dataset (https://www.cancerrxgene.org/). In addition, the expression data of the PRISM cell lines was downloaded from the DepMap portal (https://depmap.org/portal/) while the GDSC cell lines from the GDSC dataset. Both projects provided the area under the dose–response curve (area under the curve-AUC) values as a measure of drug sensitivity, and a lower AUC value indicates increased sensitivity to treatment. All compounds with more than 20% missing data were excluded from this study. Then K-nearest neighbor (k-NN) imputation was applied to impute the missing AUC values. Finally, sensitivity data for 981 compounds over 22 BRCA cell lines from PRISM and 262 compounds over 46 BRCA cell lines from GDSC were included in this research.
All statistical analyses were performed in GraphPad Prism software and R statistical software (v4.0.0, R Core Team, R Foundation for Statistical Computing, Vienna, Austria). Before the statistical test, the normality test and the homogeneity test of variance were performed. Comparison of a continuous variable in two or multiple groups was performed using a parametric test (Student’s t-test or analysis of variance) if the variable was normally distributed or a nonparametric test (Wilcoxon rank-sum test or Kruskal–Wallis test). Correlation between two continuous variables was measured by either Pearson’s r correlation or Spearman’s rank-order correlation. The hazard ratio (HR) was estimated using a Cox regression model. When survival analysis was carried out using Kaplan–Meier methods, the patients with OS missing or OS time < 30 days had been deleted to reduce statistical bias, and the log-rank test was used to determine the statistical significance of differences. For all statistical analyses, a two-tailed P < 0.05 was considered statistically significant. All the analyses process had been uploaded to the github website in a R script version (https://github.com/Hannisal/M3C.git).
Widespread changes of m3C associated genes expression in multiple cancer types
The differential expression analysis of METTL2A, METTL2B METTL6 and METTL8 was conducted between tumor tissue and normal tissue in 24 different cancer types. The m3C associated genes showed aberrant expression in 17 cancer types for METTL2A (Fig. 1a), 15 cancer types for METTL2B (Fig. 1b), 14 cancer types for METTL6 (Fig. 1c) and 13 cancer types for METTL8 (Fig. 1d). As shown in Fig. 1, the expression levels of all m3C associated genes were higher than the corresponding normal tissues in the tumor tissue of Cervical squamous cell carcinoma and endocervical adenocarcinoma, Cholangiocarcinoma, Colon adenocarcinoma, Esophageal carcinoma, Head and Neck squamous cell carcinoma, Liver hepatocellular carcinoma (LIHC), Lung adenocarcinoma (LUAD), Lung squamous cell carcinoma and Stomach adenocarcinoma (STAD). In contrast, there was down-regulation of these four genes in Kidney Chromophobe (KICH), Kidney renal clear cell carcinoma (KIRC), and Thyroid carcinoma Compared to METTL8, the other three genes (METTL2A, METTL2B and METTL6) exhibited a more homogenized co-expression in various cancer types.
The prognostic role of m3C associated genes in multiple types of cancer
To study the prognostic role of m3C associated genes in multiple types of cancer, the Kaplan–Meier methods and Cox regression were utilized in the TCGA cohort. The Cox proportional hazard ratio and Kaplan–Meier curve assay results were consistent in terms of the prognostic trend in multiple cancers (Fig. 2, Table 1). In LIHC, METTL2A, METTL2B and METTL6 high expression was associated with poor prognosis. Contrarily, METTL2A and METTL2B low expression were related to poor prognosis in KIRC. Respectively, METTL2A high expression was related to poor prognosis in Adrenocortical carcinoma, Breast invasive carcinoma (BRCA), KICH, Pancreatic adenocarcinoma (PAAD) and Uveal Melanoma (UVM), while low expression was related to poor prognosis in Esophageal carcinoma and Rectum adenocarcinoma (READ). METTL6 high expression was associated with poor prognosis in LIHC, LUAD and STAD. METTL8 high expression was related to poor prognosis in Bladder Urothelial Carcinoma (BLCA), Kidney renal papillary cell carcinoma, Brain Lower Grade Glioma, LUAD and Thyroid carcinoma. Among the 15 cancer types, there was at least one m3C associated gene that was significant related to the prognosis. We carried out the PINSPlus cluster analysis to find out different subtypes of cancer and presented their survival profiles. The results showed that abnormal expression of m3C associated genes may be used as candidate markers for molecular subtypes of BLCA, KICH, LIHC and UVM (Fig. 3).
The variation on genetic alteration status of m3C associated genes in multiple types of cancer
We observed the genetic alteration status of m3C associated genes in different cancer types of the TCGA cohorts. As shown in Fig. 4, The top three cancer types with genetic alteration of METTL2A is BRCA(7%), Mesothelioma (6%), and Lymphoid Neoplasm Diffuse Large B-cell Lymphoma (4%)(Fig. 4a). The top three cancer types with genetic alteration of METTL2B is Uterine Corpus Endometrial Carcinoma (5%), Skin Cutaneous Melanoma (2%), and Ovarian serous cystadenocarcinoma (4%) (Fig. 4b). The frequency of METTL6 alteration is the highest in BLCA (4.5%), Uterine Corpus Endometrial Carcinoma (3%), and KIRC (2.5%) (Fig. 4c). METTL8 show lower alteration frequency of 3% in Ovarian serous cystadenocarcinoma, 3% in Cholangiocarcinoma, and 2% in Lung squamous cell carcinoma (Fig. 4d). Interestingly, the main genetic alteration type of METTL2A is amplification in BRCA (Fig. 4e). Together, these results reveal the conserved heterogeneity of m3C associated genes in multiple cancers, showing that m3C associated gene disorders could play an important role in different cancers.
METTL2A protein expression was upregulated in BRCA tissues
The significant expression change, the role in prognosis and amplification on genetic alteration status of METTL2A in BRCA made it attractive. Next to determine the protein expression level of METTL2A in BRCA tissues, we used IHC staining and CPTAC database. The IHC result showed that METTL2A protein expression was higher in tumor tissues compared to adjacent noncancerous tissues (Fig. 5a, b). Consistently, the protein level of METTL2A was significantly higher in BRCA tumor tissues compared with normal tissues in CPTAC database (Fig. 5c). Altogether, these results suggest that METTL2A might be an oncogene and could play an important role in the progression of BRCA.
Univariate and multivariate analysis of METTL2A in BRCA
The univariate and multivariate analysis were ultilized to investigated the independent prognostic value of METTL2A in the BRCA cohort (Table 2). Based the 719 BRCA patients with expression data and corresponding clinical profiles from TCGA cohort, the patients with clinical features as male, TX, NX, MX, Indeterminate, Equivocal and OS.time < 30 were removed for the cox proportional hazard regression analyses. Finally, in 681 BRCA patients, the Cox proportional hazard regression model was used to evaluate the impact of METTL2A expression and other clinical feature factors on survival. Univariate analysis showed that age (HR, 2.03; 95% CI 1.34–3.07; P = 8.00E-04), T stage (HR, 1.43; 95% CI 1.11–1.84; P = 0.0052), N stage (HR, 1.52; 95% CI 1.22–1.9; P = 2.00E-04), M stage (HR, 5.79; 95% CI 2.99–11.21; P < 0.001), and METTL2A expression (HR, 1.4; 95% CI 1–1.97; P = 0.0496) were important predictors of survival. In addition, multivariate analysis results showed that the high expression of METTL2A was an important independent predictor of poor overall survival (HR, 1.44; 95% CI 1–2.06; P = 0.0482).
The association between METTL2A expression and clinical parameters of BRCA
To further identify the role of METTL2A in BRCA, we further explored METTL2A expression and clinical features of BRCA using GEO data (Additional file 1: Fig. S1). From the GSE3744 data analysis, we observed METTL2A was significantly overexpressed in breast cancer tissues. In GEO BRCA cohorts including GSE1456 (n = 159), GSE4922 (n = 242) and GSE6532 (n = 238), the K-M plot revealed that BRCA patients with METTL2A high expression exhibited poor prognosis, which was consistent with the TCGA BRCA cohort.
Considering the bias coming with gender, only female BRCA patients in TCGA cohort (n = 711) were included for further analysis. Table 3 summarizes the correlation between METTL2A expression level and various clinical features in BRCA patients. METTL2A high expression was more prevalent in patients with advanced T stage and high grade (Fig. 6a–c). There was no significant relationship between METTL2A expression and age, N and M stage (Additional file 2: Fig. S2a–d). The expression of METTL2A was higher in HER2 positive and ER positive breast cancer patients (Fig. 6d, e). But there was no significant difference in the PR subgroup (Additional file 2: Fig. S2e–g). In addition, p53 mutation-positive patients showed higher METTL2A expression (Fig. 6f). Taken together, our results reveal that higher METTL2A expression was associated with higher T stage, higher grade, P53 mutation positive, HER2 positive or ER positive status in BRCA patients.
Identification of METTL2A-related signaling pathways and potential therapeutic agents in BRCA
Based on the expression of METTL2A, GSEA analysis was performed in the TCGA cohort. The results revealed that eight tumor hallmarks were enriched in the METTL2A high-expression subgroup, which included UNFOLDED_PROTEIN_RESPONSE, MYC_TARGETS_V1, G2M_CHECKPOINT, E2F_TARGETS, DNA_REPAIR, MYC_TARGETS_V2, MTORC1_SIGNALING AND MITOTIC_SPINDLE (Fig. 7a–h). Meanwhile, the genes with r > 0.3 and P < 0.05 in TCGA were included for GO enrichment analysis to study METTL2A related GO pathway in breast tumor tissue. The results showed that several cell cycle-related pathways were enriched, such as chromosome segregation, DNA biosynthetic process, spindle, chromosomal region, catalytic activity, acting on RNA, tubulin binding to name a few (Fig. 7i). The consistency of our results indicates that breast cancer cells with high METTL2A expression could be active in proliferation and DNA damage response pathways.
The drug sensitivity analysis was performed using GDSC and PRISM-derived drug response data. First, differential drug response analysis between METTL2A high expression (top decile) and low expression (bottom decile) groups was conducted to identify compounds with lower estimated AUC values in the METTL2A high expression group (log2FC < − 0.10 and P < 0.05) (Fig. 8a–d). The therapeutic effect of 13 agents was identified in GDSC and PRISM, respectively. It was shown that Trifluridine, PD407824 and Taselisib were potential therapeutic agents for METTL2A high expression BRCA patients. In contrast, patients with METTL2A high expression demonstrated resistance to pelitinib, LY2606368, cetrimonium, nitarsone, refametinib, ICL1100013, trametinib, ispinesib, Genentech and SCH772984. Trifluridine replaces thymine in DNA replication and directly mixes into the DNA double strands . PD407824 is a CHK1 and WEE1 inhibitor  and Taselisib is a potent PI3K inhibitor . Together with the results from METTL2A associated gene network, we speculated these drugs may be as effective treatment of BRCA patients with the over-expression of METTL2A.
M3C modification was first found in Saccharomyces cerevisiae . In S. cerevisiae, Trm140 and Trm141 catalyze m3C formation in tRNAThr and tRNASer, respectively . Four Trm140 and Trm141 homologs have been discovered in humans: METTL2A, METTL2B, METTL6, and METTL8. METTL2A and METTL2B were shown to modify m3C at position 32 of tRNAThr isoacceptors and tRNAArg(CCU) . Here, we found significant correlation between m3C associated genes and tumor malignancy in transcription and genetic alteration in a variety of tumor types. Consistent with the former researches [12, 15], the negative relationship between the expression of METTL6 and the prognostic of LUAD and LIHC had been verified in the TCGA cohort. In addition, the four m3C associated genes shared the same expression pattern among several cancer types (containing BRCA, LIHC, LUAD, PAAD, and READ), which implied that they all function as oncogenes in these cancers. As described above, METTL8 had been proved to be an oncogene in BRCA with its direct m3C modification on the mRNA of ARID1A , while METTL6 shown a tendency of amplification in malignant breast cancer patients . Interestingly, compared with METTL8 and METTL6, differential expression of METTL2A in BRCA tumor tissues and adjacent normal tissues was more significantly.
Expectedly, METTL2A high expression correlated with poorer survival outcomes of BRCA patients in the TCGA and multiple GEO cohorts. Based on the genetic alteration analysis of m3C associated genes in BRCA, it was convinced that there was a higher genetic alteration of METTL2A with the rate of amplification up to 7% which was about seven times higher than METTL6 . Therefore, it was further proved that METTL2A played a more important oncogene role in BRCA than the other m3C associated genes. Furthermore, we uncovered that the expression level of METTL2A was significantly positively related to T stage, grade, P53 mutation status, HER2, and ER positive status. Using IHC and CPTAC, we found that the protein level of METTL2A was up-regulated in BRCA tumor tissues. These results suggest that METTL2A might be an oncogene that plays an important role in BRCA progression. However, larger sample sizes and further research would be necessary for validation and expanding on our research.
The related Hallmark signaling pathway of METTL2A in BRCA was analyzed by the GSEA and GO database. The results showed that G2M_CHECKPOINT, E2F_TARGETS, DNA_REPAIR, and proliferation and DNA damage response pathways were correlated with the high expression of METTL2A. These results suggested that METTL2A may affect cell cycle and proliferation of cells to promote the progression of BRCA. In addition, Trifluridine, PD407824 and Taselisib were considered to be potential therapeutic agents for METTL2A high expression BRCA patients. Taselisib is currently in Phase III trials for postmenopausal women with estrogen receptor-positive (ER +) breast cancer  and non-small cell lung cancer (NSCLC) . Taselisib is a potent PI3K inhibitor that inhibits the activity of PIK3CA, PIK3CB and PIK3CG . Taselisib may be a potential treatment strategy for METTL2A high expression ER-positive BRCA patients. Furthermore, we identified another compound, Trifluridine, which interferes with DNA synthesis and inhibits cell proliferation [44, 45]. PD407824 is a CHK1 and WEE1 inhibitor which associated with METTL2A regulated gene network of DNA synthesis and cell proliferation in BRCA tumors. Taken together, METTL2A maybe a potential new therapeutic target for BRCA treatment. Nevertheless, future in-depth research and clinical studies are warranted to substantiate our findings.
This study expands our knowledge of m3C modification-associated genes in pan-cancer. Importantly, we found that the expression level of METTL2A might be a biomarker for the prognosis of BRCA tumors. Three drugs including Taselisib, Trifluridine, and PD407824 showed good predictive potential in BRCA-positive patients with METTL2A overexpression.
Availability of data and materials
All the datasets used in our analysis are publicly available and all web links are described in the “Methods” section. All codes and processed data were available on Github (https://github.com/Hannisal/M3C.git).
The cancer genome atlas
Differential drug response analysis
Gene expression omnibus
Robust multi-array average
Copy number variation
Clinical proteomic tumor analysis consortium
Gene ontology enrichment analysis
Gene set enrichment analysis
Boccaletto P, Machnicka MA, Purta E, Piatkowski P, Baginski B, Wirecki TK, de Crecy-Lagard V, Ross R, Limbach PA, Kotter A, et al. MODOMICS: a database of RNA modification pathways. 2017 update. Nucleic Acids Res. 2018;46:D303–7.
Huang H, Weng H, Chen J. m6A modification in coding and non-coding RNAs: roles and therapeutic implications in cancer. Cancer Cell. 2020;37:270–88.
Zhang C, Jia G. Reversible RNA Modification N1-methyladenosine (m1A) in mRNA and tRNA. Genomics Proteomics Bioinformatics. 2018;16:155–61.
Nombela P, Miguel-Lopez B, Blanco S. The role of m6A, m5C and Psi RNA modifications in cancer: novel therapeutic opportunities. Mol Cancer. 2021;20:18.
Hall RH. Isolation of 3-methyluridine and 3-methylcytidine from solubleribonucleic acid. Biochem Biophys Res Commun. 1963;12:361–4.
Arimbasseri AG, Iben J, Wei FY, Rijal K, Tomizawa K, Hafner M, Maraia RJ. Evolving specificity of tRNA 3-methyl-cytidine-32 (m3C32) modification: a subset of tRNAsSer requires N6-isopentenylation of A37. RNA. 2016;22:1400–10.
Han L, Marcus E, D’Silva S, Phizicky EM. S-cerevisiae Trm140 has two recognition modes for 3-methylcytidine modification of the anticodon loop of tRNA substrates. RNA. 2017;23:406–19.
Xu L, Liu X, Sheng N, Oo KS, Liang J, Chionh YH, Xu J, Ye F, Gao YG, Dedon PC, Fu XY. Three distinct 3-methylcytidine (m3C) methyltransferases modify tRNA and mRNA in mice and humans. J Biol Chem. 2017;292:14695–703.
Mao XL, Li ZH, Huang MH, Wang JT, Zhou JB, Li QR, Xu H, Wang XJ, Zhou XL. Mutually exclusive substrate selection strategy by human m3C RNA transferases METTL2A and METTL6. Nucleic Acids Res. 2021;49:8309–23.
Schöller E, Marks J, Marchand V, Bruckmann A, Powell CA, Reichold M, Mutti CD, Dettmer K, Feederle R, Hüttelmaier S, et al. Balancing of mitochondrial translation through METTL8-mediated m3C modification of mitochondrial tRNAs. Mol Cell. 2021. https://doi.org/10.1016/j.molcel.2021.10.018.
Huang MH, Peng GX, Mao XL, Wang JT, Zhou JB, Zhang JH, Chen M, Wang ED, Zhou XL. Molecular basis for human mitochondrial tRNA m3C modification by alternatively spliced METTL8. Nucleic Acids Res. 2022;50:4012–28.
Tan XL, Moyer AM, Fau-Fridley BL, Fridley BL, Fau-Schaid DJ, Schaid DJ, Fau-Niu N, Niu N, Fau-Batzler AJ, Batzler AJ, Fau-Jenkins GD, Jenkins GD, Fau-Abo RP, Abo RP, Fau-Li L, Li L, Fau-Cunningham JM, Cunningham JM, Fau-Sun Z, et al. Genetic variation predicting cisplatin cytotoxicity associated with overall survival in lung cancer patients receiving platinum-based chemotherapy. Clin Cancer Res. 2011. https://doi.org/10.1158/1078-0432.CCR-11-1133.
Gatza ML, Silva GO, Parker JS, Fan C, Perou CM. An integrated genomics approach identifies drivers of proliferation in luminal-subtype human breast cancer. Nat Genet. 2014;46:1051–9.
Abdel Ghafar MA-O, El-Rashidy MA, Gharib F, Al-Ashmawy GM. Impact of XRCC1 genetic variants on its tissue expression and breast cancer risk: a case-control study. Environ Mol Mutagen. 2021. https://doi.org/10.1002/em.22456.
Ignatova VV, Kaiser S, Ho JSY, Bing XY, Stolz P, Tan YX, Lee CL, Gay FPH, Lastres PR, Gerlini R, et al. METTL6 is a tRNA m3C methyltransferase that regulates pluripotency and tumor cell growth. Sci Adv. 2020;6:4551.
Abdel Ghafar MA-O, Elkhouly RA, Elnaggar MH, Mabrouk MM, Darwish SA, Younis RL, Elkholy RA-O. Utility of serum neuropilin-1 and angiopoietin-2 as markers of hepatocellular carcinoma. J Investig Med. 2021. https://doi.org/10.1136/jim-2020-001744.
Lee SA, Lee KA-O, Kim H, Cho JA-O. METTL8 mRNA methyltransferase enhances cancer cell migration via direct binding to ARID1A. Int J Mol Sci. 2021. https://doi.org/10.3390/ijms22115432.
Zhang LH, Zhang XY, Hu T, Chen XY, Li JJ, Raida M, Sun N, Luo Y, Gao X. The SUMOylated METTL8 Induces R-loop and tumorigenesis via m3C. iScience. 2020;23:100968.
Abdel Ghafar MT, Soliman NA. Metadherin (AEG-1/MTDH/LYRIC) expression: significance in malignancy and crucial role in colorectal cancer. Adv Clin Chem. 2022;106:235–80.
Goldman MJ, Craft B, Hastie M, Repecka K, McDade F, Kamath A, Banerjee A, Luo Y, Rogers D, Brooks AN, et al. Visualizing and interpreting cancer genomics data via the Xena platform. Nat Biotechnol. 2020;38:675–8.
Weinstein JN, Fau-Collisson EA, Collisson EA, Fau-Mills GB, Mills GB, Fau-Shaw KRM, Shaw KR, Fau-Ozenberger BA, Ozenberger BA, Fau-Ellrott K, Ellrott K, Fau-Shmulevich I, Shmulevich I, Fau-Sander C, Sander C, Fau-Stuart JM, Stuart JM. The Cancer genome atlas pan-cancer analysis project. Nat Genet. 2013. https://doi.org/10.1038/ng.2764.
Sayers EW, Beck J, Bolton EE, Bourexis D, Brister JR, Canese K, Comeau DC, Funk K, Kim S, Klimke W, et al. Database resources of the national center for biotechnology information. Nucleic Acids Res. 2021. https://doi.org/10.1093/nar/gkaa892.
Tang Z, Kang B, Li C, Chen T, Zhang Z. GEPIA2: an enhanced web server for large-scale expression profiling and interactive analysis. Nucleic Acids Res. 2019;47:W556–60.
Gao JJ, Aksoy BA, Dogrusoz U, Dresdner G, Gross B, Sumer SO, Sun YC, Jacobsen A, Sinha R, Larsson E, et al. Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Sci Signal. 2013. https://doi.org/10.1126/scisignal.2004088.
Cerami E, Gao J, Dogrusoz U, Gross BE, Sumer SO, Aksoy BA, Jacobsen A, Byrne CJ, Heuer ML, Larsson E, et al. The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. Cancer Discov. 2012;2:401–4.
Nguyen H, Shrestha S, Draghici S, Nguyen T. PINSPlus: a tool for tumor subtype discovery in integrated genomic data. Bioinformatics. 2019;35:2843–6.
Nguyen H, Tran D, Tran B, Roy M, Cassell A, Dascalu S, Draghici S, Nguyen T. SMRT: randomized data transformation for cancer subtyping and big data analysis. Front Oncol. 2021;11:725133.
Alimonti A, Carracedo A, Clohessy JG, Trotman LC, Nardella C, Egia A, Salmena L, Sampieri K, Haveman WJ, Brogi E, et al. Subtle variations in Pten dose determine cancer susceptibility. Nat Genet. 2010;42:454–8.
Richardson AL, Wang ZC, De Nicolo A, Lu X, Brown M, Miron A, Liao X, Iglehart JD, Livingston DM, Ganesan S. X chromosomal abnormalities in basal-like human breast cancer. Cancer Cell. 2006;9:121–32.
Hall P, Ploner A, Bjohle J, Huang F, Lin CY, Liu ET, Miller LD, Nordgren H, Pawitan Y, Shaw P, et al. Hormone-replacement therapy influences gene expression profiles and is associated with breast-cancer prognosis: a cohort study. BMC Med. 2006;4:16.
Pawitan Y, Bjohle J, Amler L, Borg AL, Egyhazi S, Hall P, Han X, Holmberg L, Huang F, Klaar S, et al. Gene expression profiling spares early breast cancer patients from adjuvant therapy: derived and validated in two population-based cohorts. Breast Cancer Res. 2005;7:R953-964.
Ivshina AV, George J, Senko O, Mow B, Putti TC, Smeds J, Lindahl T, Pawitan Y, Hall P, Nordgren H, et al. Genetic reclassification of histologic grade delineates new clinical subtypes of breast cancer. Cancer Res. 2006;66:10292–301.
Loi S, Haibe-Kains B, Fau-Desmedt C, Desmedt C, Fau-Lallemand F, Lallemand F, Fau-Tutt AM, Tutt AM, Fau-Gillet C, Gillet C, Fau-Ellis P, Ellis P, Fau-Harris A, Harris A, Fau-Bergh J, Bergh J, Fau-Foekens JA, Foekens JA, Fau-Klijn JGM, et al. Definition of clinically distinct molecular subtypes in estrogen receptor-positive breast carcinomas through genomic grade. J Clin Oncol. 2007. https://doi.org/10.1200/JCO.2006.07.1522.
Barrett T, Wilhite SE, Ledoux P, Evangelista C, Kim IF, Tomashevsky M, Marshall KA, Phillippy KH, Sherman PM, Holko M, et al. NCBI GEO: archive for functional genomics data sets-update. Nucleic Acids Res. 2013;41:D991–5.
Chandrashekar DS, Bashel B, Balasubramanya SAH, Creighton CJ, Ponce-Rodriguez I, Chakravarthi B, Varambally S. UALCAN: a portal for facilitating tumor subgroup gene expression and survival analyses. Neoplasia. 2017;19:649–58.
Chen FJ, Chandrashekar DS, Varambally S, Creighton CJ. Pan-cancer molecular subtypes revealed by mass-spectrometry-based proteomic characterization of more than 500 human cancers. Nat Commun. 2019;10:5679.
Yu GC, Wang LG, Han YY, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. Om J Integr Biol. 2012;16:284–7.
Bordonaro R, Calvo A, Auriemma A, Hollebecque A, Rubovszky G, Saunders MP, Pápai Z, Prager G, Stein A, André T, et al. Trifluridine/tipiracil in combination with oxaliplatin and either bevacizumab or nivolumab in metastatic colorectal cancer: a dose–expansion, phase I study. ESMO Open. 2021. https://doi.org/10.1016/j.esmoop.2021.100270.
Arora S, Bisanz KM, Peralta LA, Basu GD, Choudhary A, Tibes R, Azorsa DO. RNAi screening of the kinome identifies modulators of cisplatin response in ovarian cancer cells. Gynecol Oncol. 2010;118:220–7.
Song KW, Edgar KA, Hanan EA-O, Hafner MA-O, Oeh J, Merchant M, Sampath D, Nannini MA, Hong R, Phu L, et al. RTK-dependent inducible degradation of mutant PI3Kα drives GDC-0077 (inavolisib) efficacy. Cancer Discov. 2022. https://doi.org/10.1158/2159-8290.CD-21-0072.
Dent S, Cortes J, Im YH, Dieras V, Harbeck N, Krop IE, Wilson TR, Cui N, Schimmoller F, Hsu JY, et al. Phase III randomized study of taselisib or placebo with fulvestrant in estrogen receptor-positive, PIK3CA-mutant, HER2-negative, advanced breast cancer: the SANDPIPER trial. Ann Oncol. 2021;32:197–207.
Langer CJ, Redman MW, Wade JL 3rd, Aggarwal C, Bradley JD, Crawford J, Stella PJ, Knapp MH, Miao J, Minichiello K, et al. SWOG S1400B (NCT02785913), a phase II study of GDC-0032 (Taselisib) for previously treated PI3K-positive patients with stage IV squamous cell lung cancer (Lung-MAP Sub-Study). J Thorac Oncol. 2019;14:1839–46.
Ding X, Faber K, Shi Y, McKnight J, Dorshorst D, Ware JA, Dean B. Validation and determination of taselisib, a beta-sparing phosphoinositide 3-kinase (PI3K) inhibitor, in human plasma by LC–MS/MS. J Pharm Biomed Anal. 2016;126:117–23.
Kawazoe A, Ando T, Hosaka H, Fujita J, Koeda K, Nishikawa K, Amagai K, Fujitani K, Ogata K, Watanabe K, et al. Safety and activity of trifluridine/tipiracil and ramucirumab in previously treated advanced gastric cancer: an open-label, single-arm, phase 2 trial. Lancet Gastroenterol Hepatol. 2021;6:209–17.
Habib EM, Nosiar NA, Eid MA, Taha AM, Sherief DE, Hassan AE, Abdel Ghafar MT. MiR-150 expression in chronic myeloid leukemia: relation to imatinib response. Lab Med. 2022;53:58–64.
This study was supported by grants from the National Natural Science Foundation of China (81972773 to Z.L.), the Scientific Research Project of Hunan Provincial Health Commission (20200216 to J.L.), the Shenzhen Science and Technology Innovation Commission (JCYJ20190809103203711 to P.X.), the Changsha Municipal Natural Science Foundation (kq2014211 to L.F.) and the Postgraduate Scientific Research Innovation Project of Hunan Province (No. CX20200259 to S.W.).
Ethics approval and consent to participate
The studies involving human participants were reviewed and approved by the medical ethics committee of the First Hospital of Changsha, China. The patients/participants provided their written informed consent to participate in this study.
The authors declare no competing financial interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
ShuaiWang is the first author.
Fig. S1 Aberrant expression of METTL2A in BRCA from GEO cohorts. a, b The expression alteration of METTL2A probes between tumor (red) and normal (blue) tissue samples in BRCA from GSE3744. c The Overall survival curve of METTL2A in GSE1456. d The Disease-Free survival curve of METTL2A in GSE4922. e, f The Relapse-Free survival curves of METTL2A in e GSE1456 and f GSE6532.
Fig. S2 Relationship between METTL2A expression and clinicopathological parameters in BRCA form GEO cohorts. Distribution of METTL2A expression stratified by a, b Age, c N stage, d M stage, and e-g PR status.
About this article
Cite this article
Wang, S., Li, H., Liu, J. et al. Integrative analysis of m3C associated genes reveals METTL2A as a potential oncogene in breast Cancer. J Transl Med 20, 476 (2022). https://doi.org/10.1186/s12967-022-03683-2