1. Introduction
Skin cancer is among the most common cancers worldwide, with keratinocyte skin cancers, basal cell carcinoma, squamous cell carcinoma, and malignant melanoma being the primary types [
1]. Among these, malignant melanoma is highly aggressive and deadly [
2]. Melanomas exhibit high mutation rates, the BRAF oncogene being frequently mutated and contributing to rapid tumour growth and treatment resistance [
2,
3]. The newly diagnosed cases of melanoma have been rising steadily worldwide, with projections suggesting an increase of more than 50% between 2020 and 2040 [
4]. Although treatments such as the BRAF inhibitor Vemurafenib and immune checkpoint blockers bring hope to melanoma patients, primary lack of response and acquired treatment resistance remain urgent challenges [
5,
6,
7].
Conjugated linolenic acids (CLnAs) are a group of isomers of α-linolenic acid (ALA) characterized by at least two conjugated double bonds. These fatty acids are predominantly found in the seeds of a limited number of plants [
8]. For example, punicic acid (PunA, C18:3 c9, t11, c13), α-eleostearic acid (α-ESA, C18:3 c9, t11, t13) and jacaric acid (JA, C18:3 c8, t10, c12) account, respectively, for up to 80% of pomegranate seed oil, approximately 50% of
Ricinodendron heudelotii seed oil [
9] and approximately 35% of the blue jacaranda seed oil [
10]. CLnAs have shown strong anti-carcinogenic potential in multiple experimental models [
11]. An early study reported that tung oil, which has a high content of α-ESA (approximately 80% of total fatty acids), exerts intense cytotoxic effects on DLD-1 colorectal adenocarcinoma cells, HepG2 hepatoma cells and A549 lung carcinoma cells [
12]. Subsequent research confirmed the strong cytotoxicity of α-ESA in several colon cancer cell lines [
13]. Similarly, JA significantly inhibits leukaemia cell proliferation [
14] and was recently shown to suppress breast cancer cell proliferation [
15]. PunA also exhibits cytotoxicity towards HCT-116 colorectal and FaDu hypopharyngeal cancer cell lines [
16], as well as various prostate cancer cell lines [
17]. Interestingly, PunA, α-ESA, and JA are highly toxic to proliferative Caco-2 colorectal cancer cells, while the toxicity is significantly decreased when Caco-2 cells have differentiated into epithelial cells forming a functional intestinal barrier [
18]. Furthermore, α-ESA does not impair the growth of normal human liver cells, while effectively inhibiting the development of breast cancer cells [
19]. Collectively, these findings suggest that CLnA toxicity is preferentially directed toward cancer cells without affecting normal cells. Despite extensive evidence of CLnA cytotoxicity in various cancers, their effects on melanoma cells remain largely unexplored.
CLnAs were initially claimed to induce apoptosis in cancer cells [
19,
20]. However, a growing body of research has indicated that CLnAs induce cancer cell death through ferroptosis. Ferroptosis is an iron-dependent form of non-apoptotic regulated cell death caused by the unrestrained accumulation of lipid hydroperoxides in cell membranes, leading to lethal membrane damage [
21]. An early study already revealed that α-ESA is cytotoxic to human monocytic leukaemia cells through a mechanism involving lipid peroxidation [
22]. More recently, α-ESA was shown to induce
bona fide ferroptosis in triple negative breast cancer cells [
23], β-ESA in fibrosarcoma HT-1080, brain neuroblastoma SK-N-SH, and clear-cell renal carcinoma 786-O cells [
24], PunA in hypopharyngeal, colorectal and prostate cancer cells [
16,
17], and JA in both triple-negative and luminal A breast cancer cell lines. Yet, whether this ferroptotic mechanism extends to melanoma cells remains to be elucidated.
Cancer cells exhibit an elevated need for fatty acids, which are essential for energy production and the synthesis of new cell membranes to sustain cell proliferation [
25,
26]. The introduction of high amounts of peroxidable polyunsaturated fatty acids (PUFAs) is expected to enhance lipid hydroperoxide production in cancer cells, thereby predisposing them to ferroptosis [
27]. Acyl-CoA synthetase long-chain family member 4 (ACSL4) catalyses the activation of PUFAs into PUFA-CoAs, and functions as a ferroptosis inducer [
28,
29]. PUFA-CoAs can be incorporated into phospholipids (PLs) within lipid membrane structures, through the activity of multiple lysophosphatidylcholine acyltransferases (LPCATs) [
30]. In contrast, antioxidant enzymes that function at the lipid membrane act as ferroptosis inhibitors. This is particularly the case for glutathione peroxidase 4 (GPX4), a monomeric selenoenzyme that reduces toxic phospholipid hydroperoxides (PLOOH) to non-toxic phospholipid alcohols (PLOH), using the reduced form of glutathione as a co-substrate. Through its action, GPX4 prevents the accumulation of hydroperoxides in the membrane [
27,
31]. ACSL4 and GPX4 have thus emerged as valuable targets in cancer therapy [
32,
33]. Accordingly, these enzymes were selected for mechanistic studies aiming at evaluating the involvement of the ferroptosis process in the cytotoxic effects of CLnAs [
34].
This study employed three melanoma cell lines, two of human origin (A375, WM266.4) and one of zebrafish origin (ZMEL1). Including both human and zebrafish melanoma cells helped to overcome the limitations of single-species model and enhanced the translational relevance of our findings. The A375 melanoma cell line, derived from a primary tumour, known for its high angiogenic and metastatic potential, exhibits rapid tumour growth and invasive behaviour [
35], making it a well-established model for studying melanoma progression and therapeutic responses. Similarly, the WM266.4 cell line, originating from a metastatic lymph node, is highly metastatic and frequently used to investigate aggressive melanoma behaviour and effectiveness of anti-cancer treatments. In parallel, the zebrafish melanoma cell line ZMEL1 offers a distinct advantage: its gene expression profile closely mirrors that of human melanoma cell lines while providing an opportunity to explore ferroptosis in a non-mammalian context [
2,
36,
37]. Beyond this, the use of ZMEL1 is an essential step toward using zebrafish as an
in vivo model for melanoma research.
We first evaluated the toxicity of CLnAs on these melanoma cell lines. To elucidate the underlying cell death mechanism, we examined the impact of combining CLnAs with ferroptosis, necroptosis or apoptosis inhibitors. In addition, we used chemical inhibitors of ACSL4 and GPX4 to evaluate their roles in the CLnA-induced cell death process. Given the lower sensitivity of ZMEL1 cells to CLnAs, this zebrafish-derived line was further used to investigate CLnA incorporation and metabolic processing in melanoma cells. ZMEL1 cells were also used to assess the expression levels of the acsl4 and gpx4 genes in the presence of CLnAs. To our knowledge, this is the first study exploring the anti-melanoma potential of CLnAs, providing a new direction for melanoma treatment.
2. Materials and Methods
2.1 Cell culture
The human melanoma cell lines A375 and WM266.4 were kindly gifted by Professor Bénédicte Jordan (LDRI, UCLouvain, Belgium). The zebrafish melanoma cell line ZMEL1 was kindly provided by the Memorial Sloan Kettering Cancer Center (USA). Human melanoma cells were cultured at 37°C with 5% CO2, in RPMI-1640 medium (21875034, Gibco), supplemented with 10% foetal bovine serum (F7524-500, Merck) and 5% penicillin-streptomycin (15140122, Gibco). ZMEL1 cells were maintained at 28.5°C with 5% CO2, in Dulbecco’s Modified Eagle’s Medium (DMEM) (L0102-500, VWR), supplemented with 10% foetal bovine serum, 5% penicillin-streptomycin and 5% Glutamax (35050-038, Gibco). Cells were regularly monitored for mycoplasma contamination.
2.2 Cell viability test
Impact of fatty acids on melanoma cell viability. Before being tested on cells, all fatty acids were conjugated to bovine serum albumin (BSA, A7030-100G, Sigma) in phosphate-buffered saline (PBS, P4417-100TAB, Sigma) to achieve a FA:BSA ratio of 4:1 (w/w). The resulting stock solutions were then diluted in culture medium to different working concentrations. Melanoma cells were harvested from the flasks using 0.25% (w/v) Trypsin (15090046, Gibco). Human melanoma cells and zebrafish melanoma cells were then seeded into 96-well plates at an initial density of 1×10
4 and 4×10
4 cells per well, respectively. Following a 24-hour adhesion period, cells were cultured with medium containing different CLnA isomers, namely JA, PunA, α-ESA, and β-eleostearic acid (β-ESA, C18:3 t9, t11, t13), at different concentrations (5, 10, 20, 40 and 80 µM) for 72 h. Control treatments included either oleic acid (OLA), a monounsaturated fatty acid, or ALA, the non-conjugated counterpart of CLnAs, at different concentrations (5, 10, 20, 40 and 80 µM), as well as a negative control without any additional fatty acid. The treatment concentration and time were determined according to preliminary experiments and previous research [
16]. After 72 h of fatty acid treatment, the culture medium was removed, and 100 µL of PrestoBlue (12083745, Fisher Scientific) solution, diluted in PBS at a ratio of 1:9 (v/v), was added to each well. After 1 h of incubation at the cell culture temperature, cell viability was assessed using a Fluoroskan Ascent FL fluorometer (Thermo Scientific) at 530/584 nm (excitation/emission) according to the manufacturer’s instructions.
Assessment of cell death mechanisms underlying CLnA toxicity in melanoma cells. The three cell lines were seeded in 96-well plates. Human melanoma cells (A375 and WM266.4) were exposed to 5 µM PunA or JA, while ZMEL1 cells were treated with 20 µM PunA or JA, each in combination with increasing doses of ferroptosis inhibitors, either ferrostatin-1 (fer-1, SML0583, Sigma), α-tocopherol (α-T, 258024, Sigma), or deferoxamine mesylate (DFOM, D9533, Sigma), as well as with the necroptosis inhibitor necrostatin-1(nec-1, 480065, Sigma) and the apoptosis inhibitor ZVAD-FMK (ZVAD, S7023, Selleck Chemicals). Absolute ethanol was used as vehicle for α-T, while DMSO was used as vehicle for fer-1, DFOM, nec-1 and ZVAD. Additional 96-well plates were subjected to either 2.5 µM PunA (A375 and WM266.4 cells) or 10 µM PunA (ZMEL1 cells), in combination with an increasing dose of the ACSL4 inhibitor PRGL493 (HY-139180, MedChemExpress). Similarly, 0.625 µM PunA was used for the human melanoma cells and 5 µM PunA for ZMEL1 cells, to study the effects of the combination with increasing doses of the GPX4 inhibitors RSL3 (S8155, Selleck Chemicals) or ML210 (S0788, Selleck Chemicals). DMSO was used as vehicle for PRGL493, RSL3 and ML210.
2.3 Kinetics of ZMEL1 cell viability loss
To determine the optimal treatment concentration and time point to harvest ZMEL1 cells for fatty acid uptake studies and mRNA expression analyses, we assessed cell viability over time under different treatment conditions. ZMEL1 cells were seeded at an initial density of 4×104 cells per well into 96-well plates and were exposed to three concentrations (i.e. 10, 20 or 40 µM) of PunA or JA for either 0, 12, 24, 48 or 72 h followed by a viability test.
2.4 Fatty acid uptake and metabolic processing in ZMEL1 cells
Owing to their lower sensitivity to CLnAs in comparison to the human melanoma cells under investigation, ZMEL1 cells were used to evaluate the uptake and distribution of PunA, JA, or ALA in different fractions, including neutral lipids (NLs), free fatty acids (FFAs), and phospholipids (PLs). Cells were either untreated (control), or treated with 20 μM of PunA, JA or ALA, and were harvested at 0, 4 and 12 hours after treatment. Total lipids were extracted using the Bligh and Dyer method [
16,
38]. The internal standard, consisting of triheptadecanoin (for NLs), nonadecanoic acid (for FFAs), and 1,2-dibehenoyl-sn-glycero-3-phosphocholine (for PLs), was added to each sample to evaluate extraction efficiency. Samples were dried under nitrogen stream at 30°C and resuspended in chloroform. Resuspended samples were loaded onto solid phase extraction columns (12102089, Agilent Technologies) and NL, FFA and PL fractions were eluted with chloroform:2-propanol (2:1, v/v), diethyl ether:acetic acid (98:2, v/v) and methanol, respectively [
39,
40]. The eluted fractions were then evaporated under nitrogen stream at 30°C and methylated under alkaline conditions with 0.5 mL of 0.1 M KOH in methanol at 70°C for 1 hour, followed by acidic conditions at 70°C for 15 minutes after addition of 0.2 mL of 1.2 M HCl in methanol. Fatty acid methyl esters (FAMEs) were extracted using 1 mL of hexane. An injection standard, methyl-undecanoate (20-1100-13, Larodan), was added to each sample for accurate measurement during analysis. FAMEs were then injected and separated using a gas chromatograph (Trace 1310, Thermo Fisher Scientific) equipped with an autosampler (TriPlusAS) and a RT-2560 capillary column (biscyanopropylpolysiloxane, 100 m length, 0.25 mm internal diameter, 0.2 mm film thickness, Restek). The flow of hydrogen was used as the carrier gas at a constant pressure of 200 kPa. The temperature program for the GC was as follows: the temperature was initially set at 80°C; it was ramped up to 175°C at a rate of 25°C/min and was maintained for 25 minutes; it was then increased to 200°C at 10°C/min and was held for 20 minutes; it was then increased to 220°C at 10°C/min and was held for 5 minutes; it was then elevated to 235°C at 10°C/min with a final hold of 15 minutes, before being brought back to 80°C at a rate of 20°C/min. FAMEs were detected using a flame ionization detector (FID) set at a constant temperature of 255°C with an air flow of 350 mL/min, hydrogen flow of 35 mL/min, and nitrogen flow of 40 mL/min. A calibration standard consisting of a mixture of 43 pure methyl ester standards (Larodan and Nu-Check Prep) was utilized to identify unknown peaks based on their retention times and to quantify them using known concentrations. A PunA methyl ester standard (20-1875, Larodan) of known concentration was used to identify and quantify the PunA peak in each sample. JA was quantified on the basis of the PunA standard; the values were thus expressed in PunA equivalents. Chromatograms were processed using ChromQuest 5.0 software (Thermo Fisher Scientific).
2.5 Impact of CLnAs on acsl4 and gpx4 expression in ZMEL1 cells
For the assessment of acsl4 and gpx4 mRNA expression in ZMEL1 cells, cells were seeded into 6-well plates at the initial density of 1.13×106 per well. Once adhered to the plate, cells were treated with 20 μM of JA, PunA, α-ESA, and β-ESA, as well as with OLA and ALA or with DMEM culture medium (as a negative control) for 24 h. After the treatment, total RNA was extracted using the High Pure RNA Isolation Kit (11828665001, Roche). One microgram of RNA was reverse transcribed into cDNA using the iScript™ cDNA Synthesis Kit (1708891, Biorad) according to the manufacturer’s instructions, using an Applied Biosystems SimpliAmp Thermal Cycler (A24811, ThermoFisher). The synthesized cDNAs were amplified by real-time qPCR using the GoTaq qPCR mixture (Promega) on an Applied Biosystems StepOnePlus Real-Time PCR System (4376600, ThermoFisher) over 40 cycles. Primers used for the qPCR are listed in Supplementary Table 1. The expression levels of target genes were normalized to the expression of housekeeping genes (HKGs), namely β-actin 2 (ACTB2), β-2-microglobulin (B2M), Hypoxanthine-guanine phosphoribosyl transferase (HPRT1) and TATA-binding protein (TBP).
2.6 Statistical analysis
Data are expressed as mean ± standard error of the mean (SEM) of three independent replicates. Statistical analyses were performed with GraphPad Prism 10 software using one-way or two-way ANOVA with Dunnett’s test or Tukey multiple comparison test, when appropriate. Statistical significance relative to the control or another treatment was determined as follows: *p ≤ 0.05; **p ≤ 0.01; ***p ≤ 0.001; ****p ≤ 0.0001.
4. Discussion
Given their higher susceptibility to oxidation, as compared to their non-conjugated counterparts, CLnAs have shown promising anti-cancer effects in several cancer types. However, their potential activity against melanoma remains unexplored. Here we show for the first time that CLnAs exert significant cytotoxic effects on melanoma cells of both human (A375 and WM266.4) and zebrafish (ZMEL1) origin. Still, the tested human melanoma cell lines displayed a markedly higher sensitivity to CLnAs, as reflected by their IC50 values, which were approximately 4- to 10-fold lower than those observed for the zebrafish melanoma cell line under investigation.
To assess if CLnAs induce melanoma cell death through ferroptosis, as they do in other cancer cell models [
11], we utilized three distinct inhibitors targeting key mechanistic nodes of lipid peroxidation. Iron, as both a lipoxygenase cofactor and a Fenton reaction catalyst, was chelated using DFOM [
21,
45]. Peroxyl radical propagation was inhibited or prevented with the lipophilic antioxidant α-T and the alkoxyl radical scavenger fer-1 [
46,
47]. In contrast to inhibitors of necroptosis and apoptosis, all three ferroptosis inhibitors effectively restored melanoma cell viability, unequivocally identifying ferroptosis as the key cell death program triggered by CLnAs in melanoma cells.
The toxicity of CLnAs towards melanoma cells showed clear isomer-specific patterns. In human melanoma cell lines, PunA and JA were the most potent compounds, whereas α-ESA and β-ESA showed weaker but comparable toxic effects. In ZMEL1 cells, JA remained the most cytotoxic isomer, followed by PunA and α-ESA, while β-ESA exhibited only limited cytotoxicity. The greater cytotoxicity of JA and PunA may be linked to their c,t,c geometric configuration (8, 10, 12 for JA; 9, 11, 13 for PunA), which enhances their oxidative susceptibility [
46] and thereby their ability to induce ferroptosis. Nevertheless, the relationship between CLnA configuration and cytotoxicity appears to strongly depend on cell context. In breast cancer cells, Beatty
et al. reported potent ferroptosis induction by JA (IC50 1.8 μM) but notably weaker activity for PunA (IC50 19.4 μM). In DLD-1 human colorectal adenocarcinoma cells, JA was also found to be the most potent isomer after a 24-hour treatment, followed by α-ESA, PunA, and β-ESA [
48]. These findings indicate that the cytotoxicity differences among CLnA isomers are not restricted to melanoma cells. The disparity in sensitivity to CLnAs across studies may be linked to reported variations in the activity of enzymes involved in fatty acid incorporation into polar lipids, fatty acid oxidation, and ferroptosis across different types of cancer cells [
49,
50]. As an example, Beatty
et al. demonstrated that CLnA-induced ferroptosis in breast cancer cells mostly rely on ACSL1 [
23], whereas our results point to ACSL4 as a key contributor to ferroptosis induction in melanoma cells. Taken together, these findings highlight (i) the importance of isomer geometric configuration, which governs their intrinsic physicochemical properties, such as oxidative stability, and (ii) the influence of cell-specific enzymatic machineries, such as ACSL isoforms, in shaping cancer cell resistance to CLnAs.
The differential cytotoxicity of CLnA isomers was further substantiated by their distinct responses to ferroptosis inhibitors. In A375 and WM266.4 cells, comparable concentrations of ferroptosis inhibitors (i.e. DFOM, α-T, or fer-1) were required to restore viability in both JA- and PunA-treated cells, aligning with the similar toxicity levels of these two isomers in human melanoma cells. In contrast, ZMEL1 cells required higher inhibitor concentrations to mitigate JA-induced cell death than PunA-induced cell death, corroborating the higher cytotoxicity of JA in this zebrafish cell line. All these findings underscore the variable cytotoxicity of CLnA isomers towards melanoma cells, with JA and PunA emerging as the most potent inducers of ferroptosis in melanoma cells among the four tested CLnA isomers. This positions JA and PunA as promising candidates for future in vivo studies.
Ferroptosis is triggered by the accumulation of PUFA-derived lipid peroxides. ACSL4 indirectly facilitates their formation through the incorporation of these fatty acids into the phospholipids, while GPX4 reduces these peroxides and mitigates ferroptotic cell death [
51]. Our study demonstrates that CLnA toxicity in melanoma cells is significantly influenced by the functional activity of these two enzymes. Specifically, the inhibition of ACSL4 significantly reduced the toxicity of PunA, whereas inhibition of GPX4 markedly enhanced it.
The lipid analysis revealed that PunA and JA were effectively incorporated into PLs and NLs in ZMEL1 cells. Even though JA was more cytotoxic than PunA to the ZMEL1 cells, its accumulation in the PL fraction was lower. In contrast, both isomers extensively accumulated in the NL fraction, up to a similar extent after 12 hours of incubation. The significant incorporation of JA and PunA into triglycerides, presumably stored within lipid droplets (LDs), presents a complex mechanistic question. On one hand, LDs are often regarded as protective organelles that shield the PUFAs from lipid peroxidation [
52,
53]. Indeed, studies have shown that exceeding the buffering capacity of triglyceride storage in LDs leads to PUFA-induced ferroptosis [
40]. On the other hand, a growing body of evidence suggests that LDs can also act as pro-ferroptotic platforms [
52,
54]. For example, Beatty
et al. demonstrated that α-ESA accumulates in LDs in breast cancer cells where it may act as a potent source of lipid ROS, which can then propagate to cellular membranes [
23]. Similarly, Lange
et al. recently revealed that while LDs typically sequester PUFAs to prevent damage, the loss of the LD-localized antioxidant FSP1 transforms PUFA-rich LDs into sites of active lipid peroxidation that initiate ferroptosis [
54]. Considering melanoma cells, it remains to be determined whether CLnA cytotoxicity is mainly driven by their incorporation into NLs or whether their integration into PLs, where peroxidation can directly trigger ferroptosis, is the dominant mechanism
The greater cytotoxicity of JA to ZMEL1 cells as compared to PunA might be due to the higher susceptibility of JA to oxidation in the cellular environment [
18]. As a consequence, a small amount of JA in the membrane would be sufficient to trigger cell death. Alternatively, JA accumulation in LDs could serve as a localized hub for initiating the ferroptotic cascade, as previously proposed for α-ESA [
23]. Finally, a combination of both mechanisms cannot be ruled out based on available data.
The metabolic trafficking of CLnAs into either PLs or NLs is necessarily preceded by their activation to acyl-CoA esters, a process catalyzed by ACSL enzymes. We found that CLnA treatment led to a significant upregulation of
acsl4a expression when compared to cells treated with ALA or OLA, indicating a specific transcriptional response to these pro-ferroptotic lipids. Interestingly, a significant change was not observed when comparing to untreated control cells. This is not unprecedented, as a similar lack of significant
acsl4 transcriptional change has been reported in other models of PUFA-induced ferroptosis [
40]. This suggests that the basal enzymatic activity of ACSL4, rather than its dramatic transcriptional induction, may be the rate-limiting factor for mediating CLnA toxicity. This conclusion is supported by the fact that the inhibition of ACSL4 effectively mitigated the cytotoxic effects of CLnAs in melanoma cells, unequivocally highlighting the critical role of this enzyme in activating these fatty acids to drive ferroptosis.
Regarding
gpx4b, the dominant
gpx4 isoform in ZMEL1 cells, our results showed that its expression was not significantly altered by CLnA treatment compared to the control or to cells treated with OLA or ALA, apart from a slight upregulation induced by PunA compared to ALA. The minimal transcriptional alteration observed for
gpx4b aligns with findings by Beatty
et al. and Cuvelier
et al., who reported unaffected GPX4 protein levels in breast cancer cells treated with αESA or JA, respectively [
15,
40]. This collective evidence indicates that CLnAs neither induce ferroptosis by suppressing GPX4 expression nor trigger significant upregulation of this key antioxidant enzyme in response to the resulting lipid peroxidation. These findings actually align with the saturation of the antioxidant capacity of the GPX4 system which is ultimately overwhelmed by the tide of CLnA-derived lipid hydroperoxides, thereby leading to ferroptosis [
23,
55].
While GPX4 is a key oxidation defender, it is important to recognize that other antioxidant systems, such as the FSP1/CoQ10 axis and DHCR7, BH4-GCH1 also contribute to mitigating ferroptosis [
17,
56]. Thus, although CLnAs do not significantly alter
gpx4 expression in melanoma cells, the possibility remains that they modulate the expression levels of other players within this intricate antioxidant network.
Zebrafish has emerged as an important model in cancer research, particularly in melanomas [
57]. The consistent findings regarding cytotoxicity and the modulation of the ACSL4/GPX4 axis across both human and zebrafish melanoma cells underscore shared mechanisms, highlighting the translational relevance of the zebrafish model. Interestingly, the fluorescently labelled ZMEL1 cell line developed by Heilmann
et al. has the potential to track live formation of metastases in the transparent Casper zebrafish [
58]. Validating our findings using this zebrafish melanoma model would be a critical next step to observe tumour response to CLnA treatments.
Author Contributions
Conceptualization, Zhuo Zhang, Yvan Larondelle,Cathy Debier, Olivier Feron and Melissa Page; Data Curation, Zhuo Zhang; Formal Analysis, Zhuo Zhang; Funding acquisition, Yvan Larondelle; Investigation, Zhuo Zhang, Alice Valembois, Caroline Rosier and Renaud Bonnevie; Methodology, Zhuo Zhang, Ineke Neefs, Aurélien Warnant and Perrine Vermonden; Resources, Yvan Larondelle and Cathy Debier; Validation, Zhuo Zhang, Alice Valembois, and Caroline Rosier; Visualization, Zhuo Zhang; Writing–Original Draft Preparation, Zhuo Zhang; Writing–Review & Editing, Zhuo Zhang, Yvan Larondelle, Cathy Debier, Melissa Page, Perrine Vermonden and Olivier Feron.
Figure 1.
Effects of different fatty acids on the viability of melanoma cell lines of human (A375, WM266.4) or zebrafish (ZMEL1) origin. The A375 (a), WM266.4 (b) and ZMEL1 (c) cell lines were treated for 72 hours with varying concentrations of the following fatty acids: oleic acid (OLA), α-linolenic acid (ALA), jacaric acid (JA), punicic acid (PunA), α-eleostearic acid (α-ESA) and β-eleostearic acid (β-ESA). In addition, a kinetic study (d) was performed on ZMEL1 cells treated with PunA or JA at 10, 20, or 40 μM over 72 h.
Relative cell viability was normalized to control cells that were cultured with medium without added fatty acid, defined as 100%. Data are presented as mean ± standard error of the mean (SEM) of 3 independent experiments. Statistical significance was assessed using two-way ANOVA with Dunnett’s test, comparing different concentrations (5, 10, 20, 40, 80 μM) of fatty acid treatments against 0 μM in Figure 1a
-c and every time point against 0 h in Figure 1d
. p<0.05(*), p<0.01(**), p<0.001(***), p<0.0001(****).
Figure 1.
Effects of different fatty acids on the viability of melanoma cell lines of human (A375, WM266.4) or zebrafish (ZMEL1) origin. The A375 (a), WM266.4 (b) and ZMEL1 (c) cell lines were treated for 72 hours with varying concentrations of the following fatty acids: oleic acid (OLA), α-linolenic acid (ALA), jacaric acid (JA), punicic acid (PunA), α-eleostearic acid (α-ESA) and β-eleostearic acid (β-ESA). In addition, a kinetic study (d) was performed on ZMEL1 cells treated with PunA or JA at 10, 20, or 40 μM over 72 h.
Relative cell viability was normalized to control cells that were cultured with medium without added fatty acid, defined as 100%. Data are presented as mean ± standard error of the mean (SEM) of 3 independent experiments. Statistical significance was assessed using two-way ANOVA with Dunnett’s test, comparing different concentrations (5, 10, 20, 40, 80 μM) of fatty acid treatments against 0 μM in Figure 1a
-c and every time point against 0 h in Figure 1d
. p<0.05(*), p<0.01(**), p<0.001(***), p<0.0001(****).

Figure 2.
Effects of ferroptosis inhibitors on the CLnA-treated melanoma cells. Viability of A375, WM266.4 and ZMEL1 cells treated with punicic acid (PunA) or jacaric acid (JA) was assessed in the presence of increasing doses of ferrostatin-1 (fer-1; a-c), α-tocopherol (α-T; d-f), or deferoxamine mesylate (DFOM; g-i). Relative cell viability was normalized to control cells that were cultured with medium without any added inhibitor or CLnA, defined as 100%. For the cells cultured without CLnA, the culture medium still included the vehicles used for fer-1, α-T, and DFOM. Data are presented as mean ± standard error of the mean (SEM) of 3 independent experiments. Dose-response curves have been fitted to the data. Statistical significance was assessed using two-way ANOVA with Dunnett’s test, comparing different inhibitor concentrations against 0 μM. p<0.05(*), p<0.01(**), p<0.001(***), p<0.0001(****). Only significant differences are shown.
Figure 2.
Effects of ferroptosis inhibitors on the CLnA-treated melanoma cells. Viability of A375, WM266.4 and ZMEL1 cells treated with punicic acid (PunA) or jacaric acid (JA) was assessed in the presence of increasing doses of ferrostatin-1 (fer-1; a-c), α-tocopherol (α-T; d-f), or deferoxamine mesylate (DFOM; g-i). Relative cell viability was normalized to control cells that were cultured with medium without any added inhibitor or CLnA, defined as 100%. For the cells cultured without CLnA, the culture medium still included the vehicles used for fer-1, α-T, and DFOM. Data are presented as mean ± standard error of the mean (SEM) of 3 independent experiments. Dose-response curves have been fitted to the data. Statistical significance was assessed using two-way ANOVA with Dunnett’s test, comparing different inhibitor concentrations against 0 μM. p<0.05(*), p<0.01(**), p<0.001(***), p<0.0001(****). Only significant differences are shown.
Figure 3.
Effects of ACSL4 or GPX4 inhibitors on the PunA-induced cytotoxicity in melanoma cells. Viability of A375, WM266.4 and ZMEL1 cells treated with punicic acid (PunA) was assessed in the presence of increasing doses of PRGL493 (a-c), RSL3 (d-f) or ML210 (g-i). Relative cell viability was normalized to control cells that were cultured with medium without any added inhibitor or CLnA, defined as 100%. For the cells cultured without CLnA, the culture medium still included the vehicles used for PRGL493, RSL3 and ML210. Data are presented as mean ± standard error of the mean (SEM) of 3 independent experiments. Dose-response curves have been fitted to the data. Statistical significance was assessed by two-way ANOVA with Dunnett’s test, comparing different tested inhibitor concentrations (0.03, 0.1, 0.3, 1, 3, 10 μM) against 0 μM. p<0.05(*), p<0.01(**), p<0.001(***), p<0.0001(****). Only significant differences are shown.
Figure 3.
Effects of ACSL4 or GPX4 inhibitors on the PunA-induced cytotoxicity in melanoma cells. Viability of A375, WM266.4 and ZMEL1 cells treated with punicic acid (PunA) was assessed in the presence of increasing doses of PRGL493 (a-c), RSL3 (d-f) or ML210 (g-i). Relative cell viability was normalized to control cells that were cultured with medium without any added inhibitor or CLnA, defined as 100%. For the cells cultured without CLnA, the culture medium still included the vehicles used for PRGL493, RSL3 and ML210. Data are presented as mean ± standard error of the mean (SEM) of 3 independent experiments. Dose-response curves have been fitted to the data. Statistical significance was assessed by two-way ANOVA with Dunnett’s test, comparing different tested inhibitor concentrations (0.03, 0.1, 0.3, 1, 3, 10 μM) against 0 μM. p<0.05(*), p<0.01(**), p<0.001(***), p<0.0001(****). Only significant differences are shown.
Figure 4.
Cellular uptake and distribution of PunA, JA or ALA in ZMEL1 cells. Enrichment in punicic acid (PunA), jacaric acid (JA), and α-linolenic acid (ALA) in ZMEL1 cells treated with 20 μM of each fatty acid for 4, or 12 hours is shown in neutral lipids (NLs; a), phospholipids (PLs; b), and both fractions combined (NLs+PLs; c). All values at 0 hour were below the limit of detection. Data are presented as mean ± standard error of the mean (SEM) of 3 independent cultures and are expressed as the ratio of PunA, JA or ALA amount to the total amount of all fatty acids in each fraction (%, nmol/nmol). Statistical significance was assessed using two-way ANOVA with Tukey’s multiple comparisons. Statistical differences among PunA, JA, and ALA, at both 4 h and 12 h, are presented by the letters a, b, c and x, y, z above the bars, respectively.
Figure 4.
Cellular uptake and distribution of PunA, JA or ALA in ZMEL1 cells. Enrichment in punicic acid (PunA), jacaric acid (JA), and α-linolenic acid (ALA) in ZMEL1 cells treated with 20 μM of each fatty acid for 4, or 12 hours is shown in neutral lipids (NLs; a), phospholipids (PLs; b), and both fractions combined (NLs+PLs; c). All values at 0 hour were below the limit of detection. Data are presented as mean ± standard error of the mean (SEM) of 3 independent cultures and are expressed as the ratio of PunA, JA or ALA amount to the total amount of all fatty acids in each fraction (%, nmol/nmol). Statistical significance was assessed using two-way ANOVA with Tukey’s multiple comparisons. Statistical differences among PunA, JA, and ALA, at both 4 h and 12 h, are presented by the letters a, b, c and x, y, z above the bars, respectively.
Figure 5.
Modulation of acsl4a and gpx4b gene expression in ZMEL1 cells treated with fatty acids. The fold changes in relative quantity (RQ) were calculated for acsl4a (a), or gpx4b (b) in cells exposed to 20 μM of a tested fatty acid for 24 hours, relative to the control cells which cultured with DMEM without any added fatty acid. The following fatty acids were tested: α-linolenic acid (ALA), oleic acid (OLA), punicic acid (PunA), α-eleostearic acid (α-ESA), β-eleostearic acid (β-ESA) and jacaric acid (JA). The RQ was calculated as RQ = 2−ΔCt, the ΔCt value representing the difference in Ct value between the target gene and a panel of reference genes, which are actin beta 2, beta-2-microglobulin, hypoxanthine phosphoribosyl transferase 1 and TATA-binding protein. Data are presented as mean ± standard error of the mean (SEM) of 3 independent experiments (N=3, n=3). Significance was assessed by one-way ANOVA with Tukey’s multiple comparison test among treatments. Significant differences between the fatty acid-treated groups and the control are marked above the columns. p<0.05(*); p<0.01(**); p<0.001(***), p<0.0001(****). Only significant differences are shown.
Figure 5.
Modulation of acsl4a and gpx4b gene expression in ZMEL1 cells treated with fatty acids. The fold changes in relative quantity (RQ) were calculated for acsl4a (a), or gpx4b (b) in cells exposed to 20 μM of a tested fatty acid for 24 hours, relative to the control cells which cultured with DMEM without any added fatty acid. The following fatty acids were tested: α-linolenic acid (ALA), oleic acid (OLA), punicic acid (PunA), α-eleostearic acid (α-ESA), β-eleostearic acid (β-ESA) and jacaric acid (JA). The RQ was calculated as RQ = 2−ΔCt, the ΔCt value representing the difference in Ct value between the target gene and a panel of reference genes, which are actin beta 2, beta-2-microglobulin, hypoxanthine phosphoribosyl transferase 1 and TATA-binding protein. Data are presented as mean ± standard error of the mean (SEM) of 3 independent experiments (N=3, n=3). Significance was assessed by one-way ANOVA with Tukey’s multiple comparison test among treatments. Significant differences between the fatty acid-treated groups and the control are marked above the columns. p<0.05(*); p<0.01(**); p<0.001(***), p<0.0001(****). Only significant differences are shown.
