Report

Frequent Mutation of BAP1 in Metastasizing Uveal Melanomas

See allHide authors and affiliations

Science  03 Dec 2010:
Vol. 330, Issue 6009, pp. 1410-1413
DOI: 10.1126/science.1194472

An Eye on Metastasis

Despite the considerable progress being made in elucidating the cell biology of metastasis, little is known about the genetic alterations that promote metastasis of human tumors, the cause of most cancer deaths. A potentially important clue now emerges from the work of Harbour et al. (p. 1410, published online 4 November), who used an exome-sequencing approach to search for genetic mutations in uveal melanomas, an eye cancer associated with a high rate of fatal metastasis. Remarkably, over 80% of tumor samples with a high metastatic risk had inactivating somatic mutations in the gene encoding BAP1 (BRCA1-associated protein 1), a nuclear protein involved in controlling protein degradation. Thus, in this tumor type, mutational inactivation of BAP1 may be a key event in the acquisition of metastatic competence.

Abstract

Metastasis is a defining feature of malignant tumors and is the most common cause of cancer-related death, yet the genetics of metastasis are poorly understood. We used exome capture coupled with massively parallel sequencing to search for metastasis-related mutations in highly metastatic uveal melanomas of the eye. Inactivating somatic mutations were identified in the gene encoding BRCA1-associated protein 1 (BAP1) on chromosome 3p21.1 in 26 of 31 (84%) metastasizing tumors, including 15 mutations causing premature protein termination and 5 affecting its ubiquitin carboxyl-terminal hydrolase domain. One tumor harbored a frameshift mutation that was germline in origin, thus representing a susceptibility allele. These findings implicate loss of BAP1 in uveal melanoma metastasis and suggest that the BAP1 pathway may be a valuable therapeutic target.

Uveal melanoma (UM) is the most common primary cancer of the eye and has a strong propensity for fatal metastasis (1). UMs are divided into class 1 (low metastatic risk) and class 2 (high metastatic risk) based on a validated multigene clinical prognostic assay included in the TNM classification system (2, 3). However, the genetic basis of metastasis remains unclear. Oncogenic mutations in the Gαq stimulatory subunit GNAQ are common in UM (4), but these mutations occur early in tumorigenesis and are not correlated with molecular class or metastasis (5, 6). On the other hand, class 2 tumors are strongly associated with monosomy 3 (7), suggesting that loss of one copy of chromosome 3 may unmask a mutant gene on the remaining copy that promotes metastasis.

Using exome capture followed by massively parallel sequencing (810), we analyzed two class 2 tumors that were monosomic for chromosome 3 (MM 56 and MM 70) and matching normal DNA from peripheral blood lymphocytes. Both tumors contained inactivating mutations in BAP1, located at chromosome 3p21.1 (Fig. 1A) (10). MM 56 contained a C/G to T/A transition that created a premature termination codon (p.W196X). MM 70 contained a deletion of 11 base pairs in exon 11, leading to a frameshift and premature termination of the BAP1 protein (p.Q322fsX100). The matched normal DNA samples did not contain these mutations, indicating that they were likely to be somatic in origin. No gene on chromosome 3 other than BAP1 contained deleterious somatic mutations that were present in both tumors (table S3).

Fig. 1

Inactivating mutations in BAP1 occur frequently in uveal melanomas. (A) Sanger sequence traces of MM 056 and MM 070 at the sites of the mutations. Location of mutated base in MM 056 and the start of the deletion of MM 070 are indicated (arrows). The noncoding BAP1 strand is shown for MM 070. (B) Map of BAP1 gene and location of BAP1 mutations. BAP1 contains 17 exons (shaded boxes) that encode a 728–amino acid protein. Introns are not to scale. Mutations are shown below the gene figure as indicated. The UCH domain [amino acids (aa) 1 to 188] and UCH37-like domain (ULD) (aa 635 to 693) are indicated (12, 13). The critical Q, C, H, and D residues of the active site (Gln85, Cys91, His169, and Asp184) are indicated with asterisks. The catalytic cysteine is indicated with a circle. Also shown are the NHNY consensus sequence for interaction with HCFC1 (aa 363 to 365, exon 11), nuclear localization signals (NLS) at aa 656 to 661 (exon 15) and aa 717 to 722 (exon 17), the BARD1 binding domain within the region bounded by aa 182 to 240 (13), and the BRCA1 binding domain within aa 598 to 729 (11). (C) Location of BAP1 missense mutations in the UCH domain aligned to the crystal structure of UCH-L3 (21). Three-dimensional structure of UCH-L3 was visualized with MMDB software (22). The small molecule near C91G, H169Q, and S172R represents a suicide inhibitor, illustrating the critical location of these mutations for catalytic activity. (D) BAP1 mRNA levels measured by quantitative RT-PCR in 9 nonmetastasizing class 1 UMs and 28 metastasizing class 2 UMs. (E) Relationship between BAP1 mRNA levels (measured by quantitative RT-PCR) and type of BAP1 mutation in 9 UMs with nonsense and other truncating mutations, 10 UMs with missense mutations (together with small in-frame deletions, splice acceptor, and stop codon read-through mutations), and 4 class 2 UMs in which no BAP1 mutations were detected. Single-letter abbreviations for the amino acid residues are as follows: A, Ala; C, Cys; D, Asp; E, Glu; F, Phe; G, Gly; H, His; I, Ile; K, Lys; L, Leu; M, Met; N, Asn; P, Pro; Q, Gln; R, Arg; S, Ser; T, Thr; V, Val; W, Trp; and Y, Tyr.

BAP1 encodes a nuclear ubiquitin carboxy-terminal hydrolase (UCH), one of several classes of deubiquitinating enzymes (11). In addition to the UCH catalytic domain, BAP1 contains a UCH37-like domain (ULD) (12), binding domains for BRCA1 and BARD1, which form a tumor suppressor heterodimeric complex (13), and a binding domain for HCFC1, which interacts with histone-modifying complexes during cell division (12, 14, 15). BAP1 also interacts with ASXL1 to form the Polycomb group repressive deubiquitinase complex (PR-DUB), which is involved in stem cell pluripotency and other developmental processes (16, 17). BAP1 exhibits tumor suppressor activity in cancer cells (11, 13), and BAP1 mutations have been reported in a small number of breast and lung cancer samples (11, 18).

To further investigate BAP1, genomic DNA from 29 additional class 2 UMs and 26 class 1 UMs were subjected to Sanger resequencing of all BAP1 exons. Altogether, BAP1 mutations were identified in 26 of 31 (84%) class 2 tumors, including 13 out-of-frame deletions and two nonsense mutations leading to premature protein termination, six missense mutations, four in-frame deletions, and one mutation predicted to produce an abnormally extended BAP1 polypeptide (Fig. 1, A to C) (10). Three of the missense mutations affected catalytic residues of the UCH active site (C91 and H169), two occurred elsewhere in the UCH domain, and one affected the ULD (Fig. 1, B and C). All BAP1 missense mutations and in-frame deletions affected phylogenetically conserved amino acids (fig. S1). Only one of 26 class 1 tumors contained a BAP1 mutation (NB101) (10). This case may represent a transition state in which the tumor has sustained a BAP1 mutation but has not yet converted to class 2, suggesting that BAP1 mutations may precede the emergence of the class 2 signature. Somatic BAP1 mutations were also detected in two of three metastatic tumors (10).

One copy of chromosome 3 was missing in all 17 BAP1-mutant class 2 tumors for which cytogenetic data were available, consistent with chromosome 3 loss uncovering recessive BAP1 mutations (10). Normal DNA from 20 patients with BAP1-mutant class 2 primary tumors and the two with metastatic tumors was available and did not contain a BAP1 mutation, indicating that the mutations were somatic in origin. However, we detected one germline mutation (p.E402fsX2; c.1318-1319insA) in the patient with tumor MM 087 (table S2), suggesting that germline alterations in BAP1 can predispose to UM. GNAQ mutation status was available in 15 cases. GNAQ mutations were present in 4 of 9 BAP1 mutant tumors and 3 of 6 BAP1 wild-type tumors, indicating that there was no correlation between GNAQ and BAP1 mutation status.

Quantitative reverse transcription polymerase chain reaction (RT-PCR) showed that BAP1 mRNA levels were significantly lower in class 2 tumors compared with class 1 tumors (P < 0.0001) (Fig. 1D). Truncating mutations were associated with significantly lower mRNA levels than missense mutations (P = 0.001) (Fig. 1E), consistent with nonsense-mediated mRNA decay in the former group. Class 2 tumors in which BAP1 mutations were not identified expressed very low levels of BAP1 mRNA (Fig. 1E).

To determine whether the low BAP1 mRNA levels in class 2 tumors without detectable BAP1 mutations may be explained by DNA methylation, we performed a preliminary analysis of DNA methylation of BAP1. This did not reveal a convincing difference between class 1 and class 2 tumors (10). However, analysis of the BAP1 promoter was limited by an unusually complex CpG island that will require further work to resolve. Thus, we cannot rule out a role for methylation in class 2 tumors in which BAP1 mutations were not found. However, with almost 85% of class 2 tumors harboring mutations, we do not expect that methylation will be a major mechanism of BAP1 inactivation. An alternative explanation is that these tumors may contain very large deletions of the BAP1 locus not detectable by our sequencing method.

Immunofluorescence revealed abundant nuclear BAP1 protein in two class 1 tumors but virtually none in four BAP1 mutant class 2 tumors (fig. S3). This was expected for the two tumors with mutations expected to cause premature protein terminations (MM 091 and MM 100), but it was surprising for the two tumors with missense mutations (MM 071 and MM 135) and suggests that these mutations lead to protein instability.

RNA interference (RNAi)–mediated knock down of BAP1 in 92.1 UM cells, which did not harbor a detectable BAP1 mutation, recapitulated many characteristics of the dedifferentiated class 2 UM phenotype (19). Cells transfected with control small interfering RNA (siRNA) exhibited typical melanocytic morphology, including dendritic projections and cytoplasmic melanosomes (Fig. 2A), whereas cells transfected with BAP1 siRNA lost these features, developed a rounded epithelioid morphology, and grew as multicellular nonadherent spheroids, strikingly similar to the features of class 2 clinical biopsy samples (Fig. 2A). Microarray gene expression profiling of 92.1 UM cells transfected with control versus BAP1 siRNA showed that most of the top genes that discriminate between class 1 and class 2 tumors shifted in the class 2 direction in BAP1-depleted cells compared with control cells (fig. S4). Similarly, depletion of BAP1 shifted the gene expression profile of the multi-gene clinical prognostic assay toward the class 2 signature (Fig. 2B). BAP1 depletion caused a reduction in mRNA levels of neural crest migration genes (ROBO1), melanocyte differentiation genes (CTNNB1, EDNRB, and SOX10), and other genes that are down-regulated in class 2 tumors (LMCD1 and LTA4H) (19). In contrast, BAP1 depletion caused an increase in mRNA levels of CDH1 and the proto-oncogene KIT, which are highly expressed in class 2 tumors (20). Similar results were seen in other UM cell lines and with an independent BAP1 siRNA (10).

Fig. 2

UM cells depleted of BAP1 acquire properties that are typical of metastasizing class 2 tumor cells. (A) Phase-contrast photomicrographs of 92.1 uveal melanoma cells transfected with BAP1 or control siRNA at the indicated days. Bottom panels show representative examples of class 1 and class 2 uveal melanoma cells obtained from patient biopsy samples (Papanicolaou stain). Scale bars, 10 μm. (B) 92.1 cells transfected with BAP1 siRNA and evaluated after 5 days. BAP1 protein levels were efficiently depleted to less than 95% of control levels (see Western blot). Upper panel depicts principal component analysis to show the effect of BAP1 knockdown on gene expression signature. The small spheres represent the training set of known class 1 (blue) and class 2 (red) tumors. Large spheres represent the control-transfected (gray) and BAP1 siRNA–transfected (red) cells. The lower panel depicts mRNA levels measured by quantitative RT-PCR of a panel of melanocyte lineage genes, presented as relative change (fold change) in BAP1 siRNA/control siRNA–transfected cells. Results are representative of three independent experiments.

GNAQ mutations occur early in UM and are not sufficient for malignant transformation (4), but they may create a dependency of the tumor cells on constitutive GNAQ activity. In contrast, BAP1 mutations occur later in UM progression and coincide with the onset of metastatic behavior. Thus, simultaneous targeting of both genetic alterations might have synergistic therapeutic effects. One potential strategy to counteract the effects of BAP1 mutation would be to inhibit the RING1 ubiquitinating activity that normally opposes the deubiquitinating activity of BAP1 (16). Our findings strongly implicate mutational inactivation of BAP1 as a key event in the acquisition of metastatic competence in UM, and they expand the role of BAP1 and other deubiquitinating enzymes as potential therapeutic targets in cancer.

Supporting Online Material

www.sciencemag.org/cgi/content/full/science.1194472/DC1

Materials and Methods

Figs. S1 to S5

Tables S1 to S3

References

References and Notes

  1. Materials and methods are available as supporting material on Science Online.
  2. This work was supported by grants to J.W.H. from the National Cancer Institute (R01 CA125970), Barnes-Jewish Hospital Foundation, Kling Family Foundation, Tumori Foundation, Horncrest Foundation, and a Research to Prevent Blindness David F. Weeks Professorship, and by awards to the Department of Ophthalmology and Visual Sciences at Washington University from a Research to Prevent Blindness, Inc., unrestricted grant, and the NIH Vision Core Grant P30 EY02687c. E.D.O.R. was supported by NIH National Institute of Arthritis and Musculoskeletal and Skin Diseases training grant AR007279-31A1. We thank J. Hoisington-Lopez from the Center for Genome Sciences (Washington University School of Medicine) for running Solexa paired end sequences and M. Lovett for comments on the manuscript. J.W.H. and Washington University may receive income based on a license of related technology by the university to Castle Biosciences, Inc.
View Abstract

Navigate This Article