Skip to main content
Download PDF
Download ePub

Lipidome characterisation and sex-specific differences in type 1 and type 2 diabetes mellitus

Cardiovascular Diabetology volume 23, Article number: 109 (2024) Cite this article

Abstract

Background

In this study, we evaluated the lipidome alterations caused by type 1 diabetes (T1D) and type 2 diabetes (T2D), by determining lipids significantly associated with diabetes overall and in both sexes, and lipids associated with the glycaemic state.

Methods

An untargeted lipidomic analysis was performed to measure the lipid profiles of 360 subjects (91 T1D, 91 T2D, 74 with prediabetes and 104 controls (CT)) without cardiovascular and/or chronic kidney disease. Ultra-high performance liquid chromatography-electrospray ionization mass spectrometry (UHPLC-ESI-MS) was conducted in two ion modes (positive and negative). We used multiple linear regression models to (1) assess the association between each lipid feature and each condition, (2) determine sex-specific differences related to diabetes, and (3) identify lipids associated with the glycaemic state by considering the prediabetes stage. The models were adjusted by sex, age, hypertension, dyslipidaemia, body mass index, glucose, smoking, systolic blood pressure, triglycerides, HDL cholesterol, LDL cholesterol, alternate Mediterranean diet score (aMED) and estimated glomerular filtration rate (eGFR); diabetes duration and glycated haemoglobin (HbA1c) were also included in the comparison between T1D and T2D.

Results

A total of 54 unique lipid subspecies from 15 unique lipid classes were annotated. Lysophosphatidylcholines (LPC) and ceramides (Cer) showed opposite effects in subjects with T1D and subjects with T2D, LPCs being mainly up-regulated in T1D and down-regulated in T2D, and Cer being up-regulated in T2D and down-regulated in T1D. Also, Phosphatidylcholines were clearly down-regulated in subjects with T1D. Regarding sex-specific differences, ceramides and phosphatidylcholines exhibited important diabetes-associated differences due to sex. Concerning the glycaemic state, we found a gradual increase of a panel of 1-deoxyceramides from normoglycemia to prediabetes to T2D.

Conclusions

Our findings revealed an extensive disruption of lipid metabolism in both T1D and T2D. Additionally, we found sex-specific lipidome changes associated with diabetes, and lipids associated with the glycaemic state that can be linked to previously described molecular mechanisms in diabetes.

Introduction

Diabetes mellitus (DM) is characterized by chronic hyperglycaemia that leads to heterogenous disturbances of metabolism [1] and its continuing rise is a major concern in society [2]. DM is one of the main risk factors for cardiovascular disease and other conditions [3]. Therefore, a better understanding of diabetes pathophysiology has become a subject of major interest in research.

Lipid disruption has been associated with many human diseases, leading to the rise in the relevance of lipidomics, an emerging field involving the study of lipids and factors that interact with lipids [4]. The association between lipids and DM has been widely recognised. Advanced lipoprotein analyses have shown a reduction in serum concentrations of triglycerides (TG), cholesterol, and apolipoprotein (Apo)B-containing lipoproteins when comparing subjects at the onset of T1D and after achieving optimal glycaemic control [5]. In addition, epidemiological studies have shown a close relationship between low-density lipoprotein cholesterol (LDL-cholesterol) and high-density lipoprotein cholesterol (HDL-cholesterol) concentrations in T2D. However, the complexity of the associations between diabetes and lipid metabolites is underestimated in these studies since lipoproteins contain a great variety of lipids that remain unanalysed, highlighting the crucial role lipidomics can play [6].

The association between lipid species and T1D still needs to be fully understood. Most of the studies have focused on biomarker discovery for T1D risk during childhood [7]. One study reported several lipid species significantly altered in subjects at the onset of T1D and after achieving glycaemic control [8]. On the other hand, the number of studies focused on lipidomic changes associated with T2D risk is higher [9,10,11,12]. Lipidome differences between normoglycemic, prediabetic and T2D subjects have also been described [13], as well as lipid species significantly associated with T2D complications, such as diabetic retinopathy [14], diabetic neuropathy [15] or diabetic nephropathy [16].

Nevertheless, to the best of our knowledge, there is a lack of studies comparing the lipidome of subjects with T1D and T2D. Moreover, some studies have shown evidence that the risk of diabetes complications differs between the sexes [17]; however, the underlying mechanisms behind these sex-specific differences are poorly understood.

In the present study, we extensively investigated the serum lipidome of T1D, T2D, and non-diabetic subjects through an untargeted lipidomics analysis using Ultra High-Performance Liquid Chromatography-Mass Spectrometry (UHPLC-MS). Our objectives were to identify lipid subspecies that differ between (1) subjects with T1D and with T2D, (2) subjects with T1D and non-diabetic subjects, (3) subjects with T2D and non-diabetic subjects, (4) to determine sex-specific differences in each of the above-mentioned comparisons, and (5) compare the lipidome between normoglycemia, prediabetes and T2D.

Methods

Participants

In this study, 536 participants including 156 with T1D, 159 with T2D, and 221 without diabetes and matched by sex and BMI, were selected from previous cohorts, at the University Hospitals Arnau de Vilanova (Lleida, Spain), Germans Trias i Pujol (Badalona, Spain), Clinic (Barcelona, Spain), and the Primary Care Center Mollerussa (Lleida, Spain) [18,19,20,21] (Additional File1 - Figure S1). The inclusion criteria for all groups were: aged between 20 and 85 years, the absence of established chronic kidney disease (defined as calculated glomerular filtration rate < 60 mL/min and/or urine albumin/creatinine ratio > 299 mg/g), and absence of known clinical cardiovascular events or associated revascularization procedures, including coronary heart disease, cerebrovascular disease, or peripheral vascular disease (including the diagnosis of diabetic foot disease).

Age, sex, tobacco exposure and pharmacological treatment were recorded. Diabetes duration was acquired from the medical records. Subjects were considered to have hypertension or dyslipidaemia if they were under anti-hypertensive or lipid-lowering treatment, respectively. Anthropometric data, weight, height, waist circumference, BMI, and blood pressure were obtained using standard methods. The standard biochemical analysis included glucose and glycated hemoglobin (HbA1c), lipid profile, and estimated glomerular filtration rate calculated according to the Chronic Kidney Disease Epidemiology Collaboration (CKD-EPI) equation [22]. The dietary pattern was assessed using the alternate Mediterranean Diet score (aMED), as described previously [23]. This score includes monounsaturated-to-saturated fat ratio, legumes, vegetables, nuts, fruits, nuts, cereals, fish, meat and wine. This score ranges from 0 to 9, with higher scores indicative of a higher adherence to the MedDiet.

Blood samples were collected in the fasting state and blood tests were conducted using standard laboratory methods [18]. Urine tests were performed in subjects with diabetes following standard laboratory methods. Subjects with normoglycemia and prediabetes were classified using the American Diabetes Association criteria (glycated haemoglobin (HbA1c) < 5.7% or fasting plasma glucose ≤ 100 mg/dL for normoglycemia; and HbA1c between 5.7% and < 6.5% or fasting plasma glucose between 101 mg/dL and < 126 mg/dL for prediabetes) [24].

Blood samples for the lipidomic analyses were collected in the fasting state with EDTA tubes, processed immediately after extraction, and stored at − 80 °C at the biobanks of the participant centres until determination.

From the 536 samples, 23 were discarded due to technical problems, and for the present study, the lipid profiles of 360 participants, 91 with T1D, 91 with T2D and 178 without diabetes were selected from the full cohort (Additional File 1 - Figure S1). Serum samples for all the selected participants were collected by University Hospital Arnau of Vilanova (Lleida, Spain) and the Mollerussa Primary Health care area (Lleida, Spain), as part of previously published studies [18,19,20]. The acquisition, processing and storage of these samples were performed in the same facility, and the geographic location of the subjects was the same, thus avoiding any variability due to sample origin that might be present in the full cohort (Additional File 1 - Figure S1).

Sample preparation

Due to the large number of samples in the full cohort, samples were randomly assigned to one of 6 batches. To reduce the impact of technical factors, the sample order within each batch was randomized before sample preparation, and then again prior to measurement of the lipid profile by UHPLC-ESI-MS/MS. All serum samples were defrosted on ice, and each sample was aliquoted (50 µL) to create a pooled quality control (QC) representative of all samples in the study. The pooled QC was vortexed, further aliquoted (50 µL), and stored at -80 °C until the analysis of each of the 6 batches of QC samples. Lipid extraction was performed by mixing 50 µL of biological sample or QC with 150 µL isopropanol (LC-MS grade), vortexed for 20 s, and centrifuged at 22,000 g for 20 min at 4 °C. 120 µL of the supernatant was transferred to a low recovery vial and transferred to the LC sample manager at 4 °C.

Ultra-high-performance liquid chromatography-mass spectrometry

Samples were maintained at 4 °C and analysed by applying UHPLC-MS methods using a Dionex UltiMate 3000 Rapid Separation LC system (Thermo Fisher Scientific, MA, USA) coupled with a heated electrospray Q Exactive Focus mass spectrometer (Thermo Fisher Scientific, MA, USA). Non-polar extracts were analysed on a Hypersil GOLD column (100 × 2.1 mm, 1.9 μm; Thermo Fisher Scientific, MA, USA). Mobile phase A consisted of 10 mM ammonium formate and 0.1% formic acid in 60% acetonitrile/water and mobile phase B consisted of 10 mM ammonium formate and 0.1% formic acid in 90% propan-2-ol/water. Flow rate was set for 0.40 mL/min with the following gradient: t = 0.0, 20% B; t = 0.5, 20% B, t = 8.5, 100% B; t = 9.5, 100% B; t = 11.5, 20% B; t = 14.0, 20% B, all changes were linear with curve = 5. The column temperature was set to 55 °C and the injection volume was 2µL. Data were acquired in positive and negative ionization mode separately within the mass range of 150–2000 m/z at resolution 70,000 (FWHM at m/z 200). Ion source parameters were: sheath gas = 50 arbitrary units, Aux gas = 13 arbitrary units, sweep gas 3 arbitrary units, spray voltage 3.5 kV (positive ion mode) and 3.1 kV (negative ion mode), Capillary temp = 263 °C, and Aux gas heater = 425 °C. Data dependent MS2 in ‘Discover mode’ was applied for the MS/MS spectral acquisition with the following settings: resolution at 17,500 (FWHM at m/z 200), isolation width 3.0 m/z, stepped normalised collision energy at 20, 50 and 80%. Spectra were acquired at three mass ranges 200–400 m/z, 400–700 m/z and 700–1500 m/z on the pooled QC samples. Thermo ExactiveTune (2.8 SP1 build 2806) software was used to control the instrument in both cases, with data acquired in profile mode. Quality control samples were acquired in both profile and dependent scan mode at the start of the run (i.e., 7 QCs MS1 only, 3 QCs with MS2) and then every seventh injection with two QC samples at the end of the analytical batch. Preparation blank samples were analysed between QCs 5 and 6 and at the end of the analytical batch.

Mass spectrometry raw data processing

Raw data acquired in each analytical batch were converted from the instrument-specific format to a mzML file format using the open access ProteoWizard (version 3.0.11417) msconvert tool [25]. Deconvolution was performed by the R package XCMS (version 1.46.0, running in the Galaxy workflow environment) [26]. Isotopologue Parameter Optimization (IPO - version 1.0.0) [27] was used to optimise the XCMS peak picking parameters. A data matrix of metabolite features (m/z-retention time pairs) versus samples was constructed with peak areas provided.

Assessment of data quality and peak matrix filtering

The first five QCs for each batch were used to equilibrate the analytical system and therefore subsequently removed before the data was processed and analysed. Data matrices were corrected for run-order drift in intensity for each lipid feature separately using the Quality Control-Robust Spline Correction (QC-RSC) algorithm [28] in the R environment using the pmp package [29]. Principal Component Analysis (PCA) was used to identify and remove (PCs 1 and 2, Hotelling T2p < 0.05) suspected outlier (QC) samples within each batch to ensure robust correction. Blank samples at the start and end of a run were used to remove features from non-biological origins. Any feature with an average QC intensity less than 20 times the average intensity of the blanks was removed. Any sample with > 50% missing values was excluded from further analysis. Metabolite features with RSD > 30% and present in less than 90% of the QC samples were deleted from the dataset. Features with a < 50% detection rate over all samples were also removed. All data preparation steps were undertaken in R using the structToolbox package [30, 31].

Statistical analysis

After removing the observations with missing values in the variables included in the models (Additional File 1 - Figure S1), the clinical data of participants was summarised as mean (standard deviation) for continuous variables and as frequency (percentage) for categorical data, using the compareGroups R package [32].

Analysis of the UHPLC-ESI-MS/MS data was conducted in the R environment [33]. Prior to the statistical analysis, Probabilistic Quotient Normalization (PQN) [34], using the mean of the QC samples as a reference, was applied. Data were log-transformed to reduce skewness.

Multiple linear regression models were used to assess the association between each metabolite and T1D, T2D, and non-diabetic controls (CT). Three different comparisons were performed: T1D against T2D, T1D against CT, and T2D against CT. The models were adjusted by sex, age, hypertension, dyslipidaemia, body mass index, glucose, smoking, systolic blood pressure, TG, HDL cholesterol, LDL cholesterol, alternate Mediterranean diet score and estimated glomerular filtration rate (eGFR); diabetes duration and glycated haemoglobin (HbA1c) were also included in the T1D and T2D comparison.

Regarding the sex-specific differences related to diabetes, the same linear models were used, but an interaction term between diabetes and sex (basal level: men) was added. Using this configuration, the p-value and regressor (βMen) associated with the diabetes variable were assigned to men. The effect of diabetes in women (βWomen) was computed by summing the regressor of the diabetes variable and the interaction between diabetes and sex. The standard error of the effect was computed, the t-value was obtained by dividing βWomen by its standard error, and the p-value associated with diabetes in women was computed using this t-value.

Concerning the lipid alterations associated with the glycaemic state, we considered the following categories: normoglycemia, prediabetes and T2D. The models were adjusted using the confounding factors mentioned above, and a numeric variable defining the glycaemic state: 0, normoglycemia; 1, prediabetes; 2, T2D.

In all analyses, False Discovery Rate (FDR) correction was performed, and a corrected p-value of < 0.05 was considered significant. For each comparison, we present significant p-values that correspond to lipids when not considering the interaction effect between diabetes and sex (all subjects), p-values for men, p-values for women, and p-values for the glycaemic state.

A description of each analysis is shown in Additional File 1 (Table S1). LipidSearch was used to annotate lipid species. Annotations with grades A or B were mapped to XCMS-detected features based on an absolute ppm error less than 5 and an absolute retention time tolerance of less than 5 s.

Results

Clinical and biological parameters

The baseline characteristics for each group according to diabetes status are shown in Table 1. Subjects with T1D had longer diabetes duration, and higher HDL-cholesterol in comparison with subjects with T2D. On the other hand, subjects with T2D were older, had a higher BMI, higher HbA1c, and higher frequency of hypertension and dyslipidaemia. In Additional File 1 (Table S2), the baseline characteristics for each comparison are shown.

Table 1 Baseline characteristics of subjects with type 1, type 2 diabetes, prediabetes and normoglycemia
Full size table

Data are mean (SD) for continuous variables and number (%) for categorical variables. For continuous variables, the p-values are obtained using a student’s t-test and for categorical variables, a chi-squared test. BMI, body mass index; DM, diabetes mellitus; HbA1c, glycated haemoglobin; sBP, systolic blood pressure; dBP, diastolic blood pressure; TG, triglycerides; HDL, high-density lipoprotein; LDL, low-density lipoprotein; aMED, alternate Mediterranean diet score; eGFR, estimated glomerular filtration rate.

Lipid markers of diabetes mellitus

Table 2 shows the overall number of lipidomic features that were significantly different in the comparisons (corrected p-value lower than 0.05), and of those, the number that were annotated using LipidSearch. When comparing T1D and T2D subjects, 30 lipid species (positive and negative ion modes combined) obtained a corrected p-value lower than 0.05 both in males and females (all subjects), 16 lipids were found in females and 17 in males. In the comparison between T1D and CT, 35 lipids were statistically significant in all subjects, 11 lipids in women only, and 16 in men only. Finally, when comparing T2D vs. CT, 15 lipids were significantly different in all subjects, 10 in women only, and 10 in men only. Overall, 54 unique lipid species from 15 classes were determined as significant features across all comparisons.

UpSet plots depicting the number of unique and shared (i.e., intersections) significant lipidomic features in the different comparisons in positive and negative acquisition modes are shown in Additional File 1 (Figures S2 and S3), respectively.

Table 2 Summary of the number of significant lipidomic features in each comparison
Full size table
Fig. 1
figure 1

Figure 1 shows the lipid classes that differ in the comparisons according to diabetes status. When comparing both diabetic conditions, lysophosphatidylcholines (LPC) and ceramides (Cer) were more importantly altered than other lipid classes (Fig. 1A). In a similar way, LPCs, phosphatidylcholines (PC), phosphatidylethanolamines (PE) and TGs were especially altered in T1D (Fig. 1B), as well as Cer in T2D (Fig. 1C). Additional File 2 (Table S3) reports the mass-to-charge ratio (mz) and retention time (rt) for each lipid ion significantly associated with one of these conditions in at least one of the analyses. Moreover, the range of corrected p-values, ionization mode and the list of analyses corresponding to the significant corrected p-value is also shown in Additional File 2 (Table S3)

Full size image

Figure 1. Manhattan plots of the minus logarithm of the corrected p-values (y axis) for each of the lipid classes (x axis) obtained in the analysis of A) T1D vs. T2D, B) T1D vs. CT and C) T2D vs. CT. Corrected p-values are shown for the features that had been annotated using LipidSearch and fulfilled the quality criteria described in the Methods section. The dashed line indicates the threshold of significance (\(0.05\)). AcCa, acylcarnitine; Cer, ceramide; ChE, cholesterol esther; Co, coenzyme; DG, diacylglycerol, dMePE, dimethylphosphatidylethanolamine; Hex1Cer, hexosylceramides; Hex2Cer, dihexosylceramides; Hex3Cer, trihexosylceramides; LdMePE, lysodimethylphosphatidylethanolamine; LPA, Lysophosphatidic acid; LPC, lysophosphatidylcholine; LPE, lysophosphatidylethanolamine; PC, phosphatidylcholine; PE, phosphatidylethanolamine; PI, phosphatidylinositol; SM, sphingomyelin; ST, sterol; StE, Stigmasteryl ester; TG, triacylglycerol; ZyE, Zymosteryl ester

In Figure 2, the fold-changes and significance level of the lipids significantly altered in at least one of the analyses are shown.

Fig. 2
figure 2

Fold-change values and statistical significance obtained for each lipid determined as being significantly different in at least one of the nine analyses conducted. Statistical significance is indicated using asterisks: corrected p-value (p) < 0.05 (*), p < 0.01 (**), p < 0.001 (***), p < 0.0001 (****). Bars in the right side of each panel indicate a positive fold-change value, while bars in the left indicate a negative one. A positive fold-change indicates that the lipid is increased in the first group (e.g. in T1D vs. T2D, TG(18:1_18:1_18:2) is significantly increased in T1D with respect to T2D). The different colours in the background of the plot show the different lipid classes. In the left, the names of the lipid subspecies are shown. The nomenclatures of the type 16:1e indicate that the fatty acid (FA) of the glycerophospholipid is linked to the glycerol moiety by an ether bond, therefore, the mentioned glycerophospholipid is an ether-glycerophospholipid

Full size image

In general, TGs, sphingomyelins (SMs), PCs, diacylglycerols (DGs) and ceramide lipids were down-regulated in subjects with T1D, while phosphatidylethanolamines (PE) lysophosphatidylethanolamines (LPEs) and LPCs were mainly up-regulated when compared to CT. On the other hand, some LPCs were down-regulated in T2D subjects and ceramide lipids were mostly up-regulated with respect to controls. Results comparing T1D and T2D support the observed opposing effects seen in comparisons against the control group. A generalized opposite disruption of ceramide lipids and LPCs in subjects with T1D and subjects with T2D could be detected. Contrary to LPCs, ceramide lipids were mainly upregulated in T2D with respect to T1D. Specifically, the ceramide lipids increased in T2D were 1-deoxyceramides (i.e., Cer(m18:1_22:0), Cer(m18:0_22:0), Cer(m18:1_23:0), Cer(m18:0_23:0), Cer(m18:0_24:1) and Cer(m18:0_24:0)) (Fig. 2).

Several sex-specific lipidomic differences were detected, as shown in Fig. 2. These alterations are further illustrated in Figs. 3 and 4 for T2D and controls comparison and T1D and controls comparison, respectively.

Fig. 3
figure 3

Boxplots of lipids that are significantly associated with T2D in men or women. The nomenclatures of the type 16:1e indicate that the fatty acid (FA) of the glycerophospholipid is linked to the glycerol moiety by an ether bond, therefore, the mentioned glycerophospholipid is an ether-glycerophospholipid

Full size image
Fig. 4
figure 4

Boxplots of 20 selected lipids significantly associated with T1D in men or women

Full size image

In general, our results showed that PC levels were higher in normoglycemic women than their male counterparts.

Figure 5 shows lipids significantly associated with a numeric variable that describes glycaemic stag, by considering the prediabetes stage. Although Co(Q10) was not statistically significant (q-value = 0.07), the boxplot shown in Fig. 5 shows a gradual decrease through the glycaemic progression.

Fig. 5
figure 5

Boxplots of lipids significantly associated with the numeric variable defined as 0 for normoglycemia, 1 for prediabetes and 2 for T2D

Full size image

Our results show a consistent alteration of ceramides, revealing a gradual increase of these lipids in the stage of prediabetes to T2D.

Discussion

Through an untargeted lipidomic serum profiling approach, we investigated the lipidomic alterations in type 1 and type 2 diabetes; the sex-specific differences in these diseases and the glycaemic state by considering the prediabetes stage. The results reported in the present study have revealed 54 lipid species belonging to 15 different lipid classes potentially implicated in well-known mechanisms involved in type 1 and type 2 diabetes.

Lysophosphatidylcholine acyltransferase (LCAT) activity

Our results have revealed a panel of LPCs significantly increased in T1D, specifically LPC(26:0), LPC(20:1), LPC(18:1) and LPC(18:2). These findings are in agreement with other studies, where LPC(18:3) has been consistently found to be positively associated with T1D or in those at risk [7]. The mechanisms behind this increase are unclear; however, SMs have been described as inhibitors of LCAT, enzyme responsible of LPCs and cholesteryl esters synthesis from PCs [35]. In our study, SM(d18:2_22:2) was significantly reduced in T1D subjects, which could be a potential explanation for an enhanced activity of LCAT and a subsequent increase of LPC levels. Furthermore, we observed a general decrease of a set of PC species in subjects with T1D. This aligns with some studies that showed a negative association between PC levels and T1D risk [36, 37], and concurs with the suspected enhancement of LCAT activity. The potential consequences of these changes should also be considered. Serum levels of LPC(18:1) and LPC(18:2) were found to be significantly higher in children with T1D after a diabetes ketoacidosis or a hypoglycaemic episode [38], suggesting an association between certain diabetes complications and LPCs.

Although non-significant, our results showed a reduction of LPC levels in T2D compared to controls, and a significant reduction of LPC(20:1) and LPC(18:2) when compared to subjects with T1D. Moreover, PC(18:0_22:5) and PC(18:0_20:3) were significantly increased in subjects with T2D versus controls. Our results are aligned with other studies that showed a positive association between PCs and T2D risk [39], and suggested a reduction of LCAT activity in T2D [40, 41] and a negative correlation between LCAT activity and HbA1c [40,41,42]. A potential mechanism behind this relationship is the non-enzymatic glycation of apolipoprotein A-I in HDL particles due to hyperglycemia [40, 42].

Polyunsaturated fatty acids (PUFAs)

Our findings also showed a general down-regulation of TG and DG containing PUFAs in subjects with T1D when compared to controls and subjects with T2D. The decrease of oleic (16:1) and palmitoleic (18:1) acids, as well as a reduction in the synthesis of omega-6 and omega-3 PUFAs in T1D has been previously described. This reduction is caused by decreased activity of Δ6 desaturase and the stearoyl-CoA desaturases (SCD). mRNA transcription of these desaturases is activated by the protein SREBP-1c, the expression and activation of which is modulated by insulin and therefore, reduced in T1D [43]. This mechanism could be an explanation for the significant reduction of PUFA-containing TGs and DGs in subjects with T1D. This becomes even more feasible if we consider that all the TGs and DGs found to be significantly reduced in our study, contain at least one of the previously mentioned fatty acids (16:1, 18:2, 18:3 and 20:5).

Ether lipids

Peroxisomal defects causing decreased levels of serum ether lipids have been associated with neurodegenerative diseases, cancer, obesity, hypertension [44] and T2D [45]. Plasmalogen PEs containing PUFA have been associated with lower risk of T2D [46]. However, in the present study, PE(20:0p_18:2) was significantly increased in subjects with T1D and T2D when compared to controls, and PE(16:0p_18:1) and PE(16:1e_18:1) were significantly increased in subjects with T1D. On the other hand, previous studies have revealed reduced levels of ether PCs in lean subjects with T2D [45] and in HDL particles of subjects with T2D [9], as well as inverse associations between ether PCs and T2D risk [39]. Moreover, ether PCs have been positively associated with longevity [47]. Although our results do not show a significant association between ether PCs and diabetes, we have found lower levels of ether LPCs (LPC(18:1e) and LPC(18:2e)) in subjects with T2D, supporting the previously hypothesized peroxisomal defect associated with T2D.

Phosphatidylethanolamines and Lysophosphatidylethanolamines

Regarding phosphatidylethanolamines (PEs), we showed an increase of PE(18:0_20:3) and PE(16:0_18:2) in subjects with T2D. Previous work has reported increased serum levels of PEs in T2D obese subjects. The mechanism proposed to explain this alteration was a higher abundance of PEs in VLDL particles compared to HDL particles, and a relative increase of VLDL particles in obese T2D subjects [45]. Furthermore, another study showed enrichment of PE(38:5), PE(38:6), and PE(40:7) in HDL particles in patients with T2D and CHD [9]. This could explain the increase in PEs observed in the present study, since the lipid species observed in lipidomics mainly stem from circulant lipoproteins. Additionally, PEs have been described as modulators of inflammation and apoptosis [47]. In line with this, our results revealed two PE species, PE(18:0_18:1) and PE(18:1_18:2), significantly increased in subjects with T1D, and interestingly one of their metabolic products LPE(18:1) also increased in T1D. On the other hand, PE(18:0_22:6) was significantly decreased in subjects with T1D, as well as one of its plausible products, LPE(22:6).

Sphingolipids

The direct inhibition of the insulin-signalling pathway caused by sphingolipids has been widely described. Ceramides accumulation interferes in the insulin-stimulated activation of protein kinase B (Akt/PKB), which decreases glucose uptake in skeletal muscles and activates gluconeogenesis and glycogenolysis in the liver. On the other hand, several studies support the theory that elevated intracellular levels of sphingolipids may hinder mitochondrial respiratory chain activity, thus causing alterations in mitochondrial metabolism [48]. Our results revealed an increase of a set of 1-deoxyceramides, Cer(m18:1_22:0), Cer(m18:0_22:0), Cer(m18:1_23:0), Cer(m18:0_23:0), Cer(m18:0_24:1), Cer(m18:0_24:0), in subjects with T2D compared to subjects with T1D and controls. Moreover, this significant association is maintained when comparing normoglycemia, prediabetes and T2D, showing a gradual increase with the glycaemic state. Related to this, it is quite relevant that we found coenzyme Q10, Co(Q10), to be significantly decreased in subjects with T2D compared to controls. Co(Q10) has a key role in the electron transport chain of the mitochondria and its deficiency in subjects with T2D has been previously described [49]. This concurs with our results and is probably related to the above-mentioned mitochondrial dysfunction. Even though this molecule has not been found significantly associated in the prediabetes analysis (q-value = 0.07), the corrected p-value is close to significance. The plot of its progression has been added to the main manuscript, and it is possible to see a non-linear gradual decrease from normoglycemia to T2D.

Sex-specific metabolic changes in type 1 and type 2 diabetes

It has been shown that women have a steeper age-related increase of ceramide levels [50]. The loss of oestrogens during and after menopause has been proposed as the main mechanism behind this pattern [50, 51], but other processes have been proposed, such as the differences in sex steroids or the higher levels of oxidative stress in post-menopausal women [51]. Moreover, pre-menopausal women have better cardiovascular health and CVD outcomes than men, but this tendency changes during and after menopause. Ceramides might have a key role in this process, due to the strong relationship between oestrogens and sphingolipid metabolism and the association of ceramides with apoptosis, oxidative stress, inflammation and endothelial dysfunction [52]. Further, menopause has been associated with an increased risk of T2D [53]. In our study, the age of the female subgroup with T2D was 57.3 (SD: 10.6) years. Therefore, we might assume that in a large proportion of the subjects, menopause was playing a role in the lipidomic differences observed. We found one 1-deoxyceramide, Cer(m18:1_20:0), significantly associated with T2D in women but not in men, and in general, the fold-changes and the significance level of the significant 1-deoxyceramides in T2D were higher in women than in men (Fig. 2). Diabetes has been shown to attenuate the protective effect of the female sex in the development of cardiac diseases and nephropathy [54]. The specific lipids that differ between sexes found in the present study could explain the greater impact of T2D complications in post-menopausal women.

Our results also revealed a greater T1D-associated alteration of ceramide metabolism in men, specifically, Cer(d18:1_20:0) and Cer(d18:1_18:0) were significantly decreased only in men. Reduced levels of very long chain ceramide species have been associated with the development of macroalbuminuria [55], while male sex has been reported as a risk factor for the development of macroalbuminuria associated with T1D [56]. Moreover, we revealed a panel of LPCs significantly increased only in men. It has been shown that LPA and LPCs accumulate in the kidney and promote renal inflammation and tubulo-interstitial fibrosis in diabetic rodent models. Six species of LPAs and LPCs were found to be significantly enriched in the urine, but not in plasma, of people with T2D with nephropathy [57]. The mentioned sex-related lipidic differences could help to explain the worse prognosis of T1D-related diabetic nephropathy in men compared to women.

Strengths and limitations

Our study has several strengths, such as the large number of covariables used in our linear models to minimize confounding, the untargeted approach that allows for a more comprehensive characterisation of the lipidome in people with diabetes, and the consideration of sex-specific lipid differences associated with diabetes. There are also some limitations. First, our findings have not been validated in an independent cohort, and secondly, the observational nature of our study does not allow us to make causal inference. Therefore, further research is required to assess diabetes progression and its complications. Interestingly, our study shows the need to investigate this matter in a sex-specific manner.

Conclusions

In conclusion, we detected a panel of lipids associated with T1D and T2D, sex-specific differences in lipid metabolism disruption associated with diabetes and lipids associated with the glycaemic state, by considering the prediabetes stage. A large part of the lipids reported in this study have previously been linked to T1D, T2D and/or their complications in the literature, thus confirming their role in diabetes. Regarding sex-specific differences, we reported several lipid species associated with T2D only in women that have been previously related to menopause. This could help explain an unfavourable prognosis of T2D in women of older age compared to their male counterparts. In a similar way, we have shown a set of lipids associated with T1D only in men that have been previously linked to diabetic nephropathy, potentially explaining the worse prognosis of diabetic nephropathy in men. Our findings point to the need of establishing sex-specific strategies in the management and research on diabetes mellitus and its associated comorbidities and suggest the importance of lipidomics in advancing personalized medicine.

Data availability

No datasets were generated or analysed during the current study.

References

  1. Kerner AW, Brückel J, Definition. Classification and diagnosis of diabetes Mellitus. Exp Clin Endocrinol Diabetes. 2014;122:384–6. https://doi.org/10.1055/s-0034-1366278.

    Article  CAS  PubMed  Google Scholar 

  2. Saeedi P, Petersohn I, Salpea P, Malanda B, Karuranga S, Unwin N et al. Global and regional diabetes prevalence estimates for 2019 and projections for 2030 and 2045: Results from the International Diabetes Federation Diabetes Atlas, 9th edition. 2019. https://doi.org/10.1016/j.diabres.2019.107843.

  3. Rojo-Martínez G, Valdés S, Soriguer F, Vendrell J, Urrutia I, Pérez V, et al. Incidence of diabetes mellitus in Spain as results of the nation-wide cohort di@bet.es study. Sci Rep. 2020;10. https://doi.org/10.1038/s41598-020-59643-7.

  4. Wenk MR. The emerging field of lipidomics. Nat Rev Drug Discov 2005.

  5. Castelblanco E, Hernández M, Ortega E, Amigó N, Real J, Granado-Casas M, et al. Outstanding improvement of the advanced lipoprotein profile in subjects with new-onset type 1 diabetes mellitus after achieving optimal glycemic control. Diabetes Res Clin Pract. 2021;182. https://doi.org/10.1016/j.diabres.2021.109145.

  6. Kim EJ, Ramachandran R, Wierzbicki AS. Lipidomics in diabetes. Curr Opin Endocrinol Diabetes Obes. 2022;29:124–30. https://doi.org/10.1097/MED.0000000000000704.

    Article  CAS  PubMed  Google Scholar 

  7. Suvitaival T. Lipidomic abnormalities during the pathogenesis of type 1 diabetes: a quantitative review. Curr Diab Rep. 2020;20. https://doi.org/10.1007/s11892-020-01326-8.

  8. Julve J, Genua I, Quifer-Rada P, Yanes Ó, Barranco-Altirriba M, Hernández M, et al. Circulating metabolomic and lipidomic changes in subjects with new-onset type 1 diabetes after optimization of glycemic control. Diabetes Res Clin Pract. 2023;197. https://doi.org/10.1016/j.diabres.2023.110578.

  9. Cardner M, Yalcinkaya M, Goetze S, Luca E, Balaz M, Hunjadi M, et al. Structure-function relationships of HDL in diabetes and coronary heart disease. JCI Insight. 2020;5. https://doi.org/10.1172/jci.insight.131491.

  10. Barber MN, Risis S, Yang C, Meikle PJ, Staples M, Febbraio MA, et al. Plasma lysophosphatidylcholine levels are reduced in obesity and type 2 diabetes. PLoS ONE. 2012;7. https://doi.org/10.1371/journal.pone.0041456.

  11. Rauschert S, Uhl O, Koletzko B, Kirchberg F, Mori TA, Huang RC, et al. Lipidomics reveals associations of phospholipids with obesity and insulin resistance in young adults. J Clin Endocrinol Metab. 2016;101:871–9. https://doi.org/10.1210/jc.2015-3525.

    Article  CAS  PubMed  Google Scholar 

  12. Diamanti K, Cavalli M, Pan G, Pereira MJ, Kumar C, Skrtic S, et al. Intra- and inter-individual metabolic profiling highlights carnitine and lysophosphatidylcholine pathways as key molecular defects in type 2 diabetes. Sci Rep. 2019;9. https://doi.org/10.1038/s41598-019-45906-5.

  13. Zhong H, Fang C, Fan Y, Lu Y, Wen B, Ren H et al. Lipidomic profiling reveals distinct differences in plasma lipid composition in healthy, prediabetic, and type 2 diabetic individuals. Gigascience 2017;6. https://doi.org/10.1093/gigascience/gix036.

  14. Ding C, Wang N, Wang Z, Yue W, Li B, Zeng J, et al. Integrated Analysis of Metabolomics and Lipidomics in plasma of T2DM patients with Diabetic Retinopathy. Pharmaceutics. 2022;14. https://doi.org/10.3390/pharmaceutics14122751.

  15. Afshinnia F, Reynolds EL, Rajendiran TM, Soni T, Byun J, Savelieff MG, et al. Serum lipidomic determinants of human diabetic neuropathy in type 2 diabetes. Ann Clin Transl Neurol. 2022;9:1392–404. https://doi.org/10.1002/acn3.51639.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Xu T, Xu X, Zhang L, Zhang K, Wei Q, Zhu L, et al. Lipidomics reveals serum specific lipid alterations in Diabetic Nephropathy. Front Endocrinol (Lausanne). 2021;12. https://doi.org/10.3389/fendo.2021.781417.

  17. Ohkuma T, Peters SAE, Woodward M. Sex differences in the association between diabetes and cancer: a systematic review and meta-analysis of 121 cohorts including 20 million individuals and one million events. Diabetologia. 2018. https://doi.org/10.1007/s00125-018-4664-5.

  18. Alonso N, Traveset A, Rubinat E, Ortega E, Alcubierre N, Sanahuja J, et al. Type 2 diabetes-associated carotid plaque burden is increased in patients with retinopathy compared to those without retinopathy. Cardiovasc Diabetol. 2015;14:1–9. https://doi.org/10.1186/s12933-015-0196-1.

    Article  Google Scholar 

  19. Carbonell M, Castelblanco E, Valldeperas X, Betriu À, Traveset A, Granado-Casas M, et al. Diabetic retinopathy is associated with the presence and burden of subclinical carotid atherosclerosis in type 1 diabetes. Cardiovasc Diabetol. 2018;17:1–10. https://doi.org/10.1186/s12933-018-0706-z.

    Article  CAS  Google Scholar 

  20. Vilanova MB, Falguera M, Marsal JR, Rubinat E, Alcubierre N, Catelblanco E, et al. Prevalence, clinical features and risk assessment of pre-diabetes in Spain: the prospective Mollerussa cohort study. BMJ Open. 2017;7. https://doi.org/10.1136/bmjopen-2016-015158.

  21. Catalan M, Herreras Z, Pinyol M, Sala-Vila A, Amor AJ, de Groot E, et al. Prevalence by sex of preclinical carotid atherosclerosis in newly diagnosed type 2 diabetes. Nutr Metabolism Cardiovasc Dis. 2015;25:742–8. https://doi.org/10.1016/j.numecd.2015.04.009.

    Article  CAS  Google Scholar 

  22. Levey AS, Stevens LA, Schmid CH, Zhang Y, Castro AF III, Feldman HI, et al. A New equation to Estimate glomerular filtration rate. Ann Intern Med. 2009;9:604–12.

    Article  Google Scholar 

  23. Falguera M, Castelblanco E, Rojo-López MI, Vilanova MB, Real J, Alcubierre N, et al. Mediterranean diet and healthy eating in subjects with prediabetes from the mollerussa prospective observational cohort study. Nutrients. 2021;13:1–10. https://doi.org/10.3390/nu13010252.

    Article  Google Scholar 

  24. American Diabetes Association. Classification and diagnosis of diabetes. Diabetes Care. 2017;40:S11–24. https://doi.org/10.2337/dc17-S005.

    Article  Google Scholar 

  25. Chambers MC, et al. A cross-platform toolkit for Mass Spectrometry and Proteomics. Businesslive. 2012;30:1. https://doi.org/10.1038/nbt.2377.A.

    Article  Google Scholar 

  26. Smith CA, Want EJ, O’Maille G, Abagyan R, Siuzdak G. XCMS: Processing mass spectrometry data for metabolite profiling using nonlinear peak alignment, matching, and identification. Anal Chem. 2006;78:779–87. https://doi.org/10.1021/ac051437y.

    Article  CAS  PubMed  Google Scholar 

  27. Libiseller G, Dvorzak M, Kleb U, Gander E, Eisenberg T, Madeo F, et al. IPO: a tool for automated optimization of XCMS parameters. BMC Bioinformatics. 2015;16:1–10. https://doi.org/10.1186/s12859-015-0562-8.

    Article  CAS  Google Scholar 

  28. Kirwan JA, Broadhurst DI, Davidson RL, Viant MR. Characterising and correcting batch variation in an automated direct infusion mass spectrometry (DIMS) metabolomics workflow. Anal Bioanal Chem. 2013;405:5147–57. https://doi.org/10.1007/s00216-013-6856-7.

    Article  CAS  PubMed  Google Scholar 

  29. Jankevics A, Lloyd GR, Weber RJM. pmp: Peak Matrix Processing and signal batch correction for metabolomics datasets 2021.

  30. Lloyd GR, Weber RJM. structToolbox: Data processing & analysis tools for Metabolomics and other omics 2022.

  31. Lloyd GR, Jankevics A, Weber RJM. Struct: an R/Bioconductor-based framework for standardized metabolomics data analysis and beyond. Bioinformatics. 2020;36:5551–2. https://doi.org/10.1093/bioinformatics/btaa1031.

    Article  CAS  PubMed Central  Google Scholar 

  32. Subirana I, Héctor C, Cresib S, Vila J. Building Bivariate Tables: The compareGroups Package for R. vol. 57. 2014.

  33. R Core Team. R: A Language and Environment for Statistical Computing 2021.

  34. Dieterle F, Ross A, Schlotterbeck G, Senn H. Probabilistic quotient normalization as robust method to account for dilution of complex biological mixtures. Application in1H NMR metabonomics. Anal Chem. 2006;78:4281–90. https://doi.org/10.1021/ac051632c.

    Article  CAS  PubMed  Google Scholar 

  35. Rousset X, Vaisman B, Amar M, Sethi AA, Remaley AT. Lecithin:Cholesterol Acyltransferase: From Biochemistry to Role in Cardiovascular Disease. n.d.

  36. La Torre D, Seppänen-Laakso T, Larsson HE, Hyötyläinen T, Ivarsson SA, Lernmark Å, et al. Decreased cord-blood phospholipids in young age-at-onset type 1 diabetes. Diabetes. 2013;62:3951–6. https://doi.org/10.2337/db13-0215.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Orešič M, Simell S, Sysi-Aho M, Näntö-Salonen K, Seppänen-Laakso T, Parikka V, et al. Dysregulation of lipid and amino acid metabolism precedes islet autoimmunity in children who later progress to type 1 diabetes. J Exp Med. 2008;205:2975. https://doi.org/10.1084/JEM.20081800.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Małachowska B, Pietrowska K, Młynarski W, Szadkowska A, Krętowski A, Ciborowski M, et al. Multi-timepoint metabolic fingerprinting of a Post-episode period of hypoglycemia and ketoacidosis among children with type 1 diabetes. Front Mol Biosci. 2022;9. https://doi.org/10.3389/fmolb.2022.869116.

  39. Floegel A, Stefan N, Yu Z, Mühlenbruch K, Drogan D, Joost H-G et al. Identification of serum metabolites Associated with risk of type 2 diabetes using a targeted Metabolomic Approach n.d. https://doi.org/10.2337/db12-0495.

  40. Rajab A, Rabizadeh S, Bitaraf M, Ghanei A, Esteghamati A, Nakhjavani M, et al. Lecithin-cholesterol acyltransferase (LCAT) activity is lower in patients with type 2 diabetes in the presence of metabolic syndrome. Bioact Compd Health Dis. 2020;3:66–73. https://doi.org/10.31989/bchd.v3i4.690.

    Article  Google Scholar 

  41. Ghanei A, Esfahanian F, Esteghamati A, Nakhjavani M. Lecithin cholesterol acyltransferase activity is decreased in type 2 diabetes mellitus. 2007.

  42. Nakhjavani M, Esteghamati A, Esfahanian F, Ghanei A, Rashidi A, Hashemi S. HbA1c negatively correlates with LCAT activity in type 2 diabetes. Diabetes Res Clin Pract. 2008;81:38–41. https://doi.org/10.1016/j.diabres.2008.01.018.

    Article  CAS  PubMed  Google Scholar 

  43. Brenner RR. Antagonism between type 1 and type 2 diabetes in unsaturated fatty acid biosynthesis. Future Lipidol. 2006;1:631–40. https://doi.org/10.2217/17460875.1.5.631.

    Article  CAS  Google Scholar 

  44. Dean JM, Lodhi IJ. Structural and functional roles of ether lipids. Protein Cell. 2018;9:196–206. https://doi.org/10.1007/s13238-017-0423-5.

    Article  CAS  PubMed  Google Scholar 

  45. Hannich JT, Loizides-Mangold U, Sinturel F, Harayama T, Vandereycken B, Saini C, et al. Ether lipids, sphingolipids and toxic 1-deoxyceramides as hallmarks for lean and obese type 2 diabetic patients. Acta Physiol. 2021;232. https://doi.org/10.1111/apha.13610.

  46. Lu J, ManLam S, Wan Q, Shi L, Huo Y, Chen L, et al. High-coverage targeted lipidomics reveals novel serum lipid predictors and lipid pathway dysregulation antecedent to type 2 diabetes onset in normoglycemic Chinese adults. Diabetes Care. 2019;42:2117–26. https://doi.org/10.2337/dc19-0100.

    Article  CAS  PubMed  Google Scholar 

  47. Gonzalez-Covarrubias V. Lipidomics in longevity and healthy aging n.d. https://doi.org/10.1007/s10522-013-9450-7.

  48. Roszczyc-Owsiejczuk K, Zabielski P. Sphingolipids as a culprit of mitochondrial dysfunction in insulin resistance and type 2 diabetes. Front Endocrinol (Lausanne). 2021;12. https://doi.org/10.3389/fendo.2021.635175.

  49. Shen Q, Pierce JD. Supplementation of coenzyme Q10 among patients with type 2 diabetes mellitus. Healthc (Switzerland). 2015;3:296–309. https://doi.org/10.3390/healthcare3020296.

    Article  Google Scholar 

  50. Mielke MM, Venkata V, Bandaru R, Han D, An Y, Resnick SM et al. Demographic and clinical variables affecting mid-to late-life trajectories of plasma ceramide and dihydroceramide species n.d. https://doi.org/10.1111/acel.12369.

  51. Muilwijk M, Callender N, Goorden S, Vaz FM, van Valkengoed IGM. Sex differences in the association of sphingolipids with age in Dutch and south Asian Surinamese living in Amsterdam, the Netherlands. Biol Sex Differ 2021. https://doi.org/10.1186/s13293-020-00353-0.

  52. Arosio B, Corbi G, Davinelli S, Giordano V, Liccardo D, Rapacciuolo A, et al. Sex differences in Cardiovascular diseases: a matter of Estrogens, Ceramides, and sphingosine 1-Phosphate. Int J Mol Sci. 2022;23. https://doi.org/10.3390/ijms23074009.

  53. Slopien R, Wender-Ozegowska E, Rogowicz-Frontczak A, Meczekalski B, Zozulinska-Ziolkiewicz D, Jaremek JD, et al. Menopause and diabetes: EMAS clinical guide. Maturitas. 2018;117:6–10. https://doi.org/10.1016/J.MATURITAS.2018.08.009.

    Article  CAS  PubMed  Google Scholar 

  54. Kautzky-Willer A, Harreiter J, Pacini G. Sex and gender differences in risk, pathophysiology and complications of type 2 Diabetes Mellitus 2016. https://doi.org/10.1210/er.2015-1137.

  55. Klein RL, Hammad SM, Baker NL, Hunt KJ, Al Gadban MM, Cleary PA, et al. Decreased plasma levels of select very long chain ceramide species are associated with the development of nephropathy in type 1 diabetes. Metabolism. 2014;63:1287–95. https://doi.org/10.1016/j.metabol.2014.07.001.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Raile K, Galler A, Hofer S, Herbst A, Dunstheimer D, Busch P et al. Diabetic Nephropathy in 27,805 children, adolescents, and adults with type 1 diabetes effect of diabetes duration, A1C, hypertension, dyslipidemia, diabetes onset, and sex 2007. https://doi.org/10.2337/dc07-0282.

  57. Saulnier-Blache J-S, Feigerlova E, Halimi JM, Gourdy P, Roussel R, Guerci B et al. Urinary lysophopholipids are increased in diabetic patients with nephropathy 2017. https://doi.org/10.1016/j.jdiacomp.2017.04.024.

Download references

Acknowledgements

We want to particularly acknowledge the participants, and the IGTP-HUGTP (B.0000643) and IRBLleida Biobank (B.0000682) integrated in the PLATAFORMA BIOBANCOS (PT20/00050 and PT20/00021, respectively). The authors are also thankful to the COST Action EpiLipidNET, CA19105 - Pan-European Network in Lipidomics and EpiLipidomics. The first author gratefully acknowledges the Universitat Politècnica de Catalunya and Banco Santander for the financial support of her predoctoral grant. We would like to thank Joana Rossell and Minerva Granado-Casas for her feedback during the revision process. The authors acknowledge Amanda Prowse (Lochside Medical Communications Ltd.) for support in editing the paper.

Funding

This work was funded by Spanish Ministry of Health, Instituto de Salud Carlos III (Madrid, Spain) grants PI15/0625 (to DM and EC), PI17/01362 (to NA), PI18/0328 (to DM), FEDER “Una manera de hacer Europa”, and by Fundació La Marató de TV3 2016 (303/C/2016) (201602.30.31) (to NA). Spanish Ministry of Economy and Competitiveness, grant PID2021-122952OB-I00 funded by AEI 10.13039/501100011033 and by ERDF A way of making Europe, Instituto de Salud Carlos III (grant AC22/00035); and the CERCA Programme / Generalitat de Catalunya (to AP). This research was supported by CIBER-Consorcio Centro de Investigación Biomèdica en Red-CIBERDEM (leading group CB15/00071) and Biomedical Research Networking Center in Bioengineering, Biomaterials and Nanomedicine (CIBER-BBN), Instituto de Salud Carlos III, Ministry of Science and Innovation. Institut de Recerca de l’Hospital de la Santa Creu i Sant Pau is accredited by the Generalitat de Catalunya as Centre de Recerca de Catalunya (CERCA). B2SLab is certified as 2021 SGR 01052.

Author information

Authors and Affiliations

  1. Department of Endocrinology & Nutrition, Hospital de la Santa Creu i Sant Pau, Barcelona, Spain

    Maria Barranco-Altirriba & Didac Mauricio

  2. Departament d’Enginyeria de Sistemes, Automàtica i Informàtica Industrial, Universitat Politècnica de Catalunya, B2SLab, Barcelona, Spain

    Maria Barranco-Altirriba & Alexandre Perera-Lluna

  3. Networking Biomedical Research Centre in the subject area of Bioengineering, Biomaterials and Nanomedicine, CIBER-BBN), Madrid, Spain

    Maria Barranco-Altirriba & Alexandre Perera-Lluna

  4. Institut de Recerca Sant Joan de Déu, Esplugues de Llobregat, Barcelona, Spain

    Maria Barranco-Altirriba & Alexandre Perera-Lluna

  5. CIBER de Diabetes y Enfermedades Metabólicas Asociadas, Instituto de Salud Carlos III, Barcelona, Spain

    Núria Alonso, Oscar Yanes, Jordi Capellades, Josep Franch-Nadal & Didac Mauricio

  6. Servicio de Endocrinología y Nutrición, Hospital Universitario e Instituto de Investigación en Ciencias de la Salud Germans Trias i Pujol, Badalona, Barcelona, Spain

    Núria Alonso

  7. Departament de Medicina, Universitat Autònoma de Barcelona (UAB), Barcelona, Spain

    Núria Alonso

  8. School of Biosciences, University of Birmingham, Edgbaston, Birmingham, B15 2TT, UK

    Ralf J. M. Weber, Gavin R. Lloyd, Catherine Winder & Warwick B Dunn

  9. Phenome Centre Birmingham, University of Birmingham, Edgbaston, Birmingham, B15 2TT, UK

    Ralf J. M. Weber, Gavin R. Lloyd, Catherine Winder & Warwick B Dunn

  10. Institute of Metabolism and Systems Research, University of Birmingham, Edgbaston, Birmingham, B15 2TT, UK

    Ralf J. M. Weber & Gavin R. Lloyd

  11. Department of Endocrinology & Nutrition, Hospital Universitari Arnau de Vilanova, Institut de Recerca Biomèdica de Lleida (IRBLleida), Lleida, Spain

    Marta Hernandez

  12. Department of Electronic Engineering, Universitat Rovira i Virgili, IISPV, Tarragona, Spain

    Oscar Yanes

  13. Institut Investigació Sanitària Pere Virgili (IISPV), Tarragona, Spain

    Jordi Capellades

  14. Biomolecular Mass Spectrometry and Proteomics, Bijvoet Center for Biomolecular Research, Utrecht Institute for Pharmaceutical Sciences, University of Utrecht, Utrecht, The Netherlands

    Andris Jankevics

  15. Netherlands Proteomics Center, Utrecht, The Netherlands

    Andris Jankevics

  16. Centre for Metabolomics Research, Department of Biochemistry, Cell and Systems Biology, Institute of Systems, Molecular and Integrative Biology, University of Liverpool, L69 7ZB, Liverpool, UK

    Catherine Winder & Warwick B Dunn

  17. Institut d’Investigació Biomèdica, Centre Atenció Primària Cervera, Gerència d’Atenció Primària, Universitat de Lleida, Institut Català de la Salut, Lleida, Spain

    Mireia Falguera

  18. DAP-Cat Group, Unitat de Suport a La Recerca Barcelona Ciutat, Institut Universitari d’Investigació en Atenció Primària Jordi Gol, Barcelona, Spain

    Josep Franch-Nadal

  19. Division of Endocrinology, Metabolism and Lipid Research, Department of Internal Medicine, Washington University School of Medicine, 63110, St. Louis, MO, USA

    Esmeralda Castelblanco

  20. Unitat de Suport a la Recerca Barcelona, Institut Universitari d’Investigació en Atenció Primària Jordi Gol i Gurina, 08007, Barcelona, Spain

    Esmeralda Castelblanco

  21. Institut d’Investigació Biomèdica Sant Pau (IR Sant Pau), 08041, Barcelona, Spain

    Didac Mauricio

  22. Faculty of Medicine, University of Vic, Vic, Spain

    Didac Mauricio

Authors
  1. Maria Barranco-Altirriba
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Núria Alonso
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Ralf J. M. Weber
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Gavin R. Lloyd
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Marta Hernandez
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Oscar Yanes
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Jordi Capellades
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Andris Jankevics
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Catherine Winder
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Mireia Falguera
    View author publications

    You can also search for this author in PubMed Google Scholar

  11. Josep Franch-Nadal
    View author publications

    You can also search for this author in PubMed Google Scholar

  12. Warwick B Dunn
    View author publications

    You can also search for this author in PubMed Google Scholar

  13. Alexandre Perera-Lluna
    View author publications

    You can also search for this author in PubMed Google Scholar

  14. Esmeralda Castelblanco
    View author publications

    You can also search for this author in PubMed Google Scholar

  15. Didac Mauricio
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

DM, NA and EC conceived and designed the study; WD, CW, MH, MF, JF, OY, JC participated in data acquisition, RW, GL and AJ participated in lipidomic data processing; MB and AP contributed to statistical analysis and the design of figures and tables; MB and EC drafted the manuscript; DM, AP, EC, MB, RW and GL contributed to expert review, data interpretation and literature review; DM, AP and EC supervised the study; DM, EC, NA contributed with funding acquisition. All authors have reviewed and agree to the published version of the manuscript. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Corresponding authors

Correspondence to Esmeralda Castelblanco or Didac Mauricio.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the local ethics committee of the University Hospital Germans Trias i Pujol (PI-15-147), following the principles of the Declaration of Helsinki. All participants provided written informed consent.

Consent for publication

Not applicable.

Competing interests

Prof. Mauricio is a co-author of this study and an Editorial Board member of the Cardiovascular diabetology journal. He was not involved in handling this manuscript during the submission and the review processes.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

12933_2024_2202_MOESM1_ESM.pdf

Additional File 1: Figure S1. Flow chart of the participants recruitment. Table S1. Description of the analyses performed. Table S2. Baseline characteristics for each comparison. Figure S2. Upset plot of the positive acquisition mode results. Figure S3. Upset plot of the negative acquisition mode results.

12933_2024_2202_MOESM2_ESM.xlsx

Additional File 2: Table S3. Mass-to-charge ratio (mz) and retention time (rt) for each lipid ion significantly associated with T1D and/or T2D in at least one of the analyses.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Barranco-Altirriba, M., Alonso, N., Weber, R.J.M. et al. Lipidome characterisation and sex-specific differences in type 1 and type 2 diabetes mellitus. Cardiovasc Diabetol 23, 109 (2024). https://doi.org/10.1186/s12933-024-02202-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12933-024-02202-5

Keywords

Cardiovascular Diabetology

ISSN: 1475-2840