Next Article in Journal
Variability in Breast Cancer Biomarker Assessment and the Effect on Oncological Treatment Decisions: A Nationwide 5-Year Population-Based Study
Next Article in Special Issue
Roles of microRNAs in Regulating Cancer Stemness in Head and Neck Cancers
Previous Article in Journal
Epithelial to Mesenchymal Transition Regulates Surface PD-L1 via CMTM6 and CMTM7 Induction in Breast Cancer
Previous Article in Special Issue
Tumor Suppressor miR-584-5p Inhibits Migration and Invasion in Smoking Related Non-Small Cell Lung Cancer Cells by Targeting YKT6
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Relevance of 2′-O-Methylation and Pseudouridylation for the Malignant Melanoma

by
Simon Jasinski-Bergner
1,*,
Juliane Blümke
1,
Claudia Wickenhauser
2 and
Barbara Seliger
1,3
1
Institute for Medical Immunology, Faculty of Medicine, Martin Luther University Halle-Wittenberg, Magdeburger Straße 2, 06108 Halle (Saale), Germany
2
Institute for Pathology, Faculty of Medicine, Martin Luther University Halle-Wittenberg, 06108 Halle (Saale), Germany
3
Fraunhofer Institute for Cell Therapy and Immunology, 04103 Leipzig, Germany
*
Author to whom correspondence should be addressed.
Cancers 2021, 13(5), 1167; https://doi.org/10.3390/cancers13051167
Submission received: 12 February 2021 / Accepted: 4 March 2021 / Published: 9 March 2021
(This article belongs to the Special Issue MicroRNA and Cancer)

Abstract

:

Simple Summary

This study investigates the expression, the histological localization, and the influence of the factors involved in 2′-O-methylation and pseudouridylation on prognostic relevant markers, proliferation markers, overall survival, molecular immune surveillance and evasion mechanisms within the malignant melanoma. Statistically significant positive correlations to the expression of markers involved in cell proliferation were observed. The upregulation of the RNA modifying factors was of prognostic relevance in this tumor disease with a negative impact on the overall survival of melanoma patients. Furthermore, the factors involved in 2′-O-methylation and pseudouridylation were statistically significant negative correlated to the expression of human leukocyte antigen class I genes as well as of components of the antigen processing machinery.

Abstract

The two RNA modifications 2′-O-methylation and pseudouridylation occur on several RNA species including ribosomal RNAs leading to an increased translation as well as cell proliferation associated with distinct functions. Using malignant melanoma (MM) as a model system the proteins mediating these RNA modifications were for the first time analyzed by different bioinformatics tools and public available databases regarding their expression and histological localization. Next to this, the impact of these RNA-modifying factors on prognostic relevant processes and marker genes of malignant melanoma was investigated and correlated to immune surveillance and evasion strategies. The RNA modifying factors exerted statistically significant positive correlations to the expression of genes involved in cell proliferation and were statistically significant negative correlated to the expression of human leukocyte antigen class I genes as well as of components of the antigen processing machinery in malignant melanoma. Upregulation of the RNA modifying proteins was of prognostic relevance in this tumor disease with a negative impact on the overall survival of melanoma patients. Furthermore, the expression of known oncogenic miRs, which are induced in malignant melanoma, directly correlated to the expression of factors involved in these two RNA modifications.

1. Introduction

Malignant melanoma (MM) refers to a neoplasm of melanocytes comprising of neural crest-derived cells located in particular in the stratum basal of the skin’s epidermis and in the uvea of the eye. In contrast to keratinocytes, melanocytes are not linked to the basal lamina by hemidesmosomes and do not have desmosomes to neighboring keratinocytes. Instead, melanocytes as well as adjacent keratinocytes express, e.g., E-cadherin on their cell surface [1], which is an important factor contributing to the cell migration and invasion of melanoma cells after malignant transformation of melanocytes. Due to this migratory ability, growth characteristics as well as other properties including resistance to radiation and certain chemotherapeutics [2], MM represents the deadliest type of skin cancer with increasing incidences [3]. UV light exposition leading to DNA damages and oxidative stress in melanocytes [4], nevus number, pigmentation characteristics as well as genetic mutations are known predispositions for the processes leading to malignant transformation [5]. Concerning the hereditary predisposition, it is known that between 5–10% of melanoma cases are familial [6] and mainly induced by germ line mutations in tumor suppressor genes involved in proliferation control and DNA repair of cells.
Benign nevi defined as growth arrested, clonal neoplasms of melanocytes initiated by well-defined oncogenic mutations in the mitogen-activated protein kinase pathway, such as the BRAFV600E-activating mutation [7], may develop dysplastic features and are classified as MM precursor lesions.
An essential impact on the cell proliferation per se and exemplarily referred to MM is confined by ribosomes. An average mammalian cell has between 5 and 10 million ribosomes, which are large multi-subunit ribonucleoprotein complexes and are required for protein synthesis. Therefore, the formation of ribosomes and their linked translational efficiency (TE) directly influence cell proliferation, metabolism and vitality—parameters that are modulated upon malignant transformation [8] and indicative for neoplastic progression [9]. Thus, an enhanced TE is directly proportional to the ribosomal density, since the number of ribosomes on a coding transcript is linked to the efficacy of translation [10,11,12]. Indeed, an increase of ribosomes has been described in various cancers suggesting that the number and modifications of ribosomes drive tumorigenesis [13].
The size of the nucleoli as major location of the ribosomal assembly represents one important cytomorphologic parameter for the determination of malignancy in melanocytic lesions and has prognostic relevance [14]. Interestingly, expression of MYC as well as mutations in the tumor suppressor genes TP53 and RB1 enhanced ribosomal biogenesis [14,15].
Next to processes of transcription, the splicing and assembly of the eukaryotic 80 S ribosomes RNA modifications including 2′-O-methylation as well as the pseudouridylation occur in nucleoli, which involve a class of non-coding small RNAs termed sno-RNAs. Both processes represent the most common modifications of rRNAs and are important steps for their maturation as well as stabilization. In mammalian cells, the 5.8S rRNA, the 18S and the 28S rRNA have in total > 100 2′-O-methylations and 95 pseudouridinylations [16,17]. In addition to the rRNAs, the class of spliceosomal snoRNAs contain both modifications underlining the influence of snoRNAs for the splicing process [18]. The altered splicing pattern after malignant transformation involve mRNAs of many cancer related and tumor biological important genes [19]. In addition, transfer RNA (tRNA) molecules, microRNAs (miRs) and messenger RNAs (mRNAs) can be modified by 2′-O-methylation [20].
For the 2′-O-methylation, a methyl residue is added to the ribose backbone, whereby fibrillarin (FBL) is the methyltransferase using S-adenosyl methionine (SAM) as a methyl donor (Figure 1A). FBL acts in a complex with nucleolar protein (NOP) 56, NOP58, 15.5K (SNU13) and the guide RNA [21]. The introduction of the methyl residue leads to steric alterations and increases the hydrophobicity thereby protecting RNA molecules from nucleolytic attacks [22]. The sum of all 2′-O-methylations within a RNA molecule can also affect the secondary structure and therefore possible interactions including RNA/RNA interactions, RNA-protein interactions [23,24] as well as mRNA splicing, stability and translation [25].
Many sites of the pseudouridylation are evolutionarily conserved [26,27]. In humans, the uridine is converted to pseudouridine in rRNA, sno/scaRNA and snRNA by dyskerin (DKC1), which is a component of a complex consisting of one H/ACA snoRNA and four core proteins, namely GAR1 ribonucleoprotein (GAR1), NHP2 ribonucleoprotein (NHP2), NOP10 ribonucleoprotein (NOP10) and DKC1 itself [28]. The substitution of uridine with pseudouridine introduces a novel H bond donor on the non-Watson Crick site of the nucleotide affecting the secondary structure and consequently structure related interactions [29] (Figure 1B). Pseudouridylation is involved in many important processes of gene expression including spliceosomal small nuclear ribonucleoprotein (snRNP) biogenesis, efficiency of pre-mRNA splicing and translation fidelity [30]. Another human pseudouridine synthase is PUS10 involved in the miR maturation, which results after its depletion in a reduced expression of mature miRs [31]. Furthermore, the 2′-O-methylation and the pseudouridylation have also an impact on the binding of RNA-binding proteins (RBPs) and are required for an appropriate mRNA splicing [25,30], which is altered in various tumor entities including MM [32,33,34]. In this study, the expression of key molecules involved in 2′-O-methylation and pseudouridylation is addressed in more detail using MM as a model.

2. Results

2.1. Expression Pattern and Localization of RNA-Modifying Proteins in the Skin

According to the information provided by The Human Protein Atlas, the protein expression involved in the pseudouridylation including DKC1, GAR1, NHP2 and NOP10 were only localized in the nucleus, while the staining of the proteins involved in 2′-O-methylation, namely FBL, NOP56, NOP58 as well as 15.5K (SNU13) revealed a more heterogenic localization summarized in Table 1. FBL, NOP56 and NOP58 were located in the nucleus and with the exception of NOP56 strongly expressed in the epidermis and in the highly proliferation active cells of the stratum germinativum. In contrast, the staining of 15.5K revealed a cytoplasmic and membranous expression almost exclusively in cells of the epidermis layers of the skin (Table 1, Figure 2).
The correlations of factors involved in 2′-O-methylation and pseudouridylation with melanoma relevant proliferation markers, prognostic markers, and with genes responsible for immune surveillance as well as for immune evasion and positive as well as negative regulations were seen as indicated and visualized in the summarizing heatmap in Figure 3. The respective R and p values are listed within Table 2, Table 3 and Table 4. The results will be addressed in more details in the following paragraphs.

2.2. Correlation of RNA-Modifying Proteins with Tumor Cell Proliferation

To address whether the expression of factors involved in 2′-O-methylation and pseudouridylation correlates with the pathological increased proliferation rates in MM, the expression data of these RNA-modifying factors were analyzed in The Cancer Genome Atlas (TCGA) data sets consisting of 214 samples from MM patients provided by the R2 data base. The following proliferation markers known to play a role in MM were investigated: Ki67 (MKI67), PCNA, cyclin A (CCNA1), cyclin B (CCNB1), MCM2, MCM4, and mitosin (CENPF) [35,36,37,38].
As summarized in Table 2, positive correlations between both RNA modifying factors and the proliferation markers analyzed were found. The effect was the strongest for MCM4, which correlated statistically significant to all factors involved into 2′-O-methylation and pseudouridylation with the exception of NOP10. In addition, PCNA, cyclin B and mitosin showed strong correlations, while the expression of CCNA1 did not correlate statistically significant with any of these factors.
Based on the strong positive correlation between the factors involved in 2′-O-methylation and pseudouridylation with most of the clinical relevant proliferation markers, a possible correlation between these factors and prognostic marker genes relevant for MM was evaluated such as melanoma antigen recognized by T cells (MART)1, S100 calcium binding protein B (S100B), S100A4, S100A9, melanocyte inducing transcription factor (MITF), matrix metallopeptidase 2 (MMP2), nucleoside diphosphate kinase 1 (NM23), cluster of differentiation (CD) 44, premelanosome protein (PMEL), and BCL2 apoptosis regulator (BCL2) using the same TCGA data sets [39,40,41,42,43] (Table 3). Next to these marker genes the invasion depth termed Breslow’s depth is of prognostic value. The average Breslow score of the 214 MM patients analyzed was 2.5 mm. All other clinical parameters of this MM patient cohort (n = 214) of the analyzed microarray data can be obtained from the original literature Jönsson et al., 2015 [44].
The two markers MART1 and MITF, but also NM23 showed statistically significant positive correlations to the expression of RNA modifying proteins. In contrast, the S100 family members S100B, S100A4, S100A9 and MMP2 exhibiting a prognostic potential in MM exerted a statistically significant negative correlation, whereas CD44 did not show any correlation and PMEL as well as BCL2 were only weakly positivey correlated.
To address whether the factors involved in 2′-O-methylation and pseudouridylation are of central importance for prognosis in MM, the same TCGA data set was used for the generation of Kaplan–Meier plots. As shown in Figure 4A–H, a direct correlation of the expression level of these factors with the overall survival (OS) of MM patients was detected and statistically significant for FBL, NOP58, and GAR1. In contrast, low NOP10 expression levels predicted a statistically non-significant correlation with a worse patients’ outcome (Figure 4). However, it is noteworthy that the use of TCGA data cannot replace extensive protein based analyses of ex vivo MM specimens.

2.3. Correlation of the Expression of RNA Modifying Factors with Immune Modulatory Genes

Due to the increased implementation of immunotherapies for the treatment of hematopoietic and solid tumors, it was analyzed whether the factors involved in 2′-O-methylation and pseudouridylation directly or indirectly affect transcripts participating in the immune surveillance of tumor cells. Therefore, their expression patterns were correlated to molecules involved in the immune recognition/evasion of malignant and/or virus transformed cells from immune effector cells. Using the same melanoma data set, an impressive statistically significant negative correlation between immunological relevant molecules, in particular the human leucocyte antigen (HLA) class Ia and Ib and components of the HLA class I antigen processing machinery (APM), with the factors involved in 2′-O-methylation and pseudouridylation was detected (Table 4). This might also explain the negative correlation of these factors with the OS of MM patients, since their impaired expression was associated with a reduced anti-tumoral immune cell response [45] (Figure 4). The divergent statistically significant positive correlation pattern of NOP10 to molecules involved in immune surveillance was also associated with a statistically significant positive correlation to the OS in MM patients (Figure 4H). However, further studies are needed to investigate whether the statistically significant negative correlated immunological relevant genes are directly negatively regulated upon a reduced activity/expression of the factors involved 2′-O-methylation and pseudouridylation.

2.4. Correlation of miR Expression with RNA-Modifying Factors

The processes of 2′-O-methylation and pseudouridylation also occur in other RNA species with a strong impact to diverse molecular processes of malignant transformation and thus are of clinical relevance. These other RNA species include microRNAs (miRs), which are non-coding single stranded RNAs with an approximately length of ~22 nt [46], binding specifically and preferentially, but not exclusively, within the 3′- untranslated region (UTR) to their target mRNAs sequence [40,47] thereby causing a translational repression and mRNA decay [48].
Some miRs can be classified into oncogenic, tumor suppressive and/or immune modulatory miRs [48]. In addition, human viruses including herpes viruses encode for viral miRs, which also affect cancer related cellular processes [49].
Several oncogenic miRs have been reported to be frequently overexpressed in MM. Fattore and co-authors (2017) even grouped the most representative deregulated miRs in melanoma with regard to the different steps of tumor progression. These potentially oncogenic miRs with putative diagnostic and/or prognostic values include miR-9, miR-10a, miR-10b, miR-17-5p, miR-18a, miR-21, miR-26b, miR-92a, miR-221, miR-222, miR-126, miR-145, miR146, miR-182, miR-514, miR-520d and miR-527 [50].
The pseudouridylation of different human RNA species is mediated by different pseudouridine synthases, such as pseudouridine synthase (PUS) 1, TruB pseudouridine synthase 2 (TRUB2), TRUB1, PUS3, PUS4, RNA pseudouridine synthase D3 (RPUSD3), RPUSD4, PUS7, pseudouridine synthase 7 like (PUS7L), PUS10 and DKC1 [28]. From these candidates, TRUB1 has been recently reported to modulate the stability of hsa-let-7 [51]. However, it is noteworthy that the actual list of enzymes involved in 2′-O-methylation and pseudouridylation of miRs might be incomplete.
Thus far, only little is known, which of these enzymes are able to modify miRs and which alterations in the respective miR-mediated functions and/or miR stability are caused by such modifications.
To determine whether the expression of known oncogenic miRs in MM might be correlated with the expression of enzymes involved into 2′-O-methylation and pseudouridylation, the TCGA data of 214 human MM samples were further analyzed for miR expression. As shown in Table 5, there existed no evidence for a global impact of the induction of oncogenic miRs, while miR-21 was statistically significant negatively correlated to the expression of the investigated enzymes in MM.
Analysis of the TCGA data from 214 human melanoma tumors provided by the R2 database (https://hgserver1.amc.nl/ (accessed on 12 February 2021)) for correlation with enzymes known to be involved in miR 2′-O-methylation and pseudouridylation with the amount of increased and/or stabilized oncogenic miRs reported to be induced in MM.

3. Discussion

Recently, the aberrant expression pattern of DKC1, NHP2, and NOP10 in several cancer entities has been reviewed [52]. Due to their biological functions, the localization of these factors was almost completely limited to the nucleus, since 2′-O-methylation and pseudouridylation occur within the nucleoli. Furthermore, an increase within the highly proliferational active cells of the epidermal stratum germinativum was demonstrated, which reflects the involvement and increased number of ribosomes during proliferation. This was further underlined by statistically significant positive correlations to proliferation markers suggesting that some factors involved in 2′-O-methylation and pseudouridylation might also be suitable markers for proliferation themselves.
However, this study did not answer whether in the case of positive correlations between the factors involved in 2′-O-methylation and pseudouridylation and, e.g., proliferation markers these positive correlations are an indirect result due to an enhanced proliferation per se mediated by increased translational efficacy caused by modifications within ribosomal RNAs, or whether the mRNAs of the positively correlated proliferation markers are directly modified by these factors or by both mechanisms in parallel.
The positive correlations described in this study are based upon TCGA data, not taken posttranscriptional mechanisms of gene regulation into account. Therefore, in depth protein-based studies applying human melanoma tissue specimens are required to proof this hypothesis.
In regard to the observed statistically significant negative and coordinated correlations to molecules involved in immune surveillance/evasion it is noteworthy that (i) HLA Ia and Ib genes, HLA class related MICs and major APM components are located within the major histocompatibility complex (MHC) locus on chromosome 6p21 and that (ii) the adaptive immunity phylogenetically occurs with the jawed vertebrates [53] and, e.g., HLA class Ib genes are evolutionary even younger, which might have an impact on the observed negative correlations.
Due to the limitation of transcriptome based microarray data sets it is necessary to underline that also indirect effects may cause such down regulations. The increased proliferation itself might decrease the antigen processing and presentation efficacy in melanoma cells. In turn, the impaired antigen processing and presentation is a frequent immune evasion strategy of tumors [54], which might be a prerequisite for an increased proliferation rate of tumor cells. Thus far, the sequence of these events in MM has not been studied.
RNA modifications on miRs affecting their stability or decay might also indirectly contribute to the regulation of cancer related processes, even independently of the translational processes linked to the number of ribosomes. This opens a new regulatory dimension for RNA modifying enzymes (Figure 5), which has to be explored in more detail.
Thus far, it is controversially discussed whether mature mammalian miRs contain 2′-O-methylations as reported for plants and Drosophila [55]. In Drosophila, 2′-O-methylation of miRs occurred age dependent and its inhibition resulted in an accelerated neurodegeneration and shorter life span [56]. Despite discrepancies regarding the missing 2′-O-methylation in mammalian miRs, the human miR-21-5p has been shown to contain a 3′-terminal 2′-O-methylation, which enhances the stability of this oncogenic miR in lung cancer patients. Interestingly, HENMT1 was identified to act as methyltransferase [57].
Indirect effects concerning the positive or negative correlations of factors involved in 2′-O-methylation and pseudouridylation with miRs might occur. Several miRs are processed from introns after the splicing, while for the process of splicing the snoRNAs as huge group of the snRNAs in involved. These snoRNAs harbor themselves 2′-O-methylations and pseudouridylations contributing to the stability and functionality of the snoRNAs.
In this study, the impact and relevance of factors involved in 2′-O-methylation and pseudouridylation of different RNA species related to the processes of tumor formation and progression were summarized in MM as a tumor model. Using different bioinformatics tools statistically significant positive correlations between proliferation and prognostic marker genes relevant for the MM with the factors involved into 2′-O-methylation and pseudouridylation were described for the first time.
The impact of these molecules on other RNA species including miRs has recently been investigated, but needs to be further studied in detail. This will lead to the identification of proteins involved in such miR modifications, which might have an impact on the function of such modified miRs. However, the existence of miR 2′-O-methylation and pseudouridylation in general must not necessarily involve tumor relevant miRs and processes.

4. Materials and Methods

Immunohistochemistry

The frequency and localization of protein levels involved in the two RNA modification processes, namely FBL, NOP56, NOP58, 15.5K (SNU13), DKC1, GAR1, NHP2, and NOP10, were analyzed by evaluation of immunohistochemistry data of healthy normal skin sections and melanoma specimen provided by the free available online data base. The Human Protein Atlas [58,59,60].
Bioinformatic analyses of gene expression data and correlation with patients’ overall survival.
The factors involved in RNA 2′-O-methylation and pseudouridylation were correlated with genes encoding for above mentioned ribosomal RNAs carrying multiple of such modifications. Unfortunately, the applied Microarrays (Illumina Human HT-12V4.0 Chips) of the 214 MM patients published by Jönsson and co-authors 2015 [44] contained only probes against the human RN5S9 transcript encoding the 5S ribosomal RNA. The factors showed a positive correlation to the RN5S9, which was statistically significant for NOP56, NOP58 and NHP2.
The correlation coefficients reflecting R values and the respective p values were calculated by the R2 database (https://hgserver1.amc.nl/ (accessed on 12 February 2021)). For visualization these R values were presented in a heatmap generated by GraphPad Prism 8 (GraphPad Software, San Diego, CA, USA).
Correlation of the expression pattern of genes involved into 2′-O-methylation and pseudouridylation with preselected genes involved in cell proliferation, prognosis and immune recognition/evasion in MM is based upon TCGA data of 214 melanoma specimens provided by the R2 database (https://hgserver1.amc.nl/ (accessed on 12 February 2021)). Statistically significant positive or negative correlations (p < 0.05) are highlighted in red or green. The same TCGA data set was applied for the generation of the respective Kaplan–Meier Plots.

5. Conclusions

Thus, a link between 2′-O-methylation and pseudouridylation to cell proliferation, host immunity and oncogenic miRs exists in MM suggesting that both RNA modifications and factors involved in this process represent suitable targets for tumor therapy and putative novel prognostic markers.

Author Contributions

S.J.-B. designed the study. J.B. and S.J.-B. performed the bioinformatics analyses. B.S., C.W., and S.J.-B. wrote the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the German-Israeli Foundation (GIF; I-37-414.11-2016) and Werner Jackstädt Foundation. We acknowledge the financial support of the Open Access Publication Fund of the Martin-Luther-University Halle-Wittenberg.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

APM, antigen processing machinery; BCL2, BCL2 apoptosis regulator; CCNA1, cyclin A1; CD44, cluster of differentiation 44; CENPF, centromere protein F; EMT, epithelial to mesenchymal transition; GAR1, GAR1 ribonucleoprotein; HLA, human leukocyte antigen; lncRNA, long non-coding RNA; MAPK, mitogen-activated protein kinase; MART-1, Melanoma Antigen Recognized by T-cells; MCM2, minichromosome maintenance complex component 2; MCM4, minichromosome maintenance complex component 4; MHC, major histocompatibility complex; MKI67, marker of proliferation Ki-67; MITF, melanocyte inducing transcription factor; MM, malignant melanoma; MMP2, matrix metallopeptidase 2; miR, microRNA; mRNA, messenger RNA; n.d., not determined; NHP2, NHP2 ribonucleoprotein; NM23, NME/NM23 nucleoside diphosphate kinase 1; NOP, nucleolar protein; PCNA, proliferating cell nuclear antigen; PMEL; premelanosome protein; PUS1, pseudouridine synthase 1; PUS7L, pseudouridine synthase 7 like; RBP, RNA-binding protein; RPUSD3, RNA pseudouridine synthase D3; rRNA, ribosomal RNA; S, Svedberg; S100B, S100 calcium binding protein B; scaRNA, small cajal body RNA; sncRNA, small non-coding RNA; snRNA, small nuclear RNA; SNU13, small nuclear ribonucleoprotein 13; TCGA, The Cancer Genome Atlas; TE, translational efficiency; tRNA, transfer RNA; UV, ultraviolet.

References

  1. Haass, N.K.; Herlyn, M. Normal human melanocyte homeostasis as a paradigm for understanding melanoma. J. Investig. Dermatol. Symp. Proc. 2005, 10, 153–163. [Google Scholar] [CrossRef] [Green Version]
  2. Laikova, K.V.; Oberemok, V.V.; Krasnodubets, A.M.; Gal’chinsky, N.V.; Useinov, R.Z.; Novikov, I.A.; Temirova, Z.Z.; Gorlov, M.V.; Shved, N.A.; Kumeiko, V.V.; et al. Advances in the Understanding of Skin Cancer: Ultraviolet Radiation, Mutations, and Antisense Oligonucleotides as Anticancer Drugs. Molecules 2019, 24, 1516. [Google Scholar] [CrossRef] [Green Version]
  3. Bray, F.; Ren, J.S.; Masuyer, E.; Ferlay, J. Global estimates of cancer prevalence for 27 sites in the adult population in 2008. Int. J. Cancer 2013, 132, 1133–1145. [Google Scholar] [CrossRef] [PubMed]
  4. Sample, A.; He, Y.Y. Mechanisms and prevention of UV-induced melanoma. Photodermatol. Photoimmunol. Photomed. 2018, 34, 13–24. [Google Scholar] [CrossRef] [PubMed]
  5. Hawkes, J.E.; Truong, A.; Meyer, L.J. Genetic predisposition to melanoma. Semin. Oncol. 2016, 43, 591–597. [Google Scholar] [CrossRef]
  6. Debniak, T. Familial malignant melanoma—Overview. Hered. Cancer Clin. Pract. 2004, 2, 123–129. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Damsky, W.E.; Bosenberg, M. Melanocytic nevi and melanoma: Unraveling a complex relationship. Oncogene 2017, 36, 5771–5792. [Google Scholar] [CrossRef] [Green Version]
  8. Stepinski, D. The nucleolus, an ally, and an enemy of cancer cells. Histochem. Cell Biol. 2018, 150, 607–629. [Google Scholar] [CrossRef]
  9. Pierard, G.E. Cell proliferation in cutaneous malignant melanoma: Relationship with neoplastic progression. ISRN Dermatol. 2012, 2012, 828146. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Szavits-Nossan, J.; Ciandrini, L. Inferring efficiency of translation initiation and elongation from ribosome profiling. Nucleic Acids Res. 2020, 48, 9478–9490. [Google Scholar] [CrossRef]
  11. von der Haar, T. Mathematical and Computational Modelling of Ribosomal Movement and Protein Synthesis: An overview. Comput. Struct. Biotechnol. J. 2012, 1, e201204002. [Google Scholar] [CrossRef] [Green Version]
  12. Szavits-Nossan, J.; Ciandrini, L. Accurate measures of translational efficiency and traffic using ribosome profiling. bioRxiv 2019. [Google Scholar] [CrossRef] [Green Version]
  13. Pelletier, J.; Thomas, G.; Volarevic, S. Ribosome biogenesis in cancer: New players and therapeutic avenues. Nat. Rev. Cancer 2018, 18, 51–63. [Google Scholar] [CrossRef] [PubMed]
  14. Penzo, M.; Montanaro, L.; Trere, D.; Derenzini, M. The Ribosome Biogenesis-Cancer Connection. Cells 2019, 8, 55. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Brighenti, E.; Calabrese, C.; Liguori, G.; Giannone, F.A.; Trere, D.; Montanaro, L.; Derenzini, M. Interleukin 6 downregulates p53 expression and activity by stimulating ribosome biogenesis: A new pathway connecting inflammation to cancer. Oncogene 2014, 33, 4396–4406. [Google Scholar] [CrossRef] [Green Version]
  16. Maden, B.E. The numerous modified nucleotides in eukaryotic ribosomal RNA. Prog. Nucleic Acid Res. Mol. Biol. 1990, 39, 241–303. [Google Scholar] [CrossRef]
  17. Ofengand, J.; Bakin, A.; Wrzesinski, J.; Nurse, K.; Lane, B.G. The pseudouridine residues of ribosomal RNA. Biochem. Cell Biol. 1995, 73, 915–924. [Google Scholar] [CrossRef]
  18. Jady, B.E.; Kiss, T. A small nucleolar guide RNA functions both in 2′-O-ribose methylation and pseudouridylation of the U5 spliceosomal RNA. EMBO J. 2001, 20, 541–551. [Google Scholar] [CrossRef] [PubMed]
  19. Silipo, M.; Gautrey, H.; Tyson-Capper, A. Deregulation of splicing factors and breast cancer development. J. Mol. Cell Biol. 2015, 7, 388–401. [Google Scholar] [CrossRef] [Green Version]
  20. Somme, J.; Van Laer, B.; Roovers, M.; Steyaert, J.; Versees, W.; Droogmans, L. Characterization of two homologous 2′-O-methyltransferases showing different specificities for their tRNA substrates. RNA 2014, 20, 1257–1271. [Google Scholar] [CrossRef] [Green Version]
  21. Peng, Y.; Yu, G.; Tian, S.; Li, H. Co-expression and co-purification of archaeal and eukaryal box C/D RNPs. PLoS ONE 2014, 9, e103096. [Google Scholar] [CrossRef]
  22. Sproat, B.S.; Lamond, A.I.; Beijer, B.; Neuner, P.; Ryder, U. Highly efficient chemical synthesis of 2′-O-methyloligoribonucleotides and tetrabiotinylated derivatives; novel probes that are resistant to degradation by RNA or DNA specific nucleases. Nucleic Acids Res. 1989, 17, 3373–3386. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Inoue, H.; Hayase, Y.; Imura, A.; Iwai, S.; Miura, K.; Ohtsuka, E. Synthesis and hybridization studies on two complementary nona(2′-O-methyl)ribonucleotides. Nucleic Acids Res. 1987, 15, 6131–6148. [Google Scholar] [CrossRef]
  24. Lacoux, C.; Di Marino, D.; Boyl, P.P.; Zalfa, F.; Yan, B.; Ciotti, M.T.; Falconi, M.; Urlaub, H.; Achsel, T.; Mougin, A.; et al. BC1-FMRP interaction is modulated by 2′-O-methylation: RNA-binding activity of the tudor domain and translational regulation at synapses. Nucleic Acids Res. 2012, 40, 4086–4096. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Elliott, B.A.; Ho, H.T.; Ranganathan, S.V.; Vangaveti, S.; Ilkayeva, O.; Abou Assi, H.; Choi, A.K.; Agris, P.F.; Holley, C.L. Modification of messenger RNA by 2′-O-methylation regulates gene expression in vivo. Nat. Commun. 2019, 10, 3401. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Li, X.; Ma, S.; Yi, C. Pseudouridine: The fifth RNA nucleotide with renewed interests. Curr. Opin. Chem. Biol. 2016, 33, 108–116. [Google Scholar] [CrossRef] [PubMed]
  27. van der Feltz, C.; DeHaven, A.C.; Hoskins, A.A. Stress-induced Pseudouridylation Alters the Structural Equilibrium of Yeast U2 snRNA Stem II. J. Mol. Biol. 2018, 430, 524–536. [Google Scholar] [CrossRef]
  28. Penzo, M.; Guerrieri, A.N.; Zacchini, F.; Trere, D.; Montanaro, L. RNA Pseudouridylation in Physiology and Medicine: For Better and for Worse. Genes 2017, 8, 301. [Google Scholar] [CrossRef] [Green Version]
  29. Ge, J.; Yu, Y.T. RNA pseudouridylation: New insights into an old modification. Trends Biochem. Sci. 2013, 38, 210–218. [Google Scholar] [CrossRef] [Green Version]
  30. Zhao, Y.; Dunker, W.; Yu, Y.T.; Karijolich, J. The Role of Noncoding RNA Pseudouridylation in Nuclear Gene Expression Events. Front. Bioeng. Biotechnol. 2018, 6, 8. [Google Scholar] [CrossRef] [Green Version]
  31. Song, J.; Zhuang, Y.; Zhu, C.; Meng, H.; Lu, B.; Xie, B.; Peng, J.; Li, M.; Yi, C. Differential roles of human PUS10 in miRNA processing and tRNA pseudouridylation. Nat. Chem. Biol. 2020, 16, 160–169. [Google Scholar] [CrossRef] [PubMed]
  32. Chen, J.; Weiss, W.A. Alternative splicing in cancer: Implications for biology and therapy. Oncogene 2015, 34, 1–14. [Google Scholar] [CrossRef] [PubMed]
  33. Ma, F.C.; He, R.Q.; Lin, P.; Zhong, J.C.; Ma, J.; Yang, H.; Hu, X.H.; Chen, G. Profiling of prognostic alternative splicing in melanoma. Oncol. Lett. 2019, 18, 1081–1088. [Google Scholar] [CrossRef] [Green Version]
  34. Singh, B.; Eyras, E. The role of alternative splicing in cancer. Transcription 2017, 8, 91–98. [Google Scholar] [CrossRef] [Green Version]
  35. Guzinska-Ustymowicz, K.; Pryczynicz, A.; Kemona, A.; Czyzewska, J. Correlation between proliferation markers: PCNA, Ki-67, MCM-2 and antiapoptotic protein Bcl-2 in colorectal cancer. Anticancer Res. 2009, 29, 3049–3052. [Google Scholar]
  36. Ladstein, R.G.; Bachmann, I.M.; Straume, O.; Akslen, L.A. Ki-67 expression is superior to mitotic count and novel proliferation markers PHH3, MCM4 and mitosin as a prognostic factor in thick cutaneous melanoma. BMC Cancer 2010, 10, 140. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Ohsie, S.J.; Sarantopoulos, G.P.; Cochran, A.J.; Binder, S.W. Immunohistochemical characteristics of melanoma. J. Cutan. Pathol. 2008, 35, 433–444. [Google Scholar] [CrossRef]
  38. Reddy, V.B.; Gattuso, P.; Aranha, G.; Carson, H.J. Cell proliferation markers in predicting metastases in malignant melanoma. J. Cutan. Pathol. 1995, 22, 248–251. [Google Scholar] [CrossRef]
  39. Bresnick, A.R.; Weber, D.J.; Zimmer, D.B. S100 proteins in cancer. Nat. Rev. Cancer 2015, 15, 96–109. [Google Scholar] [CrossRef] [Green Version]
  40. Fander, J.; Kielstein, H.; Buttner, M.; Koelblinger, P.; Dummer, R.; Bauer, M.; Handke, D.; Wickenhauser, C.; Seliger, B.; Jasinski-Bergner, S. Characterizing CD44 regulatory microRNAs as putative therapeutic agents in human melanoma. Oncotarget 2019, 10, 6509–6525. [Google Scholar] [CrossRef] [Green Version]
  41. Ilmonen, S.; Hernberg, M.; Pyrhonen, S.; Tarkkanen, J.; Asko-Seljavaara, S. Ki-67, Bcl-2 and p53 expression in primary and metastatic melanoma. Melanoma Res. 2005, 15, 375–381. [Google Scholar] [CrossRef] [PubMed]
  42. Weinstein, D.; Leininger, J.; Hamby, C.; Safai, B. Diagnostic and prognostic biomarkers in melanoma. J. Clin. Aesthet. Dermatol. 2014, 7, 13–24. [Google Scholar] [PubMed]
  43. Kycler, W.; Grodecka-Gazdecka, S.; Breborowicz, J.; Filas, V.; Teresiak, M. Prognostic factors in melanoma. Rep. Pract. Oncol. Radiother. 2006, 11, 39–48. [Google Scholar] [CrossRef] [Green Version]
  44. Cirenajwis, H.; Ekedahl, H.; Lauss, M.; Harbst, K.; Carneiro, A.; Enoksson, J.; Rosengren, F.; Werner-Hartman, L.; Torngren, T.; Kvist, A.; et al. Molecular stratification of metastatic melanoma using gene expression profiling: Prediction of survival outcome and benefit from molecular targeted therapy. Oncotarget 2015, 6, 12297–12309. [Google Scholar] [CrossRef] [Green Version]
  45. Quandt, D.; Fiedler, E.; Boettcher, D.; Marsch, W.; Seliger, B. B7-h4 expression in human melanoma: Its association with patients’ survival and antitumor immune response. Clin. Cancer Res. 2011, 17, 3100–3111. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Starega-Roslan, J.; Krol, J.; Koscianska, E.; Kozlowski, P.; Szlachcic, W.J.; Sobczak, K.; Krzyzosiak, W.J. Structural basis of microRNA length variety. Nucleic Acids Res. 2011, 39, 257–268. [Google Scholar] [CrossRef] [Green Version]
  47. Friedrich, M.; Vaxevanis, C.K.; Biehl, K.; Mueller, A.; Seliger, B. Targeting the coding sequence: Opposing roles in regulating classical and non-classical MHC class I molecules by miR-16 and miR-744. J. Immunother. Cancer 2020, 8. [Google Scholar] [CrossRef] [PubMed]
  48. Jasinski-Bergner, S.; Mandelboim, O.; Seliger, B. The role of microRNAs in the control of innate immune response in cancer. J. Natl. Cancer Inst. 2014, 106. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Jasinski-Bergner, S.; Mandelboim, O.; Seliger, B. Molecular mechanisms of human herpes viruses inferring with host immune surveillance. J. Immunother. Cancer 2020, 8. [Google Scholar] [CrossRef]
  50. Fattore, L.; Costantini, S.; Malpicci, D.; Ruggiero, C.F.; Ascierto, P.A.; Croce, C.M.; Mancini, R.; Ciliberto, G. MicroRNAs in melanoma development and resistance to target therapy. Oncotarget 2017, 8, 22262–22278. [Google Scholar] [CrossRef] [Green Version]
  51. Kurimoto, R.; Chiba, T.; Ito, Y.; Matsushima, T.; Yano, Y.; Miyata, K.; Yashiro, Y.; Suzuki, T.; Tomita, K.; Asahara, H. The tRNA pseudouridine synthase TruB1 regulates the maturation of let-7 miRNA. EMBO J. 2020, 39, e104708. [Google Scholar] [CrossRef]
  52. Penzo, M.; Montanaro, L. Turning Uridines around: Role of rRNA Pseudouridylation in Ribosome Biogenesis and Ribosomal Function. Biomolecules 2018, 8, 38. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Hirano, M.; Das, S.; Guo, P.; Cooper, M.D. The evolution of adaptive immunity in vertebrates. Adv. Immunol. 2011, 109, 125–157. [Google Scholar] [CrossRef]
  54. Bukur, J.; Jasinski, S.; Seliger, B. The role of classical and non-classical HLA class I antigens in human tumors. Semin. Cancer Biol. 2012, 22, 350–358. [Google Scholar] [CrossRef] [PubMed]
  55. Motorin, Y.; Marchand, V. Detection and Analysis of RNA Ribose 2′-O-Methylations: Challenges and Solutions. Genes 2018, 9, 642. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Abe, M.; Naqvi, A.; Hendriks, G.J.; Feltzin, V.; Zhu, Y.; Grigoriev, A.; Bonini, N.M. Impact of age-associated increase in 2′-O-methylation of miRNAs on aging and neurodegeneration in Drosophila. Genes Dev. 2014, 28, 44–57. [Google Scholar] [CrossRef] [Green Version]
  57. Liang, H.; Jiao, Z.; Rong, W.; Qu, S.; Liao, Z.; Sun, X.; Wei, Y.; Zhao, Q.; Wang, J.; Liu, Y.; et al. 3′-Terminal 2′-O-methylation of lung cancer miR-21-5p enhances its stability and association with Argonaute 2. Nucleic Acids Res. 2020, 48, 7027–7040. [Google Scholar] [CrossRef]
  58. Ponten, F.; Jirstrom, K.; Uhlen, M. The Human Protein Atlas--a tool for pathology. J. Pathol. 2008, 216, 387–393. [Google Scholar] [CrossRef]
  59. Thul, P.J.; Akesson, L.; Wiking, M.; Mahdessian, D.; Geladaki, A.; Ait Blal, H.; Alm, T.; Asplund, A.; Bjork, L.; Breckels, L.M.; et al. A subcellular map of the human proteome. Science 2017, 356. [Google Scholar] [CrossRef]
  60. Uhlen, M.; Fagerberg, L.; Hallstrom, B.M.; Lindskog, C.; Oksvold, P.; Mardinoglu, A.; Sivertsson, A.; Kampf, C.; Sjostedt, E.; Asplund, A.; et al. Proteomics. Tissue-based map of the human proteome. Science 2015, 347, 1260419. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Molecular RNA modifications on RNAs involving snoRNAs. (A) 2′-O-methylation and (B) pseudouridylation.
Figure 1. Molecular RNA modifications on RNAs involving snoRNAs. (A) 2′-O-methylation and (B) pseudouridylation.
Cancers 13 01167 g001
Figure 2. Representative immunohistochemical stainings of selected factors involved in 2′-O-methylation and pseudouridylation in human skin and melanoma specimen.
Figure 2. Representative immunohistochemical stainings of selected factors involved in 2′-O-methylation and pseudouridylation in human skin and melanoma specimen.
Cancers 13 01167 g002
Figure 3. Summarizing visualization as a heatmap of the positive and negative correlations of factors involved in 2′-O-methylation and pseudouridylation with melanoma relevant markers of cell proliferation, prognostic markers, and genes involved in immune surveillance as well as immune evasion.
Figure 3. Summarizing visualization as a heatmap of the positive and negative correlations of factors involved in 2′-O-methylation and pseudouridylation with melanoma relevant markers of cell proliferation, prognostic markers, and genes involved in immune surveillance as well as immune evasion.
Cancers 13 01167 g003
Figure 4. Impact of factors involved in 2′-O-methylation and pseudouridylation on the overall survival in MM. (AH): Kaplan–Meier Plots were generated with TCGA data of 214 melanoma specimens provided by the R2 database (https://hgserver1.amc.nl/ (accessed on 12 February 2021)).
Figure 4. Impact of factors involved in 2′-O-methylation and pseudouridylation on the overall survival in MM. (AH): Kaplan–Meier Plots were generated with TCGA data of 214 melanoma specimens provided by the R2 database (https://hgserver1.amc.nl/ (accessed on 12 February 2021)).
Cancers 13 01167 g004
Figure 5. Impact of the RNA modifications analyzed on tumor relevant processes.
Figure 5. Impact of the RNA modifications analyzed on tumor relevant processes.
Cancers 13 01167 g005
Table 1. Localization of the proteins involved into 2′-O-methylation and pseudouridylation in skin and MM specimen based upon the information provided by The Human Protein Atlas.
Table 1. Localization of the proteins involved into 2′-O-methylation and pseudouridylation in skin and MM specimen based upon the information provided by The Human Protein Atlas.
RNA ModificationFactorLocalization within SkinLocalization within MM
factors involved into 2′-O-methylationFBLnuclear, extra strong within the cells of the epidermis especially in the stratum germinativumnuclear
NOP56nuclearnuclear
NOP58nuclear strongest expression within the cells of the epidermal layers, especially within the stratum germinativumcytoplasmic and nuclear
15.5K (SNU13)cytoplasmic and membranous restricted to the cells of the epidermisnuclear
factors involved into pseudouridylationDKC1nuclearnuclear
GAR1nuclearnuclear
NHP2nuclearnuclear
NHP10nuclearnuclear
Table 2. Correlation of the factors involved in 2′-O-methylation and pseudouridylation with known proliferation markers in MM.
Table 2. Correlation of the factors involved in 2′-O-methylation and pseudouridylation with known proliferation markers in MM.
Correlated ExpressionFBLNOP56NOP58SNU13DKC1GAR1NHP2NOP10
MKI67R = 0.081R = 0.216R = 0.235R = 0.050R = 0.185R = 0.154R = −0.133R = 0.011
p = 0.241p = 1.47 × 10−3p = 5.22 × 10−4p = 0.468p = 6.50 × 10−3p = 0.024p = 0.053p = 0.874
PCNAR = 0.152R = 0.581R = 0.448R = 0.106R = 0.463R = 0.238R = 0.372R = 0.258
p = 0.026p = 9.59 × 10−21p = 5.58 × 10−12p = 0.122p = 9.06 × 10−13p = 4.55 × 10−4p = 1.92 × 10−8p = 1.36 × 10−4
CCNA1R = −0.046R = 0.065R = 0.131R = −0.079R = 0.060R = −0.063R = 0.066R = 0.172
p = 0.499p = 0.341p = 0.056p = 0.250p = 0.383p = 0.356p = 0.339p = 0.011
CCNB1R = 0.265R = 0.439R = 0.492R = 0.136R = 0.398R = 0.423R = 0.463R = 0.165
p = 8.62 × 10−5p = 1.78 × 10−11p = 1.88 × 10−14p = 0.047p = 1.50 × 10−9p = 1.10 × 10−10p = 9.44 × 10−13p = 0.015
MCM2R = 0.236R = 0.251R = 0.302R = 0.166R = 0.325R = 0.233R = 0.123R = 0.031
p = 4.94 × 10−4p = 2.07 × 10−4p = 6.94 × 10−6p = 0.015p = 1.16 × 10−6p = 5.80 × 10−4p = 0.072p = 0.656
MCM4R = 0.334R = 0.546R = 0.466R = 0.262R = 0.554R = 0.482R = 0.377R = 0.057
p = 5.60 × 10−7p = 4.78 × 10−18p = 6.15 × 10−13p = 1.04 × 10−4p = 1.34 × 10−18p = 7.93 × 10−14p = 1.29 × 10−8p = 0.403
CENPFR = 0.363R = 0.428R = 0.495R = 0.162R = 0.509R = 0.515R = 0.198R = 0.040
p = 4.69 × 10−8p = 6.01 × 10−11p = 1.33 × 10−14p = 0.018p = 1.62 × 10−15p = 6.59 × 10−16p = 3.61 × 10−3p = 0.556
Based on TCGA data of 214 human MM samples provided by the R2 database (https://hgserver1.amc.nl/ (accessed on 12 February 2021)), statistically significant positive correlations (p < 0.05) of proteins in 2′-O-methylation and pseudouridylation to proliferation markers of MM were found and highlighted in red. Background color: One group of factors belong to the brighter grey. The 2nd group to the darker grey.
Table 3. Correlation of the factors involved in 2′-O-methylation and pseudouridylation with known prognostic markers in MM.
Table 3. Correlation of the factors involved in 2′-O-methylation and pseudouridylation with known prognostic markers in MM.
Correlated ExpressionFBLNOP56NOP58SNU13DKC1GAR1NHP2NOP10
MART1R = 0.009 R = 0.304 R = −0.024 R = 0.134 R = 0.265 R = 0.203 R = 0.259 R = 0.115
p = 0.893p = 6.12 × 10−6p = 0.728p = 0.050p = 8.59 × 10−5p = 2.92 × 10−3p = 1.25 × 10−4p = 0.094
S100BR = −0.248 R = −0.038 R = −0.152 R = −0.085 R = −0.090 R = −0.067 R = 0.133 R = 0.021
p = 2.53 × 10−4p = 0.577p = 0.026p = 0.218p = 0.191p = 0.327p = 0.053p = 0.757
S100A4R = −0.325 R = −0.414 R = −0.223 R = −0.240 R = −0.342 R = −0.274 R = −0.258 R = 0.042
p = 1.20 × 10−6p = 2.79 × 10−10p = 1.01 × 10−3p = 3.96 × 10−4p = 2.98 × 10−7p = 4.75 × 10−5p = 1.33 × 10−4p = 0.539
S100A9R = −0.234 R = −0.192 R = −0.212 R = −0.185 R = −0.295R = −0.234R = −0.147 R = 0.104
p = 5.57 × 10−4p = 4.92 × 10−3p = 1.81 × 10−3p = 6.72 × 10−3p = 1.13 × 10−5p = 5.58 × 10−4p = 0.031p = 0.130
MITFR = 0.076 R = 0.388 R = 0.050 R = 0.124 R = 0.345 R = 0.256 R = 0.304 R = 0.124
p = 0.271p = 4.05 × 10−9p = 0.468p = 0.070p = 2.24 × 10−7p = 1.51 × 10−4p = 5.90 × 10−6p = 0.070
MMP2R = −0.129R = −0.321 R = −0.271 R = −0.190 R = −0.226 R = −0.040 R = −0.309 R = −0.050
p = 0.059p = 1.67 × 10−6p = 5.96 × 10−5p = 5.20 × 10−3p = 8.54 × 10−4p = 0.561p = 4.10 × 10−6p = 0.465
NM23R = −0.005 R = 0.534 R = 0.424 R = 0.053 R = 0.442 R = 0.425 R = 0.531 R = 0.321
p = 0.944p = 3.56 × 10−17p = 9.77 × 10−11p = 0.445p = 1.21 × 10−11p = 8.80 × 10−11p = 5.60 × 10−17p = 1.66 × 10−6
CD44R = −0.091 R = −0.026 R = −0.174 R = 0.034 R = 0.101 R = 0.012 R = −0.040 R = −0.040
p = 0.184p = 0.705p = 0.011p = 0.618p = 0.139p = 0.863p = 0.557p = 0.556
PMELR = −0.033 R = 0.281 R = −0.061 R = 0.029 R = 0.171 R = 0.126 R = 0.160 R = 0.179
p = 0.630p = 2.96 × 10−5p = 0.371p = 0.670p = 0.012p = 0.066p = 0.019p = 8.58 × 10−3
BCL2R = 0.217 R = −0.059 R = −0.072 R = 0.183 R = 0.008 R = 0.031 R = −0.101 R = −0.235
p = 1.40 × 10−3p = 0.392p = 0.292p = 7.21 × 10−3p = 0.907p = 0.652p = 0.140p = 5.15 × 10−4
Based on TCGA data of 214 MM samples provided by the R2 database (https://hgserver1.amc.nl/ (accessed on 12 February 2021)) statistically significant positive (red) or negative (green) correlations (p < 0.05) of proteins in 2′-O-methylation and pseudouridylation to prognostic markers of MM are highlighted. Background color: One group of factors belong to the brighter grey. The 2nd group to the darker grey.
Table 4. Correlation of the expression of factors involved in 2′-O-methylation and pseudouridylation with preselected molecules involved into immune recognition/evasion.
Table 4. Correlation of the expression of factors involved in 2′-O-methylation and pseudouridylation with preselected molecules involved into immune recognition/evasion.
Correlated ExpressionFBLNOP56NOP58SNU13DKC1GAR1NHP2NOP10
(NHP2L1)
molecules contributing to recognition by immune effector cellsMICAR = −0.341 R = 0.100 R = −0.152 R = −0.128 R = −0.112 R = −0.189 R = 0.121 R = 0.245
p = 3.10 × 10−7p = 0.146p = 0.026p = 0.061p = 0.103p = 5.59 × 10−3p = 0.079p = 2.88 × 10−4
MICBR = −0.241 R = −0.119 R = −0.069 R = −0.014 R = −0.149 R = −0.304 R = −0.106 R = 0.175
p = 3.72 × 10−4p = 0.082p = 0.314p = 0.834p = 0.029p = 5.75 × 10−6p = 0.121p = 0.010
ULBP1R = 0.109 R = 0.020 R = 0.109 R = 0.011 R = 0.067 R = 0.056 R = 0.077 R = 0.029
p = 0.112p = 0.767p = 0.112p = 0.877p = 0.328p = 0.412p = 0.264p = 0.677
ULBP2R = 0.047 R = 0.078 R = 0.107 R = −0.020 R = 0.096 R = 0.120 R = 0.093 R = −0.023
p = 0.496p = 0.255p = 0.117p = 0.766p = 0.162p = 0.080p = 0.176p = 0.742
ULBP3R = −0.091 R = 0.028 R = 0.005 R = −0.008 R = 0.036 R = 0.058 R = 0.060 R = 0.098
p = 0.185p = 0.688p = 0.947p = 0.912p = 0.603p = 0.398p = 0.380p = 0.154
ULBP4R = 0.048 R = −0.104 R = −0.171 R = 0.202 R = −0.015 R = −0.008 R = −0.024 R = −0.107
p = 0.488p = 0.129p = 0.012p = 2.93 × 10−3p = 0.825p = 0.906p = 0.729p = 0.117
ULBP5R = −0.011 R = −0.007 R = 0.018 R = 0.022 R = −0.001 R = 0.069 R = −0.022 R = 0.024
p = 0.867p = 0.915p = 0.791p = 0.744p = 1.000p = 0.313p = 0.752p = 0.730
ULBP6R = −0.005 R = 0.020 R = −0.132 R = −0.029 R = −0.093 R = −0.138 R = 0.025 R = −0.004
p = 0.940p = 0.772p = 0.055p = 0.674p = 0.175p = 0.044p = 0.715p = 0.951
HLA-AR = −0.364 R = −0.261 R = −0.370 R = −0.054 R = −0.195 R = −0.238R = −0.229 R = 0.110
p = 4.16 × 10−8p = 1.13 × 10−4p = 2.47 × 10−8p = 0.431p = 4.12 × 10−3p = 4.55 × 10−4p = 7.55 × 10−4p = 0.109
HLA-BR = −0.369 R = −0.353 R = −0.342 R = −0.024 R = −0.305 R = −0.365 R = −0.218 R = 0.114
p = 2.66 × 10−8p = 1.16 × 10−7p = 2.81 × 10−7p = 0.723p = 5.60 × 10−6p = 3.92 × 10−8p = 1.30 × 10−3p = 0.096
HLA-CR = −0.351 R = −0.264 R = −0.323 R = −0.126 R = −0.298 R = −0.360 R = −0.048 R = 0.145
p = 1.29 × 10−7p = 9.52 × 10−5p = 1.35 × 10−6p = 0.066p = 9.08 × 10−6p = 5.96 × 10−8p = 0.481p = 0.034
B2MR = −0.407R = −0.281 R = −0.224 R = −0.150 R = −0.198 R = −0.318 R = −0.178 R = 0.314
p = 6.26 × 10−10p = 2.97 × 10−5p = 9.60 × 10−4p = 0.028p = 3.67 × 10−3p = 2.04 × 10−6p = 9.12 × 10−3p = 2.73 × 10−6
TAP1R = −0.251 R = −0.074 R = −0.176 R = 0.067 R = −0.114 R = −0.154 R = 0.018 R = 0.189
p = 2.05 × 10−4p = 0.278p = 9.93 × 10−3p = 0.331p = 0.097p = 0.024p = 0.788p = 5.58 × 10−3
TAP2R = −0.102 R = −0.042 R = −0.064 R = 0.144 R = 0.027 R = −0.031 R = −0.081 R = 0.034
p = 0.137p = 0.538p = 0.355p = 0.035p = 0.695p = 0.655p = 0.241p = 0.624
TAPBPR = −0.210 R = −0.309 R = −0.374 R = 0.070 R = −0.156 R = −0.219 R = −0.328 R = −0.088
p = 1.99 × 10−3p = 3.94 × 10−6p = 1.70 × 10−8p = 0.311p = 0.022p = 1.25 × 10−3p = 9.18 × 10−7p = 0.19
LMP2R = −0.342 R = −0.214 R = −0.209 R = −0.038 R = −0.262 R = −0.338 R = −0.014 R = 0.190
p = 2.96 × 10−7p = 1.60 × 10−3p = 2.07 × 10−3p = 0.577p = 1.03 × 10−4p = 4.19 × 10−7p = 0.833p = 5.41 × 10−3
LMP7R = −0.299 R = −0.159 R = −0.233 R = 0.055 R = −0.164 R = −0.180 R = 0.092 R = 0.165
p = 8.42 × 10−6p = 0.020p = 5.83 × 10−4p = 0.420p = 0.016p = 8.24 × 10−3p = 0.179p = 0.015
LMP10R = −0.092 R = −0.148 R = −0.205 R = 0.008 R = −0.180 R = −0.245 R = −0.057 R = 0.138
p = 0.180p = 0.030p = 2.55 × 10−3p = 0.909p = 8.33 × 10−3p = 2.99 × 10−4p = 0.406p = 0.044
ERAAPR = −0.208 R = −0.141 R = −0.122 R = −0.104 R = −0.161 R = −0.197 R = −0.006 R = 0.036
p = 2.24 × 10−3p = 0.040p = 0.075p = 0.131p = 0.019p = 3.73 × 10−3p = 0.932p = 0.602
ERP57R = 0.004 R = 0.077 R = 0.014 R = −0.031 R = 0.081 R = −0.017 R = 0.012 R = −0.073
p = 0.951p = 0.262p = 0.841p = 0.648p = 0.239p = 0.804p = 0.859p = 0.285
CALRR = −0.169 R = 0.073 R = −0.102 R = −0.241 R = −0.010 R = −0.072 R = −0.191 R = 0.204
p = 0.013p = 0.288p = 0.138p = 3.64 × 10−4p = 0.888p = 0.296p = 5.11 × 10−3p = 2.74 × 10−3
CANXR = −0.303 R = −0.430 R = −0.332 R = −0.147 R = −0.300 R = −0.370 R = −0.320 R = −0.022
p = 6.34 × 10−6p = 4.64 × 10−11p = 6.62 × 10−7p = 0.031p = 8.14 × 10−6p = 2.48 × 10−8p = 1.79 × 10−6p = 0.744
molecules contributing to immune evasionPD-L1
(B7-H1)
R = −0.107 R = −0.051 R = −0.004 R = −0.103 R = 0.010 R = −0.062 R = 0.003 R = 0.041
p = 0.117p = 0.460p = 0.954p = 0.132p = 0.888p = 0.363p = 0.966p = 0.548
PD-L2
(PDCD1LG2)
R = −0.151 R = −0.150 R = −0.092 R = −0.096 R = −0.114 R = −0.215 R = 0.014 R = 0.068
p = 0.027p = 0.028p = 0.182p = 0.163p = 0.097p = 1.59 × 10−3p = 0.838p = 0.319
HLA-GR = −0.230 R = −0.190 R = −0.345 R = −0.008 R = −0.217 R = −0.280 R = −0.111 R = 0.081
p = 6.90 × 10−4p = 5.32 × 10−3p = 2.20 × 10−7p = 0.905p = 1.39 × 10−3p = 3.19 × 10−5p = 0.104p = 0.236
HLA-ER = −0.321 R = −0.290 R = −0.418 R = −0.045 R = −0.296 R = −0.380 R = −0.236 R = −0.013
p = 1.59 × 10−6p = 1.65 × 10−5p = 1.95 × 10−10p = 0.509p = 1.06 × 10−5p = 9.50 × 10−9p = 4.92 × 10−4p = 0.845
HLA-FR = −0.323 R = −0.274 R = −0.329 R = −0.005 R = −0.267 R = −0.318 R = −0.182 R = 0.098
p = 1.41 × 10−6p = 4.71 × 10−5p = 8.52 × 10−7p = 0.944p = 7.79 × 10−5p = 1.98 × 10−6p = 7.47 × 10−3p = 0.154
CD155
(PVR)
R = 0.029 R = 0.436 R = 0.098 R = 0.011 R = 0.227 R = 0.263 R = 0.125 R = 0.148
p = 0.678p = 2.41 × 10−11p = 0.151p = 0.874p = 8.28 × 10−4p = 1.00 × 10−4p = 0.067p = 0.030
B7-H4
(VTCN1)
R = −0.003 R = −0.059 R = −0.016 R = −0.048 R = −0.031 R = 0.034 R = 0.062 R = −0.086
p = 0.963p = 0.388p = 0.814p = 0.482p = 0.657p = 0.624p = 0.369p = 0.210
Bioinformatics analysis of TCGA data of 214 melanoma specimens provided by the R2 database (https://hgserver1.amc.nl/ (accessed on 12 February 2021)) was performed. Statistically significant positive (red) or negative (green) correlations (p < 0.05) are highlighted. Background color: One group of factors belong to the brighter grey. The 2nd group to the darker grey.
Table 5. Correlation of miR expression with enzymes known to be involved in 2′-O-methylation or pseudouridylation of miRs.
Table 5. Correlation of miR expression with enzymes known to be involved in 2′-O-methylation or pseudouridylation of miRs.
Induced miRs in MM with Diagnostic/Prognostic
Relevance
So Far in Literature Mentioned Enzymes with
Putative Role for 2′-O-Methylation of miRs
So Far in Literature Mentioned Enzymes with Putative Role for Pseudouridylation of miRs
FBLHENMT1DKC1TRUB1
miR-9-5pR = 0.007 R = 0.050 R = −0.026 R = 0.132
p = 0.916p = 0.469p = 0.708p = 0.053
miR-10aR = 0.020 R = −0.036 R = −0.004 R = −0.047
p = 0.767p = 0.603p = 0.957p = 0.492
miR-10bR = −0.012 R = 0.009 R = −0.064 R = −0.200
p = 0.857p = 0.896p = 0.351p = 3.33 × 10−3
miR-17-5pn.d.n.d.n.d.n.d.
miR-18an.d.n.d.n.d.n.d.
miR-21R = −0.168 R = 0.017 R = −0.304 R = −0.139
p = 0.014p = 0.800p = 5.88 × 10−6p = 0.042
miR-26bR = 0.191 R = −0.008 R = 0.112 R = −0.091
p = 4.96 × 10−3p = 0.906p = 0.101p = 0.184
miR-92aR = 0.106 R = −0.055 R = 0.155 R = 0.059
p = 0.122p = 0.424p = 0.023p = 0.390
miR-221R = −0.029 R = 0.100 R = 0.016 R = 0.064
p = 0.676p = 0.144p = 0.815p = 0.352
miR-222R = −0.015 R = 0.083 R = 0.083 R = −0.073
p = 0.825p = 0.227p = 0.224p = 0.285
miR-126R = 0.088 R = −0.025 R = −0.034 R = −0.091
p = 0.198p = 0.711p = 0.620p = 0.187
miR-145n.d.n.d.n.d.n.d.
miR146aR = −0.015 R = −0.081 R = 0.003 R = 0.098
p = 0.822p = 0.235p = 0.968p = 0.154
miR-182R = 0.008 R = −0.068 R = −0.045 R = 0.071
p = 0.905p = 0.322p = 0.510p = 0.299
miR-514R = −0.110 R = 0.065 R = 0.011 R = 0.052
p = 0.109p = 0.340p = 0.874p = 0.446
miR-520dn.d.n.d.n.d.n.d.
miR-527n.d.n.d.n.d.n.d.
Statistically significant positive or negative correlations (p < 0.05) are highlighted in red or green. Background color: One group of factors belong to the brighter grey. The 2nd group to the darker grey.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Jasinski-Bergner, S.; Blümke, J.; Wickenhauser, C.; Seliger, B. Relevance of 2′-O-Methylation and Pseudouridylation for the Malignant Melanoma. Cancers 2021, 13, 1167. https://doi.org/10.3390/cancers13051167

AMA Style

Jasinski-Bergner S, Blümke J, Wickenhauser C, Seliger B. Relevance of 2′-O-Methylation and Pseudouridylation for the Malignant Melanoma. Cancers. 2021; 13(5):1167. https://doi.org/10.3390/cancers13051167

Chicago/Turabian Style

Jasinski-Bergner, Simon, Juliane Blümke, Claudia Wickenhauser, and Barbara Seliger. 2021. "Relevance of 2′-O-Methylation and Pseudouridylation for the Malignant Melanoma" Cancers 13, no. 5: 1167. https://doi.org/10.3390/cancers13051167

APA Style

Jasinski-Bergner, S., Blümke, J., Wickenhauser, C., & Seliger, B. (2021). Relevance of 2′-O-Methylation and Pseudouridylation for the Malignant Melanoma. Cancers, 13(5), 1167. https://doi.org/10.3390/cancers13051167

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop