Skip to main content
Download PDF

GRASLND regulates melanoma cell progression by targeting the miR-218-5p/STAM2 axis

Journal of Translational Medicine volume 22, Article number: 684 (2024) Cite this article

Abstract

Background

Increasing evidence suggests that long noncoding RNAs (lncRNAs) play important regulatory roles in biological processes and are dysregulated in numerous tumors. The lncRNA GRASLND functions as an oncogene in many cancers, but its role in skin cutaneous melanoma (SKCM) requires further investigation.

Methods

SiRNA transfection, wound − healing and transwell assays were performed to evaluate the effect of GRASLND on cellular function.

Results

The present study demonstrated that GRASLND expression is increased in SKCM tissues and cell lines. The high expression of GRASLND was correlated with poor prognosis and immunotherapy outcomes. Knockdown of GRASLND significantly inhibited cell migration and invasion. In addition, we found that miR-218-5p directly binds to its binding site on GRASLND, and GRASLND and miR-218-5p demonstrate mutual inhibition. Furthermore, the miR-218-5p inhibitor partially eliminated the knockdown of GRASLND and inhibited its expression. We also demonstrated that GRASLND acts as a miR-218-5p sponge that positively regulates STAM2 expression in SKCM cells.

Conclusion

In summary, these data suggest that GRASLND functions by regulating miR-218-5p/STAM2 expression, suggesting an important role for the lncRNA‒miRNA–mRNA functional network and a new potential therapeutic target for SKCM.

Introduction

Skin melanoma (SKCM) is a common skin tumor causing ~ 55,500 deaths annually [1]. Melanoma accounts for approximately 1.7% of all newly diagnosed primary malignant cancers worldwide, and melanoma-related deaths account for ~ 0.7% of all cancer-related deaths [1]. In recent years, the incidence of SKCM has increased annually, imposing a substantial social and economic burden [2]. The main current treatment for SKCM is surgery, but even with surgical resection, the prognosis remains poor. The 5-year survival rates for patients with limited, regional, and metastatic SKCM were 98%, 64%, and 23%, respectively [3]. This poor survival rate is partially because the disease is challenging to diagnose in the early stage and less curable in the late stage. Therefore, screening for potential molecular diagnostic markers and therapeutic targets is essential for predicting patient prognosis and developing personalized treatments.

Long noncoding RNAs (lncRNAs) are a class of RNA molecules greater than 200 nucleotides in length that are transcribed from genes but not translated into proteins. Initially, lncRNAs were thought to be “junk sequences” that did not perform important functions. Recent studies have shown that lncRNAs are involved in tumorigenesis by maintaining proliferative signals, activating invasion and metastasis, and inducing angiogenesis (Table 1) [4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22]. LncRNAs also play important roles in SKCM. For example, silencing the lncRNA SAMMSON inhibits melanoma survival by disrupting important p32-mediated mitochondrial functions [23]. In addition, the lncRNA SLNCR1 increases melanoma invasion by transcriptionally upregulating matrix metalloproteinase 9 (MMP9) through cooperation with brain-specific homologous box protein 3a (Brn3a) and the androgen receptor (AR) [24]. Iyer et al. [25] identified 339 melanoma-associated lncRNAs; however, the effects of most lncRNAs on melanoma are currently unknown.

Table 1 The role of lncRNAs in various cancers
Full size table

In this study, we screened a series of differentially expressed lncRNAs between SKCM and healthy skin tissues via gene chip and The Cancer Genome Atlas Program (TCGA )RNA-seq data integration analyses, constructed a model via machine learning algorithms, and further downscaled the model to obtain the key lncRNA GRASLND. Subsequently, we performed cellular experiments and utilized multiorganomics data mining to explore the functions, mechanisms, and clinical biomarkers of the lncRNAs. Our findings suggest that GRASLND is a potential marker for prognostic and immunotherapeutic response determination in SKCM patients, and an in-depth study of its functional mechanism may provide new targets and therapeutic strategies for melanoma immunotherapy.

Methods and materials

Tissue sample collection

Melanoma (n = 16) and paraneoplastic normal (n = 12) tissue specimens were obtained from patients with cutaneous melanoma who underwent tumor resection at the Third Xiangya Hospital of Central South University. The histologic type of our selected case was nevoid melanoma. The study was approved by the Ethics Committee of the Third Xiangya Hospital of Central South University, and informed consent was obtained from all the participants.

Cell culture and siRNA transfection

Two melanoma cell lines, A2058 and SK-MEL-28, were purchased from the American Type Culture Collection (ATCC). A2058 and SK-MEL-28 cells were cultured in Dulbecco’s modified Eagle’s medium (DMEM)(Gibco, USA) and Roswell Park Memorial Institute 1640 medium(RPMI-1640)(Gibco, USA), respectively. All cells were cultured in a medium containing 10% fetal bovine serum (FBS) (Gibco, USA) and 1% 100× penicillin-streptomycin (Beyotime Biotechnology, China) and incubated at 37 °C under 5% CO2 until adherent growth. Hieff Trans® in vitro siRNA/miRNA Transfection Reagent Transfection Kit (Yeasen, China) was used in strict accordance with the instructions provided. The GRASLND siRNA and nontargeting negative control (NC) used in this study were obtained from Gemma Genetics with the following sequences: si-918, GGCAAGAUAAAUGACAAUAAATT, and si-1208, GUGGGAGUAGCAUCCACAUAATT.

qRT‒PCR

Total RNA was extracted from the cells using TRIzol Reagent (Yeason, China) and reverse transcribed using 5× Hifair® One Step RT SuperMix (Yeason, China). cDNA was used as a template for amplification using SYBR Green PCR mix reagent (Yeason, China). The primers used were synthesized by Guangzhou Aiki Biotechnology Co. The sequences of primers used were as follows: GRASLND, forward, 5′-AGGATTCAGGGGATGCACAG-3′ and reverse, 5′-TGGGCTGAAGATGAGACGTT-3′; GAPDH, forward, 5’-TGCACCACCAACTGCTTAGC-3’; and reverse 5’-GGCATGGACTGTGGTCATGAG-3’. GAPDH was used as a reference to normalize GRASLND expression. The relative expression of GRASLND was calculated using the 2-ΔΔCt method. Each quantitation of mRNA levels represents data from three independent experiments.

Transwell migration

After target gene silencing or overexpression, the cells were grouped and inoculated in 8 μm wells in the upper chamber of a Transwell plate (Corning, NY). The upper chamber of the Transwell system was supplemented with 200 µl of serum-free DMEM, and the lower chamber was supplemented with 600 µl of 10% FBS DMEM. After 48 h of incubation at 37 °C with 5% CO2, each well was washed twice with PBS. The migrated cells were then fixed with fixative and stained with 0.1% crystal violet. Six fields of view were randomly photographed under a microscope, and the migrated cells were counted.

Wound-healing assay

The monolayer was then scratched with a sterile 20 µl plastic pipette tip, and the detached cells were removed with PBS. The cells were cultured in low-serum DMEM or RPMI 1640 medium to inhibit cell proliferation, and images were captured using an inverted microscope (ZEISS, Axio Observer A1) at the indicated times (0, 12, 24, and 48 h).

Fluorescence in situ hybridization (FISH)

The probes used in this study were designed and synthesized by RiboBio (Guangzhou, China). Fluorescence in situ hybridization (FISH) was performed using a FISH kit (RiboBio, China) following the manufacturer’s instructions. 6 × 104 cells were seeded 24-well plate with cover glass and cultured to 60-70% density. The cells were fixed using paraformaldehyde and permeabilized with triton x-100. Subsequently, 20 µl of prehybridization solution was added to each well, and then 200 µl of hybridization solution containing the probe was added. The mixture was hybridized overnight at 37 °C in the dark. The hybridization solution was removed, and the cells were washed, keyed, sealed with an anti-fluorescence quencher, and photographed using a laser confocal microscope (Zeiss, Germany).

Dual luciferase assay

GRASLND-WT (wild-type), GRASLND-MUT (mutant), STAM2-WT (wild-type), and STAM2-MUT (mutant) melanoma cells were cultured in 6-well plates and transfected with the miR-218-5p mimic or control siRNA. 48 h after transfection, we assessed the luciferase activity of the firefly and sea pansy (Renilla) strains using the Dual-Luciferase Reporter Assay System (Promega). Relative luciferase activity was calculated from the ratio of firefly/sea pansy (Renilla) luciferase activity and normalized to that of the control. The firefly luciferase activity of each transfected well was normalized to the Renilla luciferase activity. Each assay was performed in triplicate.

Protein isolation and Western blotting

Total protein was extracted from melanoma cells using RIPA extraction reagent (Shanghai Biyuntian Company, China) and a protease inhibitor mixture (Shanghai Biyuntian Company, China). Proteins were separated by 10% SDS‒PAGE and then electroblotted onto a polyvinylidene fluoride (PVDF) membrane (Millipore). The membranes were blocked with 5% skim milk for 1 h at room temperature and incubated overnight at 4 °C with various specific primary antibodies. The next day, the blots were washed and incubated with horseradish peroxidase-conjugated secondary antibodies (Millipore) for 2 h. The protein bands were detected using the ECL Chemiluminescent Substrate Reagent (Biosharp, China). β-Actin served as a loading control.

Statistical analysis

Statistical analyses were performed using R language and GraphPad Prism 8.0.2 software. Student’s t-test was used to determine the significance of differences between the two groups. Significant differences between multiple groups were calculated using analysis of variance (ANOVA). Survival curves were generated via Kaplan‒Meier analysis and the log-rank test. Correlation analysis was performed using Pearson’s correlation coefficient. p < 0.05 was considered to indicate significance.

Results

GRASLND is highly expressed in melanoma

In this study, we first obtained the differentially expressed RNAs from the GSE15605 dataset via the GEO2R online tool. Subsequently, we obtained the differentially expressed RNAs from the TCGA dataset using the GEPIA2(http://gepia2.cancer-pku.cn/#index) online tool. Next, we analyzed the intersection of these two sets of differentially expressed RNAs and used the Ensembl database to obtain the shared differentially expressed RNAs annotated with lncRNAs. Through this series of analyses, we found 187 differentially expressed RNAs, 17 of which were lncRNAs (Fig. 1A). Based on these 17 lncRNAs, we used The least absolute shrinkage and selection operator (LASSO) regression analysis to construct a risk score model consisting of 7 lncRNAs (LINC00520, BANCR, GAPLINC, TFAP2A-AS1, LINC01234, GRASLND, and CTXND1) (Fig. 1B, C). Kaplan‒Meier curves revealed a significant difference in survival between the high-risk and low-risk patients (p < 0.05), while the 1–3–5 year receiver operating characteristic (ROC) curves of the risk model and the area under the curve (AUC) analysis indicated the strong predictive ability of the model (Fig. 1D-F). These results confirmed the high accuracy and reliability of the lncRNA-based risk score model we constructed for predicting the survival status of SKCM patients.

Fig. 1
figure 1

Construction of the marker lncRNA prognostic model. (A) Seventeen lncRNAs differentially expressed in the GSE15605 dataset and TCGA cohort. (B) Graph of the LASSO coefficient distribution of marker lncRNAs in GSE15605. (C) Graph of the tenfold cross-validation error rate. (D) Distribution of risk scores (upper panel), survival statuses (middle panel), and expression characteristics of 7 prognosis-related lncRNAs in SKCM. (E) Survival curves comparing the overall survival probability of melanoma patients in the low-risk versus high-risk groups. (F) ROC curves predicting the risk of death at 1, 3, and 5 years in melanoma patients. (G) One-way Cox regression analysis of prognostically relevant lncRNAs in 7 patients

Full size image

To further evaluate the clinical significance of these seven lncRNAs in cutaneous melanoma, we performed a one-way Cox regression analysis. The results showed that the hazard ratios (HRs) of BANCR, CTXND1, GAPLINC, LINC00520, LINC01234, TFAP2A-AS1, and GRASLND were 0.9546, 0.96241, 1.21251, 0.87747, 1.23661, 1.0247, and 1.09077, respectively (Fig. 1G). Combined with the p values, the trends of two of these molecules, GRASLND and TFAP2A-AS1, were in good agreement with the model, and we subsequently analyzed these molecules in depth.

High expression of GRASLND in melanoma is significantly associated with poor prognosis

To further validate the two molecules screened in the prescreening, we first assessed the expression level of GRASLND in melanoma and normal tissues by qRT–PCR, and the results showed that its expression was significantly upregulated in SKCM (Fig. 2A). In addition, we analyzed the relationship between GRASLND in the TCGA SKCM cohort and the clinical characteristics of the patients. GRASLND expression was significantly greater in the group of patients older than 65 years than in the group of patients younger than 65 years (Fig. 2B). However, the expression of GRASLND has no difference for any of the clinical characteristics, including sex, Clark stage, clinical stage, T stage, M stage, and N stage (Fig. 2C-H).

Fig. 2
figure 2

Prognostic value of GRASLND in melanoma patients. (A) Differential expression of GRASLND in normal and SKCM tumor tissues measured by qRT–PCR. (B-H) Differential expression of GRASLND between different age groups, sexes, Clark staging periods, clinical staging periods, T staging periods, M staging periods, and N staging periods. (I, J) TCGA data were analyzed to evaluate the correlation between the expression of GRASLND and OS/DSS. (K, L) Correlations of GRASLND expression with OS/RFS in different datasets

Full size image

To further explore the relationship between GRASLND expression levels and survival outcomes in patients with SKCM, we first analyzed the relationship between GRASLND expression and the corresponding prognostic information using SKCM cohort data from the TCGA database. The overall survival (OS) analysis revealed that patients with high GRASLND expression had lower survival rates (Fig. 2I). However, disease-specific survival (DSS) analysis showed that survival was not significant in patients with high GRASLND expression levels compared to those with low expression (Fig. 2J). These results suggest that high GRASLND expression levels are strongly associated with poor prognosis in melanoma patients. Finally, we used two external datasets, GSE190113 and GSE33713, to analyze the correlation between the expression level of GRASLND and OS and relapse-free survival (RFS) outcomes. Consistent with the aforementioned TCGA SKCM cohort data, in the GSE190113 dataset, the survival analyses revealed that patients with high GRASLND expression had poor OS (Fig. 2K). The GSE33713 dataset analysis revealed that patients with high GRASLND expression also exhibited poor RFS (Fig. 2L). These results further support the association between GRASLND expression and poor prognosis.

The expression of GRASLND negatively correlates with immune infiltration and immune molecules expression

To further investigate the relationship between GRASLND expression and immune cell infiltration, we calculated the correlation between GRASLND expression and the abundance of various immune cells (e.g., B cells, CD4 + T cells, CD8 + T cells, endothelial cells, macrophages, and NK cells). The results showed a negative correlation between GRASLND expression and CD8 + T-cell abundance. This finding suggested that GRASLND expression may also affect the infiltration of this immune cell subset (Fig. 3A). We further investigated the relationship between GRASLND expression and the immune checkpoint molecules expression(e.g., LAG3, PDCD1, and TIGIT). We compared the expression levels of immune checkpoint molecules in the high GRASLND expression (G1) and low GRASLND expression(G2) groups and used independent two-sample t-tests to assess the differences between the two groups. We found that LAG3 expression was significantly greater in the G1 subgroup than in the G2 subgroup, whereas the expression of PDCD1 and TIGIT was significantly lower in the G2 subgroup than in the G1 subgroup. These results suggest that GRASLND expression may impact the expression of immune checkpoint molecules, thereby modulating the function of immune cells in the tumor microenvironment (Fig. 3B). Furthermore, we found that patients with high GRASLND expression had lower OS and PFS rates (Fig. 3C-D). In conclusion, our analysis revealed a negative correlation between GRASLND expression immune cell infiltration, and immune checkpoint molecule expression.

Fig. 3
figure 3

Correlations of GRASLND expression with immune infiltration and immunotherapy outcomes. (A) Correlation analysis of GRASLND expression with the abundances of different immune cells. (B) Comparison of the expression levels of immune checkpoint molecules in the GRASLND high expression (G1) and low expression (G2) groups. (C, D) Correlation between GRASLND expression and OS/PFS outcomes

Full size image

GRASLND significantly promotes SKCM cell migration

To determine the cell types expressing GRASLND, we extracted the expression levels of GRASLND in cell subpopulation species from two datasets, GSE72056 and GSE115978_aPD1. The results showed that GRASLND was expressed in eight cell types in both datasets: B cells, CD4 + T cells, CD8 + T cells, endothelial cells, fibroblasts, malignant cells, monocytes/macrophages, and value-added T cells (Fig. 4A). The localization of GRASLND in A2058 and SK-MEL-28 was examined by FISH. The results showed that GRASLND was located mainly in the nucleus and partially in the cytoplasm (Fig. 4B). We used two siRNA(si-918 and si-1208) to knock down GRASLND in SKCM cells (Fig. 4C). Silencing of GRASLND decreased the wound healing rate, indicating that silencing of GRASLND inhibited cell migratory capacity (Fig. 4D). In the Transwell assays, GRASLND silencing significantly reduced the number of migrating cells (Fig. 4E). We also constructed a GRASLND overexpression model by lentiviral transduction in A2058 cells (Fig. 4F). In agreement with the findings of the previous experiments, GRASLND overexpression significantly enhanced the migratory ability of the cells (Fig. 4G, H). These findings suggest that GRASLND plays an important role in promoting the migration of SKCM cells.

Fig. 4
figure 4

Construction of the marker lncRNA prognostic model. (A) Expression levels of GRASLND in various cells in the GSE72056 and GSE115978_aPD1 single-cell transcriptome datasets. (B) Fluorescence in situ hybridization showed that GRASLND was expressed both in the nucleus and cytoplasm of A2058 and SK-MEL-28 cells. (C) Transfection of siRNA or NC targeting GRASLND and measurement of GRASLND expression levels in A2058 and SK-MEL-28 cells by qRT‒PCR. (D) The wound healing assay showed that the knockdown of GRASLND inhibited the migration ability of A2058 cells and SK-MEL-28 cells. (E) Migration assays showed that GRASLND knockdown significantly inhibited the migration of A2058 and SK-MEL-28 cells. (F) After overexpressing GRASLND, the expression level of GRASLND in A2058 cells was measured by qRT‒PCR. (G) The wound healing assay showed that the migration ability of A2058 cells was enhanced after overexpression of GRASLND. (H) Migration assays showed that overexpression of GRASLND significantly enhanced the migration ability of A2058 cells

Full size image

GRASLND functions as a sponge of miR-218-5p

After identifying GRASLND as an oncogene in SKCM, our next goal was to determine the underlying mechanisms involved. Considering the cytoplasmic localization of GRASLND in SKCM cells, we used online tools such as TargetScan 7.2 (http://www.targetscan.org/vert_72/) and the GSE database to predict the potential binding sites of miRNAs to GRASLND (Fig. 5A-D). We detected elevated miR-218-5p expression after transfection of the NC or miR-218-5p mimics (Fig. 5F). Subsequently, we specifically explored the relationship between miR-218-5p and GRASLND.

Fig. 5
figure 5

Flow chart of the screening of miRNAs and target genes associated with GRASLND. (A, C) Volcano plots of differentially expressed genes. (B, D) Venn diagram of differential gene expression. (E) Flowchart showing target gene annotation. (F) Transfection of mimics or NC targeting miR-218-5p and measurement of miR-218-5p expression levels in cells by qRT‒PCR. (G, H) A dual-luciferase reporter gene assay was used to measure luciferase activity in SK-MEL-28 cells

Full size image

We mutated the binding sequence of miR-218-5p to GRASLND and cotransfected SK-MEL-28 cells with reporter vectors carrying the wild-type GRASLND sequence, effectively reducing luciferase activity. However, mutation of the miR-218-5p binding site eliminated the inhibitory effect of the miR-218-5p mimic on the GRASLND- and STAM2-driven luciferase activities (Fig. 5G, H).

GRASLND positively regulates STAM2 expression by sponging miR-218-5p

To further explore whether GRASLND acts as a miRNA sponge to positively regulate mRNA expression in a competing endogenous RNA (ceRNA)-dependent manner, we selected SKCM cells for in-vitro experiments. First, we determined that the expression level of GRASLND was significantly reduced after transfection with shGRASLND (Fig. 6A), the expression level of miR-218-5p was also significantly decreased after transfection with the miR-218-5p inhibitor (Fig. 6B), and the expression level of STAM2 was significantly increased after transfection with STAM2 OE (Fig. 6C). In SKCM cells, miR-218-5p attenuated the effect of GRASLND. shGRASLND in A2058 and SK-MEL-28 cells resulted in significantly greater miR-218-5p expression and significantly lower STAM2 mRNA expression (Fig. 6D-F). A2058 and SK-MEL-28 cells were transfected with shNC, shGRASLND, or shGRASLND in combination with the miR-218-5p inhibitor or with the STAM2 OE to investigate the effect of miR-218-5p inhibition on STAM2 expression. shGRASLND reduced STAM2 expression, while the miR-218-5p inhibitor partially inhibited this effect (Fig. 6D-F). Given that GRASLND can inhibit miR-218-5p expression, we investigated whether miR-218-5p affects the GRASLND knockdown-mediated inhibition of SKCM cell progression. As expected, we observed that the downregulation of miR-218-5p expression abrogated the shGRASLND-induced reduction in A2058 cell migration and invasion (Fig. 6G-J). Compared with those of cells transfected with shGRASLND alone, the migration and invasion of A2058 cells were significantly greater after cotransfection with the miR-218-5p inhibitor and shGRASLND (Fig. 6G-J). Similarly, the migratory and invasive abilities of A2058 cells were significantly greater after cotransfection of STAM2 OE and shGRASLND than after transfection with shGRASLND alone (Fig. 6G-J). These results suggest that the downregulation of GRASLND expression inhibits SKCM invasion and migration by negatively regulating miR-218-5p expression and that GRASLND can act as a miR-218-5p sponge to positively regulate STAM2 expression.

Fig. 6
figure 6

GRASLND affects the migration and invasion of SKCM cells by regulating miR-218-5p expression. (A) After silencing GRASLND, qRT‒PCR was performed to measure GRASLND expression. (B) After transfection of cells with the miR-218-5p inhibitor or NC, qRT‒PCR was performed to measure miR-218-5p expression. (C) After overexpressing GRASLND, qRT‒PCR was performed to measure GRASLND expression. (D, E) The RNA and protein expression levels of GRASLND, miR-218-5p, and STAM2 in A2058 cells and SK-MEL-28 cells were measured by qRT‒PCR and Western blotting after the transfection of NC or shGRASLND and after the cotransfection of shGRASLND with the miR-218-5p inhibitor. (F) Grayscale statistics of the protein expression of STAM2 in A2058 cells and SK-MEL-28 cells after transfection with NC or shGRASLND and cotransfection of shGRASLND with the miR-218-5p inhibitor. (G-J) Effect on the migration and invasion of A2058 cells after transfection with NC or shGRASLND or cotransfection of shGRASLND with the miR-218-5p inhibitor

Full size image

Discussion

Skim cutaneous melanoma is one of the most life-threatening skin cancers worldwide and one of the most common malignant tumors [26, 27]. However, the molecular mechanisms underlying its malignant behavior are poorly understood. LncRNAs are involved in various cellular processes and have been implicated in the development and progression of malignant tumors [28,29,30,31]. In the present study, we screened for a key lncRNA, GRASLND, which is associated with SKCM immunotherapy, promotes SKCM metastasis, and is significantly correlated with SKCM malignant behavior, disease-specific survival, OS, and PFS rate. In addition, GRASLND expression is associated with immune infiltration and immunotherapy outcomes in SCKM patients. Our findings will potentially lead to the discovery of new candidate markers for the immunization against and treatment of SKCM.

Little is known about GRASLND, which is located on human chromosome 2p25.1 and plays a key role in the malignant progression of tumors [32,33,34]. Increased expression of GRASLND in gastric cancer tissues has been reported to promote tumor metastasis and proliferation by targeting the miR-30c-2-3p/LOX axis, suggesting a strong correlation between GRASLND, poor prognosis, and advanced disease in SKCM patients [33]. In addition, GRASLND promotes mesenchymal stem cell chondrogenesis by inhibiting the interferon type II signaling pathway [35]. Furthermore, GRASLND promotes bladder cancer progression through the miR-455-5p/SOX11 axis [36]. In our study, we revealed increased expression of GRASLND in SKCM tissues, suggesting a strong correlation between GRASLND, poor prognosis, and immunotherapy in SKCM patients. In addition, functional experiments confirmed the promotional effect of GRASLND on SKCM migration and invasion. We also found mutual inhibition of GRASLND and miR-218-5p expression, suggesting that GRASLND has a miR-218-5p-mediated tumor suppressor effect on SKCM cells. We also demonstrated that GRASLND acts as a miR-218-5p sponge to positively regulate STAM2 expression in SKCM cells. Thus, our results demonstrate the oncogenic role of GRASLND in the malignant process of SKCM.

The ceRNA mode is a characteristic way for lncRNAs to participate in gene regulation, whereby lncRNAs competitively bind to miRNAs to eliminate the endogenous repressive effects of miRNAs on their target transcripts [37, 38]. In the present study, GRASLND was suggested to act as a ceRNA by binding to miR-218-5p. It has been shown that miR-218-5p is associated with the regulation of proliferation and migration through the function of specific target genes [39,40,41]. For example, miR-218-5p in the endometrial microenvironment prevents the migration of ectopic endometrial stromal cells by inhibiting LASP1 [42]. miR-218-5p regulates skin and hair follicle development by targeting SFRP2 through the Wnt/β-catenin signaling pathway [43]. miR-218-5p targets LHFPL3 to regulate human glioma cell proliferation, migration, and epithelial–mesenchymal transition [44]. Similarly, we found that downregulation of GRASLND significantly increased miR-218-5p levels. Subsequently, by predicting potential miR-218-5p targets via TargetScan, miRTarBase, and miRDB, we screened STAM2 as a target of miR-218-5p. STAM2 is a phosphotyrosine protein that is a member of the endosomal sorting complex required for transport (ESCRT-0) [45, 46]. Kaymaz et al. revealed that STAM2 is a functional target of miR-2278, which acts as a tumor suppressor by inhibiting leukemia cell proliferation and inducing apoptosis [47]. STAM2 knockdown has been reported to inhibit proliferation, migration, and invasion by affecting the JAK2/STAT3 signaling pathway in gastric cancer [46]. Our rescue experiments further confirmed that overexpression of STAM2 reversed the effects of miR-218-5p on cells in vitro. Moreover, most lncRNAs function by interacting with various proteins [48]. Similarly, in SKCM, considering its localization, GRASLND may interact with STAM2 proteins and alter gene expression or participate in gene regulation as a cofactor of transcription factors and chromatin modifiers [49, 50]. These data suggest that miR-218-5p regulates STAM2 expression in SKCM, thus regulating SKCM cell function.

In conclusion, our study highlights the value of GRASLND as a prognostic marker or therapeutic target for SKCM. GRASLND predicts poor prognosis in SKCM patients and is correlated with tumor metastasis. In addition, GRASLND promoted SKCM invasion and migration by regulating STAM2 expression through the competitive uptake of miR-2218-5p. One of the shortcomings of this study is that no animal experiments were performed, and in vivo experiments are essential; therefore, we will perform these experiments in the future.

Data Availability

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

References

  1. Schadendorf D, van Akkooi ACJ, Berking C, Griewank KG, Gutzmer R, Hauschild A, et al. Melanoma Lancet. 2018;392(10151):971–84.

    Article  PubMed  Google Scholar 

  2. Carr S, Smith C, Wernberg J. Epidemiology and risk factors of Melanoma. Surg Clin North Am. 2020;100(1).

  3. Rebecca VW, Somasundaram R, Herlyn M. Pre-clinical modeling of cutaneous melanoma. Nat Commun. 2020;11(1):2858.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Rafiee A, Riazi-Rad F, Havaskary M, Nuri F. Long noncoding RNAs: regulation, function and cancer. Biotechnol Genet Eng Rev. 2018;34(2):153–80.

    Article  CAS  PubMed  Google Scholar 

  5. Luo W, Wang J, Xu W, Ma C, Wan F, Huang Y, et al. LncRNA RP11-89 facilitates tumorigenesis and ferroptosis resistance through PROM2-activated iron export by sponging mir-129-5p in bladder cancer. Cell Death Dis. 2021;12(11):1043.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Li Y, Shi B, Dong F, Zhu X, Liu B, Liu Y. LncRNA KCNQ1OT1 facilitates the progression of bladder cancer by targeting MiR-218-5p/HS3ST3B1. Cancer Gene Ther. 2021;28(3–4):212–20.

    Article  CAS  PubMed  Google Scholar 

  7. Lou Y, Yu Y, Xu X, Zhou S, Shen H, Fan T, et al. Long non-coding RNA LUCAT1 promotes tumourigenesis by inhibiting ANXA2 phosphorylation in hepatocellular carcinoma. J Cell Mol Med. 2019;23(3):1873–84.

    Article  CAS  PubMed  Google Scholar 

  8. Yuan K, Lan J, Xu L, Feng X, Liao H, Xie K, et al. Long noncoding RNA TLNC1 promotes the growth and metastasis of liver cancer via inhibition of p53 signaling. Mol Cancer. 2022;21(1):105.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Fang D, Ou X, Sun K, Zhou X, Li Y, Shi P, et al. m6A modification-mediated lncRNA TP53TG1 inhibits gastric cancer progression by regulating CIP2A stability. Cancer Sci. 2022;113(12):4135–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Zhu Y, Huang C, Zhang C, Zhou Y, Zhao E, Zhang Y, et al. LncRNA MIR200CHG inhibits EMT in gastric cancer by stabilizing miR-200c from target-directed miRNA degradation. Nat Commun. 2023;14(1):8141.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Ouyang S, Zhou X, Chen Z, Wang M, Zheng X, Xie M. LncRNA BCAR4, targeting to miR-665/STAT3 signaling, maintains cancer stem cells stemness and promotes tumorigenicity in colorectal cancer. Cancer Cell Int. 2019;19:72.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Peng CL, Zhao XJ, Wei CC, Wu JW. LncRNA HOTAIR promotes colon cancer development by down-regulating miRNA-34a. Eur Rev Med Pharmacol Sci. 2019;23(13):5752–61.

    PubMed  Google Scholar 

  13. Feng J, Zhang J, Li Y, Cheng W, Liu Y, Chen Z, et al. Inhibition of lncRNA PCAT19 promotes breast cancer proliferation. Cancer Med. 2023;12(10):11971–82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Zeng H, Hou Y, Zhou X, Lang L, Luo H, Sun Y, et al. Cancer-associated fibroblasts facilitate premetastatic niche formation through lncRNA SNHG5-mediated angiogenesis and vascular permeability in breast cancer. Theranostics. 2022;12(17):7351–70.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Yang Q, Wang M, Xu J, Yu D, Li Y, Chen Y, et al. LINC02159 promotes non-small cell lung cancer progression via ALYREF/YAP1 signaling. Mol Cancer. 2023;22(1):122.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Chen Y-L, Liu Y-N, Lin Y-T, Tsai M-F, Wu S-G, Chang T-H, et al. LncRNA SLCO4A1-AS1 suppresses lung cancer progression by sequestering the TOX4-NTSR1 signaling axis. J Biomed Sci. 2023;30(1):80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Wen S, Wei Y, Zen C, Xiong W, Niu Y, Zhao Y. Long non-coding RNA NEAT1 promotes bone metastasis of prostate cancer through N6-methyladenosine. Mol Cancer. 2020;19(1):171.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Huang W, Chen Q, Lu Y, Kong Z, Wan X, Huang Y et al. Androgen-responsive oncogenic lncRNA RP11-1023L17.1 enhances c-Myc protein Stability in prostate Cancer. Int J Mol Sci. 2022;23(20).

  19. Mao T-L, Fan M-H, Dlamini N, Liu C-L. LncRNA MALAT1 facilitates ovarian Cancer progression through promoting Chemoresistance and Invasiveness in the Tumor Microenvironment. Int J Mol Sci. 2021;22:19.

    Article  Google Scholar 

  20. Wang C, Wang J, Shen X, Li M, Yue Y, Cheng X, et al. LncRNA SPOCD1-AS from ovarian cancer extracellular vesicles remodels mesothelial cells to promote peritoneal metastasis via interacting with G3BP1. J Exp Clin Cancer Res. 2021;40(1):101.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Luo A, Lan X, Qiu Q, Zhou Q, Li J, Wu M, et al. LncRNA SFTA1P promotes cervical cancer progression by interaction with PTBP1 to facilitate TPM4 mRNA degradation. Cell Death Dis. 2022;13(11):936.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Fan Y, Sheng W, Meng Y, Cao Y, Li R. LncRNA PTENP1 inhibits cervical cancer progression by suppressing miR-106b. Artif Cells Nanomed Biotechnol. 2020;48(1):393–407.

    Article  CAS  PubMed  Google Scholar 

  23. Leucci E, Vendramin R, Spinazzi M, Laurette P, Fiers M, Wouters J, et al. Melanoma addiction to the long non-coding RNA SAMMSON. Nature. 2016;531(7595):518–22.

    Article  CAS  PubMed  Google Scholar 

  24. Schmidt K, Joyce CE, Buquicchio F, Brown A, Ritz J, Distel RJ, et al. The lncRNA SLNCR1 mediates Melanoma Invasion through a conserved SRA1-like region. Cell Rep. 2016;15(9):2025–37.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Iyer MK, Niknafs YS, Malik R, Singhal U, Sahu A, Hosono Y, et al. The landscape of long noncoding RNAs in the human transcriptome. Nat Genet. 2015;47(3):199–208.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Pietrobono S, Anichini G, Sala C, Manetti F, Almada LL, Pepe S, et al. ST3GAL1 is a target of the SOX2-GLI1 transcriptional complex and promotes melanoma metastasis through AXL. Nat Commun. 2020;11(1):5865.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Dimitriou F, Long GV, Menzies AM. Novel adjuvant options for cutaneous melanoma. Ann Oncol. 2021;32(7):854–65.

    Article  CAS  PubMed  Google Scholar 

  28. Gibb EA, Brown CJ, Lam WL. The functional role of long non-coding RNA in human carcinomas. Mol Cancer. 2011;10:38.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Sanchez Calle A, Kawamura Y, Yamamoto Y, Takeshita F, Ochiya T. Emerging roles of long non-coding RNA in cancer. Cancer Sci. 2018;109(7):2093–100.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Yu J, Liu Y, Guo C, Zhang S, Gong Z, Tang Y, et al. Upregulated long non-coding RNA LINC00152 expression is associated with progression and poor prognosis of tongue squamous cell carcinoma. J Cancer. 2017;8(4):523–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Zhu Y, He D, Bo H, Liu Z, Xiao M, Xiang L, et al. The MRVI1-AS1/ATF3 signaling loop sensitizes nasopharyngeal cancer cells to paclitaxel by regulating the Hippo-TAZ pathway. Oncogene. 2019;38(32):6065–81.

    Article  CAS  PubMed  Google Scholar 

  32. Tong X, Yang Z, Wang Q, Zhang D. RNF144A-AS1 promotes the development of glioma cells by targeting miR-665/HMGA1 axis. Neurosci Lett. 2021;765:136259.

    Article  CAS  PubMed  Google Scholar 

  33. Li Z, Shi L, Li X, Wang X, Wang H, Liu Y. RNF144A-AS1, a TGF-β1- and hypoxia-inducible gene that promotes tumor metastasis and proliferation via targeting the miR-30c-2-3p/LOX axis in gastric cancer. Cell Biosci. 2021;11(1):177.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. LncRNA RNF144A-AS1. Promotes bladder Cancer Progression via RNF144A-AS1/miR-455-5p/SOX11 Axis [Retraction]. Onco Targets Ther. 2022;15:681–2.

    Article  Google Scholar 

  35. non-coding RNA GRASLND enhances chondrogenesis via suppression of the interferon type II signaling pathway.pdf>.

  36. Bi H, Shang Z, Jia C, Wu J, Cui B, Wang Q, et al. LncRNA RNF144A-AS1 promotes bladder Cancer Progression via RNF144A-AS1/miR-455-5p/SOX11 Axis. Onco Targets Ther. 2020;13:11277–88.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Nie W, Ge H-j, Yang X-q, Sun X, Huang H, Tao X et al. LncRNA-UCA1 exerts oncogenic functions in non-small cell lung cancer by targeting miR-193a-3p. Cancer Lett. 2016;371(1).

  38. Liu B, Li J, Cairns MJ. Identifying miRNAs, targets and functions. Brief Bioinform. 2014;15(1).

  39. Yamamoto N, Kinoshita T, Nohata N, Itesako T, Yoshino H, Enokida H, et al. Tumor suppressive microRNA-218 inhibits cancer cell migration and invasion by targeting focal adhesion pathways in cervical squamous cell carcinoma. Int J Oncol. 2013;42(5):1523–32.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Liu Z, Mao L, Wang L, Zhang H, Hu X. miR–218 functions as a tumor suppressor gene in cervical cancer. Mol Med Rep. 2020;21(1):209–19.

    CAS  PubMed  Google Scholar 

  41. Kogo R, How C, Chaudary N, Bruce J, Shi W, Hill RP, et al. The microRNA-218 ~ survivin axis regulates migration, invasion, and lymph node metastasis in cervical cancer. Oncotarget. 2015;6(2):1090–100.

    Article  PubMed  Google Scholar 

  42. Zhang Z, Wang Y, Zeng L, Yu K, Wang Y, Luo Y, et al. Mir-218-5p in endometrial microenvironment prevents the migration of ectopic endometrial stromal cells by inhibiting LASP1. Reprod Biol Endocrinol. 2022;20(1):64.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Zhao B, Chen Y, Yang N, Chen Q, Bao Z, Liu M, et al. Mir-218-5p regulates skin and hair follicle development through Wnt/β-catenin signaling pathway by targeting SFRP2. J Cell Physiol. 2019;234(11):20329–41.

    Article  CAS  PubMed  Google Scholar 

  44. Li Z, Qian R, Zhang J, Shi X. MiR-218-5p targets LHFPL3 to regulate proliferation, migration, and epithelial-mesenchymal transitions of human glioma cells. Biosci Rep. 2019;39(3).

  45. Mizuno E, Kawahata K, Okamoto A, Kitamura N, Komada M. Association with hrs is required for the early endosomal localization, stability, and function of STAM. J Biochem. 2004;135(3):385–96.

    Article  CAS  PubMed  Google Scholar 

  46. Yang Y, Zhang Q, Liang J, Yang M, Wang Z, Tang D, et al. STAM2 knockdown inhibits proliferation, migration, and invasion by affecting the JAK2/STAT3 signaling pathway in gastric cancer. Acta Biochim Biophys Sin (Shanghai). 2021;53(6):697–706.

    Article  CAS  PubMed  Google Scholar 

  47. Kaymaz BT, Günel NS, Ceyhan M, Çetintaş VB, Özel B, Yandım MK, et al. Revealing genome-wide mRNA and microRNA expression patterns in leukemic cells highlighted hsa-miR-2278 as a tumor suppressor for regain of chemotherapeutic imatinib response due to targeting STAT5A. Tumour Biol. 2015;36(10):7915–27.

    Article  CAS  PubMed  Google Scholar 

  48. Ferrè F, Colantoni A, Helmer-Citterich M. Revealing protein-lncRNA interaction. Brief Bioinform. 2016;17(1):106–16.

    Article  PubMed  Google Scholar 

  49. Zhu P, He F, Hou Y, Tu G, Li Q, Jin T, et al. A novel hypoxic long noncoding RNA KB-1980E6.3 maintains breast cancer stem cell stemness via interacting with IGF2BP1 to facilitate c-Myc mRNA stability. Oncogene. 2021;40(9):1609–27.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Sun Q, Hao Q, Prasanth KV. Nuclear Long Noncoding RNAs: key regulators of Gene expression. Trends Genet. 2018;34(2):142–57.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This study was funded by the National Natural Science Foundation of China (82071101), and Guangdong Basic and Applied Basic Research Foundation (2021A1515011142, 2023A1515012343, 2022A1515220099), Supported by 2020 Li Ka Shing Foundation Cross-Disciplinary Research Grant (2020LKSFG18B, 2020LKSFG02E), and Guangdong University Innovation Team Project (2021KCXTD047), Provincial science and technology innovation strategy special project funding program (210714106901245, STKJ202209067, STKJ2023004).

Author information

Author notes

  1. Aiwei Ma, Wenqi Shi, Liyun Chen and Zijian Huang contributed equally to this work.

Authors and Affiliations

  1. Department of Plastic Surgery and Burns Center, Second Affiliated Hospital, Shantou University Medical College, Shantou, Guangdong, 515051, China

    Aiwei Ma, Wenqi Shi, Zijian Huang, Yiwen Zhang, Zixuan Tang, Wenshi Jiang, Mengjing Xu, Wancong Zhang & Shijie Tang

  2. Plastic Surgery Institute of Shantou University Medical College, Shantou, Guangdong, 515051, China

    Aiwei Ma, Wenqi Shi, Liyun Chen, Zijian Huang, Yiwen Zhang, Zixuan Tang, Wenshi Jiang, Mengjing Xu, Wancong Zhang & Shijie Tang

  3. Shantou Plastic Surgery Clinical Research Center, Shantou, Guangdong, 515051, China

    Aiwei Ma, Wenqi Shi, Liyun Chen, Zijian Huang, Yiwen Zhang, Zixuan Tang, Wenshi Jiang, Mengjing Xu, Wancong Zhang & Shijie Tang

  4. Research Center of Translational Medicine, Second Affiliated Hospital of Shantou University Medical College, Shantou, Guangdong, 515051, China

    Liyun Chen

  5. Department of Plastic and Reconstructive Surgery, Central South University Third Xiangya Hospital, Changsha, China

    Jianda Zhou

Authors
  1. Aiwei Ma
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Wenqi Shi
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Liyun Chen
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Zijian Huang
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Yiwen Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Zixuan Tang
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Wenshi Jiang
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Mengjing Xu
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Jianda Zhou
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Wancong Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

  11. Shijie Tang
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

Aiwei Ma contributed to data analysis and drafted the manuscript. Wancong Zhang, Jianda Zhou and Shijie Tang contributed to data verification. Zijian Huang contributed to data acquisition and data analysis. Aiwei Ma, Wenqi Shi, Yiwen Zhang, Zixuan Tang, Wenshi Jiang and Mengjing Xu conducted the experiment. Wenqi Shi and Liyun Chen contributed to the study design and revision of the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Wancong Zhang or Shijie Tang.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

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

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 http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) 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

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Ma, A., Shi, W., Chen, L. et al. GRASLND regulates melanoma cell progression by targeting the miR-218-5p/STAM2 axis. J Transl Med 22, 684 (2024). https://doi.org/10.1186/s12967-024-05397-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12967-024-05397-z

Keywords

Journal of Translational Medicine

ISSN: 1479-5876

Contact us