Review and cross-validation of gene expression signatures and melanoma prognosis

Identification of prognostic mRNAs in metastatic cutaneous melanoma

Cutaneous melanoma is the most common cause of skin cancer-related deaths worldwide. There is an urgent need to identify prognostic biomarkers to facilitate decision-making for treatment of metastatic cutaneous melanoma. Gene expression microarrays and RNA-seq technology have recently improved or changed current prognostic and therapeutic strategies for several cancers. However, according to the current melanoma staging system, prognosis is almost entirely dependent on clinicopathological features. To identify novel prognostic biomarkers, we investigated gene expression and clinical data for patients with cutaneous melanoma from three cohorts of The Cancer Genome Atlas and Gene Expression Omnibus. Kaplan-Meier survival analysis using median values of each gene as cutoff value revealed that nine genes (ABCC3, CAPS2, CCR6, CDCA8, CLU, DPF1, PTK2B, SATB1, and SYNE1) were statistically significant prognostic biomarkers of metastatic cutaneous melanoma in all three independent cohorts. Low expression of two genes (CDCA8 and DPF1) and high expression of seven genes (ABCC3, CAPS2, CCR6, CLU, PTK2B, SATB1, and SYNE) were significantly associated with positive metastatic cutaneous melanoma prognoses. In conclusion, we suggest nine novel prognostic biomarkers for cutaneous metastatic melanoma.

Survival between synchronous and non-synchronous multiple primary cutaneous melanomas-a SEER database analysis

Background

There is no criterion to distinguish synchronous and non-synchronous multiple primary cutaneous melanomas (MPMs). This study aimed to distinguish synchronous and non-synchronous MPMs and compare the survivals of them using the Surveillance, Epidemiology, and End Results database.

Methods

Synchronous and non-synchronous MPMs were distinguished by fitting the double log transformed distribution of the time interval between the first and second primary cutaneous melanomas (TIFtS) through a piecewise linear regression. The overall and melanoma-specific survivals were compared by the Kaplan-Meier method and Cox proportional hazard model through modeling the occurrence of synchronous MPMs as a time-dependent variable.

Results

The distribution of TIFtS was composed by three power-law distributions. According to its first inflection point, synchronous MPMs were defined as tumors that occurred within 2 months. The Kaplain-Meier plot revealed a significant inferior survival for synchronous MPMs than non-synchronous MPMs (P < 0.0001), and the occurrence of synchronous MPM was a risk factor for overall survival of cutaneous melanoma (CM) (hazard ratio: 2.213; (95% CI [2.087-2.346]); P < 0.0001).

Conclusions

This study provided data analysis evidences for using 2 months to distinguish synchronous MPMs and non-synchronous MPMs. Furthermore, the occurrence of synchronous MPM was a risk factor for prognosis of patients with CM.

Local recurrence and metastasis in patients with malignant melanomas after surgery: A single-center analysis of 202 patients in South Korea

Malignant melanoma (MM) is a lethal skin cancer in Western countries. Although the incidence is low in Asians compared to that in Caucasians, it is increasing. However, literature regarding risk factors for prognosis of MM patients who have undergone surgical excision in Asian is limited. This study aimed to investigate the predictive factors for local recurrence and metastasis in MM patients who underwent surgical treatment at a single tertiary-level hospital in Korea. Patients who underwent surgery for MM at our institution between January 1998 and December 2014 were analyzed. We retrospectively investigated risk factors for local recurrence and metastasis after surgery. In cases with distant metastasis, tumor thickness (adjusted Hazard Ratio (HR), 6.139; 95% confidence interval (CI), 2.152 to 17.509; P = 0.001) and increased mitotic number [(0-1/mm² vs 2-6/mm²: adjusted HR, 4.483; 95% CI, 1.233 to 16.303; P = 0.023); (0-1/mm² vs > 6/mm²: adjusted HR, 10.316; 95% CI, 2.871 to 37.063; P < 0.001)] were associated with risk in multivariate analysis. Regarding local recurrence, tumor thickness (T4 [≥4mm] vs T1) was found to be a significant risk factor (adjusted HR, 8.461; 95% CI, 2.514 to 28.474; P = 0.001). Our data revealed tumor thickness and increased mitotic count were significant risk factors for local recurrence and distant metastasis in Korean patients with MM after surgery.

Gene expression profiling in aggressive digital papillary adenocarcinoma sheds light on the architecture of a rare sweat gland carcinoma

BACKGROUND

Sweat gland carcinomas are rare cutaneous adnexal malignancies. Aggressive digital papillary adenocarcinoma (ADPA) represents a very rare subentity, thought to arise almost exclusively from the sweat glands of the fingers and toes. The aetiology of sweat gland carcinomas and ADPA is largely unknown. ADPAs are most likely driven by somatic mutations. However, somatic mutation patterns are largely unexplored, creating barriers to the development of effective therapeutic approaches to the treatment of ADPA.

OBJECTIVES

To investigate the transcriptome profile of ADPA using a sample of eight formalin-fixed, paraffin-embedded tissue samples of ADPA and healthy control tissue.

METHODS

Transcriptome profiling was performed using the Affymetrix PrimeView Human Gene Expression Microarray and findings were validated via reverse transcription of RNA and real-time quantitative polymerase chain reaction.

RESULTS

Transcriptome analyses showed increased tumour expression of 2266 genes, with significant involvement of cell cycle, ribosomal and crucial cancer pathways. Our results point to tumour overexpression of FGFR2 (P = 0·001).

CONCLUSIONS

The results indicate the involvement of crucial oncogenic driver pathways, highlighting cell cycle and ribosomal pathways in the aetiology of ADPA. Suggested tumour overexpression of FGFR2 raises the hope that targeting the fibroblast growth factor (FGF)/FGF receptor axis might be a promising treatment for ADPA and probably for the overall group of sweat gland carcinomas.

Guidelines of care for the management of primary cutaneous melanoma

The incidence of primary cutaneous melanoma continues to increase each year. Melanoma accounts for the majority of skin cancer-related deaths, but treatment is usually curative following early detection of disease. In this American Academy of Dermatology clinical practice guideline, updated treatment recommendations are provided for patients with primary cutaneous melanoma (American Joint Committee on Cancer stages 0-IIC and pathologic stage III by virtue of a positive sentinel lymph node biopsy). Biopsy techniques for a lesion that is clinically suggestive of melanoma are reviewed, as are recommendations for the histopathologic interpretation of cutaneous melanoma. The use of laboratory, molecular, and imaging tests is examined in the initial work-up of patients with newly diagnosed melanoma and for follow-up of asymptomatic patients. With regard to treatment of primary cutaneous melanoma, recommendations for surgical margins and the concepts of staged excision (including Mohs micrographic surgery) and nonsurgical treatments for melanoma in situ, lentigo maligna type (including topical imiquimod and radiation therapy), are updated. The role of sentinel lymph node biopsy as a staging technique for cutaneous melanoma is described, with recommendations for its use in clinical practice. Finally, current data regarding pregnancy and melanoma, genetic testing for familial melanoma, and management of dermatologic toxicities related to novel targeted agents and immunotherapies for patients with advanced disease are summarized.

Altered integrin expression patterns shown by microarray in human cutaneous melanoma

A large variety of molecular pathways in melanoma progression suggests that no individual molecular alteration is crucial in itself. Our aim was to define the molecular alterations underlying metastasis formation. Gene expression profiling was performed using microarray and qRT-PCR to define alterations between matched primary and metastatic melanoma cell lines. These data were integrated with publicly available unmatched tissue data. The invasiveness of cell lines was determined by Matrigel invasion assays and invasive clones from primary melanoma-derived cell lines were also selected. Two metastatic cell line models were created: the regional lymph node WM983A-WM983A-WM983B and the distant lung WM793B-WM793B-1205Lu metastatic models. The majority of metastasis genes were downregulated and enriched in adhesion and ITGA6-B4 pathways. Upregulation of immune pathways was characteristic of distant metastases, whereas increased Rap1 signaling was specific for regional (sub)cutaneous metastases. qRT-PCR analysis of selected integrins (A2, A3, A4, A9, B5, B8, A6, B1, and B3) highlighted the possible importance of ITGA3/4 and B8 in the metastatic process, distinguishing regional and distant metastases. We identified functionally relevant gene clusters that influenced metastasis formation. Our data provide further evidence that integrin expression patterns may be important in distant metastasis formation.

Cross-platform comparison of independent datasets identifies an immune signature associated with improved survival in metastatic melanoma

Platform and study differences in prognostic signatures from metastatic melanoma (MM) gene expression reports often hinder consensus arrival. We performed survival/outcome-based pairwise comparisons of three independent MM gene expression profiles using the threshold-free algorithm rank-rank hypergeometric overlap analysis (RRHO). We found statistically significant overlap for genes overexpressed in favorable outcome (FO) groups, but no overlap for poor outcome (PO) groups. This “favorable outcome signature” (FOS) of 228 genes coinciding on all three overlapping gene lists showed immune function predominated in FO MM. Surprisingly, specific cell signature-enrichment analysis showed B cell-associated genes enriched in FO MM, along with T cell-associated genes. Higher levels of B and T cells (p<0.05) and their relative proximity (p<0.05) were detected in FO-to-PO tumor comparisons from an independent MM patients cohort. Finally, expression of FOS in two independent Stage III MM tumor datasets correctly predicted clinical outcome in 12/14 and 44/70 patients using a weighted gene voting classifier (area under the curve values 0.96 and 0.75, respectively). This RRHO-based, cross-study analysis emphasizes the RRHO approach power, confirms T cells relevance for prolonged MM survival, supports a favorable role for B cells in anti-melanoma immunity, and suggests B cells potential as means of intervention in melanoma treatment.

Recent advances in melanoma research via "omics" platforms

Melanoma has a high mortality rate and metastatic melanoma is highly resistant to conventional therapies. "Omics" fields such as proteomics and microRNA and exosome studies have provided new knowledge to complement the information generated by genomic studies. This work aimed to review the current status of biomarker discovery for melanoma through multi-"omics" platforms. A few sets of novel microRNAs and proteins are described, some of them with important implications in suppressing melanoma at different stages. Upregulation of genes involved in angiogenesis, immunosuppressive factors, modification of stroma, capture of melanoma cells in lymph nodes and factors responsible for tumour cell recruitment have been identified in exosomes, among molecules with other functions. A remarkable series of proteins involved in epithelial-mesenchymal/mesenchymal-epithelial transitions, inflammation, motility, proliferation and progression processes, centrosome amplification, aneuploidy, inhibition of CD8+ effector T-cells, and metastasis in general were identified. Genomic and protein-protein interactions or metabolome levels were not analysed. Proteomics tools such as Orbitrap shotgun mass spectrometry or deep mining proteomic analysis utilizing high-resolution reversed phase nanoseparation in combination with mass spectrometry are also discussed. The application of these tools together with bioinformatics approaches applied to the clinical setting will enable the implementation of personalized medicine in the near future.

A Well-Controlled Experimental System to Study Interactions of Cytotoxic T Lymphocytes with Tumor Cells

While T cell-based immunotherapies are steadily improving, there are still many patients who progress, despite T cell-infiltrated tumors. Emerging evidence suggests that T cells themselves may provoke immune escape of cancer cells. Here, we describe a well-controlled co-culture system for studying the dynamic T cell - cancer cell interplay, using human melanoma as a model. We explain starting material, controls, and culture parameters to establish reproducible and comparable cultures with highly heterogeneous tumor cells. Low passage melanoma cell lines and melanoma-specific CD8+ T cell clones generated from patient blood were cultured together for up to 3 days. Living melanoma cells were isolated from the co-culture system by fluorescence-activated cell sorting. We demonstrate that the characterization of isolated melanoma cells is feasible using flow cytometry for protein expression analysis as well as an Agilent whole human genome microarray and the NanoString technology for differential gene expression analysis. In addition, we identify five genes (ALG12, GUSB, RPLP0, KRBA2, and ADAT2) that are stably expressed in melanoma cells independent of the presence of T cells or the T cell-derived cytokines IFNγ and TNFα. These genes are essential for correct normalization of gene expression data by NanoString. Further to the characterization of melanoma cells after exposure to CTLs, this experimental system might be suitable to answer a series of questions, including how the affinity of CTLs for their target antigen influences the melanoma cell response and whether CTL-induced gene expression changes in melanoma cells are reversible. Taken together, our human T cell - melanoma cell culture system is well suited to characterize immune-related mechanisms in cancer cells.

Clinical and biological determinants of melanoma progression: Should all be considered for clinical management?

Cutaneous melanoma is a heterogeneous disease affecting the regulation of multiple genes and proteins that contribute to melanoma progression. Survival for patients with locally invasive disease varies greatly, even within tumour stages based on current prognostic criteria. This has prompted investigations into the value of additional clinical or biological parameters predicting survival. In particular, the improved knowledge of tumour biology has fed the hope that the outcome may be predicted at the molecular level. The prognostic value of numerous potential biomarkers has therefore been evaluated in protein and gene expression studies, and genomic associations with melanoma prognosis are beginning to emerge. These potential biomarkers interrogate key tumour and host processes important for tumour development and progression, such as proliferation, invasion and migration through epithelial mesenchymal transition or the host immune or vascular responses. This research may allow more individualised information on prognosis if the challenges regarding the quality and validation of studies are overcome.

Biomarkers in melanoma

Malignant melanoma remains the skin cancer with the highest number of mortalities worldwide. While early diagnosis and complete surgical excision remain the best possibility for curing disease, prognosis at the stage of metastasis is still poor. Recent years have brought about considerable advances in terms of understanding the pathogenesis of melanoma and treating advanced disease. The discovery of activating BRAF mutations in around 50% of tumors has led to the introduction of targeted therapies downregulating BRAF signaling output. These have been further refined as combination therapies, which by targeting multiple targets have further improved the clinical outcome. A comparable, potentially even superior therapeutic alternative has been the introduction of immunotherapeutic approaches, including PD-1 and CTLA-4 checkpoint blockade therapies. Despite all genetic knowledge acquired in recent years, a clearly applicable prognostic signature of clinical value has not been established. General prognostic assessment of cutaneous melanoma remains based on clinical and pathological criteria (most importantly tumor thickness). The main challenges lying ahead are to establish a reliable prognostic test effectively determining which tumors will metastasize. Additionally establishing biomarkers which will allow patients to be stratified according to the most promising systemic therapy (immunotherapies and/or BRAF inhibitor therapies) is of utmost importance for patients with metastasized disease. Identifying serum biomarkers enabling disease to be monitored as well as determining tumor properties (i.e. resistance) would also be of great value. While initial results have proven promising, there remains much work to be done.

Identification, Review, and Systematic Cross-Validation of microRNA Prognostic Signatures in Metastatic Melanoma

In metastatic melanoma, it is vital to identify and validate biomarkers of prognosis. Previous studies have systematically evaluated protein biomarkers or mRNA-based expression signatures. No such analyses have been applied to microRNA (miRNA)-based prognostic signatures. As a first step, we identified two prognostic miRNA signatures from publicly available data sets (Gene Expression Omnibus/The Cancer Genome Atlas) of global miRNA expression profiling information. A 12-miRNA signature predicted longer survival after surgery for resection of American Joint Committee on Cancer stage III disease (>4 years, no sign of relapse) and outperformed American Joint Committee on Cancer standard-of-care prognostic markers in leave-one-out cross-validation analysis (error rates 34% and 38%, respectively). A similar 15-miRNA biomarker derived from The Cancer Genome Atlas miRNA-seq data performed slightly worse (39%) than these current biomarkers. Both signatures were then assessed for replication in two independent data sets and subjected to systematic cross-validation together with the three other miRNA-based prognostic signatures proposed in the literature to date. Five miRNAs (miR-142-5p, miR-150-5p, miR-342-3p, miR-155-5p, and miR-146b-5p) were reproducibly associated with patient outcome and have the greatest potential for application in the clinic. Our extensive validation approach highlighted among multiple independent cohorts the translational potential and limitations of miRNA signatures, and pointed to future directions in the analysis of this emerging class of markers.

Molecular markers to complement sentinel node status in predicting survival in patients with high-risk locally invasive melanoma

Sentinel lymph node status is a major prognostic marker in locally invasive cutaneous melanoma. However, this procedure is not always feasible, requires advanced logistics and carries rare but significant morbidity. Previous studies have linked markers of tumour biology to patient survival. In this study, we aimed to combine the predictive value of established biomarkers in addition to clinical parameters as indicators of survival in addition to or instead of sentinel node biopsy in a cohort of high-risk melanoma patients. Patients with locally invasive melanomas undergoing sentinel lymph node biopsy were ascertained and prospectively followed. Information on mortality was validated through the National Death Index. Immunohistochemistry was used to analyse proteins previously reported to be associated with melanoma survival, namely Ki67, p16 and CD163. Evaluation and multivariate analyses according to REMARK criteria were used to generate models to predict disease-free and melanoma-specific survival. A total of 189 patients with available archival material of their primary tumour were analysed. Our study sample was representative of the entire cohort (N = 559). Average Breslow thickness was 2.5 mm. Thirty-two (17%) patients in the study sample died from melanoma during the follow-up period. A prognostic score was developed and was strongly predictive of survival, independent of sentinel node status. The score allowed classification of risk of melanoma death in sentinel node-negative patients. Combining clinicopathological factors and established biomarkers allows prediction of outcome in locally invasive melanoma and might be implemented in addition to or in cases when sentinel node biopsy cannot be performed.

Germline determinants of clinical outcome of cutaneous melanoma

Cutaneous melanoma (CM) is the most lethal form of skin cancer. Despite the constant increase in melanoma incidence, which is in part due to incremental advances in early diagnostic modalities, mortality rates have not improved over the last decade and for advanced stages remain steadily high. While conventional prognostic biomarkers currently in use find significant utility for predicting overall general survival probabilities, they are not sensitive enough for a more personalized clinical assessment on an individual level. In recent years, the advent of genomic technologies has brought the promise of identification of germline DNA alterations that may associate with CM outcomes and hence represent novel biomarkers for clinical utilization. This review attempts to summarize the current state of knowledge of germline genetic factors studied for their impact on melanoma clinical outcomes. We also discuss ongoing problems and hurdles in validating such surrogates, and we also project future directions in discovery of more powerful germline genetic factors with clinical utility in melanoma prognostication.

Prognostic and predictive biomarkers in melanoma: an update

Malignant melanoma is one of the most aggressive cancers. Several new therapeutic strategies that focus on immuno- and/or targeted therapy have been developed, which have entered clinical trials or already been approved. This review provides an update on prognostic and predictive biomarkers in melanoma that may be used to improve the clinical management of patients. Prognostic markers include conventional histopathological characteristics, chromosomal aberrations, gene expression patterns and miRNA profiles. There is a trend towards multi-marker assays and whole-genome molecular screening methods to determine the prognosis of individual patients. Predictive biomarkers, including targeted components of signal transduction, developmental or transcriptional pathways, can be used to determine patient response towards a particular treatment or combination thereof. The rapid evolution of sequencing technologies and multi-marker screening will change the spectrum of patients who become candidates for therapeutic agents, and in addition create new ethical and regulatory challenges.

Lymphatic and blood vasculature in primary cutaneous melanomas of the scalp and neck

BACKGROUND

Scalp/neck melanomas have a poor prognosis, possibly because of a rich vascular supply that prompts tumor cells' dissemination.

METHODS

We compared the accuracy of immunohistochemical (IHC) staining with morphology for the identification of lymphovascular invasion in 156 scalp/neck melanomas. We then analyzed the association of vessel invasion and density with pathological features and survival.

RESULTS

IHC-detected lymphatic vessel invasion (LVI) and blood vessel invasion (BVI) were identified in 34.6% and 13.5% of cases, respectively. IHC increased the LVI/BVI detection compared to morphology (40.4% vs 16.6%; p < .001). The degree of peritumoral and intratumoral blood vessel density (BVD) was greater than lymphatic vessel density (LVD). Ulceration was the only factor independently associated with intratumoral (p = .029) and peritumoral (p = .047) BVD. Tumor thickness was the only independent predictor of survival (p = .002).

CONCLUSION

IHC allows accurate assessment of lymphovascular invasion in scalp/neck melanomas. In these tumors, we observed a high incidence of BVI, which deserves further investigations.

Pembrolizumab joins the anti-PD-1 armamentarium in the treatment of melanoma

Pembrolizumab (MK-3475) is a monoclonal antibody that binds to the PD-1 receptor on T cells and prevents binding to its ligands PD-L1 and PD-L2. Blocking this receptor frees T cells from the inhibitory effects of PD-L1 and allows them to mediate antitumor effects against cancer cells. In a large Phase I study of 411 patients with melanoma, high durable response rates over a range of doses and schedules have been shown with very little toxicity. A Phase III study of pembrolizumab comparing two schedules of administration with the current standard treatment with the anti-CTLA-4 monoclonal antibody is in progress. Combinations with other checkpoint inhibitors as well as other anticancer agents are also being evaluated. Approval of pembrolizumab for the treatment of melanoma is expected.

Lymph-Node Ratio in Patients with Cutaneous Melanoma: A Multi-Institution Prognostic Study

BACKGROUND

Lymph node ratio (LNR)-the number of metastatic lymph nodes (LNs) over the number of excised LNs after lymphadenectomy-is a prognostic factor for many solid tumors, but controversies still exist for skin melanoma. We investigated the prognostic relevance of LNR in melanoma patients and formulated a proposal for considering the LNR in the current American Joint Committee on Cancer (AJCC) N staging system.

METHODS

Retrospective data of 2,526 melanoma patients with LN metastasis from nine Italian institutions were collected in a multicenter database. The prognostic value of the LNR (categorized as A, ≤0.1; B, 0.11-0.25; and C, >0.25) was assessed by multivariable survival analysis.

RESULTS

LNR was a significant independent prognostic factor for melanoma-specific survival (LNR B vs. A: hazard ratio [HR] 1.47, 95 % CI 1.16-1.87, p = 0.002; LNR C vs. A: HR 1.84, 95 % CI 1.29-2.61, p = 0.001). The LNR had prognostic value in patients with AJCC N1a (one positive LN after sentinel LN biopsy [SLNB], HR 2.33, 95 % CI 1.49-3.63, p < 0.001) and N2a (two to three positive LNs after SLNB, HR 1.62, 95 % CI 1.09-2.40, p = 0.016) substages, but not in those with N1b (one clinically positive LN, p = 0.765), N2b (two to three clinically positive LNs, p = 0.165), and N3 (≥ four positive LNs, p = 0.084) substages.

CONCLUSION

The LNR is a prognostic factor in melanoma patients with one (AJCC N1a) and two to three (AJCC N2a) positive LNs after SLNB. This easy-to-obtain parameter should be considered for the staging of melanoma patients with LN metastasis, along with the number of positive LNs.

Principles for the organization of gene-sets

A gene-set, an important concept in microarray expression analysis and systems biology, is a collection of genes and/or their products (i.e. proteins) that have some features in common. There are many different ways to construct gene-sets, but a systematic organization of these ways is lacking. Gene-sets are mainly organized ad hoc in current public-domain databases, with group header names often determined by practical reasons (such as the types of technology in obtaining the gene-sets or a balanced number of gene-sets under a header). Here we aim at providing a gene-set organization principle according to the level at which genes are connected: homology, physical map proximity, chemical interaction, biological, and phenotypic-medical levels. We also distinguish two types of connections between genes: actual connection versus sharing of a label. Actual connections denote direct biological interactions, whereas shared label connection denotes shared membership in a group. Some extensions of the framework are also addressed such as overlapping of gene-sets, modules, and the incorporation of other non-protein-coding entities such as microRNAs.

Independent replication of a melanoma subtype gene signature and evaluation of its prognostic value and biological correlates in a population cohort

Development and validation of robust molecular biomarkers has so far been limited in melanoma research. In this paper we used a large population-based cohort to replicate two published gene signatures for melanoma classification. We assessed the signatures prognostic value and explored their biological significance by correlating them with factors known to be associated with survival (vitamin D) or etiological routes (nevi, sun sensitivity and telomere length). Genomewide microarray gene expressions were profiled in 300 archived tumors (224 primaries, 76 secondaries). The two gene signatures classified up to 96% of our samples and showed strong correlation with melanoma specific survival (P=3 x 10(-4)), Breslow thickness (P=5 x 10(-10)), ulceration (P=9.x10-8) and mitotic rate (P=3 x 10(-7)), adding prognostic value over AJCC stage (adjusted hazard ratio 1.79, 95%CI 1.13-2.83), as previously reported. Furthermore, molecular subtypes were associated with season-adjusted serum vitamin D at diagnosis (P=0.04) and genetically predicted telomere length (P=0.03). Specifically, molecular high-grade tumors were more frequent in patients with lower vitamin D levels whereas high immune tumors came from patients with predicted shorter telomeres. Our data confirm the utility of molecular biomarkers in melanoma prognostic estimation using tiny archived specimens and shed light on biological mechanisms likely to impact on cancer initiation and progression.

One-step extrapolation of the prediction performance of a gene signature derived from a small study

OBJECTIVE

Microarray-related studies often involve a very large number of genes and small sample size. Cross-validating or bootstrapping is therefore imperative to obtain a fair assessment of the prediction/classification performance of a gene signature. A deficiency of these methods is the reduced training sample size because of the partition process in cross-validation and sampling with replacement in bootstrapping. To address this problem, we aim to obtain a prediction performance estimate that strikes a good balance between bias and variance and has a small root mean squared error.

METHODS

We propose to make a one-step extrapolation from the fitted learning curve to estimate the prediction/classification performance of the model trained by all the samples.

RESULTS

Simulation studies show that the method strikes a good balance between bias and variance and has a small root mean squared error. Three microarray data sets are used for demonstration.

CONCLUSIONS

Our method is advocated to estimate the prediction performance of a gene signature derived from a small study.

Vascular dysfunction and increased metastasis of B16F10 melanomas in Shb deficient mice as compared with their wild type counterparts

Background

Shb is a signaling protein downstream of vascular endothelial growth factor receptor-2 and Shb deficiency has been found to restrict tumor angiogenesis. The present study was performed in order to assess metastasis in Shb deficiency using B16F10 melanoma cells.

Methods

B16F10 melanoma cells were inoculated subcutaneously on wild type or Shb +/− mice. Primary tumors were resected and lung metastasis determined after tumor relapse. Lung metastasis was also assessed after bone marrow transplantation of wild type bone marrow to Shb +/− recipients and Shb +/− bone marrow to wild type recipients. Primary tumors were subject to immunofluorescence staining for CD31, VE-cadherin, desmin and CD8, RNA isolation and isolation of vascular fragments for further RNA isolation. RNA was used for real-time RT-PCR and microarray analysis.

Results

Numbers of lung metastases were increased in Shb +/− or −/− mice and this coincided with reduced pericyte coverage and increased vascular permeability. Gene expression profiling of vascular fragments isolated from primary tumors and total tumor RNA revealed decreased expression of different markers for cytotoxic T cells in tumors grown on Shb +/− mice, suggesting that vascular aberrations caused altered immune responses.

Conclusions

It is concluded that a unique combinatorial response of increased vascular permeability and reduced recruitment of cytotoxic CD8+ cells occurs as a consequence of Shb deficiency in B16F10 melanomas. These changes may promote tumor cell intravasation and metastasis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12885-015-1269-y) contains supplementary material, which is available to authorized users.

Molecular pathology of cutaneous melanoma

Cutaneous melanoma is associated with strong prognostic phenotypic features, such as gender, Breslow's thickness and ulceration, although the biological significance of these variables is largely unknown. It is likely that these features are surrogates of important biological events rather than directly promoting cutaneous melanoma progression. In this article, we address the molecular mechanisms that drive these phenotypic changes. Furthermore, we present a comprehensive overview of recurrent genetic abnormalities, both germline and somatic, in relation to cutaneous melanoma subtypes, ultraviolet exposure and anatomical localization, as well as pre-existing and targeted therapy-induced mutations that may contribute to resistance. The increasing knowledge of critically important oncogenes and tumor-suppressor genes is promoting a transition in melanoma diagnosis, in which single-gene testing will be replaced by multiplex and multidimensional analyses that combine classical histopathological characteristics with the molecular profile for the prognostication and selection of melanoma therapy.

How anti-PD1 treatments are changing the management of melanoma

The introduction of immunotherapy based on the blockade of the PD1/PD-L1 checkpoints has been associated with high response rates and durable remissions of disease in patients with metastatic melanoma, to the extent that it is now considered the standard of care for a wide range of patients, irrespective of their BRAF or NRAS mutation status. In addition, more frequent follow-up of patients who are at high risk of recurrence after surgical treatment appears to be justified, as does neoadjuvant treatments in order to render patients treatable by surgery. The limitations of this treatment include failure of some patients to respond, a low rate of complete responses and relapses of the disease during treatment. New initiatives in order to overcome these limitations include the identification of biomarkers for the selection responders and evaluations of treatment combinations that will increase responses and their durability. The latter includes combinations with antibodies against other checkpoints on T cells and cotreatments with inhibitors of resistance pathways in melanoma.

Downregulation of reversion-inducing cysteine-rich protein with Kazal motifs in malignant melanoma: inverse correlation with membrane-type 1-matrix metalloproteinase and tissue inhibitor of metalloproteinase 2

The invasive phenotype of many tumors is associated with an imbalance between the matrix metalloproteinases (MMPs) and their inhibitors, tissue inhibitors of metalloproteinases (TIMPs), and the membrane-anchored reversion-inducing cysteine-rich protein with Kazal motifs (RECK). RECK inhibits MMP-2, MMP-9, and MT1-MMP, and has been linked to patient survival and better prognosis in several types of tumors. However, despite the wide implication of these MMPs in melanoma establishment and progression, the role of RECK in this type of tumor is still unknown. Here, we analyzed the expression of RECK, TIMP1, TIMP2, TIMP3, MT1MMP, MMP2, and MMP9 in two publicly available melanoma microarray datasets and in a panel of human melanoma cell lines. We found that RECK is downregulated in malignant melanoma, accompanied by upregulation of MT1MMP and TIMP2. In both datasets, we observed that the group of samples displaying higher RECK levels show lower median expression levels of MT1MMP and TIMP2 and higher levels of TIMP3. When tested in a sample-wise manner, these correlations were statistically significant. Inverse correlations between RECK, MT1MMP, and TIMP2 were verified in a panel of human melanoma cell lines and in a further reduced model that includes a pair of matched primary tumor-derived and metastasis-derived cell lines. Taken together, our data indicate a consistent correlation between RECK, MT1MMP, and TIMP2 across different models of clinical samples and cell lines and suggest evidence of the potential use of this subset of genes as a gene signature for diagnosing melanoma.

The immune-related role of BRAF in melanoma

Background

The existence of a dichotomy between immunologically active and quiescent tumor phenotypes has been recently recognized in several types of cancer. The activation of a Th1 type of immune signature has been shown to confer better prognosis and likelihood to respond to immunotherapy. However, whether such dichotomy depends on the genetic make‐up of individual cancers is not known yet. BRAF and NRAS mutations are commonly acquired during melanoma progression. Here we explored the role of BRAF and NRAS mutations in influencing the immune phenotype based on a classification previously identified by our group.

Methods

One‐hundred‐thirteen melanoma metastases underwent microarray analysis and BRAF and NRAS genotyping. Allele‐specific PCR was also performed in order to exclude low‐frequency mutations.

Results

Comparison between BRAF and NRAS mutant versus wild type samples identified mostly constituents or regulators of MAPK and related pathways. When testing gene lists discriminative of BRAF, NRAS and MAPK alterations, we found that 112 BRAF‐specific transcripts were able to distinguish the two immune‐related phenotypes already described in melanoma, with the poor phenotype associated mostly with BRAF mutation. Noteworthy, such association was stronger in samples displaying low BRAF mRNA expression. However, when testing NRAS mutations, we were not able to find the same association.

Conclusion

This study suggests that BRAF mutation‐related specific transcripts associate with a poor phenotype in melanoma and provide a nest for further investigation.

BRAF and NRAS status was assessed in 113 melanoma metastases by Sanger sequencing and high sensitive allele‐specific PCR.

The expression of BRAF‐specific genes categorized the metastases in two divergent groups.

The mutant group associated with a poor phenotype.

The association between BRAF mutation and the poor phenotype was stronger in samples displaying low BRAF mRNA expression.

Functional interpretation of BRAF expression‐discriminative genes revealed pathways related to an unfavorable phenotype.

Molecular subtyping of metastatic melanoma based on cell ganglioside metabolism profiles

BACKGROUND

In addition to alterations concerning the expression of oncogenes and onco-suppressors, melanoma is characterized by the presence of distinctive gangliosides (sialic acid carrying glycosphingolipids). Gangliosides strongly control cell surface dynamics and signaling; therefore, it could be assumed that these alterations are linked to modifications of cell behavior acquired by the tumor. On these bases, this work investigated the correlations between melanoma cell ganglioside metabolism profiles and the biological features of the tumor and the survival of patients.

METHODS

Melanoma cell lines were established from surgical specimens of AJCC stage III and IV melanoma patients. Sphingolipid analysis was carried out on melanoma cell lines and melanocytes through cell metabolic labeling employing [3-3H]sphingosine and by FACS. N-glycolyl GM3 was identified employing the 14 F7 antibody. Gene expression was assayed by Real Time PCR. Cell invasiveness was assayed through a Matrigel invasion assay; cell proliferation was determined through the soft agar assay, MTT, and [3H] thymidine incorporation. Statistical analysis was performed using XLSTAT software for melanoma hierarchical clustering based on ganglioside profile, the Kaplan-Meier method, the log-rank (Mantel-Cox) test, and the Mantel-Haenszel test for survival analysis.

RESULTS

Based on the ganglioside profiles, through a hierarchical clustering, we classified melanoma cells isolated from patients into three clusters: 1) cluster 1, characterized by high content of GM3, mainly in the form of N-glycolyl GM3, and GD3; 2) cluster 2, characterized by the appearance of complex gangliosides and by a low content of GM3; 3) cluster 3, which showed an intermediate phenotype between cluster 1 and cluster 3. Moreover, our data demonstrated that: a) a correlation could be traced between patients' survival and clusters based on ganglioside profiles, with cluster 1 showing the worst survival; b) the expression of several enzymes (sialidase NEU3, GM2 and GM1 synthases) involved in ganglioside metabolism was associated with patients' survival; c) melanoma clusters showed different malignant features such as growth in soft agar, invasiveness, expression of anti-apoptotic proteins.

CONCLUSIONS

Ganglioside profile and metabolism is strictly interconnected with melanoma aggressiveness. Therefore, the profiling of melanoma gangliosides and enzymes involved in their metabolism could represent a useful prognostic and diagnostic tool.

High expression of the mismatch repair protein MSH6 is associated with poor patient survival in melanoma

OBJECTIVES

The outcome of patients with primary melanoma (PM) cannot be completely explained based on currently adopted clinical-histopathologic criteria. In this study, we evaluated the potential prognostic value of mismatch repair protein expression in PMs.

METHODS

We examined the immunohistochemical staining of mismatch repair proteins in 18 benign nevi and 101 stage I to III PMs and investigated their association with tumor clinicopathologic variables and melanoma mortality.

RESULTS

Expression of MSH2, MLH1, and PMS2 was high in benign nevi and reduced in a subset of PMs. Conversely, MSH6 expression was absent or extremely low in benign nevi and increased in a subset of PMs. In the multivariate analysis, including sex, age, Breslow thickness, and ulceration, high MSH6 expression in PMs (ie, immunostaining in >20% of tumor cells) was significantly associated with an increased risk of melanoma mortality (relative risk, 3.76; 95% confidence interval, 1.12-12.70).

CONCLUSIONS

MSH6 protein expression can be a valuable marker to improve prognosis assessment in PMs.

Skinomics: past, present and future for diagnostic microarray studies in dermatology

Easily accessible, skin was among the first targets analyzed using 'omics' and dermatology embraced the approaches very early. Microarrays have been used to define disease markers, identify transcriptional changes and even trace the course of treatment. Melanoma and psoriasis have been explored using microarrays. Particularly noteworthy is the multinational mapping of psoriasis susceptibility loci. The transcriptional changes in psoriasis have been identified using hundreds of biopsies. Epidermal keratinocytes have been studied because they respond to UV light, infections, inflammatory and immunomodulating cytokines, toxins and so on. Epidermal differentiation genes are being characterized and are expressed in human epidermal stem cells. Exciting discoveries defining human skin microbiomes have opened a new field of research with great medical potential. Specific to dermatology, the non-invasive skin sampling for microarray studies, using tape stripping, has been developed; it promises to advance dermatology toward 'omics' techniques directly applicable to the personalized medicine of the future.

Immune cells in the melanoma microenvironment hold information for prediction of the risk of recurrence and response to treatment

Melanoma is an immunogenic tumor and immunotherapy treatment has established an increase in disease-free and overall survival in melanoma patients. However, a complex network of immunosuppressive mechanisms has been demonstrated to occur at the tumor site and in locoregional immune districts, such as sentinel lymph nodes (SLNs). The interplay between tumor cells and the local microenvironment leads to a tumor-driven shaping of the immune response that results in a heterogeneous cellular and molecular composition of tumor infiltrating lymphocytes (TILs). Several studies have reported the potential prognostic value of TILs infiltrating primary tumors and the association of 'immune signature' in SLNs and in melanoma metastases with prognosis and responsiveness to immunotherapeutic approaches. However, a systematic and deeper characterization of the local immunological status of TILs and SLNs is still required to refine melanoma stage classification.

Evolving role of tumor antigens for future melanoma therapies

ABSTRACT: 

Human tumor rejection antigens recognized by T lymphocytes were first defined in the early 1990s and the identification of shared tumor-restricted antigens sparked hopes for the development of a therapeutic vaccination to treat cancer, including melanoma. Despite decades of intense preclinical and clinical research, the success of anticancer vaccines based on these antigens has been limited. While melanoma is a highly immunogenic tumor, the ability to prime immunity with vaccines has not generally translated into objective disease regression. However, with the development of small molecules targeting oncogenic proteins, such as V600-mutated BRAF, and immune checkpoint inhibitors with demonstrable long-lasting clinical benefit, new opportunities for antigen-targeted directed therapies are emerging.

Isolation and molecular characterization of circulating melanoma cells

Melanoma is an invasive malignancy with a high frequency of blood-borne metastases, but circulating tumor cells (CTCs) have not been readily isolated. We adapted microfluidic CTC capture to a tamoxifen-driven B-RAF/PTEN mouse melanoma model. CTCs were detected in all tumor-bearing mice and rapidly declined after B-RAF inhibitor treatment. CTCs were shed early from localized tumors, and a short course of B-RAF inhibition following surgical resection was sufficient to dramatically suppress distant metastases. The large number of CTCs in melanoma-bearing mice enabled a comparison of RNA-sequencing profiles with matched primary tumors. A mouse melanoma CTC-derived signature correlated with invasiveness and cellular motility in human melanoma. CTCs were detected in smaller numbers in patients with metastatic melanoma and declined with successful B-RAF-targeted therapy. Together, the capture and molecular characterization of CTCs provide insight into the hematogenous spread of melanoma.

Genetic alterations and personalized medicine in melanoma: progress and future prospects

High-throughput sequencing technologies are providing new insights into the genetic alterations involved in melanomagenesis. It appears likely that most genetic events important in the pathogenesis of melanoma will be discovered over the next few years. Genetic analysis is also increasingly being used to direct patient care. In parallel with the discovery of new genes and the elucidation of molecular pathways important in the development of melanoma, therapies targeting these pathways are becoming available. In other words, the age of personalized medicine has arrived, characterized by molecular profiling of melanoma to identify the relevant genetic alterations and the abnormal signaling mechanisms involved, followed by selection of optimal, individualized therapies. In this review, we summarize the key genetic alterations in melanoma and the development of targeted agents against melanomas bearing specific mutations. These developments in melanoma serve as a model for the implementation of personalized medicine for patients with all cancers.

Disturbed protein-protein interaction networks in metastatic melanoma are associated with worse prognosis and increased functional mutation burden

For disseminated melanoma, new prognostic biomarkers and therapeutic targets are urgently needed. The organization of protein-protein interaction networks was assessed via the transcriptomes of four independent studies of metastatic melanoma and related to clinical outcome and MAP-kinase pathway mutations (BRAF/NRAS). We also examined patient outcome-related differences in a predicted network of microRNAs and their targets. The 32 hub genes with the most reproducible survival-related disturbances in co-expression with their protein partner genes included oncogenes and tumor suppressors, previously known correlates of prognosis, and other proteins not previously associated with melanoma outcome. Notably, this network-based gene set could classify patients according to clinical outcomes with 67-80% accuracy among cohorts. Reproducibly disturbed networks were also more likely to have a higher functional mutation burden than would be expected by chance. The disturbed regions of networks are therefore markers of clinically relevant, selectable tumor evolution in melanoma which may carry driver mutations.

Transcriptome profiling identifies HMGA2 as a biomarker of melanoma progression and prognosis

The genetic alterations contributing to melanoma pathogenesis are incompletely defined, and few independent prognostic features have been identified beyond the clinicopathological characteristics of the primary tumor. We used transcriptome profiling of 46 primary melanomas, 12 melanoma metastases, and 16 normal skin (N) samples to find genes associated with melanoma development and progression. Results were confirmed using immunohistochemistry and real-time PCR and replicated in an independent set of 330 melanomas using AQUA analysis of tissue microarray (TMA). Transcriptome profiling revealed that transcription factor HMGA2, previously unrecognized in melanoma pathogenesis, is significantly upregulated in primary melanoma and metastases (P-values=1.2 × 10(-7) and 9 × 10(-5)) compared with N. HMGA2 overexpression is associated with BRAF/NRAS mutations (P=0.0002). Cox proportional hazard regression model and log-rank test showed that HMGA2 is independently associated with disease-free survival (hazard ratio (HR)=6.3, 95% confidence interval (CI)=1.8-22.3, P=0.004), overall survival (OS) (stratified log-rank P=0.008), and distant metastases-free survival (HR=6.4, 95% CI=1.4-29.7, P=0.018) after adjusting for American Joint Committee on Cancer (AJCC) stage and age at diagnosis. Survival analysis in an independent replication TMA of 330 melanomas confirmed the association of HMGA2 expression with OS (P=0.0211). Our study implicates HMGA2 in melanoma progression and demonstrates that HMGA2 overexpression can serve as an independent predictor of survival in melanoma.

A new 12-gene diagnostic biomarker signature of melanoma revealed by integrated microarray analysis

Genome-wide microarray technology has facilitated the systematic discovery of diagnostic biomarkers of cancers and other pathologies. However, meta-analyses of published arrays often uncover significant inconsistencies that hinder advances in clinical practice. Here we present an integrated microarray analysis framework, based on a genome-wide relative significance (GWRS) and genome-wide global significance (GWGS) model. When applied to five microarray datasets on melanoma published between 2000 and 2011, this method revealed a new signature of 200 genes. When these were linked to so-called ‘melanoma driver’ genes involved in MAPK, Ca²⁺, and WNT signaling pathways we were able to produce a new 12-gene diagnostic biomarker signature for melanoma (i.e., EGFR, FGFR2, FGFR3, IL8, PTPRF, TNC, CXCL13, COL11A1, CHP2, SHC4, PPP2R2C, and WNT4). We have begun to experimentally validate a subset of these genes involved in MAPK signaling at the protein level, including CXCL13, COL11A1, PTPRF and SHC4 and found these to be over-expressed in metastatic and primary melanoma cells in vitro and in situ compared to melanocytes cultured from healthy skin epidermis and normal healthy human skin. While SHC4 has been reported previously to be associated to melanoma, this is the first time CXCL13, COL11A1, and PTPRF have been associated with melanoma on experimental validation. Our computational evaluation indicates that this 12-gene biomarker signature achieves excellent diagnostic power in distinguishing metastatic melanoma from normal skin and benign nevus. Further experimental validation of the role of these 12 genes in a new signaling network may provide new insights into the underlying biological mechanisms driving the progression of melanoma.

4-Bromo-2-(piperidin-1-yl)thiazol-5-yl-phenyl methanone (12b) inhibits Na+/K(+)-ATPase and Ras oncogene activity in cancer cells

The in vitro growth inhibitory activity of 26 thiazoles (including 4-halogeno-2,5-disubtituted-1,3-thiazoles) and 5 thienothiazoles was assessed on a panel of 6 human cancer cell lines, including glioma cell lines. (4-Chloro-2-(piperidin-1-yl)thiazol-5-yl)(phenyl)methanone (12a) and (4-bromo-2-(piperidin-1-yl)thiazol-5-yl)(phenyl)methanone (12b) displayed ~10 times greater in vitro growth inhibitory activity than perillyl alcohol (POH), which therapeutically benefits glioma patients through the inhibition of both alpha-1 Na(+)/K(+)-ATPase (NAK) and Ras oncogene activity. The in vitro cytostatic activities (as revealed by quantitative videomicroscopy) displayed by 12a and 12b were independent of the intrinsic resistance to pro-apoptotic stimuli associated with cancer cells. Compounds 12a and 12b displayed relatively similar inhibitory activities on purified guinea pig brain preparations that mainly express NAK alpha-2 and alpha-3 subunits, whereas only compound 12b was efficacious against purified guinea pig kidney preparations that mainly express the NAK alpha-1 subunit, which is also expressed in gliomas, melanomas and non-small-cell lung cancers NSCLCs.

An attempt at a molecular prediction of metastasis in patients with primary cutaneous melanoma

Background

Current prognostic clinical and morphological parameters are insufficient to accurately predict metastasis in individual melanoma patients. Several studies have described gene expression signatures to predict survival or metastasis of primary melanoma patients, however the reproducibility among these studies is disappointingly low.

Methodology/Principal Findings

We followed extended REMARK/Gould Rothberg criteria to identify gene sets predictive for metastasis in patients with primary cutaneous melanoma. For class comparison, gene expression data from 116 patients with clinical stage I/II (no metastasis) and 72 with III/IV primary melanoma (with metastasis) at time of first diagnosis were used. Significance analysis of microarrays identified the top 50 differentially expressed genes. In an independent data set from a second cohort of 28 primary melanoma patients, these genes were analyzed by multivariate Cox regression analysis and leave-one-out cross validation for association with development of metastatic disease. In a multivariate Cox regression analysis, expression of the genes Ena/vasodilator-stimulated phosphoprotein-like (EVL) and CD24 antigen gave the best predictive value (p = 0.001; p = 0.017, respectively). A multivariate Cox proportional hazards model revealed these genes as a potential independent predictor, which may possibly add (both p = 0.01) to the predictive value of the most important morphological indicator, Breslow depth.

Conclusion/Significance

Combination of molecular with morphological information may potentially enable an improved prediction of metastasis in primary melanoma patients. A strength of the gene expression set is the small number of genes, which should allow easy reevaluation in independent data sets and adequately designed clinical trials.

BRAF mutation, NRAS mutation, and the absence of an immune-related expressed gene profile predict poor outcome in patients with stage III melanoma

Prediction of outcome for melanoma patients with surgically resected macroscopic nodal metastases is very imprecise. We performed a comprehensive clinico-pathologic assessment of fresh-frozen macroscopic nodal metastases and the preceding primary melanoma, somatic mutation profiling, and gene expression profiling to identify determinants of outcome in 79 melanoma patients. In addition to disease stage <II at initial presentation, the following clinical and pathologic factors were independent predictors of improved outcome (odds ratios for survival >4 years, 90% confidence interval): the presence of a nodular component in the primary melanoma (6.8, 0.6-76.0), and small cell size (11.1, 0.8-100.0) or low pigmentation (3.0, 0.8-100.0) in the nodal metastases. Absence of BRAF mutation (20.0, 1.0-1000.0) or NRAS mutation (16.7, 0.6-1000.0) were both favorable prognostic factors. A 46-gene expression signature with strong overrepresentation of immune response genes was predictive of better survival (10.9, 0.4-325.6); in the full cohort, median survival was >100 months in those with the signature, but 10 months in those without. This relationship was validated in two previously published independent stage III melanoma data sets. We conclude that the presence of BRAF mutation, NRAS mutation, and the absence of an immune-related expressed gene profile predict poor outcome in melanoma patients with macroscopic stage III disease.