Hsa-let-7c-5p, hsa-miR-130b-3p, and hsa-miR-142-3p as Novel miRNA Biomarkers for Melanoma Progression

Wu, Xuerui; Wang, Yu; Chen, Chen; Xue, Yadong; Zheng, Shuyun; Cai, Limin

doi:https://doi.org/10.1155/2022/5671562

Genetics Research

On this page

Abstract Introduction Materials and Methods Results Discussion Conclusions Data Availability Ethical Approval Disclosure Conflicts of Interest Authors’ Contributions Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2022 | Article ID 5671562 | https://doi.org/10.1155/2022/5671562

Hsa-let-7c-5p, hsa-miR-130b-3p, and hsa-miR-142-3p as Novel miRNA Biomarkers for Melanoma Progression

Xuerui Wu,¹Yu Wang,¹Chen Chen,¹Yadong Xue,¹Shuyun Zheng,¹and Limin Cai¹

Academic Editor: Muhammad Babar Khawar

Received08 Apr 2022

Revised31 May 2022

Accepted15 Jun 2022

Published07 Jul 2022

Abstract

This study aimed to screen miRNA biomarkers for melanoma progression. Raw melanoma data were downloaded from the Gene Expression Omnibus (GSE34460, GSE35579, GSE18509, and GSE24996) and the Cancer Genome Atlas (TCGA). Then, all differentially expressed miRNAs (DEmiRNAs) between benign vs. primary, metastatic vs. benign, and metastatic vs. primary groups were obtained in the GSE34460 and GSE35579 datasets, and the miRNAs related to disease progression were further screened. Then, the miRNA-gene network was constructed, followed by enrichment, survival, and cluster analyses. Differentially expressed genes (DEGs), tumor-infiltrating immune cells, and tumor mutation burden (TMB) between subtypes were analyzed. miRNAs were verified in the GSE18509 and GSE24996 datasets. A total of 132 and 209 DEmiRNAs were obtained in the GSE34460 and GSE35579 datasets, respectively, and 27 DEmiRNAs related to disease progression were screened. hsa-miR-106b-5p, hsa-miR-27b-3p, and hsa-miR-141-3p had a higher degree and were regulated by numerous genes in the miRNA-gene network. Moreover, four miRNAs were associated with prognosis: hsa-let-7c-5p, hsa-miR-130b-3p, hsa-miR-142-3p, and hsa-miR-509-3p. Furthermore, the bidirectional hierarchical clustering of 27 miRNAs was classified into three subtypes, and TMB and four types of immune cells, including activated dendritic cells, naïve CD4 T cells, M1 macrophages, and plasma cells, showed significant differences among the three subtypes. The expression levels of most miRNAs in the GSE18509 and GSE24996 datasets were consistent with those in the training dataset. These miRNAs, including hsa-let-7c-5p, hsa-miR-130b-3p, and hsa-miR-142-3p, and activated dendritic cells, naïve CD4 T cells, M1 macrophages, and plasma cells may play vital roles in the pathogenesis of melanoma.

1. Introduction

Melanoma is a highly aggressive malignant tumor originating from neural crest melanocytes [1, 2]. Melanoma cells can infiltrate the tissues, lymph, and blood vessels. Melanoma accounts for 3% of all malignant tumors of the skin [3], with its incidence rate increasing rapidly worldwide [4, 5]. In 2017, the incidence rate of melanoma in China was 0.9/10 million by age [6]. Melanoma is not sensitive to radiotherapy, chemotherapy, and biological immunotherapy, and its invasiveness and early metastasis make it the deadliest among skin malignant tumors, with only a 6% five-year survival rate of stage IV patients [7]. Thus, it is necessary to screen for reliable and effective biomarkers for the diagnosis and prognosis of melanoma.

MicroRNAs (miRNAs) are gene regulatory factors with excellent functions, ranging in size from 19 to 25 nucleotides [8]. miRNAs participate in the posttranscriptional control of gene expression by binding to the 3′-untranslated region (3′-UTR) of the target mRNA and then inhibits or degrades the target gene [9]. This molecular-level imbalance in biological processes is a potential mechanism for the development of many diseases. Numerous miRNAs have been reportedly involved in melanoma development. Lu et al. found that five miRNAs, including miR-25, miR-204, miR-211, miR-510, and miR-513c, could be used as prognostic biomarkers for patients with cutaneous melanoma cancer [10]. Zhu et al. revealed that miRNA-378a-3p, miRNA-23b-3p, and miRNA-200c-3p are associated with the occurrence of skin cutaneous melanoma [11]. Li et al. suggested that miR-300 and miR-629 are associated with melanoma [12]. However, these studies screened the miRNAs related to melanoma based on a single dataset with a small sample size, and few studies have identified the miRNAs involved in the progression of melanoma. Therefore, miRNA biomarkers for melanoma progression should be identified based on multiple datasets.

This study aimed to screen miRNA biomarkers for melanoma progression. Using the data downloaded from the Gene Expression Omnibus (GEO) and the Cancer Genome Atlas (TCGA) databases, miRNAs related to melanoma progression were obtained, the miRNA-gene network was constructed, and prognosis-related miRNAs were screened. Moreover, cluster analysis was conducted based on miRNAs related to melanoma progression, and the differentially expressed genes (DEGs), tumor-infiltrating immune cells, and tumor mutation burden (TMB) between subtypes were analyzed. This study provided a powerful basis for identifying novel therapeutic targets for melanoma patients. A flowchart of this study is shown in Figure 1.

2. Materials and Methods

2.1. Data Collection

Four datasets that included some patients with melanoma were included from the GEO database [13], including GSE34460 (8 benign melanoma, 9 primary melanoma, and 4 metastatic melanoma, detection platform: GPL15019 Agilent-031181 Unrestricted_Human_miRNA_V16.0_Microarray 030840 (miRBase release 14.0 miRNA ID version)), GSE35579 (11 benign melanoma, 20 primary melanoma, and 21 metastatic melanoma, detection platform: GPL15183 CRUK/Melton lab-Human melanoma-71-v2-microRNA expression profiling), GSE18509 (8 benign melanoma, 8 metastatic melanoma, detection platform: GPL9081 Agilent-016436 Human miRNA Microarray 1.0 G4472A (miRNA ID version)), and GSE24996 (8 benign melanoma, 15 primary melanoma, and 8 metastatic melanoma, detection platform: GPL6955 Agilent-016436 Human miRNA Microarray 1.0 (feature number version)). Moreover, the miRNA mature strand expression RNAseq data (log2 (RPM+1)), gene expression RNAseq data (log2 (norm_count+1)), and the survival information of the corresponding samples (overall survival (OS), disease-specific survival (DSS), and progression-free interval (PFI)) of skin melanoma (SKCM) were obtained from TCGA of UCSC-Xena [14]. The annotated MAF file processed by the MuTect software of somatic mutation data (SNPs and small INDELs) of TCGA-SKCM was obtained from the Genomic Data Commons (GDC) database (https://portal.gdc.cancer.gov) for subsequent mutation correlation analysis.

2.2. Data Processing

For the GSE34460, GSE18509, and GSE24996 datasets, the raw data were preprocessed using limma [15] in the R package. The data preprocessing processes included data reading, RMA background adjustment, quantile normalization, and normalization. The expression matrix of probes was obtained, and the probes were annotated to miRNA ID according to the platform annotation information; then, the miRNA ID was converted using miRBase website [16], and the obtained miRNAs were used for the subsequent analysis. For the GSE35579 dataset, the processed probe expression matrix was downloaded, which was log-transformed, and the probes were annotated to miRNA ID according to the platform annotation information. Then, the miRNA ID was converted using miRBase website, and the obtained miRNAs were used for the subsequent analysis. For the TCGA-SKCM dataset, samples with an OS. time > 0 were selected for subsequent cluster analysis.

2.3. Differential Expression Analysis

Based on the GSE34460 and GSE35579 datasets, the typical Bayesian test provided by the R package limma package was used to screen the differentially expressed miRNAs (DEmiRNAs) between benign vs. primary, metastatic vs. benign, and metastatic vs. primary, with a cutoff value of and |logFC| > 0.5. Then, all overlapping DEmiRNAs screened between benign vs. primary, metastatic vs. benign, and metastatic vs. primary were obtained in each dataset, and the common DEmiRNAs were obtained from the two datasets. Based on the common DEmiRNAs, the expressed values of the common DEmiRNAs in each sample of the GSE34460 and GSE35579 datasets were extracted. For the same group of samples, the average of each miRNA was calculated as the average level of the miRNA in the sample group. After comparing the expression levels of DEmiRNAs in the benign, primary, and metastatic groups, the DEmiRNAs were divided into four types: expression continuous upregulation (up-up), expression continuous downregulation (down-down), upregulated first and then downregulated (up-down), and downregulated first and then upregulated (down-up). The miRNAs that overlapped in the same types were considered miRNAs related to disease progression for subsequent analysis.

2.4. Prediction of Target Genes of miRNA and Enrichment Analysis

Based on the six databases, miRWalk, Microt4, miRanda, miRDB, RNA22, and TargetScan, the miRNA-gene interaction pairs were screened using the miRWalk 2.0 tool [17]. The miRNA-gene interaction pairs overlapped in the six databases were identified. Moreover, the “Skin Cutaneous Melanoma” was used as the keyword to screen the disease-related genes in the GeneCards database [18]. The disease-related genes were intersected with target genes of miRNA that overlapped in the six databases, and the miRNA-gene interaction pairs were further obtained, followed by visualization using Cytoscape [19]. Additionally, clusterProfiler [20] in the R package was used to perform enrichment analysis of the genes. The value was corrected using Benjamini and Hochberg (BH), and the adjust was regarded as a statistically significant difference.

2.5. Survival Analysis

Based on the miRNAs related to disease progression, the expression values of the miRNAs in each sample were extracted from the miRNA expression data in TCGA-SKCM samples. The patients were then allocated to high and low-miRNA expression groups according to the median expression value, and the Kaplan–Meier curve method was used to evaluate the association between miRNAs and prognosis.

2.6. Cluster Analysis of miRNAs

Based on the miRNAs related to disease progression and the expression value of miRNAs in each sample in TCGA-SKCM, bidirectional hierarchical clustering was conducted based on the centered Pearson correlation algorithm using pheatmap [21] to identify different subtypes. The Kaplan–Meier curve method was used to evaluate the association between subtypes and prognosis.

2.7. Identification of DEGs between Subtypes

The gene expression matrix corresponding to each subtype sample was extracted, and the DEGs were identified between the two subtypes using limma in the R package. The value was corrected using BH, and the adjust and |logFC| > 0.5 were set as the cutoff values. Moreover, the DAVID tool [22] was used to perform enrichment analysis on the common DEGs with a threshold value of and count ≥ 2.

2.8. Tumor-Infiltrating Immune Cells

Based on the gene expression level of samples in each subtype, the CIBERSORT algorithm [23] was utilized to evaluate the proportion of immune cells in subtypes, and the differences in the proportion of immune cells in subtypes were analyzed using the Kruskal–Wallis test with a threshold value of .

2.9. TMB Analysis

The Oncoplot function in Maftools [24] was used to visualize the top 20 genes with the highest mutation frequency in each subtype based on the genetic mutation information in each sample. The TMB value was calculated, and the differences in TMB values in subtypes were analyzed using the t-test with a threshold value of .

2.10. Validation Analysis

To verify that the miRNAs were related to melanoma, the GSE18509 and GSE24996 datasets were used as the validation datasets. The expression values of miRNAs in each sample in the GSE18509 and GSE24996 datasets were extracted, and the differences in expression levels between benign and metastatic groups were analyzed using the t-test.

3. Results

3.1. Identification of DEmiRNAs

Based on the GSE34460 and GSE35579 datasets, DEmiRNAs between benign vs. primary, metastatic vs. benign, and metastatic vs. primary were identified. As given in Table 1, 80, 107, and 24 DEmiRNAs were obtained between benign vs. primary, metastatic vs. benign, and metastatic vs. primary in the GSE34460 dataset, respectively, and 138, 127, and 78 DEmiRNAs were screened between benign vs. primary, metastatic vs. benign, and metastatic vs. primary in the GSE35579 dataset, respectively. Thus, a total of 132 and 209 overlapping DEmiRNAs were obtained in the GSE34460 and GSE35579 datasets, respectively. Then, totally 61 common DEmiRNAs were screened (Figure 2(a)). A total of 27 DEmiRNAs that overlapped in the same types were considered as miRNAs related to disease progression, including 18 and 9 miRNAs in down-down- and up-up types, respectively (Figure 2(b)). Heatmaps of the 27 DEmiRNAs in the GSE34460 and GSE35579 datasets are shown in Figures 2(c) and 2(d), respectively.

(a)

(b)

(c)

(d)

3.2. Prediction of Target Genes of miRNA and Enrichment Analysis

A total of 2330 miRNA-gene interaction pairs containing 26 miRNAs and 1490 genes were screened using miRWalk 2.0. Then, the disease-related genes screened in the GeneCards database were intersected with target genes of miRNA, and a total of 604 miRNA-gene interaction pairs were further obtained, including 26 miRNAs and 377 genes. The miRNA-gene network is shown in Figure 3(a). hsa-miR-106b-5p, hsa-miR-27b-3p, and hsa-miR-141-3p had higher degree and were regulated by numerous genes (Table 2). The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis showed that the target genes of 14 miRNAs were involved in the MAPK and PI3K-Akt signaling pathways, and the top5 enriched pathways are shown in Figure 3(b).

(a)

(b)

3.3. Survival Analysis

To evaluate the relationship between the 27 miRNAs and prognosis, survival analysis was conducted using the Kaplan–Meier curve method, and the results showed that four miRNAs were associated with OS (hsa-let-7c-5p, hsa-miR-130b-3p, hsa-miR-142-3p, and hsa-miR-509-3p) (all , Figure 4(a)), and three miRNAs were correlated with DSS (hsa-let-7c-5p, hsa-miR-130b-3p, and hsa-miR-142-3p) (all , Figure 4(b)); however, there were no miRNAs related to PFI.

(a)

(b)

3.4. Cluster Analysis of miRNAs

Bidirectional hierarchical clustering was conducted based on the expression value of the 27 miRNAs in each sample in TCGA-SKCM. As shown in Figure 5(a), the bidirectional hierarchical clustering of the 27 miRNAs was classified into three subtypes (clusters 1, 2, and 3), and most miRNAs were highly expressed in cluster 2; however, miRNAs showed low expression in clusters 1 and 3. The principal component analysis (PCA) of three subtypes is shown in Figure 5(b). Moreover, survival analysis showed that clusters 1 and 2 were associated with poor OS, and cluster 3 correlated with better OS (Figure 5(c)). Moreover, the Kaplan–Meier curve was drawn between the two clusters, and the results showed that clusters 3 and 1 and clusters 3 and 2 had significant differences (all , Figures 5(d) and 5(e)). There was no significant difference between clusters 1 and 2 (, Figure 5(f)).

(a)

(b)

(c)

(d)

(e)

(f)

3.5. Identification of DEGs between Subtypes

A total of 1846, 679, and 2257 DEGs were identified between clusters 1 and 2, clusters 1 and 3, and clusters 2 and 3, respectively, and 47 common DEGs were identified (Figure 6(a)). PCA of the 47 common DEGs in clusters 1, 2, and 3 is shown in Figure 6(b). Moreover, enrichment analysis was conducted on 47 common DEGs. These were enriched in 13 gene ontology (GO)-biology process (BP) (Figure 6(c)) and three KEGG pathways (osteoclast differentiation, amoebiasis, and rheumatoid arthritis) (Figure 6(d)).

(a)

(b)

(c)

(d)

Figure 6

Identification of differentially expressed genes (DEGs) between subtypes. (a) The heatmap of 47 common DEGs. (b) The principal component analysis (PCA) of 47 common DEGs. (c) The gene ontology (GO) enrichment analysis. Length of bar chart represents the number of genes enriched in each term. The color of the ball represents the value. (d) The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. The size of the ball represents the number of genes enriched in each term. The color of the ball represents the value.

3.6. Tumor-Infiltrating Immune Cells

The CIBERSORT algorithm was used to evaluate the proportion of subtypes of immune cells, and the differences in these proportions were analyzed using the Kruskal–Wallis test. The results showed significant differences in four subtypes of immune cells, including activated dendritic cells (), naïve CD4 T cells (), M1 macrophages (), and plasma cells (). A box diagram of the proportion of 22 subtypes of immune cells is shown in Figure 7.

3.7. TMB Analysis

The waterfall plot of the top 20 genes with the highest mutation frequency in each subtype is shown in Figure 8(a), which revealed that the gene mutation frequency in cluster 2 was low, the gene mutation frequency in cluster 1 was high, and TTN, MUC16, BRAF, and DNAH5 in the three subtypes had high mutation frequencies. In addition, the differences in TMB values in subtypes were analyzed, and the TMB value in cluster 2 was significantly lower than that in the other two subtypes () (Figure 8(b)).

(a)

(b)

(c)

(d)

(e)

(f)

3.8. Validation Analysis

The GSE18509 and GSE24996 datasets were used as validation datasets to verify that the 27 miRNAs were related to melanoma. As none of the probes in the GSE18509 dataset matched hsa-miR-424-3p and hsa-miR-23b-5p, 25 miRNAs were verified. The heatmap showed that 25 miRNAs were well distinguished between samples from the benign and metastatic groups (Figure 8(c)), suggesting that 25 miRNAs were related to melanoma. The differences in expression levels between benign and metastatic groups were evaluated using the t-test, and the box diagram showed that the expression levels of most miRNAs were significantly different, and the regulation trend was consistent with the results of the training dataset (Figure 8(d)). Moreover, the expression of 27 DEmiRNAs was analyzed using the GSE24996 dataset. Because the two probes in the GSE24996 dataset were not matched, 25 miRNAs were verified. The results showed that the expression of 22 DEmiRNAs was significantly different among the benign, primary, and metastatic groups (Figures 8(e) and 8(f)), suggesting that the above results were reliable.

4. Discussion

miRNAs are aberrantly expressed in numerous cancers [25], including melanoma. This study aimed to screen miRNA biomarkers of melanoma using data downloaded from public databases. hsa-miR-106b-5p, hsa-miR-27b-3p, and hsa-miR-141-3p had a higher degree and were regulated by numerous genes in the miRNA-gene network, and the target genes of miRNAs were involved in the MAPK and PI3K-Akt signaling pathways. Moreover, four miRNAs were associated with prognosis: hsa-let-7c-5p, hsa-miR-130b-3p, hsa-miR-142-3p, and hsa-miR-509-3p. Moreover, the bidirectional hierarchical clustering of the 27 miRNAs was classified into three subtypes, and TMB and four subtypes of immune cells, including activated dendritic cells, naïve CD4 T cells, M1 macrophages, and plasma cells, showed significant differences among the three subtypes.

Sixty-one common DEmiRNAs were screened between benign vs. primary, metastatic vs. benign, and metastatic vs. primary in the GSE34460 and GSE35579 datasets. Among these, 27 DEmiRNAs were related to disease progression. Moreover, the miRNA-gene network was constructed, and hsa-miR-106b-5p, hsa-miR-27b-3p, and hsa-miR-141-3p had a higher degree and were regulated by numerous genes. Chen et al. found that hsa-miR-106-5p promotes the cell cycle progression of malignant melanoma by targeting PTEN [26]. Bonnu et al. revealed that hsa-miR-106b-5p regulates prostate cancer by interfering with the endoplasmic reticulum stress repair pathway and reducing the expression of tumor suppressor genes involved in many biological processes [27]. Liu et al. suggested that hsa-miR-106b-5p promotes colorectal cancer cell migration and invasion [28]. Gencia et al. found that hsa-miR-27b-3p showed altered expression in all three melanoma types [29]. Miao et al. indicated that hsa-miR-27b-3p inhibits glioma development by targeting YAP1 [30]. Verrando et al. found that transnonachlor reduces the level of hsa-miR-141-3p in human melanocytes in vitro and promotes the characteristics of melanoma cells [31]. In addition, the KEGG pathway enrichment analysis showed that the target genes of 14 miRNAs were involved in the MAPK and PI3K-Akt signaling pathways. Chen et al. illustrated that FARP1 promotes cell proliferation by regulating the MAPK signaling pathway in cutaneous melanoma [32]. Wu et al. found that lncRNA OR3A4 facilitates the invasion and migration of melanoma cells through the PI3K/Akt signaling pathway [33]. Thus, miRNAs may be involved in the progression of melanoma via the MAPK and PI3K-Akt signaling pathways.

Three miRNAs, hsa-let-7c-5p, hsa-miR-130b-3p, and hsa-miR-142-3p, were related to both OS and DFF, and hsa-miR-509-3p was associated with OS. Liu et al. found that hsa-let-7c-5p was associated with melanoma metastasis [34]. Fu et al. suggested that hsa-let-7c-5p restrains breast cancer cell proliferation and facilitates apoptosis by targeting ERCC6 [35]. Murria et al. identified miRNAs associated with melanoma survival, and the results suggested poorer survival in melanomas with hsa-miR-130b-3p overexpression [36]. Tembe et al. indicated that hsa-miR-142-3p is related with the prognosis of patients with metastatic melanoma [37]. Li et al. found that GPC6 is an early biomarker for metastatic progression of melanoma, which can be regulated by hsa-miR-509-3p [38]. Patil et al. revealed that hsa-miR-509-3p inhibits the migration, invasion, and proliferation of osteosarcoma cells [39]. Our results were in line with those reported in previous studies, suggesting that hsa-let-7c-5p, hsa-miR-130b-3p, hsa-miR-142-3p, and hsa-miR-509-3p are associated with the prognosis of patients with melanoma.

Tumor-infiltrating immune cells play a vital role in promoting or inhibiting tumor growth [40]. The bidirectional hierarchical clustering of the 27 miRNAs was classified into three subtypes, and four subtypes of immune cells, including activated dendritic cells, naïve CD4 T cells, M1 macrophages, and plasma cells, showed significant differences among the three subtypes. Among the immune cells involved in cancer immunity, properly activated plasmacytoid dendritic cells play a vital role in building a bridge between inherent and adaptive immune responses and directly eliminating cancer cells [41]. Dendritic cells with potent antigen-presenting abilities have long been considered a critical factor in antitumor immunity [42]. Su et al. found that the abundance of naïve CD4⁺ T cells is closely related to Tregs, suggesting that the prognosis of breast cancer patients is poor [43]. Tumor-associated macrophages could be used as treatment targets in oncology [44, 45]. Plasma cells are the main effectors of adaptive immunity and are responsible for producing antibodies to protect the body against pathogens. Numerous studies have reported that plasma cells are related to the pathogenesis of cancer, such as triple-negative breast cancer [46], lung cancer [47], and colorectal cancer [48]. Moreover, TMB could be considered an immunotherapy biomarker for cancer [49–51]. In this study, the TMB value in cluster 2 was significantly lower than that in the other two subtypes, suggesting that these miRNAs might play important roles in the pathogenesis of melanoma by involving tumor-infiltrating immune cells and TMB.

However, this study has certain limitations. First, the data analyzed in this study were downloaded from public databases, and external validation is required. Moreover, the sample size should be expanded in future studies. Relevant experiments must be performed to verify the miRNAs identified in our bioinformatics analyses.

5. Conclusions

In summary, this study identified numerous novel miRNA biomarkers, including hsa-miR-106b-5p, hsa-miR-27b-3p, and hsa-miR-141-3p, which play important roles in melanoma and might provide new insights into the pathogenesis of melanoma.

Data Availability

The data supporting the findings of this study are available from the GEO database (GSE34460, GSE35579, GSE18509, and GSE24996 datasets).

Ethical Approval

The animal experiments in this study were approved by the Ethics Committee of the First Affiliated Hospital of Harbin Medical University.

Disclosure

Xuerui Wu and Yu Wang are the co-first authors.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Authors’ Contributions

XW, LC, and SZ conceptualized and designed the study. XW, CC, YX, and YW acquired data. LC, XW, and SZ analyzed and interpreted data, XW, CC, YX, and YW performed statistical analysis. LC obtained fund and revised the manuscript for important intellectual content. XW and SZ drafted the manuscript. All authors have read and approved the final manuscript. Xuerui Wu and Yu Wang contributed equally to this work.

Acknowledgments

This work was supported by the Research Innovation Fund of the First Affiliated Hospital of Harbin Medical University (2019M07).

References

G. Brancaccio, T. Russo, A. Lallas, E. Moscarella, M. Agozzino, and G. Argenziano, “Melanoma: clinical and dermoscopic diagnosis,” Giornale italiano di dermatologia e venereologia: Organo ufficiale Societa italiana di dermatologia e sifilografia, vol. 152, no. 3, pp. 213–223, 2017.
View at: Publisher Site | Google Scholar
C. Longo and G. Pellacani, “Melanomas,” Dermatologic Clinics, vol. 34, no. 4, pp. 411–419, 2016.
View at: Publisher Site | Google Scholar
K. Eddy and S. Chen, “Overcoming immune evasion in melanoma,” International Journal of Molecular Sciences, vol. 21, no. 23, p. 8984, 2020.
View at: Publisher Site | Google Scholar
K. Namikawa and N. Yamazaki, “Targeted therapy and immunotherapy for melanoma in Japan,” Current Treatment Options in Oncology, vol. 20, no. 1, p. 7, 2019.
View at: Google Scholar
S. S. Yde, P. Sjoegren, M. Heje, and L. B. Stolle, “Mucosal melanoma: a literature review,” Current Oncology Reports, vol. 20, no. 3, p. 28, 2018.
View at: Publisher Site | Google Scholar
Y. Wu, Y. Wang, L. Wang, P. Yin, Y. Lin, and M. Zhou, “Burden of melanoma in China, 1990-2017: findings from the 2017 global burden of disease study,” International Journal of Cancer, vol. 147, no. 3, pp. 692–701, 2020.
View at: Publisher Site | Google Scholar
J. Weber, M. Mandala, M. Del Vecchio et al., “Adjuvant nivolumab versus ipilimumab in resected stage III or IV melanoma,” New England Journal of Medicine, vol. 377, no. 19, pp. 1824–1835, 2017.
View at: Publisher Site | Google Scholar
T. X. Lu and M. E. Rothenberg, “MicroRNA,” The Journal of Allergy and Clinical Immunology, vol. 141, no. 4, pp. 1202–1207, 2018.
View at: Publisher Site | Google Scholar
I. Cannell, Y. Kong, and M. Bushell, “How do microRNAs regulate gene expression?” Biochemical Society Transactions, vol. 36, no. 6, pp. 1224–1231, 2008.
View at: Publisher Site | Google Scholar
T. Lu, S. Chen, L. Qu, Y. Wang, H. Chen, and C. He, “Identification of a five-miRNA signature predicting survival in cutaneous melanoma cancer patients,” PeerJ, vol. 7, Article ID e7831, 2019.
View at: Publisher Site | Google Scholar
J. Zhu, J. Deng, L. Zhang et al., “Reconstruction of lncRNA-miRNA-mRNA network based on competitive endogenous RNA reveals functional lncRNAs in skin cutaneous melanoma,” BMC Cancer, vol. 20, no. 1, p. 927, 2020.
View at: Publisher Site | Google Scholar
S. Huang, C. Zou, Y. Tang et al., “miR-582-3p and miR-582-5p suppress prostate cancer metastasis to bone by repressing TGF-β signaling,” Molecular Therapy—Nucleic Acids, vol. 16, pp. 91–104, 2019.
View at: Publisher Site | Google Scholar
T. Barrett, D. B. Troup, S. E. Wilhite et al., “NCBI GEO: mining tens of millions of expression profiles--database and tools update,” Nucleic Acids Research, vol. 35, pp. D760–D765, 2007.
View at: Publisher Site | Google Scholar
M. Goldman, B. Craft, A. N. Brooks, J. Zhu, and D. Haussler, “The UCSC xena platform for cancer genomics data visualization and interpretation,” 2018, https://www.researchgate.net/publication/325249266_The_UCSC_Xena_Platform_for_cancer_genomics_data_visualization_and_interpretation.
View at: Google Scholar
G. K. Smyth, “Limma: linear models for microarray data,” Bioinformatics and computational biology solutions using R and Bioconductor, Springer, New York, NY, USA, 2013.
View at: Google Scholar
K. Ana, M. Birgaoanu, and S. Griffiths-Jones, “miRBase: from microRNA sequences to function,” Nucleic acids research, vol. 47, 2019.
View at: Google Scholar
H. Dweep and N. Gretz, “miRWalk2.0: a comprehensive atlas of microRNA-target interactions,” Nature Methods, vol. 12, no. 8, p. 697, 2015.
View at: Publisher Site | Google Scholar
G. Stelzer, N. Rosen, I. Plaschkes et al., “The GeneCards suite: from gene data mining to disease Genome sequence analyses,” Current Protocols in Bioinformatics, vol. 54, no. 1, 2016.
View at: Publisher Site | Google Scholar
P. Shannon, A. Markiel, O. Ozier et al., “Cytoscape: a software environment for integrated models of biomolecular interaction networks,” Genome Research, vol. 13, no. 11, pp. 2498–2504, 2003.
View at: Publisher Site | Google Scholar
G. Yu, L. G. Wang, Y. Han, and Q. Y. He, “clusterProfiler: an R package for comparing biological themes among gene clusters,” Omics-a Journal of Integrative Biology, vol. 16, no. 5, pp. 284–287, 2012.
View at: Publisher Site | Google Scholar
A. Johnson, “(R) pheatmap—clustering on some columns only,” 2016, https://stat.ethz.ch/pipermail/r-help/2016-October/442514.html.
View at: Google Scholar
D. W. Huang, B. T. Sherman, and R. A. Lempicki, “Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources,” Nature Protocols, vol. 4, no. 1, pp. 44–57, 2009.
View at: Publisher Site | Google Scholar
P. Charoentong, F. Finotello, M. Angelova et al., “Pan-cancer immunogenomic analyses reveal genotype-immunophenotype relationships and predictors of response to checkpoint blockade,” Cell Reports, vol. 18, no. 1, pp. 248–262, 2017.
View at: Publisher Site | Google Scholar
A. Mayakonda, D. C. Lin, Y. Assenov, C. Plass, and H. P. Koeffler, “Maftools: efficient and comprehensive analysis of somatic variants in cancer,” Genome Research, vol. 28, no. 11, pp. 1747–1756, 2018.
View at: Publisher Site | Google Scholar
Y. S. Lee and A. Dutta, “MicroRNAs in cancer,” Annual Review of Pathology: Mechanisms of Disease, vol. 4, no. 1, pp. 199–227, 2009.
View at: Publisher Site | Google Scholar
X. Chen, P. Chen, S. S. Chen, T. Ma, G. Shi, and Y. Zhou, “miR-106b-5p promotes cell cycle progression of malignant melanoma by targeting PTEN,” Oncology Reports, vol. 39, no. 1, pp. 331–337, 2018.
View at: Google Scholar
C. H. Bonnu, A. N. Ramadhani, R. B. Saputro et al., “The potential of hsa-mir-106b-5p as liquid biomarker in prostate cancer patients in Indonesia,” Asian Pacific Journal of Cancer Prevention, vol. 22, no. 3, pp. 837–842, 2021.
View at: Publisher Site | Google Scholar
S. Liu, G. An, Q. Cao, T. Li, X. Jia, and L. Lei, “The miR-106b/NR2F2-AS1/PLEKHO2 Axis regulates migration and invasion of colorectal cancer through the MAPK pathway,” International Journal of Molecular Sciences, vol. 22, no. 11, p. 5877, 2012.
View at: Publisher Site | Google Scholar
I. Gencia, F. Baderca, S. Avram et al., “A preliminary study of microRNA expression in different types of primary melanoma,” Bosnian Journal of Basic Medical Sciences, vol. 20, no. 2, pp. 197–208, 2020.
View at: Publisher Site | Google Scholar
W. Miao, N. Li, B. Gu, G. Yi, Z. Su, and H. Cheng, “MiR-27b-3p suppresses glioma development via targeting YAP1,” Biochemistry and Cell Biology, vol. 98, no. 4, pp. 466–473, 2020.
View at: Publisher Site | Google Scholar
P. Verrando, M. Capovilla, and R. Rahmani, “Trans-nonachlor decreases miR-141-3p levels in human melanocytes in vitro promoting melanoma cell characteristics and shows a multigenerational impact on miR-8 levels in Drosophila,” Toxicology, vol. 368-369, pp. 129–141, 2016.
View at: Publisher Site | Google Scholar
Z. H. Chen and L. H. Wang, “FARP1 facilitates cell proliferation through modulating MAPK signaling pathway in cutaneous melanoma,” The American Journal of Dermatopathology, vol. 41, no. 12, pp. 908–913, 2019.
View at: Publisher Site | Google Scholar
J. Wu, M. Y. Zhou, X. P. Yu, Y. Wu, and P. L. Xie, “Long noncoding RNA OR3A4 promotes the migration and invasion of melanoma through the PI3K/AKT signaling pathway,” European Review for Medical and Pharmacological Sciences, vol. 24, no. 21, Article ID 10917, 2020.
View at: Publisher Site | Google Scholar
X. Y. Liu, Q. S. Tang, H. C. Chen, X. L. Jiang, and H. Fang, “Lentiviral miR30-based RNA interference against heparanase suppresses melanoma metastasis with lower liver and lung toxicity,” International Journal of Biological Sciences, vol. 9, no. 6, pp. 564–577, 2013.
View at: Publisher Site | Google Scholar
X. Fu, X. Mao, Y. Wang, X. Ding, and Y. Li, “Let-7c-5p inhibits cell proliferation and induces cell apoptosis by targeting ERCC6 in breast cancer,” Oncology Reports, vol. 38, no. 3, pp. 1851–1856, 2017.
View at: Publisher Site | Google Scholar
R. Murria Estal, B. de Unamuno Bustos, G. Pérez Simó et al., “MicroRNAs expression associated with aggressive clinicopathological features and poor prognosis in primary cutaneous melanomas,” Melanoma Research, vol. 31, no. 1, pp. 18–26, 2021.
View at: Publisher Site | Google Scholar
V. Tembe, S. J. Schramm, M. S. Stark et al., “MicroRNA and mRNA expression profiling in metastatic melanoma reveal associations with BRAF mutation and patient prognosis,” Pigment cell & melanoma research, vol. 28, no. 3, pp. 254–266, 2015.
View at: Publisher Site | Google Scholar
K. A. Lisowska, M. Pindel, K. Pietruczuk et al., “The influence of a single hemodialysis procedure on human T lymphocytes,” Scientific Reports, vol. 9, no. 1, p. 5041, 2019.
View at: Publisher Site | Google Scholar
S. L. Patil, A. Palat, Y. Pan et al., “MicroRNA-509-3p inhibits cellular migration, invasion, and proliferation, and sensitizes osteosarcoma to cisplatin,” Scientific Reports, vol. 9, no. 1, Article ID 19089, 2019.
View at: Publisher Site | Google Scholar
S. S. Wang, W. Liu, D. Ly, H. Xu, L. Qu, and L. Zhang, “Tumor-infiltrating B cells: their role and application in anti-tumor immunity in lung cancer,” Cellular and Molecular Immunology, vol. 16, no. 1, pp. 6–18, 2019.
View at: Publisher Site | Google Scholar
M. Monti, F. Consoli, R. Vescovi, M. Bugatti, and W. Vermi, “Human plasmacytoid dendritic cells and cutaneous melanoma,” Cells, vol. 9, no. 2, p. 417, 2020.
View at: Publisher Site | Google Scholar
F. Veglia and D. I. Gabrilovich, “Dendritic cells in cancer: the role revisited,” Current Opinion in Immunology, vol. 45, pp. 43–51, 2017.
View at: Publisher Site | Google Scholar
S. Su, J. Liao, J. Liu et al., “Blocking the recruitment of naive CD4 (+) T cells reverses immunosuppression in breast cancer,” Cell Research, vol. 27, no. 4, pp. 461–482, 2017.
View at: Publisher Site | Google Scholar
A. Mantovani, F. Marchesi, A. Malesci, L. Laghi, and P. Allavena, “Tumour-associated macrophages as treatment targets in oncology,” Nature Reviews Clinical Oncology, vol. 14, no. 7, pp. 399–416, 2017.
View at: Publisher Site | Google Scholar
Y. Xia, L. Rao, H. Yao, Z. Wang, P. Ning, and X. Chen, “Engineering macrophages for cancer immunotherapy and drug delivery,” Advanced Materials, vol. 32, no. 40, Article ID e2002054, 2020.
View at: Publisher Site | Google Scholar
I. Bar, I. Theate, S. Haussy et al., “MiR-210 is overexpressed in tumor-infiltrating plasma cells in triple-negative breast cancer,” Journal of Histochemistry and Cytochemistry, vol. 68, no. 1, pp. 25–32, 2020.
View at: Publisher Site | Google Scholar
M. Fujimoto, A. Yoshizawa, S. Sumiyoshi et al., “Stromal plasma cells expressing immunoglobulin G4 subclass in non-small cell lung cancer,” Human Pathology, vol. 44, no. 8, pp. 1569–1576, 2013.
View at: Publisher Site | Google Scholar
J. Berntsson, B. Nodin, J. Eberhard, P. Micke, and K. Jirström, “Prognostic impact of tumour-infiltrating B cells and plasma cells in colorectal cancer,” International Journal of Cancer, vol. 139, no. 5, pp. 1129–1139, 2016.
View at: Publisher Site | Google Scholar
T. A. Chan, M. Yarchoan, E. Jaffee et al., “Development of tumor mutation burden as an immunotherapy biomarker: utility for the oncology clinic,” Annals of Oncology, vol. 30, no. 1, pp. 44–56, 2019.
View at: Publisher Site | Google Scholar
B. S. Jang, W. Han, and I. A. Kim, “Tumor mutation burden, immune checkpoint crosstalk and radiosensitivity in single-cell RNA sequencing data of breast cancer,” Radiotherapy & Oncology, vol. 142, pp. 202–209, 2020.
View at: Publisher Site | Google Scholar
R. Barroso-Sousa, E. Jain, O. Cohen et al., “Prevalence and mutational determinants of high tumor mutation burden in breast cancer,” Annals of Oncology, vol. 31, no. 3, pp. 387–394, 2020.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2022 Xuerui Wu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

303

Downloads

543

Citations