Activities of multiple cancer-related pathways are associated with BRAF mutation and predict the resistance to BRAF/MEK inhibitors in melanoma cells

Drug resistance is a major obstacle in the targeted therapy of melanoma using BRAF/MEK inhibitors. This study was to identify BRAF V600E-associated oncogenic pathways that predict resistance of BRAF-mutated melanoma to BRAF/MEK inhibitors. We took in silico approaches to analyze the activities of 24 cancer-related pathways in melanoma cells and identify those whose activation was associated with BRAF V600E and used the support vector machine (SVM) algorithm to predict the resistance of BRAF-mutated melanoma cells to BRAF/MEK inhibitors. We then experimentally confirmed the in silico findings. In a microarray gene expression dataset of 63 melanoma cell lines, we found that activation of multiple oncogenic pathways preferentially occurred in BRAF-mutated melanoma cells. This finding was reproduced in 5 additional independent melanoma datasets. Further analysis of 46 melanoma cell lines that harbored BRAF mutation showed that 7 pathways, including TNFα, EGFR, IFNα, hypoxia, IFNγ, STAT3, and MYC, were significantly differently expressed in AZD6244-resistant compared with responsive melanoma cells. A SVM classifier built on this 7-pathway activation pattern correctly predicted the response of 10 BRAF-mutated melanoma cell lines to the MEK inhibitor AZD6244 in our experiments. We experimentally showed that TNFα, EGFR, IFNα, and IFNγ pathway activities were also upregulated in melanoma cell A375 compared with its sub-line DRO, while DRO was much more sensitive to AZD6244 than A375. In conclusion, we have identified specific oncogenic pathways preferentially activated in BRAF-mutated melanoma cells and a pathway pattern that predicts resistance of BRAF-mutated melanoma to BRAF/MEK inhibitors, providing novel clinical implications for melanoma therapy.