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.     

A simple but highly effective approach to evaluate the prognostic performance of gene expression signatures

Background

Highly parallel analysis of gene expression has recently been used to identify gene sets or ‘signatures’ to improve patient diagnosis and risk stratification. Once a signature is generated, traditional statistical testing is used to evaluate its prognostic performance. However, due to the dimensionality of microarrays, this can lead to false interpretation of these signatures.

Principal Findings

A method was developed to test batches of a user-specified number of randomly chosen signatures in patient microarray datasets. The percentage of random generated signatures yielding prognostic value was assessed using ROC analysis by calculating the area under the curve (AUC) in six public available cancer patient microarray datasets. We found that a signature consisting of randomly selected genes has an average 10% chance of reaching significance when assessed in a single dataset, but can range from 1% to ∼40% depending on the dataset in question. Increasing the number of validation datasets markedly reduces this number.

Conclusions

We have shown that the use of an arbitrary cut-off value for evaluation of signature significance is not suitable for this type of research, but should be defined for each dataset separately. Our method can be used to establish and evaluate signature performance of any derived gene signature in a dataset by comparing its performance to thousands of randomly generated signatures. It will be of most interest for cases where few data are available and testing in multiple datasets is limited.     

BRAF and RKIP are significantly decreased in cutaneous squamous cell carcinoma

Background. Actinic keratosis (AK) is a well-established pre-cancerous skin lesion that has the potential to progress to squamous cell carcinoma (SCC). However, little is known about the implication of BRAF and RKIP expression, or about the incidence of BRAF mutations in the formation of these cutaneous diseases. The RAS oncogene has been proposed to significantly contribute to skin cancer development. Moreover, numerous BRAF mutations have been detected in melanoma biopsy specimens and cell lines.

Objectives. This study aimed to measure the mRNA levels of the genes BRAF and RKIP in AK, as well as their possible implication in the progress of AK to SCC. All biopsy specimens were also screened for BRAF mutations within exons 11 and 15.

Patients and methods. Expression levels of the genes BRAF and RKIP were examined in 16 AKs and 12 SCCs by RT-qPCR. A novel allele-specific qPCR method, in combination with direct DNA sequencing, was performed in order to inspect the frequency of the V600E mutation in exon 15, as well as to examine the mutation status of the gene within exon 11.

Results. Significant down-regulation was noted for both genes in SCC, compared to normal tissue (for BRAF, P=0.002; for RKIP, P     

Concordance of somatic mutational profile in multiple primary melanomas

This study aimed to determine the frequency and concordance of BRAF and NRAS mutation in tumours arising in patients with multiple primary melanoma (MPM ). Patients with MPM managed at one of three tertiary referral centres in Melbourne, Australia, from 2010 to 2015 were included. Incident and subsequent melanomas underwent mutation testing. Cohen's kappa (κ) coefficient assessed agreement between incident and subsequent primary melanomas for both BRAF and NRAS mutation status (mutant versus wild‐type). Mutation testing of at least two primary tumours from 64 patients was conducted. There was poor agreement for both BRAF and NRAS mutation status between incident and subsequent melanomas (κ = 0.10, 95% CI ?0.10 to 0.42; κ = 0.06, 95% CI ?0.10 to 0.57, respectively). In view of the low concordance in BRAF mutation status between incident and subsequent melanomas, mutational analysis of metastatic tissue, rather than of a primary melanoma, in patients with MPM should be used to guide targeted therapy.     

Identification of PLX4032‐resistance mechanisms and implications for novel RAF inhibitors

BRAF inhibitors improve melanoma patient survival, but resistance invariably develops. Here we report the discovery of a novel BRAF mutation that confers resistance to PLX4032 employing whole‐exome sequencing of drug‐resistant BRAF^V600K melanoma cells. We further describe a new screening approach, a genome‐wide piggyBac mutagenesis screen that revealed clinically relevant aberrations (N‐terminal BRAF truncations and CRAF overexpression). The novel BRAF mutation, a Leu505 to His substitution (BRAF^L505H), is the first resistance‐conferring second‐site mutation identified in BRAF mutant cells. The mutation replaces a small nonpolar amino acid at the BRAF‐PLX4032 interface with a larger polar residue. Moreover, we show that BRAF^L505H, found in human prostate cancer, is itself a MAPK‐activating, PLX4032‐resistant oncogenic mutation. Lastly, we demonstrate that the PLX4032‐resistant melanoma cells are sensitive to novel, next‐generation BRAF inhibitors, especially the ‘paradox‐blocker’ PLX8394, supporting its use in clinical trials for treatment of melanoma patients with BRAF‐mutations.     

BRAF initiating mutations in the papillary thyroid carcinoma

Among genetic alterations most important for the initiation of papillary thyroid carcinoma (PTC) is mutation T1799A in the BRAF gene which is the most frequent event (54.5%) in this type of thyroid cancer. It is seen in all stages, from microcarcinoma through clinically overt disease to anaplastic cancer. It has been shown that BRAF mutation is correlated with PTC histotype. It is identified most frequently in classical PTC and in tall cell variant. Moreover, BRAF mutation is described more often in older patients, whereas in young patients RET/PTC rearrangements dominate. In PTC cases with BRAF mutation V600E the prognosis is poorer, with more cancer invasiveness, metastasis and recurrence. The presence of BRAF mutation is related to the specific gene expression signature, different than in cancer cases showing RET/PTC rearrangement or no known initiating mutation.     

BRAF mutation and RASSF1A expression in thyroid carcinoma of southern Italy

Aim of this work is to provide a detailed comparison of clinical‐pathologic features between well‐differentiated and poorly differentiated tumors according to their BRAF and RASSF1A status. We analyzed RASSF1A methylation by MSP and BRAF mutation by LCRT‐PCR with LightMix® kit BRAF V600E in neoplastic thyroid tissues. Immunohistochemical evaluation of RASSF1A expression was also performed by standard automated LSAB‐HRP technique. An overall higher degree of RASSF1A over‐expression than normal thyroid parenchyma surrounding tumors (P < 0.05) has been found in all malignant well‐differentiated lesions. Moreover, statistically significant higher levels of RASSF1A expression were observed in differentiated cancers associated to an inflammatory autoimmune background (P = 0.01). Amplifiable DNA for LC PCR with LightMix® kit BRAF V600E was obtained in nine PTCs, four FVPTCs, five ATCs, and one control. The V600E mutation was found in 13 of 18 (72%) tumors. BRAF was mutated in 6 of 9 (66%) classical PTC, in 2 of 4 (50%) follicular variant PTC and in all ACs (100%). The overall frequency of RASSF1A promoter methylation observed was 20.5% (9 cases out 44). Hypermethylation of RASSF1A in primary tumors was variable according to histotypes ranging from100% (5/5) in ACs to only 12.5% (4/32) in PTCs. We show a correlation between RASSF1A methylation status and RASSF1A protein expression. Finally, we conclude that BRAF V600E mutation and RASSF1A methylation were pathogenetic event restricted to a subgroup of PTC/FVPTCs in early stage and to clinically aggressive ATCs. J. Cell. Biochem. 114: 1174–1182, 2013. © 2012 Wiley Periodicals, Inc.     

Mutation load in melanoma is affected by MC1R genotype

Whole‐genome sequencing of matched germline and tumour pairs in a well‐characterized cohort of melanoma patients allowed investigation of associations between melanoma body site, age at melanoma onset and MC1R variant status with overall mutation burden and specific base pair changes observed in the corresponding melanoma. We observed statistically significant associations between mutation burden in melanoma and body site, age at onset and MC1R genotype, for both ultraviolet radiation (UVR) signature changes (C>T and CC>TT) and non‐UVR base pair substitutions, as well as with overall variant load.     

Analysis of BRAF and NRAS Mutation Status in Advanced Melanoma Patients Treated with Anti-CTLA-4 Antibodies: Association with Overall Survival?

Ipilimumab and tremelimumab are human monoclonal antibodies (Abs) against cytotoxic T-lymphocyte antigen-4 (CTLA-4). Ipilimumab was the first agent to show a statistically significant benefit in overall survival in advanced melanoma patients. Currently, there is no proven association between the BRAFV600 mutation and the disease control rate in response to ipilimumab. This analysis was carried out to assess if BRAFV600 and NRAS mutation status affects the clinical outcome of anti-CTLA-4-treated melanoma patients. This is a retrospective multi-center analysis of 101 patients, with confirmed BRAF and NRAS mutation status, treated with anti-CTLA-4 antibodies from December 2006 until August 2012. The median overall survival, defined from the treatment start date with the anti-CTLA-4. Abs-treatment to death or till last follow up, of BRAFV600 or NRAS mutant patients (n = 62) was 10.12 months (95% CI 6.78–13.2) compared to 8.26 months (95% CI 6.02–19.9) in BRAFV600/NRASwt subpopulation (n = 39) (p = 0.67). The median OS of NRAS mutated patients (n = 24) was 12.1 months and although was prolonged compared to the median OS of BRAF mutated patients (n = 38, mOS = 8.03 months) or BRAFV600/NRASwt patients (n = 39, mOS = 8.26 months) the difference didn’t reach statistical significance (p = 0.56). 69 patients were able to complete 4 cycles of anti-CTLA-4 treatment. Of the 24 patients treated with selected BRAF- or MEK-inhibitors, 16 patients received anti-CTLA 4 Abs following either a BRAF or MEK inhibitor with only 8 of them being able to finish 4 cycles of treatment. Based on our results, there is no difference in the median OS in patients treated with anti-CTLA-4 Abs implying that the BRAF/NRAS mutation status alone is not sufficient to predict the outcome of patients treated with anti-CTLA-4 Abs.     

Prognostic significance of tumor iNOS and COX-2 in stage III malignant cutaneous melanoma

Purpose  New prognostic markers are needed for malignant melanoma. Inducible nitric oxide synthase (iNOS) and cyclooxygenase type 2 (COX-2) have been described to correlate with progression of melanoma. Moreover, activating mutations in BRAF/NRAS oncogenes are often detected in melanoma. The BRAF/NRAS mutation status and expression of COX-2 and iNOS were examined to compare their prognostic value for overall survival (OS) in stage III malignant cutaneous melanoma. Experimental design  The expression of iNOS and COX-2 in metastatic lymph nodes from 21 rapidly progressing (OS from date of diagnosis of stage III disease ≤14 months) and 17 slowly progressing (OS ≥60 months) stage III cutaneous melanoma patients was examined by immunohistochemistry. The presence of BRAF/NRAS mutations was analyzed using direct DNA sequencing. χ² exact trend test and logistic regression analysis were used for statistical analysis. Results  Both iNOS (P = 0.002) and COX-2 (P = 0.048) alone significantly predicted OS. The BRAF/NRAS mutation status did not significantly differ between patient groups, although iNOS significantly (P = 0.013) correlated with BRAF mutation frequency. Furthermore, the odds ratio (OR) with respect to OS of iNOS (OR = 10.4) was higher than that of COX-2 (OR = 5.6) and was stable in the multivariate analysis of OS together with disease stage IIIB/C, ulceration, number of metastatic lymph nodes, and Breslow tumor thickness. Conclusion  Our data show that iNOS is an independent and stronger prognostic factor for OS in stage III malignant cutaneous melanoma than COX-2.     

Assessment of association between BRAF-V600E mutation status in melanomas and clinical response to ipilimumab

Ipilimumab, a fully human monoclonal antibody against cytotoxic T lymphocyte antigen-4, has demonstrated significant improvement in overall survival in previously treated advanced melanoma patients. The BRAF inhibitor, vemurafenib, has shown up to 78% objective response rates in melanoma patients harboring the BRAF-V600E mutation but not in patients lacking the mutation. As an immune potentiator, the mechanism of action of ipilimumab may not be dependent of the activity of the BRAF pathway. To test this, we investigated whether the clinical activity of ipilimumab would be affected by the BRAF-V600E mutation status of the tumors. Thus, this retrospective analysis was carried using a set of tumor biopsies from a completed phase II clinical trial. CA184004 was a randomized, double-blind, multicenter trial of 82 previously treated or untreated patients with unresectable stage III/IV melanoma. Patients received ipilimumab 3 or 10 mg/kg every 3 weeks for four doses followed by maintenance dosing in eligible patients. The BRAF-V600E mutation status for 80 patients was determined in tumor biopsies by PCR-based assays. Data on disease control were available for 69 patients with evaluated BRAF-V600E mutation status. Rates of objective responses and stable disease in patients with BRAF-V600E mutation positive tumors (30%) were comparable to those in patients with the wild-type gene (~33%). Eleven patients displayed Durable Disease Control (DDC) of which 55% had BRAF-V600E mutation positive tumors and 45% did not. In the 48 patients showing no DDC, the mutation frequency was 50%. In this study, no association between BRAF-V600E mutation status of melanoma tumors and DDC after treatment with ipilimumab was detected.     

Neuroblastoma rat sarcoma mutated melanoma: That’s what we got so far

Neuroblastoma rat sarcoma (NRAS) mutation, occurring in about 20%–30% of cutaneous melanomas, leads to activation of RAS‐RAF‐MAPK cascade and represents a clear distinct clinicopathological entity in melanoma. In contrast with BRAF mutant melanoma, no specific target therapies are available outside the setting of clinical trials. In the field of immunoncology, the predictive role of NRAS mutation with respect to checkpoint inhibitors treatment has not clearly established and deserves further investigation. At present, the standard treatment is the same as for BRAF wild type melanoma. Ongoing trials are exploring novel combination strategies among patients with advanced NRAS mutant melanoma.     

A five‐microRNA signature for individualized prognosis evaluation and radiotherapy guidance in patients with diffuse lower‐grade glioma

Despite the prognostic value of IDH and other gene mutations found in diffuse glioma, markers that judge individual prognosis of patients with diffuse lower‐grade glioma (LGG) are still lacking. This study aims to develop an expression‐based microRNA signature to provide survival and radiotherapeutic response prediction for LGG patients. MicroRNA expression profiles and relevant clinical information of LGG patients were downloaded from The Cancer Genome Atlas (TCGA; the training group) and the Chinese Glioma Genome Atlas (CGGA; the test group). Cox regression analysis, random survival forests‐variable hunting (RSFVH) screening and receiver operating characteristic (ROC) were used to identify the prognostic microRNA signature. ROC and TimeROC curves were plotted to compare the predictive ability of IDH mutation and the signature. Stratification analysis was conducted in patients with radiotherapy information. Gene ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were performed to explore the biological function of the signature. We identified a five‐microRNA signature that can classify patients into low‐risk or high‐risk group with significantly different survival in the training and test datasets (P < 0.001). The five‐microRNA signature was proved to be superior to IDH mutation in survival prediction (AUCtraining = 0.688 vs 0.607). Stratification analysis found the signature could further divide patients after radiotherapy into two risk groups. GO and KEGG analyses revealed that microRNAs from the prognostic signature were mainly enriched in cancer‐associated pathways. The newly discovered five‐microRNA signature could predict survival and radiotherapeutic response of LGG patients based on individual microRNA expression.     

RB-pathway disruption in breast cancer: Differential association with disease subtypes,disease-specific prognosis and therapeutic response

In breast cancer, inactivation of the RB tumor suppressor gene is believed to occur via multiple mechanisms to facilitate tumorigenesis. However, the prognostic and predictive value of RB status in disease-specific clinical outcomes has remained uncertain. We investigated RB pathway deregulation in the context of both ER-positive and ER-negative disease using combined microarray datasets encompassing over 900 breast cancer patient samples. Disease-specific characteristics of RB pathway deregulation were investigated in this dataset by evaluating correlation among pathway genes as well as differential expression across patient tumor populations defined by ER status. Survival analysis among these breast cancer samples demonstrates that the RB-loss signature is associated with poor disease outcome within several independent cohorts. Within the ER-negative subpopulation, the RB-loss signature is associated with improved response to chemotherapy and longer relapse-free survival. Additionally, while individual genes in the RB target signature closely reproduce its prognostic value, they also serve to predict and monitor response to therapeutic compounds, such as the cytostatic agent PD-0332991. These results indicate that the RB-loss signature expression is associated with poor outcome in breast cancer, but predicts improved response to chemotherapy based on data in ER-negative populations. While the RB-loss signature, as a whole, demonstrates prognostic and predictive utility, a small subset of markers could be sufficient to stratify patients based on RB function and inform the selection of appropriate therapeutic regimens.Key words: RB, breast cancer, microarray, proliferation, cytostatics