Uveal melanoma (UM) is a highly metastatic cancer that, in contrast to cutaneous melanoma, is largely unresponsive to checkpoint immunotherapy. Here, we interrogate the tumor microenvironment at single-cell resolution using scRNA-seq of 59,915 tumor and non-neoplastic cells from 8 primary and 3 metastatic samples. Tumor cells reveal novel subclonal genomic complexity and transcriptional states. Tumor-infiltrating immune cells comprise a previously unrecognized diversity of cell types, including CD8⁺ T cells predominantly expressing the checkpoint marker LAG3, rather than PD1 or CTLA4. V(D)J analysis shows clonally expanded T cells, indicating that they are capable of mounting an immune response. An indolent liver metastasis from a class 1B UM is infiltrated with clonally expanded plasma cells, indicative of antibody-mediated immunity. This complex ecosystem of tumor and immune cells provides new insights into UM biology, and LAG3 is identified as a potential candidate for immune checkpoint blockade in patients with high risk UM.

PURPOSE: Validated and accurate prognostic testing is critical for precision medicine in uveal melanoma (UM). Our aims were to (1) prospectively validate an integrated prognostic classifier combining a 15-gene expression profile (15-GEP) and PRAME RNA expression and (2) identify clinical variables that enhance the prognostic accuracy of the 15-GEP/ PRAME classifier. MATERIALS AND METHODS: This study included 1,577 patients with UM of the choroid and/or ciliary body who were enrolled in the Collaborative Ocular Oncology Group Study Number 2 (COOG2) and prospectively monitored across 26 North American centers. Test results for 15-GEP (class 1 or class 2) and PRAME expression status (negative or positive) were available for all patients. The primary end point was metastasis-free survival (MFS). RESULTS: 15-GEP was class 1 in 1,082 (68.6%) and class 2 in 495 (31.4%) patients. PRAME status was negative in 1,106 (70.1%) and positive in 471 (29.9%) patients. Five-year MFS was 95.6% (95% CI, 93.9 to 97.4) for class 1/ PRAME (–), 80.6% (95% CI, 73.9 to 87.9) for class 1/ PRAME (+), 58.3% (95% CI, 51.1 to 66.4) for class 2/ PRAME (–), and 44.8% (95% CI, 37.9 to 52.8) for class 2/ PRAME (+). By multivariable Cox proportional hazards analysis, 15-GEP was the most important independent predictor of MFS (hazard ratio [HR], 5.95 [95% CI, 4.43 to 7.99]; P < .001), followed by PRAME status (HR, 1.82 [95% CI, 1.42 to 2.33]; P < .001). The only clinical variable demonstrating additional prognostic value was tumor diameter. CONCLUSION: In the largest prospective multicenter prognostic biomarker study performed to date in UM to our knowledge, the COOG2 study validated the superior prognostic accuracy of the integrated 15-GEP/ PRAME classifier over 15-GEP alone and clinical prognostic variables. Tumor diameter was found to be the only clinical variable to provide additional prognostic information. This prognostic classifier provides an advanced resource for risk-adjusted metastatic surveillance and adjuvant trial stratification in patients with UM.

Purpose: The strong association between BAP1 mutations and metastasizing Class 2 uveal melanoma (UM) suggests that epigenetic alterations may play a significant role in tumor progression. Thus, we characterized the impact of BAP1 loss on the DNA methylome in UM. Experimental Design: Global DNA methylation was analyzed in 47 Class 1 and 45 Class 2 primary UMs and in UM cells engineered to inducibly deplete BAP1. RNA-Seq was analyzed in 80 UM samples and engineered UM cells. Results: Hypermethylation on chromosome 3 correlated with downregulated gene expression at several loci, including 3p21, where BAP1 is located. Gene set analysis of hypermethylated and downregulated genes identified axon guidance and melanogenesis as deregulated pathways, with several of these genes located on chromosome 3. A novel hypermethylated site within the BAP1 locus was found in all Class 2 tumors, suggesting that BAP1 itself is epigenetically regulated. Highly differentially methylated probes were orthogonally validated using bisulfite sequencing, and they successfully distinguished Class 1 and Class 2 tumors in 100% of cases. In functional validation experiments, BAP1 knockdown in UM cells induced methylomic repatterning similar to UM tumors, enriched for genes involved in axon guidance, melanogenesis, and development. Conclusions: This study, coupled with previous work, suggests that the initial event in the divergence of Class 2 UM from Class 1 UM is loss of one copy of chromosome 3, followed by mutation of BAP1 on the remaining copy of chromosome 3, leading to the methylomic repatterning profile characteristic of Class 2 UMs.

Uveal melanoma (UM) is the most common primary cancer of the eye and is associated with a high rate of metastatic death. UM can be stratified into two main classes based on metastatic risk, with class 1 UM having a low metastatic risk and class 2 UM having a high metastatic risk. Class 2 UM have a distinctive genomic, transcriptomic, histopathologic, and clinical phenotype characterized by biallelic inactivation of the BAP1 tumor-suppressor gene, an immune-suppressive microenvironment enriched for M2-polarized macrophages, and poor response to checkpoint-inhibitor immunotherapy. To identify potential mechanistic links between BAP1 loss and immune suppression in class 2 UM, we performed an integrated analysis of UM samples, as well as genetically engineered UM cell lines and uveal melanocytes (UMC). Using RNA sequencing (RNA-seq), we found that the most highly upregulated gene associated with BAP1 loss across these datasets was PROS1, which encodes a ligand that triggers phosphorylation and activation of the immunosuppressive macrophage receptor MERTK. The inverse association between BAP1 and PROS1 in class 2 UM was confirmed by single-cell RNA-seq, which also revealed that MERTK was upregulated in CD163+ macrophages in class 2 UM. Using ChIP-seq, BAP1 knockdown in UM cells resulted in an accumulation of H3K27ac at the PROS1 locus, suggesting epigenetic regulation of PROS1 by BAP1. Phosphorylation of MERTK in RAW 264.7 monocyte–macrophage cells was increased upon coculture with BAP1^−/− UMCs, and this phosphorylation was blocked by depletion of PROS1 in the UMCs. These findings were corroborated by multicolor immunohistochemistry, where class 2/BAP1-mutant UMs demonstrated increased PROS1 expression in tumor cells and increased MERTK phosphorylation in CD163+ macrophages compared with class 1/BAP1-wildtype UMs. Taken together, these findings provide a mechanistic link between BAP1 loss and the suppression of the tumor immune microenvironment in class 2 UMs, and they implicate the PROS1–MERTK pathway as a potential target for immunotherapy in UM.

PRAME is a CUL2 ubiquitin ligase subunit that is normally expressed in the testis but becomes aberrantly overexpressed in many cancer types in association with aneuploidy and metastasis. Here, we show that PRAME is expressed predominantly in spermatogonia around the time of meiotic crossing-over in coordination with genes mediating DNA double strand break repair. Expression of PRAME in somatic cells upregulates pathways involved in meiosis, chromosome segregation and DNA repair, and it leads to increased DNA double strand breaks, telomere dysfunction and aneuploidy in neoplastic and non-neoplastic cells. This effect is mediated at least in part by ubiquitination of SMC1A and altered cohesin function. PRAME expression renders cells susceptible to inhibition of PARP1/2, suggesting increased dependence on alternative base excision repair pathways. These findings reveal a distinct oncogenic function of PRAME that can be targeted therapeutically in cancer.

Retinoblastoma (Rb) is a deadly childhood eye cancer that is classically initiated by inactivation of the RB1 tumor suppressor. Clinical management continues to rely on nonspecific chemotherapeutic agents that are associated with treatment resistance and toxicity. Here, we analyzed 103 whole exomes, 20 whole transcriptomes, 5 single-cell transcriptomes, and 4 whole genomes from primary Rb tumors to identify previously unknown Rb dependencies. Several recurrent genomic aberrations implicate estrogen-related receptor gamma (ESRRG) in Rb pathogenesis. RB1 directly interacts with and inhibits ESRRG, and RB1 loss uncouples ESRRG from negative regulation. ESRRG regulates genes involved in retinogenesis and oxygen metabolism in Rb cells. ESRRG is preferentially expressed in hypoxic Rb cells in vivo. Depletion or inhibition of ESRRG causes marked Rb cell death, which is exacerbated in hypoxia. These findings reveal a previously unidentified dependency of Rb cells on ESRRG, and they implicate ESRRG as a potential therapeutic vulnerability in Rb.