Esophageal adenocarcinoma is a cancer with rising incidence and poor survival. Most such cancers arise in a specialized intestinal metaplastic epithelium, which is diagnostic of Barrett's esophagus. In a genome-wide association study, we compared esophageal adenocarcinoma cases (n = 2,390) and individuals with precancerous Barrett's esophagus (n = 3,175) with 10,120 controls in 2 phases. For the combined case group, we identified three new associations. The first is at 19p13 (rs10419226: P = 3.6 × 10(-10)) in CRTC1 (encoding CREB-regulated transcription coactivator), whose aberrant activation has been associated with oncogenic activity. A second is at 9q22 (rs11789015: P = 1.0 × 10(-9)) in BARX1, which encodes a transcription factor important in esophageal specification. A third is at 3p14 (rs2687201: P = 5.5 × 10(-9)) near the transcription factor FOXP1, which regulates esophageal development. We also refine a previously reported association with Barrett's esophagus near the putative tumor suppressor gene FOXF1 at 16q24 and extend our findings to now include esophageal adenocarcinoma.

Deleterious germline variants in CDKN2A account for around 40% of familial melanoma cases, and rare variants in CDK4, BRCA2, BAP1 and the promoter of TERT have also been linked to the disease. Here we set out to identify new high-penetrance susceptibility genes by sequencing 184 melanoma cases from 105 pedigrees recruited in the UK, The Netherlands and Australia that were negative for variants in known predisposition genes. We identified families where melanoma cosegregates with loss-of-function variants in the protection of telomeres 1 gene (POT1), with a proportion of family members presenting with an early age of onset and multiple primary tumors. We show that these variants either affect POT1 mRNA splicing or alter key residues in the highly conserved oligonucleotide/oligosaccharide-binding (OB) domains of POT1, disrupting protein-telomere binding and leading to increased telomere length. These findings suggest that POT1 variants predispose to melanoma formation via a direct effect on telomeres.

At least 17 genomic regions are established as harboring melanoma susceptibility variants, in most instances with genome-wide levels of significance and replication in independent samples. Based on genome-wide single nucleotide polymorphism (SNP) data augmented by imputation to the 1,000 Genomes reference panel, we have fine mapped these regions in over 5,000 individuals with melanoma (mainly from the GenoMEL consortium) and over 7,000 ethnically matched controls. A penalized regression approach was used to discover those SNP markers that most parsimoniously explain the observed association in each genomic region. For the majority of the regions, the signal is best explained by a single SNP, which sometimes, as in the tyrosinase region, is a known functional variant. However in five regions the explanation is more complex. At the CDKN2A locus, for example, there is strong evidence that not only multiple SNPs but also multiple genes are involved. Our results illustrate the variability in the biology underlying genome-wide susceptibility loci and make steps toward accounting for some of the "missing heritability."

Truncating germline mutations in the tumor suppressor gene BRCA-1 associated protein-1 (BAP1) have been reported in families predisposed to developing a wide range of different cancer types including uveal melanoma and cutaneous melanoma. There has also been an association between amelanotic tumor development and germline BAP1 mutation suggesting a possible phenotypic characteristic of BAP1 mutation carriers. Though there have been many types of cancer associated with germline BAP1 mutation, the full spectrum of disease association is yet to be ascertained. Here we describe a Danish family with predominantly uveal melanoma but also a range of other tumor types including lung, neuroendocrine, stomach, and breast cancer; as well as pigmented skin lesions. Whole-exome sequencing identified a BAP1 splice mutation located at c.581-2A>G, which leads to a premature truncation of BAP1 in an individual with uveal melanoma. This mutation was carried by several other family members with melanoma or various cancers. The finding expands on the growing profile of BAP1 as an important uveal and cutaneous melanoma tumor suppressor gene and implicates its involvement in the development of lung, and stomach cancer.

We conducted a multi-stage genome-wide association study of natural hair color in more than 10,000 men and women of European ancestry from the United States and Australia. An initial analysis of 528,173 single nucleotide polymorphisms (SNPs) genotyped on 2,287 women identified IRF4 and SLC24A4 as loci highly associated with hair color, along with three other regions encompassing known pigmentation genes. We confirmed these associations in 7,028 individuals from three additional studies. Across these four studies, SLC24A4 rs12896399 and IRF4 rs12203592 showed strong associations with hair color, with p = 6.0x10(-62) and p = 7.46x10(-127), respectively. The IRF4 SNP was also associated with skin color (p = 6.2x10(-14)), eye color (p = 6.1x10(-13)), and skin tanning response to sunlight (p = 3.9x10(-89)). A multivariable analysis pooling data from the initial GWAS and an additional 1,440 individuals suggested that the association between rs12203592 and hair color was independent of rs1540771, a SNP between the IRF4 and EXOC2 genes previously found to be associated with hair color. After adjustment for rs12203592, the association between rs1540771 and hair color was not significant (p = 0.52). One variant in the MATP gene was associated with hair color. A variant in the HERC2 gene upstream of the OCA2 gene showed the strongest and independent association with hair color compared with other SNPs in this region, including three previously reported SNPs. The signals detected in a region around the MC1R gene were explained by MC1R red hair color alleles. Our results suggest that the IRF4 and SLC24A4 loci are associated with human hair color and skin pigmentation.

We concurrently examine the whole genome, transcriptome, methylome, and immune cell infiltrates in baseline tumors from 77 patients with advanced cutaneous melanoma treated with anti-PD-1 with or without anti-CTLA-4. We show that high tumor mutation burden (TMB), neoantigen load, expression of IFNγ-related genes, programmed death ligand expression, low PSMB8 methylation (therefore high expression), and T cells in the tumor microenvironment are associated with response to immunotherapy. No specific mutation correlates with therapy response. A multivariable model combining the TMB and IFNγ-related gene expression robustly predicts response (89% sensitivity, 53% specificity, area under the curve [AUC], 0.84); tumors with high TMB and a high IFNγ signature show the best response to immunotherapy. This model validates in an independent cohort (80% sensitivity, 59% specificity, AUC, 0.79). Except for a JAK3 loss-of-function mutation, for patients who did not respond as predicted there is no obvious biological mechanism that clearly explained their outlier status, consistent with intratumor and intertumor heterogeneity in response to immunotherapy.

Acral melanoma (AM) tumors arise on the palms, soles, fingers, toes, and nailbeds. A comprehensive systematic meta-analysis of AM genomic aberrations has not been conducted to date. A literature review was carried out to identify studies sequencing AM. Whole-genome/exome data from 181 samples were identified. Targeted panel sequencing data from MSK-IMPACT were included as a validation cohort (n = 92), and studies using targeted hot spot sequencing were also collated for BRAF (n = 26 studies), NRAS (n = 21), and KIT (n = 32). Statistical analysis indicated BRAF, NRAS, PTEN, TYRP1, and KIT as significantly mutated genes. Frequent copy-number aberrations were also found for important cancer genes, such as CDKN2A, KIT, MDM2, CCND1, CDK4, and PAK1, among others. Mapping genomic alterations within the context of the hallmarks of cancer identified four components frequently altered, including (i) sustained proliferative signaling and (ii) evading growth suppression, (iii) genome instability and mutation, and (iv) enabling replicative immortality. This analysis provides the largest analysis of genomic aberrations in AM in the literature to date and highlights pathways that may be therapeutically targetable.

We previously identified miR-4731-5p (miR-4731) as a melanoma-enriched microRNA following comparison of melanoma with other cell lines from solid malignancies. Additionally, miR-4731 has been found in serum from melanoma patients and expressed less abundantly in metastatic melanoma tissues from stage IV patients relative to stage III patients. As miR-4731 has no known function, we used biotin-labelled miRNA duplex pull-down to identify binding targets of miR-4731 in three melanoma cell lines (HT144, MM96L and MM253). Using the miRanda miRNA binding algorithm, all pulled-down transcripts common to the three cell lines (n=1092) had potential to be targets of miR-4731 and gene-set enrichment analysis of these (via STRING v9.1) highlighted significantly associated genes related to the 'cell cycle' pathway and the 'melanosome'. Following miR-4731 overexpression, a selection (n=81) of pull-down transcripts underwent validation using a custom qRT-PCR array. These data revealed that miR-4731 regulates multiple genes associated with the cell cycle (e.g. CCNA2, ORC5L, and PCNA) and the melanosome (e.g. RAB7A, CTSD, and GNA13). Furthermore, members of the synovial sarcoma X breakpoint family (SSX) (melanoma growth promoters) were also down-regulated (e.g. SSX2, SSX4, and SSX4B) as a result of miR-4731 overexpression. Moreover, this down-regulation of mRNA expression resulted in ablation or reduction of SSX4 protein, which, in keeping with previous studies, resulted in loss of 2D colony formation. We therefore speculate that loss of miR-4731 expression in stage IV patient tumours supports melanoma growth by, in part; reducing its regulatory control of SSX expression levels.

Next generation sequencing of uveal melanoma (UM) samples has identified a number of recurrent oncogenic or loss-of-function mutations in key driver genes including: GNAQ, GNA11, EIF1AX, SF3B1 and BAP1. To search for additional driver mutations in this tumor type we carried out whole-genome or whole-exome sequencing of 28 tumors or primary cell lines. These samples have a low mutation burden, with a mean of 10.6 protein changing mutations per sample (range 0 to 53). As expected for these sun-shielded melanomas the mutation spectrum was not consistent with an ultraviolet radiation signature, instead, a BRCA mutation signature predominated. In addition to mutations in the known UM driver genes, we found a recurrent mutation in PLCB4 (c.G1888T, p.D630Y, NM_000933), which was validated using Sanger sequencing. The identical mutation was also found in published UM sequence data (1 of 56 tumors), supporting its role as a novel driver mutation in UM. PLCB4 p.D630Y mutations are mutually exclusive with mutations in GNA11 and GNAQ, consistent with PLCB4 being the canonical downstream target of the former gene products. Taken together these data suggest that the PLCB4 hotspot mutation is similarly a gain-of-function mutation leading to activation of the same signaling pathway, promoting UM tumorigenesis.

The epigenetic modifier EZH2 is part of the polycomb repressive complex that suppresses gene expression via histone methylation. Activating mutations in EZH2 are found in a subset of melanoma that contributes to disease progression by inactivating tumor suppressor genes. In this study we have targeted EZH2 with a specific inhibitor (GSK126) or depleted EZH2 protein by stable shRNA knockdown. We show that inhibition of EZH2 has potent effects on the growth of both wild-type and EZH2 mutant human melanoma in vitro particularly in cell lines harboring the EZH2Y646 activating mutation. This was associated with cell cycle arrest, reduced proliferative capacity in both 2D and 3D culture systems, and induction of apoptosis. The latter was caspase independent and mediated by the release of apoptosis inducing factor (AIFM1) from mitochondria. Gene expression arrays showed that several well characterized tumor suppressor genes were reactivated by EZH2 inhibition. This included activating transcription factor 3 (ATF3) that was validated as an EZH2 target gene by ChIP-qPCR. These results emphasize a critical role for EZH2 in the proliferation and viability of melanoma and highlight the potential for targeted therapy against EZH2 in treatment of patients with melanoma.

We performed a genome-wide association study of melanoma in a discovery cohort of 2,168 Australian individuals with melanoma and 4,387 control individuals. In this discovery phase, we confirm several previously characterized melanoma-associated loci at MC1R, ASIP and MTAP-CDKN2A. We selected variants at nine loci for replication in three independent case-control studies (Europe: 2,804 subjects with melanoma, 7,618 control subjects; United States 1: 1,804 subjects with melanoma, 1,026 control subjects; United States 2: 585 subjects with melanoma, 6,500 control subjects). The combined meta-analysis of all case-control studies identified a new susceptibility locus at 1q21.3 (rs7412746, P = 9.0 × 10(-11), OR in combined replication cohorts of 0.89 (95% CI 0.85-0.95)). We also show evidence suggesting that melanoma associates with 1q42.12 (rs3219090, P = 9.3 × 10(-8)). The associated variants at the 1q21.3 locus span a region with ten genes, and plausible candidate genes for melanoma susceptibility include ARNT and SETDB1. Variants at the 1q21.3 locus do not seem to be associated with human pigmentation or measures of nevus density.

RAF inhibitor therapy yields significant reductions in tumour burden in the majority of V600E-positive melanoma patients; however, resistance occurs within 2-18 months. Here we demonstrate that the mixed lineage kinases (MLK1-4) are MEK kinases that reactivate the MEK/ERK pathway in the presence of RAF inhibitors. Expression of MLK1-4 mediates resistance to RAF inhibitors and promotes survival in V600E-positive melanoma cell lines. Furthermore, we observe upregulation of the MLKs in 9 of 21 melanoma patients with acquired drug resistance. Consistent with this observation, MLKs promote resistance to RAF inhibitors in mouse models and contribute to acquired resistance in a cell line model. Lastly, we observe that a majority of MLK1 mutations identified in patients are gain-of-function mutations. In summary, our data demonstrate a role for MLKs as direct activators of the MEK/ERK pathway with implications for melanomagenesis and resistance to RAF inhibitors.

Analysis of 501 melanoma exomes identified RASA2, encoding a RasGAP, as a tumor-suppressor gene mutated in 5% of melanomas. Recurrent loss-of-function mutations in RASA2 were found to increase RAS activation, melanoma cell growth and migration. RASA2 expression was lost in ≥30% of human melanomas and was associated with reduced patient survival. These findings identify RASA2 inactivation as a melanoma driver and highlight the importance of RasGAPs in cancer.

MicroRNAs (miRNAs) are a family of small, non-coding RNA species functioning as negative regulators of multiple target genes including tumour suppressor genes and oncogenes. Many miRNA gene loci are located within cancer-associated genomic regions. To identify potential new amplified oncogenic and/or deleted tumour suppressing miRNAs in lung cancer, we inferred miRNA gene dosage from high dimensional arrayCGH data. From miRBase v9.0 (http://microrna.sanger.ac.uk), 474 human miRNA genes were physically mapped to regions of chromosomal loss or gain identified from a high-resolution genome-wide arrayCGH study of 132 primary non-small cell lung cancers (NSCLCs) (a training set of 60 squamous cell carcinomas and 72 adenocarcinomas). MiRNAs were selected as candidates if their immediately flanking probes or host gene were deleted or amplified in at least 25% of primary tumours using both Analysis of Copy Errors algorithm and fold change (≥ ± 1.2) analyses. Using these criteria, 97 miRNAs mapped to regions of aberrant copy number. Analysis of three independent published lung cancer arrayCGH datasets confirmed that 22 of these miRNA loci showed directionally concordant copy number variation. MiR-218, encoded on 4p15.31 and 5q35.1 within two host genes (SLIT2 and SLIT3), in a region of copy number loss, was selected as a priority candidate for follow-up as it is reported as underexpressed in lung cancer. We confirmed decreased expression of mature miR-218 and its host genes by qRT-PCR in 39 NSCLCs relative to normal lung tissue. This downregulation of miR-218 was found to be associated with a history of cigarette smoking, but not human papilloma virus. Thus, we show for the first time that putative lung cancer-associated miRNAs can be identified from genome-wide arrayCGH datasets using a bioinformatics mapping approach, and report that miR-218 is a strong candidate tumour suppressing miRNA potentially involved in lung cancer.

Barrett's esophagus (BE) is the metaplastic replacement of squamous with columnar epithelium in the esophagus, as a result of reflux. It is the major risk factor for the development of esophageal adenocarcinoma (EAC). Methylation of CpG dinucleotides of normally unmethylated genes is associated with silencing of their expression, and is common in EAC. This study was designed to determine at what stage, in the progression from BE to EAC, methylation of key genes occurs.

Knowledge of key drivers and therapeutic targets in mucosal melanoma is limited due to the paucity of comprehensive mutation data on this rare tumor type. To better understand the genomic landscape of mucosal melanoma, here we describe whole genome sequencing analysis of 67 tumors and validation of driver gene mutations by exome sequencing of 45 tumors. Tumors have a low point mutation burden and high numbers of structural variants, including recurrent structural rearrangements targeting TERT, CDK4 and MDM2. Significantly mutated genes are NRAS, BRAF, NF1, KIT, SF3B1, TP53, SPRED1, ATRX, HLA-A and CHD8. SF3B1 mutations occur more commonly in female genital and anorectal melanomas and CTNNB1 mutations implicate a role for WNT signaling defects in the genesis of some mucosal melanomas. TERT aberrations and ATRX mutations are associated with alterations in telomere length. Mutation profiles of the majority of mucosal melanomas suggest potential susceptibility to CDK4/6 and/or MEK inhibitors.

Background: Little is known about the prognostic significance of somatically mutated genes in metastatic melanoma (MM). We have employed a combined clinical and bioinformatics approach on tumor samples from cutaneous melanoma (SKCM) as part of The Cancer Genome Atlas project (TCGA) to identify mutated genes with potential clinical relevance. Methods: After limiting our DNA sequencing analysis to MM samples (n = 356) and to the CANCER CENSUS gene list, we filtered out mutations with low functional significance (snpEFF). We performed Cox analysis on 53 genes that were mutated in ≥3% of samples, and had ≥50% difference in incidence of mutations in deceased subjects versus alive subjects. Results: Four genes were potentially prognostic [RAC1, FGFR1, CARD11, CIITA; false discovery rate (FDR) < 0.2]. We identified 18 additional genes (e.g., SPEN, PDGFRB, GNAS, MAP2K1, EGFR, TSC2) that were less likely to have prognostic value (FDR < 0.4). Most somatic mutations in these 22 genes were infrequent (< 10%), associated with high somatic mutation burden, and were evenly distributed across all exons, except for RAC1 and MAP2K1. Mutations in only 9 of these 22 genes were also identified by RNA sequencing in >75% of the samples that exhibited corresponding DNA mutations. The low frequency, UV signature type and RNA expression of the 22 genes in MM samples were confirmed in a separate multi-institution validation cohort (n = 413). An underpowered analysis within a subset of this validation cohort with available patient follow-up (n = 224) showed that somatic mutations in SPEN and RAC1 reached borderline prognostic significance [log-rank favorable (p = 0.09) and adverse (p = 0.07), respectively]. Somatic mutations in SPEN, and to a lesser extent RAC1, were not associated with definite gene copy number or RNA expression alterations. High (>2+) nuclear plus cytoplasmic expression intensity for SPEN was associated with longer melanoma-specific overall survival (OS) compared to lower (≤ 2+) nuclear intensity (p = 0.048). We conclude that expressed somatic mutations in infrequently mutated genes beyond the well-characterized ones (e.g., BRAF, RAS, CDKN2A, PTEN, TP53), such as RAC1 and SPEN, may have prognostic significance in MM.

In general, melanoma can be considered as a UV-driven disease with an aggressive metastatic course and high mutational load, with only few tumors (acral, mucosal, and uveal melanomas) not induced by sunlight and possessing a lower mutational load. The most commonly activated pathway in melanoma is the mitogen-activated protein kinase (MAPK) pathway. However, the prognostic significance of mutational stratification is unclear and needs further investigation. Here, in silico we combined mutation data from 162 melanomas subjected to targeted deep sequencing with mutation data from three published studies. Tumors from 870 patients were grouped according to BRAF, RAS, NF1 mutation or triple-wild-type status and correlated with tumor and patient characteristics. We found that the NF1-mutated subtype had a higher mutational burden and strongest UV mutation signature. Searching for co-occurring mutated genes revealed the RASopathy genes PTPN11 and RASA2, as well as another RAS domain-containing gene RASSF2 enriched in the NF1 subtype after adjustment for mutational burden. We found that a larger proportion of the NF1-mutant tumors were from males and with older age at diagnosis. Importantly, we found an increased risk of death from melanoma (disease-specific survival, DSS; HR, 1.9; 95% CI, 1.21-3.10; P = 0.046) and poor overall survival (OS; HR, 2.0; 95% CI, 1.28-2.98; P = 0.01) in the NF1 subtype, which remained significant after adjustment for age, gender, and lesion type (DSS P = 0.03, OS P = 0.06, respectively). Melanoma genomic subtypes display different biological and clinical characteristics. The poor outcome observed in the NF1 subtype highlights the need for improved characterization of this group.

Circulating tumor DNA (ctDNA) may serve as a surrogate to tissue biopsy for noninvasive identification of mutations across multiple genetic loci and for disease monitoring in melanoma. In this study, we compared the mutation profiles of tumor biopsies and plasma ctDNA from metastatic melanoma patients using custom sequencing panels targeting 30 melanoma-associated genes. Somatic mutations were identified in 20 of 24 melanoma biopsies, and 16 of 20 (70%) matched-patient plasmas had detectable ctDNA. In a subgroup of seven patients for whom matching tumor tissue and plasma were sequenced, 80% of the mutations found in tumor tissue were also detected in ctDNA. However, TERT promoter mutations were only detected by ddPCR, and promoter mutations were consistently found at lower concentrations than other driver mutations in longitudinal samples. In vitro experiments revealed that mutations in promoter regions of TERT and DPH3 are underrepresented in ctDNA. While the results underscore the utility of using ctDNA as an alternative to tissue biopsy for genetic profiling and surveillance of the disease, our study highlights the underrepresentation of promoter mutations in ctDNA and its potential impact on quantitative liquid biopsy applications.

Uveal melanoma (UM) is the most common intraocular tumour in adults and despite surgical or radiation treatment of primary tumours, ~50% of patients progress to metastatic disease. Therapeutic options for metastatic UM are limited, with clinical trials having little impact. Here we perform whole-genome sequencing (WGS) of 103 UM from all sites of the uveal tract (choroid, ciliary body, iris). While most UM have low tumour mutation burden (TMB), two subsets with high TMB are seen; one driven by germline MBD4 mutation, and another by ultraviolet radiation (UVR) exposure, which is restricted to iris UM. All but one tumour have a known UM driver gene mutation (GNAQ, GNA11, BAP1, PLCB4, CYSLTR2, SF3B1, EIF1AX). We identify three other significantly mutated genes (TP53, RPL5 and CENPE).