Aberrant RNA Splicing Events Driven by Mutations of RNA-Binding Proteins as Indicators for Skin Cutaneous Melanoma Prognosis

Identification of Alternative Splicing in Proteomes of Human Melanoma Cell Lines without RNA Sequencing Data

Alternative splicing is one of the main regulation pathways in living cells beyond simple changes in the level of protein expression. Most of the approaches proposed in proteomics for the identification of specific splicing isoforms require a preliminary deep transcriptomic analysis of the sample under study, which is not always available, especially in the case of the re-analysis of previously acquired data. Herein, we developed new algorithms for the identification and validation of protein splice isoforms in proteomic data in the absence of RNA sequencing of the samples under study. The bioinformatic approaches were tested on the results of proteome analysis of human melanoma cell lines, obtained earlier by high-resolution liquid chromatography and mass spectrometry (LC-MS). A search for alternative splicing events for each of the cell lines studied was performed against the database generated from all known transcripts (RefSeq) and the one composed of peptide sequences, which included all biologically possible combinations of exons. The identifications were filtered using the prediction of both retention times and relative intensities of fragment ions in the corresponding mass spectra. The fragmentation mass spectra corresponding to the discovered alternative splicing events were additionally examined for artifacts. Selected splicing events were further validated at the mRNA level by quantitative PCR.

Expression profile and prognostic values of LSM family in skin cutaneous melanoma

Background

The like-Smith (LSM) family is a group of RNA-binding proteins involved in RNA metabolism. However, their involvement in tumors, particularly skin cutaneous melanoma (SKCM), is not fully understood. In this study, we focused on the expression profiles and prognostic values of the LSM family in SKCM.

Methods

Raw data were downloaded from The Cancer Genome Atlas. The expression profile and prognostic value of LSM genes in SKCM were explored using the GEPIA, cBioPortal, and HPA databases. Protein–protein and gene–gene interaction analyses were performed using STRING and GeneMANIA. Enrichment and Cox regression analysis were conducted using R software. The TISIDB database was used to explore the relationship between LSMs and immunomodulators. Receiver operating characteristic curves and nomogram models were constructed to validate prognostic values.

Results

mRNA and protein expression levels of LSM2, LSM4, and LSM12 were significantly elevated in SKCM. The upregulated mRNA expression of LSM2 (p = 0.0013) and LSM4 (p = 0.0043) was significantly correlated with poor overall survival in patients with SKCM, whereas only LSM2 (p = 0.049) overexpression was markedly associated with worse disease-free survival. LSM2 overexpression was an independent risk factor (p = 0.013) and was confirmed to have a high prognostic value in SKCM using the receiver operating characteristic curve (AUC = 0.942) and nomogram models. All LSM genes were identified as genomic mutations, whereas alteration of LSM2 (p = 0.0153) significantly affected the overall survival in patients with SKCM. Significant correlations were observed between LSM family expression, immune cell infiltration, and immunomodulator. Furthermore, function and pathway enrichment analysis showed that the LSM family was mainly RNA binding proteins and involved in RNA splicing and degradation.

Conclusion

Expression profiles and prognostic values of LSM in SKCM were inconsistent. Among the LSM family, only LSM2 may serve as a potential poor prognosticator and immunotherapeutic target of SKCM.

Lq-based robust analytics on ultrahigh and high dimensional data

Ultrahigh and high dimensional data are common in regression analysis for various fields, such as omics data, finance, and biological engineering. In addition to the problem of dimension, the data might also be contaminated. There are two main types of contamination: outliers and model misspecification. We develop an unique method that takes into account the ultrahigh or high dimensional issues and both types of contamination. In this article, we propose a framework for feature screening and selection based on the minimum Lq-likelihood estimation (MLqE), which accounts for the model misspecification contamination issue and has also been shown to be robust to outliers. In numerical analysis, we explore the robustness of this framework under different outliers and model misspecification scenarios. To examine the performance of this framework, we conduct real data analysis using the skin cutaneous melanoma data. When comparing with traditional screening and feature selection methods, the proposed method shows superiority in both variable identification effectiveness and parameter estimation accuracy.

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.