We have extensively studied the phenotypic heterogeneity of patient-derived melanoma cells. Here, whole-exome sequencing revealed novel variants of genes associated with the MAPK, NOTCH, Hippo, cell-cycle, senescence, and ubiquitindependent pathways, which could contribute to the observed phenotypic diversity between cell lines. Focusing on mutations in the MAPK pathway-associated genes, we found BRAF (BRAF V600E ) and RAS subtypes, including NRAS Q61R and the rare HRAS Q61R variant, and additional alterations potentially leading to different ERK1/2 activity. Both RAS Q61R cell lines harbored a MEK1 P124S variant and exerted a low level of phospho-MEK1/2. Activity of the MAPK pathway was further attenuated in NRAS Q61R /MEK P124S cells by trametinib, and this effect was also shown in HRAS Q61R / MEK P124S melanoma cells. The observed variability in doubling time might be a consequence of diverse MAPK and PI3K/AKT pathway activities, but not exclusively, as a senescence program was also executed to different extent in distinct melanoma cell lines. Low percentages of senescent cells might result from mutations in CDKN2A, E2F3, and EZH2, and a high c-MYC expression. Vemurafenib and trametinib induced senescence concomitantly with c-MYC downregulation and irrespectively of CDKN2A mutation, but the EZH2 S412C variant might limit senescence induction. Damaging alterations in Hippo pathway-associated genes were accompanied with variability in the phosphorylation level of YAP1/TAZ and CTGF expression. Our study also suggests opposite activity of NOTCH2 F1209V and NOTCH2 N2002S variants. Additionally, we found a novel FBXW7 V418M variant that retained its function in melanoma cells. The obtained molecular data might be further exploited in genotype-phenotype relationship studies and in identifying novel biomarkers and therapies for melanomas.