Skip to main content

miR-146a-5p impairs melanoma resistance to kinase inhibitors by targeting COX2 and regulating NFkB-mediated inflammatory mediators



Targeted therapy with BRAF and MEK inhibitors has improved the survival of patients with BRAF-mutated metastatic melanoma, but most patients relapse upon the onset of drug resistance induced by mechanisms including genetic and epigenetic events. Among the epigenetic alterations, microRNA perturbation is associated with the development of kinase inhibitor resistance. Here, we identified and studied the role of miR-146a-5p dysregulation in melanoma drug resistance.


The miR-146a-5p-regulated NFkB signaling network was identified in drug-resistant cell lines and melanoma tumor samples by expression profiling and knock-in and knock-out studies. A bioinformatic data analysis identified COX2 as a central gene regulated by miR-146a-5p and NFkB. The effects of miR-146a-5p/COX2 manipulation were studied in vitro in cell lines and with 3D cultures of treatment-resistant tumor explants from patients progressing during therapy.


miR-146a-5p expression was inversely correlated with drug sensitivity and COX2 expression and was reduced in BRAF and MEK inhibitor-resistant melanoma cells and tissues. Forced miR-146a-5p expression reduced COX2 activity and significantly increased drug sensitivity by hampering prosurvival NFkB signaling, leading to reduced proliferation and enhanced apoptosis. Similar effects were obtained by inhibiting COX2 by celecoxib, a clinically approved COX2 inhibitor.


Deregulation of the miR-146a-5p/COX2 axis occurs in the development of melanoma resistance to targeted drugs in melanoma patients. This finding reveals novel targets for more effective combination treatment.

Video Abstract

Graphical Abstract


Treatment with BRAF and MEK inhibitors (BRAF/MEKi) has significantly improved the outcome of metastatic cutaneous melanoma patients with BRAF-mutated tumors. However, in the majority of patients, the clinical benefits are limited by the emergence of drug resistance [1]. Although BRAF/MEKi-resistant patients can benefit from different therapeutic opportunities, such as immune checkpoint blockade, strategies aimed at impairing the onset of resistance by drug combinations remain an urgent unmet clinical need.

The complex network of genomic modifications driving tumor evolution under pressure from drug treatment leads to MAPK signaling reactivation and/or AKT pathway sustainment and strongly conditions the tumor microenvironment (TME) [2]. Despite treatment, epigenetic changes governing drug-induced melanoma cell plasticity and disease progression activate protumorigenic signaling pathways that include microRNAs (miR) as regulators [3]. miR act in a cell context-dependent manner by playing key regulatory roles in cell survival and proliferation signaling pathways, contributing to the acquisition of BRAF/MEKi resistance [4, 5]. Our previous work highlighted the link between the acquisition of an HIF1/NFkB-directed proinflammatory phenotype and the dysregulation of a set of miR involved in prosurvival signaling pathways in BRAFi-resistant melanoma cell lines [6]. The acquisition of an HIF1/NFkB-directed proinflammatory phenotype was characterized by the induced production and release of several factors, including TNFα, CCL5, IL6, VEGFA, CXCL8, PGE2 and CCL2, which had a prosurvival effect on melanoma cells [6].

miR-146a-5p (miR-146a) is transcriptionally regulated by NFkB. It negatively controls NFkB-associated signaling pathways and release of proinflammatory cytokines through downregulating IRAK1 and TRAF6, and by negatively regulating expression of TNFα [7,8,9,10,11]. In melanoma, antithetic roles for miR-146a as both a tumor suppressor and tumor promotor were described [12,13,14,15], and its genetic polymorphism was associated with melanoma risk [16]. In particular, miR-146a overexpression was reported to promote melanoma cell growth by downregulating its direct targets NUMB and LFNG, which are involved in the NOTCH/PTEN/AKT pathway [12, 13]. Other investigators reported miR-146a involvement in the regulation of melanoma cell migration and invasive ability via targeting SMAD4 [15], ITGAV and ROCK1 [13]. Nonetheless, in other tumor types, such as breast, pancreatic and gastric cancer, miR-146a was reported to suppress cell proliferation, invasion and metastasis by repressing EGFR expression [7]. In lung cancer, miR-146a overexpression was shown to restore cisplatin sensitivity by targeting cyclin J and ATG12, while in liver cancer, miR-146a expression increased the radiosensitivity of tumor cells by downregulating RPA3 expression [7]. Of note, by negatively regulating EGFR expression, miR-146a could enhance the antitumor effects of drugs targeting this receptor [7].

Among the proinflammatory factors involved in tumor progression, cyclooxygenase-2 (COX2) and its product prostaglandin E2 (PGE2) have been found to be upregulated in many cancer types, favoring cell growth and survival and contributing to the generation of an immune suppressive TME [17]. Prostaglandin production depends on the activity of two isoforms of the COX enzyme, COX1 and COX2. While COX1 is constitutively expressed, COX2 is induced by different factors, such as growth factors and proinflammatory cytokines, including CCL2. Celecoxib is a nonsteroidal anti-inflammatory drug that specifically inhibits COX2, commonly used for the treatment of rheumatoid arthritis [18]. This drug is currently being repurposed in a clinical trial evaluating its treatment and prevention of cancer [18].

Notably, low miR-146a levels have been found in many tumor types, including melanoma [7], and were associated with increased COX2 levels in lung cancer [19]. However, the implication of the miR-146a/COX2 axis and its potential inhibition by celecoxib has not yet been investigated in the context of BRAF/MEKi resistance. Here, we identified the miR-146a/COX2 axis as playing a critical role in the mechanisms sustaining melanoma resistance; thus, this axis represents a druggable signaling pathway for improving the effects of BRAF/MEKi therapy in melanoma patients.


Cell cultures and reagents

The melanoma cell lines used in this study were described previously [6] and were periodically checked for mycoplasma and authentication by STR analysis (Gene Print 10 System, Promega, Fitchburg, WI, USA). Cell Counting kit 8 (CCK8, Sigma) was used to determine cell proliferation, LDH cytotoxicity assay kit (Pierce, Waltham, MA, USA) for LDH release, and Caspase-Glo 8 and 3/7 assays for detecting activated caspases (Promega, Fitchburg, WI, USA). The Prostaglandin E2 ELISA Kit (Cayman Chemical, Ann Arbor, MI, USA) was used for the measurement of PGE2 from culture supernatant of melanoma cells and of 3D cultures. Vemurafenib (PLX4032), Dabrafenib (BRAFi) and Trametinib (MEKi) were purchased (Selleck, Houston, TX, USA) and used at 3 μM for PLX4032 and 3 nM for BRAF/MEKi; recombinant sTRAIL (Vinci Biochem, Florence, Italy) was used at 50 ng/mL and the COX2 inhibitors celecoxib and NS398 (Sigma) at 50 μM.

For cells transfection of miR mimic, inhibitor and control oligos (50 nM), PTGS2 siRNA or siRNA control (100 nM) (Thermo Fisher Scientific, Waltham, MA, USA), Metafectene was used (Biontex, Munich, Germany). Transfection efficiency was monitored by qRT-PCR for miR-146a and by western blot for COX2 protein expression at 72 h. For Luciferase Assay, pLightSwitch-3’UTR-PTGS2 or the pLightSwitch empty vector (1 μg) (Active Motif, Carlsbad, CA, USA) were cotrasfected with either miR-146a mimic or negative control for 48 h and the LightSwitch Luciferase assay system (Active Motif) was used to analyze the luminescence production. All experiments were performed at 72 h after treatment, seeding 1,5 × 104 cells/well in 96 wells plates.

Western blot analysis and antibody array

Total protein was extracted from transfected cells using RIPA buffer (Thermo Fisher Scientific) and quantified using BCA protein assay (Bio-Rad, Hercules, CA, USA). Cell lysates (30 μg) were resolved on 4–12% SDS-PAGE and transferred to nitrocellulose membranes (Thermo Fisher Scientific); membranes were incubated with specific antibodies (supplementary Table S1). Peroxidase-conjugated secondary antibodies anti-mouse immunoglobulin and anti-rabbit immunoglobulin G were used. Antibody Array kit for Human Apoptosis Signaling (Raybiotech, Peachtree Corners, GA, USA) was used to test protein lysates. Chemiluminescence was measured by Uvitec Imaging System (Cleaver Scientific, Cambridge, UK) and quantified by Nine-Alliance software.

RNA extraction and qRT-PCR analysis

RNA was extracted from melanoma cells and specimens with the NucleoSpin RNA isolation kit (Macherey Nagel, Bethlehem, PA, USA). Extracted RNAs were quantified spectrophotometrically, and the absorbance ratio at 260/280 and 260/230 were measured to assure RNA quality and purity. qRT-PCR analysis was performed using Thermo Fisher Scientific reagents for gene transcripts and Exiqon reagents for miR (supplementary Table S1). The endogenous controls used for normalization were ACTB and RPL13A for genes and U6 snRNA for miR. qRT-PCR was carried out in triplicate and run on the ABI Prism 7900HT or on the QuantStudio 7 Flex instruments and analysis was performed using SDS software, version 2.2.2 and with QuantStudio 6 and 7 Flex software. The results are presented as 2−ΔCt ± SD for direct comparisons.

miR and gene expression profiling data analysis

miR expression profiles were generated using NanoString nCounter Human v2 miR Expression Assay that contains 800 human endogenous miR. Raw data were normalized and log2-transformed using the NanoStringNorm R package [20], setting the required parameters to the following values: Probe.Correction.Factor to “adjust”; CodeCount to “sum”; Background to “mean”; SampleContent to “top.geo.mean” and OtherNorm to “none”. Raw and processed miR expression data are available at the NCBI Gene Expression Omnibus with accession date GSE141314. Raw miR expression data for GSE67635 [21] and GSE68841 [6] were downloaded from GEO and pre-processed using limma [22] with “normexp” method for background correction, quantile normalization and log2 transformation. Replicated probes were collapsed calculating the average expression and for each miR the probe with highest variance across samples was selected. Differential expression analysis was carried out using limma [22] (paired design for NanoString data). P-values were corrected for multiple testing using the Benjamini-Hochberg false discovery rate (FDR) method. A nominal p-value < 0.05 was applied to select differentially expressed miR for NanoString data and an FDR < 0.05 for GSE67635 and GSE68841.

Gene expression profiles of metastatic tumor samples were generated using HumanHT-12 WG-DASL V4.0 R2 expression beadchip (Illumina, San Diego, CA, USA). RNA labeling, processing and hybridization were performed according the manufacturer’s standard protocols. Microarrays were scanned with Illumina BeadArray Reader and raw expression data were obtained using Illumina BeadStudio v3.3.8 and processed using the lumi Bioconductor package [23] as previously described [6]. The data were deposited in GEO with accession number GSE141484. Functional analyses were generated through the use of Ingenuity pathway Analysis (IPA) (Qiagen, Hilden, Germany).

3D cultures of tumor explants in RCCS bioreactor

3D Cultures of Melanoma Tissue explants were set in Rotary Cell Culture System Bioreactor (RCCS) (Synthecon Inc., Houston TX, USA) as described by Ferrarini [24]. Tumor tissue cubes were obtained by a 3 mm biopsy puncher from tumor specimens obtained from the pathologist, and cultured in duplicate in the bioreactor chambers for 3 days, in culture medium or in the presence of BRAF/MEKi (3 nM) and/or miR-146a mimic or negative control synthethic oligos (50 nM). When recovered, samples were halved, a part was fixed in formalin and paraffin-embedded for histopathology, and a part was fixed in RNAlater solution (Thermo Fisher Scientific) for RNA analysis. Quantification of secreted cytokines in the 3D culture supernatant was carried by Cytometric Bead Array (CBA, BD Biosciences, San Jose, CA, USA) using the BD software FCAP Array v3.0.1.


Melanoma sections were stained by specific antibodies (supplementary Table S1) for ki67, cleaved Caspase 3, phospho-ERK, and COX2 after antigen retrieval performed by heating in a pressure cooker with 0,5 mM EDTA pH 8 for 15 or 20 min and using a peroxidase-labelled polymer (UltraVision Quanto Detection System HRP Polymer, Thermo Fisher Scientific) and brown DAB as a chromogen (Dako Agilent, Glostrup, Denmark). Sections were scanned using the Aperio ScanScope XT systems (Aperio Technologies, Leica Microsystems).

Statistical analysis

Statistical analyses were performed using Prism software v.5 and v.8 (La Jolla, CA, USA). Comparisons between continuous variables in two groups were performed using an unpaired two tailed Student’s t-test or Mann-Whitney U test. For comparisons involving more than two groups, one-way ANOVA was used, followed by Bonferroni correction. For comparisons between two dose-response curves, two-way ANOVA followed by Bonferroni correction was used. The correlation between linear variables was calculated using Pearson or Spearman’s correlation coefficients. Drug interaction was evaluated by a standard approach that allows a value to be assigned to a drug combination (interaction index); the interaction index was calculated as the ratio between expected cell growth and observed cell growth. If the experimentally measured effect of the drug combination was equal to, higher than or lower than 1, the combination was considered to be additive, synergistic or antagonistic, respectively. Data are presented as the mean ± SD.


miR-146a is regulated in BRAFi-resistant melanoma cell lines

To expand the analysis of the miR involved in establishing a resistant phenotype, we analyzed the miR expression profiles of six melanoma cell lines which developed acquired resistance to vemurafenib (PLX4032) and cross resistance to dabrafenib (BRAFi), trametinib (MEKi) and the combination of dabrafenib and trametinib (BRAF/MEKi) and their parental drug-sensitive counterparts (supplementary Fig. S1A). We identified 22 differentially expressed miR between the resistant and sensitive cells; of these miR, 6 were upregulated, and 16 were downregulated in the resistant cells (supplementary Table S2). To strengthen our findings and focus on miR which could be relevant for resistance, acquired and intrinsic, we intersected these 22 miR to those previously identified in our laboratory in a similar profiling setting (GSE68841; 6) and to those differentially expressed between melanoma cell lines sensitive (cluster 1) or intrinsically resistant to the BRAFi PLX4720 (cluster 2) (GSE67635, 21). Three miR were commonly modulated in all three datasets (Fig. 1a). Two of them, miR-100-5p and miR-125b-5p, upregulated in resistant cells, have been previously investigated [6]. The third miR, miR-146a, was downregulated in resistant cells and its down-regulation was indeed confirmed in each of the six BRAFi-resistant melanoma cell lines of the current study compared to their sensitive counterparts (Fig. 1b). First evidence of its association with drug resistance was given by the inverse correlation of miR-146a and the IC50 values of BRAF/MEKi used in combination or as single agents (Fig. 1c and supplementary Fig. S1B).

Fig. 1

Downregulation of miR-146a is associated with resistance to BRAF/MEKi in melanoma cell lines and its manipulation affects drug response and TRAIL-induced apoptosis. a Venn diagram illustrating common miR associated with BRAF/MEKi resistance in GSE141314, GSE67635 and GSE68841 datasets. b miR-146a shows significant downregulation in melanoma resistant cells (R) compared with their sensitive counterparts (S). logFC: log2 Fold Change. P-value obtained from differential expression analysis performed with limma package. c Inverse correlation between miR-146a expression levels and IC50 values of combined BRAF/MEKi in melanoma cell lines (Spearman analysis). d Forced expression of miR-146a (+m-miR-146a) increases the effects of BRAF/MEKi treatment in LM16R, LM69 and LM70 cell lines, as shown by reduced cell growth and increased cell cytotoxicity and apoptosis, evaluated by CCK8, LDH and caspase 8 and 3/7 activity. LM69 and LM70 short term melanoma cultures were generated from treatment resistant melanoma lesions surgically excised from two patients. Based on qRT-PCR results, upon transfection miR-146a levels were up to 70-fold higher (range 76–1223) than the levels detectable in cells transfected with control oligos. Data are plotted as fold increase compared to cells transfected with mimic negative control. *: p < 0.05, **: p < 0.01, ***: p < 0.0001 by Student’s unpaired t test. e Forced expression of miR-146a increased the effect of treatment with sTRAIL. f Positive correlation between the levels of miR-146a expression and of activated caspase 3/7 induced by sTRAIL treatment (Spearman analysis). RLU: Relative Light Unit. Results shown are representative of 2 experiments performed in triplicate (d, e)

Modulation of miR-146a affects the BRAF/MEKi drug response in melanoma cells

To assess whether miR-146a expression levels affect the drug sensitivity of melanoma cells, we transiently manipulated its expression in gain- and loss-of-function assays and tested the effects on kinase inhibitor treatment. miR-146a hyperexpression in cell lines with acquired or intrinsic resistance significantly reduced cell growth, increased cytotoxicity and apoptosis induced by drug treatment compared with control conditions (Fig. 1d and supplementary Fig. S1CE). Conversely, following miR-146a inhibitor transfection, the effect of drug treatment on cell cytotoxicity and apoptosis was reduced in sensitive LM16 cells (supplementary Fig. S1D).

Interestingly, we observed that resistance to BRAF/MEKi was extended to apoptosis determined by exposure to recombinant soluble TRAIL (sTRAIL), a functional mediator of cytotoxic immune cells (supplementary Fig. S1A). miR-146a overexpression in drug-resistant cells increased their sensitivity to sTRAIL-induced apoptosis (Fig. 1e, supplementary Fig. S1E), while miR-146a inhibition was protective (supplementary Fig. S1D). miR-146a expression levels were negatively correlated with sTRAIL-induced caspase 3/7 levels (Fig. 1f), further supporting the functional contribution of miR-146a to regulating cell apoptosis.

Forced expression of miR-146a reduces drug resistance via AKT/ERK and NFkB signaling cascade

BRAFi resistance is associated with a proinflammatory phenotype, characterized by the endogenous upregulation of the HIF1/NFkB-directed gene expression of cyto-chemokines, growth factors and other inflammatory mediators [6]. Because miR-146a is subject to NFkB regulation [8], we tested its involvement in drug resistance-induced NFkB modulation. miR-146a ectopic expression decreased NFkB activity (Fig. 2a), acting as a negative feedback loop regulator [9]. As expected, the prosurvival AKT/mTOR/ERK pathway and MAPK signaling pathway were downregulated, as indicated by the reduced pAKT, AKT, pp70S6K and pERK levels (Fig. 2a). The link between miR-146a downregulation and melanoma cell growth is sustained by the results of IPA analysis showing that putative miR-146a target genes, upregulated in GSE68841 [6], associate to “Cellular Growth and Proliferation” and to “Cellular Movement” Network functions.

Fig. 2

miR-146a regulates NFkB-related signaling pathways associated with resistance and directly targets PTGS2. a Immunoblot showing reduced levels of NFkB, AKT, p70S6K, and ERK in melanoma cells transfected with synthetic miR-146a (m-miR-146a) compared to mimic negative control transfectants (left). Reduction of phosphorylated proteins in miR-146a transfectants calculated by signal quantification (right). b Reduced expression of miR-146a target genes in transfectants (m-miR-146a: miR-146a mimic; S: mimic scrambled control). c Expression levels of PTGS2 and PTGES genes in resistant cell lines (R) compared to their sensitive counterparts (S). d Positive correlation between PTGS2 expression and IC50 values of BRAF/MEKi in melanoma cell lines (Spearman analysis). e Higher levels of COX2 protein expression in resistant cells (R) compared to their sensitive counterparts (S). Specific signals were quantified and expressed as the ratio of COX2/actin intensity. f Reduction of COX2 protein expression in LM16R resistant cells transfected with miR-146a mimic (m-miR-146a) and increased levels in LM36 parental sensitive cells transfected with miR-146a inhibitor (i-miR-146a). Changes in COX2 expression were calculated by signal quantification. P values were calculated by Student’s unpaired t test in A and F, and by Mann-Whitney U test in B. *: p < 0.05, **: p < 0.01, ***: p < 0.0001

To identify genes regulated by miR-146a via direct targeting or NFkB, we first extracted those differentially expressed in the GSE68841 dataset [6] from the list of 1667 NFkB-transcribed genes reported in the database built by Yang [25]. Twenty-four out of 527 common genes were validated direct target genes of miR-146a according to the miRWalk 2.0 database (supplementary Table S3) [26]. They included CCL5, CXCL8, EGFR and PTGS2 genes related to the proinflammatory and invasive signatures of resistant cells [6, 27]. In addition, among the genes indirectly regulated by miR-146a through NFkB activity, we found CCL2, IL6, and VEGFA, which contribute to the proinflammatory phenotype, CD274 (PDL1), which promotes melanoma cell proliferation [28], and BCL2L1 and MCL1 antiapoptotic genes. All these genes displayed reduced expression upon forced miR-146a expression in melanoma cells (Fig. 2b).

miR-146a impairs drug resistance by regulating the expression of PTGS2 gene encoding COX2.

Among the predicted miR-146a gene targets, the PTGS2 gene encoding COX2 attracted our attention because COX2 overexpression was reported to be associated with tumor progression and poor prognosis in melanoma [29,30,31]. COX2 transcript and protein expression levels were higher in several cells with acquired resistance compared to their sensitive counterparts and positively correlated with BRAF/MEKi resistance, revealing a direct association between COX2 and resistance (Fig. 2cde and supplementary Fig. S2B). Resistant cell lines also showed upregulation of PTGES, the gene encoding microsomal prostaglandin E2 synthase 1 (mPGES1), which is functionally coupled with COX2 in PGE2 production (Fig. 2c).

By luciferase reporter assays, we detected a significant decrease in luciferase activity upon melanoma cells transfection with the pLightSwitch-3’UTR-PTGS2 vector together with the miR-146a mimic when compared with controls, confirming a role for miR-146a in the post-transcriptional regulation of PTGS2 (supplementary Fig. S2A). Finally, forced miR-146a expression in resistant cells resulted in a decrease in COX2 protein levels, while inhibition of miR-146a in sensitive cells induced an increase in COX2 levels (Fig. 2f). These findings further confirm the involvement of miR-146a in COX2 regulation.

Targeting COX2 restores melanoma sensitivity to BRAF/MEKi

The above results prompted us to investigate whether targeting COX2 may impact BRAF/MEKi resistance. COX2 knockdown by small interfering RNA (siCOX2) reduced COX2 expression levels and enhanced melanoma cell sensitivity to BRAF/MEKi (Fig. 3a and supplementary Fig. S2CD). Pharmacologic inhibition of COX2 activity by the selective inhibitor celecoxib and BRAF/MEKi significantly reduced cell growth and increased cell cytotoxicity (Fig. 3b and supplementary Fig. S2F). The combined treatment resulted in an additive inhibition of cell growth (interaction index = 1), a result also confirmed by dose-response assays (supplementary Fig. S2FG).

Fig. 3

miR-146a spoils resistance of melanoma cells by repressing COX2 expression. a Silencing of COX2 by siRNA transfection (siCOX2) enhanced cell growth inhibition and increased cell cytotoxicity and apoptosis upon treatment with BRAF/MEKi. b Combined treatment with celecoxib increased cell growth inhibition and cytotoxicity by BRAF/MEKi. : interaction index = 1. c Increased apoptotic signaling after combination treatment with celecoxib and BRAF/MEKi. Phosphorylated and cleaved apoptotic factors showing increased expression in cells treated with the drug combination compared to BRAF/MEKi treated cells (Fold change > 1.3) by apoptosis signaling antibody array analysis. d Reduction of PGE2 release in culture media by treatment with the COX2 inhibitors celecoxib and NS398, alone or in combination with BRAF/MEKi. All experiments were carried with the LM47R cell line. P values were calculated by Student’s unpaired t test in B, and by one-way ANOVA followed by Bonferroni correction in D. *: p < 0.05, **: p < 0.01, ***: p < 0.0001

A screening of apoptosis-related proteins by antibody arrays showed that combined treatment significantly increased the phosphorylation of the proapoptotic proteins BAD and caspase 7 through activating SMAD2, JNK and p38/MAPK apoptosis-related signaling pathways. Moreover, induction of ATM signaling resulted in non canonical NFkB activation [32], determining the phosphorylation of components of the pathways involved in DNA damage-induced apoptosis, such as PARP1, p53 and TAK1 (Fig. 3c).

One consequence of COX2 overexpression is increased PGE2 production, which has multiple protumoral effects, including apoptosis resistance [17]. According to the high expression levels of COX2 and PTGES, resistant cell lines indeed secreted high levels of PGE2, which were significantly inhibited upon cell exposure to the COX2 inhibitors celecoxib and NS398, independent of BRAF/MEKi treatment (Fig. 3d and supplementary Fig. S2E). Taken together, these data indicate that the enhanced apoptotic effects induced by BRAF/MEKi in the presence of COX2 inhibition result from diminished PGE2 production, which reduces prosurvival signaling and activates apoptotic pathways and DNA damage signal transducers.

BRAF/MEKi resistance affects the miR-146a/COX2 axis in melanoma patients

To determine whether the miR-146a downregulation and COX2 upregulation observed in BRAF/MEKi-resistant melanoma cell lines in vitro also occurred in vivo in melanoma patients, we analyzed tumor samples surgically excised for local treatment from patients undergoing BRAF/MEKi therapy. Tissue immunostaining confirmed higher COX2 expression in tumor lesions progressing during BRAF/MEKi therapy compared to matched pretreatment lesions in 3 out of 7 tested cases (Fig. 4a and supplementary Table S4). High PTGS2 and PTGES and low miR-146a expression levels were detected in treatment-resistant tumors when compared to tumors from untreated patients (Fig. 4b). Despite the different numbers of treated and untreated samples analyzed and the diverse levels of gene expression detected in the tumor samples, significantly higher PTGS2 and PTGES was evident in tumors from treated patients. For miR-146a, although several untreated tumors show high expression levels, while on the contrary low levels are detected in treated specimens, no significant differences can be observed between untreated and treated samples.

Fig. 4

Treatment resistant tumors display miR-146a/COX2 axis deregulation. a COX2 immunostaining in tumors excised from patients progressing during BRAF/MEKi therapy (Post) compared to pre-therapy matched melanoma lesions (Pre). Scale bar: 10 μm. b Lower miR-146a and higher PTGS2 and PTGES expression levels in tumors progressing in patients during BRAF/MEKi therapy (Treated, n = 12) compared to tumors from unmatched untreated patients (Untreated, n = 31). *: p < 0.05, **: p < 0.01 by Mann-Whitney U test. c Heatmap showing the gene expression pattern of PTGS2, CXCL8, PTGES, IL6, VEGFA, CD274, MCL1 in progressing tumors (Post-treatment) compared to matched pre-therapy lesions (Pre-treatment) from seven metastatic melanoma patients receiving BRAFi therapy

Similar to what was observed in cell lines, miR-146a expression levels were inversely correlated with the levels of its target genes EGFR, CCL2 and BCL2L1 in melanoma metastases (supplementary Fig. S3). We then analyzed an in-house gene profiling dataset of seven matched metastatic tumor samples obtained from a cohort of patients before therapy and after the onset of treatment resistance for expression values of PTGS2 and of a number of additional direct or indirect miR-146a targets. In comparison to the matched pre-therapy sample, resistant lesions from five of the patients (Pts 1, 2, 3, 5 and 6) coordinately displayed higher expression of PTGS2. The remaining two patients (Pts 4 and 7) had instead lower expression of PTGS2 in their resistant sample. Different numbers of additional direct or indirect miR-146a targets, five in Pt 1, four in Pts 5 and 6, three in Pt 2 and two in Pt 3, displayed an expression trend similar to PTGS2. Among them, VEGFA, CD274, MCL1 were upregulated in most of the resistant lesions, whereas PTGES increase was never observed (Fig. 4c). Although an heterogenous pattern in term of gene expression levels was clearly displayed among patients, deregulation of the miR-146a/COX2 axis occurs in a subset of melanoma patients and is associated with the development of BRAF/MEKi drug resistance.

Functional relevance of miR-146a/COX2 axis manipulation in BRAF/MEKi resistance

To investigate the effects of miR-146a overexpression and COX2 inhibition on melanoma tumors, we performed 3D experiments with tumor explants using a short-term dynamic culture system in a bioreactor [24]. The preservation of TME architecture and cellularity by this culture system was confirmed in preliminary experiments (supplementary Fig. S4). Compared to control conditions, forced miR-146a expression led to increased miR-146a levels accompanied by decreased PTGS2 transcription and downregulation of most of the studied direct and indirect miR-146a target genes (Fig. 5ac and supplementary Fig. S5A).

Fig. 5

miR-146a overexpression reduces COX2 and increases drug sensitivity in 3D cultures from BRAF/MEKi-resistant tumors. a Overexpression of miR-146a upon transfection of specific mimic (m-miR-146a) in 3D tumor explants (left) downregulates COX2 expression compared to scrambled-transfected control (S) (right). b Reduced COX2 immunostaining in 3D cultures upon miR-146a forced expression (m-miR-146a) compared to mimic scrambled control (S). c Regulation of PTGS2 and miR-146a target genes (upper panel), and decreased release of CCL2, IL6, IL8 and VEGFA (lower panel) upon treatment with miR-146a mimic (m-miR-146a) and BRAF/MEKi in 3D cultures of resistant tumors in comparison to scrambled-transfected control (S). AU: arbitrary units. d PGE2 release in culture media from 3D tumor explants upon celecoxib treatment. **: p < 0.01 by Student’s unpaired t test. e Effects of treatment with BRAF/MEKi combined with transfection of miR-146a mimic (m-miR-146a) or of scrambled control (S) in resistant tumors: reduction of COX2 and of Ki67 positive cells and increase of cleaved caspase 3 immunostaining. Scale bars: 10 μm

In a set of 3D cultures from the lesions of 5 treatment-resistant metastatic patients, forced miR-146a upregulation led to reduced COX2 protein and gene levels and decreased PTGES, EGFR, CXCL8, CCL5, CCL2, IL6 and VEGFA, CD274, MCL1 and BCL2L1 levels (Fig. 5bc). These results indicated that downregulated miR-146a expression occurring through the acquisition of resistance to BRAF/MEKi-activated tumor-derived components contributes to melanoma survival and the generation of a tumor-promoting immunosuppressive TME; this finding was also confirmed by the reduced CCL2, IL6, IL8 and VEGFA secretion in the culture supernatants (Fig. 5c). Moreover, these observations led us to test the relevance of COX2 inhibition on PGE2 production in 3D tumor explants from a BRAF/MEKi-resistant patient and a treatment naïve patient. In both tumors, the COX2 inhibitor celecoxib significantly reduced PGE2 release, demonstrating that COX2 activity regulates PGE2 production also in tumor tissue (Fig. 5d). Consistently, a tissue analysis revealed that forced miR-146a expression remarkably reduced tumor cell proliferation in 3D cultured explants upon BRAF/MEKi treatment. These effects were observed in 3D cultures of BRAF-mutated melanoma tumors from naïve patients (supplementary Fig. S5B) and confirmed in tumors from patients undergoing treatment (Fig. 5e).

Our results strongly support a central role for the miR-146a/COX2 axis in modulating key components of tumor growth and survival associated with the onset of BRAF/MEKi resistance.


In the present study, we depict for the first time a role of miR-146a in controlling BRAF/MEKi resistance in melanoma cells, through a mechanism involving COX2. Chronic exposure to BRAFi induces miR-146a decrease in melanoma cells and increases expression of COX2, which boosts tumor cell survival to BRAFi resistance. Interestingly, we found that the increased BRAF/MEKi sensitivity induced by forced miR-146a expression in BRAF/MEKi-resistant cells were achieved by regulating COX2 expression. Our findings are consistent with previous observations demonstrating an inverse relationship between miR-146a and COX2 in lung cancer cells [19]. The analysis of tumors progressing in patients during BRAF/MEKi treatment revealed that miR-146a downregulation and the associated COX2 upregulation are not unique to in vitro models of BRAF/MEKi resistance as also occur in tumors in vivo. Due to the multiplicity of molecular alterations that characterize resistant melanoma, these results point to an important role of the miR-146a/COX2 axis in resistance. Moreover, the inhibition of COX2 by forced miR-146a expression in drug-resistant melanoma lesions is associated with the downregulation of immunosuppressive cyto-chemokines contributing to melanoma survival.

Decreased miR-146a expression was previously associated with drug resistance in lung, cervical, and hepatocellular cancer cells, and its overexpression suppressed cell growth and migration and improved drug sensitivity by inducing apoptosis [7]. Consistent with these findings, our study shows that miR-146a is significantly downregulated in six melanoma cell lines with acquired resistance and poorly expressed in two lines derived from drug-resistant tumors, and that its overexpression can restore sensitivity to BRAF/MEKi treatment by affecting both cell proliferation and apoptosis although with a different effect extent in the studied cell lines. Several studies have shown that decreased miR-146a expression leads to the overexpression of NFkB-mediated inflammatory factors contributing to immunosuppression and disease progression, and that miR-146a replacement can prevent the inflammatory state [8, 33]. Activation of NFkB signaling, along with low expression and activity of MITF, represent hallmarks of melanoma resistance to targeted therapy, predictive of poor prognosis for patients treated with BRAF/MEKi. NFkB-high transcriptional state could be present prior to therapy or could be induced by the acquisition of drug resistance. TNFα is a strong NFkB agonist, able to induce transcriptional changes promoting drug resistance [34]. In this view, we speculate that the downregulation of miR-146a could be considered an epigenetic mechanism contributing to the maintenance of the NFkB-high transcriptional state associated to resistance, because of its negative regulation of NFkB signaling and TNFα expression. In line with these findings, our experiments demonstrated that miR-146a ectopic expression in resistant cells modulated the levels of NFkB activity, downregulated the AKT/mTOR/ERK and MAPK signaling pathways, and reduced the expression of NFkB-induced mediators involved in cell growth and survival, including COX2 and EGFR, the immunosuppressive molecules IL8, CCL5, CCL2, IL6, VEGFA and PDL1, and the antiapoptotic genes BCL2L1 and MCL1.

miR-146a is involved in controlling the inflammatory response of innate immune system cells, particularly monocytes/macrophages [35]. Several studies have reported that miR-146a is associated with the negative regulation of immune activation and cancer-related immunosuppression [36]. For instance, miR-146a knock-out mice with melanoma survived longer [37], and miR-146a favors immunosuppression by increasing the regulatory T cell population in colorectal cancer [38]. In addition, miR-146a, one of the most studied myeloid miR [39], has been reported to serve as negative feedback modulator in the TLR4-mediated activation of NFkB–related genes [8] and promote M2 polarization in both humans and mice [40]. As we also reported, high miR-146a levels were associated with the induction of myeloid suppressor cells (MDSCs) and resistance to immunotherapy in melanoma patients [41]. Notably, miR-146a upregulation with concomitant increased COX2, PDL1, VEGFA, CCL2, IL6, IL8 and MCL1 expression characterized the MDSCs induced in vitro by melanoma extracellular vesicles from monocytes from healthy donors, suggesting cell type-specific regulation of the miR-146a/COX2 axis.

COX2 is an inducible enzyme essential in the biosynthesis of PGE2 through the enzyme mPGES1. Overexpression of COX2 is reportedly associated with a dismal prognosis in several tumor types [29, 30, 42]. Constitutive COX2 expression is triggered by pathways activated by oncogenic stimuli and cyto-chemokines, such as MAPK, PI3K/AKT and NFkB that are frequently hyperexpressed in most of resistant tumors [17]. Upregulation of the COX2/PGE2 axis favors proliferation, angiogenesis, invasion, apoptosis resistance, and the activation of immunosuppressive cells contributing to tumor progression and therapy resistance [43]. We showed that COX2 and PTGES are both overexpressed in several BRAF/MEKi-resistant melanoma cells and tumor tissues, and that COX2 knockdown significantly enhanced drug effects in melanoma cell. The COX2 selective inhibitor celecoxib alone and in combination with BRAF/MEKi has shown tumor inhibitory effects in preclinical melanoma studies [44, 45]. Clinical trials have also reported the effect of COX2 inhibitors, alone or combined with other treatment modalities, in cancer patients [17, 18]. However, there are no clinical reports investigating the treatment effect of COX2 inhibitors on BRAF/MEKi melanoma resistance. Moreover, it has been reported that COX2 can modulate PDL1 expression, and that celecoxib reduces PDL1 tumor expression in vitro and favors cytotoxic T cell responses [46, 47]. Other studies have described that COX2/mPGES1/PGE2 signaling regulates PDL1 and IDO in tumor-associated macrophages and MDSC cells, thus inducing an immunosuppressive phenotype [48,49,50]. In line with these reports, our data showed a coordinated expression of COX2 and PDL1 in resistant melanoma lesions and in melanoma cells upon forced miR-146a expression.

Our data show that celecoxib treatment reduced PGE2 production in resistant cells and in 3D cultures of resistant tumors, and its combination with BRAF/MEKi increased the drugs’ effect. Celecoxib-induced cell death is characterized by the activation of intrinsic and extrinsic apoptosis pathways and of NFkB signaling cascades induced by DNA damage, in line with published literature [18, 32, 51]. Our results support the targeting of COX2 by direct inhibition with selective drugs or by miR-146a delivery represents an option to improve the response to targeted therapy in melanoma patients.


We here report a tumor-suppressive role of miR-146a in BRAF/MEKi-resistant melanoma and demonstrate that the manipulation of miR-146a/COX2/PGE2 axis can restore BRAF/MEKi chemosensitivity (Fig. 6).

Fig. 6

Schematic representation of the miR-146a/COX2 axis associated with BRAF/MEKi resistance

Availability of data and materials

All data generated or analyzed during this study are included either in this article or in the supplementary information files.



Actin Beta


BRAF inhibitor, MEK inhibitor


Cell Counting Kit 8


C-C Motif Chemokine Ligand 2


C-C Motif Chemokine Ligand 5






C-X-C Motif Chemokine Ligand 8


Epidermal Growth Factor Receptor


Extracellular Signal–Regulated Kinase


False Discovery Rate


Hypoxia-Inducible Factor 1


Indoleamine 2,3-Dioxygenase 1


Interleukin 6


Ingenuity Pathway Analysis


Interleukin 1 Receptor Associated Kinase 1


Integrin Subunit Alpha V


Lactate Dehydrogenase


Mitogen Activated Protein Kinase


Myeloid Cell Leukemia 1


Myeloid Derived Suppressor Cells


Melanocyte Inducing Transcription Factor


Microsomal Prostaglandin E Synthase-1


Mammalian Target Of Rapamycin


Nuclear Factor kappa-light-chain-enhancer of activated B cells


Programmed Death-Ligand 1


Prostaglandin E2


Phosphatidylinositol-4,5-Bisphosphate 3-Kinase


Phosphatase and Tensin Homolog


Prostaglandin E Synthase


Prostaglandin-Endoperoxide Synthase 2


Quantitative Real-Time Polymerase Chain Reaction


Rotary Cell Culture System Bioreactor


Rho Associated Coiled-Coil Containing Protein Kinase 1


Ribosomal Protein L13a


Small interfering RNA COX2


Short Tandem Repeat


Soluble TNF-related apoptosis-inducing ligand


Toll Like Receptor 4


Tumor Microenvironment


Tumor Necrosis Factor alpha


TNF Receptor Associated Factor 6


Vascular Endothelial Growth Factor A


  1. 1.

    Kozar I, Margue C, Rothengatter S, Haan C, Kreis S. Many ways to resistance: how melanoma cells evade targeted therapies. Biochim Biophys Acta Rev Cancer. 2019;1871(2):313–22.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  2. 2.

    Hugo W, Shi H, Sun L, Piva M, Song C, Kong X, Moriceau G, Hong A, Dahlman KB, Johnson DB, Sosman JA, Ribas A, Lo RS. Non-genomic and immune evolution of melanoma acquiring MAPKi resistance. Cell. 2015;162(6):1271–85.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  3. 3.

    Sarkar D, Leung EY, Baguley BC, Finlay GJ, Askarian-Amiri ME. Epigenetic regulation in human melanoma: past and future. Epigenetics. 2015;10(2):103–21.

    PubMed  PubMed Central  Article  Google Scholar 

  4. 4.

    Diaz-Martinez M, Benito-Jardon L, Teixido J. New insights in melanoma resistance to BRAF inhibitors: a role for microRNAs. Oncotarget. 2018;9(83):35374–5.

    PubMed  PubMed Central  Article  Google Scholar 

  5. 5.

    Fattore L, Costantini S, Malpicci D, Ruggiero CF, Ascierto PA, Croce CM, Mancini R, Ciliberto G. MicroRNAs in melanoma development and resistance to target therapy. Oncotarget. 2017;8(13):22262–78.

    PubMed  PubMed Central  Article  Google Scholar 

  6. 6.

    Vergani E, Di Guardo L, Dugo M, Rigoletto S, Tragni G, Ruggeri R, Perrone F, Tamborini E, Gloghini A, Arienti F, Vergani B, Deho P, De Cecco L, Vallacchi V, Frati P, Shahaj E, Villa A, Santinami M, De Braud F, Rivoltini L, Rodolfo M. Overcoming melanoma resistance to vemurafenib by targeting CCL2-induced miR-34a, miR-100 and miR-125b. Oncotarget. 2016;7(4):4428–41.

    PubMed  Article  PubMed Central  Google Scholar 

  7. 7.

    Iacona JR. Lutz CS: miR-146a-5p: expression, regulation, and functions in cancer. Wiley Interdiscip Rev RNA. 2019;10(4):e15338.

    Google Scholar 

  8. 8.

    Taganov KD, Boldin MP, Chang KJ, Baltimore D. NF-kappaB-dependent induction of microRNA miR-146, an inhibitor targeted to signaling proteins of innate immune responses. Proc Natl Acad Sci U S A. 2006;103(33):12481–6.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  9. 9.

    Paik JH, Jang JY, Jeon YK, Kim WY, Kim TM, Heo DS, Kim CW. MicroRNA-146a downregulates NFkappaB activity via targeting TRAF6 and functions as a tumor suppressor having strong prognostic implications in NK/T cell lymphoma. Clin Cancer Res. 2011;17(14):4761–71.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  10. 10.

    Ma X, Becker Buscaglia LE, Barker JR, Li Y. MicroRNAs in NF-kappaB signaling. J Mol Cell Biol. 2011;3(3):159–66.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  11. 11.

    El Gazzar M, Church A, Liu T, McCall CE. MicroRNA-146a regulates both transcription silencing and translation disruption of TNF-α during TLR4-induced gene reprogramming. J Leukoc Biol. 2011;90(3):509–19.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  12. 12.

    Forloni M, Dogra SK, Dong Y, Conte D Jr, Ou J, Zhu LJ, Deng A, Mahalingam M, Green MR. Wajapeyee N: miR-146a promotes the initiation and progression of melanoma by activating Notch signaling. Elife. 2014;3:e01460.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  13. 13.

    Raimo M, Orso F, Grassi E, Cimino D, Penna E, De Pitta C, Stadler MB, Primo L, Calautti E, Quaglino P, Provero P. Taverna D: miR-146a exerts differential effects on melanoma growth and Metastatization. Mol Cancer Res. 2016;14(6):548–62.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  14. 14.

    Anelli V, Ordas A, Kneitz S, Sagredo LM, Gourain V, Schartl M, Meijer AH, Mione M. Ras-induced miR-146a and 193a target Jmjd6 to regulate melanoma progression. Front Genet. 2018;9:675.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  15. 15.

    Pu W, Shang Y, Shao Q. Yuan X: miR-146a promotes cell migration and invasion in melanoma by directly targeting SMAD4. Oncol Lett. 2018;15(5):7111–7.

    PubMed  PubMed Central  Google Scholar 

  16. 16.

    Gomez-Lira M, Ferronato S, Orlandi E, Dal Molin A, Malerba G, Frigerio S, Rodolfo M, Romanelli MG. Association of microRNA 146a polymorphism rs2910164 and the risk of melanoma in an Italian population. Exp Dermatol. 2015;24(10):794–5.

    PubMed  Article  PubMed Central  Google Scholar 

  17. 17.

    Goradel HN, Najafi M, Salehi E, Farhood B, Mortezaee K. Cyclooxygenase-2 in cancer: a review. J Cell Physiol. 2019;234(5):5683–99.

    Article  CAS  Google Scholar 

  18. 18.

    Toloczko-Iwaniuk N, Dziemianczyk-Pakiela D, Nowaszewska BK, Celinska-Janowicz K, Miltyk W. Celecoxib in Cancer therapy and prevention - review. Curr Drug Targets. 2019;20(3):302–15.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  19. 19.

    Cornett AL, Carol S. Regulation of COX-2 expression by miR-146a in lung cancer cells. RNA. 2014;20(9):1419–30.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  20. 20.

    Waggott DM. NanoStringNorm: normalize NanoString miRNA and mRNA data. R Package Version. 2017;1(2):1

    Google Scholar 

  21. 21.

    Rambow F, Job B, Petit V, Gesbert F, Delmas V, Seberg H, Meurice G, Van Otterloo E, Dessen P, Robert C, Gautheret D, Cornell RA, Sarasin A, Larue L. New functional signatures for understanding melanoma biology from tumor cell lineage-specific analysis. Cell Rep. 2015;13(4):840–53.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  22. 22.

    Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W. Smyth GK: limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  23. 23.

    Du P, Kibbe WA. Lin SM: lumi: a pipeline for processing Illumina microarray. Bioinformatics. 2008;24(13):1547–8.

    CAS  PubMed  Article  Google Scholar 

  24. 24.

    Ferrarini M, Steimberg N, Ponzoni M, Belloni D, Berenzi A, Girlanda S, Caligaris-Cappio F, Mazzoleni G, Ferrero E. Ex-vivo dynamic 3-D culture of human tissues in the RCCS bioreactor allows the study of multiple myeloma biology and response to therapy. PLoS One. 2013;8(8):e71613.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  25. 25.

    Yang Y, Wu J, Wang J. A database and functional annotation of NF-κB target genes. Int J Clin Exp Med. 2016;9(5):7986–95.

    CAS  Google Scholar 

  26. 26.

    miRWalk 2.0. Accessed 12 December 2019.

  27. 27.

    Dugo M, Nicolini G, Tragni G, Bersani I, Tomassetti A, Colonna V, Del Vecchio M, De Braud F, Canevari S, Anichini A, Sensi M. A melanoma subtype with intrinsic resistance to BRAF inhibition identified by receptor tyrosine kinases gene-driven classification. Oncotarget. 2015;6(7):5118–33.

    PubMed  PubMed Central  Article  Google Scholar 

  28. 28.

    Clark C, Gupta H, Sareddy G, Pandeswara S, Lao S, Yuan B, Drerup JM, Padron A, Conejo-Garcia J, Murthy K, Liu Y, Turk MJ, Thedieck K, Hurez V, Li R, Vadlamudi R, Curiel TJ. Tumor-intrinsic PD-L1 signals regulate cell growth, pathogenesis, and autophagy in ovarian Cancer and melanoma. Cancer Res. 2016;76(23):6964–74.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  29. 29.

    Becker MR, Siegelin MD, Rompel R, Enk AH, Gaiser T. COX-2 expression in malignant melanoma: a novel prognostic marker? Melanoma Res. 2009;19(1):8–16.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  30. 30.

    Panza E, De Cicco P, Ercolano G, Armogida C, Scognamiglio G, Anniciello AM, Botti G, Cirino G, Ianaro A. Differential expression of cyclooxygenase-2 in metastatic melanoma affects progression free survival. Oncotarget. 2016;7(35):57077–85.

    PubMed  PubMed Central  Article  Google Scholar 

  31. 31.

    Soares CD, Borges CF, Sena-Filho M, Almeida OP, Stelini RF, Cintra ML, Graner E, Zecchin KG, Jorge J. Prognostic significance of cyclooxygenase 2 and phosphorylated Akt1 overexpression in primary nonmetastatic and metastatic cutaneous melanomas. Melanoma Res. 2017;27(5):448–56.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  32. 32.

    Wang W, Mani AM, Wu ZH. DNA damage-induced nuclear factor-kappa B activation and its roles in cancer progression. J Cancer Metastasis Treat. 2017;3:45–59.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  33. 33.

    Perry MM, Moschos SA, Williams AE, Shepherd NJ, Larner-Svensson HM, Lindsay MA. Rapid changes in microRNA-146a expression negatively regulate the IL-1beta-induced inflammatory response in human lung alveolar epithelial cells. J Immunol. 2008;180:5689–98.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  34. 34.

    Konieczkowski DJ, Johannessen CM, Abudayyeh O, Kim JW, Cooper ZA, Piris A, Frederick DT, Barzily-Rokni M, Straussman R, Haq R, Fisher DE, Mesirov JP, Hahn WC, Flaherty KT, Wargo JA, Tamayo P, Garraway LA. A melanoma cell state distinction influences sensitivity to MAPK pathway inhibitors. Cancer Discov. 2014;4(7):816–27.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  35. 35.

    Labbaye C, Testa U. The emerging role of MIR-146A in the control of hematopoiesis, immune function and cancer. J Hematol Oncol. 2012;5:13–8722-5-13.

    Article  CAS  Google Scholar 

  36. 36.

    Essandoh K, Li Y, Huo J, Fan GC. MiRNA-mediated macrophage polarization and its potential role in the regulation of inflammatory response. Shock. 2016;46(2):122–31.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  37. 37.

    Mastroianni J, Stickel N, Andrlova H, Hanke K, Melchinger W, Duquesne S, Schmidt D, Falk M, Andrieux G, Pfeifer D, Dierbach H, Schmitt-Graeff A, Meiss F, Boerries M. Zeiser R: miR-146a controls immune response in the melanoma microenvironment. Cancer Res. 2019;79(1):183–95.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  38. 38.

    Khorrami S, Zavaran Hosseini A, Mowla SJ, Soleimani M, Rakhshani N, Malekzadeh R. MicroRNA-146a induces immune suppression and drug-resistant colorectal cancer cells. Tumour Biol. 2017;39(5):1010428317698365.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  39. 39.

    O’Connell RM, Zhao JL, Rao DS. MicroRNA function in myeloid biology. Blood. 2011;118(11):2960–9.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  40. 40.

    Lee HM, Kim TS, Jo EK. MiR-146 and miR-125 in the regulation of innate immunity and inflammation. BMB Rep. 2016;49(6):311–8.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  41. 41.

    Huber V, Vallacchi V, Fleming V, Hu X, Cova A, Dugo M, Shahaj E, Sulsenti R, Vergani E, Filipazzi P, De Laurentiis A, Lalli L, Di Guardo L, Patuzzo R, Vergani B, Casiraghi E, Cossa M, Gualeni A, Bollati V, Arienti F, De Braud F, Mariani L, Villa A, Altevogt P, Umansky V, Rodolfo M, Rivoltini L. Tumor-derived microRNAs induce myeloid suppressor cells and predict immunotherapy resistance in melanoma. J Clin Invest. 2018;128(12):5505–16.

    PubMed  PubMed Central  Article  Google Scholar 

  42. 42.

    Wang ZM, Liu J, Liu HB, Ye M, Zhang YF, Yang DS. Abnormal COX2 protein expression may be correlated with poor prognosis in Oral Cancer: a meta-analysis. Biomed Res Int. 2014;364207.

  43. 43.

    Tong D, Liu Q, Wang LA, Xie Q, Pang J, Huang Y, Wang L, Liu G, Zhang D, Lan W, Jiang J. The roles of the COX2/PGE2/EP axis in therapeutic resistance. Cancer Metastasis Rev. 2018;37(2–3):355–68.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  44. 44.

    Gowda R, Kardos G, Sharma A, Singh S, Robertson GP. Nanoparticle-based Celecoxib and Plumbagin for the synergistic treatment of melanoma. Mol Cancer Ther. 2017;16(3):440–52.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  45. 45.

    Escuin-Ordinas H, Atefi M, Fu Y, Cass A, Ng C, Huang RR, Yashar S, Comin-Anduix B, Avramis E, Cochran AJ, Marais R, Lo RS, Graeber TG, Herschman HR, Ribas A. COX-2 inhibition prevents the appearance of cutaneous squamous cell carcinomas accelerated by BRAF inhibitors. Mol Oncol. 2014;8(2):250–60.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  46. 46.

    Markosyan N, Chen EP, Evans RA, Ndong V, Vonderheide RH, Smyth EM. Mammary carcinoma cell derived cyclooxygenase 2 suppresses tumor immune surveillance by enhancing intratumoral immune checkpoint activity. Breast Cancer Res. 2013;15(5):R75.

    PubMed  PubMed Central  Article  Google Scholar 

  47. 47.

    Botti G, Fratangelo F, Cerrone M, Liguori G, Cantile M, Anniciello AM, Scala S, D'Alterio C, Trimarco C, Ianaro A, Cirino G, Caraco C, Colombino M, Palmieri G, Pepe S, Ascierto PA, Sabbatino F, Scognamiglio G. COX-2 expression positively correlates with PD-L1 expression in human melanoma cells. J Transl Med. 2017;15(1):46.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  48. 48.

    Prima V, Kaliberova LN, Kaliberov S, Curiel DT, Kusmartsev S. COX2/mPGES1/PGE2 pathway regulates PD-L1 expression in tumor-associated macrophages and myeloid-derived suppressor cells. Proc Natl Acad Sci U S A. 2017;114(5):1117–22.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  49. 49.

    Hennequart M, Pilotte L, Cane S, Hoffmann D, Stroobant V, Plaen E, Van den Eynde BJ. Constitutive IDO1 expression in human tumors is driven by Cyclooxygenase-2 and mediates intrinsic immune resistance. Cancer Immunol Res. 2017;5(8):695–709.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  50. 50.

    Kim SH, Roszik J, Cho SN, Ogata D, Milton DR, Peng W, Menter DG, Ekmekcioglu S, Grimm EA. The COX2 effector microsomal PGE2 synthase 1 is a regulator of immunosuppression in cutaneous melanoma. Clin Cancer Res. 2019;25(5):1650–63.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  51. 51.

    Jendrossek V. Targeting apoptosis pathways by Celecoxib in cancer. Cancer Lett. 2013;332(2):313–24.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

Download references


We thank Martina Filugelli, Paola Deho and Agata Cova for their technical contribution.


This study was funded by CARIPLO Foundation (2015–0911 to VV), by AIRC (IG 17462 to MR) and AIRC Special Program Innovative Tools for Cancer Risk Assessment and early Diagnosis 5 × 1000 (no. 12162 to LR and MSe). ES is a PhD Student of Open University Milton Keynes, UK.

Author information




EV, MR and VV conceived the study and designed the experiments. EV, SF, ES, BVe, MF and MGL performed the experiments. EV, MC, BVa, MSe, BEL and EF interpreted the data. MD, CG, MSe and LL performed the bioinformatics and statistical analysis. MDV, LDG, FG, IM, MSa, ET provided the clinical samples. EV, VV and MR wrote the manuscript. LR, MSe and VH revised the manuscript. The author(s) read and approved the final manuscript.

Corresponding author

Correspondence to Viviana Vallacchi.

Ethics declarations

Ethics approval and consent to participate

The studies with human samples were approved by the Institutional Review Board and by an Independent Ethics Committee, and patient samples were included upon informed consent.

Consent for publication

All authors agree with the content of the manuscript.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Supplementary information

Additional file 1: Supplementary Table S1.

List of used antibodies and molecular probes. Supplementary Table S2. Differentially expressed miR between cell lines with acquired resistance to vemurafenib and their sensitive counterparts. Supplementary Table S3. Genes directly targeted by miR-146a among the NFkB-regulated genes differentially expressed in LM16R/LM16 dataset (GSE68841). Supplementary Table S4. Characteristics of patients treated with vemurafenib and their matched pre- and post-treatment melanoma metastasis (GSE141484).

Additional file 2: Supplementary Fig. S1.

Altered miR-146a expression influences BRAF/MEKi sensitivity and apoptosis in melanoma. A) Boxplots representing IC50 values of 6 matched PLX4032-resistant (R) and sensitive (S) melanoma cell lines to PLX4032 (vemurafenib), BRAFi (dabrafenib), MEKi (trametinib), to the combined treatment BRAF/MEKi, and to sTRAIL-induced apoptosis. B) Inverse correlation between miR-146a expression levels and IC50 values of PLX4032, BRAFi and MEKi in melanoma cell lines (Spearman analysis). C) Forced expression of miR-146a (+m-miR-146a) in LM16R cell line and in LM69 and LM70 short term cultures increased the effects of PLX4032 treatment as shown by reduced cell growth and increased cell cytotoxicity and apoptosis, evaluated by CCK8, LDH and caspase 8 and 3/7 activity. D) Inhibition of miR-146a expression (i-miR-146a) in LM16 cells increased cell proliferation and decreased the release of LDH and the apoptosis rate as evaluated by CCK8, LDH and caspase 3/7 activity. E) Overexpression of miR-146a (m-miR-146a) upon PLX4032, BRAF/MEKi and sTRAIL treatments in LM47R cells increased cell cytotoxicity and apoptosis, as evaluated by LDH and caspase 8 and 3/7 activity. Data are plotted compared to scrambled control. *: p < 0.05, **: p < 0.01, ***: p < 0.0001 by Student’s unpaired t test. Supplementary Fig. S2. COX2 inhibition increases sensitivity to PLX4032 and to BRAF/MEKi and reduces PGE2 release. A) Relative luciferase activity after co-transfection of PTGS2 3’UTR luciferase reporter vector or control vector with miR-146a mimic or mimic negative control. Experiment was performed in LM16R cells. B) Inverse correlation between miR-146a expression levels and IC50 values of PLX4032 in melanoma cell lines (Spearman analysis). C) Western blot analyses showing downregulation of COX2 after transfection with siRNA against COX2 (siCOX2) compared to scrambled control (S). COX2 protein levels were downregulated to 15%, as determined by quantification of the signal and expressed as the ratio of COX2/Actin intensity. Experiment was performed with LM47R cells. D) Silencing of COX2 by siRNA transfection enhanced cell growth inhibition and increased cell cytotoxicity and apoptosis upon treatment with PLX4032 compared to scrambled control. E) PGE2 release in culture supernatants following 24 h treatment with the COX2 inhibitors celecoxib and NS398, alone or in combination with BRAF/MEKi. F) Combined treatment with celecoxib increased cell growth inhibition and cytotoxicity by PLX4032. : interaction index = 1. G) Dose-dependent effect of celecoxib alone (5, 10, 20, 40, 50 and 80 μM) and combined with PLX4032 (3 μM) on cell growth. : interaction index = 1 at the dose of 50 μM of celecoxib combined with PLX4032. P values were calculated by Student’s unpaired t test in A, D and F, by one-way ANOVA followed by Bonferroni correction in E, and by two-way ANOVA followed by Bonferroni correction in G. *: p < 0.05, **: p < 0.01, ***: p < 0.0001. Supplementary Fig. S3. miR-146a expression levels are inversely correlated with levels of target genes. Analysis of correlation between miR-146a expression and its target genes EGFR, CCL2 and BCL2L1 in metastatic melanoma specimens. The rp and p values resulting from Pearson analysis are shown. Supplementary Fig. S4. Dynamic culture in bioreactor preserves TME architecture of 3D tumor explants. Representative histological images of H&E-stained 3D tumor sections after 3 days culture in bioreactor. The marked areas in the left panels are shown at higher magnification in the right panels to display the preserved cellularity and tumor tissue histo-architecture. Scale bar 80 μM. Supplementary Fig. S5. miR-146a overexpression increases drug sensitivity in 3D cultures of tumor explants. A) Modulation of direct and indirect miR-146a target genes upon transfection of miR-146a mimic (m-miR-146a) compared to scrambled-transfected control (S) in 3D cultures of BRAF-mutated melanoma tumors from naïve patients. AU: arbitrary units. B) Decreased staining for Ki67 and pERK proliferating cells and increased Caspase 3 positive cells in 3D cultures upon treatment with BRAF/MEKi and transfection with miR-146a (m-miR-146a) or scrambled control (S). Scale bar 10 um.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Vergani, E., Dugo, M., Cossa, M. et al. miR-146a-5p impairs melanoma resistance to kinase inhibitors by targeting COX2 and regulating NFkB-mediated inflammatory mediators. Cell Commun Signal 18, 156 (2020).

Download citation


  • Melanoma resistance
  • BRAF/MEK inhibitors
  • microRNA
  • miR-146a-5p
  • COX2