Tumor cell invasion is one of the key processes during cancer progression, leading to life-threatening metastatic lesions in melanoma. As methylation of cancer-related genes plays a fundamental role during tumorigenesis and may lead to cellular plasticity which promotes invasion, our aim was to identify novel epigenetic markers on selected invasive melanoma cells. Using Illumina BeadChip assays and Affymetrix Human Gene 1.0 microarrays, we explored the DNA methylation landscape of selected invasive melanoma cells and examined the impact of DNA methylation on gene expression patterns. Our data revealed predominantly hypermethylated genes in the invasive cells affecting the neural crest differentiation pathway and regulation of the actin cytoskeleton. Integrative analysis of the methylation and gene expression profiles resulted in a cohort of hypermethylated genes (IL12RB2, LYPD6B, CHL1, SLC9A3, BAALC, FAM213A, SORCS1, GPR158, FBN1 and ADORA2B) with decreased expression. On the other hand, hypermethylation in the gene body of the EGFR and RBP4 genes was positively correlated with overexpression of the genes. We identified several methylation changes that can have role during melanoma progression, including hypermethylation of the promoter regions of the ARHGAP22 and NAV2 genes that are commonly altered in locally invasive primary melanomas as well as during metastasis. Interestingly, the down-regulation of the methylcytosine dioxygenase TET2 gene, which regulates DNA methylation, was associated with hypermethylated promoter region of the gene. This can probably lead to the observed global hypermethylation pattern of invasive cells and might be one of the key changes during the development of malignant melanoma cells. K E Y W O R D S cell invasion, DNA methylation, gene body, malignant melanoma, TET2

Selective inhibition of the mutant BRAF protein is a highly promising therapeutic approach for melanoma patients carrying the BRAF V600E mutation. Despite the remarkable clinical response, most patients develop resistance and experience tumour regrowth. To clarify the molecular background of BRAF inhibitor resistance, we generated four drug-resistant melanoma cell lines from paired primary/metastatic cell lines using a vemurafenib analogue PLX4720. Three of the resistant cell lines showed decreased proliferation after drug withdrawal, but the proliferation of one cell line (WM278RES) increased notably. Furthermore, we observed opposite phenomena in which a ‘drug holiday’ could not only be beneficial but also contribute to tumour progression. Using genomic and proteomic approaches, we found significantly different alterations between the sensitive and resistant cell lines, some of which have not been reported previously. In addition to several other changes, copy number gains were observed in all resistant cell lines on 8q24.11–q24.12 and 8q21.2. Gene expression analysis showed that most genes upregulated in the resistant cell lines were associated with cell motility and angiogenesis. Increased expression of six proteins (ANGPLT4, EGFR, Endoglin, FGF2, SerpinE1 and VCAM-1) and decreased expression of two proteins (osteopontin and survivin) were observed consistently in all resistant cell lines. In summary, we identified new genomic alterations and characterized the protein expression patterns associated with the resistant phenotype. Although several proteins have been shown to be associated with BRAF resistance, our study is the first to describe the association of VCAM-1 and osteopontin with BRAF resistance.

A large variety of molecular pathways in melanoma progression suggests that no individual molecular alteration is crucial in itself. Our aim was to define the molecular alterations underlying metastasis formation. Gene expression profiling was performed using microarray and qRT-PCR to define alterations between matched primary and metastatic melanoma cell lines. These data were integrated with publicly available unmatched tissue data. The invasiveness of cell lines was determined by Matrigel invasion assays and invasive clones from primary melanoma-derived cell lines were also selected. Two metastatic cell line models were created: the regional lymph node WM983A-WM983A-WM983B and the distant lung WM793B-WM793B-1205Lu metastatic models. The majority of metastasis genes were downregulated and enriched in adhesion and ITGA6-B4 pathways. Upregulation of immune pathways was characteristic of distant metastases, whereas increased Rap1 signaling was specific for regional (sub)cutaneous metastases. qRT-PCR analysis of selected integrins (A2, A3, A4, A9, B5, B8, A6, B1, and B3) highlighted the possible importance of ITGA3/4 and B8 in the metastatic process, distinguishing regional and distant metastases. We identified functionally relevant gene clusters that influenced metastasis formation. Our data provide further evidence that integrin expression patterns may be important in distant metastasis formation.

Glyphosate is the most commonly used herbicide around the world, which led to its accumulation in the environment and consequent ubiquitous human exposure. Glyphosate is marketed in numerous glyphosate-based herbicide formulations (GBHs) that include co-formulants to enhance herbicidal effect of the active ingredient, but are declared as inert substances. However, these other ingredients can have biologic activity on their own and may interact with the glyphosate in synergistic toxicity. In this study, we focused to compare the cytogenetic effect of the active ingredient glyphosate and three marketed GBHs (Roundup Mega, Fozat 480, and Glyfos) by investigating cytotoxicity with fluorescent co-labeling and WST-1 cell viability assay as well as genotoxicity with cytokinesis block micronucleus assay in isolated human mononuclear white blood cells. Glyphosate had no notable cytotoxic activity over the tested concentration range (0–10,000 μM), whereas all the selected GBHs induced significant cell death from 1,000 μM regardless of metabolic activation (S9). Micronucleus (MN) formation induced by glyphosate and its formulations at sub-cytotoxic concentrations (0–100 μM) exhibited a diverse pattern. Glyphosate caused statistically significant increase of MN frequency at the highest concentration (100 μM) after 20-h exposure. Contrarily, Roundup Mega exerted a significant genotoxic effect at 100 μM both after 4- and 20-h exposures; moreover, Glyfos and Fozat 480 also resulted in a statistically significant increase of MN frequency from the concentration of 10 μM after 4-h and 20-h treatment, respectively. The presence of S9 had no effect on MN formation induced by either glyphosate or GBHs. The differences observed in the cytotoxic and genotoxic pattern between the active principle and formulations confirm the previous concept that the presence of co-formulants in the formulations or the interaction of them with the active ingredient is responsible for the increased toxicity of herbicide products, and draw attention to the fact that GBHs are still currently in use, the toxicity of which rivals that of POEA-containing formulations (e.g., Glyfos) already banned in Europe. Hence, it is advisable to subject them to further comprehensive toxicological screening to assess the true health risks of exposed individuals, and to reconsider their free availability to any users.

Osteopontin (OPN) is a multifunctional glycoprotein that physiologically interacts with different types of integrins. It is considered to be a possible prognostic biomarker in certain tumor types; however, various splicing isoforms exist, which have not been investigated in melanoma. We aimed to define the relative expression pattern of five OPN isoforms and clarify the prognostic significance of the splice variants in melanoma. We also aimed to investigate the expression pattern of eight integrins in the same tumors. Gene expression analyses revealed that the relative expression of OPNa, OPNb, and OPNc is significantly higher in metastatic tumors compared to primary lesions (p < 0.01), whereas the expression of OPN4 and OPN5 was low in both. The more aggressive nodular melanomas had higher expression levels compared to the superficial spreading subtype (p ≤ 0.05). The relative expression of the eight tested integrins was low, with only the expression of ITGB3 being detectable in nodular melanoma (Medianlog2 = 1.274). A positive correlation was found between Breslow thickness and the expression of OPNc variant, whereby thicker tumors (>4 mm) had significantly higher expression (p ≤ 0.05). The Breslow thickness was negatively correlated with the expression of OPN4, and similarly with ITGA2. OPNc also exhibited significant positive correlation with the presence of metastasis. Our data show that high expression of OPNa, OPNb, and especially OPNc and low expression of OPN4 and ITGA2 are associated with an advanced stage of tumor progression and poor prognosis in melanoma.