- Open Access
MicroRNA expression profiling predicts clinical outcome of carboplatin/paclitaxel-based therapy in metastatic melanoma treated on the ECOG-ACRIN trial E2603
Clinical Epigeneticsvolume 7, Article number: 58 (2015)
Carboplatin/paclitaxel (CP), with or without sorafenib, result in objective response rates of 18–20 % in unselected chemotherapy-naïve patients. Molecular predictors of survival and response to CP-based chemotherapy in metastatic melanoma (MM) are critical to improving the therapeutic index.
Intergroup trial E2603 randomized MM patients to CP with or without sorafenib. Expression data were collected from pre-treatment formalin-fixed paraffin-embedded (FFPE) tumor tissues from 115 of 823 patients enrolled on E2603. The selected patients were balanced across treatment arms, BRAF status, and clinical outcome. We generated data using Nanostring array (microRNA (miRNA) expression) and DNA-mediated annealing, selection, extension and ligation (DASL)/Illumina microarrays (HT12 v4) (mRNA expression) with protocols optimized for FFPE samples. Integrative computational analysis was performed using a novel Tree-guided Recursive Cluster Selection (T-ReCS)  algorithm to select the most informative features/genes, followed by TargetScan miRNA target prediction (Human v6.2) and mirConnX  for network inference.
T-ReCS identified PLXNB1 as negatively associated with progression-free survival (PFS) and miR-659-3p as the primary miRNA associated positively with PFS. miR-659-3p was differentially expressed based on PFS but not based on treatment arm, BRAF or NRAS status. Dichotomized by median PFS (less vs greater than 4 months), miR-659-3p expression was significantly different. High miR-659-3p expression distinguished patients with responsive disease (complete or partial response) from patients with stable disease. miR-659-3p predicted gene targets include NFIX, which is a transcription factor known to interact with c-Jun and AP-1 in the context of developmental processes and disease.
This novel integrative analysis implicates miR-659-3p as a candidate predictive biomarker for MM patients treated with platinum-based chemotherapy and may serve to improve patient selection.
The incidence of melanoma has increased progressively over the last 30 years with an estimated 76,100 new cases in 2014 and a projected 9700 deaths (http://seer.cancer.gov/statfacts/html/melan.html). The prognosis of metastatic disease is historically poor, with 5-year survival rates of less than 6 % and median survival of 6 to 9 months [3, 4]. Until 2010, alkylating agents were the mainstay of treatment but associated with response rates less than 10–15 % [5, 6]. Platinum-based chemotherapy with carboplatin and paclitaxel (CP) is associated with response rates of 18–20 % in chemotherapy-naïve patients and 11–12 % in the second-line setting [7, 8]. Significant therapeutic advances have been made in recent years with FDA approval of the monoclonal antibody to CTLA-4, ipilimumab, which results in durable responses and improved overall survival (OS) in patients with metastatic melanoma (MM) compared with peptide vaccine or dacarbazine . Antibody therapies against PD-1 and PD-L1 have also been associated with durable responses in MM, and the first agent targeting this pathway was granted FDA approval in September 2014 [10–12]. Finally, the selective targeting of BRAF activating mutations, which occur in about half of melanomas, with vemurafenib, dabrafenib, and trametinib, is associated with marked clinical activity and improved OS compared with chemotherapy [13, 14].
Despite the recent successes of immunotherapy and targeted therapy, chemotherapy continues to be a significant modality utilized worldwide in MM, especially after progression on immunotherapy and/or antitumor targeted agents. Recalcitrance of melanoma to chemotherapy is a significant issue due to limited response rates and lack of predictive biomarkers of chemotherapy benefit. Identification of biomarkers of response and survival for patients treated with chemotherapy is critical to allow more refined application of this modality and improving treatment outcomes. MicroRNAs (miRNAs) are endogenous non-coding RNAs approximately 22 nucleotides in length that regulate gene expression at the post-transcriptional level and play an important role in the regulation of tumor suppressor genes, oncogenes, and genes involved in the epigenetic machinery and cancer cell metastasis . miR-15a, miR-16, and the let-7 family function predominately as tumor suppressors, while the miR-17-92 cluster functions as oncogenes [16–20]. Differential expression of four miRNAs (miR-205, miR-23b, miR-146a, and miR-155) distinguishes melanoma cell lines from nevi and primary or metastatic melanoma . Differential miRNA expression also distinguishes among common melanoma histologic subtypes (superficial spreading melanoma vs nodular melanoma) and benign nevi compared with either primary or metastatic melanoma [21–24]. Of greatest importance, miRNAs have been shown to have prognostic implications in melanoma. Downregulation of miR-125b in patients with stage T2 primary cutaneous melanoma is associated with the presence of micrometastatic disease to the sentinel lymph node . Analysis of primary and metastatic melanoma specimens shows that overexpression of an 18-miRNA signature is significantly correlated with longer survival (defined as a post-recurrence survival of more than 18 months) , indicating the value of miRNAs as candidate prognostic biomarkers. To date, there have been no studies examining miRNAs as candidate predictive biomarkers of response to chemotherapy in metastatic melanoma.
The Eastern Cooperative Oncology Group (ECOG) Study, E2603, was a phase III randomized controlled clinical trial evaluating the addition of the unselective multi-tyrosine kinase inhibitor, sorafenib, to platinum-based chemotherapy in patients with previously untreated MM. This phase III trial was the first to examine targeted therapy in melanoma but failed to meet its primary endpoint of demonstrating an OS benefit in MM patients treated with CP alone vs CP with sorafenib. The response rates and progression-free survival (PFS) were similar in either treatment arm. In the present study, our primary objective was to identify molecular predictors of survival and response to CP-based chemotherapy in the MM patients treated on E2603. We applied a novel method for cluster selection  to miRNA and gene expression data in this well-defined patient population in an effort to identify miRNAs as candidate biomarkers of prognosis and response to platinum-based chemotherapy. We then performed novel integrative analysis to construct a proposed genome-wide regulatory network, which lends insight into the role of our candidate miRNA biomarkers in the landscape of MM signaling. To our knowledge, our study is the first to utilize miRNAs as predictive biomarkers for therapy response in melanoma.
Results and discussion
High-quality miRNA and mRNA expression data were generated from legacy samples from 115 MM patients treated on E2603. The clinicopathologic features for the patients included in this analysis are summarized in Table 1. The median PFS for this cohort of patients was 4.4 months and the median OS was 9.8 months.
PLXNB1 and RAD23 are associated with shortened PFS in metastatic melanoma
Among the single genes selected with Tree-guided Recursive Cluster Selection (T-ReCS), PLXNB1, a known tumor suppressor of melanoma, had the most significant negative association with PFS (p = 0.0001; Fig. 1a) [27, 28]. T-ReCS found gene RAD23 in the same group as PLXNB1. RAD23, a component of the nucleotide excision repair mechanism, was also negatively associated with PFS. In the other cluster, T-ReCS found BCMO1, which encodes a key enzyme in beta-carotene metabolism to vitamin A and is thus important for skin protection. BCMO1 was part of a very significant group variable, which included LOC389936, a shRNA construct with no known association with melanoma.
miR-659-3p is associated with improved PFS in metastatic melanoma
T-ReCS identified miR-659-3p as the primary miRNA associated positively with PFS (p = 0.008). The miR-659-3p cluster also included miR-219-3p and miR-519-5p (Fig. 1b), which had less significant or no association with PFS (p = 0.036 and p = 0.12, respectively). We examined the Cox regression coefficients of this selected cluster to interrogate our results. Indeed, all of these miRNAs were positively associated with PFS. miR-659-3p was not differentially expressed based on the treatment arm (CP with sorafenib vs CP without sorafenib; p = 0.6), BRAF mutation status (mutant vs wild-type; p = 0.34) or NRAS mutation status (mutant vs wild-type; p = 0.55). We dichotomized PFS less than 4 months vs greater than 4 months, which is the median PFS of patients treated with CP with or without sorafenib on E2603. Higher miR-659-3p expression remained prognostic of longer PFS (p = 0.03). Further partitioning of the dataset into terciles or quantiles did not yield any significant results, presumably due to sample size. We also noticed that higher miR-659-3p expression was generally observed in patients who achieved an objective response (complete or partial response) as opposed to patients with stable disease by Response Evaluation Criteria In Solid Tumors (RECIST) (p = 0.04.)
miR-659-3p gene targets
miR-659-3p gene targets were inferred from TargetScan Human v6.2 and from the mirConnX prior information file . The list of predicted targets included genes relevant to the MAP kinase pathway (BRAF and NRAS), genes relevant to the PI3K/MTOR, FGFR, and RAC1 pathways, as well as genes involved in DNA repair (Table 2; Fig. 2b). However, the expression of most of those genes was not significantly anti-correlated with miR-659-3p expression. This might be explained by the fact that each individual miRNA target usually decreases the expression of the host gene by a small amount. However, we did find a significant anti-correlation (p < 0.01) between miR-659-3p and its predicted targets COG6, ENOX2, HCFC2, NFIX, SDC1, and TSC22D2 (Fig. 2a). NFIX (nuclear factor 1 X-type), a member of the NF-1 family, is particularly important since it has been implicated in many developmental and biological processes and diseases, and the role of the NF-1 pathway is emerging as a potential therapeutic target [30–33]. Studies have also shown that it interacts with the oncogenes c-Jun and AP-1 .
miR-331-3p distinguishes OS in metastatic melanoma
Running T-ReCS with OS as the outcome variable, we found miR-331-3p as the primary miRNA positively associated with OS, while its cluster also included miR-93-5p. miR-331-3p was not differentially expressed based on the treatment arm (CP with sorafenib vs CP without sorafenib; p = 0.706), BRAF mutation status (mutant vs wild-type; p = 0.1909) or NRAS mutation status (mutant vs wild-type; p = 0.52).
We have utilized a novel integrative analysis of miRNA and gene expression data to identify molecular biomarkers of PFS and response and found miR-659-3p expression to be positively associated with longer PFS, distinguishing patients with disease that is responsive to CP-based therapy. Characterization of the functional significance of miR-659-3p has been limited to the neuroscience literature, where it has been shown to bind a common genetic variant of and regulate the expression of the gene progranulin, which may significantly increase the risk for frontotemporal dementia . The functional significance of miR-659-3p in cancer and melanoma is as of yet unknown. Of particular significance in our study is the identification of gene targets predicted for miR-659-3p and the construction of a proposed regulatory network which may lend insight into the mechanistic role of miR-659-3p in response to cytotoxic chemotherapy in melanoma.
Gene expression analysis showed PLXNB1 to be associated with foreshortened PFS. PLXNB1 has previously been shown to be lost in MM and in deeply invasive primary tumors and is strongly inhibited by MAP kinase signaling in melanoma cells and melanocytes [27, 28]. In preclinical models, plexin B1 suppressed proliferation, enhanced migration, stimulated Akt activation and, of particular interest, has been reported to render melanoma cell lines resistant to cisplatin-induced apoptosis . Plexin B1 has been shown to inhibit c-Met in response to its ligand hepatocyte growth factor (HGF) and is predicted to be a classic tumor suppressor protein in melanomas in which progression is c-Met dependent . RAD23 was also identified in our model as associated with shorter PFS and is a component of the protein complex that specifically complements the nucleotide excision repair defect of xeroderma pigmentosum group C (XP-c) cell extracts in vitro . Furthermore, RAD23 has been shown to be downregulated in glutathione S-transferase M1 (GSTM1)-null melanomas in patients reporting histories of sunburns . RAD23 has been implicated as a biomarker of response to epigenetic therapy, in particular HDAC inhibition, through a regulatory circuit that involves RAD23B and HDAC6 [38, 39].
Predicted gene targets of miR-659-3p selected in our analysis are genes important to the MAP kinase signaling pathway, in particular BRAF and NRAS, and genes relevant to the PI3K/MTOR, FGFR, and RAC1 pathways, as well as genes involved in DNA repair. BRAF, a member of the RAF kinase family of serine/threonine protein kinases, is mutated in 40–50 % of melanomas and results most frequently from a valine (V) to glutamic acid (E) substitution at residue 600 (BRAF V600E), rendering mutant BRAF protein which no longer requires dimerization for its activity, and is highly predictive of therapeutic response to the specific BRAF inhibitors, vemurafenib and dabrafenib. Mutations in NRAS occur in 15–20 % of melanomas and result in constitutive activation of the NRAS protein and enhanced MAP kinase signaling [40–42]. FGFR1 signaling has been implicated in melanoma progression ; introduction of antisense oligonucleotides targeted toward FGFR1 into metastatic cell lines results in decreased proliferation and signs of differentiation [44, 45]; and injection of an antisense FGFR1 construct into primary and metastatic melanomas grown in nude mice results in inhibition of tumor growth and induction of apoptosis [46, 47]. RAC1, a member of the Rho family of GTPases, is involved in cellular proliferation, survival, motility, and invasiveness . The RAC1 mutation resulting in the substitution of proline for serine at residue 29 is the most frequently mutated amino acid in MM, after V600 mutations in BRAF and the glutamine (Q) substitution at position 61 in NRAS. RAC1 mutations are observed in 4 % of all melanomas and 9 % of sun-exposed melanomas and may occur either independently or in conjunction with BRAF and NRAS mutations, suggesting growth advantages distinct from activation of the MAP kinase signaling pathway [49, 50]. RAC1 and other members of the Rho family of GTPases have been implicated as key modulators of sensitivity to DNA-damaging agents; in particular, chemical inhibition of RAC1 has been shown to sensitize drug-resistant ovarian cancer cell lines to cisplatin [51–53].
While the association between miR-659-3p and these predicted genes targets is intriguing, it is important to note that our gene expression data do not show strong correlations with miR-659-3p and at this point, these findings are not confirmed. These findings do, however, suggest biological relevance of miR-659-3p in the context of known oncogenic drivers in melanoma, an association which requires further investigation in future studies. We did find a significant anti-correlation between miR-659-3p and a small set of genes. NFIX, which has been implicated in many developmental and biological processes and diseases, is one of them and has been shown to interact with known oncogenes. This interaction may represent a clue as to functional relevance of miR-659-3p and requires further validation.
Our study identified an association between miR-659-3p and the clinical activity of CP-based therapy implicating miR-659-3p as a candidate predictive biomarker in patients with metastatic melanoma. This analysis did not include untreated MM patients, and therefore, we are unable to evaluate the prognostic value of miR-659-3p, i.e., we are unable to assess if the association between miR-659-3p and PFS persists in the absence of therapy. These data were generated from formalin-fixed paraffin-embedded (FFPE) legacy samples from a phase III multicenter clinical trial that represents a limitation of this study. Recovery of RNA from archival FFPE tissues can be challenging due to strand breakage and cross-linking by formalin fixation, which may inhibit polymerase in PCR-based molecular assays. We overcame this challenge by utilizing technology optimized for FFPE samples. Specifically, the miRNA expression assay used in this study allows for direct digital detection and counting of miRNAs in a single reaction without the need for an amplification step. Another limitation is that this bio-analysis represents a subset of the 823 patients enrolled on E2603. Archival tissue was not compulsory in E2603, and therefore, genotyping information is available for only 179 patients treated on this trial. It is notable that the patients included in this analysis were similar to the overall E2603 population and the population of patients for whom genotyping information is available [7, 54].
This novel integrative analysis suggests miR-659-3p as a candidate biomarker of response in MM patients treated with platinum-based chemotherapy and may serve to improve patient selection. To our knowledge, our study is the first to identify a miRNA species as a predictor of therapy outcome in metastatic melanoma. While our proposed regulatory network embeds miR-659-3p in the landscape of melanoma signaling, functional studies are paramount to the validation of this miRNA as a mediator of melanoma signaling. The value of the proposed regulatory network is to inform upon pathways which should be interrogated as we continue to elucidate the role of miR-659-3p as a predictive biomarker. Independent validation of the predictive value of miR-659-3p is currently ongoing in an independent cohort of MM patients treated with carboplatin and paclitaxel.
Patients and study design
E2603 enrolled patients with pathologically confirmed advanced unresectable or metastatic melanoma without brain involvement, who had no previous treatment with chemotherapeutic or MAPK pathway targeting agents. The study design has been previously described . Key inclusion criteria were the presence of measurable disease by RECIST version 1.0, at least 18 years of age, ECOG performance status of 0 or 1, and satisfactory baseline organ function. Pre-treatment specimen collection was conducted as part of the clinical trial and analyzed with institutional review board approval. From the 823 patients enrolled on the E2603 clinical trial, 620 archival pre-treatment samples were available for analysis. One hundred and twenty pre-treatment formalin fixed paraffin-embedded (FFPE) tumor tissues from 115 patients were selected to ensure a balanced representation of patients based on treatment arm, BRAF status, and clinical outcome (PFS less than 4 months vs greater than 4 months) were included in this analysis.
Treatment and assessment of response
The patients were randomized in a double-blind fashion to receive CP, either with or without sorafenib. Paclitaxel was administered at 225 mg/m2 over 3 h followed by carboplatin AUC 6 IV over 30 min on day 1 of each 21-day cycle during cycles 1–4. During cycles 5 through 10, paclitaxel was administered at 175 mg/m2 IV over 3 h followed by carboplatin AUC 5 IV over 30 min. Sorafenib 400 mg or placebo was administered twice daily by mouth on days 2 through 19 of each 21-day cycle. Upon completion of 10 cycles of chemotherapy, sorafenib 400 mg twice daily by mouth or placebo was administered daily until disease progression or intolerable toxicity. Tumor response was assessed every 2 cycles during cycles 1 through 10 and then every 3 cycles. Response was defined by RECIST version 1.0.
Genotyping was performed using a custom iPlex Sequenom panel interrogating 74 mutations in 13 genes in the Perelman School of Medicine Molecular Profiling Facility, as previously described .
Total RNA was isolated from melanoma tissues using the PerfectPure RNA Tissue kit (5Prime Inc., MD, USA) at the University of Pittsburgh Cancer Institute (UPCI). RNA was quantified using Ribogreen RNA quantitation Kit (Molecular Probes, Eugene OR). RNA quality was also evaluated by RNA integrity number using the Agilent Bioanalyzer.
miRNA microarray expression profiling and gene expression analysis
Profiling of miRNA levels was performed at the UPCI utilizing technology optimized for FFPE samples: Nanostring technology (Ncounter Human v2 miRNA Expression Assay) for miRNA analysis of 800 human miRNAs was derived from miRBase v.18. Total RNA was used as input for nCounter miRNA sample preparation reactions. All sample preparation and processing were performed according to the manufacturer’s protocol. Following ligation, sample preparation reactions were purified and diluted, and hybridization reactions were performed by incubation at 65 °C for a minimum of 18 h. Hybridized probes were purified and counted on the nCounter Prep Station and Digital Analyzer (NanoString) following the manufacturer’s instructions. For each assay, a high-density scan (600 fields of view) was performed. Whole-genome gene expression analysis was also carried out using DNA-mediated annealing, selection, extension, and ligation (DASL)/Illumina HT-12 v4 Expression BeadChip, a technology optimized for FFPE samples, at the UPCI according to the manufacturer’s protocol.
Statistical analysis for mRNA and miRNA biomarkers predictive of survival
We performed feature selection at the cluster level utilizing a novel algorithm (T-ReCS) . Briefly, starting from a standard feature selection method against a target variable (in our case, PFS or OS), T-ReCS expands the single features into clusters of features according to two statistical criteria that maintain the prediction accuracy. The clusters are defined dynamically as groups of similarly expressed genes. T-ReCS works in two steps. First, it partitions the feature space into a tree structure with the root being the set of all variables (mRNA or miRNA genes, respectively), while each leaf is a single variable. Each internal node in the tree represents a cluster of variables, its descendants, whose values (expression) are similar. The lower in the tree a node is, the more similar the patterns of its members are. For our implementation, we use the tree produced by our Recursive K-means Spectral Clustering (ReKS), which was previously developed for heterogeneous clinical data . The second step of T-ReCS is feature selection operating on the tree structure iteratively. Starting with a set of variables (tree leafs) selected with any feature selection method, it ascends the tree one level at a time and statistically assesses whether the group variables (internal nodes) retain the same predictive properties. For statistical assessment, T-ReCS uses two conditional independence tests with p value thresholds determined by cross-validation. For survival target variables, T-ReCS uses survival MMPC (sMMPC)  for single-variable selection. The association between selected genes and survival was tested using Cox regression (Cox proportional hazards model) .
Regulatory associations between predictive mRNA and miRNA features
A regulatory network was constructed from the mRNA and miRNA genes and their clusters that were selected as characteristic of PFS. The mRNA targets of miRNAs were predicted using TargetScan Human v6.2 , while transcription factor (TF)→gene and TF→miRNA regulatory interactions were obtained from the prior (static) network of mirConnX . mirConnX is a user-friendly web interface for inferring, displaying, and parsing mRNA and miRNA gene regulatory networks. mirConnX infers regulatory networks by combining expression correlation values (dynamic network) with regulatory interactions obtained from the literature and from computational predictions (static network).
Huang GT, Tsamardinos I, Raghu V, Kaminski N, Benos PV. T-ReCS: stable selection of dynamically formed groups of features with application to prediction of clinical outcomes. In: Pacific Symposium in Biocomputing. USA: Big Island, Hawai'i; 2015
Huang GT, Athanassiou C, Benos PV. mirConnX: condition-specific mRNA-microRNA network integrator. Nucleic Acids Res. 2011;39:W416–23.
Cummins D, Cummins JM, Pantle H, Silverman MA, Leonard AL. Cutaneous malignant melanoma. Mayo Clin Proc. 2006;81:500–7.
Barth A, Wanek LA, Morton DL. Prognostic factors in 1,521 melanoma patients with distant metastases. J Am Coll Surg. 1995;181:193–201.
Middleton MR, Grob JJ, Aaronson N, Fierlbeck G, Tilgen W, Seiter S, et al. Randomized phase III study of temozolomide versus dacarbazine in the treatment of patients with advanced metastatic malignant melanoma. J Clin Oncol. 2000;18:158–66.
Patel PM, Suciu S, Mortier L, Kruit WH, Robert C, Schadendorf D, et al. Extended schedule, escalated dose temozolomide versus dacarbazine in stage IV melanoma: final results of a randomised phase III study (EORTC 18032). Eur J Cancer. 2011;47:1476–83.
Flaherty KT LS, Schuchter LM, Flaherty LE, Wright JJ, Leming PD, Kirkwood JM. Final results of E2603: a double-blind, randomized phase III trial comparing carboplatin/paclitaxel with or without sorafenib in metastatic melanoma. J Clin Oncol. 2010;28:Abstract 8511.
Hauschild A, Agarwala SS, Trefzer U, Hogg D, Robert C, Hersey P, et al. Results of a phase III, randomized, placebo-controlled study of sorafenib in combination with carboplatin and paclitaxel as second-line treatment in patients with unresectable stage III or stage IV melanoma. J Clin Oncol. 2009;27:2823–30.
Hodi FS, O'Day SJ, McDermott DF, Weber RW, Sosman JA, Haanen JB, et al. Improved survival with ipilimumab in patients with metastatic melanoma. N Engl J Med. 2010;363:711–23.
Hamid O, Robert C, Daud A, Hodi FS, Hwu WJ, Kefford R, et al. Safety and tumor responses with lambrolizumab (anti-PD-1) in melanoma. N Engl J Med. 2013;369:134–44.
Topalian SL, Hodi FS, Brahmer JR, Gettinger SN, Smith DC, McDermott DF, et al. Safety, activity, and immune correlates of anti-PD-1 antibody in cancer. N Engl J Med. 2012;366:2443–54.
Brahmer JR, Tykodi SS, Chow LQ, Hwu WJ, Topalian SL, Hwu P, et al. Safety and activity of anti-PD-L1 antibody in patients with advanced cancer. N Engl J Med. 2012;366:2455–65.
Chapman PB, Hauschild A, Robert C, Haanen JB, Ascierto P, Larkin J, et al. Improved survival with vemurafenib in melanoma with BRAF V600E mutation. N Engl J Med. 2011;364:2507–16.
Hauschild A, Grob JJ, Demidov LV, Jouary T, Gutzmer R, Millward M, et al. Dabrafenib in BRAF-mutated metastatic melanoma: a multicentre, open-label, phase 3 randomised controlled trial. Lancet. 2012;380:358–65.
Davalos V, Esteller M. MicroRNAs and cancer epigenetics: a macrorevolution. Curr Opin Oncol. 2010;22:35–45.
Calin GA, Sevignani C, Dumitru CD, Hyslop T, Noch E, Yendamuri S, et al. Human microRNA genes are frequently located at fragile sites and genomic regions involved in cancers. Proc Natl Acad Sci U S A. 2004;101:2999–3004.
Johnson SM, Grosshans H, Shingara J, Byrom M, Jarvis R, Cheng A, et al. RAS is regulated by the let-7 microRNA family. Cell. 2005;120:635–47.
Calin GA, Dumitru CD, Shimizu M, Bichi R, Zupo S, Noch E, et al. Frequent deletions and down-regulation of micro-RNA genes miR15 and miR16 at 13q14 in chronic lymphocytic leukemia. Proc Natl Acad Sci U S A. 2002;99:15524–9.
Tagawa H, Seto M. A microRNA cluster as a target of genomic amplification in malignant lymphoma. Leukemia. 2005;19:2013–6.
Hayashita Y, Osada H, Tatematsu Y, Yamada H, Yanagisawa K, Tomida S, et al. A polycistronic microRNA cluster, miR-17-92, is overexpressed in human lung cancers and enhances cell proliferation. Cancer Res. 2005;65:9628–32.
Philippidou D, Schmitt M, Moser D, Margue C, Nazarov PV, Muller A, et al. Signatures of microRNAs and selected microRNA target genes in human melanoma. Cancer Res. 2010;70:4163–73.
Poliseno L, Haimovic A, Segura MF, Hanniford D, Christos PJ, Darvishian F, et al. Histology-specific microRNA alterations in melanoma. J Invest Dermatol. 2012;132:1860–8.
Xu Y, Brenn T, Brown ER, Doherty V, Melton DW. Differential expression of microRNAs during melanoma progression: miR-200c, miR-205 and miR-211 are downregulated in melanoma and act as tumour suppressors. Br J Cancer. 2012;106:553–61.
Kozubek J, Ma Z, Fleming E, Duggan T, Wu R, Shin DG, et al. In-depth characterization of microRNA transcriptome in melanoma. PLoS One. 2013;8, e72699.
Glud M, Rossing M, Hother C, Holst L, Hastrup N, Nielsen FC, et al. Downregulation of miR-125b in metastatic cutaneous malignant melanoma. Melanoma Res. 2010;20:479–84.
Segura MF, Belitskaya-Levy I, Rose AE, Zakrzewski J, Gaziel A, Hanniford D, et al. Melanoma microRNA signature predicts post-recurrence survival. Clin Cancer Res. 2010;16:1577–86.
Argast GM, Croy CH, Couts KL, Zhang Z, Litman E, Chan DC, et al. Plexin B1 is repressed by oncogenic B-Raf signaling and functions as a tumor suppressor in melanoma cells. Oncogene. 2009;28:2697–709.
Stevens L, McClelland L, Fricke A, Williamson M, Kuo I, Scott G. Plexin B1 suppresses c-Met in melanoma: a role for plexin B1 as a tumor-suppressor protein through regulation of c-Met. J Invest Dermatol. 2010;130:1636–45.
Huang GT, Cunningham KI, Benos PV, Chennubhotla CS. Spectral clustering strategies for heterogeneous disease expression data. Pac Symp Biocomput. 2013;212–223.
Holmfeldt P, Pardieck J, Saulsberry AC, Nandakumar SK, Finkelstein D, Gray JT, et al. Nfix is a novel regulator of murine hematopoietic stem and progenitor cell survival. Blood. 2013;122:2987–96.
Martynoga B, Mateo JL, Zhou B, Andersen J, Achimastou A, Urban N, et al. Epigenomic enhancer annotation reveals a key role for NFIX in neural stem cell quiescence. Genes Dev. 2013;27:1769–86.
Dixon C, Harvey TJ, Smith AG, Gronostajski RM, Bailey TL, Piper M. Nuclear factor one X regulates Bobby sox during development of the mouse forebrain. Cell Mol Neurobiol. 2013;33:867–73.
Pistocchi A, Gaudenzi G, Foglia E, Monteverde S, Moreno-Fortuny A, Pianca A, et al. Conserved and divergent functions of Nfix in skeletal muscle development during vertebrate evolution. Development. 2013;140:1528–36.
Ravichandran V, Sabath BF, Jensen PN, Houff SA, Major EO. Interactions between c-Jun, nuclear factor 1, and JC virus promoter sequences: implications for viral tropism. J Virol. 2006;80:10506–13.
Rademakers R, Eriksen JL, Baker M, Robinson T, Ahmed Z, Lincoln SJ, et al. Common variation in the miR-659 binding-site of GRN is a major risk factor for TDP43-positive frontotemporal dementia. Hum Mol Genet. 2008;17:3631–42.
Budden T, Bowden NA. The role of altered nucleotide excision repair and UVB-induced DNA damage in melanomagenesis. Int J Mol Sci. 2013;14:1132–51.
Steinberg ML, Hubbard K, Utti C, Clas B, Hwang BJ, Hill HZ, et al. Patterns of persistent DNA damage associated with sun exposure and the glutathione S-transferase M1 genotype in melanoma patients. Photochem Photobiol. 2009;85:379–86.
New M, Olzscha H, Liu G, Khan O, Stimson L, McGouran J, et al. A regulatory circuit that involves HR23B and HDAC6 governs the biological response to HDAC inhibitors. Cell Death Differ. 2013;20:1306–16.
Khan O, Fotheringham S, Wood V, Stimson L, Zhang C, Pezzella F, et al. HR23B is a biomarker for tumor sensitivity to HDAC inhibitor-based therapy. Proc Natl Acad Sci U S A. 2010;107:6532–7.
Lovly CM, Dahlman KB, Fohn LE, Su Z, Dias-Santagata D, Hicks DJ, et al. Routine multiplex mutational profiling of melanomas enables enrollment in genotype-driven therapeutic trials. PLoS One. 2012;7, e35309.
Long GV, Menzies AM, Nagrial AM, Haydu LE, Hamilton AL, Mann GJ, et al. Prognostic and clinicopathologic associations of oncogenic BRAF in metastatic melanoma. J Clin Oncol. 2011;29:1239–46.
Curtin JA, Fridlyand J, Kageshita T, Patel HN, Busam KJ, Kutzner H, et al. Distinct sets of genetic alterations in melanoma. N Engl J Med. 2005;353:2135–47.
Gartside MG, Chen H, Ibrahimi OA, Byron SA, Curtis AV, Wellens CL, et al. Loss-of-function fibroblast growth factor receptor-2 mutations in melanoma. Mol Cancer Res. 2009;7:41–54.
Becker D, Lee PL, Rodeck U, Herlyn M. Inhibition of the fibroblast growth factor receptor 1 (FGFR-1) gene in human melanocytes and malignant melanomas leads to inhibition of proliferation and signs indicative of differentiation. Oncogene. 1992;7:2303–13.
Xerri L, Battyani Z, Grob JJ, Parc P, Hassoun J, Bonerandi JJ, et al. Expression of FGF1 and FGFR1 in human melanoma tissues. Melanoma Res. 1996;6:223–30.
Wang Y, Becker D. Antisense targeting of basic fibroblast growth factor and fibroblast growth factor receptor-1 in human melanomas blocks intratumoral angiogenesis and tumor growth. Nat Med. 1997;3:887–93.
Valesky M, Spang AJ, Fisher GW, Farkas DL, Becker D. Noninvasive dynamic fluorescence imaging of human melanomas reveals that targeted inhibition of bFGF or FGFR-1 in melanoma cells blocks tumor growth by apoptosis. Mol Med. 2002;8:103–12.
Mack NA, Whalley HJ, Castillo-Lluva S, Malliri A. The diverse roles of Rac signaling in tumorigenesis. Cell cycle (Georgetown, Tex). 2011;10:1571–81.
Krauthammer M, Kong Y, Ha BH, Evans P, Bacchiocchi A, McCusker JP, et al. Exome sequencing identifies recurrent somatic RAC1 mutations in melanoma. Nat Genet. 2012;44:1006–14.
Hodis E, Watson IR, Kryukov GV, Arold ST, Imielinski M, Theurillat JP, et al. A landscape of driver mutations in melanoma. Cell. 2012;150:251–63.
Ho H, Aruri J, Kapadia R, Mehr H, White MA, Ganesan AK. RhoJ regulates melanoma chemoresistance by suppressing pathways that sense DNA damage. Cancer Res. 2012;72:5516–28.
Bid HK, Roberts RD, Manchanda PK, Houghton PJ. RAC1: an emerging therapeutic option for targeting cancer angiogenesis and metastasis. Mol Cancer Ther. 2013;12:1925–34.
Goncharov K, Zeineldin R. Sensitizing resistant ovarian cancer to chemotherapy through inhibition of small GTPases [abstract]. In: Proceedings of the 103rd Annual Meeting of the American Association for Cancer Research. Chicago: AACR; 2012. p. Abstract nr 4227 2012.
Wilson MA, Zhao F, Letrero R, D'Andrea K, Rimm DL, Kirkwood JM, et al. Correlation of somatic mutations and clinical outcome in melanoma patients treated with carboplatin, paclitaxel, and sorafenib. Clin Cancer Res. 2014;20:3328–37.
Lagani V, Tsamardinos I. Biomarker signature identification in “omics” data with multi-class outcome. Computational Structural Biotechnol J. 2013;6
Cox DR. Regression models and life-tables. J R Stat Soc B. 1972;34:187–220.
Lewis BP, Shih IH, Jones-Rhoades MW, Bartel DP, Burge CB. Prediction of mammalian microRNA targets. Cell. 2003;115:787–98.
This study was conducted by the ECOG-ACRIN Cancer Research Group (Robert L. Comis, MD and Mitchell D. Schnall, MD, PhD, Group Co-Chairs) and supported in part by Public Health Service Grants CA180820, CA180794, CA180844, CA15488, and CA18086; the UPCI Skin Cancer SPORE (P50 CA121973-02) with the UPCI Cancer Biomarkers Facility (P30CA047904); the National Institutes of Health grant (R01LM012087) and the Conquer Cancer Foundation of ASCO Young Investigator Award and the National Cancer Institute, National Institutes of Health, and the Department of Health and Human Services. Its content is solely the responsibility of the authors and does not necessarily represent the official views of the National Cancer Institute, National Institutes of Health, and the Department of Health and Human Services.
Liza C. Villaruz, Grace Huang, Marjorie Romkes, Shama Buch, Tomoko Nukui, Keith Flaherty, Sandra Lee, Melissa A. Wilson, Katherine L. Nathanson, Panayiotis V. Beno, and Hussein A. Tawbi declare that they have no competing interests. John Kirkwood reports personal fees from BMS, Merck, GSK, Celgene, Vical , and Ziopharm.
LCV is responsible for the concept, design, and conduct of the study and manuscript preparation. GH performed the bioinformatics analysis of the mRNA and microRNA expression data in the lab of PVB. MR who performed the gene expression and microRNA profiling. JMK is the Chairman of the Melanoma Committee of ECOG-ACRIN. SCB performed the gene expression and microRNA profiling in MR’s lab. TN performed the gene expression and microRNA profiling in MR’s lab. KTF is the principal investigator of ECOG-ACRIN Cancer Research Group intergroup trial E2603. SJL is the ECOG-ACRIN biostatistician responsible for the E2603 clinical data. MAW performed mutation profiling in the lab of KLN. KLN performed mutation profiling. PVB performed the bioinformatics analysis of the mRNA and microRNA expression data. HAT provided oversight and mentorship in the concept, design, and conduct of the study and manuscript preparation. All authors read and approved the final manuscript.
Liza C. Villaruz and Grace Huang contributed equally to this work.