- Research article
- Open Access
Reconstruction of lncRNA-miRNA-mRNA network based on competitive endogenous RNA reveals functional lncRNAs in skin cutaneous melanoma
BMC Cancer volume 20, Article number: 927 (2020)
Human skin cutaneous melanoma is the most common and dangerous skin tumour, but its pathogenesis is still unclear. Although some progress has been made in genetic research, no molecular indicators related to the treatment and prognosis of melanoma have been found. In various diseases, dysregulation of lncRNA is common, but its role has not been fully elucidated. In recent years, the birth of the “competitive endogenous RNA” theory has promoted our understanding of lncRNAs.
To identify the key lncRNAs in melanoma, we reconstructed a global triple network based on the “competitive endogenous RNA” theory. Gene Ontology and KEGG pathway analysis were performed using DAVID (Database for Annotation, Visualization, and Integration Discovery). Our findings were validated through qRT-PCR assays. Moreover, to determine whether the identified hub gene signature is capable of predicting the survival of cutaneous melanoma patients, a multivariate Cox regression model was performed.
According to the “competitive endogenous RNA” theory, 898 differentially expressed mRNAs, 53 differentially expressed lncRNAs and 16 differentially expressed miRNAs were selected to reconstruct the competitive endogenous RNA network. MALAT1, LINC00943, and LINC00261 were selected as hub genes and are responsible for the tumorigenesis and prognosis of cutaneous melanoma.
MALAT1, LINC00943, and LINC00261 may be closely related to tumorigenesis in cutaneous melanoma. In addition, MALAT1 and LINC00943 may be independent risk factors for the prognosis of patients with this condition and might become predictive molecules for the long-term treatment of melanoma and potential therapeutic targets.
Human skin cutaneous melanoma (SKCM) is the most common and dangerous type of skin tumour [1, 2]. Worldwide, approximately 232,000 (1.7%) cases of cutaneous melanoma are reported among all newly diagnosed primary malignant cancers, and this disease results in approximately 55,500 cancer deaths (0.7% of all cancer deaths) [1, 3]. The incidence of melanoma in Australia, New Zealand, Norway, Sweden, the UK, and the USA from 1982 to 2011 has shown increases of approximately 3% annually and will further increase until 2022 . In 2015, there were 3.1 million people with melanoma, resulting in 59,800 deaths . Nevertheless, 95,710 cases of melanoma in situ will be newly diagnosed in 2020 . The high incidence and high mortality of melanoma indicate that researchers must further study this disease. Although some achievements have been made in the genetic research of melanoma, markers related to diagnosis and treatment are needed.
Tumorigenesis often results from aberrant transcriptomes, including aberrant levels of coding RNA and noncoding RNA [6,7,8]. It has been proven that lncRNAs have various effects, including regulation of gene transcription, post-transcriptional regulation and epigenetic regulation [9,10,11,12]. In addition, dysregulation of lncRNAs has been observed in various diseases [13,14,15,16]. Unfortunately, the functions of lncRNAs are more difficult to identify than those of coding RNAs. Until now, only a few lncRNAs have been identified as crucial factors in the tumorigenesis and development of melanoma, including ZNNT1, THOR and SAMMSON [14, 15, 17]. Thus, how to locate them and define their functions is a challenge of current research.
The effect of miRNAs on malignancies has been verified in many ways. Studies have suggested that lncRNAs can regulate miRNA abundance by binding and sequestering them . Thus, we aimed to study the function of lncRNAs by studying the interactions among lncRNAs, mRNAs and miRNAs. In 2011, the competitive endogenous RNA (ceRNA) hypothesis proposed a novel regulatory mechanism between noncoding RNA and coding RNA [19,20,21]. This theory indicated that any RNA transcript harbouring miRNA-response elements (MREs) can sequester miRNAs from other targets sharing the same MREs and thereby regulate their expression [19,20,21]. That is, the RNA transcripts that can be cross regulated by each other can be biologically predicted according to their common MREs [20, 22]. Evidence has shown that ceRNAs exist in several species and contexts and might play an important role in various biological processes, such as tumorigenesis . Systematic analysis of the ceRNA network has been performed in multiple tumours, such as gastric cancer, bladder cancer, and ovarian cancer, contributing to a better understanding of tumorigenesis and facilitating the development of lncRNA-directed diagnostics and therapeutics against this disease [23,24,25]. Unfortunately, however, such functional interactions have not yet been elucidated in melanoma.
In this study, we used bioinformatics methods to construct the ceRNA network of cutaneous melanoma and to identify the key lncRNAs involved in melanomagenesis. Through the reconstruction of a ceRNA network, we identified and verified that the key ceRNA molecules play a crucial role in the tumorigenesis and prognosis of SKCM. (Work flow was shown in Fig. 1).
Human melanoma miRNA expression data were downloaded from the NCBI GEO database (GEO (http://www.ncbi.nlm.nih.gov/geo) , including GSE24996, GSE35579, and GSE62372, which are array-based datasets. The GSE24996 dataset consists of 8 benign nevus tissue samples and 23 primary melanoma tissue samples. The GSE35579 dataset consists of 11 benign nevus tissue samples and 20 primary melanoma tissue samples. The GSE62372 dataset consists of 9 benign nevus tissue samples and 92 primary melanoma tissue samples. mRNA and lncRNA expression data were also downloaded from the NCBI GEO database (GSE112509), which is a sequence-based dataset. The GSE112509 dataset consists of 23 benign nevus tissue samples and 57 primary melanoma tissue samples.
Identification of DEMis, DELs and DEMs
For identification of the differentially expressed miRNAs (DEMis) between primary melanoma and benign nevus samples, “R” (version 3.4.2, https://www.r-project.org/)  was used with the “limma” package after normalization . For identification of the differentially expressed lncRNAs (DELs) and mRNAs (DEMs) between primary melanoma and benign nevus samples, “R” (version 3.4.2, https://www.r-project.org/)  was used with the “DESeq2” package . The DEMis, DELs and DEMs were selected according to |log2FC| > 1 and adjusted P-value < 0.05.
Prediction of target lncRNAs and mRNAs
For prediction of the target lncRNAs and mRNAs through DEMis, starBase (starbase.sysu.edu.cn) was used in our study . Multiple lncRNA/mRNA-predicting programmes (PITA, RNA22, miRmap, DIANA-microT, miRanda, PicTar and TargetScan) were used in starBase . For accuracy, only when the target mRNA was predicted in at least four predicted programmes on starBase would it be chosen as the predicted target mRNA. Then, these predicted target lncRNAs and mRNAs were merged with DEMs and DELs, respectively.
Reconstruction of the ceRNA network
The ceRNA network was reconstructed based on ceRNA theory  and as follows: (1) Expression correlation between DELs and DEMs was evaluated using the Pearson correlation coefficient (PCC). The DEL-DEM pairs with PCC > 0.4 and P-value < 0.01 were considered coexpressed lncRNA-mRNA pairs. (2) Both lncRNAs and mRNAs in the pairs were negatively correlated with their common miRNAs. (3) The ceRNA network was reconstructed and visualized using Cytoscape (version 3.7.1, https://cytoscape.org/) [31, 32].
Functional enrichment analysis
For functional enrichment, Gene Ontology (GO) biological process (BP), cell component (CC), molecular function (MF) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of mRNAs in the ceRNA network were performed using DAVID (version 6.8, https://david.ncifcrf.gov/) [33, 34].
Hub gene selection and reconstruction of key ceRNA subnetworks
To reconstruct our key ceRNA subnetwork, we first selected hub genes according to the node degrees of the ceRNA network we reconstructed above by calculating the number of lncRNA-miRNA and miRNA-mRNA pairs. For these key lncRNAs, GO-BP, GO-CC, GO-MF and KEGG pathway annotation were performed according to their first mRNA neighbours by using DAVID (version 6.8, https://david.ncifcrf.gov/) [33, 34].
Sample selection for qRT-PCR validation
To validate findings in the ceRNA network, we selected the top three hub genes to determine their expression in cutaneous melanoma and skin tissues. Twelve patients with cutaneous melanoma and three healthy patients were included in this study. The study protocol was approved by the Ethics Committee of The First Affiliated Hospital, Sun Yat-sen University. All patients provided written informed consent in compliance with the code of ethics of the World Medical Association (Declaration of Helsinki). The eligible patients for this study had to meet the following criteria: (1) histologically confirmed as melanoma; (2) received no radiotherapy, chemotherapy or biotherapy before surgery. The exclusion criteria were as follows: (1) previous malignancies; (2) concomitant malignancies; (3) serious active infection; and (4) pregnancy or lactation.
Eligible cutaneous melanoma patients were from The First Affiliated Hospital, Sun Yat-sen University (Guangzhou, Guangdong, China) or the Cancer Center of Guangzhou Medical University (Guangzhou, Guangdong, China). Each tumour sample was matched with adjacent apparently normal tissues removed during the same operation. Frozen sections were made from these tissues and examined by at least three pathologists. The clinicopathological features of twelve skin cutaneous melanoma patients (51.67 ± 14.57 years old) for qRT-PCR validation are shown in Table 1. Three healthy patients from The First Affiliated Hospital, Sun Yat-sen University (Guangzhou, Guangdong, China) were included in this study. These patients were scheduled to undergo split-thickness skin grafting due to deep partial burn wounds. Each normal skin sample was obtained from the donor site. All the samples were frozen immediately after the operation and were stored in liquid nitrogen until RNA isolation.
RNA isolation and qRT-PCR
Total RNA was extracted from all fresh-frozen samples using TRIzol reagent (Invitrogen, USA). The OD value (260/280) of all RNA extracted samples was greater than 1.8. For each replicate, complementary DNA (cDNA) was synthesized from 2 μg RNA using the GoScript Reverse Transcription System (Promega, USA). The qRT-PCR comprised 10 μl of GoTaq qPCR Master Mix (2×) (Promega, USA), 2 μl of diluted cDNA template (1:10) and 10 μM of each primer contributing to a total volume of 20 μl. Reactions were run in an ABI 7500 real-time PCR system (Applied Biosystems, USA) under the following conditions: 95 °C for 10 mins and 40 cycles of 95 °C for 15 s and 60 °C for 60 s. Melting curves were derived for every reaction to ensure a single product. Relative gene expression was evaluated according to the ddCT method, using the human GAPDH gene as an endogenous control for RNA load and gene expression in the analysis. All experiments were performed in triplicate. GraphPad Prism 8 (GraphPad Software, USA) was used to output figures.
The primers were as follows: MALAT1 Fw.: GACGAGTTGTGCTGCGAT; MALAT1 Rev.: TTCTGTGTTATGCCTGGTTA; LINC00943 Fw.: GGATTGGATTGTGGATTGC; LINC00943 Rev.: CAGGTCTCAGTTCAGTGTT; LINC00261 Fw.: CTTCTTGACCACATCTTACAC; LINC00261 Rev.: GGACCATTGCCTCTTGATT; GAPDH Fw: GAGAGGGAAATCGTGCGTGAC; GAPDH Rev.: CATCTGCTGGAAGGTGGACA.
Multivariate cox regression model for survival analysis
To carry out a multivariate Cox regression analysis for survival analysis of patients with MALAT1, LINC00943, and LINC00261 CNV-deficient cutaneous melanoma, we first used the UCSC genome browser (http://genome.ucsc.edu/index.html) to determine the number and region of exons of MALAT1, LINC00943, and LINC00261. All information belongs to the hg19 database (Table 2). A total of 537 SKCM patients were from the Skin Cutaneous Melanoma (TCGA, PanCancer Atlas, https://gdc.cancer.gov/about-data/publications/pancanatlas)  and Metastatic Melanoma (DFCI, Science 2015, https://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?study_id=phs000452.v2.p1) [36,37,38] datasets. Raw data were downloaded from cBioPortal (http://www.cbioportal.org/) . By further analysing the copy number variation (CNV) data of these 537 patients, we determined whether each melanoma sample had deletions of these exons. Seg. means ≤ − 0.3 were considered CNV deficiency, others were considered without CNV deficiency (see https://docs.gdc.cancer.gov/Data/Bioinformatics_Pipelines/CNV_Pipeline/, and CNV and patient information are shown in Supplementary Table 1).
To determine which factors should be included in the multivariate Cox regression model, we first performed the univariate Cox regression model for survival analysis. Factors that were statistically significant (p < 0.05) in the univariate Cox regression model were included in the multivariate Cox regression model, and the multivariate Cox regression model for survival analysis was performed. SPSS 22.0 was used for the analysis of the Cox regression model.
Identification of DEMs, DELs and DEMis and reconstruction of the lncRNA-miRNA-mRNA (ceRNA) network
After standardization of the GEO datasets, 56, 70 and 34 DEMis between benign nevus tissues and primary melanoma tissues were identified in GSE24996, GSE35579 and GSE62372, respectively (Supplementary Table 2, Fig. 2a-f). The candidate 18 miRNAs were shared in at least two datasets (Fig. 3a): hsa-miRNA-378a-3p, hsa-miRNA-23b-3p, hsa-miRNA-140-3p, hsa-miRNA-99a-5p, hsa-miRNA-100-5p, hsa-miRNA-204-5p, hsa-miRNA-211-5p, hsa-miRNA-205-5p, hsa-miRNA-224-5p, hsa-miRNA-200b-3p, hsa-miRNA-200c-3p, hsa-miRNA-125b-5p, hsa-miRNA-149-5p, hsa-miRNA-21-5p, hsa-miRNA-20b-5p, hsa-miRNA-424-5p, hsa-miRNA-203a-3p and hsa-miRNA-1826. According to method 2.3, 2361 mRNAs and 277 lncRNAs were predicted using these miRNAs. We ruled out two of these 18 DEMis, hsa-miRNA-203a-3p and hsa-miRNA-1826, because no predicted gene was found in starBase according to method 2.3. In addition, 5953 DEMs and 665 DELs between benign nevus tissues and primary melanoma tissues were identified in GSE112509 (Fig. 2g and h). As a result, a total of 898 DEMs and 53 DELs were selected for further analysis according to method 2.3 (Fig. 3b and c). Finally, 898 DEMs, 53 DELs and 16 DEMis were selected for further reconstruction of the lncRNA-miRNA-mRNA (ceRNA) network.
The lncRNA-miRNA-mRNA (ceRNA) network, consisting of 53 lncRNA nodes, 16 miRNA nodes, 898 mRNA nodes and 609 edges, was reconstructed and visualized using Cytoscape (Fig. 4a).
KEGG pathway and GO enrichment analysis of lncRNAs based on the ceRNA network
We used DAVID to analyse the biological classification of DEMs according to method 2.5. The results of the top 15 significant GO terms and KEGG pathways are shown in Table 3 and Fig. 4b-e. Sixty pathways were significantly enriched through KEGG pathway analysis, including the PI3K-Akt signalling pathway, focal adhesion, proteoglycans in cancer, pathway in cancer and, most importantly, melanomagenesis. The results of GO-BP analysis revealed 172 enriched terms, particularly in the regulation of transcription, such as positive regulation of transcription from the RNA polymerase II promoter, positive regulation of transcription (DNA-templated), and transcription from the RNA polymerase II promoter.
Hub gene selection
According to the node degree in the ceRNA network, we found that three lncRNAs, MALAT1, LINC00943, and LINC00261, had the highest number of lncRNA-miRNA and miRNA-mRNA pairs, suggesting that these three lncRNAs could be chosen as hub nodes, and the results are shown in Table 4. Therefore, these three lncRNAs might play an essential role in melanomagenesis and might be considered key lncRNAs.
Reconstruction of the MALAT1/LINC00943/LINC00261-miRNA-mRNA subnetworks
MALAT, LINC00943, LINC00261 and their paired miRNAs and mRNAs were used to reconstruct key ceRNA subnetworks. The MALAT1 ceRNA network consists of 1 lncRNA node, 9 miRNA nodes, 158 mRNA nodes and 209 edges, as shown in Fig. 5a. The LINC00943 ceRNA network consists of 1 lncRNA node, 7 miRNA nodes, 182 mRNA nodes and 209 edges, as shown in Fig. 6a. The LINC00261 ceRNA network consists of 1 lncRNA node, 5 miRNA nodes, 123 mRNA nodes and 163 edges, as shown in Fig. 7a. The results of functional analysis revealed that 75 GO-BP, 21 GO-CC, 15 GO-MF and 20 pathways were enriched in the MALAT1-miRNA-mRNA subnetwork; 67 GO-BP, 14 GO-CC, 17 GO-MF and 13 pathways were enriched in the LINC00943-miRNA-mRNA subnetwork; and 42 GO-BP, 7 GO-CC, 10 GO-MF and 7 pathways were enriched in the LINC00261-miRNA-mRNA subnetwork. The results of the top 10 significant GO terms and KEGG pathways of these three lncRNAs are shown in Fig. 5b-e, Fig. 6b-e, Fig. 7b-e, and Tables 5, 6, 7.
Expression of MALAT1, LINC00943 and LINC00261 is higher in tumour tissues
To confirm the expression of MALAT1, LINC00943 and LINC00261 in melanoma tissues, we evaluated the MALAT1, LINC00943 and LINC00261 expression levels in the cancer tissues from 12 melanoma patients (see Table 1) and 3 healthy tissues via qRT-PCR, as shown in Fig. 8. The results showed that the expression of MALAT1, LINC00943 and LINC00261 was significantly higher in the tumour tissues than in the healthy tissues (p = 0.0243, p = 0.0005, p < 0.0001, respectively). Additionally, the expression of MALAT1, LINC00943 and LINC00261 was significantly higher in the tumour tissues than in the adjacent normal tissues (p = 0.0002, p < 0.0001, p < 0.0001, respectively). However, no significant difference was observed between the healthy tissues and the adjacent normal skin tissues in the expression of MALAT1, LINC00943 and LINC00261 (p = 0.366, p = 0.379, p = 0.262, respectively). The results are consistent with those discussed above. Thus, the expression of MALAT1, LINC00943 and LINC00261 is increased in melanoma and may be responsible for the tumorigenesis of melanoma.
MALAT1 and LINC00943 are independent risk factors for the prognosis of cutaneous melanoma
A univariate Cox regression model for survival analysis of age, sex and stage was performed, and the results are shown in Supplementary Table 3. Then, the multivariate Cox regression model for survival analysis of MALAT1, LINC00943, and LINC00261 was performed. The results showed that the overall survival time and disease-free survival time of the patients with MALAT1 or LINC00943 CNV deficiency were significantly lower than those without it, and the difference was significant (details are shown in Table 8 and Fig. 9a-d), suggesting that MALAT1 and LINC00943 are independent risk factors for the prognosis of cutaneous melanoma. Although the overall survival time and disease-free survival time of patients with LINC00261 deletion were lower than those without it, the difference was not significant (p = 0.535, p = 0.694) (details are shown in Table 8 and Fig. 9e- f).
In this study, three lncRNAs, MALAT1, LINC00943 and LINC00261, were identified according to the reconstructed ceRNA network. Among these key lncRNAs found in this study, MALAT1 has been demonstrated to be related to various malignant tumours [40,41,42,43,44]. Studies have confirmed that MALAT1 is a valuable prognostic marker and a promising therapeutic target in lung cancer metastasis [40, 41]. A study also suggested that MALAT1 plays an important role in tumour progression and could serve as a promising therapeutic target . Through the study of the whole-genome mutational landscape and characterization of noncoding and structural mutations in liver cancer, Fujimoto A. and colleagues discovered that MALAT1 is closely related to liver carcinogenesis.46 In addition, a study revealed a novel mechanism of MALAT1-regulated autophagy-related chemoresistance in gastric cancer . At present, it is believed that MALAT1 is mainly responsible for regulating the proliferation, migration and invasion of tumour cells. According to our findings, MALAT1 might also be a crucial factor in the tumorigenesis and development of melanoma. In this subnetwork, we found nine lncRNA-miRNA pairs: miRNA-378a-3p, miRNA-23b-3p, miRNA-224-5p, miRNA-204-5p, miRNA-205-5p, miRNA-200c-3p, miRNA-200b-3p, miRNA-149-5p, and miRNA-211-5p. Among them, MALAT1 was shown to regulate chemoresistance via miRNA-23b-3p sequestration in gastric cancer . In ovarian cancer, a study suggested that MALAT1-miRNA-211-5p may act as a key mediator in the prevention of this disease . MALAT1 is also involved in promoting renal cell carcinoma through interaction with miRNA-205-5p . Studies have confirmed that MALAT1 functions in liver and lung cancer through miRNA-204-5p [47, 48]. In addition, targeting the MALAT1/miRNA-200c-3p axis in a xenograft endometrial carcinoma model strongly inhibited tumour growth .
Moreover, studies have illustrated that these miRNAs are closely related to melanoma in several ways. miRNA-378a-3p can regulate oncogenic PARVA expression in melanoma, preventing its progression . miRNA-23b-3p was shown to be a tumour suppressor gene in melanoma . miRNA-224-5p can be regulated by E2F1 to drive EMT through TXNIP downregulation in melanoma, and it can inhibit uveal melanoma cell proliferation, migration, and invasion by targeting PIK3R3/AKT3 [52, 53]. miRNA-204-5p, known as a tumour suppressor gene in melanoma, was associated with the CDKN2A pathway and NRAS gene and contributed to BRAF inhibitor resistance [51, 54, 55].. miRNA-205-5p suppresses proliferation and induces senescence via regulation of E2F1 in melanoma [51, 56,57,58]. miRNA-200b/c-3p act as potential diagnostic and prognostic markers for melanoma [59,60,61]. Upregulation of miRNA-149-5p, directly regulated by p53, results in increased expression of Mcl-1 and resistance to apoptosis in melanoma cells . Most importantly, studies have confirmed that miRNA-211-5p plays a major role as a tumour suppressor via various targets in melanoma [51, 55, 59, 63, 64]. Moreover, MALAT1 is an independent risk factor for the prognosis of SKCM according to multivariate Cox regression model analysis. Thus, we believe that MALAT1 may contribute to the tumorigenesis and survival of SKCM.
Little is known about LINC00943. According to the LINC00943-miRNA-mRNA subnetwork, miRNA-99a-5p, miRNA-100-5p, miRNA-23b-3p, miRNA-204-5p, miRNA-224-5p, miRNA-149-5p and miRNA-125b-5p closely interacted with LINC00943. No connection between LINC00943 and these miRNAs has been discovered yet; however, these miRNAs were also demonstrated to be associated with melanoma, except miRNA-99a-5p. The links between miRNA-204-5p, miRNA-224-5p, miRNA-149-5p and melanoma are discussed above. In addition, miRNA-23b was suggested as a tumour suppressor gene.54 miRNA-100-5p and miRNA-125b-5p are associated with resistance to treatment with immune checkpoint inhibitors in melanoma . Additionally, we confirmed that LINC00943 is an independent risk factor for the prognosis of SKCM. Therefore, understanding the relationships among LINC00943, miRNAs and malignancies may provide further information for future research on melanoma and other malignancies.
Seven KEGG pathways were enriched based on the LINC00261 subnetwork. One of these pathways, the PI3K/Akt signalling pathway, has been proven to play a critical role in tumorigenesis , especially in melanoma . Additionally, a study has demonstrated that LINC00261 promotes cancer cell proliferation and metastasis in human choriocarcinoma . However, LINC00261 has shown a strong capacity in improving the chemotherapeutic response and survival of patients with oesophageal cancer . In gastric cancer, LINC00261 can suppress tumour metastasis by regulating epithelial-mesenchymal transition . Moreover, LINC00261 can block cellular proliferation by activating the DNA damage response . LINC00261 may affect the biological behaviour of different tumours in different ways. Therefore, it is essential to further explore the role of LINC00261 in different tumours. However, five miRNAs, miRNA-23b-3p, miRNA-211-5p, miRNA-205-5p, miRNA-140-3p and miRNA-125b-5p, interacted with LINC00261 according to the LINC00261-miRNA-mRNA subnetwork. Similarly, no connection between LINC00261 and these miRNAs has been discovered yet. The roles of miRNA-23b-3p, miRNA-211-5p, miRNA-205-5p, and miRNA-125b-5p in melanoma are discussed above. miRNA-140-3p was reported to be regulated by MALAT1 in uveal melanoma cells . The multivariate Cox regression model for survival suggested that LINC00261 was not a risk factor for the prognosis of SKCM, however, the median overall survival and disease-free survival time for patients with LINC00261 CNV deficiency were significantly lower than those without LINC00261 CNV deficiency (17.03 m vs 61.05 m, 13.50 vs 25.02).
Three of the 16 predicted miRNAs were not associated with MALAT1, LINC00943 and LINC00261: miRNA-21-5p, miRNA-20b-5p and miRNA-424-5p. They are closely related to SGMS1.AS1, EPB41L4A.AS1 and SNHG1 according to the ceRNA network. Little is known about miRNA-424-5p in melanoma, while studies have suggested that miRNA-20b-5p may inhibit tumour metastasis via regulation of the PAR-1 receptor in melanoma cells , and miRNA-21 may regulate melanoma cell proliferation, migration, and apoptosis through the ERK/NF-κB signalling pathway by targeting SPRY1, PDCD4 and PTEN [74, 75].
This study advances our understanding of tumorigenesis and development in cutaneous melanoma from the perspective of the ceRNA theory. In addition, MALAT1 and LINC00943 may be independent risk factors for the prognosis of patients with cutaneous melanoma and might become predictive molecules for the long-term treatment of melanoma and potential therapeutic targets. Further studies are required to validate the role of MALAT1, LINC00943 and LINC00261 in cutaneous melanoma.
Availability of data and materials
The data that support results of the present study are available from GEO datasets (including GSE24996(https://www.ncbi.nlm.nih.gov/ geo/query/acc.cgi?acc = GSE24996),GSE35579(https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE35579),GSE62372(https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE62372), and GSE112509(https://www.ncbi.nlm.nih.gov/geo/query/ acc.cgi?acc = GSE112509)), and cBioportal (http://www.cbioportal.org/), DAVID (https://david.ncifcrf.gov/), and starbase miRNA-mRNA Interactions (http://starbase.sysu.edu.cn/agoClipRNA.php?source=mRNA), and starbase miRNA-lncRNA Interactions (http://starbase.sysu.edu.cn/agoClipRNA. php?source = lncRNA) database.
Competitive endogenous RNA
Skin cutaneous melanoma
Long non-coding RNAs
- NCBI GEO:
National center for biotechnology information gene expression omnibus
Kyoto encyclopedia of genes and genomes
Database for annotation, visualization, and integration discovery
Differential expressed miRNAs
Differential expressed lncRNAs
Differential expressed mRNAs
Copy number variation
McGuire S. World Cancer report 2014. Geneva, Switzerland: World Health Organization, International Agency for Research on Cancer, WHO press, 2015. Adv Nutr. 2016;7(2):418–9.
Berwick M, Erdei E, Hay J. Melanoma epidemiology and public health. Dermatol Clin. 2009;27(2):205–14 viii.
Schadendorf D, van Akkooi ACJ, Berking C, Griewank KG, Gutzmer R, Hauschild A, Stang A, Roesch A, Ugurel S. Melanoma. Lancet. 2018;392(10151):971–84.
Disease GBD, Injury I, Prevalence C. Global, regional, and national incidence, prevalence, and years lived with disability for 310 diseases and injuries, 1990-2015: a systematic analysis for the global burden of Disease study 2015. Lancet. 2016;388(10053):1545–602.
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2020. CA Cancer J Clin. 2020;70(1):7–30.
Burrell RA, McGranahan N, Bartek J, Swanton C. The causes and consequences of genetic heterogeneity in cancer evolution. Nature. 2013;501(7467):338–45.
Lander ES, Linton LM, Birren B, Nusbaum C, Zody MC, Baldwin J, Devon K, Dewar K, Doyle M, FitzHugh W, et al. Initial sequencing and analysis of the human genome. Nature. 2001;409(6822):860–921.
Yost SE, Smith EN, Schwab RB, Bao L, Jung H, Wang X, Voest E, Pierce JP, Messer K, Parker BA, et al. Identification of high-confidence somatic mutations in whole genome sequence of formalin-fixed breast cancer specimens. Nucleic Acids Res. 2012;40(14):e107.
Goodrich JA, Kugel JF. Non-coding-RNA regulators of RNA polymerase II transcription. Nat Rev Mol Cell Biol. 2006;7(8):612–6.
Yoon JH, Abdelmohsen K, Gorospe M. Posttranscriptional gene regulation by long noncoding RNA. J Mol Biol. 2013;425(19):3723–30.
Kiefer JC. Epigenetics in development. Dev Dyn. 2007;236(4):1144–56.
Mikkelsen TS, Ku M, Jaffe DB, Issac B, Lieberman E, Giannoukos G, Alvarez P, Brockman W, Kim TK, Koche RP, et al. Genome-wide maps of chromatin state in pluripotent and lineage-committed cells. Nature. 2007;448(7153):553–60.
Joung J, Engreitz JM, Konermann S, Abudayyeh OO, Verdine VK, Aguet F, Gootenberg JS, Sanjana NE, Wright JB, Fulco CP, et al. Genome-scale activation screen identifies a lncRNA locus regulating a gene neighbourhood. Nature. 2017;548(7667):343–6.
Leucci E, Vendramin R, Spinazzi M, Laurette P, Fiers M, Wouters J, Radaelli E, Eyckerman S, Leonelli C, Vanderheyden K, et al. Melanoma addiction to the long non-coding RNA SAMMSON. Nature. 2016;531(7595):518–22.
Hosono Y, Niknafs YS, Prensner JR, Iyer MK, Dhanasekaran SM, Mehra R, Pitchiaya S, Tien J, Escara-Wilke J, Poliakov A, et al. Oncogenic role of THOR, a conserved Cancer/testis long non-coding RNA. Cell. 2017;171(7):1559–72 e1520.
Montes M, Nielsen MM, Maglieri G, Jacobsen A, Hojfeldt J, Agrawal-Singh S, Hansen K, Helin K, van de Werken HJG, Pedersen JS, et al. The lncRNA MIR31HG regulates p16(INK4A) expression to modulate senescence. Nat Commun. 2015;6:6967.
Li P, He J, Yang Z, Ge S, Zhang H, Zhong Q, Fan X. ZNNT1 long noncoding RNA induces autophagy to inhibit tumorigenesis of uveal melanoma by regulating key autophagy gene expression. Autophagy. 2020;16(7):1186–99.
Jalali S, Bhartiya D, Lalwani MK, Sivasubbu S, Scaria V. Systematic transcriptome wide analysis of lncRNA-miRNA interactions. PLoS One. 2013;8(2):e53823.
Ala U, Karreth FA, Bosia C, Pagnani A, Taulli R, Leopold V, Tay Y, Provero P, Zecchina R, Pandolfi PP. Integrated transcriptional and competitive endogenous RNA networks are cross-regulated in permissive molecular environments. Proc Natl Acad Sci U S A. 2013;110(18):7154–9.
Salmena L, Poliseno L, Tay Y, Kats L, Pandolfi PP. A ceRNA hypothesis: the Rosetta stone of a hidden RNA language? Cell. 2011;146(3):353–8.
Tay Y, Rinn J, Pandolfi PP. The multilayered complexity of ceRNA crosstalk and competition. Nature. 2014;505(7483):344–52.
Rinn JL, Chang HY. Genome regulation by long noncoding RNAs. Annu Rev Biochem. 2012;81:145–66.
Liu XH, Sun M, Nie FQ, Ge YB, Zhang EB, Yin DD, Kong R, Xia R, Lu KH, Li JH, et al. Lnc RNA HOTAIR functions as a competing endogenous RNA to regulate HER2 expression by sponging miR-331-3p in gastric cancer. Mol Cancer. 2014;13:92.
Zhong Z, Huang M, Lv M, He Y, Duan C, Zhang L, Chen J. Circular RNA MYLK as a competing endogenous RNA promotes bladder cancer progression through modulating VEGFA/VEGFR2 signaling pathway. Cancer Lett. 2017;403:305–17.
Chang L, Guo R, Yuan Z, Shi H, Zhang D. LncRNA HOTAIR regulates CCND1 and CCND2 expression by sponging miR-206 in ovarian Cancer. Cell Physiol Biochem. 2018;49(4):1289–303.
Edgar R, Domrachev M, Lash AE. Gene expression omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002;30(1):207–10.
Breuer J: R (Software). In., edn.; 2017.
Smyth GK: limma: Linear Models for Microarray Data. 2005.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.
Li JH, Liu S, Zhou H, Qu LH, Yang JH. starBase v2.0: decoding miRNA-ceRNA, miRNA-ncRNA and protein-RNA interaction networks from large-scale CLIP-Seq data. Nucleic Acids Res. 2014;42(Database issue):D92–7.
Liao Q, Liu C, Yuan X, Kang S, Miao R, Xiao H, Zhao G, Luo H, Bu D, Zhao H, et al. Large-scale prediction of long non-coding RNA functions in a coding-non-coding gene co-expression network. Nucleic Acids Res. 2011;39(9):3864–78.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.
Huang da W, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009;37(1):1–13.
Huang da W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.
Hoadley KA, Yau C, Hinoue T, Wolf DM, Lazar AJ, Drill E, Shen R, Taylor AM, Cherniack AD, Thorsson V, et al. Cell-of-origin patterns dominate the molecular classification of 10,000 tumors from 33 types of Cancer. Cell. 2018;173(2):291–304 e296.
Van Allen EM, Miao D, Schilling B, Shukla SA, Blank C, Zimmer L, Sucker A, Hillen U, Foppen MHG, Goldinger SM, et al. Genomic correlates of response to CTLA-4 blockade in metastatic melanoma. Science. 2015;350(6257):207–11.
Tryka KA, Hao L, Sturcke A, Jin Y, Wang ZY, Ziyabari L, Lee M, Popova N, Sharopova N, Kimura M, et al. NCBI's database of genotypes and phenotypes: dbGaP. Nucleic Acids Res. 2014;42(Database issue):D975–9.
Mailman MD, Feolo M, Jin Y, Kimura M, Tryka K, Bagoutdinov R, Hao L, Kiang A, Paschall J, Phan L, et al. The NCBI dbGaP database of genotypes and phenotypes. Nat Genet. 2007;39(10):1181–6.
Gao J, Aksoy BA, Dogrusoz U, Dresdner G, Gross B, Sumer SO, Sun Y, Jacobsen A, Sinha R, Larsson E, et al. Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Sci Signal. 2013;6(269):pl1.
Gutschner T, Hammerle M, Eissmann M, Hsu J, Kim Y, Hung G, Revenko A, Arun G, Stentrup M, Gross M, et al. The noncoding RNA MALAT1 is a critical regulator of the metastasis phenotype of lung cancer cells. Cancer Res. 2013;73(3):1180–9.
Gutschner T, Hammerle M, Diederichs S. MALAT1 -- a paradigm for long noncoding RNA function in cancer. J Mol Med (Berl). 2013;91(7):791–801.
Lai MC, Yang Z, Zhou L, Zhu QQ, Xie HY, Zhang F, Wu LM, Chen LM, Zheng SS. Long non-coding RNA MALAT-1 overexpression predicts tumor recurrence of hepatocellular carcinoma after liver transplantation. Med Oncol. 2012;29(3):1810–6.
Fujimoto A, Furuta M, Totoki Y, Tsunoda T, Kato M, Shiraishi Y, Tanaka H, Taniguchi H, Kawakami Y, Ueno M, et al. Whole-genome mutational landscape and characterization of noncoding and structural mutations in liver cancer. Nat Genet. 2016;48(5):500–9.
YiRen H, YingCong Y, Sunwu Y, Keqin L, Xiaochun T, Senrui C, Ende C, XiZhou L, Yanfan C. Long noncoding RNA MALAT1 regulates autophagy associated chemoresistance via miR-23b-3p sequestration in gastric cancer. Mol Cancer. 2017;16(1):174.
Tao F, Tian X, Ruan S, Shen M, Zhang Z. miR-211 sponges lncRNA MALAT1 to suppress tumor growth and progression through inhibiting PHF19 in ovarian carcinoma. FASEB J. 2018;32:fj.201800495RR.
Hirata H, Hinoda Y, Shahryari V, Deng G, Nakajima K, Tabatabai ZL, Ishii N, Dahiya R. Long noncoding RNA MALAT1 promotes aggressive renal cell carcinoma through Ezh2 and interacts with miR-205. Cancer Res. 2015;75(7):1322–31.
Li J, Wang J, Chen Y, Li S, Jin M, Wang H, Chen Z, Yu W. LncRNA MALAT1 exerts oncogenic functions in lung adenocarcinoma by targeting miR-204. Am J Cancer Res. 2016;6(5):1099–107.
Tan X, Huang Z, Li X. Long non-coding RNA MALAT1 interacts with miR-204 to modulate human Hilar Cholangiocarcinoma proliferation, migration, and invasion by targeting CXCR4. J Cell Biochem. 2017;118(11):3643–53.
Li Q, Zhang C, Chen R, Xiong H, Qiu F, Liu S, Zhang M, Wang F, Wang Y, Zhou X, et al. Disrupting MALAT1/miR-200c sponge decreases invasion and migration in endometrioid endometrial carcinoma. Cancer Lett. 2016;383(1):28–40.
Velazquez-Torres G, Shoshan E, Ivan C, Huang L, Fuentes-Mattei E, Paret H, Kim SJ, Rodriguez-Aguayo C, Xie V, Brooks D, et al. A-to-I miR-378a-3p editing can prevent melanoma progression via regulation of PARVA expression. Nat Commun. 2018;9(1):461.
Kozubek J, Ma Z, Fleming E, Duggan T, Wu R, Shin DG, Dadras SS. In-depth characterization of microRNA transcriptome in melanoma. PLoS One. 2013;8(9):e72699.
Knoll S, Furst K, Kowtharapu B, Schmitz U, Marquardt S, Wolkenhauer O, Martin H, Putzer BM. E2F1 induces miR-224/452 expression to drive EMT through TXNIP downregulation. EMBO Rep. 2014;15(12):1315–29.
Li J, Liu X, Li C, Wang W. miR-224-5p inhibits proliferation, migration, and invasion by targeting PIK3R3/AKT3 in uveal melanoma. J Cell Biochem. 2019;120(8):12412–21.
Galasso M, Morrison C, Minotti L, Corra F, Zerbinati C, Agnoletto C, Baldassari F, Fassan M, Bartolazzi A, Vecchione A, et al. Loss of miR-204 expression is a key event in melanoma. Mol Cancer. 2018;17(1):71.
Diaz-Martinez M, Benito-Jardon L, Alonso L, Koetz-Ploch L, Hernando E, Teixido J. miR-204-5p and miR-211-5p contribute to BRAF inhibitor resistance in melanoma. Cancer Res. 2018;78(4):1017–30.
Xu Y, Brenn T, Brown ER, Doherty V, Melton DW. Differential expression of microRNAs during melanoma progression: miR-200c, miR-205 and miR-211 are downregulated in melanoma and act as tumour suppressors. Br J Cancer. 2012;106(3):553–61.
Dar AA, Majid S, de Semir D, Nosrati M, Bezrookove V. Kashani-Sabet M: miRNA-205 suppresses melanoma cell proliferation and induces senescence via regulation of E2F1 protein. J Biol Chem. 2011;286(19):16606–14.
Sanchez-Sendra B, Martinez-Ciarpaglini C, Gonzalez-Munoz JF, Murgui A, Terradez L, Monteagudo C. Downregulation of intratumoral expression of miR-205, miR-200c and miR-125b in primary human cutaneous melanomas predicts shorter survival. Sci Rep. 2018;8(1):17076.
Mirzaei H, Gholamin S, Shahidsales S, Sahebkar A, Jaafari MR, Mirzaei HR, Hassanian SM, Avan A. MicroRNAs as potential diagnostic and prognostic biomarkers in melanoma. Eur J Cancer. 2016;53:25–32.
Elson-Schwab I, Lorentzen A, Marshall CJ. MicroRNA-200 family members differentially regulate morphological plasticity and mode of melanoma cell invasion. PLoS One. 2010;5(10):e13176.
Zhao H, Xing G, Wang Y, Luo Z, Liu G, Meng H. Long noncoding RNA HEIH promotes melanoma cell proliferation, migration and invasion via inhibiting miR-200b/a/429. Biosci Rep. 2017;37:BSR20170682.
Jin L, Hu WL, Jiang CC, Wang JX, Han CC, Chu P, Zhang LJ, Thorne RF, Wilmott J, Scolyer RA, et al. MicroRNA-149*, a p53-responsive microRNA, functions as an oncogenic regulator in human melanoma. Proc Natl Acad Sci U S A. 2011;108(38):15840–5.
Bell RE, Khaled M, Netanely D, Schubert S, Golan T, Buxbaum A, Janas MM, Postolsky B, Goldberg MS, Shamir R, et al. Transcription factor/microRNA axis blocks melanoma invasion program by miR-211 targeting NUAK1. J Invest Dermatol. 2014;134(2):441–51.
Levy C, Khaled M, Iliopoulos D, Janas MM, Schubert S, Pinner S, Chen PH, Li S, Fletcher AL, Yokoyama S, et al. Intronic miR-211 assumes the tumor suppressive function of its host gene in melanoma. Mol Cell. 2010;40(5):841–9.
Huber V, Vallacchi V, Fleming V, Hu X, Cova A, Dugo M, Shahaj E, Sulsenti R, Vergani E, Filipazzi P, et al. Tumor-derived microRNAs induce myeloid suppressor cells and predict immunotherapy resistance in melanoma. J Clin Invest. 2018;128(12):5505–16.
Fruman DA, Rommel C. PI3K and cancer: lessons, challenges and opportunities. Nat Rev Drug Discov. 2014;13(2):140–56.
Davies MA. The role of the PI3K-AKT pathway in melanoma. Cancer J. 2012;18(2):142–7.
Wang Y, Xue K, Guan Y, Jin Y, Liu S, Wang Y, Liu S, Wang L, Han L. Long noncoding RNA LINC00261 suppresses cell proliferation and invasion and promotes cell apoptosis in human Choriocarcinoma. Oncol Res. 2017;25(5):733–42.
Lin K, Jiang H, Zhuang SS, Qin YS, Qiu GD, She YQ, Zheng JT, Chen C, Fang L, Zhang SY. Long noncoding RNA LINC00261 induces chemosensitization to 5-fluorouracil by mediating methylation-dependent repression of DPYD in human esophageal cancer. FASEB J. 2019;33(2):1972–88.
Yu Y, Li L, Zheng Z, Chen S, Chen E, Hu Y. Long non-coding RNA linc00261 suppresses gastric cancer progression via promoting slug degradation. J Cell Mol Med. 2017;21(5):955–67.
Shahabi S, Kumaran V, Castillo J, Cong Z, Nandagopal G, Mullen DJ, Alvarado A, Correa MR, Saizan A, Goel R, et al. LINC00261 is an epigenetically regulated tumor suppressor essential for activation of the DNA damage response. Cancer Res. 2019;79(12):3050–62.
Sun L, Sun P, Zhou QY, Gao X, Han Q. Long noncoding RNA MALAT1 promotes uveal melanoma cell growth and invasion by silencing of miR-140. Am J Transl Res. 2016;8(9):3939–46.
Saleiban A, Faxalv L, Claesson K, Jonsson JI, Osman A. miR-20b regulates expression of proteinase-activated receptor-1 (PAR-1) thrombin receptor in melanoma cells. Pigment Cell Melanoma Res. 2014;27(3):431–41.
Mao XH, Chen M, Wang Y, Cui PG, Liu SB, Xu ZY. MicroRNA-21 regulates the ERK/NF-kappaB signaling pathway to affect the proliferation, migration, and apoptosis of human melanoma A375 cells by targeting SPRY1, PDCD4, and PTEN. Mol Carcinog. 2017;56(3):886–94.
Yang CH, Yue J, Pfeffer SR, Handorf CR, Pfeffer LM. MicroRNA miR-21 regulates the metastatic behavior of B16 melanoma cells. J Biol Chem. 2011;286(45):39172–8.
This manuscript is approved by all authors for publication. And we sincerely thank those who help finishing this article.
This article is funded by the Science and Technology Program of Guangzhou (201704020165) and Natural Science Foundation of Guangdong Province (2017A030313619).
Ethics approval and consent to participate
The study protocol was approved by the Ethics Committee of The First Affiliated Hospital, Sun Yat-sen University. All patients provided written informed consent in compliance with the code of ethics of World Medical Association (Declaration of Helsinki). All the data used in this article were from open assess databases, and no permission was required.
Consent for publication
Written informed consent for publication of this article was obtained from patients themselves.
No potential conflict of interest was declared by the authors.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Differentially expressed miRNAs in GSE24996、GSE35579、GSE62372.
Univariate COX regression model for survival analysis of age, sex and stage.
About this article
Cite this article
Zhu, J., Deng, J., Zhang, L. et al. Reconstruction of lncRNA-miRNA-mRNA network based on competitive endogenous RNA reveals functional lncRNAs in skin cutaneous melanoma. BMC Cancer 20, 927 (2020). https://doi.org/10.1186/s12885-020-07302-5
- Human skin cutaneous melanoma
- Competitive endogenous RNA