Skip to main content

Increased expression of YTHDF1 and HNRNPA2B1 as potent biomarkers for melanoma: a systematic analysis



The incidence and mortality of melanoma is increasing around the world. To deeply explain the mechanism insight into it, we conducted a systematic analysis to examine the levels of regulatory genes of the common RNA epigenetic modification-N6-methyladenosine (m6A) in patients with melanoma compared by the healthy.


We analyzed the expression of m6A Eraser, Writer, and Reader genes based on publicly available datasets on Oncomine and validated the results with a gene expression omnibus dataset. Hub genes were identified with Cytohubba and the frequency of copy number alterations was analyzed with the cBioPortal tool.


The results revealed the up-regulation of YTHDF1 and HNRNPA2B1 in melanoma. Combining the two genes improved the efficacy in diagnosing melanoma by about 10% compared to each gene alone. Hub genes identified with four analysis methods were compared and the overlapping genes were selected. These genes were enriched in several gene ontology terms. Genes related to p53-signaling consisted of CDK2, CDK1, RRM2, CCNB1, and CHEK1. All five genes were positively correlated with either YTHDF1 or HNRNPA2B1, suggesting that both genes may affect m6A modification by the five genes, further up-regulating their expression and facilitate their roles in inhibiting p53 to suppress tumorigenesis. We also observed major mutations in YTHDF1 and HNRNPA2B1 that led to their amplification in melanoma. Significant differences were observed in the clinical characteristics of patients with altered and unaltered m6A regulatory genes such as tumor stage and treatment response.


We, for the first time, identified a combination of m6A regulatory genes to diagnose melanoma. We also analyzed m6A-related genes more comprehensively based on systematic complete data. We found that YTHDF1 and HNRNPA2B1 were altered in melanoma and might influence the development of the disease through signaling pathways such as p53.


Melanoma is one of the fastest developing malignancies with strong aggressive ability, however no proper curative treatments exist at present [1, 2]. It is a type of epithelial malignant tumor originating from melanocytes with obviously increasing incidence and mortality around the world [1]. Based on insight into the molecular mechanism underlying the disease, serum lactate dehydrogenase is considered as a biomarker to predict the development of melanoma [3].

With the advancement of sequencing technology, related mutations such as B-Raf proto-oncogene (BRAF) and signaling pathways such as the mitogen-activated protein kinase (MAPK) pathway have been found. These discoveries have led to the emergence of targeted drug treatment for melanoma such as inhibitors of BRAF, mitogen-activated protein kinase kinase 7 (MAP2K7, also known as MEK), and mitogen-activated protein kinase 1 (MAPK1, also known as ERK) [2, 4, 5]. However, the actual pathological mechanism in melanoma remains unknown.

N6-methyladenosine (m6A) is a type of RNA epigenetic modification [6]. Originally reported only in mRNA, experts have found m6A in other types of RNA and involved in various bioprocesses with the development of high-throughput sequencing technology [7]. Similar to DNA methylation modification, m6A is also regulated by methyltransferase and demethylase [7, 8]. Previous reports showed that m6A methylation mainly interacted with three classes of proteins: firstly, m6A methyltransferase whose coding genes were called Writers and included methyltransferase-like 3 (METTL3), methyltransferase-like 14 (METTL14), and WT1 associated protein (WTAP); secondly, m6A demethylase whose coding genes were called Erasers and included alpha-ketoglutarate dependent dioxygenase (FTO) and alkB homolog 5 (ALKBH5); thirdly, proteins that could attach to the m6A methylation site in RNA and further played a specific role in physio-pathological processes with coding genes called Readers, which included YTH N6-methyladenosine RNA-binding protein 1 (YTHDF1), YTH N6-methyladenosine RNA-binding protein 2 (YTHDF2), E74-like ETS transcription factor 3 (ELF3), heterogeneous nuclear ribonucleoprotein C (HNRNPC), and heterogeneous nuclear ribonucleoprotein A2/B1 (HNRNPA2B1) [6, 9]. Recent studies have shown that m6A methylation is closely related to tumorigenesis and tumor development [10]. Vu et al. [11] reported that CD34+ hematopoietic stem cells with lower expression of METTL3 had higher levels of phosphorylated protein kinase B and that METTL3 promoted normal stem cells into acute myeloid leukemia cells. Zhou et al. [12] reported increased expression of FTO in tissue lesions of cervical squamous cell carcinoma (CSCC) and observed that patients with high FTO expression exhibited chemotherapy tolerance. However, few studies have investigated the role of m6A modification in the development or pathological process of melanoma.

In this research, our goal was to examine the expression of m6A regulatory genes in melanoma based on open biological data and to identify the related pathways. Firstly, we used the online Oncomine tool to determine whether the m6A-related genes presented above were altered in melanoma. We then extracted the RNA-seq data from the gene expression omnibus (GEO) website ( to validate the results from Oncomine and investigate the altered genes. Next, we used statistical methods to evaluate the hub genes involved in RNA processing and calculate the correlation between the hub genes and the selected target genes. Finally, we analyzed all the critical genes and identified the signaling pathways and mutation hot spots for m6A regulatory genes to uncover new therapy targets and determine the m6A-related biological mechanism in melanoma. A summary of this study is shown in Fig. 1.

Fig. 1

Flow chart of study. M, melanoma; N, normal; ROC, receiver operating characteristic curve


Ethics statement

All clinical data, copy number alterations (CNAs), mutations, RNA-seq data, and other data were abstracted from Oncomine, cBioPortal (, or the GEO platform, which are open to the public. All analyses were performed based on the data of previously published studies. Therefore, no ethical approval or patient consent was required for this study.

Overview of the expression of 10 genes in Oncomine

The Oncomine ( tool was used to determine the expression of m6A regulators in melanoma, five studies (Haqq et al. 2005; Talantov et al. 2005; TCGA 2003; Riker et al. 2008; Critchley et al. 2006) that had mRNA or DNA expression documents were selected. There were total 615 patients with melanoma included in this research, and their information was presented in Additional file 1: Table S1. We used “Differential analysis” online tool in Oncomine portal to systematically compare the DNA or transcript levels of the normal and patients with melanoma. The whole data in the five studies above were corrected and then median ranked analysis was performed, results were shown in Fig. 2.

Fig. 2

Overview of the expression of 10 m6A-related genes in melanoma from Oncomine. Comparison of over-expression (a) and under-expression (b) of different genes among the selected studies. 1 represents the study by Critchley-Thorne et al., 2007; 2 represents Haqq et al., 2005; 3 represents Riker et al., 2008; 4 represents Talantov et al., 2005; 5 represents TCGA, 2013; the rank for a gene is the median rank for the gene across each of the analysis methods and the P-value is for the median-ranked analysis

Validation of the expression of all genes by RNA-seq data from GEO

We obtained the RNA-seq data from GEO (23 normal vs 57 patients with melanoma). Patient information is shown in Additional file 1: Table S2. R package affyPLM and RColorBrewer were used to do the quality control. Robust Multi-Array Average method was used in data pre-processing, then the files from the normal and patients’ samples were merged. Package impute in R was used to normalize and correct the files and package limma was to determine the different expression files. Heat map and volcano graphs were drew by gplots package in R software. The expression files of target genes were abstracted to do further analysis.

Functional enrichment analysis of the normalized RNA expression files from GEO

Gene ontology (GO) term and Kyoto Encyclopedia of Genes and Genomes analysis were performed by ClueGO [13, 14]. Biological Process, Cellular Component, Molecular Function, Immune System Process ontologies and KEGG pathways were selected. GO tree interval was from 3 to 8, the minimal number of genes in GO term or pathway was 3. GO term/pathway network connectivity (Kappa Score) was 0.4, the statistical option was Enrichment/Depletion (Two-sided hypergeometric test) and pV correction was Bonferroni step down. The String tool ( was used to draw the general network. Cytoscape (CytoHubba) was used to predict and explore the critical genes, nodes’ score was calculated in this software and the scores of maximal clique centrality (MCC), maximum neighborhood component (MNC), edge percolated component (EPC), and degree methods of each genes were obtained to determine the key genes in this study.

Analysis of mutation diagram for each gene with cBioPortal

We used cBioPortal to analyze the frequency of mutations and CNAs in patients with melanoma. A total of 653 samples were included in this research. Detailed information for the patients is shown in Table 1. Oncoprint tool in cBioPortal was used to scan mutation frequency of each gene in every sample, samples with mutation were presented in Fig. 5. There were three studies (Snyder, MEJM; TCGA; Van Allen, Science) included in this research. Cancer Type tool was used to show the mutation types of every gene or total genes in each study. Mutations tool in cBioPortal was used to showed the mutation sites for each gene in melanoma, and the protein post-translational modification (PTM) site plugin in mutations tool was used to predict the occurrence of phosphorylation, acetylation, ubiquitination, methylation, and O-linked glycosylation, primarily between 0 and 370aa.

Table 1 Clinical information for altered and unaltered groups on cBioPortal


For the comparison of two different groups, the t test or Mann-Whitney test was used according to the data distribution. The Chi squared test was performed to compare the discrete data and the Kruskal–Wallis test was used to compare the expression of more than two groups that did not conform to the conditions of the parametric test. The relationship between two continuous variables was evaluated with the Pearson correlation coefficient. Receiver operating characteristic (ROC) curves and scatter diagrams were drawn using GraphPad Prism 6.0. IBM SPSS Statistics version 21.0 was used to calculate all statistical parameters. The significance level was P < 0.05.


Systematic comparison of m6A regulatory genes in melanoma and normal populations

We used the Oncomine database to compare the expression of m6A-related genes in patients with melanoma and healthy individuals reported in five studies. The inclusion criteria were as follows: patients involved in the studies had a definite diagnosis or histopathological diagnosis of melanoma; the samples were from patients confirmed to have clinical data; the data were collected from tissue samples not cell lines; gene expression profiles were available. We excluded the datasets if : they were without clinical information; the samples’ type was cell line or primary cell culture; their gene expression files were not provided on the Oncomine portal; they did not have samples of both patients with melanoma and the healthy. Samples were included in this study regard less of the age, race or sex. As shown in Fig. 2, YTHDF1 and HNRNPA2B1 were significantly over-expressed in patients with melanoma while ELF3 was down-regulated. The detailed P-values are shown in Table 2. Other genes such as HNRNPC were also altered to some degree. We further validated our findings using normalized RNA-seq data from GEO.

Table 2 P-values of different genes in five studies selected in Oncomine

Validation of the expression of target genes

Original RNA-seq files for 53 patients with melanoma and 27 healthy individuals were obtained from GEO. The gene expression for all participants is shown in a cluster heat map in Fig. 3a and the up-regulated (red) and down-regulated (green) genes are visually presented in volcano plots (Fig. 3b). We examined the 10 genes of interest in this dataset and found that only HNRNPA2B1, YTHDF1, ELF3, and YTHDF2 were significantly altered in patients with melanoma (Table 3, Fig. 3c, Additional file 2: Figure S1). Integrating the results from GEO and Oncomine, the expression of HNRNPA2B1, YTHDF1, and ELF3 showed the same patterns and were thus persuasive. We assessed the specificity and sensitivity of these three genes and observed that HNRNPA2B1 and YTHDF1 had higher ROC areas than YTHDF2, whereas ELF3 had a lower value (Table 3). Therefore, we selected HNRNPA2B1 and YTHDF1 for further research and combined the expression of the two genes to evaluate the combined specificity and sensitivity. The equation for the combination used in this study was X = logit (P) = ln (P/1 − P) = − 7.777 + 3.93*YTHDF1 + 1.024*HNRNPA2B1. Next, the predicted probability of melanoma was calculated with the following equation: predicted probability (P) = ex/(1 + ex) [15, 16]. The results indicated the efficiency of combining YTHDF1 and HNRNPA2B1 compared to the results for each gene individually, with an improvement of almost 10% (Fig. 3c, Table 3). The expression and ROC curves for ELF3, and YTHDF2 are presented in Additional file 2: Figure S1.

Fig. 3

Expression of YTHDF1 and HNRNPA2B1 in the complete RNA-seq dataset from GEO. a Cluster heat map and b volcano plot of all genes in healthy individuals (n = 23) and patients with melanoma (n = 57); b red points represent up-regulated genes, green represents down-regulated genes. c Expression level and ROC curves of YTHDF1 and HNRNPA2B1 and their combination. FC, fold change; FDR, false discovery rate; AUC, area under curve; ROC, receiver operating characteristic curve. ***P < 0.001, ****P < 0.0001

Table 3 The ROC curve-related parameters and results of Mann–Whitney test

m6A-associated mechanism in melanoma

We extracted the GO terms for all genes using ClueGO and selected 611 genes involved in RNA processing or the epigenetic activation processes (Fig. 4a). After filtering the unconnected genes from the network using the String tool (Fig. 4b), we chose 586 nodes to determine the important genes using CytoHubba in Cytoscape. The top 50 genes according to the MCC, MNC, EPC and degree were extracted and a Venn diagram was created using Venny tool online. Thirty-one genes overlapped in the results from the four types of analysis (Fig. 4c). The network and MCC score heat nodes are illustrated in Fig. 4d, which shows that cyclin dependent kinase 2 (CDK2), CDK1, cyclin-B1 (CCNB1), and checkpoint kinase 1 (CHEK1) were altered. The detailed ranks for MCC and the other analysis methods are shown in Table 4. Further we performed the correlation analysis of all 31 genes and YTHDF1 or HNRNPA2B1 by calculating their Pearson correlation coefficients (Table 5). Only 17 genes were ether positively or negatively correlated with YTHDF1 or HNRNPA2B1 (P < 0.05) (Table 5).

Fig. 4

m6A-associated mechanism in melanoma. a Top 10 GO terms for RNA processing. b Network of all genes identified with the a String tool. c Venn diagram of genes identified with different analysis methods to determine the key genes in melanoma. d Network and rank scores of 31 hub genes found in the results from all four methods in (c). e Heat map of genes that had a significant relationship with YTHDF1 or HNRNPA2B1. f Scatter plots of expression of YTHDF1 and HNRNPA2B1 and all genes in (e) and their associations. g Pathway or GO terms of the critical genes in melanoma and their relationship with YTHDF1 or HNRNPA2B1. MCC, maximal clique centrality; MNC, maximum neighborhood component; EPC, edge percolated component. * P < 0.05, δ P < 0.01, Ф P < 0.001, Δ P < 0.0001

Table 4 Evaluation of hub genes
Table 5 Pearson correlation coefficient for HNRNPA2B1, YTHDF1, and hub genes. Co-gene, genes included in the correlation analysis

We visualized the expression of each gene by creating a heat map (Fig. 4e) and found that all nodes were up-regulated. Figure 4f shows that the expression of HNRNPA2B1 and YTHDF1 was positively correlated (P < 0.05). CCNB1, CDK1, CHEK1, ribonucleotide reductase regulatory subunit M2 (RRM2), aurora kinase B (AURKB), cyclin A2 (CCNA2), serine/threonine-protein kinase (PLK1), dual specificity protein kinase (TTK), cell division control protein 45 (CDC45), enhancer of zeste homolog 2 (EZH2), proliferating cell nuclear antigen (PCNA), and DNA topoisomerase 2-alpha (TOP2A) were positively associated with the level of HNRNPA2B1 (P < 0.05) (Fig. 4f). In addition to CD86 and STAT1, which were negatively correlated with YTHDF1, other genes in Fig. 4f exhibited a positive relation with YTHDF1 (P < 0.05). The interactions of m6A-regulated genes measured with the cBioPortal tool are shown in Additional file 1: Table S3 and confirm the co-occurrence of YTHDF1 and HNRNPA2B1 (P = 0.001, Q = 0.013).

Next, we created a network map to visualize the connections between these genes and further explored the pathways behind them (Fig. 4g). YTHDF1 and HNRNPA2B1 interacted with targeted genes involved in bioprocesses such as the p53 signaling pathway and positive regulation of DNA replication, which might provide insight into the molecular mechanism underlying melanoma (Fig. 4g).

Mutation and CNAs of m6A-related genes in patients with melanoma

We analyzed the data for melanoma patients with CNA and mutation profiles on cBioPortal. All mutated samples are shown in Fig. 5a. Integration of the mutation profiles with the clinical information of the patients revealed statistical differences between patients with alterations and those without alterations in terms of the tumor stage, durable clinical benefits, neo-antigen load, serum lactate dehydrogenase, and treatment response (Table 1). The mutation diagram shows the top three altered genes were YTHDF1 (8%), ELF3 (7%), and HNRNPA2B1 (6%), with a high frequency of amplification.

Fig. 5

Copy number alteration (CNA) frequency for m6A-associated genes in melanoma. a Definition of the CNAs and their percentages for m6A regulators in melanoma. b The altered features in different studies. c Mutation diagram and PTM sites of YTHDF1 and HNRNPA2B1. PTM, protein post-translational modification; Glyco, glycosylation

Three studies were included in our systematic analysis. The amplification of genes was dominant in Snyder’s study, while mutation and amplification were the major modifications in the other two studies (Fig. 5b). For YTHDF1 and HNRNPA2B1, the predominant alteration was amplification, followed by mutation (Fig. 5b). Other genes are shown in Additional file 2: Figure S2a.

Missense mutation was the most common alteration in YTHDF1, covering most of the gene (Fig. 5c). Missense and truncating mutations were observed in HNRNPA2B1. The PTM tool showed post-transcriptional modifications such as phosphorylation, acetylation, ubiquitination, methylation, malonylation, and sumoylation in the gene (Fig. 5c). The mutation frequency and PTM sites for YTHDF2, FTO, WTAP, ALKBH5, METTL14, HNRNPC, ELF3, and METTL3 are shown in Additional file 2: Figure S2b.


Melanoma is a tumor that consists of malignant melanocytes, which are pigment-producing cells of neuroectodermal origin that can be found throughout the body (including in the iris, rectum, and skin) [2]. With the advancement of science and technology coupled with increased recognition of the disease, experts’ interpretation of the pathogenesis of melanoma has changed from the previous explanation of simple molecular interactions to epigenetic regulation [17, 18]. As one of the most common chemical modifications in messenger RNA (mRNA), m6A has received substantial attention in epigenetics recently [10, 18]. Last year, Yang et al. [19] found that FTO, an m6A demethylase, could promote tumorigenesis and anti-PD-1 resistance in melanoma, highlighting the role of FTO in immunotherapy for the disease. In the same year, Jia et al. [18] observed a decrease in m6A levels in ocular melanoma samples, which could indicate poor prognosis and predict tumor progression. However, few studies based on systematic data have been conducted to explore m6A-related genes in melanoma and their general mechanisms. In this study, we summarized the expression of 10 Writer, Reader, and Eraser genes of m6A for the first time using Oncomine with more than 500 melanoma samples. We found that YTHDF1, ELF3, HNRNPA2B1, YTHDF2, FTO, and HNRNPC were significantly altered (Fig. 2, Table 2). By integrating the results from Oncomine and GEO, the up-regulation of HNRNPA2B1 and YTHDF1 was identified (Fig. 3, Table 3). We further examined whether the combination of YTHDF1 and HNRNPA2B1 had higher efficacy in diagnosing melanoma by using linear regression and found the efficacy was improved from about 0.760 to 0.857 (P < 0.0001). This combination provides the first optimized method to distinguish patients with and without melanoma based on m6A regulatory genes rather than common molecular biomarkers such as serum lactate dehydrogenase.

We also analyzed the GO terms for all genes. For the terms associated with RNA processing, 518 genes were enriched in RNA metabolic processes, with positive regulation of RNA metabolism as the secondary term. The genes for both processes were up-regulated, indicating the occurrence of dynamic and active RNA metabolic processes in melanoma (Fig. 4a). YTHDF1, an m6A Reader, was found to augment EIF3C translation by binding to m6A-modified EIF3C mRNA and to further promote tumorigenesis and metastasis in ovarian cancer [20]. Han et al. [21] demonstrated that targeting YTHDF1 with intestinal stem cells in established tumors could result in tumor shrinkage and longer survival. In our research, we found genes involved in the p53 signaling pathway such as CCNB1, CDK1, CHEK1, RRM2, and CDK2 were positively correlated with either YTHDF1 or HNRNPA2B1 (Fig. 4f). This suggested that YTHDF1 or HNRNPA2B1 may interact with the related genes above and further influence the p53 signaling pathway, resulting in the development of melanoma (Fig. 4g). CCNB1 silencing activates the p53 signaling pathway, further inhibiting the proliferation of cells and promoting cell senescence in pancreatic cancer [22]. CDK1/2 inhibitors activate p53 in tumors, especially the inhibitors of CDK2, which maintains the balance of S-phase regulatory proteins and thus coordinates subsequent p53-independent G2/M checkpoint activation [23, 24]. Hala et al. confirmed the existence of CHEK1/p53 association in human colorectal cancer in vivo and demonstrated that tumors lacking p53 had higher levels of CHEK1 [25]. The RRM1/RRM2B enzyme is capable of retaining activity in hypoxia and the replication pressure it induces is one of the factors proposed to promote selection stress that results in the loss of key ingredients for the DNA damage response, including p53 [26,27,28,29]. Therefore, YTHDF or HNRNPA2B1 might up-regulate CCNB1, CDK1, CHEK1, RRM2, or CDK2 (Fig. 4f), all of which could inhibit the role of p53 in suppressing tumors and promote the development of melanoma.

GO terms were found related to the activation of the cell cycle, DNA integrity checkpoints, histone phosphorylation, DNA damage response, signal transduction by p53 class mediation resulting in cell cycle arrest, regulation of transcription involved in the G1/S transition in the mitotic cell cycle, etc. (Fig. 4g). These findings suggested that the gene-interaction chain pathway might affect the proliferation as well as the normal physiological and metabolic processes of melanocytes, eventually resulting in the development of melanoma. All the genes in the identified pathways were also positively related to YTHDF1 or HNRNPA2B1 (Fig. 4f), suggesting the expression of both genes might influence m6A modifications in mRNA genes and potentially result in melanoma heterogeneity at the cellular level. As shown in Fig. 4f–g, STAT1 and CD86 were negatively regulated by YTHDF1. Previous studies demonstrated that the loss of STAT1 activation or expression was found in malignant cells derived from various tumors, while the loss of CD86 expression resulted in a decrease in tumor-infiltrating T lymphocytes in diffuse B-cell large-cell lymphoma [30,31,32,33,34,35]. YTHDF1 might negatively regulate STAT1 and further facilitated melanoma development, but the roles of STAT1 and CD86 in the general network were not large enough to mitigate the trend of disease development. These genes together promoted the formation of melanoma, as shown in Fig. 4g. To fully explore the role of YTHDF1 or HNRNPA2B1, we assessed their mutations using the cBioPortal Tool (Fig. 5) and found the two genes were mainly amplified in melanoma. Major mutations occurred in the non-coding domain and affected the post-transcription of related proteins by influencing modifications to the proteins such as phosphorylation (Fig. 5c). In addition, the alterations of m6A-related genes influenced the patients’ clinical features such as the tumor size, tumor stage, or treatment response (Table 3), indicating their important roles in the progression of melanoma.

However, this research also had some limitations. Firstly, the heterogeneity among the five studies included in Oncomine portal could not be ignored. The heterogeneity might result from patients’ races, treatment or other factors that could have an influence on the expression levels of patients’ transcript files. Secondly, there were few studies including gene qualification documents from both the normal and melanoma sites, therefore the number of samples used to examine the results got from Oncomine was limited. Thirdly, we compared the levels of ten m6A regulatory genes between melanoma and the normal groups. Those ten genes included Writers as METTL14, Erasers as FTO, and Readers as YTHDF1, HNRNPA2B1, which were common genes that regulated m6A modification in RNA, while with the development new regulatory genes might emerge. Fourthly, the further experiments that explored the mechanism insight into the way m6A regulatory genes affected the development of melanoma needed to be performed in the future.


We systematically examined m6A-related genes in melanoma for the first time based on large datasets on publicly available databases. Several sources were used to confirm our results. Statistical analysis revealed the significance of YTHDF1 and HNRNPA2B1 and the combination of the two genes may provide a better approach to diagnose melanoma. Moreover, YTHDF1 or HNRNPA2B1 might interact with RNA processing-related genes such as CDK1/CDK2 and further promote the formation and progression of melanoma. This study may provide new targets for the treatment of melanoma and the evidence to explain the pathophysiology of the disease.

Availability of data and materials

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



B-Raf proto-oncogene


Mitogen-activated protein kinase


Also known as MAP2K7, mitogen-activated protein kinase kinase 7


Also known as MAPK1, mitogen-activated protein kinase 1




Methyltransferase like 3


Methyltransferase like 14


WT1 associated protein


FTO alpha-ketoglutarate dependent dioxygenase


alkB homolog 5


YTH N6-methyladenosine RNA binding protein 1


YTH N6-methyladenosine RNA binding protein 2


E74 like ETS transcription factor 3


Heterogeneous nuclear ribonucleoprotein C


Heterogeneous nuclear ribonucleoprotein A2/B1


Acute myeloid leukemia


Cervical squamous cell carcinoma


Gene expression omnibus


Receiver operating characteristic curve


Maximal clique centrality


Maximum neighborhood component


Edge percolated component


The cancer genome atlas


Confidence interval


Mann–Whitney test


The genes included in the correlation analysis


Fold change


False discovery rate


Area under curve


Protein post-translational modification




  1. 1.

    Nikolaou V, Stratigos AJ. Emerging trends in the epidemiology of melanoma. The British journal of dermatology. 2014;170(1):11–9.

    CAS  Article  Google Scholar 

  2. 2.

    Schadendorf D, Fisher DE, Garbe C, Gershenwald JE, Grob JJ, Halpern A, Herlyn M, Marchetti MA, McArthur G, Ribas A, et al. Melanoma. Nature reviews Disease primers. 2015;1:15003.

    Article  Google Scholar 

  3. 3.

    Egberts F, Kotthoff EM, Gerdes S, Egberts JH, Weichenthal M, Hauschild A. Comparative study of YKL-40, S-100B and LDH as monitoring tools for Stage IV melanoma. Eur J Cancer. 2012;48(5):695–702.

    CAS  Article  Google Scholar 

  4. 4.

    Flaherty KT, Hodi FS, Fisher DE. From genes to drugs: targeted strategies for melanoma. Nat Rev Cancer. 2012;12(5):349–61.

    CAS  Article  Google Scholar 

  5. 5.

    Flaherty KT, Infante JR, Daud A, Gonzalez R, Kefford RF, Sosman J, Hamid O, Schuchter L, Cebon J, Ibrahim N, et al. Combined BRAF and MEK inhibition in melanoma with BRAF V600 mutations. N Engl J Med. 2012;367(18):1694–703.

    CAS  Article  Google Scholar 

  6. 6.

    Zhao BS, Roundtree IA, He C. Post-transcriptional gene regulation by mRNA modifications. Nat Rev Mol Cell Biol. 2017;18(1):31–42.

    CAS  Article  Google Scholar 

  7. 7.

    Lee M, Kim B, Kim VN. Emerging roles of RNA modification: m(6)A and U-tail. Cell. 2014;158(5):980–7.

    CAS  Article  Google Scholar 

  8. 8.

    Zhang C, Fu J, Zhou Y. A Review in Research Progress Concerning m6A Methylation and Immunoregulation. Frontiers in immunology. 2019;10:922.

    CAS  Article  Google Scholar 

  9. 9.

    Meyer KD, Jaffrey SR. Rethinking m(6)A Readers, Writers, and Erasers. Annu Rev Cell Dev Biol. 2017;33:319–42.

    CAS  Article  Google Scholar 

  10. 10.

    Yu S, Li X, Liu S, Yang R, Liu X, Wu S. N(6)-Methyladenosine: a Novel RNA Imprint in Human Cancer. Frontiers in oncology. 2019;9:1407.

    Article  Google Scholar 

  11. 11.

    Vu LP, Pickering BF, Cheng Y, Zaccara S, Nguyen D, Minuesa G, Chou T, Chow A, Saletore Y, MacKay M, et al. The N(6)-methyladenosine (m(6)A)-forming enzyme METTL3 controls myeloid differentiation of normal hematopoietic and leukemia cells. Nat Med. 2017;23(11):1369–76.

    CAS  Article  Google Scholar 

  12. 12.

    Zhou S, Bai ZL, Xia D, Zhao ZJ, Zhao R, Wang YY, Zhe H. FTO regulates the chemo-radiotherapy resistance of cervical squamous cell carcinoma (CSCC) by targeting beta-catenin through mRNA demethylation. Mol Carcinog. 2018;57(5):590–7.

    CAS  Article  Google Scholar 

  13. 13.

    Bindea G, Mlecnik B, Hackl H, Charoentong P, Tosolini M, Kirilovsky A, Fridman WH, Pages F, Trajanoski Z, Galon J. ClueGO: a Cytoscape plug-into decipher functionally grouped gene ontology and pathway annotation networks. Bioinformatics (Oxford, England). 2009;25(8):1091–3.

    CAS  Article  Google Scholar 

  14. 14.

    Bindea G, Galon J, Mlecnik B. CluePedia Cytoscape plugin: pathway insights using integrated experimental and in silico data. Bioinformatics (Oxford, England). 2013;29(5):661–3.

    CAS  Article  Google Scholar 

  15. 15.

    Zhang ZN, Xu JJ, Fu YJ, Liu J, Jiang YJ, Cui HL, Zhao B, Sun H, He YW, Li QJ, et al. Transcriptomic analysis of peripheral blood mononuclear cells in rapid progressors in early HIV infection identifies a signature closely correlated with disease progression. Clin Chem. 2013;59(8):1175–86.

    CAS  Article  Google Scholar 

  16. 16.

    Xie Y, Zhang X, Zhao T, Li W, Xiang J. Natural CD8(+)25(+) regulatory T cell-secreted exosomes capable of suppressing cytotoxic T lymphocyte-mediated immunity against B16 melanoma. Biochem Biophys Res Commun. 2013;438(1):152–5.

    CAS  Article  Google Scholar 

  17. 17.

    Luo G, Xu W, Zhao Y, Jin S, Wang S, Liu Q, Chen X, Wang J, Dong F, Hu DN, et al. RNA m(6) A methylation regulates uveal melanoma cell proliferation, migration, and invasion by targeting c-Met. J Cell Physiol. 2020.

    Article  PubMed  Google Scholar 

  18. 18.

    Jia R, Chai P, Wang S, Sun B, Xu Y, Yang Y, Ge S, Jia R, Yang YG, Fan X. m(6)A modification suppresses ocular melanoma through modulating HINT2 mRNA translation. Molecular cancer. 2019;18(1):161.

    Article  Google Scholar 

  19. 19.

    Yang S, Wei J, Cui YH, Park G, Shah P, Deng Y, Aplin AE, Lu Z, Hwang S, He C, et al. m(6)A mRNA demethylase FTO regulates melanoma tumorigenicity and response to anti-PD-1 blockade. Nature communications. 2019;10(1):2782.

    Article  Google Scholar 

  20. 20.

    Liu T, Wei Q, Jin J, Luo Q, Liu Y, Yang Y, Cheng C, Li L, Pi J, Si Y, et al. The m6A reader YTHDF1 promotes ovarian cancer progression via augmenting EIF3C translation. Nucleic Acids Res. 2020;48:3816–31.

    Article  Google Scholar 

  21. 21.

    Han B, Yan S, Wei S, Xiang J, Liu K, Chen Z, Bai R, Sheng J, Xu Z, Gao X. YTHDF1-mediated translation amplifies Wnt-driven intestinal stemness. EMBO Rep. 2020;21:e49229.

    CAS  Article  Google Scholar 

  22. 22.

    Zhang H, Zhang X, Li X, Meng WB, Bai ZT, Rui SZ, Wang ZF, Zhou WC, Jin XD. Effect of CCNB1 silencing on cell cycle, senescence, and apoptosis through the p53 signaling pathway in pancreatic cancer. J Cell Physiol. 2018;234(1):619–31.

    Article  Google Scholar 

  23. 23.

    Payton M, Chung G, Yakowec P, Wong A, Powers D, Xiong L, Zhang N, Leal J, Bush TL, Santora V, et al. Discovery and evaluation of dual CDK1 and CDK2 inhibitors. Cancer Res. 2006;66(8):4299–308.

    CAS  Article  Google Scholar 

  24. 24.

    Chung JH, Bunz F. Cdk2 is required for p53-independent G2/M checkpoint control. PLoS Genet. 2010;6(2):e1000863.

    Article  Google Scholar 

  25. 25.

    Gali-Muhtasib H, Kuester D, Mawrin C, Bajbouj K, Diestel A, Ocker M, Habold C, Foltzer-Jourdainne C, Schoenfeld P, Peters B, et al. Thymoquinone triggers inactivation of the stress response pathway sensor CHEK1 and contributes to apoptosis in colorectal cancer cells. Cancer Res. 2008;68(14):5609–18.

    CAS  Article  Google Scholar 

  26. 26.

    Foskolou IP, Jorgensen C, Leszczynska KB, Olcina MM, Tarhonskaya H, Haisma B, D’Angiolella V, Myers WK, Domene C, Flashman E, et al. Ribonucleotide Reductase Requires Subunit Switching in Hypoxia to Maintain DNA Replication. Mol Cell. 2017;66(2):206-220 e209.

    Article  Google Scholar 

  27. 27.

    Gorgoulis VG, Vassiliou LV, Karakaidos P, Zacharatos P, Kotsinas A, Liloglou T, Venere M, Ditullio RA Jr, Kastrinakis NG, Levy B, et al. Activation of the DNA damage checkpoint and genomic instability in human precancerous lesions. Nature. 2005;434(7035):907–13.

    CAS  Article  Google Scholar 

  28. 28.

    Hammond EM, Kaufmann MR, Giaccia AJ. Oxygen sensing and the DNA-damage response. Curr Opin Cell Biol. 2007;19(6):680–4.

    CAS  Article  Google Scholar 

  29. 29.

    Kim JY, Lee JY. Targeting Tumor Adaption to Chronic Hypoxia: implications for Drug Resistance, and How It Can Be Overcome. Int J Mol Sci. 2017;18(9):1854.

    Article  Google Scholar 

  30. 30.

    Adamkova L, Souckova K, Kovarik J. Transcription protein STAT1: biology and relation to cancer. Folia Biol. 2007;53(1):1–6.

    CAS  Google Scholar 

  31. 31.

    Chen J, Zhao J, Chen L, Dong N, Ying Z, Cai Z, Ji D, Zhang Y, Dong L, Li Y, et al. STAT1 modification improves therapeutic effects of interferons on lung cancer cells. J Transl Med. 2015;13:293.

    Article  Google Scholar 

  32. 32.

    Dorritie KA, McCubrey JA, Johnson DE. STAT transcription factors in hematopoiesis and leukemogenesis: opportunities for therapeutic intervention. Leukemia. 2014;28(2):248–57.

    CAS  Article  Google Scholar 

  33. 33.

    Zhang X, Li X, Tan F, Yu N, Pei H. STAT1 Inhibits MiR-181a Expression to Suppress Colorectal Cancer Cell Proliferation Through PTEN/Akt. J Cell Biochem. 2017;118(10):3435–43.

    CAS  Article  Google Scholar 

  34. 34.

    Zhang Y, Molavi O, Su M, Lai R. The clinical and biological significance of STAT1 in esophageal squamous cell carcinoma. BMC Cancer. 2014;14:791.

    Article  Google Scholar 

  35. 35.

    Stopeck AT, Gessner A, Miller TP, Hersh EM, Johnson CS, Cui H, Frutiger Y, Grogan TM. Loss of B7.2 (CD86) and intracellular adhesion molecule 1 (CD54) expression is associated with decreased tumor-infiltrating T lymphocytes in diffuse B-cell large-cell lymphoma. Clin Cancer Res. 2000;6(10):3904–9.

    CAS  PubMed  Google Scholar 

Download references




This research is supported by grants from the National Natural Science Foundation of China (81671556, 81421002, 81670107, 81701608), the Tianjin Science Founds for Distinguished Young Scholar (17JCdJQJC45800), the Non-profit Central Research Institute Fund of Chinese Academy of Medical Sciences (2018PT32034 and 2019-RC-HL-013), China Postdoctoral Science Foundation (2017M623440).

Author information




TL: designed the study, downloaded and analyzed the data, wrote original draft, and drew the graphs; MG: analyzed the data, discussed with other authors; AD, obtained funding, reviewed and edited the manuscript; CQ, conducted systematic review, checked the results and edited the original draft. All authors read and approved the manuscript.

Corresponding authors

Correspondence to Anmei Deng or Cheng Qian.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare they have no competing conflicts of interest.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Table S1. Characteristics of patients with melanoma on Oncomine. No., number; y, year; n/a, not provided; * one was not provided; T, tumor. Table S2. Patient information for Manfred et al. 2018. No., number; NMM, nodular malignant melanoma; SSM, superficial spreading melanoma. Table S3. Interactions between the 10 target genes and their parameters. Q-value, adjusted P-value.

Additional file 2: Figure S1. Expression level and ROC curves of ELF3, and YTHDF2. The data was from the GEO dataset online. ROC, receiver operating characteristic curve. *P < 0.05, **P < 0.01. Figure S2. Altered characteristics for different genes in melanoma on cBioPortal. (a) Alteration frequency in different selected studies. (b) The mutation diagram and PTM sites for each gene. PTM, protein post-translational modification.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Li, T., Gu, M., Deng, A. et al. Increased expression of YTHDF1 and HNRNPA2B1 as potent biomarkers for melanoma: a systematic analysis. Cancer Cell Int 20, 239 (2020).

Download citation


  • Melanoma
  • YTHDF1
  • m6A
  • Systematic analysis
  • p53