Next Article in Journal
Deep Neural Network Models for Colon Cancer Screening
Next Article in Special Issue
Increased MARCKS Activity in BRAF Inhibitor-Resistant Melanoma Cells Is Essential for Their Enhanced Metastatic Behavior Independent of Elevated WNT5A and IL-6 Signaling
Previous Article in Journal
Adding Concomitant Chemotherapy to Postoperative Radiotherapy in Oral Cavity Carcinoma with Minor Risk Factors: Systematic Review of the Literature and Meta-Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Circulating miR-1246 and miR-485-3p as Promising Biomarkers of Clinical Response and Outcome in Melanoma Patients Treated with Targeted Therapy

by
Lauretta Levati
1,
Cristian Bassi
2,3,
Simona Mastroeni
4,
Laura Lupini
2,
Gian Carlo Antonini Cappellini
5,†,
Laura Bonmassar
1,
Ester Alvino
6,
Simona Caporali
1,‡,
Pedro Miguel Lacal
1,
Maria Grazia Narducci
1,
Ivan Molineris
7,
Federica De Galitiis
5,
Massimo Negrini
2,3,
Giandomenico Russo
1 and
Stefania D’Atri
1,*
1
Laboratory of Molecular Oncology, IDI-IRCCS, Via dei Monti di Creta 104, 00167 Rome, Italy
2
Department of Translational Medicine, University of Ferrara, Via Fossato di Mortara 70, 44121 Ferrara, Italy
3
LTTA Center, University of Ferrara, Via Fossato di Mortara 70, 44121 Ferrara, Italy
4
Clinical Epidemiology Unit, IDI-IRCCS, Via dei Monti di Creta 104, 00167 Rome, Italy
5
Department of Oncology and Dermatological Oncology, IDI-IRCCS, Via dei Monti di Creta 104, 00167 Rome, Italy
6
Institute of Translational Pharmacology, National Council of Research, Via Fosso del Cavaliere 100, 00133 Rome, Italy
7
Department of Life Science and System Biology, University of Turin, Via Accademia Albertina 13, 10123 Turin, Italy
*
Author to whom correspondence should be addressed.
Current address: UOC Oncologia, Interpresidio ASL RM2, Via dei Monti Tiburtini 387, 00157 Rome, Italy.
Current address: Regional Transplant Center Lazio (CRTL), San Camillo Hospital, Circonvallazione Gianicolense 87, 00152 Rome, Italy.
Cancers 2022, 14(15), 3706; https://doi.org/10.3390/cancers14153706
Submission received: 26 May 2022 / Revised: 16 July 2022 / Accepted: 24 July 2022 / Published: 29 July 2022

Abstract

:

Simple Summary

Several therapeutic options exist for patients with advanced BRAF-mutant melanoma. Biomarkers able to identify patients with refractory disease or with poor progression-free survival (PFS) and overall survival (OS) expectancy with respect to specific treatments might allow a more personalized therapeutic approach. Here, we profiled plasma miRNAs at baseline and at progression in patients treated with BRAF inhibitors (BRAFi) or BRAFi + MEKi. Selected miRNAs associated with response to therapy were subjected to validation by real-time quantitative RT-PCR. Receiver Operating Characteristics (ROC), Kaplan–Meier and univariate and multivariate Cox regression analyses were performed on miR-1246 and miR-485-3p baseline levels to assess their ability to discriminate between responding and non-responding patients and to determine their prognostic value for PFS and OS. Globally, our results suggest that circulating miR-1246 and miR-485-3p could be valuable biomarkers for identifying patients most likely to be resistant to targeted therapy or with a poor expectancy of PFS and OS. Prospective studies in a larger cohort of patients are warranted.

Abstract

Despite the significant improvements in advanced melanoma therapy, there is still a pressing need for biomarkers that can predict patient response and prognosis, and therefore support rational treatment decisions. Here, we investigated whether circulating miRNAs could be biomarkers of clinical outcomes in patients treated with targeted therapy. Using next-generation sequencing, we profiled plasma miRNAs at baseline and at progression in patients treated with BRAF inhibitors (BRAFi) or BRAFi + MEKi. Selected miRNAs associated with response to therapy were subjected to validation by real-time quantitative RT-PCR. Receiver Operating Characteristics (ROC), Kaplan–Meier and univariate and multivariate Cox regression analyses were performed on the validated miR-1246 and miR-485-3p baseline levels. The median baseline levels of miR-1246 and miR-485-3p were significantly higher and lower, respectively, in the group of patients not responding to therapy (NRs) as compared with the group of responding patients (Rs). In Rs, a trend toward an increase in miR-1246 and a decrease in miR-485-3p was observed at progression. Baseline miR-1246 level and the miR-1246/miR-485-3p ratio showed a good ability to discriminate between Rs and NRs. Poorer PFS and OS were observed in patients with unfavorable levels of at least one miRNA. In multivariate analysis, a low level of miR-485-3p and a high miR-1246/miR-485-3p ratio remained independent negative prognostic factors for PFS, while a high miR-1246/miR-485-3p ratio was associated with an increased risk of mortality, although statistical significance was not reached. Evaluation of miR-1246 and miR-485-3p baseline plasma levels might help clinicians to identify melanoma patients most likely to be unresponsive to targeted therapy or at higher risk for short-term PFS and mortality, thus improving their management.

1. Introduction

Cutaneous melanoma is an aggressive cancer whose incidence is increasing constantly all over the world. Although surgical resection is an effective treatment for early-stage melanoma, natural prognosis of patients with stage IV disease is extremely poor, with a median survival rate of less than 1 year and a 5-year survival rate of less than 10% [1].
In recent years, a remarkable advance in the therapy of metastatic melanoma has been achieved with the development of BRAF and MEK inhibitors (BRAFi, MEKi), approved for BRAF-mutant tumors, and monoclonal antibodies targeting T-lymphocyte-associated antigen 4 and programmed cell-death protein 1, approved for both BRAF-mutant and BRAF wild-type melanomas.
Monotherapy with BRAFi and MEKi yields objective response rates (ORR) of 50–60% and of 20–30%, respectively, and significantly prolongs progression-free survival (PFS) and overall survival (OS) with respect to dacarbazine chemotherapy [2,3]. Better results are achieved with the combination of BRAFi + MEKi—nowadays the standard-of-care of BRAF targeted approach—with ORR of 65–75% and PFS and OS rates of 15–22% and 30–38%, respectively, at 5 years [4].
Although BRAFi and the combination of BRAFi + MEKi induce undoubtedly high ORR, a non-negligible percentage of patients shows primary resistance. Moreover, long-term efficacy of therapy is limited by the development of drug resistance in nearly all patients [5,6]. Novel therapeutic approaches able to improve response to BRAFi and MEKi and to mitigate or overcome acquired resistance, as well as biomarkers able to better predict patients’ response to therapy and survival outcomes, are therefore urgently needed.
MicroRNAs (miRNAs) are a class of small non-coding RNAs that negatively regulate gene expression at the post-transcriptional level by base paring to—usually—the 3′untraslated region (UTR) of target mRNAs [7]. They are involved in numerous cellular processes, including apoptosis, proliferation, differentiation, and metabolism [7]. Each miRNA can target numerous transcripts, whereas different miRNAs can converge on a single mRNA [7]. Therefore, miRNAs can regulate a large fraction of protein-coding genes.
Aberrant expression of miRNAs has been demonstrated in a variety of human cancers, where miRNAs can operate as oncogenes or tumor suppressor genes [8,9]. In the clinical setting, tumor miRNA signatures have been correlated with patients’ outcome and response to therapy in several malignancies [10,11,12]. miRNAs appear, therefore, to have both diagnostic and prognostic significance and to potentially constitute novel targets and therapeutic agents for cancer treatment [12,13,14]. Tumor-derived miRNAs can also be detected in plasma/serum of cancer patients—associated with proteins, lipoproteins or included in extracellular vesicles (EVs)—and accumulating experimental evidence points out that circulating cell-free miRNAs (cf-miRNAs) can provide a non-invasive strategy for cancer diagnosis and prognosis, and for predicting and/or monitoring patients’ responses to therapy [15,16].
In melanoma, aberrant expression of miRNAs has been widely documented, and shown to affect several well-known pathways involved in the regulation of cell proliferation, invasiveness and survival, as well as to contribute to drug resistance [17,18,19,20,21,22,23,24,25,26,27,28,29]. Moreover, the potential value of single cf-miRNAs or cf-miRNA signatures as biomarkers for melanoma diagnosis, staging, risk of recurrence and survival prognosis has been highlighted by numerous studies [26,28,30,31]. On the other hand, to the best of our knowledge, the role of cf-miRNA as biomarkers of patients’ responses to targeted therapy has been addressed by a very limited number of investigations [21,32].
In the present study, we used next-generation sequencing (NGS) for profiling plasma cf-miRNAs at baseline and at progression in melanoma patients treated with BRAFi monotherapy or the combination of BRAFi + MEKi to identify cf-miRNAs associated with response to therapy. Selected cf-miRNAs were then subjected to validation by real-time quantitative RT-PCR and investigated for their potential value as biomarkers predictive of patients’ clinical outcomes.

2. Materials and Methods

2.1. Study Population

A total of 57 patients with inoperable stage IIIc or stage IV melanoma consecutively enrolled at IDI-IRCCS from 2013 to 2018, and for whom peripheral blood samples had been collected before the start of therapy, and when possible, at disease progression, were included in the study.
Patients who started therapy before October 2013 were either treated with vemurafenib (n. 4), enrolled in the therapeutic protocol for compassionate use of dabrafenib monotherapy (n. 12), or treated with dabrafenib + trametinib within the clinical trial MEK116513 (EUDRACT NUMBER 2011-006088-23) (n. 1). On October 2013, the compassionate use of dabrafenib + trametinib became possible, and 9 patients under dabrafenib monotherapy could switch to the combination. From October 2013 to June 2018, 29 additional patients were enrolled in the therapeutic protocol for compassionate use of dabrafenib + trametinib, which ensured dabrafenib and/or trametinib supply to patients until drug official approval by AIFA (Agenzia Italiana del Farmaco). Four patients starting therapy in 2014 and not eligible for compassionate use of dabrafenib + trametinib were treated with vemurafenib (n. 3) or dabrafenib (n.1) alone. Seven patients started therapy with vemurafenib + cobimetinib between 2016 and 2017. Cobimetinib was supplied to those patients for compassionate use until its approval by AIFA. Demographic and clinical–pathological features of patients are illustrated in Table 1 and Table S1.
Dabrafenib was given at the dose of 150 mg BID, vemurafenib at the dose of 960 mg BID, dabrafenib + trametinib at the dose of 150 mg BID and 2 mg/die, respectively, and vemurafenib + cobimetinib at the dose of 960 mg BID and 60 mg/die, respectively, for three weeks with one week of break. Baseline evaluation included medical history, physical examination, assessment of biochemical parameters—including serum lactate dehydrogenase (sLDH)—and radiologic tumor assessment with computer tomography (CT) or positron emission tomography scans. All patients underwent physical examination and evaluation of biochemical parameters monthly, whereas tumor response was determined with CT every three months, or less if required. Tumor response was classified as complete response (CR), partial response (PR), stable disease (SD) or progressive disease (PD) according to RECIST 1.1 criteria [33]. Patients achieving CR or PR were classified as “Responders” (Rs), whereas Pts with SD or PD as best response were classified as “Non-responders” (NRs). PFS was defined as the time from the start of therapy to the first observation of disease progression per RECIST 1.1 or death (event) or last follow-up (censored). OS was calculated from the start of therapy to death (event) or last follow-up (censored).
The study was conducted in accordance with Good Clinical Practice Guidelines and the Declaration of Helsinki. The study was also approved by the IDI-IRCCS Ethics Committee (ID #407/1, 2013) and a written informed consent was obtained from all patients.

2.2. Plasma Preparation and RNA Extraction

Blood was collected into BD Vacutainer® tubes (#367704, BD Biosciences, Plymouth, UK) and double centrifuged at 1200× g for 10 min at 4 °C, and the plasma was aliquoted and stored at −80 °C until use.
Total RNA, including miRNAs, was isolated from 400 μL of plasma samples using the Maxwell® RSC miRNA Tissue Kit AS1460 (Promega, Madison, WI, USA) in combination with the Maxwell® RSC Instrument AS4500 (Promega). Plasma volume was mixed with 200 μL chilled 1-thioglycerol/homogenization solution and 200 μL lysis solution. After adding 25 μL of Proteinase K, samples were incubated for 10 min at room temperature. To control for variation in recovery and amplification efficiency between RNA preparations, cel-miR-39 (14 μL of 1.6 × 106 copies/μL dilution) was added as spike-in exogenous control before loading samples in the Maxwell RSC Cartridge. All samples were eluted in 50 μL nuclease-free water, and extracted RNA was aliquoted and stored at −80 °C until use.

2.3. Small RNA Library Preparation and Sequencing

For small RNA sequencing (small RNA-seq), 5 μL of extracted RNA was used as starting material to generate cDNA barcoded libraries. Library preparation was performed using the Illumina TruSeq Small RNA Library Prep Kit v2 (Illumina, San Diego, CA, USA) following the manufacturer’s instructions. Briefly, specific adapters were ligated to the 5′ and 3′ ends of RNA molecules and used as templates for reverse transcription. The synthesized cDNA fragments were subsequently amplified by 15 cycles of polymerase chain reaction using unique index primers. The obtained libraries were quantified on Bioanalyzer 2100 with High Sensitivity DNA Kit (Agilent Technologies, Inc., Santa Clara, CA, USA), pooled together and run on 6% Novex TBE-PAGE gel (Thermo Fisher Scientific, Waltham, MA, USA) for size selection. cDNA fragments of 140–160 bp in size were excised from the gel and purified. A concluding Bioanalyzer 2100 run with a High Sensitivity DNA Kit (Agilent) was performed to verify the length range of the library and calculate the final working solution concentration. Small RNA-seq was carried out according to Illumina pipeline on NextSeq 500 Instrument (Illumina) with the Next Seq High Output kit v2 (75 cycles) (Illumina).

2.4. miRNA Data Analysis

Raw base-call data were demultiplexed using Illumina BaseSpace Sequence Hub and converted to FASTQ format. After a quality check with FastQC v.0.11.5 (https://www.bioinformatics.babraham.ac.uk/projects/fastqc/) (accessed on 3 August 2020, the adapter sequences were trimmed using Cutadapt v.1.6 (http://cutadapt.readthedocs.io/en/stable/index.html) (accessed on 3 August 2020) with these settings: -f fastq -O 5 -q 20 -m 16 -M 30 -a TGGAATTCTC, which also removed sequences <16 nucleotides. Reads were mapped using the STAR v.2.5.0a algorithm (https://www.ncbi.nlm.nih.gov/pubmed/23104886) (accessed on 3 August 2020) with these parameters: --outFilterMismatchNoverLmaverLmax 0.1 --outFilterMatchNmin 16 --outFilterScoreMinOverLread 0 --outFilterMatchNminOverLread 0 --alignIntronMax 1. The reference genome consisted of human miRNA sequences from the miRbase 21 database.
Raw counts from mapped reads were obtained using the htseq-count script v. 07.2 (--mode intersection-strict) from the HTSeq tools (http://www-huber.embl.de/HTSeq/doc/overview.html) (accessed on 3 August 2020). Counts were normalized using DESeq2 v.1.28.1 bioconductor package (http://bioconductor.org/packages/release/bioc/html/DESeq2.html) (accessed on 3 August 2020).
Normalized sequencing data were imported and analyzed in Genespring GX software v.14.9.1 (Agilent Technologies). cf-miRNAs with low expression levels in at least 50% of samples (expression value between 0–3 in the normalized data) were filtered out. A list of 290 cf-miRNAs remained for the subsequent analyses (Table S2). Differentially expressed (DE) cf-miRNAs were identified using a fold change ≥ 1.5 filter and a Benjamini-Hochberg [34] adjusted p-value cut-off of 0.15. An additional quality control on sequencing data was performed with miRTrace v. 1.0.1 (https://pubmed.ncbi.nlm.nih.gov/30514392/) (accessed on 22 June 2022).

2.5. Real-Time Quantitative RT-PCR (qRT-PCR)

The expression of mature hsa-miR-1246 and hsa-miR-485-3p was assayed using the TaqMan® MicroRNA Reverse Transcription Kit, the TaqMan® Universal PCR Master Mix No AmpErase® UNG and the specific TaqMan® MicroRNA Assays, all purchased from Applied Biosystems (Foster City, CA, USA). All experimental procedures were performed according to the manufacturer’s protocols. Five μL of total RNA extracted from plasma samples as described for cf-miRNA profiling, were reverse transcribed in a final volume of 15 μL, and qRT-PCR was done on an ABI PRISM 7000 Sequence Detection System (Applied Biosystems) in a final volume of 20 μL. Since no specific TaqMan® MicroRNA Assay was available for has-miR-92b-3p, the expression levels of this miRNA were investigated using the Qiagen miScript PCR System, based on Sybr Green chemistry. Total RNA (5 μL) was used for cDNA synthesis and reverse transcription was carried out using the miScript II RT kit (Qiagen, Valencia, CA, USA) following the manufacturer’s protocol with HiSpec buffer. qRT-PCR was performed in a final volume of 25 μL using miScript II SYBR green PCR kit (Qiagen) and the pre-validated miScript Primer Assays for hsa-miR-92b-3p (cat.# MS00032144) and for the spike-in control cel-miR-39 (cat.# MS00019789) (Qiagen) according to manufacturer’s instruction. The cycling program consisted of 15 min at 95 °C, followed by 40 cycles of denaturation at 94 °C for 15 s, annealing at 55 °C for 30 s, and extension at 72 °C for 30 s.
All qRT-PCR reactions were run in duplicate. The expression of cf-miRNAs relative to cel-miR-39 was determined using the formula 2−∆Ct, where ∆CT = CTmiR − CTcel-miR-39, and CT (i.e., threshold cycle) indicates the fractional cycle number at which the amount of amplified target reaches a fixed threshold.

2.6. Statistical Analysis

Receiver Operating Characteristics (ROC) analysis [35] was performed to assess the ability of specific cf-miRNAs to discriminate between Rs and NRs. Based on the highest Youden’s Index [36], the optimal cut-off value for sensitivity and specificity was estimated. ROC analysis was performed on the 29 cf-miRNAs DE at T0 between Rs and NRs according to small RNA-seq data, on miR-1246 and miR-485-3p T0 plasma levels determined by qRT-PCR and on the resulting miR-1246 to miR-485-3p (miR-1246/miR485-3p) ratio.
Only for cf-miRNAs DE at T0 between Rs and NRs according to small RNA-seq data, PFS was categorized into 3 groups (≤180 days; 181–365 days; ≥366 days). Differences in levels of cf-miRNAs between PFS groups were analyzed with the Cuzick’s test for trend across categories [37].
Plasma levels of cf-miRNAs determined by qRT-PCR were presented as median and Interquartile Range (IQR). The Mann–Whitney U test was used to compare between-group differences (i.e., T0/Rs and T0/NRs), while the Wilcoxon matched-pairs signed-rank test was used to evaluate before–after differences (i.e., TP/Rs vs. T0/Rs).
PFS and OS curves were estimated with the Kaplan–Meier method. The log-rank test was used to compare the PFS and OS curves in patient groups defined on the basis of miR-1246, miR-485-3p and miR-1246/miR-485-3p ratio categories according to AUC optimal cut-off, as well as of miR-1246 and miR-485-3p combinations.
The Cox proportional hazard model was used in the univariate and multivariate analyses to estimate crude and adjusted hazard ratios (HR) with 95% confidence intervals (95% CIs).
Statistical analyses were performed using Stata software version 15 (StataCorp LLC, College Station, TX, USA) and GraphPad prism software version 5.04 (La Jolla, CA, USA).

3. Results

3.1. Identification of cf-miRNAs Associated with Melanoma Resistance to Targeted Therapy Based on Small RNA-seq Data

Profiling of plasma cf-miRNAs by small RNA-seq was performed in a cohort of 33 patients treated with either vemurafenib or dabrafenib monotherapy (hereafter referred to as BRAFi/MONO) (n. 11), with dabrafenib for 4–8 months and then with dabrafenib + trametinib (n. 7), or with the combination of dabrafenib + trametinib or vemurafenib + cobimetinib (hereafter referred to as COMBO) (n. 15). Among the 33 patients, 28 had achieved an objective clinical response (Rs) and 5 had not (NRs) (Table S1). To identify miRNAs potentially associated with primary resistance to therapy, cf-miRNA profiling was carried out on plasma samples collected before the start of therapy (T0) from Rs and NRs, whereas to identify miRNAs potentially associated with secondary resistance, small RNA-seq was carried out on plasma samples collected from Rs at disease progression on therapy (TP). Two Rs were still in CR at the time of cf-miRNA assessment, and no TP samples were, therefore, available. Moreover, in the group of Rs, the T0 sample was not available for two patients, whereas the TP sample was not available for one patient. Overall, 31 T0 and 25 TP samples were, therefore, subjected to small RNA-seq.
Sequencing yielded 50.75 gigabases (Gb) of high-quality data with an average of 11.6 × 106 single-end reads per sample. The quality control performed on sequencing data by miRTrace (Figure S1) showed that 95.32% of the sequenced nucleotides had a Phred score greater than 30 and that 20.9% of the reads had a length compatible with the length of miRNAs (20–25 nucleotides). Among the total reads, 37.73% had the Illumina adapter removed and were long enough to map reliably (18 nucleotides or longer). In detail, 45.9% of the reads with length >18 nucleotides were identified as miRNA reads, whereas 1.1% were identified as rRNA reads, 0.6% as tRNA reads, 5.5% as Illumina sequences, and 47% as unknown. Furthermore, no contamination from miRNAs of non-human species were reported (Figure S1). Around 30% of the total reads were mapped by the STAR algorithm to the reference.
Bioinformatic and statistical analyses were initially conducted to identify cf-miRNAs DE at T0 between the group of NRs (T0/NRs) and the group of Rs (T0/Rs), as well as between TP/Rs and T0/Rs, considering all patients independently of the therapy received. The analyses identified five cf-miRNAs upregulated and twenty-four downregulated in T0/NRs vs. T0/Rs (Table 2), and four cf-miRNAs upregulated and nine downregulated in TP/Rs vs. T0/Rs (Table 3). miR-369-3p, miR-485-3p, miR-487a-3p and miR-4286 were downregulated, while miR-1246 was upregulated in both T0/NRs vs. T0/R and TP/Rs vs. T0/Rs, suggesting a strong association with drug resistance.
Taking into consideration that all NRs had been subjected to BRAFi/MONO, whereas the Rs had been treated with either BRAFi/MONO or COMBO, or initially with dabrafenib and then with dabrafenib + trametinib, bioinformatic and statistical analyses were also conducted to identify cf-miRNAs DE between T0/NRs and T0/Rs, including in the group of Rs only those patients treated with BRAFi/MONO (n = 5) and those treated with dabrafenib followed by the addition of trametinib (n = 7), since these latter patients had achieved their best response on treatment with the BRAFi alone. In this case, only three cf-miRNAs were found to be DE in T0/NRs vs. T0/Rs, namely miR-1246 and miR-92b-3p, upregulated, and miR-6837-3p, downregulated (Table S3). All three cf-miRNAs were also present in the list of those DE in T0/NRs vs. T0/Rs considering all therapies (Table 2). The subsequent analyses on T0 samples were, therefore, performed without distinguishing patients according to the treatment received.

3.2. Identification of cf-miRNAs Able to Discriminate NRs from Rs or Associated with PFS Based on Small RNA-seq Data

To select among cf-miRNAs DE in T0/NRs vs. T0/Rs those with the highest potential as predictive biomarkers of patients’ response to therapy, we first performed a ROC analysis to assess the ability of cf-miRNA to discriminate between NRs and Rs. As illustrated in Table 4, four cf-miRNAs upregulated and six downregulated showed an AUC value > 0.80, which is recommended for a predictive marker with potential clinical utility [38].
We next evaluated the association of cf-miRNAs with patient PFS. To this end, differences in cf-miRNA expression levels between three PFS categories, namely ≤180 days, 181–365 days and ≥366 days, were considered. Eight cf-miRNAs showed a statistically significant trend of expression across the three PFS categories (Table 5).

3.3. Evaluation of T0 and TP Plasma Levels of miR-1246, miR-485-3p and miR-92b-3p by qRT-PCR

Based on small RNA-seq data analysis, we selected miR-1246, miR-92b-3p and miR-485-3p for further studies. Indeed, miR-1246 displayed the best performance in discriminating between NRs and Rs, was upregulated in T0/NRs vs. T0/Rs and in TP/Rs vs. T0/Rs, suggesting a strong association with resistance to therapy. Moreover, miR-1246 was found to be upregulated in melanoma specimens as compared to normal tissue [39] as well as in melanoma cell lines resistant to the BRAFi PLX4720 [40]. Finally, higher levels of miR-1246 were detected in serum EVs isolated from melanoma patients as compared with those obtained from control subjects [41]. miR-92b-3p showed a potential value as biomarkers predictive for primary resistance, being upregulated in T0/NRs vs. T0/Rs and showing a good ability to discriminate between NRs and Rs. In addition, its expression at T0 was significantly associated with PFS. Although no data are available about miR-92b-3p dysregulation in cutaneous melanoma, this miRNA has been reported to function as an oncomiR in different types of malignancies and to promote chemoresistance [42,43,44,45,46]. miR-485-3p was downregulated in both T0/NRs vs. T0/Rs and TP/Rs vs. T0/Rs. Furthermore, its T0 levels were significantly associated with PFS, and a tumor suppressor function has been recently suggested in melanoma [47].
T0 and TP plasma levels of miR-1246, miR-485-3p and miR-92b-3p were determined by qRT-PCR in the first cohort of patients and in an additional 24 patients comprising 5 NRs and 19 Rs treated with either dabrafenib followed by the addition of trametinib or with COMBO. (Table S1). Among the 19 new Rs, seven patients were still in response at the time of cf-miRNA determination, and no TP samples were therefore available. The TP sample was also missing for one patient. Overall, 55 T0 samples (10 from NRs and 45 from Rs) and 36 TP samples (all from Rs) were analyzed by qRT-PCR.
Consistent with small RNA-seq data, statistically significant differences were observed in miR-1246 and miR-485-3p T0 levels between the group of Rs and NRs (Figure 1A). The median expression level of miR-1246 was 3.84 (IQR; 2.53–6.31) for Rs and 10.59 (IQR: 5.86–16.57) for NRs (p < 0.001), whereas the median expression level of miR-485-3p was 0.036 (IQR: 0.014–0.089) for Rs and 0.011 (IQR: 0.007–0.046) for NRs (p = 0.041). On the other hand, no differences were observed in miR-92b-3p T0 levels between Rs and NRs (Figure 1A).
For the group of Rs, the levels of the three cf-miRNAs at TP were then compared to those at T0, including in the analysis-only matched samples (i.e., 34 T0-TP pairs). The results illustrated in Figure 1B show that at progression, an increase in the level of miR-1246 and a decrease in that of miR-485-3p were detectable, even though the differences did not reach statistical significance. The median expression level of miR-1246 was 4.51 at T0 (IQR: 2.84–6.65) and 5.82 at TP (IQR: 3.45–8.91), whereas the median expression level of miR-485-3p was 0.038 at T0 (IQR: 0.012–0.091) and 0.020 at TP (IQR: 0.013–0.038). Comparable results were obtained when all samples (i.e., 45 T0 and 36 TP samples) were included in the analysis (Figure S2). No substantial differences were observed between TP and T0 samples in the case of miR-92b-3p (Figure 1B and Figure S2).

3.4. Evaluation of the Ability of miR-1246 and miR-485-3p T0 Plasma Levels to Discriminate NRs from Rs and of Their Association with PFS and OS

Based on the results of qRT-PCR data, miR-92b-3p was not further investigated, whereas T0 levels of miR-1246 and miR-485-3p determined by qRT-PCR in all patients were subjected to statistical analyses to assess their potential value as predictive biomarkers of patients’ response to therapy and their association with PFS and OS.
For ROC analysis, we considered each cf-miRNA alone and the miR-1246/miR485-3p ratio as a single parameter reflecting the levels of the two miRNAs in each patient. The optimal cut-off for sensitivity and specificity was 8.637 for miR-1246, 0.013 for miR-485-3p and 345.493 for miR-1246/miR485-3p ratio. As shown in Table 6, a good ability to discriminate between Rs and NRs was confirmed for miR-1246, even though the AUC was lower than that resulting from ROC analysis performed on miR-1246 plasma levels obtained by small RNA-seq. An AUC > 0.8 was also displayed by the miR-1246/miR485-3p ratio, whereas miR-485-3p showed a less satisfactory predictive performance.
Based on the threshold identified by Youden’s index, patients were then categorised according to T0 levels of miR-1246, miR-485-3p, their combinations and the miR-1246/miR485-3p ratio. Twelve (21.8%) out of fifty-five patients showed miR-1246 levels ≥ 8.637, whereas 16 patients (29.1%) showed miR-485-3p levels < 0.013. Combining these two cf-miRNA thresholds, group (a) included 34 (61.8%) patients with miR-1246 < 8.637 and miR-485-3p ≥ 0.013; group (b) included 14 (25.5%) patients with either miR-1246 ≥ 8.637 and miR-485-3p ≥ 0.013, or miR-1246 < 8.637 and miR-485-3p < 0.013; group (c) included 7 (12.7%) patients with miR-1246 ≥ 8.637 and miR-485-3p < 0.013. Seventeen (30.9%) out of fifty-five patients showed a miR-1246/miR-485-3p ratio ≥ 345.493.
Median PFS of all 55 patients was 8.3 months. Patients with miR-1246 < 8.637 showed a median PFS of 8.9 months, whereas patients with miR-1246 ≥ 8.637 showed a median PFS of 4.8 months (p = 0.003) (Figure 2A). Median PFS was 10.2 months for patients with miR-485-3p ≥ 0.013, compared with 5.4 months for patients with miR-485-3p < 0.013 (p < 0.001) (Figure 2B). When patients were classified according to the combinations of miR-1246 and miR-485-3p levels, median PFS was 10.0 months for group (a), 7.2 months for group (b), and 4.3 months for group (c) (log-rank test for trend, p < 0.001) (Figure 2C). Patients with an miR-1246/miR-485-3p ratio ≥ 345.493 also showed a worse PFS in comparison with those with a ratio < 345.493 (median PFS 10.6 vs. 5.3 months, p < 0.001) (Figure 2D). The percentages of patients without progression at 12 months from the start of therapy were also determined. As illustrated in Figure 2 for each analysis, the percentage of patients free of progression ranged between 0% and 20% in the group of those having unfavorable levels of miR-1246 and/or miR-485-3p or a high miR-1246/miR-485-3p ratio, whereas about 40% of the patients with favorable levels of both cf-miRNAs or of their ratio were still in response at that time.
In univariate Cox regression analysis (Table 7), an increased risk of progression was found for patients with miR-1246 ≥ 8.637 in comparison to patients with miR-1246 < 8.637, as well as for patients with miR-485-3p < 0.013 in comparison to patients with miR-485-3p ≥ 0.013. An increased risk of progression was also found for patients of groups (b) and (c) as compared with patients of group (a), and for patients with miR-1246/miR-485-3p ratio ≥ 345.493 in comparison to patients with a ratio < 345.493.
Univariate Cox regression analysis for PFS was also performed for patients’ baseline characteristics of known prognostic value, namely age, sex, disease stage, sLDH level, and Eastern Cooperative Oncology Group (ECOG) performance status. As illustrated in Table 7, male gender, stage M1c, elevated sLDH, and ECOG performance status ≥ 1 were all associated with an increased risk of progression in our cohort of patients.
Median OS, assessed for all 38 patients remaining after exclusion of those who changed therapy after disease progression, was 13.0 months. A poorer survival was observed for patients with miR-1246 levels ≥8.637 vs. <8.637 (median: 6.3 vs. 13.5 months, p = 0.013), (Figure 3A) and for patients with miR-485-3p levels <0.013 vs. ≥0.013 (median: 7.5 vs. 15.1 months, p = 0.026) (Figure 3B). Due to the limited number of patients in groups (b) and (c) out of the 38 patients available for OS analysis (9 and 7, respectively), the two groups were considered together for Kaplan–Meier estimates. A worse OS was also observed for patients of groups (b) and (c) vs. patients of group (a) (median: 7.5 vs. 20.5 months, log-rank test, p = 0.004) (Figure 3C) and for patients with miR-1246/miR-485-3p ratio ≥345.493 vs. <345.493 (median 6.8 vs. 16.2 months, p = 0.005) (Figure 3D). Univariate Cox regression analysis showed that miR-1246 levels ≥ 8.637, miR-485-3p levels < 0.013, the (b) + (c) combination, and miR-1246/miR-485-3p ratio ≥ 345.493 were all associated with increased risk of mortality (Table 8). As observed for PFS, male gender, stage M1c, elevated sLDH, and ECOG performance status ≥ 1 were also associated with increased risk of mortality (Table 8).

3.5. Evaluation of miR-1246 and miR-485-3p T0 Plasma Level Association with PFS and Mortality in Multivariate Cox Regression Analysis

We next assessed whether T0 levels of miR-1246, T0 levels of miR485-3p, or their combination or ratio were independent prognostic factors for PFS and/or OS. To this end, four separate multivariate Cox models were carried out for both PFS and mortality. Multivariate analysis included only those variables that were statistically significant in univariate analysis. Since sLDH was strongly associated with disease stage and ECOG performance status (Fisher’s exact test, p < 0.0001 and p = 0.016, respectively), to avoid multicollinearity we included in each model only sex and sLDH.
An interaction effect was observed between categories of miR-1246 and sLDH levels for PFS (HR = 0.10; 95% CI: 0.01–1.05, p = 0.055) and OS (HR = 0.04; 95% CI: 0.003–0.46, p = 0.011). Therefore, a stratified analysis by sLDH level was performed for miR-1246 levels (Table S4). After adjusting for sex, no association between miR-1246 and PFS was found in patients with elevated sLDH. Only one patient fell into the class of those with normal sLDH and miR-1246 plasma levels ≥ 8.637, and consequently, even if the multivariate analysis showed an increased risk of progression for this category of patients, the confidence interval was very wide (HR = 5.02; 95% CI: 0.53–47.5) and no statistical significance was achieved. No association between miR-1246 and OS was found for patients with either normal or elevated LDH.
After adjusting for sex and sLDH, miR-485-3p levels < 0.013 remained independently associated with a two-fold increased risk of progression (HR = 2.43; 95% CI 1.21–4.90, p = 0.013 compared with miR-485-3p levels ≥ 0.013) (Table S4). Considering the combination of miR-1246 and miR-485-3p, patients included in group (c), compared to patients in group (a), also reported a higher risk of progression, although statistical significance was not achieved (HR = 2.70; 95% CI 0.96–7.61, p = 0.060) (Table S4). An miR-1246/miR-485-3p ratio ≥ 345.493 also remained independently associated with a three-fold increased risk of progression (HR = 2.94, 95% CI: 1.47–5.88, p = 0.002) (Table S4).
After adjusting for sex and sLDH levels, only an miR-1246/miR-485-3p ratio ≥ 345.493 showed an increased risk of mortality (HR = 2.28; 95%CI: 0.90–5.80, p = 0.082), although this association was not statistically significant.
Elevated sLDH levels remained statistically significant in all multivariate models, with an increased risk for both progression and mortality. The strong association of sLDH levels with PFS and OS, prompted us to evaluate the usefulness of sLDH levels to predict response to therapy in our patient cohort. No statistically significant association was found between sLDH levels and response to therapy (Fisher’s exact test, p = 0.157).

4. Discussion

Circulating cf-miRNAs are emerging as valuable non-invasive biomarkers for cancer diagnosis and assessment of patients’ prognosis and response to therapy [15,16,26,28]. In this regard, most of the studies performed so far on circulating cf-miRNAs in melanoma have been focused on their role as potential diagnostic and prognostic biomarkers in treatment-naïve patients [30,31], while few investigations have been performed to assess the association of cf-miRNA expression patterns with clinical outcomes in patients subjected to immunotherapy [48,49] or targeted therapy [21,32].
In the present study, using an NGS approach to profile cf-miRNAs in a cohort of 33 patients, we identified a set of cf-miRNAs potentially linked to resistance to therapy with BRAFi and MEKi, as they were DE at baseline between patients who had responded to treatment and patients who had not, and/or at disease progression with respect to baseline, in patients who had initially achieved an objective clinical response. Among those cf-miRNAs, baseline levels of a more restricted set also showed a good ability in discriminating Rs and NRs and/or a significant association with the duration of response to treatment, suggesting their potential utility in the clinical setting. We considered miR-1246, miR-92b-3p, and miR485-3p interesting for further studies, and therefore determined their plasma level at T0 and TP by qRT-PCR in the initial cohort of 33 patients and in an additional 24 patients. The results of qRT-PCR assays were consistent with small RNA-seq data analysis for miR-1246 and miR-485-3p, but not for miR-92b-3p, which therefore was not further investigated.
A large body of experimental evidence indicates that miR-1246 acts as an oncomiR in most, although not all, types of cancer, and that its expression in the tumor tissue and/or in patients’ plasma/serum can have a diagnostic and prognostic value [50]. Indeed, with respect to normal adjacent tissue, overexpression of miR-1246 has been described in breast, colorectal, lung and ovarian cancer, as well as in hepatocellular and oral squamous cell carcinoma [50,51]. Moreover, increased levels of miR-1246 have been detected in plasma/serum of patients with those types of cancer as compared with healthy controls [50]. In patients with colorectal, lung, esophageal and breast cancer and hepatocellular and oral squamous cell carcinoma, worse disease-free survival and/or OS were also found to be associated with high miR-1246 expression in the tumor tissue and/or in plasma/serum [50,51]. An elevated level of circulating miR-1246 was also demonstrated to be a marker of therapy resistance in colorectal and breast cancer [52,53]. On the other hand, decreased expression of miR-1246 has been found in prostate [54] and renal cell carcinoma tissues [55], suggesting a tumor suppressor function of miR-1246 in these types of neoplasias. Notably, both upregulation [56] and downregulation [57] of miR-1246 have been reported in cervical carcinoma specimens.
Numerous functional studies have been performed in cell lines derived from different types of tumors to assess the impact of miR-1246 overexpression or silencing on proliferation, survival, invasiveness and chemoresistance, and to identify the most relevant targets involved in the observed effects. For instance, regarding the oncogenic function of miR-1246, Wang et al. [58] demonstrated that this miRNA promoted proliferation and invasiveness of colorectal cancer (CRC) cell lines and protected them from apoptosis through negative regulation of the CCNG2 gene, which encodes cycling G2, an atypical cyclin that induces cell cycle arrest and is frequently downregulated in tumors. miR-1246-mediated downregulation of cycling G2 expression was also demonstrated to increase the stemness-like properties and chemoresistance of pancreatic and oral cancer cells [59,60], as well as to enhanced proliferation, invasiveness and drug resistance of breast cancer cells [61]. In vitro studies addressing the role of miR-1246 in the crosstalk between cancer cells and the tumor microenviroment demonstrated that this miRNA can shuttle between CRC cells and fibroblasts, promoting migration of the tumor cells—via activation of the Wnt/b-catenin signaling pathway—and transdifferentiation of fibroblasts into cancer-associated fibroblasts that support tumor progression [62]. Similarly, Cooks et al. [63] demonstrated that CRC cells with specific gain of function mutations in p53 can reprogram neighboring macrophages into a tumor-promoting state by releasing miR-1246-enriched exosomes. Several other target genes of miR-1246 involved in its oncogenic activity have been experimentally identified, including, among others, GSK3B (Glycogen Synthase Kinase 3 Beta) in lung cancer cells, SPRED2 (Sprouty Related EVH1 Domain Containing 2) in CRC cells, CADM1 (Cell Adhesion Molecule 1) in hepatocellular carcinoma cells [50], NFE2L3 (NFE2-like bZIP transcription factor 3) in breast cancer cells [51].
A limited number of studies have also reported inhibitory effects of miR-1246 on tumor cell growth and/or invasiveness. In renal carcinoma cell lines, miR-1246 was shown to target CXCR4 (C-X-C Motif chemokine receptor 4) leading to impairment of proliferation and migration [55], while in prostate cancer cells, CDH2 (Cadherin 2) and VIM (Vimentin) were identified as direct miR-1246 target genes involved in its suppression of epithelial to mesenchimal transition (EMT) [54]. Controversial results have instead been reported in cervical cancer cells. Indeed, Chen and collaborators demonstrated that miR-1246 suppressed the THBS2 (Thrombospondin 2) gene and positively modulated cell proliferation, migration and invasion [64], whereas inhibitory effects of miR-1246 on cervical cancer cell invasiveness were reported by Yang et al. [57].
Compared to miR-1246, fewer investigations have addressed the role of miR-485-3p in cancer, but, with few exceptions [65,66,67], it appears to function as a tumor-suppressor miRNA. Its expression in tumor tissue was found to be downregulated in prostate cancer [68,69], glioblastoma [70], breast and colorectal cancer [71,72,73,74], osteosarcoma [75] and renal carcinoma [76]. Interestingly, in some of those studies, the reduced expression of miR-485-3p was linked to upregulation of a long non-coding RNA (LncRNA) [75] or a circular RNA [72,74,76] able to sponge the miRNA, thus leading to enhanced expression of target genes promoting tumor cell proliferation, survival and invasiveness. In the circulation, lower levels of miR-485-3p were detected in serum and serum exosome of glioblastoma patients as compared with healthy subjects [77,78]. Moreover, low pre-surgery levels of miR-485-3p were associated with shorter PFS and OS in glioma patients receiving radiotherapy plus chemotherapy after surgery [77].
As for miR-1246, functional studies have been performed in cancer cell lines to identify the molecular mechanisms involved in the biological activity of miR-485-3p. Several targets of this miRNA have been experimentally validated, and their downregulation has been implicated in the inhibitory effects exerted by the miRNA on tumor cell proliferation, survival and invasiveness. For instance, in prostate cancer cells, miR-485-3p was shown to inhibit proliferation, migration and invasion by targeting the TGFBR2 (Transforming Growth Factor Beta Receptor 2) gene, a key regulator of the TGF-β signaling pathway [69], while in glioblastoma cells, the inhibitory effects of miR-485-3p on proliferation and migration were linked to downregulation of the RNF135 (Ring Finger Protein 135) gene and impairment of the ERK1/2 signaling pathway [70]. In breast cancer cells, miR-485-3p was found to suppress migration, invasion and mitochondrial respiration by inhibiting the expression of the PGC-1A gene, which encodes a transcriptional coactivator that regulates genes involved in energy metabolism [71]. Other target genes whose negative modulation by miR-485-3p has been shown to impair tumor growth and metastasis include, among others, ZEB1 (Zinc Finger E-box Binding Homeobox 1) and BIRC5 (Baculoviral IAP Repeat Containing 5) in breast cancer cells [72,79], MELK (Maternal Embryonic Leucine Zipper Kinase) and JAK2 (Janus Kinase 2) in CRC cells [74,80], MET and AKT3, in osteosarcoma cells [75].
An oncogenic function of miR-485-3p has also been reported. Specifically, miR-485-3p was found to be upregulated in hepatocellular carcinoma, and to support tumor cell proliferation and survival in vitro and tumor growth and metastasis in mice by targeting the MAT1A gene, which encodes the α1 catalytic subunit of methionine adenosyltransferase [65]. In hepatocellular carcinoma, miR-485-3p was shown to also promote proliferation and invasiveness by downregulation of NTRK3, which codes for the neurotrophic tyrosine kinase receptor type 3 [67]. Moreover, suppression of this target gene was demonstrated to underlie the tumor-promoting activity of miR-485-3p in gastric cancer [81].
Few data are presently available regarding the expression and function of miR-1246 and miR-485-3p in melanoma. Armand-Labit et al. [82] reported that miR-1246 plasma levels were significantly higher in metastatic melanoma patients than in healthy controls, and confirmed the expression of this miRNA in melanoma metastases. Furthermore, the authors showed that plasma levels of miR-1246 in combination with those of miR-185 could differentiate patients from healthy controls with elevated accuracy. Increased levels of miR-1246 were also observed by Torii et al. [41] in serum EVs derived from melanoma patients as compared with those isolated from healthy controls. Interestingly, the authors demonstrated that miR-1246 contained in EVs derived from a highly metastatic melanoma cell line, as well as miR-1246 mimics, were able to increase resistance to 5-fluorouracil in endothelial cells. Enhanced expression of miR-1246 in melanoma specimens with respect to normal tissue was recently reported by Yu et al. [39], who in addition demonstrated that miR-1246 promoted melanoma cell proliferation, survival and invasiveness by targeting FOXA2, a gene with a tumor-suppressor function in several types of cancer, including melanoma [83,84,85,86]. Finally, miR-1246 was found to be upregulated in melanoma cell lines resistant to the BRAFi PLX4720 [40]. Regarding miR-485-3p, there is only a very recent study by Huo et al. [47], who reported over-expression of the LncRNA MIR155HG and downregulation of miR-485-3p in melanoma specimens, and demonstrated that the LncRNA acts as a molecular sponge of miR-485-3p, leading to increased expression of its target gene PSIP1. This gene is upregulated in several cancers and involved in tumor aggressiveness and chemoresistance [87,88,89,90]. Overall, those studies suggest that in melanoma, miR-1246 and miR-485-3p can have an oncogenic and a tumor-suppressor function, respectively.
In the present investigation, we show for the first time that baseline plasma levels of miR-1246 and miR-485-3p could represent valuable and non-invasive biomarkers to predict clinical response and prognosis in melanoma patients treated with BRAFi and MEKi. Indeed, miR-1246 and miR-485-3p baseline (i.e., T0) plasma levels were significantly higher and lower, respectively, in the group of NRs as compared with the group of Rs, and a trend toward an increase in miR-1246 and a decrease in miR-485-3p was observed in the latter group of patients at the development of secondary resistance (i.e., TP). Both baseline miR-1246 levels and the miR-1246/miR-485-3p ratio displayed a good ability to discriminate between Rs and NRs. Moreover, Kaplan–Meier curves showed that either miR-1246 plasma levels ≥ 8.637 or miR-485-3p plasma levels < 0.013 were associated with poorer PFS and OS, while univariate Cox regression analysis evidenced a more than two-fold higher risk of short-term progression and mortality for the two groups of patients showing those cf-miRNA levels. Accordingly, when the different combinations of miR-1246 and miR-485-3p were considered, the group of patients with favourable levels of both cf-miRNAs showed the best PFS and OS. Notably, the miR-1246/miR-485-3p ratio also appeared to be a useful prognostic parameter for PFS and OS. Indeed, an miR-1246/miR-485-3p ratio ≥ 345.493 was associated with significantly shorter PFS and OS and higher risk of progression and mortality. Altogether, our results are in agreement with previous studies that suggest an oncogenic and tumor-suppressor function in melanoma for miR-1246 and miR-485-3p, respectively.
Our findings also appear to be of clinical relevance. Indeed, both targeted therapy and immunotherapy with immune checkpoint inhibitors are possible therapeutic options for patients with BRAF-mutant melanoma, and the availability of novel biomarkers able to predict patients’ responses to therapy and to estimate their expectancy of PFS and OS may potentially help clinicians choose the optimal therapeutic protocol. The potential clinical utility of assessing baseline plasma levels of miR-1246 and miR-485-3p is also strengthened by the finding that miR-485-3p and the miR-1246/miR-485-3p ratio remained independently associated with PFS after adjusting for sex and sLDH in multivariate analysis. After adjusting for sex and sLDH, patients with a miR-1246/miR-485-3p ratio ≥ 345.493 showed a two-fold increased risk of mortality, but the statistical significance was not reached, most probably for the limited size of our cohort of patients. On the other hand, the stratified analysis by sLDH performed for miR-1246 levels did not evidence any significant association between miR-1246 and PFS or OS for patients with either normal or elevated sLDH, further highlighting that the simultaneous determination of miR-1246 and miR-485-3p can provide a better prognostication of patients’ outcomes.
LDH, which catalyzes the conversion of pyruvate to lactate, is a key enzyme of the glycolytic pathway and plays a central role in the Warburg effect, namely the elevated aerobic glycolysis that is considered a hallmark of cancer [91]. Accordingly, numerous studies in different types of cancer have demonstrated the importance of LDH in tumor growth and malignant behavior [91]. LDH can be released into blood stream upon cell damage, and in cancer patients sLDH is considered a biomarker linked to tumor burden and aggressiveness [91,92]. Elevated sLDH represents an independent predictor of poor outcome in patients with stage IV melanoma, and presently, sLDH is the only circulating biomarker incorporated by the AJCC in melanoma staging [1]. High levels of sLDH are also a predictor of inferior OS and PFS in melanoma patients receiving BRAFi with or without MEKi [93,94]. These latter findings were confirmed in our study, which also evidenced an interaction effect between miR-1246 and sLDH for PFS and OS, most likely because plasma levels of this miRNA can be an expression of tumor burden and biological aggressiveness, as reported for sLDH. It is, however, noteworthy that although baseline sLDH is a strong biomarker of PSF and OS in patients receiving targeted therapy, in contrast to circulating miR-1246 and the miR-1246/miR-485-3p ratio, it does not appear useful for predicting patient response to treatment. In this regard, the analysis of pooled data derived from three phase 3 clinical trials of dabrafenib + trametinib in melanoma patients evidenced that, although inferior with respect to patients with normal sLDH levels, the response rate was also elevated (about 50%) in patients with high sLDH levels [95]. Similarly, a response rate of 63% was reported by a large multicentric study in melanoma patients with elevated sLDH treated with targeted therapy [96]. Consistent with these findings, we did not find any association between baseline sLDH and response to therapy in our cohort of melanoma patients.
A previous study by Svedman et al. [32] evaluated the levels of 372 miRNAs in plasma EVs isolated before and during treatment from melanoma patients receiving therapy with BRAFi alone or in combination with MEKi. Twenty patients who had achieved PR or SD constituted the group of “disease control”, whereas eight patients progressing on therapy without any previous response constituted the group of NRs. The authors observed that increased levels of let-7g-5p during the course of treatment were associated with better disease control, and that patients displaying high levels of miR-497-5p during therapy had a longer PFS; but the authors did not find any association between pre-treatment EV miRNA levels and response to therapy. Several differences exist between the study of Svedman and collaborators and our investigation that could explain why the former did not identify pre-treatment miRNA levels associated with patients’ response to therapy. We performed cf-miRNA profiling using RNA extracted from total plasma samples, and therefore, in addition to miRNAs included in EVs, those associated with proteins and lipoproteins were also detected, increasing the possibility of identifying a clinically relevant cf-miRNA. Notably, miR-1246 was not present in the 372-miRNA set analyzed by Svedman and collaborators. Furthermore, we classified patients as R and NRs, including in the former group those who had achieved CR and PR, and in the latter group those who had shown SD or PD as their best clinical response, according to RECIST 1.1 criteria. This different classification of patients may have also favoured the identification of cf-miRNAs baseline levels associated with patient response to therapy.
The potential of two cf-miRNAs, namely miR-199b-5p and miR-4488, as biomarkers of melanoma patients’ resistance to targeted therapy was also investigated by Fattore et al. [21] in a cohort of 25 patients treated with vemurafenib or the combination of vemurafenib + cobimetinib or encorafenib + binimetinib. Consistent with their findings in melanoma cell lines sensitive or resistant to targeted therapy, and in melanoma specimens collected before the start of treatment and at progression, the authors found that miR-199b-5p and miR-4488 plasma levels were downregulated and upregulated, respectively, at disease progression as compared to baseline, and showed a good ability to discriminate pre-treatment samples from progression samples. According to small RNA-seq data analysis, these two miRNAs were not DE at TP vs. T0 in our cohort of Rs, but it is possible that this result depends on the filtering procedures and statistical analyses applied to small RNA-seq row data. Moreover, in contrast with the study of Fattore and collaborators, most of our patients were subjected to dabrafenib or dabrafenib + trametinib, and it cannot be excluded that the changes occurring in cf-miRNA expression patterns from pre-treatment to progression may be related to the type of BRAFi and MEKi received from the patients.
Interestingly, in a number of studies [97,98,99], miR-1246 in serum, plasma or exosomes was shown to be derived from the processing of the RNU2-1 transcript, a small nuclear RNA, and not from the MIR1246 gene through to the canonical pathway of miRNA biogenesis. Functionality of miR-1246 derived from RNU2-1 transcript was confirmed [99], and therefore, independently of its origin, miR-1246 remains an important miRNA involved in cancer. Future studies are required to ascertain whether circulating miR-1246 in our cohort of patients derives from the processing of RNU2-1 transcript.
We are aware that our study has some limitations. Treatment with BRAFi and MEKi in our cohort of melanoma patients was eterogenous, and it can not be excluded that this could have affected survival outcomes. Moreover, the study was retrospective, and included a relatively small number of patients. Although our results appear of clinical relevance, they need to be interpreted cautiously and validated in a larger prospective study.

5. Conclusions

Treatment with inhibitors of the BRAF/MEK pathways provides significant benefit in patients with BRAF-mutant advanced melanoma. However, heterogeneity in treatment response and its duration highlights the need for biomarkers that can allow a more personalized therapeutic approach.
Here, we showed that miR-1246 and miR-485-3p baseline plasma levels are associated with clinical response and prognosis in melanoma patients treated with BRAFi and MEKi. In particular, the baseline miR-1246/miR-485-3p ratio, which takes into consideration the levels of both miRNAs in each individual, appears to be a valuable biomarker to identify patients most likely to be unresponsive to targeted therapy or at higher risk for short-term PFS and mortality, and therefore is of potential clinical utility to improve, in addition to the presently well recognized prognostic factors, the management of patients with BRAF-mutant metastatic melanoma. Studies in a larger cohort of patients are warranted in order to further validate the predictive and prognostic potential of circulating miR-1246 and miR-485-3p.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cancers14153706/s1. Figure S1: miRTrace quality control report; Figure S2: Box-and-whisker diagrams of T0 and TP plasma levels of miR-1246, mir-92b-3p and miR-485-3p in melanoma patients treated with targeted therapy; Table S1: Demographic and clinical–pathological features of each melanoma patient included in the study; Table S2: cf-miRNA expression levels determined by small RNA-seq. Table S3: cf-miRNAs DE in T0 plasma samples of NRs vs. Rs subjected to BRAFi/MONO according to small RNA-seq data analysis; Table S4: Multivariate Cox regression analysis for PFS according to T0 plasma levels of miR-1246, miR-485-3p and clinical features of patients.

Author Contributions

Conceptualization, L.L. (Lauretta Levati) and S.D.; data analysis, L.L. (Lauretta Levati), C.B., S.M. and S.D.; data curation, L.L. (Lauretta Levati), C.B., S.M., G.C.A.C. and S.D.; funding acquisition, G.R. and S.D.; methodology, L.L. (Lauretta Levati), C.B., S.M., L.L. (Laura Lupini), L.B., S.C. and P.M.L.; patient treatment and follow-up, G.C.A.C. and F.D.G.; writing—original draft, L.L. (Lauretta Levati), C.B., S.M. and S.D.; writing—review & editing, L.L. (Lauretta Levati), C.B., S.M., G.C.A.C., E.A., M.G.N., I.M., F.D.G., M.N., G.R. and S.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Italian Association for Cancer Research (AIRC, grant IG 17585 to S.D.), and by the Italian Ministry of Health (grant PE-2016-02364102 to G.R.).

Institutional Review Board Statement

The study was conducted in accordance with Good Clinical Practice Guidelines and the Declaration of Helsinki. The study was approved by the IDI-IRCCS Ethics Committee (ID #407/1, 2013).

Informed Consent Statement

Written informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The small RNA-seq row data generated and analyzed during the current study are available in the NCBI Sequence Read Archive (SRA) database with the accession number PRJNA834668. All other data supporting the findings of this study are available within the article and its Supplementary Information files and from the corresponding author upon reasonable request.

Acknowledgments

The authors thank the nurses Cristina Sanandrea and Arianna Brunello (Department of Oncology and Dermatological Oncology, Istituto Dermopatico dell’Immacolata, IDI-IRCCS) for their collaboration in this study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Keung, E.Z.; Gershenwald, J.E. The Eighth Edition American Joint Committee on Cancer (AJCC) Melanoma Staging System: Implications for Melanoma Treatment and Care. Expert Rev. Anticancer Ther. 2018, 18, 775–784. [Google Scholar] [CrossRef] [PubMed]
  2. Mackenzie, A.; Flaherty, K.T.; Boycott, K.; Robert, C.; Schadendorf, D.; Mastrobuoni, S.; Chan, V.; Ruel, M. MEK Inhibition in BRAF-Mutated Melanoma. N. Engl. J. Med. 2012, 367, 1364–1365. [Google Scholar] [CrossRef] [PubMed]
  3. Queirolo, P.; Picasso, V.; Spagnolo, F. Combined BRAF and MEK Inhibition for the Treatment of BRAF-Mutated Metastatic Melanoma. Cancer Treat. Rev. 2015, 41, 519–526. [Google Scholar] [CrossRef] [PubMed]
  4. Robert, C.; Grob, J.J.; Stroyakovskiy, D.; Karaszewska, B.; Hauschild, A.; Levchenko, E.; Chiarion Sileni, V.; Schachter, J.; Garbe, C.; Bondarenko, I.; et al. Five-Year Outcomes with Dabrafenib plus Trametinib in Metastatic Melanoma. N. Engl. J. Med. 2019, 381, 626–636. [Google Scholar] [CrossRef] [PubMed]
  5. Savoia, P.; Zavattaro, E.; Cremona, O. Clinical Implications of Acquired Braf Inhibitors Resistance in Melanoma. Int. J. Mol. Sci. 2020, 21, 9730. [Google Scholar] [CrossRef] [PubMed]
  6. Patel, M.; Eckburg, A.; Gantiwala, S.; Hart, Z.; Dein, J.; Lam, K.; Puri, N. Resistance to Molecularly Targeted Therapies in Melanoma. Cancers 2021, 13, 1115. [Google Scholar] [CrossRef]
  7. Bartel, D.P. Metazoan MicroRNAs. Cell 2018, 173, 20–51. [Google Scholar] [CrossRef] [Green Version]
  8. Slack, F.J.; Chinnaiyan, A.M. The Role of Non-Coding RNAs in Oncology. Cell 2019, 179, 1033–1055. [Google Scholar] [CrossRef]
  9. Smolarz, B.; Durczyński, A.; Romanowicz, H.; Szyłło, K.; Hogendorf, P. MiRNAs in Cancer (Review of Literature). Int. J. Mol. Sci. 2022, 23, 2805. [Google Scholar] [CrossRef] [PubMed]
  10. Geretto, M.; Pulliero, A.; Rosano, C.; Zhabayeva, D.; Bersimbaev, R.; Izzotti, A. Resistance to Cancer Chemotherapeutic Drugs Is Determined by Pivotal MicroRNA Regulators. Am. J. Cancer Res. 2017, 7, 1350–1371. [Google Scholar]
  11. Corrà, F.; Agnoletto, C.; Minotti, L.; Baldassari, F.; Volinia, S. The Network of Non-Coding RNAs in Cancer Drug Resistance. Front. Oncol. 2018, 8, 327. [Google Scholar] [CrossRef] [Green Version]
  12. He, B.; Zhao, Z.; Cai, Q.; Zhang, Y.; Zhang, P.; Shi, S.; Xie, H.; Peng, X.; Yin, W.; Tao, Y.; et al. Mirna-Based Biomarkers, Therapies, and Resistance in Cancer. Int. J. Biol. Sci. 2020, 16, 2628–2647. [Google Scholar] [CrossRef] [PubMed]
  13. Rupaimoole, R.; Slack, F.J. MicroRNA Therapeutics: Towards a New Era for the Management of Cancer and Other Diseases. Nat. Rev. Drug Discov. 2017, 16, 203–221. [Google Scholar] [CrossRef] [PubMed]
  14. Winkle, M.; El-Daly, S.M.; Fabbri, M.; Calin, G.A. Noncoding RNA therapeutics—Challenges and potential solutions. Nat. Rev. Drug Discov. 2021, 20, 629–651. [Google Scholar] [CrossRef] [PubMed]
  15. Terrinoni, A.; Calabrese, C.; Basso, D.; Aita, A.; Caporali, S.; Plebani, M.; Bernardini, S. The Circulating MiRNAs as Diagnostic and Prognostic Markers. Clin. Chem. Lab. Med. 2019, 57, 932–953. [Google Scholar] [CrossRef] [Green Version]
  16. Zhou, E.; Li, Y.; Wu, F.; Guo, M.; Xu, J.; Wang, S.; Tan, Q.; Ma, P.; Song, S.; Jin, Y. Circulating Extracellular Vesicles Are Effective Biomarkers for Predicting Response to Cancer Therapy. EBioMedicine 2021, 67, 103365. [Google Scholar] [CrossRef] [PubMed]
  17. Pinto, R.; Strippoli, S.; De Summa, S.; Albano, A.; Azzariti, A.; Guida, G.; Popescu, O.; Lorusso, V.; Guida, M.; Tommasi, S. MicroRNA Expression in BRAF-Mutated and Wild-Type Metastatic Melanoma and Its Correlation with Response Duration to BRAF Inhibitors. Expert Opin. Ther. Targets 2015, 19, 1027–1035. [Google Scholar] [CrossRef]
  18. Fattore, L.; Costantini, S.; Malpicci, D.; Ruggiero, C.F.; Ascierto, P.A.; Croce, C.M.; Mancini, R.; Ciliberto, G. MicroRNAs in Melanoma Development and Resistance to Target Therapy. Oncotarget 2017, 8, 22262–22278. [Google Scholar] [CrossRef] [Green Version]
  19. Romano, G.; Kwong, L.N. MiRNAs, Melanoma and Microenvironment: An Intricate Network. Int. J. Mol. Sci. 2017, 18, 2354. [Google Scholar] [CrossRef] [Green Version]
  20. Ross, C.L.; Kaushik, S.; Valdes-Rodriguez, R.; Anvekar, R. MicroRNAs in Cutaneous Melanoma: Role as Diagnostic and Prognostic Biomarkers. J. Cell. Physiol. 2018, 233, 5133–5141. [Google Scholar] [CrossRef]
  21. Fattore, L.; Ruggiero, C.F.; Pisanu, M.E.; Liguoro, D.; Cerri, A.; Costantini, S.; Capone, F.; Acunzo, M.; Romano, G.; Nigita, G.; et al. Reprogramming MiRNAs Global Expression Orchestrates Development of Drug Resistance in BRAF Mutated Melanoma. Cell Death Differ. 2019, 26, 1267–1282. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Caporali, S.; Amaro, A.; Levati, L.; Alvino, E.; Lacal, P.M.; Mastroeni, S.; Ruffini, F.; Bonmassar, L.; Antonini Cappellini, G.C.; Felli, N.; et al. MiR-126-3p down-Regulation Contributes to Dabrafenib Acquired Resistance in Melanoma by up-Regulating ADAM9 and VEGF-A. J. Exp. Clin. Cancer Res. 2019, 38, 272. [Google Scholar] [CrossRef] [PubMed]
  23. Tupone, M.G.; D’Aguanno, S.; Di Martile, M.; Valentini, E.; Desideri, M.; Trisciuoglio, D.; Donzelli, S.; Sacconi, A.; Buglioni, S.; Ercolani, C.; et al. MicroRNA-378a-5p IS a Novel Positive Regulator of Melanoma Progression. Oncogenesis 2020, 9, 22. [Google Scholar] [CrossRef] [PubMed]
  24. Motti, M.L.; Minopoli, M.; Di Carluccio, G.; Ascierto, P.A.; Carriero, M.V. Micrornas as Key Players in Melanoma Cell Resistance to Mapk and Immune Checkpoint Inhibitors. Int. J. Mol. Sci. 2020, 21, 4544. [Google Scholar] [CrossRef]
  25. Lorusso, C.; De Summa, S.; Pinto, R.; Danza, K.; Tommasi, S. MiRNAs as Key Players in the Management of Cutaneous Melanoma. Cells 2020, 9, 415. [Google Scholar] [CrossRef] [Green Version]
  26. Varrone, F.; Caputo, E. The MiRNAs Role in Melanoma and in Its Resistance to Therapy. Int. J. Mol. Sci. 2020, 21, 878. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Sloane, R.A.S.; White, M.G.; Witt, R.G.; Banerjee, A.; Davies, M.A.; Han, G.; Burton, E.; Ajami, N.; Simon, J.M.; Bernatchez, C.; et al. Identification of MicroRNA-MRNA Networks in Melanoma and Their Association with PD-1 Checkpoint Blockade Outcomes. Cancers 2021, 13, 5301. [Google Scholar] [CrossRef]
  28. Nguyen, M.H.T.; Luo, Y.H.; Li, A.L.; Tsai, J.C.; Wu, K.L.; Chung, P.J.; Ma, N. Mirna as a Modulator of Immunotherapy and Immune Response in Melanoma. Biomolecules 2021, 11, 1648. [Google Scholar] [CrossRef] [PubMed]
  29. Ghafouri-Fard, S.; Gholipour, M.; Taheri, M. MicroRNA Signature in Melanoma: Biomarkers and Therapeutic Targets. Front. Oncol. 2021, 11, 608987. [Google Scholar] [CrossRef]
  30. Mumford, S.L.; Towler, B.P.; Pashler, A.L.; Gilleard, O.; Martin, Y.; Newbury, S.F. Circulating MicroRNA Biomarkers in Melanoma: Tools and Challenges in Personalised Medicine. Biomolecules 2018, 8, 21. [Google Scholar] [CrossRef] [Green Version]
  31. Polini, B.; Carpi, S.; Romanini, A.; Breschi, M.C.; Nieri, P.; Podestà, A. Circulating Cell-Free MicroRNAs in Cutaneous Melanoma Staging and Recurrence or Survival Prognosis. Pigment Cell Melanoma Res. 2019, 32, 486–499. [Google Scholar] [CrossRef] [PubMed]
  32. Svedman, F.C.; Lohcharoenkal, W.; Bottai, M.; Brage, S.E.; Sonkoly, E.; Hansson, J.; Pivarcsi, A.; Eriksson, H. Extracellular Microvesicle MicroRNAs as Predictive Biomarkers for Targeted Therapy in Metastastic Cutaneous Malignant Melanoma. PLoS ONE 2018, 13, e0206942. [Google Scholar] [CrossRef] [PubMed]
  33. Eisenhauer, E.A.; Therasse, P.; Bogaerts, J.; Schwartz, L.H.; Sargent, D.; Ford, R.; Dancey, J.; Arbuck, S.; Gwyther, S.; Mooney, M.; et al. New Response Evaluation Criteria in Solid Tumours: Revised RECIST Guideline (Version 1.1). Eur. J. Cancer 2009, 45, 228–247. [Google Scholar] [CrossRef] [PubMed]
  34. Benjamini, Y.; Hochberg, J. Controlling the False Discovery Rate: A Pratical and Powerful Approach to Multiple Testing. J. R. Stat. Soc. Ser. B 1995, 57, 289–300. [Google Scholar] [CrossRef]
  35. Zou, K.H.; O’Malley, A.J.; Mauri, L. Receiver-Operating Characteristic Analysis for Evaluating Diagnostic Tests and Predictive Models. Circulation 2007, 115, 654–657. [Google Scholar] [CrossRef] [Green Version]
  36. Youden, W.J. Index for Rating Diagnostic Tests. Cancer 1950, 3, 32–35. [Google Scholar] [CrossRef]
  37. Cuzick, J. A Wilcoxon-type Test for Trend. Stat. Med. 1985, 4, 87–90. [Google Scholar] [CrossRef]
  38. English, P.A.; Williams, J.A.; Martini, J.-F.; Motzer, R.J.; Valota, O.; Buller, R.E. A Case for the Use of Receiver Operating Characteristic Analysis of Potential Clinical Efficacy Biomarkers in Advanced Renal Cell Carcinoma. Futur. Oncol. 2016, 12, 175–182. [Google Scholar] [CrossRef] [Green Version]
  39. Yu, Y.; Yu, F.; Sun, P. MicroRNA-1246 Promotes Melanoma Progression through Targeting FOXA2. Onco. Targets. Ther. 2020, 13, 1245–1253. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Kim, J.-H.; Ahn, J.-H.; Lee, M. Upregulation of MicroRNA-1246 Is Associated with BRAF Inhibitor Resistance in Melanoma Cells with Mutant BRAF. Cancer Res. Treat. 2017, 49, 947–959. [Google Scholar] [CrossRef] [Green Version]
  41. Torii, C.; Maishi, N.; Kawamoto, T.; Morimoto, M.; Akiyama, K.; Yoshioka, Y.; Minami, T.; Tsumita, T.; Alam, M.T.; Ochiya, T.; et al. MiRNA-1246 in Extracellular Vesicles Secreted from Metastatic Tumor Induces Drug Resistance in Tumor Endothelial Cells. Sci. Rep. 2021, 11, 13502. [Google Scholar] [CrossRef] [PubMed]
  42. Gong, L.; Ren, M.; Lv, Z.; Yang, Y.; Wang, Z. MiR-92b-3p Promotes Colorectal Carcinoma Cell Proliferation, Invasion, and Migration by Inhibiting FBXW7 In Vitro and In Vivo. DNA Cell Biol. 2018, 37, 501–511. [Google Scholar] [CrossRef] [PubMed]
  43. Wang, G.; Cheng, B.; Jia, R.; Tan, B.; Liu, W. Altered Expression of MicroRNA-92b-3p Predicts Survival Outcomes of Patients with Prostate Cancer and Functions as an Oncogene in Tumor Progression. Oncol. Lett. 2021, 21, 4. [Google Scholar] [CrossRef] [PubMed]
  44. Du, Y.; Miao, Z.; Wang, K.; Lv, Y.; Qiu, L.; Guo, L. Expression Levels and Clinical Values of MiR-92b-3p in Breast Cancer. World J. Surg. Oncol. 2021, 19, 239. [Google Scholar] [CrossRef] [PubMed]
  45. Li, M.; Shan, W.; Hua, Y.; Chao, F.; Cui, Y.; Lv, L.; Dou, X.; Bian, X.; Zou, J.; Li, H.; et al. Exosomal MiR-92b-3p Promotes Chemoresistance of Small Cell Lung Cancer Through the PTEN/AKT Pathway. Front. Cell Dev. Biol. 2021, 9, 661602. [Google Scholar] [CrossRef] [PubMed]
  46. Zhao, F.; Yang, Z.; Gu, X.; Feng, L.; Xu, M.; Zhang, X. MiR-92b-3p Regulates Cell Cycle and Apoptosis by Targeting CDKN1C, Thereby Affecting the Sensitivity of Colorectal Cancer Cells to Chemotherapeutic Drugs. Cancers 2021, 13, 3323. [Google Scholar] [CrossRef]
  47. Huo, J.; Wang, Y.; Zhang, Y.; Wang, W.; Yang, P.; Zhao, W.; Zhang, M.; Cui, L.; Zhang, D. The LncRNA MIR155HG Is Upregulated by SP1 in Melanoma Cells and Drives Melanoma Progression via Modulating the MiR-485-3p/PSIP1 Axis. Anticancer Agents Med. Chem. 2022, 22, 152–159. [Google Scholar] [CrossRef]
  48. Bustos, M.A.; Gross, R.; Rahimzadeh, N.; Cole, H.; Tran, L.T.; Tran, K.D.; Takeshima, L.; Stern, S.L.; O’Day, S.; Hoon, D.S.B. A Pilot Study Comparing the Efficacy of Lactate Dehydrogenase Levels Versus Circulating Cell-Free MicroRNAs in Monitoring Responses to Checkpoint Inhibitor Immunotherapy in Metastatic Melanoma Patients. Cancers 2020, 12, 3361. [Google Scholar] [CrossRef]
  49. Nakahara, S.; Fukushima, S.; Okada, E.; Morinaga, J.; Kubo, Y.; Tokuzumi, A.; Matsumoto, S.; Tsuruta-Kadohisa, M.; Kimura, T.; Kuriyama, H.; et al. MicroRNAs That Predict the Effectiveness of Anti-PD-1 Therapies in Patients with Advanced Melanoma. J. Dermatol. Sci. 2020, 97, 77–79. [Google Scholar] [CrossRef]
  50. Ghafouri-Fard, S.; Khoshbakht, T.; Hussen, B.M.; Taheri, M.; Samadian, M. A Review on the Role of MiR-1246 in the Pathoetiology of Different Cancers. Front. Mol. Biosci. 2022, 8, 771835. [Google Scholar] [CrossRef]
  51. Dai, Y.C.; Pan, Y.; Quan, M.M.; Chen, Q.; Pan, Y.; Ruan, Y.Y.; Sun, J.G. MicroRNA-1246 Mediates Drug Resistance and Metastasis in Breast Cancer by Targeting NFE2L3. Front. Oncol. 2021, 11, 677168. [Google Scholar] [CrossRef] [PubMed]
  52. Jin, G.; Liu, Y.; Zhang, J.; Bian, Z.; Yao, S.; Fei, B.; Zhou, L.; Yin, Y.; Huang, Z. A Panel of Serum Exosomal MicroRNAs as Predictive Markers for Chemoresistance in Advanced Colorectal Cancer. Cancer Chemother. Pharmacol. 2019, 84, 315–325. [Google Scholar] [CrossRef] [PubMed]
  53. Zhang, Z.; Zhang, L.; Yu, G.; Sun, Z.; Wang, T.; Tian, X.; Duan, X.; Zhang, C. Exosomal MiR-1246 and MiR-155 as Predictive and Prognostic Biomarkers for Trastuzumab-Based Therapy Resistance in HER2-Positive Breast Cancer. Cancer Chemother. Pharmacol. 2020, 86, 761–772. [Google Scholar] [CrossRef]
  54. Bhagirath, D.; Yang, T.L.; Bucay, N.; Sekhon, K.; Majid, S.; Shahryari, V.; Dahiya, R.; Tanaka, Y.; Saini, S. MicroRNA-1246 Is an Exosomal Biomarker for Aggressive Prostate Cancer. Cancer Res. 2018, 78, 1833–1844. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Liu, H.T.; Fan, W.X. MiRNA-1246 Suppresses the Proliferation and Migration of Renal Cell Carcinoma through Targeting CXCR4. Eur. Rev. Med. Pharmacol. Sci. 2020, 24, 5979–5987. [Google Scholar] [CrossRef] [PubMed]
  56. Chen, J.; Yao, D.; Li, Y.; Chen, H.; He, C.; Ding, N.; Lu, Y.; Ou, T.; Zhao, S.; Li, L.; et al. Serum MicroRNA Expression Levels Can Predict Lymph Node Metastasis in Patients with Early-Stage Cervical Squamous Cell Carcinoma. Int. J. Mol. Med. 2013, 32, 557–567. [Google Scholar] [CrossRef] [Green Version]
  57. Yang, Y.; Xie, Y.J.; Xu, Q.; Chen, J.X.; Shan, N.C.; Zhang, Y. Down-Regulation of MIR-1246 in Cervical Cancer Tissues and Its Clinical Significance. Gynecol. Oncol. 2015, 138, 683–688. [Google Scholar] [CrossRef]
  58. Wang, S.; Zeng, Y.; Zhou, J.M.; Nie, S.L.; Peng, Q.; Gong, J.; Huo, J.R. MicroRNA-1246 Promotes Growth and Metastasis of Colorectal Cancer Cells Involving CCNG2 Reduction. Mol. Med. Rep. 2016, 13, 273–280. [Google Scholar] [CrossRef] [Green Version]
  59. Hasegawa, S.; Eguchi, H.; Nagano, H.; Konno, M.; Tomimaru, Y.; Wada, H.; Hama, N.; Kawamoto, K.; Kobayashi, S.; Nishida, N.; et al. MicroRNA-1246 Expression Associated with CCNG2-Mediated Chemoresistance and Stemness in Pancreatic Cancer. Br. J. Cancer 2014, 111, 1572–1580. [Google Scholar] [CrossRef] [Green Version]
  60. Lin, S.S.; Peng, C.Y.; Liao, Y.W.; Chou, M.Y.; Hsieh, P.L.; Yu, C.C. MiR-1246 Targets CCNG2 to Enhance Cancer Stemness and Chemoresistance in Oral Carcinomas. Cancers 2018, 10, 272. [Google Scholar] [CrossRef] [Green Version]
  61. Li, X.J.; Ren, Z.J.; Tang, J.H.; Yu, Q. Exosomal MicroRNA MiR-1246 Promotes Cell Proliferation, Invasion and Drug Resistance by Targeting CCNG2 in Breast Cancer. Cell. Physiol. Biochem. 2017, 44, 1741–1748. [Google Scholar] [CrossRef] [PubMed]
  62. Si, G.; Li, S.; Zheng, Q.; Zhu, S.; Zhou, C. MiR-1246 Shuttling from Fibroblasts Promotes Colorectal Cancer Cell Migratio. Neoplasma 2021, 68, 317–324. [Google Scholar] [CrossRef] [PubMed]
  63. Cooks, T.; Pateras, I.S.; Jenkins, L.M.; Patel, K.M.; Robles, A.I.; Morris, J.; Forshew, T.; Appella, E.; Gorgoulis, V.G.; Harris, C.C. Mutant P53 Cancers Reprogram Macrophages to Tumor Supporting Macrophages via Exosomal MiR-1246. Nat. Commun. 2018, 9, 771. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Chen, J.; Yao, D.; Zhao, S.; He, C.; Ding, N.; Li, L.; Long, F. MiR-1246 Promotes SiHa Cervical Cancer Cell Proliferation, Invasion, and Migration through Suppression of Its Target Gene Thrombospondin 2. Arch. Gynecol. Obstet. 2014, 290, 725–732. [Google Scholar] [CrossRef] [PubMed]
  65. Yang, H.; Cho, M.E.; Li, T.W.H.; Peng, H.; Ko, K.S.; Mato, J.M.; Lu, S.C. MicroRNAs Regulate Methionine Adenosyltransferase 1A Expression in Hepatocellular Carcinoma. J. Clin. Investig. 2013, 123, 285–298. [Google Scholar] [CrossRef]
  66. Dai, D.; Tan, Y.; Guo, L.; Tang, A.; Zhao, Y. Identification of Exosomal MiRNA Biomarkers for Diagnosis of Papillary Thyroid Cancer by Small RNA Sequencing. Eur. J. Endocrinol. 2019, 182, 111–121. [Google Scholar] [CrossRef]
  67. Wang, H.; Wang, N.; Zheng, X.; Wu, L.; Fan, C.; Li, X.; Ye, K.; Han, S. Circular RNA Hsa_circ_0009172 Suppresses Gastric Cancer by Regulation of MicroRNA-485-3p-Mediated NTRK3. Cancer Gene Ther. 2021, 28, 1312–1324. [Google Scholar] [CrossRef]
  68. Formosa, A.; Markert, E.K.; Lena, A.M.; Italiano, D.; Finazzi-Agro’, E.; Levine, A.J.; Bernardini, S.; Garabadgiu, A.V.; Melino, G.; Candi, E. MicroRNAs, MiR-154, MiR-299-5p, MiR-376a, MiR-376c, MiR-377, MiR-381, MiR-487b, MiR-485-3p, MiR-495 and MiR-654-3p, Mapped to the 14q32.31 Locus, Regulate Proliferation, Apoptosis, Migration and Invasion in Metastatic Prostate Cancer Cells. Oncogene 2014, 33, 5173–5182. [Google Scholar] [CrossRef]
  69. Chen, D.; Fan, J.; Li, X.; Jiao, Z.; Tang, G.; Guo, X.; Chen, H.; Wang, J.; Men, T. Downregulation of MiR-485-3p Promotes Proliferation, Migration and Invasion in Prostate Cancer through Activation of TGF-β Signaling. Histol. Histopathol. 2022, 37, 423–430. [Google Scholar] [CrossRef]
  70. Zhang, Y.; Sui, R.; Chen, Y.; Liang, H.; Shi, J.; Piao, H. Downregulation of MiR-485-3p Promotes Glioblastoma Cell Proliferation and Migration via Targeting RNF135. Exp. Ther. Med. 2019, 18, 475–482. [Google Scholar] [CrossRef] [Green Version]
  71. Lou, C.; Xiao, M.; Cheng, S.; Lu, X.; Jia, S.; Ren, Y.; Li, Z. MiR-485-3p and MiR-485-5p Suppress Breast Cancer Cell Metastasis by Inhibiting PGC-1α Expression. Cell Death Dis. 2016, 7, e2159. [Google Scholar] [CrossRef] [Green Version]
  72. Xie, C.; Li, J.; Xu, C.; Xiong, W.; Yuan, X. CircRNA DNA Methyltransferase 1 Silence Inhibits Breast Cancer Development by Regulating MicoRNA-485-3p/Zinc Finger E-Box Binding Homeobox 1 Axis. J. Obstet. Gynaecol. Res. 2021, 47, 1068–1081. [Google Scholar] [CrossRef] [PubMed]
  73. Taherdangkoo, K.; Kazemi, N.S.R.; Hajjari, M.R.; Tahmasebi, B.M. MiR-485-3p Suppresses Colorectal Cancer via Targeting TPX2. Bratisl. Lek. Listy 2020, 121, 302–307. [Google Scholar] [CrossRef] [PubMed]
  74. Su, S.; Lu, W.; Liu, J.; Li, L.; Liu, L.; Li, X.; Ye, D.; Yu, Z. Circ_0007031 Silencing Inhibits Cell Proliferation and Induces Cell Apoptosis via Downregulating MELK at a MiR-485-3p-Dependent Way in Colorectal Cancer. Biochem. Genet. 2022, 60, 576–597. [Google Scholar] [CrossRef] [PubMed]
  75. Wang, Q.; Liu, M.J.; Bu, J.; Deng, J.L.; Jiang, B.Y.; Jiang, L.D.; He, X.J. MiR-485-3p Regulated by MALAT1 Inhibits Osteosarcoma Glycolysis and Metastasis by Directly Suppressing c-MET and AKT3/MTOR Signalling. Life Sci. 2021, 268, 118925. [Google Scholar] [CrossRef]
  76. Lai, J.; Xin, J.; Fu, C.; Zhang, W. CircHIPK3 Promotes Proliferation and Metastasis and Inhibits Apoptosis of Renal Cancer Cells by Inhibiting MiR-485-3p. Cancer Cell Int. 2020, 20, 248. [Google Scholar] [CrossRef] [PubMed]
  77. Wang, Z.Q.; Zhang, M.Y.; Deng, M.L.; Weng, N.Q.; Wang, H.Y.; Wu, S.X. Low Serum Level of MiR-485-3p Predicts Poor Survival in Patients with Glioblastoma. PLoS ONE 2017, 12, e0184969. [Google Scholar] [CrossRef]
  78. Ebrahimkhani, S.; Vafaee, F.; Hallal, S.; Wei, H.; Lee, M.Y.T.; Young, P.E.; Satgunaseelan, L.; Beadnall, H.; Barnett, M.H.; Shivalingam, B.; et al. Deep Sequencing of Circulating Exosomal MicroRNA Allows Non-Invasive Glioblastoma Diagnosis. NPJ Precis. Oncol. 2018, 2, 28. [Google Scholar] [CrossRef] [Green Version]
  79. Wang, M.; Cai, W.R.; Meng, R.; Chi, J.R.; Li, Y.R.; Chen, A.X.; Yu, Y.; Cao, X.C. MiR-485-5p Suppresses Breast Cancer Progression and Chemosensitivity by Targeting Survivin. Biochem. Biophys. Res. Commun. 2018, 501, 48–54. [Google Scholar] [CrossRef]
  80. Wang, J.; Zhang, Y.; Song, H.; Yin, H.; Jiang, T.; Xu, Y.; Liu, L.; Wang, H.; Gao, H.; Wang, R.; et al. The Circular RNA CircSPARC Enhances the Migration and Proliferation of Colorectal Cancer by Regulating the JAK/STAT Pathway. Mol. Cancer 2021, 20, 81. [Google Scholar] [CrossRef]
  81. Xiong, D.; Sheng, Y.; Ding, S.; Chen, J.; Tan, X.; Zeng, T.; Qin, D.; Zhu, L.; Huang, A.; Tang, H. LINC00052 Regulates the Expression of NTRK3 by MiR-128 and MiR-485-3p to Strengthen HCC Cells Invasion and Migration. Oncotarget 2016, 7, 47593–47608. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  82. Armand-Labit, V.; Meyer, N.; Casanova, A.; Bonnabau, H.; Platzer, V.; Tournier, E.; Sansas, B.; Verdun, S.; Thouvenot, B.; Hilselberger, B.; et al. Identification of a Circulating MicroRNA Profile as a Biomarker of Metastatic Cutaneous Melanoma. Acta Derm. Venereol. 2016, 96, 29–34. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Song, Y.; Washington, M.K.; Crawford, H.C. Loss of FOXA1/2 Is Essential for the Epithelial-to-Mesenchymal Transition in Pancreatic Cancer. Cancer Res. 2010, 70, 2115–2125. [Google Scholar] [CrossRef] [Green Version]
  84. Tang, Y.; Shu, G.; Yuan, X.; Jing, N.; Song, J. FOXA2 Functions as a Suppressor of Tumor Metastasis by Inhibition of Epithelial-to-Mesenchymal Transition in Human Lung Cancers. Cell Res. 2011, 21, 316–326. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  85. Zhang, Z.; Yang, C.; Gao, W.; Chen, T.; Qian, T.; Hu, J.; Tan, Y. FOXA2 Attenuates the Epithelial to Mesenchymal Transition by Regulating the Transcription of E-Cadherin and ZEB2 in Human Breast Cancer. Cancer Lett. 2015, 361, 240–250. [Google Scholar] [CrossRef]
  86. Lin, J.; Zhang, D.; Fan, Y.; Chao, Y.; Chang, J.; Li, N.; Han, L.; Han, C. Regulation of Cancer Stem Cell Self-Renewal by HOXB9 Antagonizes Endoplasmic Reticulum Stress-Induced Melanoma Cell Apoptosis via the MiR-765-FOXA2 Axis. J. Investig. Dermatol. 2018, 138, 1609–1619. [Google Scholar] [CrossRef] [Green Version]
  87. Basu, A.; Rojas, H.; Banerjee, H.; Cabrera, I.B.; Perez, K.Y.; de León, M.; Casiano, C.A. Expression of the Stress Response Oncoprotein LEDGF/P75 in Human Cancer: A Study of 21 Tumor Types. PLoS ONE 2012, 7, e30132. [Google Scholar] [CrossRef] [Green Version]
  88. Singh, D.K.; Gholamalamdari, O.; Jadaliha, M.; Li, X.L.; Lin, Y.C.; Zhang, Y.; Guang, S.; Hashemikhabir, S.; Tiwari, S.; Zhu, Y.J.; et al. PSIP1/P75 Promotes Tumorigenicity in Breast Cancer Cells by Promoting the Transcription of Cell Cycle Genes. Carcinogenesis 2017, 38, 966–975. [Google Scholar] [CrossRef] [Green Version]
  89. Ríos-Colón, L.; Cajigas-Du Ross, C.K.; Basu, A.; Elix, C.; Alicea-Polanco, I.; Sanchez, T.W.; Radhakrishnan, V.; Chen, C.-S.; Casiano, C.A. Targeting the Stress Oncoprotein LEDGF/P75 to Sensitize Chemoresistant Prostate Cancer Cells to Taxanes. Oncotarget 2017, 8, 24915–24931. [Google Scholar] [CrossRef] [Green Version]
  90. Ortiz-Hernandez, G.L.; Sanchez-Hernandez, E.S.; Ochoa, P.T.; Elix, C.C.; Alkashgari, H.R.; McMullen, J.R.W.; Soto, U.; Martinez, S.R.; Diaz Osterman, C.J.; Mahler, M.; et al. The LEDGF/P75 Integrase Binding Domain Interactome Contributes to the Survival, Clonogenicity, and Tumorsphere Formation of Docetaxel-Resistant Prostate Cancer Cells. Cells 2021, 10, 2723. [Google Scholar] [CrossRef]
  91. Gallo, M.; Sapio, L.; Spina, A.; Naviglio, D.; Calogero, A.; Naviglio, S. Lactic Dehydrogenase and Cancer: An Overview. Front. Biosci.-Landmark 2015, 20, 1234–1249. [Google Scholar] [CrossRef] [Green Version]
  92. Van Wilpe, S.; Koornstra, R.; Den Brok, M.; De Groot, J.W.; Blank, C.; De Vries, J.; Gerritsen, W.; Mehra, N. Lactate Dehydrogenase: A Marker of Diminished Antitumor Immunity. Oncoimmunology 2020, 9, 1731942. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  93. Petrelli, F.; Ardito, R.; Merelli, B.; Lonati, V.; Cabiddu, M.; Seghezzi, S.; Barni, S.; Ghidini, A. Prognostic and Predictive Role of Elevated Lactate Dehydrogenase in Patients with Melanoma Treated with Immunotherapy and BRAF Inhibitors: A Systematic Review and Meta-Analysis. Melanoma Res. 2019, 29, 1–12. [Google Scholar] [CrossRef] [PubMed]
  94. Gassenmaier, M.; Lenders, M.M.; Forschner, A.; Leiter, U.; Weide, B.; Garbe, C.; Eigentler, T.K.; Wagner, N.B. Serum S100B and LDH at Baseline and During Therapy Predict the Outcome of Metastatic Melanoma Patients Treated with BRAF Inhibitors. Target. Oncol. 2021, 16, 197–205. [Google Scholar] [CrossRef]
  95. Long, G.V.; Grob, J.J.; Nathan, P.; Ribas, A.; Robert, C.; Schadendorf, D.; Lane, S.R.; Mak, C.; Legenne, P.; Flaherty, K.T.; et al. Factors Predictive of Response, Disease Progression, and Overall Survival after Dabrafenib and Trametinib Combination Treatment: A Pooled Analysis of Individual Patient Data from Randomised Trials. Lancet Oncol. 2016, 17, 1743–1754. [Google Scholar] [CrossRef]
  96. Knispel, S.; Gassenmaier, M.; Menzies, A.M.; Loquai, C.; Johnson, D.B.; Franklin, C.; Gutzmer, R.; Hassel, J.C.; Weishaupt, C.; Eigentler, T.; et al. Outcome of Melanoma Patients with Elevated LDH Treated with First-Line Targeted Therapy or PD-1-Based Immune Checkpoint Inhibition. Eur. J. Cancer 2021, 148, 61–75. [Google Scholar] [CrossRef]
  97. Mazières, J.; Catherinne, C.; Delfour, O.; Gouin, S.; Rouquette, I.; Delisle, M.B.; Prévot, G.; Escamilla, R.; Didier, A.; Persing, D.H.; et al. Alternative processing of the U2 small nuclear RNA produces a 19-22nt fragment with relevance for the detection of non-small cell lung cancer in human serum. PLoS ONE 2013, 8, e60134. [Google Scholar] [CrossRef] [Green Version]
  98. Baraniskin, A.; Nöpel-Dünnebacke, S.; Ahrens, M.; Jensen, S.G.; Zöllner, H.; Maghnouj, A.; Wos, A.; Mayerle, J.; Munding, J.; Kost, D.; et al. Circulating U2 small nuclear RNA fragments as a novel diagnostic biomarker for pancreatic and colorectal adenocarcinoma. Int. J. Cancer 2013, 132, E48–E57. [Google Scholar] [CrossRef] [Green Version]
  99. Xu, Y.F.; Hannafon, B.N.; Khatri, U.; Gin, A.; Ding, W.Q. The origin of exosomal miR-1246 in human cancer cells. RNA Biol. 2019, 6, 770–784. [Google Scholar] [CrossRef]
Figure 1. Box-and-whisker diagrams of miR-1246, mir-92b-3p and miR-485-3p levels in melanoma patients treated with targeted therapy. (A). cf-miRNA levels were measured by qRT-PCR in plasma samples obtained from 55 patients (45 Rs and 10 NRs) before the start of therapy (T0). (B). cf-miRNA levels were measured by qRT-PCR in matched plasma samples obtained from 34 patients at T0 and at progression (TP). The edges of each box represent the 75th and 25th percentile, respectively, and whiskers are defined according to Tukey method. The horizontal bar within each box indicates the median. The outliers are reported as dots. Relative expression levels of miR-485-3p and miR-92b-3p are multiplied by 100. Data were analyzed by nonparametric Mann–Whitney U test to compare differences between groups, (A) and by the Wilcoxon matched-pairs signed-rank test for the before–after differences (B). ** p < 0.01 and * p < 0.05.
Figure 1. Box-and-whisker diagrams of miR-1246, mir-92b-3p and miR-485-3p levels in melanoma patients treated with targeted therapy. (A). cf-miRNA levels were measured by qRT-PCR in plasma samples obtained from 55 patients (45 Rs and 10 NRs) before the start of therapy (T0). (B). cf-miRNA levels were measured by qRT-PCR in matched plasma samples obtained from 34 patients at T0 and at progression (TP). The edges of each box represent the 75th and 25th percentile, respectively, and whiskers are defined according to Tukey method. The horizontal bar within each box indicates the median. The outliers are reported as dots. Relative expression levels of miR-485-3p and miR-92b-3p are multiplied by 100. Data were analyzed by nonparametric Mann–Whitney U test to compare differences between groups, (A) and by the Wilcoxon matched-pairs signed-rank test for the before–after differences (B). ** p < 0.01 and * p < 0.05.
Cancers 14 03706 g001
Figure 2. Kaplan–Meier curves for progression-free survival (PFS) of melanoma patients treated with targeted therapy. Patients were stratified by T0 levels of miR-1246 (A), miR-485-3p (B), both miR-1246 and miR-485-3p (C), miR-1246/miR485-3p ratio (D). Group (a): miR-1246 < 8.637 and miR-485-3p ≥ 0.013; group (b): miR-1246 ≥ 8.637 and miR-485-3p ≥ 0.013, or miR-1246 < 8.637 and miR-485-3p < 0.013; group (c): miR-1246 ≥ 8.637 and miR-485-3p < 0.013. p-values are from log-rank test (A,B,D) and log-rank test for trend (C). The estimated percentages of PFS at 12 months from the start of therapy is also reported for each comparison.
Figure 2. Kaplan–Meier curves for progression-free survival (PFS) of melanoma patients treated with targeted therapy. Patients were stratified by T0 levels of miR-1246 (A), miR-485-3p (B), both miR-1246 and miR-485-3p (C), miR-1246/miR485-3p ratio (D). Group (a): miR-1246 < 8.637 and miR-485-3p ≥ 0.013; group (b): miR-1246 ≥ 8.637 and miR-485-3p ≥ 0.013, or miR-1246 < 8.637 and miR-485-3p < 0.013; group (c): miR-1246 ≥ 8.637 and miR-485-3p < 0.013. p-values are from log-rank test (A,B,D) and log-rank test for trend (C). The estimated percentages of PFS at 12 months from the start of therapy is also reported for each comparison.
Cancers 14 03706 g002
Figure 3. Kaplan–Meier curves for overall survival (OS) of melanoma patients treated with targeted therapy. Patients were stratified by T0 levels of miR-1246 (A), miR-485-3p (B), both miR-1246 and miR-485-3p (C), miR-1246/miR-485-3p ratio (D). Group (a): miR-1246 < 8.637 and miR-485-3p ≥ 0.013; group (b + c): miR-1246 ≥ 8.637 and miR-485-3p ≥ 0.013, or miR-1246 < 8.637 and miR-485-3p < 0.013, or miR-1246 ≥ 8.637 and miR-485-3p < 0.013. p-values are from log-rank test. The estimated percentages of OS at 12 months from the start of therapy are also reported for each comparison.
Figure 3. Kaplan–Meier curves for overall survival (OS) of melanoma patients treated with targeted therapy. Patients were stratified by T0 levels of miR-1246 (A), miR-485-3p (B), both miR-1246 and miR-485-3p (C), miR-1246/miR-485-3p ratio (D). Group (a): miR-1246 < 8.637 and miR-485-3p ≥ 0.013; group (b + c): miR-1246 ≥ 8.637 and miR-485-3p ≥ 0.013, or miR-1246 < 8.637 and miR-485-3p < 0.013, or miR-1246 ≥ 8.637 and miR-485-3p < 0.013. p-values are from log-rank test. The estimated percentages of OS at 12 months from the start of therapy are also reported for each comparison.
Cancers 14 03706 g003
Table 1. Baseline demographic and clinical characteristics of melanoma patients included in the study.
Table 1. Baseline demographic and clinical characteristics of melanoma patients included in the study.
Number of patients57
Unresectable Stage IIIc a5
Stage IV a52
  M1a8
  M1b9
  M1c35
Male37
Female20
Age (years, range)20–82
sLDH
  Normal30
  High b27
ECOG PS c
  030
  116
  211
Previous Therapy
  Yes 2
  No 55
Targeted Therapy d
  DAB4
  VEM7
  DAB→DAB + TRAM9
  DAB + TRAM30
  VEM + COBI7
Best Response e
  CR 12
  PR 35
  SD 3
  PD 7
Progression on therapy f
  Yes52
  No 5
Dead for melanoma f42
a Disease stage before the beginning of targeted therapy according to 7th Edition of the AJCC Cancer Staging Manualb >1.5× upper limit of normal values. c Eastern Cooperative Oncology Group (ECOG) Performance Status. d DAB, dabrafenib, TRAM, trametinib, VEM, vemurafenib, COBI, cobimetinib. e CR, complete response; PR, partial response; SD, stable disease; PD, progressive disease, according to RECIST 1.1 criteria. f N. of patients who had experienced disease progression at the time of data analysis.
Table 2. cf-miRNAs DE in T0 plasma samples of NRs vs. Rs according to small RNA-seq data analysis.
Table 2. cf-miRNAs DE in T0 plasma samples of NRs vs. Rs according to small RNA-seq data analysis.
miRNAp (Adj) apbNRs vs. RsMean Value Rs cMean Value NRs cFC c
miR-12460.0400.001up9.46911.7142.245
miR-125b-2-3p0.1040.021up3.2195.0551.836
miR-10b-5p0.1450.042up13.71015.0331.323
miR-378a-3p0.1000.018up8.8639.5440.681
miR-92b-3p0.0690.005up11.00311.6480.644
miR-6837-3p0.0210.0002down3.5830.976−2.607
miR-485-3p0.0720.007down4.7642.901−1.863
miR-224-5p0.1040.024down7.4245.615−1.810
miR-369-3p0.0880.009down4.8643.060−1.804
miR-766-3p0.1450.047down5.1893.419−1.771
miR-181c-3p0.0690.004down6.2044.526−1.678
miR-7-1-3p0.1000.016down3.7522.179−1.573
miR-379-3p0.1450.049down3.1991.673−1.527
miR-15b-3p0.1300.036down3.7562.270−1.486
miR-139-5p0.1250.032down5.8994.427−1.472
miR-487a-3p0.1000.018down3.5412.076−1.465
miR-42860.1040.023down3.8402.380−1.460
miR-136-3p0.1450.047down6.3054.883−1.421
miR-1260b0.1450.047down4.2792.870−1.410
miR-3605-5p0.1000.015down4.3342.935−1.399
miR-44350.1000.015down5.2403.901−1.339
miR-194-5p0.1250.031down5.0583.786−1.273
miR-4662a-5p0.1040.022down5.3244.102−1.223
miR-590-3p0.1250.033down4.8923.701−1.191
miR-342-3p0.1000.018down8.5937.719−0.873
miR-20a-5p0.0690.004down7.6346.857−0.777
miR-15b-5p0.0420.001down8.6587.883−0.775
miR-628-3p0.0720.006down7.3346.600−0.734
miR-223-3p0.1000.017down11.43010.757−0.674
ap (Adj), Benjamini–Hochberg adjusted p-value. b p, p-value. c Data are in log2 format. FC, Fold change.
Table 3. cf-miRNAs DE in TP vs. T0 plasma samples of Rs according to small RNA-seq data analysis.
Table 3. cf-miRNAs DE in TP vs. T0 plasma samples of Rs according to small RNA-seq data analysis.
miRNAp (Adj) apbTP vs. T0Mean Value T0 cMean Value TP cFC c
miR-12460.0860.021up9.46910.4270.958
miR-548b-5p0.0860.014up2.8273.6110.784
miR-144-5p0.0860.025up8.7929.5020.709
miR-17-3p0.0860.032up4.0034.5890.586
miR-42860.0030.000down3.8402.321−1.519
miR-485-3p0.0560.004down4.7643.461−1.303
miR-369-3p0.0860.029down4.8643.855−1.009
miR-99b-3p0.0860.026down3.5802.687−0.892
miR-664a-5p0.0860.020down4.3253.472−0.853
miR-487a-3p0.1050.066down3.5412.773−0.768
miR-5430.0860.038down4.5053.813−0.692
miR-128-1-5p0.0860.035down3.1682.546−0.622
miR-199b-5p0.0860.033down3.4882.885−0.603
ap (Adj), Benjamini–Hochberg adjusted p-value. b p, p-value. c Data are in log2 format. FC, Fold change.
Table 4. ROC analysis on cf-miRNAs DE in T0 plasma samples of NRs vs. Rs according to small RNA-seq data analysis.
Table 4. ROC analysis on cf-miRNAs DE in T0 plasma samples of NRs vs. Rs according to small RNA-seq data analysis.
miRNANRs vs. RsAUC (95% CI) aCut-Off bSensitivity (%)Specificity (%)
miR-1246up0.92 (0.79–0.99)10.2100.076.9
miR-92b-3pup0.88 (0.70–0.96)11.4100.069.2
miR-125b-2-3pup0.83 (0.66–0.95)4.680.080.1
miR-378a-3pup0.82 (0.66–0.95)8.9100.061.5
miR-10b-5pup0.72 (0.52–0.86)13.5100.042.3
miR-6837-3pdown0.90 (0.74–0.98)2.992.380.0
miR-15b-5pdown0.88 (0.70–0.96)8.569.2100.00
miR-4286down0.86 (0.70–0.96)3.288.580.0
miR-342-3pdown0.85 (0.66–0.95)8.469.2100.0
miR-487a-3pdown0.83 (0.66–0.95)3.365.4100.0
miR-20a-5pdown0.82 (0.66–0.95)7.557.7100.0
miR-223-3pdown0.79 (0.63–0.93)11.173.180.0
miR-485-3pdown0.79 (0.63–0.93)4.280.880.0
miR-628-3pdown0.79 (0.63–0.93)6.984.680.0
miR-379-3pdown0.78 (0.59–0.90)3.073.180.0
miR-1260bdown0.77 (0.59–0.90)4.057.7100.0
miR-181c-3pdown0.76 (0.59–0.90)6.069.280.0
miR-766-3pdown0.74 (0.55–0.88)4.876.980.0
miR-369-3pdown0.73 (0.55–0.88)3.596.260.0
miR-4662a-5pdown0.73 (0.55–0.88)5.553.9100.0
miR-7-1-3pdown0.73 (0.55–0.88)1.996.260.0
miR-139-5pdown0.72 (0.52–0.86)5.857.7100.0
miR-590-3pdown0.70 (0.52–0.86)5.238.5100.0
miR-224-5pdown0.69 (0.49–0.83)7.461.580.0
miR-15b-3pdown0.68 (0.49–0.83)3.946.2100.0
miR-3605-5pdown0.68 (0.49–0.83)4.261.580.0
miR-4435down0.68 (0.49–0.83)4.0100.060.0
miR-136-3pdown0.65 (0.45–0.81)6.357.780.0
miR-194-5pdown0.59 (0.39–0.75)3.9100.040.0
a AUC, area under curve; CI, confidence interval. b Each value is in log2 format.
Table 5. T0 plasma level of cf-miRNAs according to PFS categories.
Table 5. T0 plasma level of cf-miRNAs according to PFS categories.
PFS (Days) a
miRNAAll≤180181–365≥365ptrendb
miR-92b-3p11.1 (10.6–11.6)11.5 (11.4–11.6)11.2 (10.7–11.5)10.6 (10.5–10.7)<0.0001
miR-485-3p4.8 (3.9–5.4)4.0 (3.3–5.0)4.7 (4.2–5.3)5.6 (5.3–5.8)0.001
miR-15b-5p8.5 (8.3–8.9)8.4 (8.2–8.6)8.5 (8.2–8.9)8.8 (8.6–9.0)0.015
miR-181c-3p6.0 (5.6–6.7)5.9 (5.6–6.3)6.0 (5.4–6.9)6.4 (6.1–6.9)0.021
miR-487a-3p3.4 (2.7–4.2)3.2 (2.4–3.5)3.5 (2.5–4.1)4.6 (3.2–5.2)0.028
miR-379-3p3.2 (1.5–4.1)3.0 (1.4–3.4)3.3 (3.0–4.1)4.2 (2.6–4.8)0.028
miR-378a-3p8.8 (8.4–9.5)9.0 (8.7–9.6)9.1 (8.4–9.5)8.4 (8.2–8.8)0.031
miR-4662a-5p5.4 (4.9–5.6)5.4 (5.0–5.5)5.1 (4.7–5.6)5.8 (5.5–5.9)0.045
miR-12469.7 (8.7–10.6)10.2 (9.0–11.7)10.0 (9.2–10.6)8.9 (8.5–9.3)0.056
miR-628-3p7.3 (6.8–7.6)7.3 (6.8–7.4)7.0 (6.8–7.1)7.6 (7.4–7.8)0.059
miR-7-1-3p3.6 (3.1–4.5)3.3 (1.7–3.6)4.3 (3.4–4.8)3.9 (3.4–4.3)0.059
miR-224-5p7.6 (6.3–8.2)6.8 (5.9–7.7)8.1 (6.3–8.2)7.7 (7.2–8.5)0.063
miR-139-5p5.7 (5.4–6.5)5.6 (5.3–5.8)6.1 (4.9–7.1)6.2 (5.5–7.0)0.087
miR-342-3p8.5 (7.7–9.1)8.3 (7.6–9.1)8.5 (7.7–9.0)8.9 (8.4–9.3)0.087
miR-136-3p6.3 (5.5–7.0)6.1 (5.2–6.7)6.1 (5.2–6.8)6.9 (6.3–7.4)0.091
miR-42863.9 (3.2–4.5)3.5 (2.8–4.1)3.9 (3.5–4.4)4.3 (3.4–4.8)0.114
miR-223-3p11.3 (10.9–11.8)11.2 (10.7–11.4)11.3 (10.9–11.4)11.7 (11.2–11.9)0.123
miR-20a-5p7.5 (7.3–7.8)7.5 (7.1–7.8)7.4 (7.3–7.7)7.8 (7.5–8.1)0.178
miR-10b-5p13.8 (12.9–15.1)13.9 (13.3–14.4)14.6 (13.7–15.5)13.3 (12.5–13.8)0.203
miR-766-3p5.2 (4.8–6.0)4.9 (4.5–6.2)5.2 (4.8–5.7)5.8 (5.1–6.3)0.216
miR-125b-2-3p3.9 (2.0–4.8)4.2 (3.6–5.2)3.0 (1.8–4.5)3.5 (2.4–4.0)0.255
miR-194-5p5.0 (4.4–5.4)5.0 (4.8–5.4)5.0 (4.4–5.8)4.7 (4.3–5.2)0.282
miR-369-3p4.9 (4.2–5.4)5.1 (3.5–5.2)4.6 (4.4–5.0)5.3 (4.3–6.3)0.293
miR-44355.2 (4.6–5.7)5.1 (4.3–5.7)5.1 (5.0–5.5)5.3 (4.9–5.7)0.340
miR-1260b4.0 (3.6–5.4)3.9 (3.4–5.1)4.1 (3.4–5.4)4.3 (3.8–5.3)0.385
miR-15b-3p3.6 (2.9–4.4)3.3 (2.9–3.7)4.2 (3.4–4.8)3.3 (2.8–4.2)0.405
miR-590-3p4.8 (4.4–5.4)4.8 (4.6–5.3)4.6 (3.8–4.8)5.5 (4.8–5.9)0.434
miR-6837-3p3.5 (3.0–4.2)3.6 (1.5–4.6)3.5 (3.1–4.6)3.5 (3.1–3.6)0.581
miR-3605-5p4.2 (3.8–4.8)4.2 (2.9–4.5)4.9 (4.6–5.1)3.8 (3.6–4.2)0.807
a PFS, progression free survival. Values are expressed in terms of median. In parentheses: interquartile range. All values are in log2 format. b Cuzick’s test for trend across the three categories.
Table 6. AUC for T0 plasma levels of miR-1246, miR-485-3p and their ratio according to qRT-PCR data.
Table 6. AUC for T0 plasma levels of miR-1246, miR-485-3p and their ratio according to qRT-PCR data.
miRNAAUC (95%CI) aCut-OffSensitivity (%)Specificity (%)
miR-12460.85 (0.73–0.94)8.63770.0%88.9%
miR-485-3p0.71 (0.57–0.82)0.01382.2%70.0%
miR-1246/miR-485-3p ratio0.84 (0.71–0.92)345.49390.0%82.2%
a AUC, area under curve; CI, confidence interval.
Table 7. Univariate Cox regression analysis for disease progression according to T0 plasma levels of miR-1246, miR-485-3p and clinical features of patients.
Table 7. Univariate Cox regression analysis for disease progression according to T0 plasma levels of miR-1246, miR-485-3p and clinical features of patients.
CharacteristicsEvents a/PatientsHR b (95% CI b)pb
miR-1246
 <8.63738/431
 ≥8.63712/122.72 (1.37–5.39)0.004
miR-485-3p
 ≥0.01334/391
 <0.01316/162.86 (1.51–5.43)0.001
miR-1246/miR-485-3p ratio
 <345.49333/381
 ≥345.49317/173.26 (1.73–6.14)<0.0001
miR-1246 and miR-485-3p combinations c
 Group a 29/341
 Group b 14/142.53 (1.27–5.06)0.008
 Group c 7/74.92 (2.02–12.0)<0.0001
sLDH
 Normal25/301
 High d25/254.80 (2.41–9.57)<0.0001
Stage e
 IIIC + M1a + M1b18/221
 M1c32/332.72 (1.48–5.02)0.001
Gender
 Female16/201
 Male34/351.90 (1.04–3.49)0.037
Age, years
 ≤6028/301
 >6022/250.77 (0.44–1.35)0.360
ECOG PS f
 025/301
 ≥125/253.61 (1.93–6.76)<0.0001
a Event: disease progression. b HR, hazard ratio; CI, confidence interval; p, probability. Estimated by Cox’s regression model. c Group (a): miR-1246 < 8.637 and miR-485-3p ≥ 0.013; group (b): miR-1246 ≥ 8.637 and miR-485-3p ≥ 0.013, or miR-1246 < 8.637 and miR-485-3p < 0.013; group (c): miR-1246 ≥ 8.637 and miR-485-3p < 0.013. d >1.5× upper limit of normal values. e Disease stage before the beginning of targeted therapy according to 7th Edition of the AJCC Cancer Staging Manual. f ECOG PS, Eastern Cooperative Oncology Group Performance Status.
Table 8. Univariate Cox regression analysis for mortality according to T0 plasma levels of miR-1246, miR-485-3p and clinical features of patients.
Table 8. Univariate Cox regression analysis for mortality according to T0 plasma levels of miR-1246, miR-485-3p and clinical features of patients.
CharacteristicsEvents a/PatientsHR b (95% CI b)pb
miR-1246
 <8.63719/271
 ≥8.63710/112.66 (1.19–5.96)0.017
miR-485-3p
 ≥0.01318/261
 <0.01311/122.37 (1.08–5.17)0.031
miR-1246/miR-485-3p ratio
 <345.49317/251
 ≥345.49312/132.91 (1.34–6.30)0.007
miR-1246 and miR-485-3p combinations c
 Group a 14/221
 Group b + c15/162.98 (1.37–6.49)0.006
sLDH
 Normal11/191
 High d18/195.31 (2.18–13.0)<0.0001
Stage e
 IIIC + M1a + M1b5/111
 M1c24/274.77 (1.76–13.0)0.002
Gender
 Female8/151
 Male21/232.78 (1.22–6.37)0.015
Age, years
 ≤6017/201
 >6012/180.69 (0.33–1.45)0.324
ECOG PS f
 012/201
 ≥117/184.45 (1.90–10.5)0.001
a event: death. b HR, hazard ratio; CI, confidence interval; p, probability. Estimated by Cox’s regression model. c Group (a): miR-1246 < 8.637 and miR-485-3p ≥ 0.013; group (b): miR-1246 ≥ 8.637 and miR-485-3p ≥ 0.013, or miR-1246 < 8.637 and miR-485-3p < 0.013; group (c): miR-1246 ≥ 8.637 and miR-485-3p < 0.013. d >1.5× upper limit of normal values. e Disease stage before the beginning of targeted therapy according to 7th Edition of the AJCC Cancer Staging Manualf ECOG PS, Eastern Cooperative Oncology Group Performance Status.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Levati, L.; Bassi, C.; Mastroeni, S.; Lupini, L.; Antonini Cappellini, G.C.; Bonmassar, L.; Alvino, E.; Caporali, S.; Lacal, P.M.; Narducci, M.G.; et al. Circulating miR-1246 and miR-485-3p as Promising Biomarkers of Clinical Response and Outcome in Melanoma Patients Treated with Targeted Therapy. Cancers 2022, 14, 3706. https://doi.org/10.3390/cancers14153706

AMA Style

Levati L, Bassi C, Mastroeni S, Lupini L, Antonini Cappellini GC, Bonmassar L, Alvino E, Caporali S, Lacal PM, Narducci MG, et al. Circulating miR-1246 and miR-485-3p as Promising Biomarkers of Clinical Response and Outcome in Melanoma Patients Treated with Targeted Therapy. Cancers. 2022; 14(15):3706. https://doi.org/10.3390/cancers14153706

Chicago/Turabian Style

Levati, Lauretta, Cristian Bassi, Simona Mastroeni, Laura Lupini, Gian Carlo Antonini Cappellini, Laura Bonmassar, Ester Alvino, Simona Caporali, Pedro Miguel Lacal, Maria Grazia Narducci, and et al. 2022. "Circulating miR-1246 and miR-485-3p as Promising Biomarkers of Clinical Response and Outcome in Melanoma Patients Treated with Targeted Therapy" Cancers 14, no. 15: 3706. https://doi.org/10.3390/cancers14153706

APA Style

Levati, L., Bassi, C., Mastroeni, S., Lupini, L., Antonini Cappellini, G. C., Bonmassar, L., Alvino, E., Caporali, S., Lacal, P. M., Narducci, M. G., Molineris, I., De Galitiis, F., Negrini, M., Russo, G., & D’Atri, S. (2022). Circulating miR-1246 and miR-485-3p as Promising Biomarkers of Clinical Response and Outcome in Melanoma Patients Treated with Targeted Therapy. Cancers, 14(15), 3706. https://doi.org/10.3390/cancers14153706

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

Article Metrics

Back to TopTop