Journal of Clinical Medicine Article The Footprint of Type 1 Diabetes on Red Blood Cells: A Metabolomic and Lipidomic Study José Raul Herance 1,2, * , Andreea Ciudin 3,4 , Rubén Lamas-Domingo 5 , Carolina Aparicio-Gómez 1,2 , Cristina Hernández 2,3 , Rafael Simó 2,3 and Martina Palomino-Schätzlein 5,6, * 1 Medical Molecular Imaging Research Group, Vall d’Hebron Research Institute and Autonomous University of Barcelona, 08035 Barcelona, Spain; carolina.aparicio@vhir.org 2 CIBER-bbn (ISCIII), 28040 Madrid, Spain; cristina.hernandez@vhir.org (C.H.); rafael.simo@vhir.org (R.S.) 3 Diabetes and Metabolism Research Unit, Vall d’Hebron Research Institute, Autonomous University of Barcelona, 08035 Barcelona, Spain; aciudin@vhir.org 4 CIBERDEM (ISCIII), 28040 Madrid, Spain 5 NMR Facility, Centro de Investigación Príncipe Felipe, 46013 Valencia, Spain; ruben.lamas12@gmail.com 6 ProtoQSAR SL, CEEI (Centro Europeo de Empresas Innovadoras), Parque Tecnológico de Valencia, 46980 Valencia, Spain * Correspondence: raul.herance@vhir.org (J.R.H.); mpalomino@protoqsar.com (M.P.-S.); Tel.: +34-9-3489-3000 (ext. 4946) (J.R.H.); +34-9-6202-1811 (M.P.-S.) Abstract: The prevalence of diabetes type 1 (T1D) in the world populations is continuously growing. Although treatment methods are improving, the diagnostic is still symptom-based and sometimes far after onset of the disease. In this context, the aim of the study was the search of new biomarkers of the disease in red blood cells (RBCs), until now unexplored. The metabolomic and the lipidomic profile of RBCs from T1D patients and matched healthy controls was determined by NMR spectroscopy, and different multivariate discrimination models were built to select the metabolites and lipids that change most significantly. Relevant metabolites were further confirmed by univariate statistical analysis. Robust separation in the metabolomic and lipidomic profiles of RBCs from patients and controls was confirmed by orthogonal projection on latent structure discriminant analysis (OPLS-DA), Citation: Herance, J.R.; Ciudin, A.; random forest analysis, and significance analysis of metabolites (SAM). The main changes were Lamas-Domingo, R.; detected in the levels of amino acids, organic acids, creatine and phosphocreatine, lipid change length, Aparicio-Gómez, C.; Hernández, C.; and choline derivatives, demonstrating changes in glycolysis, BCAA metabolism, and phospholipid Simó, R.; Palomino-Schätzlein, M. metabolism. Our study proves that robust differences exist in the metabolic and lipidomic profile The Footprint of Type 1 Diabetes on of RBCs from T1D patients, in comparison with matched healthy individuals. Some changes were Red Blood Cells: A Metabolomic and similar to alterations found already in RBCs of T2D patients, but others seemed to be specific for Lipidomic Study. J. Clin. Med. 2023, type 1 diabetes. Thus, many of the metabolic differences found could be biomarker candidates for an 12, 556. https://doi.org/10.3390/ earlier diagnosis or monitoring of patients with T1D. jcm12020556 Academic Editor: Tomoaki Morioka Keywords: type 1 diabetes; NMR spectroscopy; metabolomics; lipidomics; red blood cells Received: 29 November 2022 Revised: 4 January 2023 Accepted: 7 January 2023 Published: 10 January 2023 1. Introduction Diabetes mellitus is the first non-communicable disease that the United Nations recognized as a 21st-century pandemic problem [1]. Diabetes is a group of diseases mainly characterized by high blood glucose levels. Type 1 diabetes (T1D) and type 2 diabetes (T2D) Copyright: © 2023 by the authors. both involve hyperglycemia, but are disorders with different patterns, mechanisms, clinical Licensee MDPI, Basel, Switzerland. presentations, and disease progressions. Nevertheless, once hyperglycemia occurs, patients This article is an open access article with all forms of diabetes can develop similar chronic complications, in different rates of distributed under the terms and progression and form [2]. conditions of the Creative Commons T2D accounts for about 80–90% of the cases of diabetes and is due to a progressive loss Attribution (CC BY) license (https:// of β-cell mass and insulin secretion, on the background of insulin resistance (IR), low-grade creativecommons.org/licenses/by/ 4.0/). inflammation, and oxidative stress [3]. In 2019, 463 million adults (20–79 years) were living J. Clin. Med. 2023, 12, 556. https://doi.org/10.3390/jcm12020556 https://www.mdpi.com/journal/jcm J. Clin. Med. 2023, 12, 556 2 of 14 with diabetes worldwide and this number could rise to 700 million by 2045 “www.idf.org (accessed on 6 January 2023)”. Due to the latent onset of this disease, there is also a long pre-detection period, and up to one-half of cases in the population may be undiagnosed [4]. T1D accounts for about 5–10% of all cases of diabetes and is due to the autoimmune destruction of beta cells of the pancreas, leading to a total incapacity of producing insulin [5,6]. The incidence and prevalence of T1D are increasing worldwide dramatically. While T2D starts generally at an older age (average diagnosis age 45–52 years) [7], T1D tends to affect much younger people (average diagnosis age 10–14 years) [8], children in many cases, although the age limit is not mandatory. Patients with T1D need complete substitutive treatment with insulin for their lifetime (both long-acting and rapid-acting insulin) and have been classically described to have unimpaired insulin resistance (Imran 2008). Nevertheless, recent scientific reports suggest that people with T1D also develop insulin resistance (IR), which is related to complications and comorbidities in patients [9]. Furthermore, mechanisms such as lipotoxicity, glucotoxicity, closely related to hy- perglycemia [10], and iatrogenic hyperinsulinemia [11] cause oxidative stress that can contribute to the development of IR in T1D. In oxidative stress, reactive oxygen species (ROS) can induce insulin receptor substrate serine/threonine phosphorylation, impair cellular redistribution of insulin signaling components, reduce GLUT4 gene transcription, or alter mitochondrial activity [12]. On the other hand, oxidative stress can increase the production of pro-inflammatory cytokines and induce inflammatory processes. Immune cells such as macrophages, which are the source of pro-inflammatory cytokines, can migrate to pancreatic islet cells [13]. Thus, chronic inflammation can also contribute to the develop- ment of IR in T1D. For instance, adipose cell enlargement can lead to a pro-inflammatory state and the formation of pro-inflammatory compounds, such as interleukin-6, C-reactive protein, and plasminogen activator inhibitor-1, which, in turn, contribute to the death of beta cells, modulation of beta cell regeneration processes, and IR [14,15]. Even if the main causes of IR in T1D have been explored, at present, there is a lack of data that connect the risk factors of IR in T1D with the metabolomic and lipidomic profile. Glycated hemoglobin (HbA1c) is the parameter used worldwide for the diagnosis and follow-up of the metabolic control of diabetes [16]. HbA1c is adult hemoglobin (HbA) with glucose bound to its b chain N-terminal valine, resulting from non-enzymatic glycation in erythrocytes, also known as red blood cells (RBCs) [17]. RBCs are the most abundant cell type in the human body and play essential roles in metabolism, transport, and respiration. Their composition, shape, abundance, and motility can be altered by diseases. Several studies revealed that RBCs suffer general alterations in diabetes [18]. For instance, RBCs’ aggregability increase [19] as well as alterations of RBC shape [20], deformability [21], and membrane phospholipid composition [22]. In addition, the glutathione redox system is affected [23] with decreases in reduced glutathione (GSH), GSH peroxidase and reductase activity, energetic metabolism, as well as alterations in adenosine triphosphate (ATP) levels [24]. Our group recently showed that the metabolomics profile of RBCs is generally altered in T2D disease [25], mainly affecting the redox and energetic metabolism, as well as the transport of certain metabolites necessary for physiological functions. However, at present, no information about the impact of T1D on the metabolomics profile of RBCs is known. On these bases, the aim of this study was to explore the general metabolic alteration of RBCs in patients with T1D to find fingerprints of insulin resistance and metabolic information for monitoring this type of patients. 2. Materials and Methods 2.1. Study Design and Subjects The study was designed as a cross-sectional study case-control. It was approved by the local Ethics Committee (PR(AG)234/2015) and carried out in accordance with the Declaration of Helsinki. Patients were recruited from the outpatient clinic of the Endocrinology Department of Vall d’Hebron University Hospital and divided into two J. Clin. Med. 2023, 12, 556 3 of 14 groups as follows: 20 patients with T1D matched by age, sex, and body mass index (BMI) with 20 healthy individuals. All patients signed the informed consent form before their inclusion in the study. The inclusion criteria were the following: (a) age > 18 years, (b) T1D, (c) written informed consent form. The exclusion criteria were the following: (a) any other type of diabetes, except for T1D (T2D, MODY, etc.), (b) tobacco use in the last 12 months, (c) alcohol intake of more than 30 g in males and 20 g in females, (d) hemoglobinopathies, (e) hemato- logic disease, (f) any other disease or treatment that might influence the RBCs’ metabolism (chronic kidney disease, cardiac valve prosthesis, chemotherapy, anticoagulants, etc.). All diabetic patients had positive islet autoantibodies of anti-GAD-65 (>50 IU/mL). This was a pilot study. At present, there are no available data regarding the impact of T1D on the RBCs’ metabolism; therefore, we were not able to perform a sample size calculation. Thus, the size has been selected by taking into account ethical considerations and the pilot studies in T2D [25,26]. All the subjects included in the study underwent a complete medical history review, anthropometric evaluation, and biochemical analysis. Anthropometrical parameters of weight (kg), waist (cm), and height (m) were measured just before the collecting of the blood samples. Biochemical analysis was performed using standardized and validated routine methodologies on the blood samples collected at fasting conditions at the biochemistry core facilities of Vall d’Hebron University Hospital. Plasmatic concentrations of glucose (mg/dL), glycosylated hemoglobin (HbA1c) (%DCCT), high-density lipoprotein (HDL) (mg/dL), low-density lipoprotein (LDL) (mg/dL), triglycerides (TG) (mg/dL), anti-GAD65 (IU/mL), and insulin (mU/mL) were measured in all the subjects. Body mass index (BMI) kg/m2 was calculated using the following formula: body weight (kg)/height2 (m) [27]. 2.2. Peripheral Blood Samples Blood samples were collected followed a similar protocol as in previous studies [25,26]. Briefly, peripheral blood was collected under fasting conditions, stored at 4 ◦ C, and processed within the first hour. Then, 5 mL of peripheral blood were poured into a quartz tube with 10 mL of Ficoll until 2 phases were separated by gravity, around 30 min. Then, the pellet was washed twice with 10 mL of cool Hank’s Balanced Salt Solution in a centrifuge at 200× g and 4 ◦ C for 10 min without breaks. Cell counting was performed on a Neubauer chamber and the purity was tested with flow cytometry. For storage, 200 µL of ice-cold methanol were added per 10 million cells for quenching and the samples were frozen directly at −80 ◦ C. 2.3. Extraction of Metabolites and Lipids Metabolites were extracted from RBCs as described previously [25]. Briefly, frozen samples were allowed to thaw for 5 min on ice. Then, 800 µL of chloroform at 4 ◦ C were added to each sample tube. After 10 min, the samples were homogenized, resuspended with a pipette, and transferred to a bigger plastic tube. Then, samples were placed in liquid nitrogen for 1 min and then allowed to thaw on ice for 2 min, for uniform cell membrane breakage. This step was performed two more times. Afterwards, 1250 µL of distilled water and 1250 µL of chloroform were added and the sample vortexed. Samples were then centrifuged at 13,000× g for 20 min at 4 ◦ C to separate the phases. The upper phase contained polar metabolites in a mixture of water/methanol and was lyophilized overnight. The lower chloroform phase containing the lipids was evaporated under nitrogen flux. Extracts were stored at −80 ◦ C until NMR sample preparation and analysis. For NMR sample preparation, the frozen samples were allowed to thaw for 5 min on ice. Then, 550 µL of NMR buffer (100 mM Na2 HPO4 in D2 O at pH 7.4) containing 0.1 mM deuterated trimethylsilylpropanoic acid as an internal standard were added to the polar samples and transferred to a 5 mm NMR tube. For the preparation of the lipidomics samples, 600 µL of deuterated chloroform containing 0.003% v/v tetramethylsilane was added to the lipidic samples and transferred to a 5 mm NMR tube. All samples were analyzed the same day and stored at 4 ◦ C until the analysis. J. Clin. Med. 2023, 12, 556 4 of 14 2.4. NMR Experiments NMR spectra were recorded as described previously [25] at 27 ◦ C on a Bruker AVII-600 spectrometer using a 5 mm triple resonance cryoprobe and processed using TopSpin version 3.2 software (Bruker BioSpin). Before starting sample measurement, the temperature was calibrated at 27 ◦ C with a 99.8% MeOD sample. Further, the parameters of the spectrometer were optimized to ensure an optimal resolution and sensitivity (width at half-height ≤ 1 Hz), with a standard sample containing 2 mM of sucrose, 0.5 mM of DSS (sodium trimethylsilyl- propanesulfonate), 2 mM of NaN3 in 90% H2 O, and 10% D2O. 1 H 1D NOESY NMR spectra were acquired with 256 free induction decays (FIDs). A 4s relaxation delay was incorporated between FIDs and the water signal was eliminated with presaturation in the aqueous samples. The FID values were multiplied by an exponential function with a 0.5 Hz line broadening factor for an optimal baseline correction, with 64K data points digitalized over a spectral width of 30 ppm. Quality control samples, consisting of a mixture of lactate, creatine, citric acid, phenylalanine, and glucose, were run in regular intervals to confirm the reproducibility of the measurements. 2.5. Data Analysis 1 H-NMR signals from the spectra were assigned with the help of previous data [25,28] and the spectral databases Human Metabolome Database and Biological Magnetic Resonance Bank [29,30]. In ambiguous cases, the assignment was confirmed by spiking the spectra with reference compounds. For metabolite quantification, spectra were normalized to total intensity. This allowed us to work with normalized concentration values and avoid the detection of changes related exclusively to differences in total sample amount and experimental error during the extraction process. The predominant glucose signals and the solvent signals were excluded from normalization. Optimal integration regions were defined for each metabolite, an attempt being made to select the signals without overlapping (Table S1). Integration was performed with global spectrum deconvolution in MestreNova 12 (Mestrelab Research S. L.). The datasets are available in the Zenodo repository [31]. For modelling, metabolite integration data were Pareto-scaled (each value being di- vided by the square root of the standard deviation of each variable) and mean-centered, for an easier interpretation of the data and to also take into account the variation of small signals. Orthogonal projection on latent structure discriminant analysis (OPLS-DA) and projection on latent structure (PLS) regression were performed with SIMCA-P 14.0 (Umet- rics, Sweden). OPLS-DA and PLS models were described with R2Y(cum) (representing the cumulative SS of all the y-variables explained by the extracted components) and Q2(cum) (describing the cumulative Q2 for all y-variables for the extracted components). Score plots (representing scores related to prediction Y t(1) versus scores related to the first orthogonal component t0) and S-plots (plotting the modelled covariation (p(1)) versus the modelled correlation (p(corr)) [32] were used to detect significant metabolites and lipids in the OPLS- DA models. The glucose signal, due to its high intensity and correlation to plasmatic glucose, was excluded from OPLS-DA analysis. OPLS-DA and PLS models were validated in SIMCA-P by a 100-times permutation (where the overfit of the model is measured by the intercept of the regression line of the correlation coefficient between the original y-variable and the permuted y-variable on the x-axis versus the cumulative R2 and Q2 on the y-axis) and analysis of variance testing of cross-validated predictive residuals (CV-ANOVA) [33]. CV-ANOVA is a significance test for the Q2YCV cross-validation using the F-distribution, based on an ANOVA assessment of the cross-validatory (CV) predictive residuals of the model. In PLS analysis, metabolites that were relevant for the correlation were identified by VIP (variable-importance-projection) values >1. Significance analysis of metabolites (SAM) analysis was performed with Metaboan- alyst [34] applying the siggenes package [35]. To each variable, a significance score was assigned based on its change relative to the standard deviation of repeated measurements. For a variable with scores greater than an adjustable threshold, its relative difference is compared to the distribution estimated by random permutations of the class labels. For J. Clin. Med. 2023, 12, 556 5 of 14 each threshold, a certain proportion of the variables in the permutation set will be found to be significant by chance. The proportion is used to calculate the FDR. In our analysis, a Delta value of 1.4 was defined to control FDR. Random Forest (RF) analysis was performed with Metaboanalyst [34] based on the randonForest package [36]. RF uses groups of classification trees, each of which is grown by random feature selection from a bootstrap sample at each branch. Class prediction is based on the majority vote of the ensemble. For initial tree construction, one-third of the samples are left out for validation. The OOB (out-of-bag) data are then used as a test sample to obtain an unbiased estimate of the classification error (OOB error). Variable importance is calculated by the measure of the increase of the OOB error when it is permuted. Hierarchical clustering was performed with Metaboanalyst with the hclust function in package stat and the result represented as a heatmap. Box plots and t-tests were performed with normalized concentration data in Graphpad Prism 8. Only metabolites that showed high scores in all three models (OPLS-DA, Random Forest, and SAM) were considered for univariate analysis. For the biological interpretation of the results and the identification of metabolic pathways, the Kegg Data Base and MetPA (Metaboanalyst) were used. 3. Results Table 1 summarizes the main characteristics of the subjects and controls included in the study. Table 1. Main characteristics of the cohort of T1D patients. T1D CT p Age (years) 32.3 ± 1.6 33.0 ± 2.4 0.821 Gender F/M 8/12 8/12 - HbA1c (DCCT%) 8.00 ± 0.29 5.12 ± 0.06 <0.001 Glucose (mg/dL) 168.7 ± 18.3 82.4 ± 2.2 <0.001 BMI (kg/m2 ) 24.1 ± 1.0 23.4 ± 0.8 0.604 Cholesterol (mg/dL) 187 ± 10 185 ± 7 0.892 HDL (mg/dL) 57 ± 3 60 ± 3 0.440 LDL (mg/dL) 116 ± 9 109 ± 6 0.525 Triacylglycerides 75 ± 9 84 ± 9 0.483 (mg/dL) CT = healthy control, F/M = female/male, HbA1c = glycosylated hemoglobin levels, BMI = body mass index. Data are shown as mean ± SD. p-values were calculated with the Student’s t-test. A metabolomic and a lipidomic analysis of RBCs of all individuals was performed by proton nuclear magnetic resonance (1 H-NMR). As a result, 68 different hydrophilic metabo- lites and 12 different types of lipid groups were identified and quantified (Supporting Tables S1 and S2). Three different discrimination methods were applied in order to analyze significant differences between both groups: OPLS-DA (Orthogonal Projections to Latent Structures Discriminant Analysis), Random Forest Classification, and SAM (Significance Analysis of Metabolites). Significant discrimination between healthy con- trols (CT) and T1D patients could be obtained for the three model types of the metabolomic and lipidomic profiles (Figure 1). In addition, a heatmap analysis was performed to have an overview of the main metabolic changes in the patient and control groups. As can be seen in Figure 2A, a similar tendency in metabolite levels in T1D patients is detected (yellow color), while controls have a higher dispersion (blue color). As a general tendency, T1D patients seem to have lower levels of amino acids and organic acids. Concerning lipid metabolism, T1D patients tended to have lower levels of unsaturated fatty acids (Figure 2B). This tendency was, however, less pronounced for polyunsaturated fatty acids (PUFA). Additionally, phosphatidylcholines were reduced as well as phospho- lipids in general. On the other hand, a decrease of lipid chain length could be detected in J. Clin. Med. 2023, 12, x FOR PEER REVIEW 6 of 14 to analyze significant differences between both groups: OPLS-DA (Orthogonal Projections J. Clin. Med. 2023, 12, 556 6 of 14 to Latent Structures Discriminant Analysis), Random Forest Classification, and SAM (Sig- nificance Analysis of Metabolites). Significant discrimination between healthy controls (CT) and T1D patients could be obtained for the three model types of the metabolomic T1Dlipidomic and patients, determined by the profiles (Figure 1).ratio of lipid methyl groups (Lipid CH3-, increase) versus linear lipid chains (Lipid–CH2, decrease). A OPLS-DA Metabolites B OPLS-DA Lipids CT CT T1D T1DM CT CT T1DM T1D T1DM 40 30 10 20 to[1] to[1] 10 0 1,0 0 Lipid -CH2- -10 0,8 Phospholipids -10 -20 3-methyladipate 0,8 Phosphatidycholine 0,6 -10 -5 0 5 10 -30 0,6-20 -10 0 10 20 30 t[1] t[1] 0.4 SIMCA-P+ 12.0.1 - 2020-08-29 23:56:57 (UTC+1) SIMCA-P+ 12.0.1 - 2020-08-29 23:29:27 (UTC+1) 0.4 p(corr)[1] p(corr)[1] 0.2 0.2 MUFA+PUFA Cholesterol -0.0 -0.0 Lipid CH3- 2 -0.2 Cholesterol -0.2 FA -CH2CH2COO- phosphocreatine -0.4 -0.4 creatine PUFA phosphocholine -0.6 Lipid CH3- 1 GSSG Cholesterol Phosphatidylethanolamine -0.6 -0.8 -0.30 -0.20 -0.10 -0.00 0.10 0.20 0.30 p[1] -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 p[1] C Random Forest Metabolites D Random Forest Lipids class class CT T1D error CT T1D error CT 18 2 0.1 CT 18 0 0.0 T1 2 18 0.1 T1 0 18 0.0 trees trees E SAM ANALYSIS Metabolites F SAM ANALYSIS Lipids cutlow: -2.609 cutlow: -1.896 cutup: Inf cutup: Inf p0 = 0.641 p0 = 0.370 Significant: 14 Significant: 42 False 0.57 False 3.24 FDR: 0.026 FDR: 0.029 Observed(i) values Observed(i) values Expected d(i) values Expected d(i) values Figure 1. Figure (A) Score 1. (A) Score and and S-plot S-plot from from OPLS-DA OPLS-DAanalysis analysiscomparison comparisonofofthetheaqueous aqueousRBC RBCprofile profileofof T1DM T1DM patients patients and andhealthy healthycontrols. controls.RY(cum) RY(cum)= 0.77, Q(cum) = 0.77, Q(cum)= 0.64, CV-ANOVA = 0.64, CV-ANOVA = 3.9 × 10× = 3.9 10−6 . −6. (B) Score and S-plot from OPLS-DA analysis comparison of the lipidomic RBC profile (B) Score and S-plot from OPLS-DA analysis comparison of the lipidomic RBC profile of T1DM of T1DM patients and healthy patients and controls. RY(cum)RY(cum) healthy controls. = 0.83, Q(cum) = 0.66, CV-ANOVA = 0.83, Q(cum) = 1.1× =101.1 = 0.66, CV-ANOVA −6. (C) Random Forest × 10−6 . (C) Random classification of aqueous RBC metabolites of T1D patients and healthy controls, 500 trees, 7 Forest classification of aqueous RBC metabolites of T1D patients and healthy controls, 500 trees, 7 predictors, randomness. (D) Random Forest classification of lipidic RBC metabolites of T1D patients and healthy controls, 500 trees, 5 predictors, fixed seed. (E) SAM analysis of aqueous metabolites, Delta = 1.4, FDR 0.026, significant compounds: 14. (F) SAM analysis of lipidic metabolites, Delta = 0.1, FDR 0.025, significant compounds: 10. In addition, a heatmap analysis was performed to have an overview of the main met- abolic changes in the patient and control groups. As can be seen in Figure 2A, a similar tendency in metabolite levels in T1D patients is detected (yellow color), while controls J. Clin. Med. 2023, 12, 556 have a higher dispersion (blue color). As a general tendency, T1D patients seem to have 7 of 14 lower levels of amino acids and organic acids. A METABOLITES class glutamine GSH GSSG threonine 2-3BPG ATP creatine phosphocholine phosphocreatine proline acetate GPC aspartate DMA lysine 3-methyl-2-oxovalerate isoleucine IMP sarcosine asparagine methionine alanine glutamate lactate formate pyroglutamate leucine 2-hydroxyisovalerate 3-methyladipate B LIPIDS class colesterol MUFA+PUFA PUFA phosphatidylethanol phosphatidylcholine phospholipids lipids CH3 FA –CH2CH21COO- lipids CH3 2 lipids –CH2- Figure Figure2.2.Heatmap Heatmapanalysis performed analysis withwith performed the Ward clustering the Ward algorithm clustering and ANOVA algorithm test. (A) and ANOVA test. Heatmap showing main aqueous metabolites. (B) Heatmap showing the main lipid groups. GSH = (A) Heatmap showing main aqueous metabolites. (B) Heatmap showing the main lipid groups. reduced glutathione, GSSG = oxidized glutathione, 2,3-BPG = 2,3-Bisphosphoglyceric acid, GPC = GSH = reduced glutathione, GSSG = oxidized glutathione, 2,3-BPG = 2,3-Bisphosphoglyceric acid, glycerophosphocholine, DMA = dimethylamine, IMP = inosinic acid, MUFA = monounsaturated GPC = glycerophosphocholine, DMA = dimethylamine, IMP = inosinic acid, MUFA = monoun- saturated fatty acids, PUFA = polyunsaturated fatty acids, CH3 = methyl group, FA = fatty acids, CH2 = lipid chain. The color scale refers to normalized concentration values. In order to make a robust metabolite selection, metabolites that had high scores in all three model types (Figure 1) were further submitted to univariate analysis. The resulting significant metabolites are represented in Figure 3, and can be divided into different metabolite groups. As can be seen, organic acids and metabolites related to lipid and choline metabolism seem to be the most affected. detected in T1D patients, determined by the ratio of lipid methyl groups (Lipid CH3-, increase) versus linear lipid chains (Lipid–CH2, decrease). In order to make a robust metabolite selection, metabolites that had high scores in all three model types (Figure 1) were further submitted to univariate analysis. The resulting significant metabolites are represented in Figure 3, and can be divided into different me- J. Clin. Med. 2023, 12, 556 tabolite groups. As can be seen, organic acids and metabolites related to lipid and choline 8 of 14 metabolism seem to be the most affected. ORGANIC ACIDS LIPID AND CHOLINE METABOLISM 2.0 *** 6 *** 15 8 * 300 *** 650 *** 600 1.5 6 4 10 200 550 A.U. A.U. A.U. A.U. A.U. 1.0 4 500 2 5 100 0.5 2 450 0.0 0 0 0 0 400 formate acetate lactate GPC Lipid CH3- Lipid -CH2- 2.5 *** 1.5 * 15 *** 8 * 80 *** 1.0 *** 2.0 0.8 6 60 1.0 10 1.5 0.6 A.U. A.U. A.U. A.U. A.U. A.U. 4 40 1.0 0.4 0.5 5 2 20 0.5 0.2 0.0 0.0 0 0 0 0.0 2-hydroxyisovalerate 3-methyl-2-oxovalerate 3-methyladipate phosphocholine Phosphatidycholine propyleneglycol BCAA METABOLISM NUCLEOTID METABOLISM CREATINE METABOLISM 6 15 1.5 1.0 8 40 ** * * *** ** ** 0.8 4 10 1.0 6 30 CT 0.6 A.U. A.U. A.U. A.U. A.U. A.U. 4 20 0.4 T1D 2 5 0.5 2 10 0.2 0 0 0.0 0.0 0 0 isoleucine leucine IMP sarcosine phosphocreatine creatine Figure Figure3. 3. Box plots of RBCofmetabolites Box plots that showed RBC metabolites that relevant showedchanges relevantin T1D patients. changes GPC patients. in T1D = glyc- erophosphocholine, CH3 = methyl group, CH2 = lipid chain, IMP = inosinic acid. Concentrations GPC = glycerophosphocholine, CH3 = methyl group, CH2 = lipid chain, IMP = inosinic acid. Con- were of normalized total intensity. * = p < 0.05, ** = p < 0.01, *** = p < 0,001 in Wilcoxon test. A.U. = centrations were of normalized total intensity. * = p < 0.05, ** = p < 0.01, *** = p < 0.001 in Wilcoxon arbitrary units from normalized concentration values. test. A.U. = arbitrary units from normalized concentration values. InInorder ordertotodetect detectwhich whichofofthethemetabolic metabolicchanges changeswereweredirectly directlycorrelated correlatedwith with plasma glucose levels in fasting conditions, we performed PLS regression plasma glucose levels in fasting conditions, we performed PLS regression analyses analyses of theof metabolomic and the the metabolomic andlipidomic profile the lipidomic versus profile glucose versus concentrations. glucose As a As concentrations. result, we ob- a result, we tained a good correlation of glucose with several metabolites (Figure S1A) and obtained a good correlation of glucose with several metabolites (Figure S1A) and a mod- a moderate correlation with some erate correlation withblood some lipids blood (Figure S1B). The lipids (Figure metabolites S1B). and theand The metabolites lipids the that were lipids that relevant for the correlation (VIP > 1) are listed in Tables S3 and S4. Most metabolites were relevant for the correlation (VIP > 1) are listed in Tables S3 and S4. Most metabolites that were significantly that were altered significantly in RBCs altered of T1D in RBCs patients of T1D also patients had also a correlation had a correlationwith withplasma plasma glucose levels. Only glycerophosphocholine (GPC), inosinic acid (IMP), and creatine were not relevantly correlated. Regarding lipidomics, although the correlation was much weaker, all relevant lipid groups (lipid CH3 , lipid CH2 , and phosphatidylcholine) seemed to be related to plasma glucose concentration. In order to correlate the metabolite changes with alterations of specific metabolic pathways, an enrichment pathway analysis was performed (Figure 4), showing that energy metabolism, fatty acid metabolism, and amino acid metabolism were the most affected metabolic pathways in T1D compared to healthy controls. not relevantly correlated. Regarding lipidomics, although the correlation was much weaker, all relevant lipid groups (lipid CH3, lipid CH2, and phosphatidylcholine) seemed to be related to plasma glucose concentration. In order to correlate the metabolite changes with alterations of specific metabolic pathways, an enrichment pathway analysis was performed (Figure 4), showing that en- J. Clin. Med. 2023, 12, 556 ergy metabolism, fatty acid metabolism, and amino acid metabolism were the most9 of af-14 fected metabolic pathways in T1D compared to healthy controls. Metabolite Sets Enrichment Overview Pyruvate Metabolism Glutathione Metabolism Valine, Leucine and Isoleucine Degradation Aspartate Metabolism Biotin Metabolism Pyruvaldehyde Degradation Glutamate Metabolism 6 × 10³ Phosphatidylcholine Metabolism Arginine and Proline Metabolism Glycine and Serine metabolism Ethanol Degradation Carnitine Synthesis × 101 40.4 Arachidonic Acid Metabolism Phospholid Biosynthesis Lysine Degradation Amino Sugar Metabolism Gluconeogenesis × 101 70.7 Fatty Acid Biosynthesis Retinol Metabolism Spingolipid Metabolism Methionine Metabolism Warburg Effect Purine Metabolism Enrichment Ratio Figure 4. Pathway enrichment overview of metabolomic changes in RBCs of T1D obtained Figure 4. Pathway enrichment overview of metabolomic changes in RBCs of T1D obtained with with Metaboanalyst. Metaboanalyst. 4. Discussion 4. Discussion In the present pilot study, we detected a clear difference in the metabolomic and lipidomicIn profile the present pilot study, of patients with we T1Ddetected a clear and control difference subjects, in thethat proving metabolomic T1D has aand lip- specific idomic profile of patients with T1D and control subjects, proving that T1D impact on RBC metabolism and transport. Our group has previously reported that the has a specific metabolomics profile of RBCs is generally altered in T2D disease [25], mainly affecting redox and energy metabolism, as well as the transport of certain metabolites with physiological functions. As far as we know, this is the first study that shows a significant impact of T1D on the metabolomics and lipidomics profile of RBCs. In our analysis, we detected a significant reduction of acetate in patients, as well as a tendency of lactate to decrease (Figure 2A). This could reflect a reduction in glycolysis, a simple pathway of glucose metabolism that is regulated by insulin secretion and sensi- bilization [37]. Previous studies reported increased plasmatic lactate and acetate levels in patients with diabetes [38]. Thus, the increment in lactate could be related to increased carbonic anhydrase activity in RBCs, transporting a higher amount of lactate into blood J. Clin. Med. 2023, 12, 556 10 of 14 plasma [39]. Our group reported similar findings in patients with T2D [25], which could indicate that similar processes happen in T1D and T2D. Other metabolites altered in RBCs were isoleucine and leucine, essential metabolites for the production and formation of hemoglobin and the production of RBCs [40]. This alteration is similar to that detected in T2D patients [25], where the BCCAs valine and leucine decreased. An opposite tendency was detected in plasma of diabetic patients where the concentration of BCAAs was increased in both T1D and T2D [41,42]. This inverse correlation could be related to the fact that RBCs with highly glycosylated hemoglobin, as occurs in both T1D and T2D diabetes, seem to have a lower BCAA uptake capacity. Interestingly, isoleucine was exclusively altered in RBCs of T1D patients, while significant changes in valine were only detected in T2D. The differences in the alteration of BCCAs depending of the type of diabetes could indicate a different pathophysiology associated with BCAAs in both types of diabetes, since this type of amino acids is transported by RBCs. Furthermore, valine and leucine could be RBCs’ fingerprints of T2D and T1D, respectively. In addition, the alteration of the organic acids 2-hydroxyisovalerate and 3-mehtyl-2-oxovalerate, key compounds in BCAAs metabolism, were only decreased in T1D, indicating a higher affectation of BCAAs in this type of diabetes. Previous studies reported changes in other kinds of amino acids associated with diabetes and pre-diabetes in many studies [43]. Interestingly, the changes in non-BCAA amino acids seemed to be less pronounced in RBCs of T1D patients than in T2D [25]. For instance, alanine, asparagine, glutamate, proline, threonine, and methionine seem not to be affected by T1D, while they change significantly in T2D. These altered amino acid levels may directly interfere with insulin-regulated glycogen synthesis, and the lower impact on amino acid concentration in RBCs may be related to a lower level of IR in T1D, reducing the enhancement of plasma concentration of these type of metabolites [44]. However, further studies are mandatory to clarify this assumption. Another noteworthy result from our study was the increase of creatine and phos- phocreatine in RBCs of T1D patients. Plasmatic creatine has been recently identified as a potential biomarker for mitochondria dysfunction and for T2D risk in the PREVEND study from Post et al [45]. Our results indicate that this increase could be more generalized in the organism, taking place in cells, and may be related to alterations in muscle or en- ergy metabolism where creatine plays an important role. This is further supported by the increase of inosine monophosphate (IMP) showed in RBCs of T1D patients, a nucleotide related to energy and generated in muscles during exercise [46]. Interestingly, neither creatine nor IMP showed a direct correlation with blood glucose. Thus, further studies could prove their usefulness as hyperglycemia-independent biomarkers for T1D in early stages or non-fasting conditions. Interestingly, we did not detect significant changes in the levels of reduced and oxidized glutathione (GSH and GSSG), although we detected a tendency of GSH levels to decrease, and the glutathione pathway seemed to be affected in enrichment analysis (Figure 4). Decreased GSH and increased GSSG levels have been previously detected in plasma [27,47,48] and RBCs [49,50] of T2D patients. This alteration was attributed to a decrease in the activities of the enzymes involved in GSH synthesis [49], a decrease in the transport rate of oxidized glutathione (GSSG) from RBCs [51], an enhanced sorbitol pathway [52], and the increase of oxidative stress produced by the elevated generation of ROS linked to T2D [53]. Our result was coherent with the fact that a decrease of GSH in blood of T1D patients was detected less frequently [47,54]. Thus, the glutathione pathway seems to be less affected by T1D. We also did not detect significant differences in 2,3-Bisphosphoglyceric acid (2,3-BPG) levels, which have been described for T2D [25]. 2,3-BPG regulates oxygen absorption by hemoglobin, and its alteration can be directly associated to hypoxia and related to oxidative stress [55]. Thus, our study indicates lower hypoxia in organs and tissues related to RBCs and reactive oxygen species (ROS) levels in RBCs in T1D patients, in comparison with J. Clin. Med. 2023, 12, 556 11 of 14 T2D. This result is in accordance with the lower impact on glutathione metabolism that we detected for T1D patients. Regarding the lipidomic profile of RBCs, we found very significant differences between T1D patients and healthy controls. This result was in coherence with previous works that detected alterations in lipid metabolism associated with T1D and T2D [56–58]. Specifically, we detected a very clear decrease in lipid chain length in combination with an increase in lipid methyl groups, which could be an indication for lipid peroxidation, as was observed previously in peripheral blood mononuclear cells [59]. On the other hand, the increase in methyl lipid groups may also be related to the higher lipid accumulation of T1D patients due to the alterations of insulin release from pancreatic beta cells that inhibit the glucose uptake in cells [60]. In addition, we detected changes in the phospholipid phosphatidyl- choline, as well as for the related metabolites phosphocholine and glycerophosphocholine (GPC), which are all important constituents of the cell membrane. As mentioned, diabetes negatively affects the membrane and the shape of RBCs [19,25,56], which is in coherence with these results. This alteration also was shown in RBCs of T2D patients [25]. Our study is only a first pilot study, and has several limitations that will have to be overcome in additional studies. For instance, although our study groups were sex-matched to compensate for metabolic differences between men and women, the low number of samples did not allow us to perform sex-specific studies to identify specific alterations for each sex group. Additionally, the effect of many other factors, such as age, comorbidities, or treatments, should be studied. Nevertheless, we show robust differences in the metabolic and lipidomic profile of RBCs from T1D patients that provide a valuable starting point for the identification of biomarkers for an earlier diagnosis or monitoring of patients with T1D. Several of the detected metabolic changes are similar to alterations detected previously in the RBCs of T2D patients. Further studies are needed to confirm if these alterations are also related to the presence of IR in T1D patients. Supplementary Materials: The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/jcm12020556/s1. Author Contributions: Conceptualization, M.P.-S., R.S. and J.R.H.; methodology, M.P.-S., R.L.-D., A.C., C.A.-G., C.H. and J.R.H.; software, M.P.-S., R.L.-D. and J.R.H.; formal analysis, M.P.-S., A.C. and J.R.H.; investigation, M.P.-S., R.L.-D., A.C., C.A.-G., C.H., R.S. and J.R.H. resources, M.P.-S. and J.R.H.; writing—manuscript preparation, M.P.-S., R.S. and J.R.H.; writing—original draft preparation, visualization, M.P.-S., R.S. and J.R.H.; and supervision, M.P.-S., R.S. and J.R.H. All authors have read and agreed to the published version of the manuscript. Funding: All consumables: and part of laboratory equipment and human resources of this work was supported by the Carlos III Health Institute and the European Regional Development Fund (PI16/02064 and PI20/01588) and the Agency for Management of University and Research Grants (AGAUR) of Catalonia (2017SGR1303). Nuclear Magnetic Resonance Spectrometers employed for metabolite and lipid quantification have been funded by Generalitat Valenciana and co-financed with ERDF funds (OP ERDF of Comunitat Valenciana 2014–2020). Institutional Review Board Statement: The study was performed with patients from the outpatient clinic of the Endocrinology Department of Vall d´Hebron University Hospital, approved by the local Ethics Committee (PR(AG)234/2015) and carried out in accordance with the Declaration of Helsinki. Informed Consent Statement: Informed consent was obtained from all subjects involved in the study. Data Availability Statement: The datasets generated and/or analysed during the current study are available in the Zenodo repository, “https://doi.org/10.5281/zenodo.6597902 (accessed on 2 September 2022)”. Acknowledgments: We thank the NMR facility of the CIPF for the availability of acquisition time. We thank Bernabé Palacios Naranjo for technical help with the sample preparation. Conflicts of Interest: The authors declare that they have no competing interests. J. Clin. Med. 2023, 12, 556 12 of 14 References 1. United Nations. Modalities, Format and Organization of the High-Level Meeting of the General Assembly on the Prevention and Control of Non-Communicable Diseases. Resolut. A/RES/65/238 2011, 11759, 1–4. 2. Skyler, J.S.; Bakris, G.L.; Bonifacio, E.; Darsow, T.; Eckel, R.H.; Groop, L.; Groop, P.H.; Handelsman, Y.; Insel, R.A.; Mathieu, C.; et al. Differentiation of Diabetes by Pathophysiology, Natural History, and Prognosis. Diabetes 2017, 66, 241–255. [CrossRef] [PubMed] 3. Care, D.; Suppl, S.S. Classification and Diagnosis of Diabetes: Standards of Medical Care in Diabetes-2020. Diabetes Care 2020, 43, S14–S31. [CrossRef] 4. Soriguer, F.; Goday, A.; Bosch-Comas, A.; Bordiú, E.; Calle-Pascual, A.; Carmena, R.; Casamitjana, R.; Castaño, L.; Castell, C.; Catalá, M.; et al. Prevalence of Diabetes Mellitus and Impaired Glucose Regulation in Spain: The Di@bet.Es Study. Diabetologia 2012, 55, 88–93. [CrossRef] [PubMed] 5. Patterson, C.C.; Karuranga, S.; Salpea, P.; Saeedi, P.; Dahlquist, G.; Soltesz, G.; Ogle, G.D. Worldwide Estimates of Incidence, Prevalence and Mortality of Type 1 Diabetes in Children and Adolescents: Results from the International Diabetes Federation Diabetes Atlas, 9th Edition. Diabetes Res. Clin. Pract. 2019, 157, 107842. [CrossRef] 6. Mobasseri, M.; Shirmohammadi, M.; Amiri, T.; Vahed, N.; Fard, H.H.; Ghojazadeh, M. Prevalence and Incidence of Type 1 Diabetes in the World: A Systematic Review and Meta-Analysis. Heal. Promot. Perspect. 2020, 10, 98–115. [CrossRef] [PubMed] 7. Pippitt, K.; Li, M.; City, S.L.; Gurgle, U.H.E.; College, U.; City, S.L. Diabetes Mellitus: Screening and Diagnosis. Am. Fam. Phys. 2016, 93, 103–109. 8. Maahs, D.M.; West, N.A.; Lawrence, J.M.; Mayer-Davis, E.J. Epidemiology of Type 1 Diabetes. Endocrinol. Metab. Clin. N. Am. 2010, 39, 481–497. [CrossRef] 9. Roglic, G.; Norris, S.L. Medicines for Treatment Intensification in Type 2 Diabetes and Type of Insulin in Type 1 and Type 2 Diabetes in Low-Resource Settings: Synopsis of the World Health Organization Guidelines on Second- and Third-Line Medicines and Type of Insulin for the Con. Ann. Intern. Med. 2018, 169, 394–397. [CrossRef] 10. Wolosowicz, M.; Lukaszuk, B.; Chabowski, A. The Causes of Insulin Resistance in Type 1 Diabetes Mellitus: Is There a Place for Quaternary Prevention? Int. J. Environ. Res. Public Health 2020, 17, 8651. [CrossRef] 11. Gregory, J.M.; Smith, T.J.; Slaughter, J.C.; Mason, H.R.; Hughey, C.C.; Smith, M.S.; Kandasamy, B.; Greeley, S.A.W.; Philipson, L.H.; Naylor, R.N.; et al. Iatrogenic Hyperinsulinemia, Not Hyperglycemia, Drives Insulin Resistance in Type 1 Diabetes as Revealed by Comparison with GCK-MODY (MODY2). Diabetes 2019, 68, 1565–1576. [CrossRef] 12. Rains, J.L.; Jain, S.K. Oxidative Stress, Insulin Signaling, and Diabetes. Free Radic. Biol. Med. 2011, 50, 567–575. [CrossRef] [PubMed] 13. Citro, A.; Campo, F.; Dugnani, E.; Piemonti, L. Innate Immunity Mediated Inflammation and Beta Cell Function: Neighbors or Enemies? Front. Endocrinol. 2021, 11, 606332. [CrossRef] 14. Eizirik, D.L.; Colli, M.L.; Ortis, F. The Role of Inflammation in Insulitis and B-Cell Loss in Type 1 Diabetes. Nat. Rev. Endocrinol. 2009, 5, 219–226. [CrossRef] 15. Tsalamandris, S.; Antonopoulos, A.S.; Oikonomou, E.; Papamikroulis, G.A.; Vogiatzi, G.; Papaioannou, S.; Deftereos, S.; Tousoulis, D. The Role of Inflammation in Diabetes: Current Concepts and Future Perspectives. Eur. Cardiol. Rev. 2019, 14, 50–59. [CrossRef] [PubMed] 16. American Diabetes Association. ADA-Standards-of-Medical-Care-Diabetes-Care-2021; Diabetes Care; American Diabetes Association: Arlington County, VA, USA, 2021; Volume 44. 17. Hare, M.J.L.; Shaw, J.E.; Zimmet, P.Z. Current Controversies in the Use of Haemoglobin A 1c. J. Intern. Med. 2012, 271, 227–236. [CrossRef] 18. Alsalhi, M.S.; Devanesan, S.; Alzahrani, K.E.; Alshebly, M.; Al-Qahtani, F.; Farhat, K.; Masilamani, V. Impact of Diabetes Mellitus on Human Erythrocytes: Atomic Force Microscopy and Spectral Investigations. Int. J. Environ. Res. Public Health 2018, 15, 2368. [CrossRef] 19. Zilberman-Kravits, D.; Harman-Boehm, I.; Shuster, T.; Meyerstein, N. Increased Red Cell Aggregation Is Correlated with HbA1C and Lipid Levels in Type 1 but Not Type 2 Diabetes. Clin. Hemorheol. Microcirc. 2006, 35, 463–471. 20. Piagnerelli, M.; Zouaoui Boudjeltia, K.; Brohee, D.; Vereerstraeten, A.; Piro, P.; Vincent, J.L.; Vanhaeverbeek, M. Assessment of Erythrocyte Shape by Flow Cytometry Techniques. J. Clin. Pathol. 2007, 60, 549–554. [CrossRef] [PubMed] 21. Sabo, A.; Jakovljević, V.; Stanulović, M.; Lepsanović, L.; Pejin, D. Red Blood Cell Deformability in Diabetes Mellitus: Effect of Phytomenadione. Int. J. Clin. Pharmacol. Ther. Toxicol. 1993, 31, 1–5. [PubMed] 22. Labrouche, S.; Freyburger, G.; Gin, H.; Boisseau, M.R.; Cassagne, C. Changes in Phospholipid Composition of Blood Cell Membranes (Erythrocyte, Platelet, and Polymorphonuclear) in Different Types of Diabetes - Clinical and Biological Correlations. Metabolism 1996, 45, 57–62. [CrossRef] [PubMed] 23. Dincer, Y.; Akcay, T.; Alademir, Z.; Ilkova, H. Effect of Oxidative Stress on Glutathione Pathway in Red Blood Cells from Patients with Insulin-Dependent Diabetes Mellitus. Metabolism 2002, 51, 1360–1362. [CrossRef] [PubMed] 24. Petruzzi, E.; Orlando, C.; Pinzani, P.; Sestini, R.; Del Rosso, A.; Dini, G.; Tanganelli, E.; Buggiani, A.; Pazzagli, M. Adenosine Triphosphate Release by Osmotic Shock and Hemoglobin A1C in Diabetic Subjects’ Erythrocytes. Metabolism 1994, 43, 435–440. [CrossRef] [PubMed] J. Clin. Med. 2023, 12, 556 13 of 14 25. Palomino-Schätzlein, M.; Lamas-Domingo, R.; Ciudin, A.; Gutiérrez-Carcedo, P.; Marés, R.; Aparicio-Gómez, C.; Hernández, C.; Simó, R.; Herance, J.R. A Translational in Vivo and in Vitro Metabolomic Study Reveals Altered Metabolic Pathways in Red Blood Cells of Type 2 Diabetes. J. Clin. Med. 2020, 9, 1619. [CrossRef] 26. Domingo-Ortí, I.; Lamas-Domingo, R.; Ciudin, A.; Hernández, C.; Herance, J.R.; Palomino-Schätzlein, M.; Pineda-Lucena, A. Metabolic Footprint of Aging and Obesity in Red Blood Cells. Aging (Albany. NY). 2021, 13, 4850–4880. [CrossRef] 27. Foster, D.; Karloff, H.; Shirley, K.E. How Well Does the Standard Body Mass Index or Variations with a Different Exponent Predict Human Lifespan? Obesity 2016, 24, 469–475. [CrossRef] 28. Palomino-Schätzlein, M.; Simó, R.; Hernández, C.; Ciudin, A.; Mateos-Gregorio, P.; Hernández-Mijares, A.; Pineda-Lucena, A.; Herance, J.R. Metabolic Fingerprint of Insulin Resistance in Human Polymorphonuclear Leucocytes. PLoS ONE 2018, 13, e0199351. [CrossRef] 29. Ulrich, E.L.; Akutsu, H.; Doreleijers, J.F.; Harano, Y.; Ioannidis, Y.E.; Lin, J.; Livny, M.; Mading, S.; Maziuk, D.; Miller, Z.; et al. BioMagResBank. Nucleic Acids Res. 2008, 36. [CrossRef] 30. Wishart, D.S.; Guo, A.C.; Oler, E.; Wang, F.; Anjum, A.; Peters, H.; Dizon, R.; Sayeeda, Z.; Tian, S.; Lee, B.L.; et al. HMDB 5.0: The Human Metabolome Database for 2022. Nucleic Acids Res. 2022, 50, D622–D631. [CrossRef] 31. Herance, J.R.; Ciudin, A.; Lamas-Domingo, R.; Aparicio-Gómez, C.; Hernández, C.; Simó, R.; Palomino-Schätzlein, M. Normalized NMR Integration Values from the Metabolomic and Lipidomic Analysis of Red Blood Cell Extract from Type 1 Diabetes Patients and Matched Healthy Subjects. Zenodo 2022. [CrossRef] 32. Wiklund, S.; Johansson, E.; Sjöström, L.; Mellerowicz, E.J.; Edlund, U.; Shockcor, J.P.; Gottfries, J.; Moritz, T.; Trygg, J. Visualization of GC/TOF-MS-Based Metabolomics Data for Identification of Biochemically Interesting Compounds Using OPLS Class Models. Anal. Chem. 2008, 80, 115–122. [CrossRef] [PubMed] 33. Eriksson, L.; Trygg, J.; Svante, W. CV-ANOVA for Significance Testing of PLS and OPLS®Models. J. Chemom. 2008, 22, 11–12. [CrossRef] 34. Pang, Z.; Chong, J.; Zhou, G.; De Lima Morais, D.A.; Chang, L.; Barrette, M.; Gauthier, C.; Jacques, P.É.; Li, S.; Xia, J. MetaboAnalyst 5.0: Narrowing the Gap between Raw Spectra and Functional Insights. Nucleic Acids Res. 2021, 49, W388–W396. [CrossRef] [PubMed] 35. Schwender, H.; Ickstadt, K. Empirical Bayes Analysis of Single Nucleotide Polymorphisms. BMC Bioinform. 2008, 9, 144. [CrossRef] [PubMed] 36. Liaw, A.; Wiener, M. Classication and Regression by RandomForest. R News 2002, 2, 18–22. 37. Guo, X.; Li, H.; Xu, H.; Woo, S.; Dong, H.; Lu, F.; Lange, A.J.; Wu, C. Glycolysis in the Control of Blood Glucose Homeostasis. Acta Pharm. Sin. B 2012, 2, 358–367. [CrossRef] 38. Akanji, A.O.; Humphreys, S.; Thursfield, V.; Hockaday, T.D.R. The Relationship of Plasma Acetate with Glucose and Other Blood Intermediary Metabolites in Non-Diabetic and Diabetic Subjects. Clin. Chim. Acta 1989, 185, 25–34. [CrossRef] 39. Ibrahim, S.I.; Amodu, A.D.; Ene-Ojo, A.S.; Ismaila, U.A.; Fakhruddeen, M. Effect of Hyperglycemia on Erythrocyte Carbonic Anhydrase and Lactic Acid in Type II Diabetic Subjects. J. Diabetes Mellit. 2016, 06, 158–165. [CrossRef] 40. Davis, J. What Are Leucine and Isoleucine? News Med. Life Sci. 2019, 1–4. 41. Mochida, T.; Tanaka, T.; Shiraki, Y.; Tajiri, H.; Matsumoto, S.; Shimbo, K.; Ando, T.; Nakamura, K.; Okamoto, M.; Endo, F. Time-Dependent Changes in the Plasma Amino Acid Concentration in Diabetes Mellitus. Mol. Genet. Metab. 2011, 103, 406–409. [CrossRef] 42. Holeček, M. Why Are Branched-Chain Amino Acids Increased in Starvation and Diabetes? Nutrients 2020, 12, 3087. [CrossRef] 43. Hameed, A.; Mojsak, P.; Buczynska, A.; Suleria, H.A.R.; Kretowski, A.; Ciborowski, M. Altered Metabolome of Lipids and Amino Acids Species: A Source of Early Signature Biomarkers of T2DM. J. Clin. Med. 2020, 9, 2257. [CrossRef] 44. Seibert, R.; Abbasi, F.; Hantash, F.M.; Caulfield, M.P.; Reaven, G.; Kim, S.H. Relationship between Insulin Resistance and Amino Acids in Women and Men. Physiol. Soc. 2015, 3, e12392. [CrossRef] [PubMed] 45. Post, A.; Groothof, D.; Schutten, J.C.; Flores-Guerrero, J.L.; Swarte, J.C.; Douwes, R.M.; Kema, I.P.; de Boer, R.A.; Garcia, E.; Connelly, M.A.; et al. Plasma Creatine and Incident Type 2 Diabetes in a General Population-Based Cohort: The PREVEND Study. Clin. Endocrinol. 2021, 94, 563–574. [CrossRef] 46. Tullson, P.C.; Bangsbo, J.; Hellsten, Y.; Richter, E.A. IMP Metabolism in Human Skeletal Muscle after Exhaustive Exercise. J. Appl. Physiol. 1995, 78, 146–152. [CrossRef] 47. Darmaun, D.; Smith, S.D.; Sweeten, S.; Sager, B.K.; Welch, S.; Mauras, N. Evidence for Accelerated Rates of Glutathione Utilization and Glutathione Depletion in Adolescents with Poorly Controlled Type 1 Diabetes. Diabetes 2005, 54, 190–196. [CrossRef] [PubMed] 48. Lorenzo-Medina, M.; Uranga, B.; Rus, A.; Martínez, R.; Puertas, C.; Blanco, M.D.; Casís, E.; Corcoy, R. Sex and Age Affect Agreement between Fasting Plasma Glucose and Glycosylated Hemoglobin for Diagnosis of Dysglycemia. Endocrinol. Diabetes y Nutr. 2017, 64, 345–354. [CrossRef] [PubMed] 49. Ahmed, R.G. The Physiological and Biochemical Effects of Diabetes on the Balance between Oxidative Stress and Antioxidant Defense System. Med. J. Islam. World Acad. Sci. 2005, 15, 32–42. 50. Sender, R.; Fuchs, S.; Milo, R. Revised Estimates for the Number of Human and Bacteria Cells in the Body. PLoS Biol. 2016, 14, e1002533. [CrossRef] J. Clin. Med. 2023, 12, 556 14 of 14 51. Murakami, K.; Ohtsuka, Y.; Fujiwara, Y.; Shimada, M. Impairment of Glutathione Metabolism in Erythrocytes With Diabetes Mellitus. Natl. Libr. Med. 1989, 38, 753–758. [CrossRef] 52. Ciuchi, E.; Odetti, P.; Prando, R. Relationship Between Glutathione and Sorbitol Concentrations in Erythrocytes From Diabetic Patients. Metabolism 1996, 45, 611–613. [CrossRef] 53. Sung, K.C.; Lee, M.Y.; Lee, J.Y.; Lee, S.H.; Kim, S.H.; Kim, S.H. Utility of ALT Concentration in Men and Women with Nonalcoholic Fatty Liver Disease: Cohort Study. J. Clin. Med. 2019, 8, 445. [CrossRef] 54. Sakhi, A.K.; Berg, T.J. Reduced Glutathione Concentrations Are Not Decreased in Red Blood Cells of Patients with Long Term Type 1-Diabetes. Scand. J. Clin. Lab. Investig. 2011, 71, 108–111. [CrossRef] 55. Płoszczyca, K.; Czuba, M.; Chalimoniuk, M.; Gajda, R.; Baranowski, M. Red Blood Cell 2,3-Diphosphoglycerate Decreases in Response to a 30 Km Time Trial Under Hypoxia in Cyclists. Front. Physiol. 2021, 12. [CrossRef] 56. Bucala, R.; Makita, Z.; Vega, G.; Grundy, S.; Koschinsky, T.; Cerami, A.; Vlassara, H. Modification of Low Density Lipoprotein by Advanced Glycation End Products Contributes to the Dyslipidemia of Diabetes and Renal Insufficiency. Proc. Natl. Acad. Sci. USA 1994, 91, 9441–9445. [CrossRef] [PubMed] 57. Murfitt, S.A.; Zaccone, P.; Wang, X.; Acharjee, A.; Sawyer, Y.; Koulman, A.; Roberts, L.D.; Cooke, A.; Griffin, J.L. Metabolomics and Lipidomics Study of Mouse Models of Type 1 Diabetes Highlights Divergent Metabolism in Purine and Tryptophan Metabolism Prior to Disease Onset. J. Proteome Res. 2018, 17, 946–960. [CrossRef] 58. Vergès, B. Lipid Disorders in Type 1 Diabetes. Diabetes Metab. 2009, 35, 353–360. [CrossRef] [PubMed] 59. Jones, A.R.; Coleman, E.L.; Husni, N.R.; Deeney, J.T.; Raval, F.; Steenkamp, D.; Dooms, H.; Nikolajczyk, B.S.; Corkey, B.E. Type 1 Diabetes Alters Lipid Handling and Metabolism in Human Fibroblasts and Peripheral Blood Mononuclear Cells. PLoS ONE 2017, 12, e0188474. [CrossRef] 60. Ritchie, R.H.; Zerenturk, E.J.; Prakoso, D.; Calkin, A.C. Lipid Metabolism and Its Implications for Type 1 Diabetes-Associated Cardiomyopathy. J. Mol. Endocrinol. 2017, 58, R225–R240. [CrossRef] [PubMed] 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.