Profiling of prognostic alternative splicing in melanoma

Generation of Schwann cell derived melanocytes from hPSCs identifies pro-metastatic factors in melanoma

The neural crest (NC) is highly multipotent and generates diverse lineages in the developing embryo. However, spatiotemporally distinct NC populations display differences in fate potential, such as increased gliogenic and parasympathetic potential from later migrating, nerve-associated Schwann cell precursors (SCPs). Interestingly, while melanogenic potential is shared by both early migrating NC and SCPs, differences in melanocyte identity resulting from differentiation through these temporally distinct progenitors have not been determined. Here, we leverage a human pluripotent stem cell (hPSC) model of NC temporal patterning to comprehensively characterize human NC heterogeneity, fate bias, and lineage development. We captured the transition of NC differentiation between temporally and transcriptionally distinct melanogenic progenitors and identified modules of candidate transcription factor and signaling activity associated with this transition. For the first time, we established a protocol for the directed differentiation of melanocytes from hPSCs through a SCP intermediate, termed trajectory 2 (T2) melanocytes. Leveraging an existing protocol for differentiating early NC-derived melanocytes, termed trajectory 1 (T1), we performed the first comprehensive comparison of transcriptional and functional differences between these distinct melanocyte populations, revealing differences in pigmentation and unique expression of transcription factors, ligands, receptors and surface markers. We found a significant link between the T2 melanocyte transcriptional signature and decreased survival in melanoma patients in the cancer genome atlas (TCGA). We performed an in vivo CRISPRi screen of T1 and T2 melanocyte signature genes in a human melanoma cell line and discovered several T2-specific markers that promote lung metastasis in mice. We further demonstrated that one of these factors, SNRPB, regulates the splicing of transcripts involved in metastasis relevant functions such as migration, cell adhesion and proliferation. Overall, this study identifies distinct developmental trajectories as a source of diversity in melanocytes and implicates the unique molecular signature of SCP-derived melanocytes in metastatic melanoma.

The Identification of Prognostic and Metastatic Alternative Splicing in Skin Cutaneous Melanoma

Skin cutaneous melanoma (SKCM) is a type of highly invasive cancer originated from melanocytes. It is reported that aberrant alternative splicing (AS) plays an important role in the neoplasia and metastasis of many types of cancer. Therefore, we investigated whether ASEs of pre-RNA have such an influence on the prognosis of SKCM and the related mechanism of ASEs in SKCM. The RNA-seq data and ASEs data for SKCM patients were obtained from the TCGA and TCGASpliceSeq database. The univariate Cox regression revealed 1265 overall survival-related splicing events (OS-SEs). Screened by Lasso regression, 4 OS-SEs were identified and used to construct an effective prediction model (AUC: .904), whose risk score was proved to be an independent prognostic factor. Furthermore, Kruskal-Wallis test and Mann-Whitney-Wilcoxon test showed that an aberrant splicing type of aminoacyl tRNA synthetase complex-interacting multifunctional protein 2 (AIMP2) regulated by CDC-like kinase 1 (CLK1) was associated with the metastasis and stage of SKCM. Besides, the overlapped signal pathway for AIMP2 was galactose metabolism identified by the co-expression analysis. External database validation also confirmed that AIMP2, CLK1, and the galactose metabolism were associated with the metastasis and stage of SKCM patients. ChIP-seq and ATAC-seq methods further confirmed the transcription regulation of CLK1, AIMP2, and other key genes, whose cellular expression was detected by Single Cell Sequencing. In conclusion, we proposed that CLK1-regulated AIMP2-78704-ES might play a critical role in the tumorigenesis and metastasis of SKCM via galactose metabolism. Besides, we established an effective model with MTMR14-63114-ES, URI1-48867-ES, BATF2-16724-AP, and MED22-88025-AP to predict the metastasis and prognosis of SKCM patients.

Down-regulated expression of miR-582 predicts poor prognosis and facilitates melanoma progression by targeting FOXC1

Melanoma is one of the most common malignant tumors that originate from nerve sheath melanocytes and are produced in the skin and other organs. This paper mainly studied the role of miR-582 in melanoma. The expression of miR-582 in melanoma cell lines and tissues was evaluated by real-time quantitative PCR. Kaplan-Meier curve and Cox proportional hazards model analysis were used to investigate the prognostic value of miR-582 in melanoma. Cell proliferation, invasion, and migration were analyzed using CCK-8 and transwell assays. The results showed that miR-582 was significantly decreased in melanoma cells and tissues. Downregulation of miR-582 was associated with ulceration, lymph node metastasis, TNM stage, and poor overall survival. The functional results showed that low expression of miR-582 can promote cell proliferation, migration, and invasion. FOXC1 was a direct target of miR-582. Overall, the expression of miR-582 is downregulated in melanoma tissues and cell lines. Low expression of miR-582 is associated with prognosis and progression of melanoma by targeting FOXC1. miR-582 may be a prognostic biomarker and a new therapeutic strategy for melanoma.

A Transcriptome-Wide Isoform Landscape of Melanocytic Nevi and Primary Melanomas Identifies Gene Isoforms Associated with Malignancy

Genetic splice variants have become of central interest in recent years, as they play an important role in different cancers. Little is known about splice variants in melanoma. Here, we analyzed a genome-wide transcriptomic dataset of benign melanocytic nevi and primary melanomas (n = 80) for the expression of specific splice variants. Using kallisto, a map for differentially expressed splice variants in melanoma vs. benign melanocytic nevi was generated. Among the top genes with differentially expressed splice variants were Ras-related in brain 6B (RAB6B), a member of the RAS family of GTPases, Macrophage Scavenger Receptor 1 (MSR1), Collagen Type XI Alpha 2 Chain (COLL11A2), and LY6/PLAUR Domain Containing 1 (LYPD1). The Gene Ontology terms of differentially expressed splice variants showed no enrichment for functional gene sets of melanoma vs. nevus lesions, but between type 1 (pigmentation type) and type 2 (immune response type) melanocytic lesions. A number of genes such as Checkpoint Kinase 1 (CHEK1) showed an association of mutational patterns and occurrence of splice variants in melanoma. Moreover, mutations in genes of the splicing machinery were common in both benign nevi and melanomas, suggesting a common mechanism starting early in melanoma development. Mutations in some of these genes of the splicing machinery, such as Serine and Arginine Rich Splicing Factor A3 and B3 (SF3A3, SF3B3), were significantly enriched in melanomas as compared to benign nevi. Taken together, a map of splice variants in melanoma is presented that shows a multitude of differentially expressed splice genes between benign nevi and primary melanomas. The underlying mechanisms may involve mutations in genes of the splicing machinery.

Relevance of 2'-O-Methylation and Pseudouridylation for the Malignant Melanoma

Simple Summary

This study investigates the expression, the histological localization, and the influence of the factors involved in 2′-O-methylation and pseudouridylation on prognostic relevant markers, proliferation markers, overall survival, molecular immune surveillance and evasion mechanisms within the malignant melanoma. Statistically significant positive correlations to the expression of markers involved in cell proliferation were observed. The upregulation of the RNA modifying factors was of prognostic relevance in this tumor disease with a negative impact on the overall survival of melanoma patients. Furthermore, the factors involved in 2′-O-methylation and pseudouridylation were statistically significant negative correlated to the expression of human leukocyte antigen class I genes as well as of components of the antigen processing machinery.

Abstract

The two RNA modifications 2′-O-methylation and pseudouridylation occur on several RNA species including ribosomal RNAs leading to an increased translation as well as cell proliferation associated with distinct functions. Using malignant melanoma (MM) as a model system the proteins mediating these RNA modifications were for the first time analyzed by different bioinformatics tools and public available databases regarding their expression and histological localization. Next to this, the impact of these RNA-modifying factors on prognostic relevant processes and marker genes of malignant melanoma was investigated and correlated to immune surveillance and evasion strategies. The RNA modifying factors exerted statistically significant positive correlations to the expression of genes involved in cell proliferation and were statistically significant negative correlated to the expression of human leukocyte antigen class I genes as well as of components of the antigen processing machinery in malignant melanoma. Upregulation of the RNA modifying proteins was of prognostic relevance in this tumor disease with a negative impact on the overall survival of melanoma patients. Furthermore, the expression of known oncogenic miRs, which are induced in malignant melanoma, directly correlated to the expression of factors involved in these two RNA modifications.

Prognostic alternative splicing signature reveals the landscape of immune infiltration in Pancreatic Cancer

Background: Pancreatic cancer (PC) is an aggressive cancer with worse survival in the world. Emerging evidence suggested that the imbalance of alternative splicing (AS) is a hallmark of cancer and indicated poor prognosis of patients. Genes-derived splicing events can produce neoepitopes for immunotherapy. However, the profound study of splicing profiling in PC is still elusive. We aimed to identification of novel prognostic signature across a comprehensive splicing landscape and reveal their relationship with tumor-infiltrating immune cells in pancreatic cancer microenvironment.

Methods: Based on integrated analysis of splicing profiling and clinical data, differentially splicing events were filtered out. Then, stepwise Cox regression analysis was applied to identify survival-related splicing events and construct prognostic signature. Functional enrichment analysis was performed to explore biology function. Kaplan-Meier curves and receiver operating characteristic (ROC) curves were performed to validate the predictive effect of predictive signature. We also verified the clinical value of prognostic signature under the influence of different clinical parameters. For deeper analysis, we evaluated the correlation between prognostic signature and infiltrating immune cells by CIBERSORT.

Results: According to systematic analyzing, a final six splicing events were identified and validated the good prognostic capability in entire TCGA dataset, validation set 1 and validation set 2 by Kaplan-Meier curves (P < 0.0001). The area under the curve (AUC) of ROC curves were also confirmed the high predictive efficiency of the prognostic signature in these three cohorts (AUC = 0.857, 0.895 and 0.788). In order to validate whether prognostic signature highlights a correlation between AS and immune contexture, CIBERSORT was performed to analyze the proportion of tumor-infiltrating immune cells in PC. Based on prognostic signature, we identified survival-related immune cells including CD8 T cells (P = 0.0111), activated CD4 memory T cells (P = 0.0329) and resting mast cells (P = 0.0352).

Conclusion: In conclusion, our study contribute to provide a promising prognostic signature based on six splicing events and revealed prognosis-related immune cells which indeed represented novel tumor drivers and provide potential targets for personalized therapeutic.