Abstract

This study aimed to investigate the prognostic value of HRAS mRNA expression in cutaneous melanoma. Cutaneous melanoma is an aggressive cancer with an increasing incidence. Few studies have focused on the transcriptional level of RAS isoforms (KRAS, NRAS, and HRAS) in cutaneous melanoma. To gain further insight into RAS isoforms at transcriptional level, we obtained the cutaneous melanoma data from cBioPortal and investigated the RAS mRNA expression levels in different stages of melanoma and evaluated their correlation with clinical characteristics and patients’ survival. Furthermore, we retrieved and analyzed the coexpression data and performed pathway enrichment analysis. Totally, 452 cutaneous melanoma cases were included in this study. We found that lower HRAS expression level was associated with longer patient survival. 206 genes that negatively correlated with HRAS expression were positively correlated with KRAS and NRAS expression. In contrast, no gene that positively correlated with HRAS expression was positively correlated with KRAS and NRAS expression. In conclusion, our data showed that transcriptional regulation was different for the three RAS isoforms in cutaneous melanoma. This study highlighted the prognostic value of HRAS mRNA expression and revealed that HRAS greatly differs from KRAS and NRAS at the transcriptional level.

1. Introduction

RAS proteins are small GTPases that signal downstream of cell surface receptors to regulate cell proliferation, differentiation, migration, and apoptosis [13]. In human, there are 3 ubiquitously expressed RAS genes, KRAS, NRAS, and HRAS, which encode four highly homologous RAS proteins (splicing variants KRAS4A/4B, NRAS, and HRAS). These RAS genes are the most frequently mutated genes in human cancers [4]. Mutations in the RAS genes are associated with around 30% of all human tumors [5]. In a specific type of cancer, usually only a single RAS isoform is mutated. For example, mutations in KRAS are common in lung, colon, and pancreatic cancers, those in NRAS predominate in melanoma, and HRAS mutations are commonly seen in bladder, head and neck, and skin cancers [69].

Cutaneous melanoma is an aggressive cancer with an increasing incidence [10]. Alterations in multiple signaling pathways that regulate cell proliferation and survival contribute to the tumorigenesis and progression of this disease [11]. Cutaneous melanoma frequently harbors activating mutations in NRAS (around 20%) or the RAS-regulated kinase BRAF (around 37%), suggesting that the RAS-RAF-MAPK pathway may be critical in the pathogenesis of cutaneous melanoma [1214]. Although NRAS shares high sequence homology with KRAS and HRAS and all these three RAS isoforms can be signaled through the RAS-RAF-MAPK pathway, HRAS and KRAS mutations are rare in this disease [15, 16]. Different functions of these three RAS isoforms have been reported in many literatures [1, 1719]. However, differential regulating of these RAS isoforms at transcriptional level was not well studied.

Currently, there has been remarkable progress in understanding melanoma pathogenesis, and numerous studies have now shown the correlation of BRAF and NRAS mutation status with clinical outcome and immune and targeted therapy strategies in melanoma [13, 20]. However, few studies have focused on the characteristics of RAS isoforms mRNA expression in cutaneous melanoma.

In this study, we sought to examine the differential regulation of RAS isoforms at mRNA level in cutaneous melanoma. We analyzed the correlation between RAS mRNA expression and different clinical characteristics and explored the genetic mutation data. By analyzing genes coexpressed with different RAS isoforms, we found significant difference between genes coexpressed with HRAS and those with KRAS/NRAS. We also analyzed the correlation of RAS mRNA expression with BRAF. Furthermore, we performed pathway enrichment analysis, which provided more clues about the differential regulation of RAS isoforms.

2. Materials and Methods

2.1. Data Sources and Bioinformatics

We obtained the cutaneous melanoma data from cBioPortal (http://www.cbioportal.org). The data of skin cutaneous melanoma (TCGA, Provisional) included mRNA expression data (RNA Seq V2 RSEM) of the three RAS isoforms (KRAS, NRAS, and HRAS), matching clinical information, genetic mutation data, and genes coexpressed with the three RAS isoforms. The cases selection included three criteria: (i) each patient ID corresponded to a unique sample ID; (ii) the case included available RAS isoforms (KRAS, NRAS, and HRAS) mRNA expression data; and (iii) the case included available age and gender information.

The following variables were obtained from cBioPortal (http://www.cbioportal.org): age, gender, AJCC stage, tumor thickness, ulceration, overall survival status, genetic mutation data, and genes coexpressed with the three RAS isoforms. To ensure accuracy, dual data extraction was conducted. Data were subsequently verified between the 2 reviewers, and discrepancies were settled through consensus discussion. To minimize subjective judgment and selection bias, investigators were blinded to outcomes.

2.2. Clinical Characteristics and Survival Analysis

We compared mRNA expression levels of the three RAS isoforms in different clinical subgroups. Clinical factors considered included age, gender, AJCC stage, tumor thickness (mm), and ulceration. Cases were divided into “low” and “high” groups based on the median RAS expression level.

For survival analysis, based on the median levels of RAS mRNA expression, cases were grouped into “low” and “high.” Because BRAF is one of the most important genes in melanoma, we also investigated the correlation of BRAF mRNA expression and overall survival. Cases were also divided into “low” and “high” groups based on the median BRAF expression level.

Because LCK was an important prognostic factor in cutaneous melanoma [21], we also analyzed the correlation between LCK and HRAS at transcriptional level.

2.3. mRNA Expression of RAS Isoforms according to Genetic Mutation Status

We compared RAS mRNA expression levels according to RAS mutation status and BRAF mutation status, as BRAF is one of the most frequently mutated genes in cutaneous melanoma. Additionally, we analyzed the correlation of RAS mRNA expression and BRAF mRNA expression.

2.4. Coexpression and Pathway Enrichment Analysis

We retrieved and analyzed the coexpression data using the coexpression tool in cBioPortal. Genes with Pearson’s correlation coefficients (CC) ≥ 0.3 were considered positively correlated with RAS expression. Genes with Pearson’s correlation coefficients (CC) ≤ −0.3 were considered negatively correlated with RAS expression. For the genes positively and negatively correlated with KRAS, NRAS, and HRAS expression, we analyzed the number of overlapping genes and enriched pathways.

The pathway enrichment analysis was carried out using the default parameters of Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway by the Database for Annotation Visualization and Integrated Discovery (DAVID) (https://david.ncifcrf.gov/). was used as the cut-off.

2.5. Statistical Analysis

The correlations between RAS mRNA expression and clinical characteristics including age, gender, AJCC stage, tumor thickness, and ulceration were analyzed by the Chi-square test. Overall survival was assessed by Kaplan-Meier analysis, and univariate and multivariate Cox regression analyses were used to calculate hazard ratios (HRs) and their 95% confidence intervals (95% CI). The RAS mRNA expression difference between wild-type and mutant groups was evaluated by Mann–Whitney test. In addition, Pearson’s correlation coefficients analysis was used for the correlation of LCK and HRAS, RAS isoforms, and BRAF at transcriptional level, respectively. All tests were 2-sided, and was considered statistically significant. All analysis was conducted using SPSS 17.0 (IBM Corporation, Armonk, NY, USA).

3. Results

3.1. Clinicopathological and Demographic Characteristics of Patients

A total 452 cutaneous melanoma patients were included in this study, with 280 males and 172 females. All cases had available age, gender, and mRNA expression data for the three RAS isoforms by RNA Seq V2 RSEM, and each patient ID corresponded to a unique sample ID. The data were obtained on August 22, 2016.

The median age was 58 years (range, 15~90 years), and the medium overall survival follow-up time was 36.93 months (range, 0.20~369.65 months). The baseline clinicopathological and demographic characteristics of patients were summarized in Table 1. A total of 216 patients (47.8%) died during the period of follow-up.

3.2. The Relationship between mRNA Expression of RAS Isoforms and Clinical Characteristics

The correlation between RAS mRNA expression and clinical characteristics was summarized in Table 2. There was significant correlation between HRAS mRNA expression and tumor thickness; higher HRAS mRNA level was correlated with higher tumor thickness (). There was no correlation between HRAS mRNA expression and age, gender, AJCC stage, and ulceration, respectively. KRAS () and NRAS () mRNA expressions were positively correlated with age, respectively. No other significant correlations between KRAS or NRAS mRNA expression and gender, AJCC stage, tumor thickness, or ulceration were found (Table 2).

3.3. Prognostic Value of HRAS mRNA Expression in Cutaneous Melanoma

Interestingly, lower HRAS expression level was associated with longer overall survival (Figure 1(c)). Cases were classified into “low” and “high” groups using median HRAS mRNA expression as cut-off. For total 452 cases, time to death was plotted in a Kaplan-Meier curve for those cases exhibiting HRAS mRNA expression above the median () or below the median () level of expression. As shown in Figure 1(c), the survival curves were significantly different (); those cases had lower mRNA expression of HRAS, surviving longer (Figure 1(c)). Furthermore, our results showed that there was a negative correlation between LCK and HRAS at transcriptional level (CC = −0.122; ).

Similar overall survival analyses were carried out for KRAS, NRAS, and BRAF mRNA expression. There were no significant differences in the relationships of KRAS (, Figure 1(a)), NRAS (, Figure 1(b)), and BRAF (, Figure 1(d)) mRNA expression and overall survival when total cases () were considered, respectively.

In addition, we used a Cox proportional hazard regression model to estimate the crude HRs of each clinicopathological characteristic (Table 3). Univariate Cox regression survival analysis showed that high HRAS mRNA expression was a risk factor for worse patient survival (HR: 1.532, 95% CI: 1.163~2.017, for overall survival). And univariate analysis revealed that age (), AJCC stage (), tumor thickness (), and ulceration () were all significantly associated with overall survival (Table 3). However, KRAS (), NRAS (), and BRAF () mRNA expressions did not correlate with patient survival (Table 3). Furthermore, multivariate Cox regression analyses on overall survival showed that HRAS mRNA expression (HR: 1.555, 95% CI: 1.066~2.269, ), AJCC stage (HR: 1.921, 95% CI: 1.336~2.762, ), and ulceration (HR: 1.787, 95% CI: 1.229~2.600, ) were independent prognostic factors.

3.4. mRNA Expression of RAS Isoforms according to Genetic Mutation Status

In 452 patients, 182 (40.3%) melanomas were BRAF-mutant, 95 (21.0%) were NRAS-mutant, and 180 (39.8%) were wild-type (no BRAF mutation and no NRAS mutation); only 5 (1.1%) were both BRAF-mutant and NRAS-mutant. There are very few cases of KRAS mutation (7 cases, 1.5%) or HRAS mutation (5 cases, 1.1%) (Table 4).

Compared to wild-type RAS isoforms, mutational RAS isoforms mRNA expressions were significantly higher, respectively ( for KRAS, for NRAS, and for HRAS) (Table 4). For the NRAS-mutant group, HRAS mRNA expression was significantly lower than NRAS wild-type group () (Table 4). For the group with BRAF mutation, NRAS mRNA expression was significantly lower than BRAF wild-type group () (Table 4).

Furthermore, we have analyzed the correlation of RAS mRNA expression and BRAF mRNA expression. The results showed that there was a positive correlation between KRAS and NRAS (CC = 0.250; ), between KRAS and BRAF (CC = 0.210; ), and between NRAS and BRAF (CC = 0.241; ), respectively. Notably, HRAS mRNA expression was negatively correlated with KRAS (CC = −0.192; ), NRAS (CC = −0.281; ), and BRAF (CC = −0.396; ) mRNA expressions, respectively (Figure 2).

3.5. Coexpression and Pathway Enrichment Analysis

To get a better picture about transcriptional regulation of RAS isoforms, we analyzed genes that were coexpressed with RAS isoforms. Genes with Pearson’s correlation coefficient (CC) ≥ 0.3 or ≤−0.3 were evaluated. 355 genes were positively correlated and 31 were negatively correlated with KRAS expression. 1217 genes were positively correlated and 650 were negatively correlated with NRAS expression. There were 1784 genes positively correlated and 1624 genes negatively correlated with HRAS expression (Figure 3).

In addition, we explored the top 10 enriched pathways by KEGG pathway analyses for genes positively (Figure 4(a)) and negatively (Figure 4(b)) correlated with HRAS expression, respectively. The enriched pathways of the genes positively correlated with HRAS included oxidative phosphorylation (hsa00190, ), Huntington’s disease (hsa05016, ), Parkinson’s disease (hsa05012, ), Alzheimer’s disease (hsa05010, ), and ribosome (hsa03010, ) (Figure 4(a)). And the enriched pathways of the genes negatively correlated with HRAS included ubiquitin-mediated proteolysis (hsa04120, ), RNA degradation (hsa03018, ), MAPK signaling pathway (hsa04010, ), Wnt signaling pathway (hsa04310, ), and RIG-I-like receptor signaling pathway (hsa04622, ) (Figure 4(b)).

Interestingly, 206 genes that negatively correlated with HRAS expression were positively correlated with KRAS and NRAS expression. In contrast, no gene that positively correlated with HRAS expression was positively correlated with KRAS and NRAS expression (Figures 3(a) and 3(b)).

To explore possible signal pathways, we studied the enriched pathways of the 206 genes that negatively correlated with HRAS but positively correlated with KRAS and NRAS expression (Table 5). We found that these genes were enriched in the following pathways (): chronic myeloid leukemia (hsa05220, ), prostate cancer (hsa05215, ), regulation of actin cytoskeleton (hsa04810, ), acute myeloid leukemia (hsa05221, ), cell cycle (hsa04110, ), spliceosome (hsa03040, ), pancreatic cancer (hsa05212, ), MAPK signaling pathway (hsa04010, ), and TGF-beta signaling pathway (hsa04350, ) (Table 5).

Similar analyses were carried out for genes positively correlated with HRAS expression and negatively correlated with KRAS/NRAS mRNA expression. There were 31 genes in this group, but no gene that negatively correlated with KRAS, NRAS, and HRAS was found (Figures 3(c) and 3(d)). For the 31 overlapping genes, no pathway was significantly associated with these genes.

4. Discussion

This study highlighted a novel role for the prognostic value of HRAS mRNA expression in cutaneous melanoma. To our knowledge, few studies examined the association of RAS mRNA expression with overall survival. Our data indicated that HRAS mRNA expression, but not KRAS or NRAS, was correlated with prognosis in cutaneous melanoma. Patients with higher HRAS mRNA expression levels in the tumors had poorer overall survival.

It has been reported that HRAS mRNA level has prognostic meaning in triple-negative breast cancers [22]. In our study, we found that HRAS mRNA level was correlated with the poor prognosis of the cutaneous melanoma. We also detected that increased HRAS mRNA level was correlated with high thickness. More mechanistic studies are needed to further elucidate the relationship between HRAS mRNA and the prognosis of cutaneous melanoma [23]. Furthermore, Because LCK was the most important prognostic factor in cutaneous melanoma [21], the negative correlation between HRAS and LCK at transcriptional level indicates that HRAS might be a critical molecule in cutaneous melanoma.

Whether RAS mutation is a prognostic factor in melanoma was controversial. Some studies found no prognostic impact of mutation status [24] and other studies demonstrated that BRAF-mutant or NRAS-mutant melanoma patients had poorer overall survival [2527].

Another interesting observation is that HRAS mRNA expression was negatively correlated with that of KRAS, NRAS, and BRAF, respectively. In contract, KRAS, NRAS, and BRAF mRNA expressions were positively correlated with each other. We also found that there were 206 genes that negatively correlated with HRAS expression significantly overlapped with genes positively correlated with KRAS/NRAS in cutaneous melanoma. Coexpression analysis showed that there was no gene coexpressed with all three isoforms of RAS. This observation suggests that regulation of HRAS in melanoma is different from that of KRAS and NRAS. Although those three RAS isoforms share similar downstream pathways, their transcriptional regulation may be different.

Despite a high degree of sequence homology among the RAS isoforms, different functions and subcellular localization of them have been reported [3]. Different RAS isoforms distinctly contribute to embryonic development, cancer development, cellular homeostasis, and differential coupling to canonical effector pathways [3, 28]. However, the mechanisms underlying the biological differences between the RAS isoforms remain unclear.

Our study suggests that, at transcriptional level, the three RAS genes are differently regulated. RAS isoforms have some distinct transcription factors. Transcription factor binding site analysis in the promoter region of HRAS suggests that p53, STAT3, c-Myc, NF-1, NF-1/L, NF-B, NF-Y, c-Myb, and Max can regulate HRAS transcription. For KRAS, the transcription factors include p53, AP-1, c-Jun, Elk-1, NF-B, STAT1, and PPAR-1. For NRAS, the transcription factors are p53, C/EBP-, c-Fos, AP-1, c-Jun, STAT5A, MyoD, and NF-1 (data obtained from GeneCards: the Human Gene Database; http://www.genecards.org/). There are 47 NF-1 binding sites in the HRAS gene promoter but much less for KRAS and NRAS. In KRAS gene promoter, there are 16 binding sites of transcription factor STAT1, which are the most abundant. In NRAS gene promoter, there are 11 binding sites of AP-1, which is the most abundant transcription factor (data obtained from GeneCards: the Human Gene Database; http://www.genecards.org/). In the pathogenesis of cutaneous melanoma, the differential regulation of the above-mentioned transcription factors may be the primary reason to regulate RAS isoforms mRNA expression.

On average, 16% of human cancers harbor activating mutations of RAS at amino acid 12, 13, or 61 [29]. With regard to NRAS, the most common oncogenic change (>80% of all NRAS mutations) is a point mutation leading to the substitution of leucine by glutamine at amino acid 61, with mutations at amino acids 12 and 13 occurring less frequently [8, 14]. Amino acid 61 mutations also account for the majority of HRAS mutations in melanoma, whereas most KRAS mutations are at amino acid 12 [14]. Although it is not clear why NRAS mutations are more frequent in melanoma compared to HRAS or KRAS mutations, there is evidence that NRAS is overexpressed in melanocytes relative to other RAS isoforms. It is also possible that NRAS may activate different signaling pathways from KRAS and HRAS, an idea supported by the observation that NRAS has greater transforming activity than KRAS in experimental models of melanoma [14, 15, 30].

RAS activating mutations have been viewed as functionally dominant because they are constitutively active. However, the status of the wide-type RAS allele may also play a role in tumors carrying mutant RAS genes [9, 30]. Several recent studies have shown that oncogenic and wild-type RAS play independent and nonredundant roles within the cell. Oncogenic RAS regulates basal effector signaling, whereas wild-type RAS mediates signaling downstream of activated receptor tyrosine kinases (RTKs) [30, 31]. Further studies are needed to determine the association between wild-type RAS and mutants.

5. Conclusions

In conclusion, our data showed that transcriptional regulation was different for the three RAS isoforms in cutaneous melanoma. This study highlighted the prognostic value of HRAS mRNA expression and revealed that HRAS greatly differs from KRAS and NRAS at the transcriptional level. Further understanding of the upstream signaling pathways regulating RAS gene expression will shed light on pathogenesis of cutaneous melanoma.

Conflicts of Interest

The authors declare that they have no actual or potential conflicts of interest.

Authors’ Contributions

Xiaohua Wan and Zhongwu Li carried out the studies, participated in data collection, and drafted the manuscript. Ruping Liu and Zhongwu Li performed the statistical analysis and participated in study design. All authors read and approved the final manuscript.

Acknowledgments

This study was supported by the Beijing Municipal Science and Technology Commission NOVA Program (No. 2010 B033), the Beijing Municipal Science and Technology Commission Capital Characteristic Clinical Application Research (no. Z141107002514077), the National Natural Science Foundation of China (no. 61501039), and Scientific Research Funding from Beijing Cancer Hospital (Approval no. 13-11).