Next Article in Journal
Etiology, Pathogenesis, Diagnosis, and Practical Implications of Hepatocellular Neoplasms
Next Article in Special Issue
Selective Internal Radiotherapy (SIRT) and Chemosaturation Percutaneous Hepatic Perfusion (CS-PHP) for Metastasized Uveal Melanoma: A Retrospective Comparative Study
Previous Article in Journal
Outcomes of Minimally Invasive and Robot-Assisted Esophagectomy for Esophageal Cancer
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

FOXD1 Is a Transcription Factor Important for Uveal Melanocyte Development and Associated with High-Risk Uveal Melanoma

by
Quincy C. C. van den Bosch
1,2,†,
Josephine Q. N. Nguyen
1,2,†,
Tom Brands
1,2,
Thierry P. P. van den Bosch
3,
Robert M. Verdijk
3,4,
Dion Paridaens
5,
Nicole C. Naus
1,
Annelies de Klein
2,
Emine Kiliç
1,*,‡ and
Erwin Brosens
2,*,‡,§ on behalf of the Rotterdam Ocular Melanoma Study Group
1
Department of Ophthalmology, Erasmus MC Cancer Center, Erasmus MC University Medical Center Rotterdam, 3000 CA Rotterdam, The Netherlands
2
Department of Clinical Genetics, Erasmus MC Cancer Center, Erasmus MC University Medical Center Rotterdam, 3000 CA Rotterdam, The Netherlands
3
Department of Pathology, Section Ophthalmic Pathology, Erasmus MC Cancer Institute, Erasmus MC University Medical Center Rotterdam, 3000 CA Rotterdam, The Netherlands
4
Department of Pathology, Leiden University Medical Center, 2333 ZA Leiden, The Netherlands
5
The Rotterdam Eye Hospital, 3011 BH Rotterdam, The Netherlands
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
These authors contributed equally to this work.
§
The Rotterdam Ocular Melanoma Study Group (ROMS) is a collaborative research group with members from the Rotterdam Eye Hospital, Departments of Ophthalmology, Pathology and Clinical Genetics, of the Erasmus MC, Rotterdam, The Netherlands.
Cancers 2022, 14(15), 3668; https://doi.org/10.3390/cancers14153668
Submission received: 27 May 2022 / Revised: 12 July 2022 / Accepted: 22 July 2022 / Published: 28 July 2022

Abstract

:

Simple Summary

Despite successful treatment of primary uveal melanoma (UM), metastases still occur in approximately 50% of the patients. Unfortunately, little is known about the mechanism behind metastasized UM. By reanalyzing publicly available single-cell RNA sequencing data of embryonic zebrafish larvae and validating the results with UM data, we have identified five transcription regulators of interest: ELL2, KDM5B, REXO4, RBFOX2 and FOXD1. The most significant finding is FOXD1, which is nearly exclusively expressed in high-risk UM and is associated with poor survival. FOXD1 is a novel gene which could be involved in the metastatic capability of UM. Elucidating its function and role in metastatic UM could help to understand and develop treatment for UM.

Abstract

Uveal melanoma (UM) is a deadly ocular malignancy, originating from uveal melanocytes. Although much is known regarding prognostication in UM, the exact mechanism of metastasis is mostly unknown. Metastatic tumor cells are known to express a more stem-like RNA profile which is seen often in cell-specific embryonic development to induce tumor progression. Here, we identified novel transcription regulators by reanalyzing publicly available single cell RNA sequencing experiments. We identified five transcription regulators of interest: ELL2, KDM5B, REXO4, RBFOX2 and FOXD1. Our most significant finding is FOXD1, as this gene is nearly exclusively expressed in high-risk UM and its expression is associated with a poor prognosis. Even within the BAP1-mutated UM, the expression of FOXD1 is correlated with poor survival. FOXD1 is a novel factor which could potentially be involved in the metastatic capacity of high-risk UM. Elucidating the function of FOXD1 in UM could provide insight into the malignant transformation of uveal melanocytes, especially in high-risk UM.

Graphical Abstract

‡ These authors contributed equally to this work.
§ The Rotterdam Ocular Melanoma Study Group (ROMS) is a collaborative research group with members from the Rotterdam Eye Hospital, Departments of Ophthalmology, Pathology and Clinical Genetics, of the Erasmus MC, Rotterdam, the Netherlands.

1. Introduction

Ocular melanomas account for around 5% of all melanoma with uveal melanoma (UM) being the most prevalent. Arising from the uvea, UM shows a distinct mutation burden pattern when compared to malignant melanoma of the skin with a predominant absence of a UV mutational signature [1,2]. UM has an incidence of one to nine per million in Europe, with a lower incidence in the southern countries and a higher incidence in the northern countries [3]. Successful treatment aims to preserve the eye and vision by stereotactic radiotherapy, brachytherapy or proton therapy and enucleation in case of more advanced tumors [4]. Despite successful local tumor control, UM is a highly aggressive disease where in the end 50% of all patients develop metastases [3,5]. Whilst much is known about UM prognostication and driver mutations [6,7,8,9,10,11], mechanisms behind metastatic UM are largely unknown. Apart from the recently discovered therapy drug Tebentafusp [12], therapeutic options are still limited in patients with metastatic UM, and therefore it is crucial to investigate metastatic events in UM [13]. Tumor cells frequently take advantage of transcription factors to alter their gene expression [14], specifically during events such as epithelial-to-mesenchymal transition [15]. Which transcription factor is used often relies on the cell type the tumor originates from [14]. Despite the fact that all UMs arise from the uvea, high-risk UMs harbor an altered transcriptome compared to low-risk UMs [16,17]. In addition, others have shown that high-risk UMs also have an increased expression of transcription regulators [16,17,18]. Here, we hypothesized that high-risk UM utilizes developmental transcription regulators used in the early establishment of uveal melanocytes for its metastatic potential. However, even though there is a clear general understanding how melanocyte cell fate is programmed [19], little is known about the mechanism behind the establishment of uveal melanocytes and what separates them from cutaneous and mucosal melanocytes.
All melanocytes spread throughout the body follow a similar gene expression profile, yet animal studies suggest an independent mechanism for ocular pigmentation during development. For instance, knock-out (KO) mice models of Cell Division Cycle 42 (Cdc42), the Actin Related Protein 2/3 (Arp2-3) complex, Brahma-related gene 1 (Brg1), Myosin X (Myo10) and Rac Family Small GTPase 1 (Rac1) revealed white spotting of the skin as a result of migration or differentiation defects of cutaneous melanocytes [20,21,22,23,24]. However, pigmentation of the eyes remained normal in these models. This phenomenon is also seen in non-mammalian vertebrate models. Genetic disruption of microphthalmia-associated transcription factor A (mitfa) and roy orbison (roy) in zebrafish generated the casper zebrafish that completely lacks melanocytes and iridophores. Nevertheless, despite the genetic disruption of melanocytic core genes casper zebrafish harbor normally pigmented and functional eyes [25]. Ocular pigmentation in zebrafish development has been shown to depend on orthodenticle homolog (otx) transcription factors rather than on mitfa or microphthalmia-associated transcription factor B (mitfb) [26]. Furthermore, using RNA in situ hybridization (ISH) of pigmentation genes in zebrafish revealed genetic divergence after the teleost genome duplication. For example, the human genome Premelanosome protein (PMEL) is expressed in all melanocytes as it is necessary for pigment functionality [27], whereas zebrafish express premelanosome protein a (pmela) in all melanocytes but restrict the expression of premelanosome protein b (pmelb) to the eyes [28]. In order to gain insights in transcriptional factors involved in early uveal melanocyte development, we argue that the zebrafish is a useful model as there is available embryonic vertebrate developmental data and genes to specify melanocyte location (eye versus all melanocytes).
To identify uveal melanocyte specific transcription regulators we mined the zebrafish atlas [29], as this single-cell RNA-sequencing (scRNA-seq) dataset yields whole zebrafish larvae at different embryonic stages. We have extracted melanocyte clusters and isolated the transcription regulators that are unique to pmelb positive cells. After the zebrafish to human orthologue prediction, we inspected the expression of the identified embryonic transcription regulators in healthy human melanocytes and UM. With this approach we aimed to (1) identify transcription regulators involved in early uveal melanocyte development and to (2) assess if these genes are active in high-risk UM.

2. Materials and Methods

2.1. (sc)RNA-Seq Datasets and Identification of Melanocyte-Specific Markers Utilizing Public Datasets

scRNA-seq datasets of whole zebrafish [29] and the human eye (GSE135922) [30] were downloaded from the gene expression omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/gds/, accessed on 24 April 2021) and used for gene expression analysis (Figure 1). The Cell Ranger count matrices acquired through the GEO Accession viewer were processed into count tables using the Seurat package [31] in R.
Analysis of public Single Cell datasets was performed using the Qiagen CLC Genomics Workbench (version 21.0.5), CLC Single Cell Analysis Module (Version 21.1) and Seurat (version 4.0). To import the data, the Cellranger Matrix exports (MEX) were used where possible, either directly or through conversion of other formats to MEX using Seurat. Where there was no compatible count data available, the reads were imported using the Workbench’s SRA download function. The Single Cell plugin comes with default pipelines for both count data and reads. The read workflow has three additional steps: annotating the reads with Cell and UMI information, trimming them and mapping them against a reference genome. Then, both pipelines continue with a quality control step before creating a normalized matrix of all samples in a dataset. All datasets were processed using the default settings and the genomes available in the CLC Genomics Workbench (Homo Sapiens Hg19 Ensembl version 75 & Danio Rerio GRCz11). The normalized matrices were then clustered based on 5000 highly variable genes, principal component analysis and Pearson’s correlation for measuring distance. Cluster data was grouped by the Leiden algorithm (resolution 0.1) and visualized as a UMAP. Additionally, clusters were sub-setted and exported for further analysis. Gene expression thresholding was determined by overall gene count distribution to include a multimodal binomial distribution of gene counts while excluding low abundant genes. This threshold was set to be at least 100 gene counts (Figure S1). The DRSC integrative ortholog prediction tool (DIOPT; Version 8.5) was used to identify human orthologues [32]. Human orthologues with high ranking were included for further analysis. Human orthologue genes were uploaded into the Ingenuity Pathway Analysis (IPA) software (Qiagen Bioinformatics, Spring release April 2022, QIAGEN Inc., https://digitalinsights.qiagen.com/IPA, accessed on 27 May 2022) [33]. IPA-curated knowledge was used to identify gene function to eventually perform a Core Analysis using default parameters, to predict the effects on pathways and biological functions.

2.2. Bulk RNA-Sequencing Dataset Analysis of the TCGA and ROMS Uveal Melanoma Cohort

The UM RNA-seq data of the 80 UM samples in the TCGA database was obtained from the Genomic Data Common (GDC) Portal (https://portal.gdc.cancer.gov/, accessed on 18 June 2021). We used whole-transcriptome sequencing data from 26 UM from the ROMS-cohort, which has been described previously by Smit et al. [17]. Our ethics committee and current informed consent does not allow sharing of individual patient or control genotype information in the public domain. Upon reasonable request, data access can be granted by the data access committee of the Department of Clinical Genetics of the Erasmus MC Rotterdam, the Netherlands. The mRNA library was sequenced on the Ion Proton sequencer (ThermoFisher Scientific, Waltham, MA, USA). Reads were trimmed using Cutadapt Version 3.4 [34], aligned to the human reference genome hg19 using HISAT2 version 2.1.0 [35] and aligned reads were counted using htseq-count Version 0.9.1 [36] and normalized using DESeq2. Trimmed mean per million (TMM) values were used to normalize sequencing depth across samples. For each gene, counts per million (CPM) were calculated. Whole-transcriptome samples (ROMS and TCGA) were included in the analysis when mutational status and survival data were available. Identified transcription regulators were then inspected in both cohorts with regards to their mutational status (EIF1AX-, SF3B1- or BAP1-mutation).

2.3. Immunohistochemistry (IHC) of FOXD1

Immunohistochemistry was performed with an automated, validated, and accredited staining system (Ventana Benchmark ULTRA, Ventana Medical Systems, Tucson, AZ, USA) using an ultraview Universal Alkaline Phosphatase Red Detection Kit (Ventana reference no. 760-501). Following deparaffinization and heat-induced antigen retrieval (CC1 for 32 min), the tissue samples were incubated for 32 min with the FOXD1 antibody (1:800 dilution, Abcam, Waltham, MA, USA, #129324). Counterstain was done by hematoxylin II stain for 12 min and a blue coloring reagent for 8 min according to the manufacturer’s instructions (Ventana Benchmark ULTRA, Ventana Medical Systems, Tucson, AZ, USA). For each slide, healthy kidney tissue was used as a positive control.

2.4. Statistical Analysis

Statistical analysis was performed using GraphPad Prism 9.0 (GraphPad Software, Inc., San Diego, CA, USA). Differences between mutational subtype were assessed using a one-way analysis of variance (ANOVA) with the Kruskal–Wallis test and Dunn’s multiple comparisons test. The survival analysis was generated using the Kaplan–Meier method and use of the Log-Rank test to find differences between the mutational subtypes. A p-value of <0.05 was considered statistically significant.

3. Results

3.1. Identification of Ocular Melanocytic Clusters in Whole Zebrafish Single-Cell RNA Sequencing

Zebrafish melanocytes were identified based on established core melanocytic markers: mitfa, dopachrome tautomerase (dct), tyrosinase (tyr), tyrosinase-related protein 1b (tyrp1b) and pmela. We analyzed the zebrafish atlas per developmental stage (24 h post fertilization (hpf), 48 hpf and 120 hpf, Figure 2a–c). Across three uniform manifold approximation and projections (UMAPs) we identified a total of seven melanocytic clusters that express the core factors with a sufficient amount of cells (>100 cells). Next, we separated ocular melanocytes from other melanocytes based on pmelb expression which revealed four pmelb+ clusters and three pmelb- clusters. To gain insight into ocular melanocyte development we extracted the read-counts per identified cluster (pmelb+ cluster in red, pmelb- clusters in blue, Figure 1e–g). In order to investigate the expression of a gene, we determined our threshold to have at least a total sum of 100 read counts based on the overall gene-count distribution. With this threshold we included a multimodal binomial distribution of gene counts while filtering out low-abundance genes (Figure S1). After setting our threshold, we extracted all genes expressed per melanocytic cluster and combined the different developmental stages to end up with a pmelb- and a pmelb+ gene list (Tables S1 and S2). We compared both clusters together, where core melanocytic genes are expressed in both clusters and pmelb is discriminatory between the clusters (Figure 3). A total of 279 genes were unique to pmelb- cells, 1817 genes are expressed in both clusters and 637 genes were unique to pmelb+ cells.

3.2. Orthologue Prediction, Function Annotation and Expression Validation in Human Melanocytes and UM

To validate 637 uniquely expressed genes in pmelb+ zebrafish cells in human melanocytes, we performed the DIOPT orthologue prediction [32]. A total of 478 genes were predicted to their human orthologues with high ranking scores (Table S3). Next, we annotated gene function and identified associated pathways with this gene list using IPA analysis (Table S4). We identified 49 genes involved in transcriptional regulation (10.2%), 107 enzymes (22.0%), 233 genes classified as ‘other’ (56.8%), and the remaining genes were spread out in multiple functions (Figure 4a). Pathway analysis shows association with cellular development, embryonic development and cancer. Next, we inspected the expression of all transcription regulators in healthy human melanocytes (GSE135922) [30] using the same scRNA-seq pipeline and thresholding as done before. In total we identified the expression of 21 transcription regulators in healthy melanocytes, illustrating that these genes are expressed in ocular melanocytes (Figure 4b). In addition, the expression of the 49 transcription regulators was analyzed in RNA-seq data from the The Cancer Genome Atlas (TCGA) and our Rotterdam Ocular Melanoma Study (ROMS) cohort (Figure 4c). We averaged the gene expression in CPM by combining all samples and found that also in UM the majority of transcription regulators are expressed. However, to gain insights into factors used in high-risk UM, we looked into the expression pattern per secondary driver mutation (BRCA1 associated protein 1 (BAP1), Splicing factor 3b Subunit 1 (SF3B1) and Eukaryotic Translation Initiation Factor 1A X-Linked (EIF1AX)).

3.3. RNA Expression Analysis of Transcription Regulator Genes of Interest in UM

We assessed the expression levels of transcription regulator genes of interest in UM in our own and publicly available UM patient cohorts (n = 26; ROMS-cohort, n = 56; TCGA cohort) of which the secondary driver gene status and survival was known [17]. After normalization and alignment, the CPM values of all genes were extracted and compared between UM subtypes. Gene expression was grouped into driver mutations (EIF1AX, n = 16; SF3B1, n = 27; BAP1, n = 39). Out of 49 transcription regulators we identified five genes associated with high-risk UM and poor prognosis. A negative correlation was found in Elongation Factor For RNA Polymerase II 2 (ELL2), Lysine Demethylase 5B (KDM5B) and REX4 Homolog, 3′-5′ Exonuclease (REXO4), where low expression is associated with poor prognosis (Figure 5 and Figure S2). Positive correlations were found by increased expression of RNA Binding Fox-1 Homolog 2 (RBFOX2) and Forkhead Box D1 (FOXD1) within the BAP1-mutated UM, which are associated with poor prognosis (Figure 5). Interestingly, FOXD1 is significantly expressed within the BAP1-mutated UM compared to SF3B1-mutated or EIF1AX-mutated UM (Figure 5c). Next, we assessed if expression of FOXD1 is correlated with overall survival of UM patients. In both cohorts, expression of FOXD1 (CPM ≥ 1) is significantly correlated with a poor prognosis. The median survival of FOXD1-positive UM is 26.5 months when compared to 129.6 months in the FOXD1-negative UM group (Figure 5d, p = <0.0001). Due to the striking expression of FOXD1 specifically in the BAP1-mutated UM, we investigated if overall survival based on FOXD1 is due to BAP1 co-segregation (n = 39). Surprisingly, even within the BAP1-mutated UM, the expression of FOXD1 is significantly correlated with a poor prognosis (Figure 5d, p = 0.0413). BAP1-mutated UMs with FOXD1 expression have a median survival of 26.1 months when compared to 41.7 months in BAP1-mutated FOXD1 negative UM. To validate our findings on protein levels, we next assessed FOXD1 levels on UM with immunohistochemistry.

3.4. FOXD1, ELL2, RBFOX2, KDM5B and REXO4 as and RNA-Based Biomarkers in Uveal Melanoma

As proof of principle that our identified genes are able to distinguish low-risk from high-risk UM, we extracted the 15 genes used in the gene expression profiling (GEP) test [37] and plotted them together with FOXD1, ELL2, RBFOX2, KDM5B and REXO4. The ROMS and TCGA cohort were sorted based on disomy 3 or monosomy 3 status, which, in line with published data, showed a clear low-risk UM (class 1, disomy 3) and high-risk UM (class 2, monosomy 3) based on the GEP test. Upregulation of FOXD1 and RBFOX2 clusters was similar to the commercially used upregulated genes CDH1, ECM1, HTR2B and RAB31 in class 2 UM. The GEP test also consists of downregulated genes in class 2 UM, where ELL2, KDM5B and REXO4 also cluster in a similar fashion (Figure 6).

3.5. Immunohistochemistry (IHC) of FOXD1 in Uveal Melanoma

For IHC analysis, we analyzed 30 UM samples of the ROMS cohort with confirmed mutational status (EIF1AX, n = 8; SF3B1, n = 9; BAP1, n = 9, EIF1AX+SF3B1, n = 2; wild-type (WT), n = 2) and matched RNA expression (discovery cohort). We were unable to stain the full RNA-seq cohort due to lack of tissue for RNA–protein confirmation of the same patient. All samples were stained with FOXD1 antibody (Figure 7a,b) and samples were scored as either positive or negative. The expression of FOXD1 on protein level proved to be nearly correlated with survival (Figure S3a, p = 0.1472). To strengthen this analysis, we have added 29 additional UM samples (replication cohort) with confirmed mutational status (EIF1AX, n = 9; BAP1, n = 20). The replication cohort alone also showed a near correlation between FOXD1 expression and survival (Figure S3b, p = 0.2340). However, when we combined both cohorts (discovery and replication cohorts, n = 59), a significant correlation between FOXD1 expression on protein level and poorer survival was assessed (Figure 7c, p = 0.0439). As seen in other studies [38], we found FOXD1 expression to be focused mainly in the cytoplasm.

4. Discussion

In this study, we identified multiple transcription regulators which are associated with early uveal melanocyte development and high-risk UM. We have determined uveal melanocytic clusters with a unique expression pattern when compared to cutaneous melanocytes. Expression of pmelb was unique to this melanocytic cluster when compared to other melanocytic clusters. To support our hypothesis that this cluster contains uveal melanocytes, the RNA ISH of pmelb elucidated spatial expression restricted to the eye, unlike pmela [28]. Due to the teleost duplication zebrafish underwent, the diversification of pmelb can be an interesting target to study ocular pigmentation without targeting cutaneous and/or mucosal melanocytes in zebrafish.
Downregulated expression of transcription regulators in high-risk UM could help identify potential tumor-suppressor genes, however only the downregulation of ELL2 is associated with prostate cancer and multiple myeloma based on our current knowledge [39,40,41,42]. KDM4B is a well-studied lysine-specific histone demethylase that is known to be overexpressed in multiple cancers [43,44,45,46]. Our findings are discordant with the literature, where a lower expression is associated with a poor prognosis (Figure S2). Our finding could be due to BAP1 co-segregation, creating a bias in the prognosis prediction. Unlike KDM5B, REXO4 is much less well-investigated in cancer biology. Current literature illustrates overexpression in hepatocellular carcinoma that is associated with tumor progression and immune infiltration [47,48]. Further investigations are needed to unravel the role of KDM5B and REXO4 loss in high-risk UM.
We identified the upregulation of two transcription regulators, RBFOX2 and FOXD1, in high-risk UM. RBFOX2 is known to be involved in cerebellar and heart development [49,50], and FOXD1 is involved in retina development [51,52], but to our knowledge neither have been associated with melanocyte development prior to this study. RBFOX2 is a well-known regulator of epithelial-to-mesenchymal transition [53], where the alternate splicing of pre-mRNA alters signal transduction via a wide range of genes depending on cellular context [54,55,56]. Further investigation of RBFOX2 overexpression in UM could reveal novel targets for drug intervention of high-risk UM.
Our most significant finding, FOXD1, was identified in 2–5 day-old zebrafish larvae; yet this was not found in healthy human uveal melanocytes (Figure 4b). We hypothesized that these transcription factors are lost after the establishment of uveal melanocytes. The human dataset consisted of uveal melanocytes of healthy adults (age range: 54–92 years) which doesn’t reflect embryonic development, whereas the zebrafish dataset (2–5 days post fertilization) showed the embryonic stage of vertebrate development. The differences in developmental stages could be a reason as to why FOXD1 was not expressed in the human dataset. scRNA-seq data of human fetal uveal melanocytes could confirm this hypothesis; however, these data are not available to our knowledge. We reanalyzed the data of IPSC-derived melanocytes during neural crest cell generation. FOXD1 was not expressed in these in vitro models, mimicking early melanocyte development (data not shown). In addition, we inspected the expression of FOXD1 in early fetal development and FOXD1 is only expressed in skeletal muscle, retinal progenitors, Müller glia, photoreceptor and ganglion cells [57]. In this stage, there are no melanocyte cells present in the eye. Additionally, we inspected the recently published ocular dataset using the interactive single-cell portal (https://singlecell.broadinstitute.org/single_cell; accessed on 29 June 2022) [58]. In adulthood, a few ocular melanocytes show expression of FOXD1, whereas the RPE and fibroblasts show a strong expression of FOXD1 (Figure S4). Aggressive tumors are known to have a higher degree of dedifferentiation, which has been correlated to reactivate expression of developmental transcription factors [59,60,61]. Therefore, we assessed the association between FOXD1 as a transcription regulator and its potential role in UM. Analyzing two independent datasets, we determined FOXD1 to be associated with poor prognosis in UM based on RNA expression. Interestingly, we identified FOXD1 as a biomarker that divides BAP1-mutated UM into two distinct groups. Previous studies have illustrated prognostic differences within the highly aggressive UM tumors, classifying monosomy 3 UM tumors into two groups [18]. Recently, loss of BAP1 has been shown to be associated with DNA methylation repatterning needed for UM to gain malignant potential [62]. Although there is no functional proof that FOXD1 and BAP1 interact, due to the distinct FOXD1 expression pattern in UM it might be involved in diverting from a low-risk to high-risk tumor. In line with previous published gene expression profiles, we illustrate that the identified transcription regulators follow a similar gene expression profile as seen in the GEP-test (Figure 6). To prove that FOXD1 is either a driver or bystander in high-risk UM, functional assays must be performed to illustrate its mechanism. However, in this article we provide proof that FOXD1 can reliably be used as a biomarker in UM. FOXD1 expression has been assessed in multiple cancer types and is often correlated to a poor prognosis [63,64,65,66]. In vitro and in vivo experiments elucidated FOXD1 as an oncogene involved in proliferation, dedifferentiation, migration, and radio- and chemoresistance [67,68,69,70,71]. For instance, in oral squamous cell carcinoma, FOXD1 promotes the epithelial–mesenchymal transition (EMT) by activation of SNAI2 [72]. Functional work on FOXD1 in cutaneous melanoma elucidated its capability to induce the expression of the tumor-specific isoform Rac Family Small GTPase 1B (RAC1B) to enhance melanoma migration and invasion [38]. Additionally, upon BRAF inhibitor treatment, FOXD1 was shown to directly bind the CTGF promoter to promote dedifferentiation and therapy resistance [64]. In line with published data on EMT, we assessed FOXD1 expression with known EMT inducers and repressors described in squamous cell carcinoma and cutaneous melanoma. However, unlike in squamous cell carcinoma and cutaneous melanoma, these EMT-like patterns are not clearly or significantly correlated with FOXD1 in UM (Figures S4a and S5b). It seems that EMT-like patterns of UM are different than what is currently been described in literature concerning the downstream effects of FOXD1. For instance, a loss of CDH1 has been describe to be necessary for EMT [73], while in the GEP test, CDH1 is used as an upregulated biomarker for class 2 UM.

5. Conclusions

In this study we found the expression of FOXD1 to be almost exclusively in BAP1-mutated UM. BAP1 mutations are associated with early metastatic disease [74,75], although its underlying mechanism remains unclear. Our analysis shows that most BAP1-mutated UMs with FOXD1 expression are correlated with poor survival. This correlation was also seen using immunohistochemistry, although the correlation was not as strong as on the RNA level. To understand what role FOXD1 plays in UM, it is necessary to perform functional assays. Due to the striking correlation of FOXD1 expression in BAP1-mutated UM, functional assays elucidating its function and targets could provide therapeutic targets for high-risk UM.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/cancers14153668/s1. Figure S1. Determination of read counts threshold; Figure S2. Expression and survival plots of KDM5B and REXO4 in the combined cohort; Figure S3. Survival plot of FOXD1 expression on IHC in the discovery cohort and validation cohort; Figure S4. UMAP of ocular compartments with FOXD1 expression highlighted; Figure S5. Heatmap of FOXD1 with EMT inducers and repressors; Table S1. Extracted gene list; Table S2. Gene expression data; Table S3. Full list of DIOPT analysis; Table S4. IPA function annotation.

Author Contributions

Conceptualization, E.K. and E.B.; methodology, T.B., Q.C.C.v.d.B. and J.Q.N.N.; software, T.B., Q.C.C.v.d.B. and J.Q.N.N.; validation, Q.C.C.v.d.B., J.Q.N.N. and T.B.; formal analysis, T.B., Q.C.C.v.d.B. and J.Q.N.N.; investigation, T.B., Q.C.C.v.d.B. and J.Q.N.N.; resources, E.K. and E.B.; data curation, T.B., Q.C.C.v.d.B. and J.Q.N.N.; writing—original draft preparation, Q.C.C.v.d.B. and J.Q.N.N.; writing—review and editing, E.K., T.B., T.P.P.v.d.B., R.M.V., D.P., N.C.N., A.d.K., E.B. and E.K.; supervision, A.d.K., E.K. and E.B.; project administration, E.K. and E.B.; funding acquisition, E.K. and E.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Henkes foundation (2021-04) J.Q.N.N., E.K.; Rotterdamse Blindenbelangen (HV/AB/B20200043) Q.C.C.v.d.B., E.B.; Collaborative Ophthalmic Research Rotterdam (6.2.0) J.Q.N.N., E.K.; Landelijke Stichting voor Blinden en Slechtzienden (2019-7) Q.C.C.v.d.B., E.B.; and Algemene Nederlandse Vereniging ter voorkoming van Blindheid (2019-07) Q.C.C.v.d.B., E.B.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Medical ethics committee of the Erasmus Medical Centre (OZR nr 2009-17, MEC-2009-375, 12 November 2009).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data from The Cancer Genome Atlas (TCGA) and repository NCBI GEO Datasets are publicly accessible. Our ethics committee does not allow sharing of individual patient or control genotype information in the public domain.

Acknowledgments

The results shown here are in part based upon data generated by the TCGA Research Network: https://www.cancer.gov/tcga (accessed on 27 May 2022).

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Abbreviations

Arp2-3Actin Related Protein 2/3
BAP1BRCA1 associated protein 1
Brg1Brahma-related gene 1
Cdc42Cell Division Cycle 42
CDH1Cadherin 1
CPMCounts per mean
dctDopachrome tautomerase
DIOPTDRSC integrative ortholog prediction tool
ECM1Extracellular Matrix Protein 1
EIF1AXEukaryotic Translation Initiation Factor 1A X-Linked
ELL2Elongation Factor For RNA Polymerase II 2
EMTEpithelial-mesenchymal transition
FOXD1Forkhead box D1
GEOGene expression omnibus
GEPGene expression profiling
hpfHours post fertilization
HTR2B5-Hydroxytryptamine receptor 2B
IHCImmunohistochemistry
IPAIngenuity pathway analysis
IPSCInduced pluripotent stem cells
ISHIn situ hybridisation
KDM5BLysine Demethylase 5B
KOKnock-out
mitfaMicrophthalmia-associated transcription factor A
mitfbMicrophthalmia-associated transcription factor B
Myo10Myosin X
OTXOrthodenticle Homeobox
PMELPremelanosome protein
pmelaPremelanosome protein A
pmelbPremelanosome protein B
RAB31Ras-related protein Rab-31
Rac1Rac Family Small GTPase 1
RAC1BRac Family Small GTPase 1B
RBFOX2RNA Binding Fox-1 Homolog 2
REXO4REX4 Homolog, 3′-5′ Exonuclease
ROMSRotterdam Ocular Melanoma Study
royRoy Orbison
scRNA-seqSingle-cell RNA sequencing
SF3B1Splicing factor 3b Subunit 1
SNAI2Snail Family Transcriptional Repressor 2
TCGAThe Cancer Genome Atlas
tyrp1bTyrosinase related protein 1B
UMUveal melanoma
UMAPUniform manifold approximation and projection
WTWild-type

References

  1. Pleasance, E.D.; Cheetham, R.K.; Stephens, P.J.; McBride, D.J.; Humphray, S.J.; Greenman, C.D.; Varela, I.; Lin, M.L.; Ordonez, G.R.; Bignell, G.R.; et al. A comprehensive catalogue of somatic mutations from a human cancer genome. Nature 2010, 463, 191–196. [Google Scholar] [CrossRef] [PubMed]
  2. Royer-Bertrand, B.; Torsello, M.; Rimoldi, D.; El Zaoui, I.; Cisarova, K.; Pescini-Gobert, R.; Raynaud, F.; Zografos, L.; Schalenbourg, A.; Speiser, D.; et al. Comprehensive Genetic Landscape of Uveal Melanoma by Whole-Genome Sequencing. Am. J. Hum. Genet. 2016, 99, 1190–1198. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Jager, M.J.; Shields, C.L.; Cebulla, C.M.; Abdel-Rahman, M.H.; Grossniklaus, H.E.; Stern, M.H.; Carvajal, R.D.; Belfort, R.N.; Jia, R.; Shields, J.A.; et al. Uveal melanoma. Nat. Rev. Dis. Primers 2020, 6, 24. [Google Scholar] [CrossRef] [PubMed]
  4. Singh, A.D.; Turell, M.E.; Topham, A.K. Uveal melanoma: Trends in incidence, treatment, and survival. Ophthalmology 2011, 118, 1881–1885. [Google Scholar] [CrossRef] [PubMed]
  5. Kujala, E.; Makitie, T.; Kivela, T. Very long-term prognosis of patients with malignant uveal melanoma. Investig. Ophthalmol. Vis. Sci. 2003, 44, 4651–4659. [Google Scholar] [CrossRef] [Green Version]
  6. Van Raamsdonk, C.D.; Griewank, K.G.; Crosby, M.B.; Garrido, M.C.; Vemula, S.; Wiesner, T.; Obenauf, A.C.; Wackernagel, W.; Green, G.; Bouvier, N.; et al. Mutations in GNA11 in uveal melanoma. N. Engl. J. Med. 2010, 363, 2191–2199. [Google Scholar] [CrossRef] [Green Version]
  7. Koopmans, A.E.; Verdijk, R.M.; Brouwer, R.W.; van den Bosch, T.P.; van den Berg, M.M.; Vaarwater, J.; Kockx, C.E.; Paridaens, D.; Naus, N.C.; Nellist, M.; et al. Clinical significance of immunohistochemistry for detection of BAP1 mutations in uveal melanoma. Mod. Pathol. 2014, 27, 1321–1330. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Martin, M.; Masshofer, L.; Temming, P.; Rahmann, S.; Metz, C.; Bornfeld, N.; van de Nes, J.; Klein-Hitpass, L.; Hinnebusch, A.G.; Horsthemke, B.; et al. Exome sequencing identifies recurrent somatic mutations in EIF1AX and SF3B1 in uveal melanoma with disomy 3. Nat. Genet. 2013, 45, 933–936. [Google Scholar] [CrossRef] [Green Version]
  9. Yavuzyigitoglu, S.; Koopmans, A.E.; Verdijk, R.M.; Vaarwater, J.; Eussen, B.; van Bodegom, A.; Paridaens, D.; Kilic, E.; de Klein, A.; Rotterdam Ocular Melanoma Study Group. Uveal Melanomas with SF3B1 Mutations: A Distinct Subclass Associated with Late-Onset Metastases. Ophthalmology 2016, 123, 1118–1128. [Google Scholar] [CrossRef]
  10. Johansson, P.; Aoude, L.G.; Wadt, K.; Glasson, W.J.; Warrier, S.K.; Hewitt, A.W.; Kiilgaard, J.F.; Heegaard, S.; Isaacs, T.; Franchina, M.; et al. Deep sequencing of uveal melanoma identifies a recurrent mutation in PLCB4. Oncotarget 2016, 7, 4624–4631. [Google Scholar] [CrossRef] [Green Version]
  11. Moore, A.R.; Ceraudo, E.; Sher, J.J.; Guan, Y.; Shoushtari, A.N.; Chang, M.T.; Zhang, J.Q.; Walczak, E.G.; Kazmi, M.A.; Taylor, B.S.; et al. Recurrent activating mutations of G-protein-coupled receptor CYSLTR2 in uveal melanoma. Nat. Genet. 2016, 48, 675–680. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Nathan, P.; Hassel, J.C.; Rutkowski, P.; Baurain, J.F.; Butler, M.O.; Schlaak, M.; Sullivan, R.J.; Ochsenreither, S.; Dummer, R.; Kirkwood, J.M.; et al. Overall Survival Benefit with Tebentafusp in Metastatic Uveal Melanoma. N. Engl. J. Med. 2021, 385, 1196–1206. [Google Scholar] [CrossRef] [PubMed]
  13. Jochems, A.; van der Kooij, M.K.; Fiocco, M.; Schouwenburg, M.G.; Aarts, M.J.; van Akkooi, A.C.; van den Berkmortel, F.; Blank, C.U.; van den Eertwegh, A.J.M.; Franken, M.G.; et al. Metastatic Uveal Melanoma: Treatment Strategies and Survival-Results from the Dutch Melanoma Treatment Registry. Cancers 2019, 11, 1007. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Reddy, J.; Fonseca, M.A.S.; Corona, R.I.; Nameki, R.; Segato Dezem, F.; Klein, I.A.; Chang, H.; Chaves-Moreira, D.; Afeyan, L.K.; Malta, T.M.; et al. Predicting master transcription factors from pan-cancer expression data. Sci. Adv. 2021, 7, eabf6123. [Google Scholar] [CrossRef] [PubMed]
  15. Stemmler, M.P.; Eccles, R.L.; Brabletz, S.; Brabletz, T. Non-redundant functions of EMT transcription factors. Nat. Cell Biol. 2019, 21, 102–112. [Google Scholar] [CrossRef]
  16. Smit, K.N.; Boers, R.; Vaarwater, J.; Boers, J.; Brands, T.; Mensink, H.; Verdijk, R.M.; van IJcken, W.F.J.; Gribnau, J.; de Klein, A.; et al. Genome-wide aberrant methylation in primary metastatic UM and their matched metastases. Sci. Rep. 2022, 12, 42. [Google Scholar] [CrossRef] [PubMed]
  17. Smit, K.N.; Chang, J.; Derks, K.; Vaarwater, J.; Brands, T.; Verdijk, R.M.; Wiemer, E.A.C.; Mensink, H.W.; Pothof, J.; de Klein, A.; et al. Aberrant MicroRNA Expression and Its Implications for Uveal Melanoma Metastasis. Cancers 2019, 11, 815. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Robertson, A.G.; Shih, J.; Yau, C.; Gibb, E.A.; Oba, J.; Mungall, K.L.; Hess, J.M.; Uzunangelov, V.; Walter, V.; Danilova, L.; et al. Integrative Analysis Identifies Four Molecular and Clinical Subsets in Uveal Melanoma. Cancer Cell 2017, 32, 204–220.e15. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Mort, R.L.; Jackson, I.J.; Patton, E.E. The melanocyte lineage in development and disease. Development 2015, 142, 620–632. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Woodham, E.F.; Paul, N.R.; Tyrrell, B.; Spence, H.J.; Swaminathan, K.; Scribner, M.R.; Giampazolias, E.; Hedley, A.; Clark, W.; Kage, F.; et al. Coordination by Cdc42 of Actin, Contractility, and Adhesion for Melanoblast Movement in Mouse Skin. Curr. Biol. 2017, 27, 624–637. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Papalazarou, V.; Swaminathan, K.; Jaber-Hijazi, F.; Spence, H.; Lahmann, I.; Nixon, C.; Salmeron-Sanchez, M.; Arnold, H.H.; Rottner, K.; Machesky, L.M. The Arp2/3 complex is crucial for colonisation of the mouse skin by melanoblasts. Development 2020, 147, dev194555. [Google Scholar] [CrossRef] [PubMed]
  22. Marathe, H.G.; Watkins-Chow, D.E.; Weider, M.; Hoffmann, A.; Mehta, G.; Trivedi, A.; Aras, S.; Basuroy, T.; Mehrotra, A.; Bennett, D.C.; et al. BRG1 interacts with SOX10 to establish the melanocyte lineage and to promote differentiation. Nucleic Acids Res. 2017, 45, 6442–6458. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Tokuo, H.; Bhawan, J.; Coluccio, L.M. Myosin X is required for efficient melanoblast migration and melanoma initiation and metastasis. Sci. Rep. 2018, 8, 10449. [Google Scholar] [CrossRef] [PubMed]
  24. Li, A.; Ma, Y.; Yu, X.; Mort, R.L.; Lindsay, C.R.; Stevenson, D.; Strathdee, D.; Insall, R.H.; Chernoff, J.; Snapper, S.B.; et al. Rac1 drives melanoblast organization during mouse development by orchestrating pseudopod- driven motility and cell-cycle progression. Dev. Cell 2011, 21, 722–734. [Google Scholar] [CrossRef] [Green Version]
  25. White, R.M.; Sessa, A.; Burke, C.; Bowman, T.; LeBlanc, J.; Ceol, C.; Bourque, C.; Dovey, M.; Goessling, W.; Burns, C.E.; et al. Transparent adult zebrafish as a tool for in vivo transplantation analysis. Cell Stem Cell 2008, 2, 183–189. [Google Scholar] [CrossRef] [Green Version]
  26. Lane, B.M.; Lister, J.A. Otx but not Mitf transcription factors are required for zebrafish retinal pigment epithelium development. PLoS ONE 2012, 7, e49357. [Google Scholar] [CrossRef] [Green Version]
  27. Watt, B.; van Niel, G.; Raposo, G.; Marks, M.S. PMEL: A pigment cell-specific model for functional amyloid formation. Pigment Cell Melanoma Res. 2013, 26, 300–315. [Google Scholar] [CrossRef] [Green Version]
  28. Cechmanek, P.B.; McFarlane, S. Retinal pigment epithelium expansion around the neural retina occurs in two separate phases with distinct mechanisms. Dev. Dyn. 2017, 246, 598–609. [Google Scholar] [CrossRef] [Green Version]
  29. Farnsworth, D.R.; Saunders, L.M.; Miller, A.C. A single-cell transcriptome atlas for zebrafish development. Dev. Biol. 2020, 459, 100–108. [Google Scholar] [CrossRef]
  30. Voigt, A.P.; Mulfaul, K.; Mullin, N.K.; Flamme-Wiese, M.J.; Giacalone, J.C.; Stone, E.M.; Tucker, B.A.; Scheetz, T.E.; Mullins, R.F. Single-cell transcriptomics of the human retinal pigment epithelium and choroid in health and macular degeneration. Proc. Natl. Acad. Sci. USA 2019, 116, 24100–24107. [Google Scholar] [CrossRef] [Green Version]
  31. Hao, Y.; Hao, S.; Andersen-Nissen, E.; Mauck, W.M., 3rd; Zheng, S.; Butler, A.; Lee, M.J.; Wilk, A.J.; Darby, C.; Zager, M.; et al. Integrated analysis of multimodal single-cell data. Cell 2021, 184, 3573–3587.e29. [Google Scholar] [CrossRef]
  32. Hu, Y.; Flockhart, I.; Vinayagam, A.; Bergwitz, C.; Berger, B.; Perrimon, N.; Mohr, S.E. An integrative approach to ortholog prediction for disease-focused and other functional studies. BMC Bioinform. 2011, 12, 357. [Google Scholar] [CrossRef] [Green Version]
  33. Kramer, A.; Green, J.; Pollard, J., Jr.; Tugendreich, S. Causal analysis approaches in Ingenuity Pathway Analysis. Bioinformatics 2014, 30, 523–530. [Google Scholar] [CrossRef]
  34. Martin, M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 2011, 17, 10–12. [Google Scholar] [CrossRef]
  35. Kim, D.; Langmead, B.; Salzberg, S.L. HISAT: A fast spliced aligner with low memory requirements. Nat. Methods 2015, 12, 357–360. [Google Scholar] [CrossRef] [Green Version]
  36. Anders, S.; Pyl, P.T.; Huber, W. HTSeq—A Python framework to work with high-throughput sequencing data. Bioinformatics 2015, 31, 166–169. [Google Scholar] [CrossRef] [PubMed]
  37. Harbour, J.W.; Chen, R. The DecisionDx-UM Gene Expression Profile Test Provides Risk Stratification and Individualized Patient Care in Uveal Melanoma. PLoS Curr. 2013, 5. [Google Scholar] [CrossRef] [PubMed]
  38. Wu, H.; Larribere, L.; Sun, Q.; Novak, D.; Sachindra, S.; Granados, K.; Umansky, V.; Utikal, J. Loss of neural crest-associated gene FOXD1 impairs melanoma invasion and migration via RAC1B downregulation. Int. J. Cancer 2018, 143, 2962–2972. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Wang, Z.; Pascal, L.E.; Chandran, U.R.; Chaparala, S.; Lv, S.; Ding, H.; Qi, L.; Wang, Z. ELL2 Is Required for the Growth and Survival of AR-Negative Prostate Cancer Cells. Cancer Manag. Res. 2020, 12, 4411–4427. [Google Scholar] [CrossRef] [PubMed]
  40. Wang, Z.; Zhong, M.; Song, Q.; Pascal, L.E.; Yang, Z.; Wu, Z.; Wang, K.; Wang, Z. Anti-apoptotic factor Birc3 is up-regulated by ELL2 knockdown and stimulates proliferation in LNCaP cells. Am. J. Clin. Exp. Urol. 2019, 7, 223–231. [Google Scholar] [PubMed]
  41. Ali, M.; Ajore, R.; Wihlborg, A.K.; Niroula, A.; Swaminathan, B.; Johnsson, E.; Stephens, O.W.; Morgan, G.; Meissner, T.; Turesson, I.; et al. The multiple myeloma risk allele at 5q15 lowers ELL2 expression and increases ribosomal gene expression. Nat. Commun. 2018, 9, 1649. [Google Scholar] [CrossRef] [PubMed]
  42. Qiu, X.; Pascal, L.E.; Song, Q.; Zang, Y.; Ai, J.; O’Malley, K.J.; Nelson, J.B.; Wang, Z. Physical and Functional Interactions between ELL2 and RB in the Suppression of Prostate Cancer Cell Proliferation, Migration, and Invasion. Neoplasia 2017, 19, 207–215. [Google Scholar] [CrossRef]
  43. Dai, B.; Hu, Z.; Huang, H.; Zhu, G.; Xiao, Z.; Wan, W.; Zhang, P.; Jia, W.; Zhang, L. Overexpressed KDM5B is associated with the progression of glioma and promotes glioma cell growth via downregulating p21. Biochem. Biophys. Res. Commun. 2014, 454, 221–227. [Google Scholar] [CrossRef] [PubMed]
  44. Chen, B.; Chen, H.; Lu, S.; Zhu, X.; Que, Y.; Zhang, Y.; Huang, J.; Zhang, L.; Zhang, Y.; Sun, F.; et al. KDM5B promotes tumorigenesis of Ewing sarcoma via FBXW7/CCNE1 axis. Cell Death Dis. 2022, 13, 354. [Google Scholar] [CrossRef] [PubMed]
  45. Yang, Y.; Li, W.; Wei, B.; Wu, K.; Liu, D.; Zhu, D.; Zhang, C.; Wen, F.; Fan, Y.; Zhao, S. MicroRNA let-7i Inhibits Histone Lysine Demethylase KDM5B to Halt Esophageal Cancer Progression. Mol. Ther. Nucleic Acids 2020, 22, 846–861. [Google Scholar] [CrossRef]
  46. Bamodu, O.A.; Huang, W.C.; Lee, W.H.; Wu, A.; Wang, L.S.; Hsiao, M.; Yeh, C.T.; Chao, T.Y. Aberrant KDM5B expression promotes aggressive breast cancer through MALAT1 overexpression and downregulation of hsa-miR-448. BMC Cancer 2016, 16, 160. [Google Scholar] [CrossRef] [Green Version]
  47. Ruan, Y.; Chen, W.; Gao, C.; Xu, Y.; Shi, M.; Zhou, Z.; Zhou, G. REXO4 acts as a biomarker and promotes hepatocellular carcinoma progression. J. Gastrointest. Oncol. 2021, 12, 3093–3106. [Google Scholar] [CrossRef] [PubMed]
  48. Chen, W.; Gao, C.; Shen, J.; Yao, L.; Liang, X.; Chen, Z. The expression and prognostic value of REXO4 in hepatocellular carcinoma. J. Gastrointest. Oncol. 2021, 12, 1704–1717. [Google Scholar] [CrossRef]
  49. Verma, S.K.; Deshmukh, V.; Thatcher, K.; Belanger, K.K.; Rhyner, A.M.; Meng, S.; Holcomb, R.J.; Bressan, M.; Martin, J.F.; Cooke, J.P.; et al. RBFOX2 is required for establishing RNA regulatory networks essential for heart development. Nucleic Acids Res. 2022, 50, 2270–2286. [Google Scholar] [CrossRef]
  50. Gehman, L.T.; Meera, P.; Stoilov, P.; Shiue, L.; O’Brien, J.E.; Meisler, M.H.; Ares, M., Jr.; Otis, T.S.; Black, D.L. The splicing regulator Rbfox2 is required for both cerebellar development and mature motor function. Genes Dev. 2012, 26, 445–460. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Newman, E.A.; Kim, D.W.; Wan, J.; Wang, J.; Qian, J.; Blackshaw, S. Foxd1 is required for terminal differentiation of anterior hypothalamic neuronal subtypes. Dev. Biol. 2018, 439, 102–111. [Google Scholar] [CrossRef] [PubMed]
  52. Hernandez-Bejarano, M.; Gestri, G.; Spawls, L.; Nieto-Lopez, F.; Picker, A.; Tada, M.; Brand, M.; Bovolenta, P.; Wilson, S.W.; Cavodeassi, F. Opposing Shh and Fgf signals initiate nasotemporal patterning of the zebrafish retina. Development 2015, 142, 3933–3942. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Braeutigam, C.; Rago, L.; Rolke, A.; Waldmeier, L.; Christofori, G.; Winter, J. The RNA-binding protein Rbfox2: An essential regulator of EMT-driven alternative splicing and a mediator of cellular invasion. Oncogene 2014, 33, 1082–1092. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Kim, J.H.; Jeong, K.; Li, J.; Murphy, J.M.; Vukadin, L.; Stone, J.K.; Richard, A.; Tran, J.; Gillespie, G.Y.; Flemington, E.K.; et al. SON drives oncogenic RNA splicing in glioblastoma by regulating PTBP1/PTBP2 switching and RBFOX2 activity. Nat. Commun. 2021, 12, 5551. [Google Scholar] [CrossRef] [PubMed]
  55. Naro, C.; De Musso, M.; Delle Monache, F.; Panzeri, V.; de la Grange, P.; Sette, C. The oncogenic kinase NEK2 regulates an RBFOX2-dependent pro-mesenchymal splicing program in triple-negative breast cancer cells. J. Exp. Clin. Cancer Res. 2021, 40, 397. [Google Scholar] [CrossRef]
  56. Cooper, A.M.; Nutter, C.A.; Kuyumcu-Martinez, M.N.; Wright, C.W. Alternative Splicing of the Aryl Hydrocarbon Receptor Nuclear Translocator (ARNT) Is Regulated by RBFOX2 in Lymphoid Malignancies. Mol. Cell Biol. 2022, 42, e0050321. [Google Scholar] [CrossRef]
  57. Cao, J.; O’Day, D.R.; Pliner, H.A.; Kingsley, P.D.; Deng, M.; Daza, R.M.; Zager, M.A.; Aldinger, K.A.; Blecher-Gonen, R.; Zhang, F.; et al. A human cell atlas of fetal gene expression. Science 2020, 370, eaba7721. [Google Scholar] [CrossRef] [PubMed]
  58. Gautam, P.; Hamashima, K.; Chen, Y.; Zeng, Y.; Makovoz, B.; Parikh, B.H.; Lee, H.Y.; Lau, K.A.; Su, X.; Wong, R.C.B.; et al. Multi-species single-cell transcriptomic analysis of ocular compartment regulons. Nat. Commun. 2021, 12, 5675. [Google Scholar] [CrossRef]
  59. Herreros-Villanueva, M.; Zhang, J.S.; Koenig, A.; Abel, E.V.; Smyrk, T.C.; Bamlet, W.R.; de Narvajas, A.A.; Gomez, T.S.; Simeone, D.M.; Bujanda, L.; et al. SOX2 promotes dedifferentiation and imparts stem cell-like features to pancreatic cancer cells. Oncogenesis 2013, 2, e61. [Google Scholar] [CrossRef] [PubMed]
  60. Ray, D.; Yun, Y.C.; Idris, M.; Cheng, S.; Boot, A.; Iain, T.B.H.; Rozen, S.G.; Tan, P.; Epstein, D.M. A tumor-associated splice-isoform of MAP2K7 drives dedifferentiation in MBNL1-low cancers via JNK activation. Proc. Natl. Acad. Sci. USA 2020, 117, 16391–16400. [Google Scholar] [CrossRef]
  61. Malta, T.M.; Sokolov, A.; Gentles, A.J.; Burzykowski, T.; Poisson, L.; Weinstein, J.N.; Kaminska, B.; Huelsken, J.; Omberg, L.; Gevaert, O.; et al. Machine Learning Identifies Stemness Features Associated with Oncogenic Dedifferentiation. Cell 2018, 173, 338–354.e15. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Field, M.G.; Kuznetsov, J.N.; Bussies, P.L.; Cai, L.Z.; Alawa, K.A.; Decatur, C.L.; Kurtenbach, S.; Harbour, J.W. BAP1 Loss Is Associated with DNA Methylomic Repatterning in Highly Aggressive Class 2 Uveal Melanomas. Clin. Cancer Res. 2019, 25, 5663–5673. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Li, D.; Fan, S.; Yu, F.; Zhu, X.; Song, Y.; Ye, M.; Fan, L.; Lv, Z. FOXD1 Promotes Cell Growth and Metastasis by Activation of Vimentin in NSCLC. Cell. Physiol. Biochem. 2018, 51, 2716–2731. [Google Scholar] [CrossRef] [PubMed]
  64. Sun, Q.; Novak, D.; Huser, L.; Poelchen, J.; Wu, H.; Granados, K.; Federico, A.; Liu, K.; Steinfass, T.; Vierthaler, M.; et al. FOXD1 promotes dedifferentiation and targeted therapy resistance in melanoma by regulating the expression of connective tissue growth factor. Int. J. Cancer 2021, 149, 657–674. [Google Scholar] [CrossRef]
  65. Chen, C.; Xu, Z.Q.; Zong, Y.P.; Ou, B.C.; Shen, X.H.; Feng, H.; Zheng, M.H.; Zhao, J.K.; Lu, A.G. CXCL5 induces tumor angiogenesis via enhancing the expression of FOXD1 mediated by the AKT/NF-kappaB pathway in colorectal cancer. Cell Death Dis. 2019, 10, 178. [Google Scholar] [CrossRef] [Green Version]
  66. Chang, S.; Sun, L.; Feng, G. SP1-mediated long noncoding RNA POU3F3 accelerates the cervical cancer through miR-127-5p/FOXD1. Biomed. Pharm. 2019, 117, 109133. [Google Scholar] [CrossRef]
  67. Nagel, S.; Meyer, C.; Kaufmann, M.; Drexler, H.G.; MacLeod, R.A. Deregulated FOX genes in Hodgkin lymphoma. Genes Chromosomes Cancer 2014, 53, 917–933. [Google Scholar] [CrossRef]
  68. Sun, D.S.; Guan, C.H.; Wang, W.N.; Hu, Z.T.; Zhao, Y.Q.; Jiang, X.M. LncRNA NORAD promotes proliferation, migration and angiogenesis of hepatocellular carcinoma cells through targeting miR-211-5p/FOXD1/VEGF-A axis. Microvasc. Res. 2021, 134, 104120. [Google Scholar] [CrossRef]
  69. Ma, X.L.; Shang, F.; Ni, W.; Zhu, J.; Luo, B.; Zhang, Y.Q. MicroRNA-338-5p plays a tumor suppressor role in glioma through inhibition of the MAPK-signaling pathway by binding to FOXD1. J. Cancer Res. Clin. Oncol. 2018, 144, 2351–2366. [Google Scholar] [CrossRef]
  70. Wu, Q.; Ma, J.; Wei, J.; Meng, W.; Wang, Y.; Shi, M. FOXD1-AS1 regulates FOXD1 translation and promotes gastric cancer progression and chemoresistance by activating the PI3K/AKT/mTOR pathway. Mol. Oncol. 2021, 15, 299–316. [Google Scholar] [CrossRef]
  71. Zhou, H.; Lv, Q.; Guo, Z. Transcriptomic signature predicts the distant relapse in patients with ER+ breast cancer treated with tamoxifen for five years. Mol. Med. Rep. 2018, 17, 3152–3157. [Google Scholar] [CrossRef] [Green Version]
  72. Chen, Y.; Liang, W.; Liu, K.; Shang, Z. FOXD1 promotes EMT and cell stemness of oral squamous cell carcinoma by transcriptional activation of SNAI2. Cell Biosci. 2021, 11, 154. [Google Scholar] [CrossRef] [PubMed]
  73. Serrano-Gomez, S.J.; Maziveyi, M.; Alahari, S.K. Regulation of epithelial-mesenchymal transition through epigenetic and post-translational modifications. Mol. Cancer 2016, 15, 18. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Decatur, C.L.; Ong, E.; Garg, N.; Anbunathan, H.; Bowcock, A.M.; Field, M.G.; Harbour, J.W. Driver Mutations in Uveal Melanoma: Associations With Gene Expression Profile and Patient Outcomes. JAMA Ophthalmol. 2016, 134, 728–733. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Uner, O.E.; See, T.R.O.; Szalai, E.; Grossniklaus, H.E.; Stalhammar, G. Estimation of the timing of BAP1 mutation in uveal melanoma progression. Sci. Rep. 2021, 11, 8923. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Schematic overview of the single-cell workflow from raw data files to validation of gene expression in uveal melanoma.
Figure 1. Schematic overview of the single-cell workflow from raw data files to validation of gene expression in uveal melanoma.
Cancers 14 03668 g001
Figure 2. Uniform manifold approximation and projection (UMAP) plot depicting single-cell transcriptomes from whole zebrafish larvae aged (a) 24, (b) 48 and (c) 120 h post fertilization (hpf). Each dot represents a single cell. Melanocyte cluster gene expression was isolated from pmelb-negative cells (blue) and pmelb-positive cells (red) and combined for further analysis (df).
Figure 2. Uniform manifold approximation and projection (UMAP) plot depicting single-cell transcriptomes from whole zebrafish larvae aged (a) 24, (b) 48 and (c) 120 h post fertilization (hpf). Each dot represents a single cell. Melanocyte cluster gene expression was isolated from pmelb-negative cells (blue) and pmelb-positive cells (red) and combined for further analysis (df).
Cancers 14 03668 g002
Figure 3. (a) Venn diagram of the melanocytic clusters in the zebrafish atlas. Clusters were divided into pmelb-positive and pmelb-negative melanocytes. A total of 1817 genes were expressed in both clusters, whereas 637 genes were uniquely expressed in the pmelb-positive cluster and 279 genes in the pmelb-negative cluster. (b) Total counts of melanocyte markers in the pmelb-positive and pmelb-negative clusters. An overall high total count of most melanocyte markers was assessed in both clusters, while pmelb is only expressed in one cluster.
Figure 3. (a) Venn diagram of the melanocytic clusters in the zebrafish atlas. Clusters were divided into pmelb-positive and pmelb-negative melanocytes. A total of 1817 genes were expressed in both clusters, whereas 637 genes were uniquely expressed in the pmelb-positive cluster and 279 genes in the pmelb-negative cluster. (b) Total counts of melanocyte markers in the pmelb-positive and pmelb-negative clusters. An overall high total count of most melanocyte markers was assessed in both clusters, while pmelb is only expressed in one cluster.
Cancers 14 03668 g003
Figure 4. (a) Gene function annotation from IPA, with including 49 transcription regulars. Expression of the 49 transcription regulator genes (expressed in early-developed zebrafish larvae) in (b) human healthy uveal melanocytes and (c) UM. Twenty-one genes were expressed in human healthy uveal melanocyte whereas the majority of the genes are expressed in UM. Red line = gene expression threshold of 100.
Figure 4. (a) Gene function annotation from IPA, with including 49 transcription regulars. Expression of the 49 transcription regulator genes (expressed in early-developed zebrafish larvae) in (b) human healthy uveal melanocytes and (c) UM. Twenty-one genes were expressed in human healthy uveal melanocyte whereas the majority of the genes are expressed in UM. Red line = gene expression threshold of 100.
Cancers 14 03668 g004
Figure 5. (a) Expression of ELL2 in UM, grouped per mutational subtype and (b) survival plot of expression in the combined (ROMS and TCGA) cohorts. (c) Expression of RBFOX2 in UM, grouped per mutational subtype and (d) survival plot of expression in the combined (ROMS and TCGA) cohorts. (e) Expression of FOXD1 in UM, grouped per mutational subtype and (f) survival plot of expression in the combined (ROMS and TCGA) cohorts and (g) BAP1-mutated UM within the combined cohort. ***, p < 0.001; ****, p < 0.0001.
Figure 5. (a) Expression of ELL2 in UM, grouped per mutational subtype and (b) survival plot of expression in the combined (ROMS and TCGA) cohorts. (c) Expression of RBFOX2 in UM, grouped per mutational subtype and (d) survival plot of expression in the combined (ROMS and TCGA) cohorts. (e) Expression of FOXD1 in UM, grouped per mutational subtype and (f) survival plot of expression in the combined (ROMS and TCGA) cohorts and (g) BAP1-mutated UM within the combined cohort. ***, p < 0.001; ****, p < 0.0001.
Cancers 14 03668 g005
Figure 6. Heatmap using Z−scores of the GEP test of UM, including FOXD1, RBFOX2, ELL2, KDM5B and REXO4 on UM samples. Samples were sorted based on disomy 3 or monosomy 3 status. FOXD1 and RBFOX2 clustered similarly to upregulated genes, whereas ELL2, KDM5B and REXO4 were more similarly clustered to downregulated genes.
Figure 6. Heatmap using Z−scores of the GEP test of UM, including FOXD1, RBFOX2, ELL2, KDM5B and REXO4 on UM samples. Samples were sorted based on disomy 3 or monosomy 3 status. FOXD1 and RBFOX2 clustered similarly to upregulated genes, whereas ELL2, KDM5B and REXO4 were more similarly clustered to downregulated genes.
Cancers 14 03668 g006
Figure 7. Immunohistochemical expression of FOXD1; (a) FOXD1-positive UM and (b) FOXD1-negative UM (scale 5 µm). (c) The survival plot of FOXD1 expression on IHC in the combined cohorts (discovery and replication cohorts). The expression of FOXD1 was correlated with a poorer survival (n = 59, p = 0.0439).
Figure 7. Immunohistochemical expression of FOXD1; (a) FOXD1-positive UM and (b) FOXD1-negative UM (scale 5 µm). (c) The survival plot of FOXD1 expression on IHC in the combined cohorts (discovery and replication cohorts). The expression of FOXD1 was correlated with a poorer survival (n = 59, p = 0.0439).
Cancers 14 03668 g007
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

van den Bosch, Q.C.C.; Nguyen, J.Q.N.; Brands, T.; van den Bosch, T.P.P.; Verdijk, R.M.; Paridaens, D.; Naus, N.C.; de Klein, A.; Kiliç, E.; Brosens, E., on behalf of the Rotterdam Ocular Melanoma Study Group . FOXD1 Is a Transcription Factor Important for Uveal Melanocyte Development and Associated with High-Risk Uveal Melanoma. Cancers 2022, 14, 3668. https://doi.org/10.3390/cancers14153668

AMA Style

van den Bosch QCC, Nguyen JQN, Brands T, van den Bosch TPP, Verdijk RM, Paridaens D, Naus NC, de Klein A, Kiliç E, Brosens E on behalf of the Rotterdam Ocular Melanoma Study Group . FOXD1 Is a Transcription Factor Important for Uveal Melanocyte Development and Associated with High-Risk Uveal Melanoma. Cancers. 2022; 14(15):3668. https://doi.org/10.3390/cancers14153668

Chicago/Turabian Style

van den Bosch, Quincy C. C., Josephine Q. N. Nguyen, Tom Brands, Thierry P. P. van den Bosch, Robert M. Verdijk, Dion Paridaens, Nicole C. Naus, Annelies de Klein, Emine Kiliç, and Erwin Brosens on behalf of the Rotterdam Ocular Melanoma Study Group . 2022. "FOXD1 Is a Transcription Factor Important for Uveal Melanocyte Development and Associated with High-Risk Uveal Melanoma" Cancers 14, no. 15: 3668. https://doi.org/10.3390/cancers14153668

APA Style

van den Bosch, Q. C. C., Nguyen, J. Q. N., Brands, T., van den Bosch, T. P. P., Verdijk, R. M., Paridaens, D., Naus, N. C., de Klein, A., Kiliç, E., & Brosens, E., on behalf of the Rotterdam Ocular Melanoma Study Group . (2022). FOXD1 Is a Transcription Factor Important for Uveal Melanocyte Development and Associated with High-Risk Uveal Melanoma. Cancers, 14(15), 3668. https://doi.org/10.3390/cancers14153668

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop