Next Article in Journal
Orchestrating Cellular Balance: ncRNAs and RNA Interactions at the Dominant of Autophagy Regulation in Cancer
Next Article in Special Issue
Molecular Aspects of Inflammation and Lipid Metabolism in Health and Disease: The Role of the Mitochondria
Previous Article in Journal
Factors Influencing Venous Remodeling in the Development of Varicose Veins of the Lower Limbs
Previous Article in Special Issue
Adipocyte- and Monocyte-Mediated Vicious Circle of Inflammation and Obesity (Review of Cellular and Molecular Mechanisms)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Alterations in Plasma Lipid Profiles Associated with Melanoma and Therapy Resistance

1
Clinical Biochemistry and Mass Spectrometry Laboratory, Health Sciences Department, Università degli Studi di Milano, 20122 Milan, Italy
2
Unit of Bioinformatics and Biostatistics, Department of Epidemiology and Data Science, Fondazione IRCCS Istituto Nazionale dei Tumori, 20133 Milan, Italy
3
Unit of Translational Immunology, Department of Experimental Oncology, Fondazione IRCCS Istituto Nazionale dei Tumori di Milano, 20133 Milan, Italy
4
Department of Diagnostic and Technology, Fondazione IRCCS Istituto Neurologico Carlo Besta, 20133 Milan, Italy
5
Molecular Pharmacology Unit, Department of Experimental Oncology, Fondazione IRCCS Istituto Nazionale dei Tumori, 20133 Milan, Italy
6
Laboratory of Neurobiology and UCV, Neurology IX Unit, Fondazione IRCCS Istituto Neurologico Carlo Besta, 20133 Milan, Italy
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
These authors contributed equally to this work.
Int. J. Mol. Sci. 2024, 25(3), 1558; https://doi.org/10.3390/ijms25031558
Submission received: 28 December 2023 / Revised: 17 January 2024 / Accepted: 24 January 2024 / Published: 26 January 2024

Abstract

:
Dysfunctions of lipid metabolism are associated with tumor progression and treatment resistance of cutaneous melanoma. BRAF/MEK inhibitor resistance is linked to alterations of melanoma lipid pathways. We evaluated whether a specific lipid pattern characterizes plasma from melanoma patients and their response to therapy. Plasma samples from patients and controls were analyzed for FASN and DHCR24 levels and lipidomic profiles. FASN and DHCR24 expression resulted in association with disease condition and related to plasma cholesterol and triglycerides in patients at different disease stages (n = 144) as compared to controls (n = 115). Untargeted lipidomics in plasma (n = 40) from advanced disease patients and controls revealed altered levels of different lipids, including fatty acid derivatives and sphingolipids. Targeted lipidomics identified higher levels of dihydroceramides, ceramides, sphingomyelins, ganglioside GM3, sphingosine, sphingosine-1-phosphate, and dihydrosphingosine, saturated and unsaturated fatty acids. When melanoma patients were stratified based on a long/short-term clinical response to kinase inhibitors, differences in plasma levels were shown for saturated fatty acids (FA 16:0, FA18:0) and oleic acid (FA18:1). Our results associated altered levels of selected lipid species in plasma of melanoma patients with a more favorable prognosis. Although obtained in a small cohort, these results pave the way to lipidomic profiling for melanoma patient stratification.

1. Introduction

Cutaneous melanoma displays a rising worldwide incidence, but in the past decade, its clinical management has been revolutionized with the introduction of BRAF and MEK inhibitors (BRAFi and MEKi) and of immunotherapy with antibodies to immune checkpoint inhibitors, achieving a dramatic improvement of patient overall survival. Several new drugs and therapeutic options have been studied and approved, together with drug combination strategies and adjuvant and neoadjuvant therapy settings, showing promising safety and efficacy results compared to conventional chemotherapy [1,2].
Metastatic disease progression and drug resistance are associated with metabolic reprogramming and dysregulated lipid metabolism [3,4]. Alterations of lipid regulatory pathways, including pathways controlling cholesterol (CHOL), fatty acid (FA) synthesis [5,6,7], and sphingolipids (SL) metabolism [8,9], have been associated with disease progression and resistance to kinase inhibitors (KI), both in patients’ tumors and in cell lines model systems.
We recently reported that transcriptomes of BRAFi/MEKi-resistant melanoma are characterized by dysregulated lipid metabolism pathways involving FA and CHOL [10] and determining alterations in the lipid composition, particularly in different classes of FA, cholesteryl esters (CE) and triglycerides (TG), along with modulated expression of enzymes regulating biosynthetic nodes of the lipid metabolism [11]. The effect of tackling lipid metabolism pathways in resistant cell lines was evidenced by lipid starvation, which reduced cell growth and increased sensitivity to the BRAFi PLX4032/vemurafenib [11]. Specifically, we showed that the combination of BRAFi with inhibitors of lipogenic enzymes fatty acid synthase (FASN) and 24-dehydrocholesterol reductase (DHCR24) resulted in favorable drug interactions, potentially exploitable for the treatment of BRAFi-resistant melanoma [11,12].
Alterations in lipid profiles, which elucidate disease mechanisms and represent valuable biomarkers, are expected to have a pivotal role in patients’ risk stratification as well as in precision medicine. These features allow for the setting up of medical treatments based on individual lipid profiles, namely the ‘lipidome’. Specifically, the ‘lipidome’ is dependent on two factors: the biological properties of cells, which reflect the genomic alterations in biological systems, and the environmental impact on those biological systems. By capturing the complex, multilayered landscape of the lipidome as well as tracking changes in lipid levels and composition, lipidomics can reveal distinct molecular phenotypes that indicate the physiological and pathological state of an organism [13,14,15,16].
In this explorative study, we aimed to evaluate the involvement of regulatory lipid pathways at the plasma level in melanoma and in drug resistance. We quantitatively assessed the levels of FASN and DHCR24 enzymes in melanoma patients and healthy control subjects in relation to plasma TG and CHOL and tested plasma lipid composition by untargeted/targeted lipidomic approaches. In addition, to uncover potential alterations of lipid species that specifically characterize the signature of melanoma KI-resistance, we quantitively determined the lipid species in plasma from patients showing a long or a short response to targeted therapy.

2. Results

2.1. Plasma Levels of DHCR24 and FASN Enzymes Are Altered in Melanoma Patients Compared to Control Subjects

Our previous studies showed that transcriptomes of BRAFi/MEKi-resistant melanoma cells are characterized by dysregulated lipid metabolism pathways involving FA and CHOL, determining alterations in the lipid composition, along with modulated expression of enzymes regulating biosynthetic nodes of the lipid metabolism. FASN and DHCR24 lipogenic enzymes were shown to be involved in the modulation of BRAFi response in melanoma [11,12]. Based on these findings, we explored the plasma levels of FASN and DHCR24 as potential biomarkers in melanoma patients and tested the association between the expression levels of the two enzymes and the different clinical-demographic features of melanoma patients, including blood total CHOL, high-density lipoprotein (HDL), low-density lipoprotein (LDL) and TG. Quantitative analysis was carried out in a cohort of 144 melanoma patients at different clinical stages (stage I–IV) and 115 control subjects, which comprised individuals at surgery for benign skin lesions and healthy blood donors. Table 1 reports the main clinical-demographic characteristics of the study cohort.
By considering all the tested plasma samples, the levels of DHCR24 (ANOVA p-value 0.019) and FASN (ANOVA p-value 0.011) were statistically associated with disease condition (Figure 1A,B).
In detail, by looking at individuals undergoing surgery for skin lesions suspected of melanoma, in subjects diagnosed with melanoma by pathological analysis (PT DH I/II), the plasma levels of FASN were shown to be significantly higher compared to all the other groups (Tukey p-values < 0.01, Figure 1C). Conversely, an opposite trend was shown for DHCR24 levels. The only significant association after multiple comparison adjustments was shown in patients at advanced stages (PT III/IV) vs. subjects at surgery of benign skin lesions (control DH) (Tukey p-value 0.017, Figure 1D). By looking at the logistic regression results, the odds ratio for FASN and DHCR24, as well as for other biochemical variables in a univariate fashion, showed significance for both FASN and DHCR24 (Table 2), although when adjusting for age, neither of these markers retained the statistical significance. No correlation between FASN and DHCR24 levels was observed. In addition, by looking at the pair-wise correlations of FASN and DHCR24 with the other variables, we observed a low magnitude pattern of correlation: FASN levels were positively correlated with CHOL, LDL, and TG, whereas DHCR24 levels were negatively correlated to TG (Figure S1).
We next analyzed the potential association of FASN and DHCR24 with treatment response duration in a group of advanced metastatic patients at the baseline of medical treatments resulting in a long- (long-term responders, LR) or a short- (short-term responders, SR) term clinical response (Table 1). Higher plasma levels of FASN were observed in LR patients (Figure 1E) with a borderline significance (OR: 1.571; 95%CI: 0.976; 2.528), whereas no association was observed between DHCR24 and treatment response (Figure 1F) and with the other considered plasma variables (Table S1). The same trend was observed in patients treated with BRAFi/MEKi (KI, N = 40) and in those treated with immunotherapy with immune checkpoint inhibitors (ICI, N = 40) (Figure 1G). The opposite behavior was observed for DHCR24, showing higher levels in SR patients (Figure 1H), while no significant association was observed between the treatment response duration and the other plasma variables considered.

2.2. Lipidomic Profiling of Plasma Reveals Altered Levels of Fatty Acids Derivatives and of Sphingolipids in Melanoma Patients

To evaluate potential alterations of systemic lipid metabolism in melanoma progression, we carried out an untargeted lipidomic analysis in 20 plasma samples from advanced melanoma patients at the baseline of targeted therapy with KI and in 20 sex-matched control subjects (control DH). Among the 935 lipid species recognized, case-control analysis at a single lipid level identified four lipids associated (after false discovery rate—FDR-adjustment) with the disease condition and showing lower values in patients, including glycerolipids and glycerophospholipids TG 56:0 (or TG 16:0/18:0/22:0), DG 25:4, DG 50:3, and PC O-39:10 (Figure 2A).
At the lipid family level, case-control analysis highlighted three families with a different expression pattern between cases and control, with free FA showing higher levels in patients compared to controls and hydroxylated phytosphingolipids Cer AP (alpha-hydroxy fatty acid phytoceramide) and HexCer AP (alpha-hydroxy fatty acid phytohexosylceramide) showing lower levels in patients (Figure 2B). In LR compared to SR patients, diacylglycerol (DG 32:5) and FA derivatives N acyl glycil serine (NAGly Ser) 25:0, CAR 26:0, and vitamin A fatty acid ester (VAE) 13:0 showed significantly higher values, while sphingomyelin SM 33:2 was under-expressed (Figure 2C). None of the latter was included in the list of lipids differentially expressed between cases and controls, suggesting a possible involvement of different lipids in the two scenarios.

2.3. Targeted Analysis of Fatty Acid and Sphingolipid Species Identify Quantitative Differences in Patients Responsive to Targeted Therapy

To deepen the analysis, we performed a targeted lipidomic quantification focused on FA and SL profiles. What appeared striking about the complex plasma SL profile is their diffuse increased concentrations in melanoma patients compared to control subjects (Figure 3A). The more remarkable differences were appreciated in the total concentrations of dihydroceramides (DHCer, FRD-W p-value: 0.006), ceramides (Cer, FRD-W p-value: 0.004), sphingomyelins (SM, FRD-W p-value: 0.001), and gangliosides (GM3, FRD-W p-value: <0.001). When we compared samples from patients classified according to response to targeted therapy, the differences in the SL pathway were less marked and lacked statistical significance (Figure 3B).
Relatively to sphingoid bases, a sharp difference was found in the concentration of the catabolic products sphingosine (Sph, FRD-W p-value: 0.009) and sphingosine-1-phosphate (S1P, FRD-W p-value: 0.011) as well as dihydrosphingosine (dhSPh, FRD-W p-value: <0.001), which is an early intermediate of the de-novo SL biosynthesis (Figure 4A). By contrast, no significant differences were found between SR and LR.
For the analysis of the plasma levels of free FA, we focused only on FA arising from the direct or related activity of FASN, which drives the synthesis of palmitate. The pattern found for the SL concentrations was also observed for the FA plasma content, as higher levels are found in melanoma patients (Figure 4B). This trend reaches a borderline significance after FRD adjustments for stearic FA (FA 18:0, FRD-W p-value: 0.065) and for the unsaturated derivatives oleic (FA 18:1, FRD-W p-value: 0.065), linoleic (FA 18:2, FRD-W p-value: 0.058), linolenic (FA 18:3, FRD-W p-value: 0.025) and palmitoleic acids (FA 16:1, FRD-W p-value: 0.072), while no significance was instead observed for palmitic acid (FA 16:0, FRD-W p-value: 0.236). In patients stratified based on treatment response, a trend towards different FA plasma levels was shown for the saturated FA 16 and FA 18 (W p-value: 0.014) and oleic acid (FA 18:1, W p-value: 0.11), which were higher in the group of responsive patients although no significances after FDR adjustments were reached (Figure 4C). Of note, superimposable results were obtained when lipid families were computed according to the Principal Component Approach (PCA).

3. Discussion

Several lines of evidence implicate the dysregulation of lipid metabolism in melanoma progression and therapy resistance [17]. Subcutaneous adipose tissue and released adipokines were suggested to promote melanoma growth [18]. Melanoma cells were shown to be able to take up lipids from extracellular vesicles released by stromal adipocytes or from aged fibroblasts by fatty acid transport protein (FATP) and use those to survive upon targeted therapy [19,20,21]. Dysregulation of CHOL homeostasis in melanoma cells was shown to impact proliferation, migration, and invasion, and overexpression of CHOL synthesis genes was associated with poor disease prognosis and drug resistance [10,22]. High prostaglandin E2 (PGE2) synthesis and increased FA oxidation are associated with poor prognosis in melanoma patients after targeted therapy with BRAFi/MEKi [6,21,23]. The expression of genes implicated in FA metabolism was shown to be associated with the expression of melanocyte inducing transcription factor (MITF), a main regulator of the proliferative to invasive phenotype switch in melanoma cells, and to AXL receptor, a hallmark of a de-differentiated and drug-resistant phenotype [2,7].
As reported in previous studies carried out in our lab, vemurafenib-resistant melanoma cell line models displayed increased levels of key enzymes controlling lipid metabolism [10,12]. Drug-induced cell selection importantly impacts on lipid assets, as resistant cell lines revealed an altered lipid composition, characterized by a general increase in mono-unsaturated fatty acids (MUFA) counterbalanced by a decrease in SFA in FA as well as in TG and CE. In addition, a trend for the increase in FA typical of SM (C20:0, C20:1, C20:2, C22:0, C24:0, and C24.1) was also evidenced [11]. When quantitatively evaluated by SL targeted analysis, the total amount of Cer (d18:1, FA from C14 to 24) and SM (d18:1, FA from C16 to 24) as well as Sph d18:1 was higher in resistant cells compared to their sensitive counterparts and appeared modulated by drug treatment. Several studies showed that SL deregulated metabolism results in the accumulation of metabolites endowed with a pro-tumoral action, including both simple Cer derivatives, including Sph, as well as complex SL, such as sphingomyelin and GM3 [9,24].
Lipidomics is the most powerful analytical tool for studying lipids in biological samples and their biochemical involvement in human diseases. In particular, bioactive lipids are considered novel biomarkers that may contribute to innovative medical treatments, early intervention, and improved patients’ prognoses. Therefore, the identification of lipidomic signatures may be relevant for the detection of potential diagnostic, prognostic, and predictive biomarkers of cancer [25]. In melanoma, examining the concentrations of a subset of lipids may shed light on cellular processes and reveal changes in targetable pathways [26].
Recently, the profiles of plasma lipidome in cancer patients have gained attention for potential translation from basic research to clinical application [27]. In this pilot experimental study, we evaluated the levels of lipogenic enzymes and lipid profiles in plasma from melanoma patients compared to healthy controls and their potential association with clinical response to treatment. As the levels of the lipogenic enzymes FASN and DHCR24 differed in melanoma patients with respect to controls, we proceeded further to lipidomic profiling studies. To this purpose, targeted and untargeted lipidomic approaches were applied. Untargeted lipidomics indicated that three major lipid families showed quantitative differences in patients, including phytosphingolipids (CerAP and HexCerAP) and FA, with the latter showing an enrichment (Figure 2B). A targeted approach towards complex SL, sphingoid bases, and free FA (saturated and unsaturated) families was then applied to narrow the analysis of these specific lipid species regulated in the CHOL and FA biosynthetic pathways. We found higher levels of SL DHCer, Cer, SM, HexCer, LacCer, and GM3, sphingoid bases (e.g., Sph, S1P, and dhSph) and of free FA (e.g., FA 16:1, 18:0, 18:1, 18:2, 18:3) in patients compared to controls (Figure 3 and Figure 4). Concerning the FA, higher levels of FA 16:0, 18:0, and 18:1 were evidenced in LR in comparison to SR (Figure 3 and Figure 4). Notably, with respect to SR, LR patients showed a trend for higher levels of DHCer, Cer, SM, and LacCer, while the levels of GM3 revealed an opposite trend with higher values in SR. Gangliosides, including the simplest GM3, are present in melanoma cells in high amounts and have been associated with melanoma progression [28]. The higher levels of GM3 found in melanoma patients as compared to controls and higher levels evidenced in poorly responsive patients suggest a possible different expression in melanoma.
In summary, the levels of specific SL, sphingoid bases, and saturated/unsaturated free FA defined a lipid pattern distinguishing advanced melanoma patients. This specific lipid pattern, together with the plasma level of FASN and DHCR24, stratifies advanced patients in SR and LR, suggesting a potential predictive value for treatment response. Although interesting, these results were obtained in a small sample set and require further validation studies in larger cohorts.
Our results are in agreement with data reported by Garandeu [29], showing higher plasma levels of individual long-chain Cer species C22:0, C24:0, and C24:1 in three patients classified as responders compared to 3 non-responders to BRAFi, suggesting that high levels of SL may distinguish responsive patients. While no evidence for a major role for plasma lipids has been shown in melanoma risk [30,31], emerging data indicate that blood metabolites impact immune cell functions and clinical response to immunotherapy in melanoma patients [32]. The association of pre-treatment serum lipid levels and response to anti-PD-1 therapy was also reported in advanced intrahepatic cholangiocarcinoma [33]. A high body mass index, identifying overweight and obesity conditions, has been associated with a marked survival advantage after both targeted and immunotherapy in male melanoma patients [34], highlighting a potential link between factors associated with obesity, such as free FA, insulin/IGF1, and proinflammatory cytokines, and the activation of an effective anti-tumor immune response. Evidence has been provided that the anti-tumor activity of immune cells is facilitated by lipids, which are essential fuels and metabolic components [35]. In particular, FA metabolism processes are crucially involved in the survival and function of immune cells [36]. In this light, the higher levels of FA and SL we detected in LR compared to SR may define a favorable condition promoting targeted therapy-induced anti-tumor immune response.
Overall, our analysis of plasma lipid composition by targeted and untargeted approaches focusing on phospholipids reveals for the first time a different lipid pattern in melanoma patients and its potential predictive value for treatment response, suggesting that lipidomic profiles may uncover novel prognostic biomarkers for melanoma. However, this is a small explorative study suffering from some limitations. Specifically, the limitations of the study are the small number of patients and the lack of a validation cohort to support the obtained results. In addition, the retrospective nature of the present research requires the conduction of future prospective studies for further validation of its clinical applicability, for which the development of advanced and simplified techniques for lipid measurement is also necessary.

4. Materials and Methods

4.1. Plasma Samples

The study was approved by the Independent Ethic Committee and by the Institutional Review Board (protocol INT 94/18), and blood samples from cutaneous melanoma patients and control subjects were obtained upon written informed consent. The cohort of melanoma patients included groups at different clinical stages as reported in Table 1: (a) patients at early disease stages, at surgery for cutaneous primary melanoma (stage I-II, PT DH I/II, n = 36) or for regional lymph node metastases (stage III, PT III, n = 28) and (b) advanced stage metastatic patients (stage III–IV, PT III–IV) at the baseline of medical treatment with BRAF/MEK KI (n = 40) or with immune checkpoint inhibitors (ICI) (n = 40). Short-term (SR, ≤12 months) or long-term (LR, ≥13 months) responses to medical treatment in this class of patients were defined based on clinical controls during therapy. The series of control subjects included (a) patients at surgery for benign skin lesions (CTR DH, n = 75) and (b) blood donors from the hospital blood transfusion service (CTR HD, n = 40). Lipidomic analyses were carried out on a subgroup of samples from 20 patients with stage III–IV melanoma at baseline of KI therapy and 20 controls (CTR DH) matched for sex. Plasma samples were obtained from blood collected from fasting subjects in 10 mL vacutainer K2EDTA tubes [Becton Dickinson (Franklin Lakes, NJ, USA)] and centrifuged at 2000× g for 10 min at RT within 1h. To eliminate cell debris, plasma was centrifuged again at 2000× g for 10 min at RT before storage in aliquots at −80 °C until use. Samples were analyzed for levels of triglycerides (TG), total cholesterol (CHOL), and high-density lipoprotein (HDL) cholesterol by standard clinical laboratory methods; low-density lipoprotein (LDL) cholesterol concentrations were estimated as CHOL − (HDL + TG/5).

4.2. Measurement of FASN and DHCR24 Plasma Levels

Plasma levels of FASN and DHCR24 were quantified using ELISA as instructed by the manufacturer’s protocols (FASN, NBP2-82512, Novus Biologicals, Centennial, CO, USA; DHCR24, KTE62074, Abbkine, Wuhan, China). Plasma samples were tested undiluted. The limits of detection for FASN and DHCR24 were 0.31–20 ng/mL and 15–240 pg/mL, respectively. A microplate reader (Tecan Infinite M1000, Zürich, Switzerland) was adopted to determine the absorbance at 450 nm.

4.3. Chemicals and Reagents for Lipidomics

We purchased the following chemicals from Sigma-Aldrich (St. Louis, MO, USA): acetonitrile, 2-propanol, methanol, chloroform, formic acid, ammonium acetate, ammonium formate, and dibutylhydroxytoluene (BHT) were purchased. Purified water at a Milli-Q grade (Burlington, MA, USA) was used to prepare all aqueous solutions.

4.4. Untargeted Lipidomics

Lipids from plasma (25 µL) were analyzed by LC-MS/MS after dilution with water (75 µL) and extraction by a mixture of methanol/chloroform (850 µL, 2:1 v/v). LC-MS/MS consisted of a Shimadzu UPLC coupled with a Triple TOF 6600 Sciex (Concord, Toronto, ON, Canada) [37]. Analysis was in duplicate for all samples in positive and negative electrospray ionization. Acquisition of spectra was also carried out by full-mass scan from m/z 200–1500 and top-20 data-dependent acquisition from m/z 50–1500. A value of 50 eV was fixed for declustering potential, with a collision energy of 35 ± 15 eV. A reverse-phase Acquity CSH C18 column 1.7 μm, 2.1 × 100 mm (Waters, Franklin, MA, USA) was used for chromatographic separation by a gradient between (A) water/acetonitrile (60:40) and (B) 2-propanol/acetonitrile (90:10), both containing 10 mM ammonium acetate and 0.1% of formic acid.

4.5. Lipidomics Data Processing

The freeware software MS-DIAL (ver. 4.0) was used to attain the spectra deconvolution, peak alignment, blank filtering, and sample normalization. Values of 0.01 and 0.05 Da were set for MS and MS/MS tolerance for peak profile, respectively. By matching spectra with the LipidBlast database or in-house built mass spectral library identification was achieved. The Lawless algorithm was employed to normalize the intensities of analytes, excluding those with a CV% superior to 30% in the QC pool sample. The following lipid classes were considered: acylcarnitines (ACar), cholesterol (CHOL), cholesterol esters (CE), dihydroceramides (DHCer), ceramides (Cer), hexosylceramides (HexCer), lactosylceramides (LacCer), globotriaosylceramide (Gb3), sulfatides (SULF), gangliosides (GM3), sphingomyelins (SM), vitamin A fatty acid ester (VAE), N acyl glycil serine (NAGlySer), lysophosphatidylcholines (LPC), lysophosphatidylethanolamines (LPE), phosphatidylcholines (PC), phosphatidylethanolamine (PE), phosphatidylinositoles (PI), and plasmalogens, that are ether-linked phosphatidylcholines (EtherPC) and vinyl linked phosphatidylethanolamines (EtherPE). The total number of carbon atoms and degree of unsaturation (i.e., PC 40:2) or specific acyl chains detected (i.e., PC 18:0/18:1) are considered to name lipids.

4.6. Targeted Lipidomics for Sphingolipids and Free Fatty Acids Quantification

The extraction of sphingolipids (SL) was performed with the same method used for untargeted lipidomics, with the addition of an alkaline methanolysis step with the aim of ameliorating low abundant SL species recovery. An LC-MS/MS system consisting of a LC Dionex 3000 UltiMate (ThermoFisher Scientific, Waltham, MA, USA) coupled to a tandem mass spectrometer AB Sciex 3200 QTRAP (AB Sciex, Concord, Toronto, ON, Canada) equipped with electro spray ionization TurboIonSpray™ source operating in positive mode (ESI+) was employed to inject the methanol-based supernatant. Mobile phases (A) water + 0.2% formic acid + 2 mM ammonium formate and (B) methanol + 0.2% formic acid + 1 mM ammonium formate were used for chromatography [38]. An Acquity BEH C8 column 1.7 μm, 2.1 × 100 mm (Waters, MA, USA) was employed for the complex SL (DHCer, Cer, SM, HexCer, LacCer and GM3), whereas a Cortecs C18 1.6 μm, 2.1 × 100 mm (Waters, MA, USA) for sphingoid bases (Sph, S1P, dhSph and dhS1P), as previously detailed [38]. In addition, in the case of free FA, plasma (50 µL) deproteinization was carried out with 100 µL of cold isopropanol added with 10 µL of IS (2.5 µg/mL isobuthoxyacetic acid, 10 µg/mL undecanoic acid and 20 µg/mL heptadecanoic acid). The samples were centrifuged (12,000× g, 5 min), and supernatants (100 µL) were separated in glass vials for the derivatization of the carboxylic function. An aliquot of 50 µL of 50 mM of 3-NPH, 50 µL of 50 mM of EDC, and 50 µL of pyridine (7%) was added to the latter. The agents used for derivatization were prepared in a solution of 70% methanol. Derivatization was carried out in an incubator at 37 °C for up to 30 min. The solution was diluted with 250 µL 0.5% of formic acid in isopropanol and directly injected into the above-described apparatus. Mass spectrometry operated in negative electrospray ionization by multiple reaction monitoring mode and chromatographic separation was achieved on a Restek Raptor C18 2.7 µm 2.1 × 100 mm (Bellefonte, PA, USA) using mobile phase (A) water + 0.1% formic acid and (B) acetonitrile. The following elution program (%B) was used: 0–2.5 min 10%, 2.5–31 min 10–99%, 31–35 99%, 35–35.2 99–10% maintained until 42 min; with the flow rate of 0.4 mL/min; the column and the autosampler temperature were 35 °C and 15 °C.

4.7. Statistical Analysis

The main characteristics of the study cohort were described by standard descriptive statistics, medians and ranges for continuous variables, and frequency tables for categorical variables. For continuous variables, the one-way analysis of variance was adopted to assess the relationship between the continuous variables and the different study outcomes. To adjust for group comparisons the Turkey correction was adopted. To assess the association between FASN and DHCR24 plasma levels with the disease status (cases vs. controls), logistic regression models [39] were also performed, and results were reported in terms of odds ratio (OR) and 95% confidence intervals (95%CI). The Pearson correlation coefficient (r) and its 95%CI were used to assess the strength of the association of the considered biochemical variables [40]. For untargeted lipidomic analyses, data were checked for integrity, filtered by interquartile range, and log-transformed before the statistical analysis. The Wilcoxon exact test (W) [41], followed by false discovery rate (FDR) adjustment, was adopted to highlight differentially expressed lipids in melanoma patients and controls as well as between treatment responses (LR vs. SR). The same approach was adopted for the analysis at the lipid family level. For the latter, the principal component approach (PCA) [42] was used to combine lipids of the same family. For targeted lipidomic analysis, plasma concentrations expressed as µM were considered and analyzed according to the W-test. For SL family analysis, the sum of the lipids with FA between C16 and C24 was considered, whereas that of saturated and unsaturated free FA between C16 and C18 was considered for the sphingoid bases and their phosphate form. A w-test followed by FDR adjustment was applied to determine families differentially expressed between melanoma patients and controls as well as between treatment responses. SAS software (Version 9.4, SAS Institute Inc., Cary, NC, USA), adopting a nominal level of 5%, was used for all statistical analysis, and R-software (R version 4.2.1, Foundation for Statistical Computing, Vienna, Austria) with the ggplot2 package (version 3.4.2) for graphical representations.

5. Conclusions

The identification of changes in lipid profiles, providing valuable biomarkers, and shedding light on mechanisms of diseases could have a crucial role in patients’ risk stratification and precision medicine, allowing for the selection and customization of medical treatment based on individual lipid profiles. In this pilot study, we evaluated whether a specific lipid pattern characterizes plasma from melanoma patients and their response to targeted therapy with KI. Taken together, our results showed that altered levels of selected lipid species are detectable in melanoma patients’ plasma specifically in patients with a more favorable prognosis. Although these results were obtained in a small patient cohort, they paved the way for the hypothesis that lipidomic profiles may uncover novel prognostic biomarkers for melanoma.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms25031558/s1.

Author Contributions

Conceptualization, M.D.C., C.M.C., E.V., P.V., R.G., R.P., P.P., L.G. and M.R.; methodology, E.V., E.C., M.A. and G.L.B.; software, C.M.C., M.R. and V.D.; validation, E.V., E.C. and M.A.; formal analysis, M.D.C., C.M.C. and V.D.; writing—original draft preparation, M.D.C., C.M.C. and G.L.B.; data curation, M.D.C., C.M.C. and V.D., supervision, P.V., R.G., R.P., P.P., L.G. and M.R.; funding acquisition, L.G. and M.R.; writing—review and editing, L.G. and M.R. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Italian Ministry of Health (grant RF2016-02361091 to L.G. and RCR2021-ACC Melanoma Working Group to M.R.).

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki and approved by the Independent Ethic Committee and by the Institutional Review Board (protocol INT 94/18, approval date 28 May 2018). Blood samples from cutaneous melanoma patients and control subjects were obtained upon written informed consent.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The original contributions presented in the study are included in the article; further inquiries can be directed to the corresponding author.

Acknowledgments

We thank Lorenza Di Guardo (Medical Oncology), Elena Tolomio and Aldo Bono (Surgical Day Hospital), Flavio Arienti (Transfusional Medicine) for the accrual of melanoma patients and control subjects, Felicetta Giardino and Martina Stroscia study nurses for blood samples, and Simona Frigerio and Paola Deho for sample processing and banking.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Villani, A.; Scalvenzi, M.; Micali, G.; Lacarrubba, F.; Fornaro, L.; Martora, F.; Potestio, L. Management of Advanced Invasive Melanoma: New Strategies. Adv. Ther. 2023, 40, 3381–3394. [Google Scholar] [CrossRef]
  2. Seth, R.; Agarwala, S.S.; Messersmith, H.; Alluri, K.C.; Ascierto, P.A.; Atkins, M.B.; Bollin, K.; Chacon, M.; Davis, N.; Faries, M.B.; et al. Systemic Therapy for Melanoma: ASCO Guideline Update. J. Clin. Oncol. 2023, 41, 4794–4820. [Google Scholar] [CrossRef]
  3. Falletta, P.; Goding, C.R.; Vivas-García, Y. Connecting metabolic rewiring with phenotype switching in melanoma. Front. Cell Dev. Biol. 2022, 10, 930250. [Google Scholar] [CrossRef]
  4. Pellerin, L.; Carrié, L.; Dufau, C.; Nieto, L.; Ségui, B.; Levade, T.; Riond, J.; Andrieu-Abadie, N. Lipid metabolic Reprogramming: Role in Melanoma Progression and Therapeutic Perspectives. Cancers 2020, 12, 3147. [Google Scholar] [CrossRef]
  5. Wilcock, D.J.; Badrock, A.P.; Wong, C.W.; Owen, R.; Guerin, M.; Southam, A.D.; Johnston, H.; Telfer, B.A.; Fullwood, P.; Watson, J.; et al. Oxidative stress from DGAT1 oncoprotein inhibition in melanoma suppresses tumor growth when ROS defenses are also breached. Cell Rep. 2022, 39, 110995. [Google Scholar] [CrossRef]
  6. Shen, S.; Faouzi, S.; Souquere, S.; Roy, S.; Routier, E.; Libenciuc, C.; André, F.; Pierron, G.; Scoazec, J.Y.; Robert, C. Melanoma persister cells are tolerant to BRAF/MEK Inhibitors via ACOX1-mediated Fatty Acid Oxidation. Cell Rep. 2020, 33, 108421. [Google Scholar] [CrossRef]
  7. Vivas-García, Y.; Falletta, P.; Liebing, J.; Louphrasitthiphol, P.; Feng, Y.; Chauhan, J.; Scott, D.A.; Glodde, N.; Chocarro-Calvo, A.; Bonham, S.; et al. Lineage-Restricted Regulation of SCD and Fatty Acid Saturation by MITF Controls Melanoma Phenotypic Plasticity. Mol. Cell 2020, 77, 120–137. [Google Scholar] [CrossRef]
  8. Garandeau, D.; Mrad, M.; Levade, T.; Perrotta, C.; Andrieu-Abadie, N.; Diab-Assaf, M. Dysregulation of Sphingolipid Metabolism in Melanoma: Roles in Pigmentation, Cell Survival and Tumor Progression. In Bioactive Sphingolipids in Cancer Biology and Therapy; Hannun, Y., Luberto, C., Mao, C., Obeid, L., Eds.; Springer: Cham, Switzerland, 2015; pp. 123–139. [Google Scholar] [CrossRef]
  9. Carrié, L.; Virazels, M.; Dufau, C.; Montfort, A.; Levade, T.; Ségui, B.; Andrieu-Abadie, N. New Insights into the Role of Sphingolipid Metabolism in Melanoma. Cells 2020, 9, 1967. [Google Scholar] [CrossRef] [PubMed]
  10. Vergani, E.; Busico, A.; Dugo, M.; Devecchi, A.; Valeri, B.; Cossa, M.; Di Guardo, L.; De Cecco, L.; Feltrin, E.; Valle, G.; et al. Genetic Layout of Melanoma Lesions Is Associated with BRAF/MEK-Targeted Therapy Resistance and Transcriptional Profiles. J. Investig. Dermatol. 2022, 142, 3030–3040.e5. [Google Scholar] [CrossRef] [PubMed]
  11. Vergani, E.; Beretta, G.L.; Aloisi, M.; Costantino, M.; Corno, C.; Frigerio, S.; Tinelli, S.; Dugo, M.; Accattatis, F.M.; Granata, A.; et al. Targeting of the Lipid Metabolism Impairs Resistance to BRAF Kinase Inhibitor in Melanoma. Front. Cell Dev. Biol. 2022, 10, 927118. [Google Scholar] [CrossRef] [PubMed]
  12. Stamatakos, S.; Beretta, G.L.; Vergani, E.; Dugo, M.; Corno, C.; Corna, E.; Tinelli, S.; Frigerio, S.; Ciusani, E.; Rodolfo, M.; et al. Deregulated FASN Expression in BRAF Inhibitor-Resistant Melanoma Cells Unveils New Targets for Drug Combinations. Cancers 2021, 13, 2284. [Google Scholar] [CrossRef]
  13. Van Meer, G.; Voelker, D.R.; Feigenson, G.W. Membrane lipids: Where they are and how they behave. Nat. Rev. Mol. Cell Biol. 2008, 9, 112–124. [Google Scholar] [CrossRef] [PubMed]
  14. Chiurchiù, V.; Leuti, A.; Maccarrone, M. Bioactive Lipids and Chronic Inflammation: Managing the Fire Within. Front. Immunol. 2018, 9, 38. [Google Scholar] [CrossRef]
  15. Hussain, G.; Anwar, H.; Rasul, A.; Imran, A.; Qasim, M.; Zafar, S.; Imran, M.; Kamran, S.K.S.; Aziz, N.; Razzaq, A.; et al. Lipids as biomarkers of brain disorders. Crit. Rev. Food Sci. Nutr. 2020, 60, 351–374. [Google Scholar] [CrossRef]
  16. van Kruining, D.; Luo, Q.; van Echten-Deckert, G.; Mielke, M.M.; Bowman, A.; Ellis, S.; Oliveira, T.G.; Martinez-Martinez, P. Sphingolipids as prognostic biomarkers of neurodegeneration, neuroinflammation, and psychiatric diseases and their emerging role in lipidomic investigation methods. Adv. Drug Deliv. Rev. 2020, 159, 232–244. [Google Scholar] [CrossRef] [PubMed]
  17. Neagu, M. Metabolic Traits in Cutaneous Melanoma. Front. Oncol. 2020, 10, 851. [Google Scholar] [CrossRef]
  18. Molinelli, E.; Ceccarelli, G.; Fantone, S.; Di Mercurio, E.; Gambini, D.; Maurizi, A.; Perugini, J.; Tossetta, G.; Brisigotti, V.; De Simoni, E.; et al. Melanoma and subcutaneous adipose tissue: Role of peritumoral adipokines in disease characterization and prognosis. Pigment Cell Melanoma Res. 2023, 36, 423–430. [Google Scholar] [CrossRef]
  19. Lazar, I.; Clement, E.; Dauvillier, S.; Milhas, D.; Ducoux-Petit, M.; LeGonidec, S.; Moro, C.; Soldan, V.; Dalle, S.; Balor, S.; et al. Adipocyte Exosomes Promote Melanoma Aggressiveness through Fatty Acid Oxidation: A Novel Mechanism Linking Obesity and Cancer. Cancer Res. 2016, 76, 4051–4057. [Google Scholar] [CrossRef]
  20. Zhang, M.; Di Martino, J.S.; Bowman, R.L.; Campbell, N.R.; Baksh, S.C.; Simon-Vermot, T.; Kim, I.S.; Haldeman, P.; Mondal, C.; Yong-Gonzales, V.; et al. Adipocyte-Derived Lipids Mediate Melanoma Progression via FATP Proteins. Cancer Discov. 2018, 8, 1006–1025. [Google Scholar] [CrossRef]
  21. Alicea, G.M.; Rebecca, V.W.; Goldman, A.R.; Fane, M.E.; Douglass, S.M.; Behera, R.; Webster, M.R.; Kugel, C.H., 3rd; Ecker, B.L.; Caino, M.C.; et al. Changes in Aged Fibroblast Lipid Metabolism Induce Age-Dependent Melanoma Cell Resistance to Targeted Therapy via the Fatty Acid Transporter FATP2. Cancer Discov. 2020, 10, 1282–1295. [Google Scholar] [CrossRef]
  22. Kuzu, O.F.; Noory, M.A.; Robertson, G.P. The Role of Cholesterol in Cancer. Cancer Res. 2016, 76, 2063–2070. [Google Scholar] [CrossRef]
  23. Delgado-Goñi, T.; Galobart, T.C.; Wantuch, S.; Normantaite, D.; Leach, M.O.; Whittaker, S.R.; Beloueche-Babari, M. Increased inflammatory lipid metabolism and anaplerotic mitochondrial activation follow acquired resistance to vemurafenib in BRAF-mutant melanoma cells. Br. J. Cancer 2020, 122, 72–81. [Google Scholar] [CrossRef]
  24. Bataller, M.; Sánchez-García, A.; Garcia-Mayea, Y.; Mir, C.; Rodriguez, I.; LLeonart, M.E. The Role of Sphingolipids Metabolism in Cancer Drug Resistance. Front. Oncol. 2021, 11, 807636. [Google Scholar] [CrossRef]
  25. Agatonovic-Kustrin, S.; Morton, D.W.; Smirnov, V.; Petukhov, A.; Gegechkori, V.; Kuzina, V.; Gorpinchenko, N.; Ramen-skaya, G. Analytical Strategies in Lipidomics for Discovery of Functional Biomarkers from Human Saliva. Dis. Markers 2019, 2019, 6741518. [Google Scholar] [CrossRef]
  26. Kim, H.Y.; Lee, H.; Kim, S.H.; Jin, H.; Bae, J.; Choi, H.K. Discovery of potential biomarkers in human melanoma cells with different metastatic potential by metabolic and lipidomic profiling. Sci. Rep. 2017, 7, 8864. [Google Scholar] [CrossRef]
  27. Sagini, K.; Urbanelli, L.; Buratta, S.; Emiliani, C.; Llorente, A. Lipid Biomarkers in Liquid Biopsies: Novel Opportunities for Cancer Diagnosis. Pharmaceutics 2023, 15, 437. [Google Scholar] [CrossRef]
  28. Hersey, P. Ganglioside antigens in tissue sections of skin, naevi, and melanoma—Implications for treatment of melanoma. Cancer Treat. Res. 1991, 54, 137–151. [Google Scholar] [CrossRef]
  29. Garandeau, D.; Noujarède, J.; Leclerc, J.; Imbert, C.; Garcia, V.; Bats, M.L.; Rambow, F.; Gilhodes, J.; Filleron, T.; Meyer, N.; et al. Targeting the Sphingosine 1-Phosphate Axis Exerts Potent Antitumor Activity in BRAFi-Resistant Melanomas. Mol. Cancer Ther. 2019, 18, 289–300. [Google Scholar] [CrossRef] [PubMed]
  30. Vinceti, M.; Malagoli, C.; Iacuzio, L.; Crespi, C.M.; Sieri, S.; Krogh, V.; Marmiroli, S.; Pellacani, G.; Venturelli, E. Serum Fatty acids and risk of cutaneous melanoma: A population-based case-control study. Dermatol. Res. Pract. 2013, 2013, 659394. [Google Scholar] [CrossRef]
  31. Liyanage, U.E.; Law, M.H.; Ong, J.S.; Cust, A.E.; Mann, G.J.; Ward, S.V.; Melanoma Meta-Analysis Consortium; Gharahkhani, P.; Iles, M.M.; MacGregor, S. Polyunsaturated fatty acids and risk of melanoma: A Mendelian randomisation analysis. Int. J. Cancer 2018, 143, 508–514. [Google Scholar] [CrossRef] [PubMed]
  32. Vilbert, M.; Koch, E.C.; Rose, A.A.N.; Laister, R.C.; Gray, D.; Sotov, V.; Penny, S.; Spreafico, A.; Pinto, D.M.; Butler, M.O.; et al. Analysis of the Circulating Metabolome of Patients with Cutaneous, Mucosal and Uveal Melanoma Reveals Distinct Metabolic Profiles with Implications for Response to Immunotherapy. Cancers 2023, 15, 3708. [Google Scholar] [CrossRef]
  33. Yang, Z.; Zhang, D.; Sima, X.; Fu, Y.; Zeng, H.; Hu, Z.; Hou, J.; Pan, Y.; Zhang, Y.; Zhou, Z.; et al. Levels of pretreatment serum lipids predict responses to PD-1 inhibitor treatment in advanced intrahepatic cholangiocarcinoma. Int. Immunopharmacol. 2023, 115, 109687. [Google Scholar] [CrossRef] [PubMed]
  34. McQuade, J.L.; Daniel, C.R.; Hess, K.R.; Mak, C.; Wang, D.Y.; Rai, R.R.; Park, J.J.; Haydu, L.E.; Spencer, C.; Wongchenko, M.; et al. Association of body-mass index and outcomes in patients with metastatic melanoma treated with targeted therapy, immunotherapy, or chemotherapy: A retrospective, multicohort analysis. Lancet Oncol. 2018, 19, 310–322. [Google Scholar] [CrossRef]
  35. Wang, Y.; Guo, Z.; Isah, A.D.; Chen, S.; Ren, Y.; Cai, H. Lipid metabolism and tumor immunotherapy. Front. Cell Dev. Biol. 2023, 11, 1187989. [Google Scholar] [CrossRef]
  36. Luo, Y.; Wang, H.; Liu, B.; Wei, J. Fatty acid metabolism and cancer immunotherapy. Curr. Oncol. Rep. 2022, 24, 659–670. [Google Scholar] [CrossRef] [PubMed]
  37. Dei Cas, M.; Zulueta, A.; Mingione, A.; Caretti, A.; Ghidoni, R.; Signorelli, P.; Paroni, R. An Innovative Lipidomic Workflow to Investigate the Lipid Profile in a Cystic Fibrosis Cell Line. Cells 2020, 9, 1197. [Google Scholar] [CrossRef] [PubMed]
  38. Morano, C.; Zulueta, A.; Caretti, A.; Roda, G.; Paroni, R.; Dei Cas, M. An Update on Sphingolipidomics: Is Something Still Missing? Some Considerations on the Analysis of Complex Sphingolipids and Free-Sphingoid Bases in Plasma and Red Blood Cells. Metabolites 2022, 12, 450. [Google Scholar] [CrossRef]
  39. Hosmer, D.W., Jr.; Lemeshow, S. Applied Logistic Regression, 2nd ed.; John Wiley & Sons, Inc.: New York, NY, USA, 2000; pp. 1–397. [Google Scholar]
  40. Artusi, R.; Verderio, P.; Marubini, E. Bravais-Pearson and Spearman correlation coefficients: Meaning, test of hypothesis and confidence interval. Int. J. Biol. Markers 2002, 17, 148–151. [Google Scholar] [CrossRef]
  41. Hollander, M.; Wolfe, D.A.; Chicken, E. Nonparametric Statistical Methods, 3rd ed.; John Wiley & Sons, Inc.: New York, NY, USA, 2014; pp. 1–848. [Google Scholar]
  42. Joliffe, I. Principal Component Analysis; Springer: New York, NY, USA, 1986. [Google Scholar]
Figure 1. Distribution of FASN and DHCR24 levels in melanoma patients and control subjects. (A,B) Box plots reflecting the distribution of FASN (ng/mL) and DHCR24 (pg/mL) in melanoma patients and control subjects in (C,D) when stratified in groups according to clinical stages and as DH (day-hospital) or HD (healthy donors) for controls groups. (E,F) Box plots reflecting the overall distribution of FASN and DHCR24 in melanoma patients according to short-term (SR, ≤12 months) or long-term (LR, ≥13 months) response to treatment, and (G,H) stratified by medical treatment (ICI—immune checkpoint inhibitors; KI—kinase inhibitors). Each box indicates the 25th and 75th centiles. The horizontal line inside the box indicates the median and the whiskers indicate the extreme measured values. The reported p-values refer to ANOVA significant results.
Figure 1. Distribution of FASN and DHCR24 levels in melanoma patients and control subjects. (A,B) Box plots reflecting the distribution of FASN (ng/mL) and DHCR24 (pg/mL) in melanoma patients and control subjects in (C,D) when stratified in groups according to clinical stages and as DH (day-hospital) or HD (healthy donors) for controls groups. (E,F) Box plots reflecting the overall distribution of FASN and DHCR24 in melanoma patients according to short-term (SR, ≤12 months) or long-term (LR, ≥13 months) response to treatment, and (G,H) stratified by medical treatment (ICI—immune checkpoint inhibitors; KI—kinase inhibitors). Each box indicates the 25th and 75th centiles. The horizontal line inside the box indicates the median and the whiskers indicate the extreme measured values. The reported p-values refer to ANOVA significant results.
Ijms 25 01558 g001
Figure 2. Untargeted lipidomic analysis of plasma samples from 20 melanoma patients at advanced stages (stage III–IV) and 20 sex-matched control subjects (control). Box plots reflecting the distribution of the differentially expressed lipids between cases and controls (A) at the single lipid level and (B) at the lipid family level. (C) Box plots reflecting the distribution of the differentially expressed lipids between LR and SR patients (short-term SR, ≤12 months, or long-term LR, ≥13 months). Each box indicates the 25th and 75th centiles. The horizontal line inside the box indicates the median and the whiskers indicate the extreme measured values expressed as µM. The reported p-values refer to FDR-adjusted Wilcoxon significant results.
Figure 2. Untargeted lipidomic analysis of plasma samples from 20 melanoma patients at advanced stages (stage III–IV) and 20 sex-matched control subjects (control). Box plots reflecting the distribution of the differentially expressed lipids between cases and controls (A) at the single lipid level and (B) at the lipid family level. (C) Box plots reflecting the distribution of the differentially expressed lipids between LR and SR patients (short-term SR, ≤12 months, or long-term LR, ≥13 months). Each box indicates the 25th and 75th centiles. The horizontal line inside the box indicates the median and the whiskers indicate the extreme measured values expressed as µM. The reported p-values refer to FDR-adjusted Wilcoxon significant results.
Ijms 25 01558 g002
Figure 3. Targeted lipidomic analysis of complex sphingolipids levels in plasma samples from 20 melanoma patients at advanced stages (stage III–IV) and 20 sex-matched control subjects (control). Box plots reflecting the distribution of plasma sphingolipid profiles at lipid family level (A) in cases and controls, and (B) in patients according to treatment response (short-term SR, ≤12 months, or long-term LR, ≥13 months). Each box indicates the 25th and 75th centiles. The horizontal line inside the box indicates the median and the whiskers indicate the extreme measured values, expressed as µM. The reported p-values refer to FDR-adjusted Wilcoxon significant results.
Figure 3. Targeted lipidomic analysis of complex sphingolipids levels in plasma samples from 20 melanoma patients at advanced stages (stage III–IV) and 20 sex-matched control subjects (control). Box plots reflecting the distribution of plasma sphingolipid profiles at lipid family level (A) in cases and controls, and (B) in patients according to treatment response (short-term SR, ≤12 months, or long-term LR, ≥13 months). Each box indicates the 25th and 75th centiles. The horizontal line inside the box indicates the median and the whiskers indicate the extreme measured values, expressed as µM. The reported p-values refer to FDR-adjusted Wilcoxon significant results.
Ijms 25 01558 g003
Figure 4. Targeted lipidomic analysis of sphingoid bases and their phosphate form and of saturated and unsaturated free fatty acids between C16 and C18 in plasma samples from 20 melanoma patients at advanced stages (stage III–IV) and 20 sex-matched control subjects (control). Box plots showing (A) the distribution of plasma sphingoid-bases profiles in cases and controls, (B) of levels of free FA in plasma of cases and controls, and (C) in patients according to treatment response (short-term SR, ≤12 months, or long-term LR, ≥13 months). Each box indicates the 25th and 75th centiles. The horizontal line inside the box indicates the median and the whiskers indicate the extreme measured values, expressed as µM. The reported p-values refer to FDR-adjusted Wilcoxon significant results.
Figure 4. Targeted lipidomic analysis of sphingoid bases and their phosphate form and of saturated and unsaturated free fatty acids between C16 and C18 in plasma samples from 20 melanoma patients at advanced stages (stage III–IV) and 20 sex-matched control subjects (control). Box plots showing (A) the distribution of plasma sphingoid-bases profiles in cases and controls, (B) of levels of free FA in plasma of cases and controls, and (C) in patients according to treatment response (short-term SR, ≤12 months, or long-term LR, ≥13 months). Each box indicates the 25th and 75th centiles. The horizontal line inside the box indicates the median and the whiskers indicate the extreme measured values, expressed as µM. The reported p-values refer to FDR-adjusted Wilcoxon significant results.
Ijms 25 01558 g004
Table 1. Study sample characteristics.
Table 1. Study sample characteristics.
FrequencyPercent
Sex
 F10741.31
 M15258.69
Patient’s Group
 PT DH I/II3613.9
 PT III2810.81
 PT III–IV8030.89
 control DH7528.96
 control HD4015.44
Treatment type (PT III–IV only)
 ICI (IPI/NIVO)4050
 KI (BRAF/MEK)4050
Response class (PT III–IV only)
 SR5568.75
 LR2531.25
Age
n; median (range)
256; 50 (17–85)
Triglycerides (mg/dL)
n; median (range)
258; 114 (40–528)
Total cholesterol (mg/dL)
n; median (range)
258; 194 (69–336)
LDL (mg/dL)
n; median (range)
258; 111.6 (5.8–242.8)
HDL (mg/dL)
n; median (range)
258; 53 (12–227)
DHCR24 (pg/mL)
n; median (range)
259; 24.78 (3.61–687.82)
FASN (ng/mL)
n; median (range)
258; 0.74 (0.08–3.97)
DH—day-hospital; HD—healthy donors; ICI—immune checkpoint inhibitors; KI—kinase inhibitors; SR—short-term responders (≤12 months); LR—long-term responders (≥13 months).
Table 2. Plasma levels of FASN, DHCR24, TG, CHOL, LDL, and HDL: results from univariate logistic regression models according to disease status.
Table 2. Plasma levels of FASN, DHCR24, TG, CHOL, LDL, and HDL: results from univariate logistic regression models according to disease status.
OR (95% CI) *OR (95% CI) **
DHCR24 (pg/mL)0.74 (0.56; 0.96)0.84 (0.62; 1.13)
FASN (ng/mL)1.43 (1.09; 1.88)1.24 (0.92;1.67)
TG (mg/dL)1.26 (0.97; 1.64)1.04 (0.79; 1.37)
Total CHOL (mg/dL)0.97 (0.76; 1.24)0.77 (0.58; 1.03)
LDL (mg/dL)0.99 (0.78; 1.27)0.83 (0.62; 1.10)
HDL (mg/dL)0.77 (0.58; 1.03)0.78 (0.56; 1.07)
* Odds ratio computed for a specific unit change for each variable equal to one SD; ** Odds ratio computed for a specific unit change for each variable equal to one SD and adjusted by age.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Dei Cas, M.; Ciniselli, C.M.; Vergani, E.; Ciusani, E.; Aloisi, M.; Duroni, V.; Verderio, P.; Ghidoni, R.; Paroni, R.; Perego, P.; et al. Alterations in Plasma Lipid Profiles Associated with Melanoma and Therapy Resistance. Int. J. Mol. Sci. 2024, 25, 1558. https://doi.org/10.3390/ijms25031558

AMA Style

Dei Cas M, Ciniselli CM, Vergani E, Ciusani E, Aloisi M, Duroni V, Verderio P, Ghidoni R, Paroni R, Perego P, et al. Alterations in Plasma Lipid Profiles Associated with Melanoma and Therapy Resistance. International Journal of Molecular Sciences. 2024; 25(3):1558. https://doi.org/10.3390/ijms25031558

Chicago/Turabian Style

Dei Cas, Michele, Chiara Maura Ciniselli, Elisabetta Vergani, Emilio Ciusani, Mariachiara Aloisi, Valeria Duroni, Paolo Verderio, Riccardo Ghidoni, Rita Paroni, Paola Perego, and et al. 2024. "Alterations in Plasma Lipid Profiles Associated with Melanoma and Therapy Resistance" International Journal of Molecular Sciences 25, no. 3: 1558. https://doi.org/10.3390/ijms25031558

APA Style

Dei Cas, M., Ciniselli, C. M., Vergani, E., Ciusani, E., Aloisi, M., Duroni, V., Verderio, P., Ghidoni, R., Paroni, R., Perego, P., Beretta, G. L., Gatti, L., & Rodolfo, M. (2024). Alterations in Plasma Lipid Profiles Associated with Melanoma and Therapy Resistance. International Journal of Molecular Sciences, 25(3), 1558. https://doi.org/10.3390/ijms25031558

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop