Identification of Sets of Cytoskeletal Related and Adhesion-related Coding Region Mutations in the TCGA Melanoma Dataset that Correlate with a Negative Outcome

The Driverless Triple-Wild-Type (BRAF, RAS, KIT) Cutaneous Melanoma: Whole Genome Sequencing Discoveries

The genetic makeup of the triple-wild-type melanoma (BRAF, NRAS and NF1) has been known for some time, but those studies grouped together rare histopathological versions with common ones, as well as mucosal and even uveal ones. Here we used whole genome sequencing to genetically characterize the triple-wild-type melanoma (TWM), termed here as BRAF, RAS and KIT wild type (the most frequent oncogenic drivers of skin melanoma), using the most common histological forms and excluding rare ones. All these tumors except one were clearly induced by UV based on the mutational signature. The tumor mutational burden was low in TWM, except in the NF1 mutant forms, and a relatively high frequency of elevated LOH scores suggested frequent homologue recombination deficiency, but this was only confirmed by the mutation signature in one case. Furthermore, all these TWMs were microsatellite-stabile. In this driverless setting, we revealed rare oncogenic drivers known from melanoma or other cancer types and identified rare actionable tyrosine kinase mutations in NTRK1, RET and VEGFR1. Mutations of TWM identified genes involved in antitumor immunity (negative and positive predictors of immunotherapy), Ca++ and BMP signaling. The two regressed melanomas of this cohort shared a 17-gene mutation signature, containing genes involved in antitumor immunity and several cell surface receptors. Even with this comprehensive genomic approach, a few cases remained driverless, suggesting that unrecognized drivers are hiding among passenger mutations.

Melanoma molecular subtyping and scoring model construction based on ligand-receptor pairs

Melanoma is a malignancy of melanocytes, responsible for a high percentage of skin cancer mortality. Ligand-Receptor pairs, a type of cellular communication, are essential for tumor genesis, growth, metastasis, and prognosis. Yet, the role of Ligand-Receptor pairs in melanoma has not been fully elucidated. Our research focused on the function of Ligand-Receptor pairs in melanoma prognosis. We screened 131 melanoma prognosis corresponded ligand-receptor pairs by analyzing the TCGA data of melanoma and the 2293 LR pairs retrieved from the connectomeDB2020 database. And further developed subtypes of melanoma according to the expression of these ligand-receptor pairs by Consensus Clustering. Then we using lasso cox regression and stepwise multivariate regression analysis established a ligand-receptor pairs-based scoring model for the evaluation of melanoma prognosis. Our study demonstrated that the ligand-receptor pairs are vital to the molecular heterogeneity of melanoma, and characterized three different melanoma ligand-receptor pairs subtypes. Among them, the C3 subtype showed a better prognosis, while the C1 subtype exhibited a low prognosis state. And our analysis then found out that this could be related to the differed activation and inhabitation of the cell cycle and immune-related pathways. Using lasso cox regression and stepwise multivariate regression analysis, we further identified 9 key ligand-receptor pairs and established a scoring model that effectively correlated with the prognosis, immune pathways, and therapy of melanoma, showing that the LR.score model was a trustworthy and independent biomarker for melanoma prognosis evaluation. In sum, we found that ligand-receptor pairs are significantly associated with the prognosis and therapy of melanoma. And our ligand-receptor-based scoring model showed potential for the evaluation of melanoma prognosis and immune therapy outcome prediction, which is crucial to the survival for the patients.

Somatic synonymous mutations in regulatory elements contribute to the genetic aetiology of melanoma

Background

Non-synonymous mutations altering tumor suppressor genes and oncogenes are widely studied. However, synonymous mutations, which do not alter the protein sequence, are rarely investigated in melanoma genome studies.

Methods

We explored the role of somatic synonymous mutations in melanoma samples from TCGA (The Cancer Genome Atlas). The pathogenic synonymous mutation and neutral synonymous mutation data were used to assess the significance of pathogenic synonymous mutations in melanoma likely to affect genetic regulatory elements using Fisher’s exact test. Poisson distribution probabilities of each gene were used to mine the genes with multiple potential functional synonymous mutations affecting regulatory elements.

Results

Concentrating on five types of genetic regulatory functions, we found that the mutational patterns of pathogenic synonymous mutations are mostly involved in exonic splicing regulators in near-splicing sites or inside DNase I hypersensitivity sites or non-optimal codon. Moreover, the sites of miRNA binding alteration exhibit a significantly lower rate of evolution than other sites. Finally, 12 genes were hit by recurrent potentially functional synonymous mutations, which showed statistical significance in the pathogenic mutations. Among them, nine genes (DNAH5, ADCY8, GRIN2A, KSR2, TECTA, RIMS2, XKR6, MYH1, SCN10A) have been reported to be mutated in melanoma, and other three genes (SLC9A2, CASR, SLC8A3) have a great potential to impact melanoma.

Conclusion

These findings confirm the functional consequences of somatic synonymous mutations in melanoma, emphasizing the significance of research in future studies.

Mutant cytoskeletal and ECM peptides sensitive to the ST14 protease are associated with a worse outcome for glioblastoma multiforme

We previously identified a set of the most frequently mutated cytoskeleton- and extracellular matrix-related proteins (CECMPs) in numerous cancer datasets. In this report, we used a bioinformatics approach to assess the impact of amino acid (AA) substitutions on the sensitivity of CECMPs to the ST14 protease (matriptase I), a transmembrane serine protease previously implicated in cancer development. Results indicated that AA substitutions in glioblastoma multiforme (GBM) CECMPs are skewed toward increased resistance to the ST14 protease, in comparison to the wild-type peptide sequence. Furthermore, the protease resistant AA substitutions represent relatively high binding affinities to HLA class I proteins, when assessing the binding specificities using HLA class I alleles matched to the source of the mutant AA. Moreover, samples representing AA substitutions that increased protease sensitivity also represented reduced overall and disease-free survival periods for patients with glioblastoma. To assess tumor specimen immunogenicity, we identified T-cell receptor (TCR) V(D)J recombinations in GBM exome files. The overlap between ST14 protease sensitive mutant barcodes and the TCR V(D)J recombination read positive barcodes represented significantly reduced survival.