Elaine Cheng

Expression Profiling Reveals Novel Pathways in the Transformation of Melanocytes to Melanomas

Affymetrix and spotted oligonucleotide microarrays were used to assess global differential gene expression comparing normal human melanocytes with six independent melanoma cell strains from advanced lesions. The data, validated at the protein level for selected genes, confirmed the overexpression in melanoma cells relative to normal melanocytes of several genes in the growth factor/receptor family that confer growth advantage and metastasis. In addition, novel pathways and patterns of associated expression in melanoma cells not reported before emerged, including the following: (a) activation of the NOTCH pathway; (b) increased Twist expression and altered expression of additional transcriptional regulators implicated in embryonic development and epidermal/mesenchymal transition; (c) coordinated activation of cancer/testis antigens; (d) coordinated down-regulation of several immune modulation genes, in particular in the IFN pathways; (e) down-regulation of several genes implicated in membrane trafficking events; and (f) down-regulation of growth suppressors, such as the Prader-Willi gene NECDIN, whose function was confirmed by overexpression of ectopic Flag-necdin. Validation of differential expression using melanoma tissue microarrays showed that reduced ubiquitin COOH-terminal esterase L1 in primary melanoma is associated with worse outcome and that increased expression of the basic helix-loop-helix protein Twist is associated with worse outcome. Some differentially expressed genes reside on chromosomal regions displaying common loss or gain in melanomas or are known to be regulated by CpG promoter methylation. These results provide a comprehensive view of changes in advanced melanoma relative to normal melanocytes and reveal new targets that can be used in assessing prognosis, staging, and therapy of melanoma patients.

PLX4032, a selective BRAFV600E kinase inhibitor, activates the ERK pathway and enhances cell migration and proliferation of BRAFWT melanoma cells

BRAFV600E/K is a frequent mutationally active tumor‐specific kinase in melanomas that is currently targeted for therapy by the specific inhibitor PLX4032. Our studies with melanoma tumor cells that are BRAFV600E/K and BRAFWT showed that, paradoxically, while PLX4032 inhibited ERK1/2 in the highly sensitive BRAFV600E/K, it activated the pathway in the resistant BRAFWT cells, via RAF1 activation, regardless of the status of mutations in NRAS or PTEN. The persistently active ERK1/2 triggered downstream effectors in BRAFWT melanoma cells and induced changes in the expression of a wide‐spectrum of genes associated with cell cycle control. Furthermore, PLX4032 increased the rate of proliferation of growth factor‐dependent NRAS Q61L mutant primary melanoma cells, reduced cell adherence and increased mobility of cells from advanced lesions. The results suggest that the drug can confer an advantage to BRAFWT primary and metastatic tumor cells in vivo and provide markers for monitoring clinical responses.

Incidence of the V600K mutation among melanoma patients with BRAF mutations, and potential therapeutic response to the specific BRAF inhibitor PLX4032

Activating mutations in BRAF kinase are common in melanomas. Clinical trials with PLX4032, the mutant-BRAF inhibitor, show promising preliminary results in patients selected for the presence of V600E mutation. However, activating V600K mutation is the other most common mutation, yet patients with this variant are currently excluded from the PLX4032 trials. Here we present evidence that a patient bearing the BRAF V600K mutation responded remarkably to PLX4032, suggesting that clinical trials should include all patients with activating BRAF V600E/K mutations.

Genome-wide screen of promoter methylation identifies novel markers in melanoma

DNA methylation is an important component of epigenetic modifications, which influences the transcriptional machinery aberrant in many human diseases. In this study we present the first genome-wide integrative analysis of promoter methylation and gene expression for the identification of methylation markers in melanoma. Genome-wide promoter methylation and gene expression of eight early-passage human melanoma cell strains were compared with newborn and adult melanocytes. We used linear mixed effect models (LME) in combination with a series of filters based on the localization of promoter methylation relative to the transcription start site, overall promoter CpG content, and differential gene expression to discover DNA methylation markers. This approach identified 76 markers, of which 68 were hyper- and eight hypomethylated (LME, P < 0.05). Promoter methylation and differential gene expression of five markers (COL1A2, NPM2, HSPB6, DDIT4L, MT1G) were validated by sequencing of bisulfite-modified DNA and real-time reverse transcriptase PCR, respectively. Importantly, the incidence of promoter methylation of the validated markers increased moderately in early and significantly in advanced-stage melanomas, using early-passage cell strains and snap-frozen tissues (n = 18 and n = 24, respectively) compared with normal melanocytes and nevi (n = 11 and n = 9, respectively). Our approach allows robust identification of methylation markers that can be applied to other studies involving genome-wide promoter methylation. In conclusion, this study represents the first unbiased systematic effort to determine methylation markers in melanoma and revealed several novel genes regulated by promoter methylation that were not described in cancer cells before.

Exome sequencing identifies recurrent mutations in NF1 and RASopathy genes in sun-exposed melanomas

We report on whole-exome sequencing (WES) of 213 melanomas. Our analysis established NF1, encoding a negative regulator of RAS, as the third most frequently mutated gene in melanoma, after BRAF and NRAS. Inactivating NF1 mutations were present in 46% of melanomas expressing wild-type BRAF and RAS, occurred in older patients and showed a distinct pattern of co-mutation with other RASopathy genes, particularly RASA2. Functional studies showed that NF1 suppression led to increased RAS activation in most, but not all, melanoma cases. In addition, loss of NF1 did not predict sensitivity to MEK or ERK inhibitors. The rebound pathway, as seen by the induction of phosphorylated MEK, occurred in cells both sensitive and resistant to the studied drugs. We conclude that NF1 is a key tumor suppressor lost in melanomas, and that concurrent RASopathy gene mutations may enhance its role in melanomagenesis.

Abnormal Acidification of Melanoma Cells Induces Tyrosinase Retention in the Early Secretory Pathway

In tyrosinase-positive amelanotic melanoma cells, inactive tyrosinase accumulates in the endoplasmic reticulum. Based on studies described here, we propose that aberrant vacuolar proton ATPase (V-ATPase)-mediated proton transport in melanoma cells disrupts tyrosinase trafficking through the secretory pathway. Amelanotic but not melanotic melanoma cells or normal melanocytes display elevated proton export as observed by the acidification of the extracellular medium and their ability to maintain neutral intracellular pH. Tyrosinase activity and transit through the Golgi were restored by either maintaining the melanoma cells in alkaline medium (pH 7.4–7.7) or by restricting glucose uptake. The translocation of tyrosinase out of the endoplasmic reticulum and the induction of cell pigmentation in the presence of the ionophore monensin or the specific V-ATPase inhibitors concanamycin A and bafilomycin A1 supported a role for V-ATPases in this process. Because it was previously shown that V-ATPase activity is increased in solid tumors in response to an acidified environment, the appearance of hypopigmented cells in tyrosinase-positive melanoma tumors may indicate the onset of enhanced glycolysis and extracellular acidification, conditions known to favor metastatic spread and resistance to weak base chemotherapeutic drugs.

Proper Folding and Endoplasmic Reticulum to Golgi Transport of Tyrosinase Are Induced by Its Substrates, DOPA and Tyrosine

Tyrosinase is essential for pigmentation and is a source of tumor-derived antigenic peptides and cellular immune response. Wild type tyrosinase in melanoma cells and certain albino mutants in untransformed melanocytes are targeted to proteolytic degradation by the 26 S proteasome due to retention of the misfolded protein in the endoplasmic reticulum and its subsequent retranslocation to the cytosol. Here, we demonstrate that the substrates DOPA and tyrosine induced in melanoma cells a transition of misfolded wild type tyrosinase to the native form that is resistant to proteolysis, competent to exit the endoplasmic reticulum, and able to produce melanin. Because the enzymatic activity of tyrosinase is induced by DOPA, we propose that proper folding of the wild type protein, just like mutant forms, is tightly linked to its catalytic state. Loss of pigmentation, therefore, in tyrosinase-positive melanoma cells is a consequence of tumor-induced metabolic changes that suppress tyrosinase activity and DOPA production within these cells.

Integrative Analysis of Epigenetic Modulation in Melanoma Cell Response to Decitabine: Clinical Implications

Decitabine, an epigenetic modifier that reactivates genes otherwise suppressed by DNA promoter methylation, is effective for some, but not all cancer patients, especially those with solid tumors. It is commonly recognized that to overcome resistance and improve outcome, treatment should be guided by tumor biology, which includes genotype, epigenotype, and gene expression profile. We therefore took an integrative approach to better understand melanoma cell response to clinically relevant dose of decitabine and identify complementary targets for combined therapy. We employed eight different melanoma cell strains, determined their growth, apoptotic and DNA damage responses to increasing doses of decitabine, and chose a low, clinically relevant drug dose to perform whole-genome differential gene expression, bioinformatic analysis, and protein validation studies. The data ruled out the DNA damage response, demonstrated the involvement of p21Cip1 in a p53-independent manner, identified the TGFβ pathway genes CLU and TGFBI as markers of sensitivity to decitabine and revealed an effect on histone modification as part of decitabine-induced gene expression. Mutation analysis and knockdown by siRNA implicated activated β-catenin/MITF, but not BRAF, NRAS or PTEN mutations as a source for resistance. The importance of protein stability predicted from the results was validated by the synergistic effect of Bortezomib, a proteasome inhibitor, in enhancing the growth arrest of decitabine in otherwise resistant melanoma cells. Our integrative analysis show that improved therapy can be achieved by comprehensive analysis of cancer cells, identified biomarkers for patients selection and monitoring response, as well as targets for improved combination 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 BRAFV600K 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 (BRAFL505H), 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 BRAFL505H, 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.

Detecting copy number status and uncovering subclonal markers in heterogeneous tumor biopsies

BackgroundGenomic aberrations can be used to determine cancer diagnosis and prognosis. Clinically relevant novel aberrations can be discovered using high-throughput assays such as Single Nucleotide Polymorphism (SNP) arrays and next-generation sequencing, which typically provide aggregate signals of many cells at once. However, heterogeneity of tumor subclones dramatically complicates the task of detecting aberrations.ResultsThe aggregate signal of a population of subclones can be described as a linear system of equations. We employed a measure of allelic imbalance and total amount of DNA to characterize each locus by the copy number status (gain, loss or neither) of the strongest subclonal component. We designed simulated data to compare our measure to existing approaches and we analyzed SNP-arrays from 30 melanoma samples and transcriptome sequencing (RNA-Seq) from one melanoma sample.We showed that any system describing aggregate subclonal signals is underdetermined, leading to non-unique solutions for the exact copy number profile of subclones. For this reason, our illustrative measure was more robust than existing Hidden Markov Model (HMM) based tools in inferring the aberration status, as indicated by tests on simulated data. This higher robustness contributed in identifying numerous aberrations in several loci of melanoma samples. We validated the heterogeneity and aberration status within single biopsies by fluorescent in situ hybridization of four affected and transcriptionally up-regulated genes E2F8, ETV4, EZH2 and FAM84B in 11 melanoma cell lines. Heterogeneity was further demonstrated in the analysis of allelic imbalance changes along single exons from melanoma RNA-Seq.ConclusionsThese studies demonstrate how subclonal heterogeneity, prevalent in tumor samples, is reflected in aggregate signals measured by high-throughput techniques. Our proposed approach yields high robustness in detecting copy number alterations using high-throughput technologies and has the potential to identify specific subclonal markers from next-generation sequencing data.