Differences in global gene expression in melanoma cell lines with and without homozygous deletion of the CDKN2A locus genes

JUN dependency in distinct early and late BRAF inhibition adaptation states of melanoma

A prominent mechanism of acquired resistance to BRAF inhibitors in BRAF (V600) -mutant melanoma is associated with the upregulation of receptor tyrosine kinases. Evidences suggested that this resistance mechanism is part of a more complex cellular adaptation process. Using an integrative strategy, we found this mechanism to invoke extensive transcriptomic, (phospho-) proteomic and phenotypic alterations that accompany a cellular transition to a de-differentiated, mesenchymal and invasive state. Even short-term BRAF-inhibitor exposure leads to an early adaptive, differentiation state change-characterized by a slow-cycling, persistent state. The early persistent state is distinct from the late proliferative, resistant state. However, both differentiation states share common signaling alterations including JUN upregulation. Motivated by the similarities, we found that co-targeting of BRAF and JUN is synergistic in killing fully resistant cells; and when used up-front, co-targeting substantially impairs the formation of the persistent subpopulation. We confirmed that JUN upregulation is a common response to BRAF inhibitor treatment in clinically treated patient tumors. Our findings demonstrate that events shared between early- and late-adaptation states provide candidate up-front co-treatment targets.

GON4L Drives Cancer Growth through a YY1-Androgen Receptor-CD24 Axis

In principle, the inhibition of candidate gain-of-function genes defined through genomic analyses of large patient cohorts offers an attractive therapeutic strategy. In this study, we focused on changes in expression of CD24, a well-validated clinical biomarker of poor prognosis and a driver of tumor growth and metastasis, as a benchmark to assess functional relevance. Through this approach, we identified GON4L as a regulator of CD24 from screening a pooled shRNA library of 176 candidate gain-of-function genes. GON4L depletion reduced CD24 expression in human bladder cancer cells and blocked cell proliferation in vitro and tumor xenograft growth in vivo Mechanistically, GON4L interacted with transcription factor YY1, promoting its association with the androgen receptor to drive CD24 expression and cell growth. In clinical bladder cancer specimens, expression of GON4L, YY1, and CD24 was elevated compared with normal bladder urothelium. This pathway is biologically relevant in other cancer types as well, where CD24 and the androgen receptor are clinically prognostic, given that silencing of GON4L and YY1 suppressed CD24 expression and growth of human lung, prostate, and breast cancer cells. Overall, our results define GON4L as a novel driver of cancer growth, offering new biomarker and therapeutic opportunities. Cancer Res; 76(17); 5175-85. ©2016 AACR.

Immunolabeling for p16, WT1, and Fli-1 in the assignment of growth phase for cutaneous melanomas

Distinction between radial growth phase (RGP) and vertical growth phase (VGP) in cutaneous melanomas is prognostically significant. Despite established morphological criteria, molecular markers to separate RGP and VGP have not been well established. The goal of this study was to investigate associations of p16, WT1, and Fli-1 with RGP-to-VGP progression, by immunohistochemistry. The p16 is a tumor suppressor, whereas WT1 and Fli-1 are transcriptional activators. The authors hypothesized that entry into VGP would be associated with decreased p16 and increased WT1 and Fli-1. Paraffin sections from 18 RGP and 15 VGP melanomas were immunostained with well-characterized antibodies to p16, WT1, and Fli-1. Melanoma growth phases were determined using precodified morphological attributes. In RGP melanomas, p16 was expressed in 15 of 18 (83%), WT1 in 17 of 17 (100%), and Fli-1 at least focally in 6 of 18 (33%). The deep dermal component of VGP melanomas stained positively for Fli-1 in 9 of 14 (64%), strongly for WT1 in 10 of 14 (71%), and strongly for p16 in only 2 of 15 (13%). Observed patterns of WT1 immunopositivity did not support the authors' hypothesis; it is not likely to be a good indicator of VGP. On the other hand, Fli-1 staining trended toward more positive deep tumor compartment staining and p16 to weaker staining in the deep compartment. At present, application of histological criteria remains the best method for assignment of growth phase in melanomas; however, p16 and possibly Fli-1 immunostains may serve as useful adjuncts in morphologically indeterminate cases.

Modulation of NKp30- and NKp46-mediated natural killer cell responses by poxviral hemagglutinin

Natural killer (NK) cells are an important element in the immune defense against the orthopox family members vaccinia virus (VV) and ectromelia virus (ECTV). NK cells are regulated through inhibitory and activating signaling receptors, the latter involving NKG2D and the natural cytotoxicity receptors (NCR), NKp46, NKp44 and NKp30. Here we report that VV infection results in an upregulation of ligand structures for NKp30 and NKp46 on infected cells, whereas the binding of NKp44 and NKG2D was not significantly affected. Likewise, infection with ectromelia virus (ECTV), the mousepox agent, enhanced binding of NKp30 and, to a lesser extent, NKp46. The hemagglutinin (HA) molecules from VV and ECTV, which are known virulence factors, were identified as novel ligands for NKp30 and NKp46. Using NK cells with selectively silenced NCR expression and NCR-CD3ζ reporter cells, we observed that HA present on the surface of VV-infected cells, or in the form of recombinant soluble protein, was able to block NKp30-triggered activation, whereas it stimulated the activation through NKp46. The net effect of this complex influence on NK cell activity resulted in a decreased NK lysis susceptibility of infected cells at late time points of VV infection when HA was expression was pronounced. We conclude that poxviral HA represents a conserved ligand of NCR, exerting a novel immune escape mechanism through its blocking effect on NKp30-mediated activation at a late stage of infection.

Expression and Function of Kruppel Like-Factors (KLF) in Carcinogenesis

Krüppel-like factor (KLF) family members share a three C2H2 zinc finger DNA binding domain, and are involved in cell proliferation and differentiation control in normal as in pathological situations. Studies over the past several years support a significant role for this family of transcription factors in carcinogenesis. KLFs can both activate and repress genes that participate in cell-cycle regulation. Among them, many up-regulated genes are inhibitors of proliferation, whereas genes that promote cell proliferation are repressed. However, several studies do present KLFs as positive regulator of cell proliferation. KLFs can be deregulated in multiple cancers either by loss of heterozygosity (LOH), somatic mutation or transcriptional silencing by promoter hypermethylation. Accordingly, KLF expression was shown to mediate growth inhibition when ectopically expressed in multiple cancer-derived cell lines through the inhibition of a number of key oncogenic signaling pathways, and to revert the tumorogenic phenotype in vivo. Taken together, these observations suggest that KLFs act as tumor suppressor. However, in some occasion, KLFs could act as tumor promoters, depending on “cellular context”. Thus, this review will discuss the roles and the functions of KLF family members in carcinogenesis, with a special focus on cancers from epithelial origin.

Somatic alterations in the melanoma genome: a high-resolution array-based comparative genomic hybridization study

We performed DNA microarray-based comparative genomic hybridization to identify somatic alterations specific to melanoma genome in 60 human cell lines from metastasized melanoma and from 44 corresponding peripheral blood mononuclear cells. Our data showed gross but nonrandom somatic changes specific to the tumor genome. Although the CDKN2A (78%) and PTEN (70%) loci were the major targets of mono-allelic and bi-allelic deletions, amplifications affected loci with BRAF (53%) and NRAS (12%) as well as EGFR (52%), MITF (40%), NOTCH2 (35%), CCND1 (18%), MDM2 (18%), CCNE1 (10%), and CDK4 (8%). The amplified loci carried additional genes, many of which could potentially play a role in melanoma. Distinct patterns of copy number changes showed that alterations in CDKN2A tended to be more clustered in cell lines with mutations in the BRAF and NRAS genes; the PTEN locus was targeted mainly in conjunction with BRAF mutations. Amplification of CCND1, CDK4, and other loci was significantly increased in cell lines without BRAF-NRAS mutations and so was the loss of chromosome arms 13q and 16q. Our data suggest involvement of distinct genetic pathways that are driven either through oncogenic BRAF and NRAS mutations complemented by aberrations in the CDKN2A and PTEN genes or involve amplification of oncogenic genomic loci and loss of 13q and 16q. It also emerges that each tumor besides being affected by major and most common somatic genetic alterations also acquires additional genetic alterations that could be crucial in determining response to small molecular inhibitors that are being currently pursued.

Activation of natural killer cells by newcastle disease virus hemagglutinin-neuraminidase

The avian paramyxovirus Newcastle disease virus (NDV) selectively replicates in tumor cells and is known to stimulate T-cell-, macrophage-, and NK cell-mediated responses. The mechanisms of NK cell activation by NDV are poorly understood so far. We studied the expression of ligand structures for activating NK cell receptors on NDV-infected tumor cells. Upon infection with the nonlytic NDV strain Ulster and the lytic strain MTH-68/H, human carcinoma and melanoma cells showed enhanced expression of ligands for the natural cytotoxicity receptors NKp44 and NKp46, but not NKp30. Ligands for the activating receptor NKG2D were partially downregulated. Soluble NKp44-Fc and NKp46-Fc, but not NKp30-Fc, chimeric proteins bound specifically to NDV-infected tumor cells and to NDV particle-coated plates. Hemagglutinin-neuraminidase (HN) of the virus serves as a ligand structure for NKp44 and NKp46, as indicated by the blockade of binding to NDV-infected cells and viral particles in the presence of anti-HN antibodies and by binding to cells transfected with HN cDNA. Consistent with the recognition of sialic acid moieties by the viral lectin HN, the binding of NKp44-Fc and NKp46-Fc was lost after desialylation. NKp44- and NKp46-CD3zeta lacZ-inducible reporter cells were activated by NDV-infected cells. NDV-infected tumor cells stimulated NK cells to produce increased amounts of the effector lymphokines gamma interferon and tumor necrosis factor alpha. Primary NK cells and the NK line NK-92 lysed NDV-infected tumor cells with enhanced efficiency, an effect that was eliminated by the treatment of target cells with the neuraminidase inhibitor Neu5Ac2en. These results suggest that direct activation of NK cells contributes to the antitumor effects of NDV.

Identification of ARHGEF17, DENND2D, FGFR3, and RB1 mutations in melanoma by inhibition of nonsense-mediated mRNA decay

Gene identification by nonsense-mediated mRNA decay inhibition (GINI) has proven to be a strategy for genome-wide discovery of genes containing inactivating mutations in colon and prostate cancers. Here, we present the first study of inhibition of the nonsense-mediated mRNA decay (NMD) pathway in melanoma. We used a combination of emetine and actinomycin D treatment to stabilize mRNAs containing premature termination codons (PTCs), followed by microarray analysis and sequencing to identify novel tumor suppressor genes (TSGs) in a panel of 12 melanoma cell lines. Stringent analysis of the array data was used to select 35 candidate genes for sequencing. Of these, 4 (11%) were found to carry PTCs, including ARHGEF17, DENND2D, FGFR3, and RB1. While RB1 mutations have previously been described in melanoma, the other three genes represent potentially novel melanoma; TSGs. ARHGEF17 showed a G1865A mutation leading to W622X in a cell line derived from a mucosal melanoma; in RB1 a C1411T base change resulting in Q471X was discovered in a cell line derived from an acral melanoma; and the FGFR3 and DENND2D genes had intronic insertions leading to PTCs in cell lines derived from superficially spreading melanomas. We conclude that although the false positive rate is high, most likely due to the lack of DNA mismatch repair gene defects, the GINI protocol is one approach to discover novel TSGs in melanoma.

Kruppel-like factor 4 induces p27Kip1 expression in and suppresses the growth and metastasis of human pancreatic cancer cells

The zinc finger transcription factor Krüppel-like factor 4 (KLF4) has been implicated in both tumor suppression and progression. However, its function in pancreatic cancer has not been well characterized. Here, we show that pancreatic cancer cell lines expressed various levels of KLF4 RNA and protein. Ectopic expression of KLF4 by FG and BxPC-3 pancreatic cancer cells resulted in cell cycle arrest and marked inhibition of cell growth in vitro and attenuation of tumor growth and metastasis in an orthotopic mouse model. Overexpression of KLF4 also led to significant induction of p27(Kip1) expression, at both the RNA and protein levels, in a dose- and time-dependent manner, indicating that KLF4 transcriptionally regulates the expression of p27(Kip1). Chromatin immunoprecipitation assays consistently showed that KLF4 protein physically interacts with the p27(Kip1) promoter. Promoter deletion and point mutation analyses indicated that a region between nucleotides -435 and -60 of the p27(Kip1) promoter and intact of the three KLF4-binding sites within that region were required for the full induction of p27(Kip1) promoter activity by KLF4. Our findings suggest that KLF4 transactivates p27(Kip1) expression and inhibits the growth and metastasis of human pancreatic cancer.

Gene expression signature associated with BRAF mutations in human primary cutaneous melanomas

With the aim to correlate BRAF mutation status with gene expression in human primary cutaneous melanomas, and thus to get more insight on the consequences of BRAF mutation on cell biology, we analyzed all expression data obtained in melanomas from which DNA was extracted from the same tissue slides that were used for the expression study. A cohort of 69 frozen primary melanoma whose oligonucleotide micro-array expression data were available, were genotyped for BRAF and NRAS genes. The expression data from these melanomas were re-analyzed according to BRAF mutational status. A set of 250 probes representing 209 genes that were significantly (raw P< or =0.001) associated with BRAF mutation status was identified and 17 of these were previously shown to be implicated in cutaneous melanoma progression or pigmentation pathway-associated genes driven by the microphthalmia transcription factor (MITF). The list of 34 top probes contained no more than 1% of false discoveries with a probability of 0.95. Among the genes that differentiated most strongly between BRAF mutated and non-mutated melanomas, there were those involved in melanoma immune response such as MAGE-D2, CD63, and HSP70. These findings support the immunogenicity of BRAF(V600E), eliciting patients T-cell responses in various in vitro assays. The genes whose expression is associated with BRAF mutations are not simply restricted to the MAPK/ERK signaling but also converge to enhanced immune responsiveness, cell motility and melanosomes processing involved in the adaptative UV response.

Gene expression profiling in melanoma identifies novel downstream effectors of p14ARF

p14ARF is inactivated by deletions/mutations in many cancer types and can suppress cell growth by both p53-dependent and p53-independent mechanisms. To identify novel downstream effectors of p14ARF, we used gene expression profiling as a primary screening tool to select candidates for follow up validation studies using in vitro cell-based assays. Gene expression profiles of a panel of 35 melanoma cell lines with either wild-type (n = 12) or mutant (n = 23) p14ARF were compared to identify genes associated with inactivation of p14ARF. Analysis of the microarray data identified 1,316 probe sets that were significantly (p < 0.01) differentially expressed between the p14ARF wild-type and mutant cell lines. Pathway analysis of these genes showed an overrepresentation of many receptor-mediated signal transduction pathways, e.g. TGFbeta, EGF, HGF, PDGF, MAPK, Wnt and integrin pathways. A number of components of these pathways, including FLRT3, RUNX2, MIG-6 and SMURF2 were confirmed as downstream targets of p14ARF using p14ARF-inducible cell lines and RNAi. We propose that regulation of these genes may contribute to melanoma development when p14ARF function is lost.

Broad spectrum identification of SUMO substrates in melanoma cells

Like phosphorylation, protein sumoylation likely represents a dynamic PTM to alter protein function in support of cell regulatory systems. The broad-spectrum impact of transient or chronic engagement of signal transduction cascades on protein sumoylation has not been explored. Here, we find that epidermal growth factor (EGF) stimulation evokes a rapid alteration in small ubiquitin modifier (SUMO) target selection, while oncogene expression alters steady-state SUMO-protein profiles. A proteomic SUMO target analysis in melanoma cells identified proteins involved in cellular signaling, growth control, and neural differentiation.

Confirmation of a BRAF mutation-associated gene expression signature in melanoma

Mutations in the BRAF oncogene occur in the majority of melanomas, leading to the activation of the mitogen-activated protein kinase pathway and the transcription of downstream effectors. As BRAF and its effectors could be good melanoma therapy targets, defining the repertoire of genes that are differentially regulated because of BRAF mutational activation is an important objective. Towards this goal, we and others have attempted to determine whether a BRAF mutation-associated gene expression profile exists. Results have been mixed, with some groups reporting a BRAF-signature and another group not. Here we resolve this issue and confirm that while gene-by-gene correlations fail to reveal a specific gene(s) whose expression correlates with BRAF status, a BRAF signature can be distinguished by analysis of global expression patterns. Specifically, we have here applied support vector machine (SVM) analysis to Affymetrix microarray data from a panel of 63 melanoma cell lines. SVMs found a BRAF signature in training samples and predicted BRAF mutation status with high accuracy (AUC=0.840) in the remaining samples. We verified this is a generalized BRAF signature by repeating the analysis in three published microarray datasets, and again found that SVMs predicted BRAF mutation well (Philadelphia: AUC=0.788; Zurich: AUC=0.688; Mannheim: AUC=0.686). An ensemble of 300 SVMs trained on our data also predicted BRAF mutation status in two of the three published datasets (Philadelphia AUC=0.778; Zurich AUC=0.719; Mannheim AUC=0.564). Taken together, these data support the existence of a BRAF mutation-specific expression signature.