We revisited the role of PRAME in a cohort of 64 patients with UM that had been part of prior publications, to which we now added a longer follow-up, and added the analysis of the TCGA cohort of 80 cases. We additionally studied the literature about the use of PRAME as prognostic marker, and we reviewed the current evidence supporting the use of PRAME as a target for immunotherapy. We confirmed the correlation of PRAME expression with negative prognostic factors, as has previously been reported
70,91: in the total LUMC cohort, the LUMC disomy 3 cohort, and in the TCGA total cohort,
PRAME expression was correlated with a higher LBD (see
Table 2,
70 Supplementary Table S4), as also reported by Schefler and Field.
69,91 To add to PRAME's negative prognostic profile, in both the LUMC and the TCGA cohort,
PRAME expression was higher in darker tumors (see
Table 2,
Fig. 2). We did not find a significant correlation of PRAME with chromosome 3 status or BAP1 expression in either cohort, and we saw a borderline correlation with 6p gain in the TCGA disomy 3 subcohort and a matched near-significant correlation in the LUMC disomy 3 subcohort. Field et al. had previously reported an association between high PRAME expression and the presence of 6p gain, among other chromosomal alterations, such as 6q loss, 8q gain and 16q loss, in both class 1 and class 2 UM.
69 At the moment, however, because we only found this association in one of the subcohorts, it is hard to say if 6p gain itself has a role or if it is part of the genomic instability postulated by Harbour et al.
90 More consistent correlations were found with the presence of an inflammatory phenotype, with enrichment in inflammation-related pathways (see
Table 3,
Figs. 5,
6) and with gain of chromosome 8q (
Table 2,
70 Supplementary Tables S4–
S7,
S11). Our differential expression analysis highlighted 8q24.3 as strictly related to
PRAME overexpression, because most of the upregulated genes are located in this region. Interestingly, this specific region has been reported as relevant in the progression of several types of cancers, and in UM in particular.
96 Speicher et al. studied 11 patients with UM with comparative genomic hybridization (CGH) and reported gain of chromosome 8 as the most frequent copy number alteration and the segment 8q24→qter as the smallest represented segment.
158 Similarly, Anbunathan et al. analyzed data from 182 patients with UM and reported gain of 8q and loss of heterozygosity in chromosome 3 to be present in >50% of cases and identified 8q24.3 as one of the focal copy number regions associated with metastasis and worse survival.
159 A literature search into the region 8q24 highlighted some interesting genes that have been related to UM progression and prognosis. C-MYC (MYC Proto-Oncogene, BHLH Transcription Factor), a known proto-oncogene located on 8q24.1, has been studied in UM and its expression has been correlated with a high proliferation, the presence of monosomy 3 or BAP1 mutation, but its association with survival is not unanimous.
160–165 DDEF1 (ArfGAP with SH3 domain, ankyrin repeat and PH domain 1), located on 8q24.21, was identified as the gene, the expression of which correlated best with the number of 8q copies in a study involving 25 patients with primary UM, was significantly higher in class 2 UM compared to class 1 UM and its overexpression increased cell motility in the Mel202 UM cell line.
166 The current study, however, found no significant increase in C-MYC in UM with 8q gain. Protein tyrosine kinase 2 (PTK2) on 8q24.3 has been shown to be amplified in breast cancer
167 and hepatocellular carcinoma,
168 and to be among the upregulated genes that exhibited hypomethylation in a study involving 182 UM cases.
159 Protein tyrosine phosphatase 4A3 (PTP4A3), located on 8q.24.3, as mentioned previously, showed a significant association with early metastases and caused an increase in cell migration and invasiveness in UM cell lines.
117 Moreover, PTP4A3 was among the genes that showed upregulation and hypomethylation in the study by Anbunathan et al.
159 Considering this information from the literature and our data, we think a good candidate that is worth studying further, both in the context of PRAME and of UM with 8q gain in general, might be PTP4A3. In our cohort, only PTP4A3 was among the most differentially expressed genes in PRAME-positive UM (adjusted
P value = 0.026, logFC = 0.819; see
Fig. 7,
Supplementary Table S11), whereas MYC had an adjusted
P value > 0.05 and DDEF1 and PTK2 had an adjusted
P value of 0.049 and 0.041, respectively, and a lower logFC (0.643 and 0.461, respectively).