Metabolomics analysis of type 2 diabetes remission identifies 12 metabolites with predictive capacity: a CORDIOPREV clinical trial study
BMC Medicine volume 20, Article number: 373 (2022)
Type 2 diabetes mellitus (T2DM) is one of the most widely spread diseases, affecting around 90% of the patients with diabetes. Metabolomics has proven useful in diabetes research discovering new biomarkers to assist in therapeutical studies and elucidating pathways of interest. However, this technique has not yet been applied to a cohort of patients that have remitted from T2DM.
All patients with a newly diagnosed T2DM at baseline (n = 190) were included. An untargeted metabolomics approach was employed to identify metabolic differences between individuals who remitted (RE), and those who did not (non-RE) from T2DM, during a 5-year study of dietary intervention. The biostatistical pipeline consisted of an orthogonal projection on the latent structure discriminant analysis (O-PLS DA), a generalized linear model (GLM), a receiver operating characteristic (ROC), a DeLong test, a Cox regression, and pathway analyses.
The model identified a significant increase in 12 metabolites in the non-RE group compared to the RE group. Cox proportional hazard models, calculated using these 12 metabolites, showed that patients in the high-score tercile had significantly (p-value < 0.001) higher remission probabilities (Hazard Ratio, HR, high versus low = 2.70) than those in the lowest tercile. The predictive power of these metabolites was further studied using GLMs and ROCs. The area under the curve (AUC) of the clinical variables alone is 0.61, but this increases up to 0.72 if the 12 metabolites are considered. A DeLong test shows that this difference is statistically significant (p-value = 0.01).
Our study identified 12 endogenous metabolites with the potential to predict T2DM remission following a dietary intervention. These metabolites, combined with clinical variables, can be used to provide, in clinical practice, a more precise therapy.
T2DM is a metabolic disorder widely identified by a generalized hyperglycaemia and insulin resistance . Nowadays, approximately 422 million adults are diagnosed with diabetes ; T2DM is the most prevalent form of the disease, affecting 90% of the patients (circa 380 million worldwide). Common co-morbidities associated with T2DM include cardiovascular diseases, blindness, nerve damage, and kidney failure [3,4,5]. The co-occurrence of coronary heart disease (CHD) along with T2DM markedly increases the risk of macrovascular complications and mortality. Indeed, macrovascular events represent circa 80% of all deaths in these patients . This current scenario urges us to find new approaches to diagnose and treat these patients. Metabolic profiling, or metabolomics, allows for the characterization of hundreds of compounds (i.e. metabolites) facilitating functional information of the metabolism at the time when the sample is taken [6,7,8] and has proven useful to understand the metabolism in the context of (i) diagnostic and prognostic biomarker discovery, (ii) therapy research, and (iii) pathways determination within T2DM .
Recent evidence has shown that T2DM is reversible when tackled in an early phase following two different strategies: calorie restriction and bariatric surgery. The former is linked to weight loss, gut permeability, and reduction in inflammatory and endotoxemia biomarkers . The latter leads to normalizing plasma glucose levels and a significant weight loss [10,11,12].
Here, we analyse the metabolomics modulations differences between individuals where T2DM has remitted (RE) and those who did not recover and consequently remained as diabetics (non-RE) after 5 years of dietary intervention. The dietary intervention consisted of two different types of diets: low-fat diet (LF) and Mediterranean diet (MED).
Aim and objectives
The aim is to characterize the metabolic profile of non-RE and RE patients in serum samples at the baseline to identify biomarkers of interest to assist in the diagnosis, monitoring, and treatment of the disease. These biomarkers can be particularly useful to predict who will recover from T2DM following a dietary intervention.
The specific objectives of our study are to analyse the metabolomics differences before the dietary intervention and identify biomarkers of interest to assist in the prediction of T2DM recovery.
Study design and participants
This study was developed within the framework of the CORDIOPREV (CORonary Diet Intervention with Olive oil and cardiovascular PREVention) study, registered at Clinicaltrials.gov (number NCT00924937). This study is an ongoing controlled, single-blind, and randomized trial, with 1002 CHD patients. The trial protocol and subsequent revisions were approved by the Reina Sofia University Hospital Ethics Committee, following the Helsinki Declaration and good clinical practices. All patients signed a written informed consent to participate in the study.
Patients’ recruitment took place between November 2009 and February 2012, mostly at Reina Sofia University Hospital, Córdoba, Andalusia, Spain, with contributions from other hospitals in Córdoba and Jaen areas, in Andalusia, Spain. Complete details of the study methods, rationale, inclusion criteria, cardiovascular risk factors, and baseline characteristics are found elsewhere . In brief, eligible participants, in the age range of 20 to 75 years, had established CHD with no clinical events in the previous 6 months. They all had at least a 5-year life expectancy and no other concurrent major diseases and were willing to participate in a long-term monitoring study .
In this work, 183 patients, from the CORDIOPREV study (https://www.cordioprev.es/index.php/es/) diagnosed with diabetes, underwent a dietary interventional study where participants were offered either an LF or MED diet for 5 years (Fig. 1). In our study, blood samples were taken during fasting (time 0) and 120 min after a glucose boost.
Oral glucose tolerance test
The patients underwent an OGTT at the baseline and, once a year, every year during the dietary intervention. Before the test, patients had fasted (from food/drugs) for 12 h and were asked to refrain from smoking and alcohol intake during the preceding 7 days. They were also asked to avoid strenuous physical activity, a day before the test. At 8:00 A.M., patients were admitted to the laboratory to perform the oral glucose tolerance test (OGTT) (75 g flavoured glucose load, Trutol 75; Custom Laboratories, Baltimore, MD, USA). Blood samples were taken at times corresponding to 0, 30, 60, 90, and 120 min to determine the glucose and insulin concentrations . The insulin sensitivity index (ISI) was calculated from the OGTT using the following formula: ISI = 10.000 ÷ √([fasting plasma glucose X fasting plasma insulin] × [mean glucose in OGTT × mean insulin in OGTT]) . HOMA-IR was calculated as described by Song et al. . Beta-cell function was calculated using the disposition index (DI) as follows: DI = ISI × [AUC30 min insulin/AUC30 min glucose], where AUC30 min is the area under the curve between baseline and that at 30 min of the OGTT for insulin (pmol/l) and glucose (mmol/l) measurements, calculated by the trapezoidal method . The indices used to determine tissue-specific insulin resistance (IR) were the hepatic insulin resistance index (HIRI, fasting plasma insulin × fasting plasma glucose) and the muscle insulin sensitivity index (MISI, (dG/dt)/mean of plasma insulin) . Insulinogenic Index (IGI) was calculated by measuring plasma insulin at 30 min − fasting plasma insulin (mU/L)/(plasma glucose at 30 min − fasting plasma glucose(mg/dL) .
The adipose tissue (AT) insulin resistance index (Adipo-IR) was determined according to the formula: Adipo-IR = fasting plasma NEFA (mM) × fasting plasma insulin (pmol/L), which has been found as a suitable and useful method in clinical practice to estimate AT insulin sensitivity .
Randomization and masking
The process of randomization has been reported elsewhere . Briefly, this is based on the following variables: sex (male, female), age (under and over 60 years old), and previous myocardial infarction (yes, no). Eight different groups were created to represent all the possible combinations of the above factors. Therefore, eight different blocks were created to assign the diets (bloc randomization). Dietitians were the only members of the intervention team to be aware of the dietary group of each participant.
The participants were randomized to consume two diets: the Med diet or an LF diet . The LF diet consists of < 30% total fat (< 10% saturated fat, 12–14% MUFA fat, and 6–8% PUFA fat), 15% protein, and a minimum of 55% carbohydrates. The Med diet consists of a minimum of 35% of calories as fat (22% MUFA fat, 6% PUFA fat, and < 10% saturated fat), 15% proteins, and a maximum of 50% carbohydrates . Neither energy restriction, nor physical activity was specifically encouraged. In both diets, the cholesterol content was adjusted to < 300 mg/d.
The Mediterranean and low-fat diets were designed to provide a wide variety of foods, including vegetables, fruit, cereals, potatoes, legumes, dairy products, meat, and fish. The participants in both intervention groups received the same intensive dietary counselling. The nutritionists administered personalized individual advice every 6 months. In addition, quarterly group education sessions were held with up to 20 participants per session; separate group sessions were performed every 3 months, and dietary counselling by phone was done every 2 months . At the beginning of the study, and every 6 months afterwards, each patient had a face-to-face interview with a nutritionist to complete a 137-item semi-quantitative food frequency questionnaire (validated in Spain ). The dietary evaluation was calculated by the 14-item Med Diet Adherence Screener, which was used for measuring adherence to the Med diet . Moreover, a 9-item dietary adherence screener was used to measure adherence to the LF diet guidelines. A more detailed report on dietary adherence has been published recently by our research group .
Diabetes remission criteria
Remission required the following: (i) the absence of glucose-lowering treatment and was defined by levels of HbA1c < 6·5%, (ii) a fasting plasma glucose < 126 mg/dl, and (iii) a 2-h plasma glucose in the 75 g OGTT < 200 mg/dl maintained for at least 2 years. This agrees with the American Diabetes Association (ADA) diagnosis criteria .
Plasma samples (100 μL) were immersed in bath ice and treated with 300 μL of 3:1 (v/v) methanol–acetonitrile (MeOH–ACN). The treated samples were vortexed for 2 min and subsequently cooled at − 20 °C for 3 min. Centrifugation was carried out for 15 min at 4 °C and 13,800 × g in a thermostatic centrifuge Thermo Sorvall Legend Micro 21 R from Thermo (Thermo Fisher Scientific, Bremen, Germany), and the supernatant phase was isolated. This phase was dried by evaporation and reconstituted with 60 μL of 3:1 (v/v) MeOH–ACN. All samples were processed in a 1200 Series LC system (Agilent Technologies, Waldbronn, Germany) coupled to an Agilent 6530 high-resolution QTOF mass spectrometer equipped with a dual electrospray ionization source.
LC–QTOF MS/MS analysis
A Poroshell 120 EC-C18 column (50 mm × 2.1 mm i.d., 2.7 μm particle size, from Agilent), kept at 25 °C, was used to carry out the chromatographic division. The mobile phases consisted of (A) 0.1% formic acid in deionized water and (B) 0.1% formic acid in acetonitrile. The protocol used for the elution consisted of 0–2 min, 5% B; 2–11 min and the percentage of mobile phase B was modified from 0 to 100%. The final percentage was held for 6 min. Five minutes post-run was included to equilibrate the column. The flow rate was maintained at 0.4 mL/min. The injected sample volume was 5.0 μL and the injector needle was washed 10 times with 70% methanol between injections. Therefore, the needle seat was flushed for 15 s at a flow rate of 4 mL/min, with 70% methanol, to avoid cross-contamination between samples. The autosampler was kept at 4 °C to increase sample stability. The settings of the electrospray ionization source, which was operated in negative and positive ionization modes, were as follows: capillary voltage ± 3.5 kV, Q1 voltage 130 V, N2 pressure in the nebulizer 35 psi; N2 flow rate and temperature as drying gas 10 L min–1 and 325 °C, respectively. MS/MS data were acquired in both polarities, using the centroid mode at a rate of 2.5 spectra s–1 in extended dynamic range mode (2 GHz). Accurate mass spectra in the MS scan were acquired in the m/z range 40–1100 and the MS/MS mode in the m/z range 30–1100. The instrument gave a typical resolution of 18,000 full width at half maximum (FWHM) at m/z 118.0862 and 35,000 FWHM at m/z 922.0098. The instrument was calibrated and tuned as recommended by the manufacturer. To assure the desired mass resolution, continuous internal calibration was performed during analyses by using the signals at m/z 121.0509 (protonated purine) and m/z 922.0098 [protonated hexakis(1H,1H,3H-tetrafluoropropoxy) phosphazine or HP-921] in the positive ion mode, while in the negative ion mode, ions with m/z 119.0362 (proton abstracted purine) and m/z 966.0007 (formate adduct of HP-921) were used. The collision energy was set at 20 V for the whole run. The analytical samples were injected in auto MS/MS acquisition mode to obtain fragmentation information from a maximum of two precursors selected per cycle with an exclusion window of 0.1 min after 2 consecutive selections of the same precursor.
The MassHunter Workstation software (version B7.00 Qualitative Analysis, Agilent Technologies, Santa Clara, CA, USA) was used to process all the data obtained by LC–QTOF in data-dependent acquisition MS/MS mode. Treatment of raw data files started with the extraction of potential molecular features (MFs) with the suited algorithm included in the software. For this aim, the extraction algorithm considered all ions exceeding 500 counts for both polarities with a maximum charge state of 2 for the obtained chromatograms. The count cut-off value was established considering the chromatographic background noise. Additionally, only MFs defined by two or more ions were considered, with a tolerance for the isotopic distribution of 0.0025 m/z for peak spacing tolerance, plus 7.0 ppm in mass accuracy. Only the following potential ions and adducts were considered in positive (H + , Na + , K + , NH4 +) and negative ionization (H − , HCOO + , Cl +) modes. Furthermore, a potential neutral loss by dehydration was also included to identify features corresponding to the same potential metabolite.
Identification of metabolites was supported on MS and MS/MS information by using METLIN MS and MS/MS databases (http://metlin.scripps.edu), the Human Metabolome Database (HMDB, 3.6 version), and the LIPID MAPS website ((http://www.lipidmaps.org); in all cases, the MFs obtained in the previous step were used. A database with all identified metabolites was used to perform a targeted compound extraction analysis using a tolerance window of 0.8 min and 5 ppm mass accuracy. This step was performed with Profinder Analysis (version B8.00, Agilent Technologies, Santa Clara, CA, USA). A table with the peak area of all identified compounds in the different samples injected was obtained as a result.
Metabolites showing in at least 80% of the samples were selected for further analysis. To allow predictive modelling, imputation was carried out when needed, substituting missing values by half the smallest value of the appropriate metabolite.
LC–MS data (polar and apolar) was imported into Matlab (R2015a, Mathworks UK) and analysed using the statistic toolbox and algorithms from Korrigan Toolbox version 0.1 (Korrigan Sciences Ltd, UK). In Matlab, matrices were log 10 normalized. The biostatistical pipeline for the multivariate statistical analysis considered a preliminary unsupervised principal component analysis (PCA), followed by a supervised pairwise O-PLS DA [24, 25], which identifies the specific modulations driven by the appropriate predictor (i.e. individuals who returned from T2DM versus those who did not).
To assess the predictive power of the O-PLS DA models, R2 (the explained variance) was calculated. This parameter evaluates the model maximizing variance given by the endogenous variables. The Q2, or goodness of prediction, assesses the predictive relevance of the model and is based on a matrix partition technique that ignores part of the data (in our case a seventh part each time), estimates the model parameters, and predicts the omitted parts using the estimates obtained previously. Q2 greater than 0 means the model has predictive value. In addition, the overfitting of the model (the difference between R2Y and Q2Y) was also considered, and only models with less than 50% overfitting were further considered. Model parameters and associated metabolites were reported and used for a Cox proportional hazard model, GLM, and ROC calculations in R (version 4.0.5 (2021–03-31, https://www.r-project.org) using the packages ‘caret’ and ‘pROC’. Unadjusted Cox proportional-hazard models calculated the hazard ratio (HR) of every metabolite previously identified in the O-PLS DA model within a 95% confidence interval (CI). This unadjusted Cox allowed the identification of the betas for every metabolite. This information was used to calculate the patient’s likelihood of recovering from diabetes by running a Cox analysis adjusted for sex, age, body mass index (BMI), HDL, triglycerides, and intensity of statin therapy based on the tertiles resulting from the multiplication of the betas previously obtained by the abundance of each metabolite for every patient. Finally, generalized linear models were calculated for (i) all the clinical variables (i.e. sex, age, BMI, HDL, triglycerides, and intensity statin therapy), (ii) all the metabolites, (iii) the glycated haemoglobin, (iv) the clinical variables and the 12 metabolites, and (v) the clinical variables and the glycated haemoglobin. ROC analyses were carried out for these three models, and AUC, sensitivity, specificity, accuracy, and threshold were estimated for the models. Finally, DeLong analysis was used to compare whether the AUCs of these models were or not significantly different between them.
BMI, waist circumference, body weight, glucose, glycated haemoglobin (HbA1c), insulin, HIRI, and homeostatic model assessment of insulin resistance (HOMA-IR) were statistically significantly higher at baseline in the non-RE group than in the RE group. Conversely, ISI and DI values were statistically significantly lower at baseline in the RE group than in the non-RE group (p < 0.05) (Additional file 1: Table S1).
O-PLS DA results from the comparison between individuals who remitted from T2DM and those who did not
The O-PLS DA was examined based on R2Y, Q2Y, and overfit parameters obtained in every pairwise comparison (see material and methods). The O-PLS DA analysis identified differences between RE and DM individuals during fasting (R2Y = 0.1420, Q2Y = 0.0018), but not after a glucose overdose (Fig. 2). DM individuals had higher levels of sphingosine (d18:2), docosenamide, oxo-tricosanoic acid, tetracosahexaenoic acid, ketodeoxycholic acid, stearoylcarnitine, diglyceride (33:4), creatine, tridecanoic acid, monoacylglycerol (22:6), dihydroxycholesterol, and biliverdin. These metabolites, as well as a ranking with the degree of statistically significant association of each metabolite, are presented in further detail in Additional file 2: Fig. S1.
Cox proportional hazard models
Adjusted and unadjusted Cox proportional hazard models were calculated for every metabolite previously identified in the O-PLS DA model. These metabolites were grouped in terciles for the Cox analysis. The results are shown in Table 1. The unadjusted beta values were used to calculate the patient risk score (beta × metabolite abundance) for all the patients and metabolites. This patient risk score grouped by terciles was used to carry out a new Cox model adjusted for sex, age, BMI, HDL, TG, and intensity statin therapy (see Fig. 3). Results from the Cox analysis of the patient’s risk score showed that individuals in the high tercile (luckily to recover from T2DM) had an HR of 2.70 compared to individuals in the low tercile; this difference was statistically significant (p-value = 0.003). Moreover, we also grouped the patients by ascending terciles of glycated haemoglobin, to carry out a Cox model adjusted for sex, age, BMI, HDL, TG, and intensity statin therapy (Additional file 2: Fig. S2). Results from the Cox analysis of the glycated haemoglobin showed that individuals in the low and medium tercile had an HR of 3.80 and 4.31 respectively, compared to individuals in the High tercile (both p-value = 0.001).
Generalized linear models and receiving operating characteristics
GLMs and ROCs were run for (i) clinical variables alone, (ii) the 12 metabolites of interest, (iii) glycated haemoglobin, (iv) the clinical variables and the 12 metabolites, and (v) the clinical variables and the glycated haemoglobin. AUC, sensitivity, specificity, accuracy, and threshold were calculated in all cases. The results are shown in Tables 2 and 3 and Fig. 4. The AUC for the clinical variables was 0.610, for the metabolites was 0.701, for the glycated haemoglobin 0.618, for the combination of clinical variables and the 12 metabolites 0.721, and for the combination of the clinical variables and the glycated haemoglobin 0.667. A DeLong test was carried out comparing the AUC from the GLM of the clinical variables alone with the GLM with the clinical variables and the 12 metabolites or the GLM with the clinical variables and the glycated haemoglobin (Fig. 4). There were not significant differences between the model with the clinical variables alone and the addition of the glycated haemoglobin. However, the p-value (0.01265) resulting from the DeLong test comparing the model with the clinical variables alone and the clinical variables with the 12 metabolites indicated that these two models were significantly different and showed that the addition of the metabolites significantly improved the prediction capacity of the model.
Our study identified 12 plasma metabolites by O-PLS DA differing between RE and non-RE patients at the baseline of the study; these metabolites were further used to build a score to assess regression probability. This score was significantly associated with a higher probability of T2DM remission. These 12 metabolites, together with the clinical variables previously described, significantly improved the T2DM remission prediction power of the model. However, the prediction capacity of the model with the clinical variables did not significantly improve when the glycated haemoglobin was added. In practice, this could be achieved by analysing a plasma sample from the patients and determining the concentrations of these molecules.
Lifestyle modifications, including the implementation of healthy diets, result in a beneficial effect on T2DM prevention and management . Recent studies indicated that it was possible to induce T2DM remission by weight loss with calorie restriction interventions . Previous studies showed that bariatric surgeries can revert T2DM regulating the glucose in plasma before a significant weight reduction is obtained [9,10,11].
The importance of the early identification of diabetic patients with a probability of achieving T2DM remission lies in the ability of the β cells to recover long-term functionality after T2DM diagnosis, and before an irreversible stage of β cell dysfunction . Moreover, an effective and efficient therapeutic action focused on disease remission is especially important in individuals with co-occurring acute myocardial infarction and T2DM, who have a higher risk of developing a new cardiovascular event than those without T2DM .
In terms of predicting T2DM remission, a variety of scores and variables were previously used to identify subjects with a probability of remission. Classical clinical variables have shown a reduced prediction capacity [29, 30]. Conversely, we have previously reported how the use of miRNAs or gut microbiota proxies improve the predictive power of the clinical variables alone, improving the estimation of the T2DM remission probability induced by the dietary intervention [31,32,33]. However, to the best of our knowledge, the plasma metabolite profile has not been applied to such end.
This study has shown a metabolic profile associated with T2DM remission in CHD patients. Here, TOF/LS-MS metabolomics at baseline was used to identify which newly diagnosed T2DM patients will benefit from a dietary intervention (a Med or a low-fat diet) to induce the remission of T2DM, with the main difference between the diets in the amount and type of dietary sources of fat and the amount of carbohydrates . Both diets were ethically appropriate for this profile of patients and no further energy restriction or physical activity was implemented.
This profile was characterized by a reduction in 12 metabolites across RE patients identified by O-PLS DA. This model also showed which metabolites found have the strongest influence. These metabolites ranked with a different degree of significance in the model as shown in Additional file 2: Fig. S1. These metabolites, according to the results from the DeLong test, significantly improved the predictive capacity of the clinical variables alone. These metabolites are included in several insulin-related pathways such as sphingolipid metabolism or alpha-linolenic acid and linoleic acid metabolism. Moreover, taken together, these metabolites may be linked to the lipid alterations associated with T2DM. In fact, the derangement of the lipid metabolism is a common complication in T2DM due to the inadequate functioning of key enzymes and pathways as well as the insulin resistance prevalent in these patients . In addition, the association between dyslipidaemia and atherosclerosis is well established, and the composition of lipid particles in diabetic dyslipidaemia has a stronger atherogenic impact on the disease compared to other kinds of dyslipidaemia [35, 36]. Thus, it would be also expected the relationship of these metabolites with lipid alterations considering that the individuals included in this study were CHD patients. Moreover, these disruptions in the metabolite profile in non-RE patients seem to be additional to the dyslipidaemia associated with CHD, which includes hypertriglyceridemia, hypercholesterolemia, and elevated LDL cholesterol . This present study brings insights into the metabolomics modulations occurring during dyslipidaemia.
The O-PLS DA model showed that out of the 12 metabolites identified, tetracosahexaenoic acid (THA), oxo-tricosanoic acid, dihydrocholesterol, and tridecanoic acid were the most discriminant between RE and non-RE; hence, high levels were observed in non-RE patients. In line with this, it is worth mentioning that THA and tridecanoic acid are rare fatty acids. Current literature is contradictory about the beneficial or detrimental effects of rare fatty acids in CHD or T2DM patients [38,39,40]. In our study, higher levels of THA and tridecanoic acid were observed in non-RE patients. This could be because our population comprises patients with established CHD in which may exist certain alterations in their lipid species profile patients due to the cardiovascular disease.
Dihydroxycholesterol is a metabolite involved in the primary bile acid biosynthesis pathway and derived from cholesterol. This metabolite was identified as an intermediate in C21-Steroid hormone metabolism . Previous studies have shown that plasma concentrations of cholesterol oxidation (ChOx) products (including different forms of hydroxycholesterol) are elevated in DM1 and DM2 patients compared to age-matched subjects without diabetes . In our study, significantly higher levels of dihydroxycholesterol were identified across those individuals who did not return from T2DM with the dietary intervention, which suggests that the lipid metabolism of those individuals is more impaired in non-RE than in RE patients. The relationship between cholesterol oxidized metabolites and lipid metabolism lies in the fact that high levels of glucose seem to promote lipidic accumulation via the DNA CpG methylation in the DNA methyltransferase-1 (DNMT1) promoter region, where cholesterol oxidized metabolites act as mediators . The high levels of glucose increase a nuclear form of oxidized cholesterol (25-hydroxycholesterol) which activates DNMT1 regulating the expression of several genes implicated in the intracellular lipid metabolism in the liver. This process results in the hypermethylation of genes directly implicated in carbohydrate and lipid metabolism of PI3K, cAMP, insulin, insulin secretion, and diabetic and non-alcoholic fatty liver disease signalling pathways. This suggests that the high hydroxycholesterol plasma levels found in our study, in the non-RE group, might be altering the expression of these genes towards a deleterious gene expression profile. This epigenetic regulation of hepatic cell metabolism seems to also have relevance in non-alcoholic, fatty liver disease, and metabolic syndrome [44, 45].
Hepatic lipogenesis can be suppressed by downregulating the gene SREBP1c (trigged by n3 supplementation), implicated in fatty acid biosynthesis. These inhibitory effects in hepatic lipogenesis are associated with a reduction in plasma levels of THA . In our study, we have observed high levels of THA amongst the non-RE group, which suggests a lack of inhibitory feedback in this pathway. This could explain the metabolic resilience to respond to the dietary intervention in these patients, which could be due to the methylation of the promoter region DNMT1. Previous studies from Nagao and colleagues (2003) showed that the depletion in THA was associated with omega 3-PUFAs supplementation. This suggests that the omega 3-PUFAs ingested in the dietary intervention are effective only in patients with low THA levels at baseline. Therefore, the high levels of THA amongst non-RE patients could predispose them towards a higher degree of metabolic impedance, preventing them to respond to the healthy diet administrated which contains 6–8% of total calories as PUFAs. However, this hypothesis needs further evidence.
We also found elevated tridecanoic acid levels in the DM group, a metabolite previously associated with hypoglycaemia, fatty liver, and cardiomyopathies, that are implicated in functions such as oxidation, cell death, and insulin resistance . The latter is especially important as high tridecanoic aid levels were found in the non-RE group, which included the patients who did not reduce their insulin resistance after the ingestion of the healthy diet intervention.
Predicting which patients with CHD can recover from T2DM is crucial since patients with co-occurring CHD and T2DM have a considerably higher risk of developing a new cardiovascular event than those without T2DM. Furthermore, some of the metabolites might be linked to other physiological processes that remain unclear and future research in this field is needed.
In summary, our study provides new plasma biomarkers to predict, in combination with clinical variables, the dietary remission capacity of T2DM patients with co-occurring CHD. Indeed, the addition of the 12 metabolites identified significantly improved the prediction power of the clinical variables alone.
It is also important to mention the limitations of this study. Firstly, this research is based on a long-term, well-controlled dietary intervention, which despite ensuring the quality of the study, may not reflect the level of compliance in a free-living population. The second limitation is that the remission of T2DM was not the primary endpoint of the CORDIOPREV trial, although it was a secondary objective of this study. However, there are no reasons to believe that this randomization would not have worked in such a large subset of participants, taking into account that the baseline characteristics in the groups of patients analysed in the current study according to the diet were similar.
This study showed the association of a specific metabolic profile in plasma with T2DM remission in patients with CHD in a dietary intervention. These metabolites, combined with clinical variables, could be used to provide in clinical practice more precise therapeutical advice. This new approach will allow the possibility to discern between newly diagnosed T2DM luckily to remit from diabetes following a dietary intervention, from those who will need a more exhaustive treatment such us anti-glycaemic drugs. That would improve the management of these patients and represent a personalized medicine approach. Moreover, our results suggest that lipid metabolism is implicated in the probability to remit from T2DM. The detailed pathway that allows for high glucose (or high dietary glucose) to produce lipids accumulation is not fully understood yet, but early downregulation of liver lipogenesis seems determinant for the metabolism to recover from T2DM. Further investigations should interrogate the potential methylation of the promoter region DNMT1, to unravel whether epigenetic changes may influence the capacity to return from T2DM.
Availability of data and materials
The data that support the findings of this study are available from the corresponding author upon reasonable request.
American Diabetes Association
Adipose tissue insulin resistance index
Area under the curve
Body mass index
Coronary heart disease
Generalized linear model
Hepatic insulin resistance index
Human Metabolome Database
Homeostatic model assessment of insulin resistance
Insulin sensitivity index
Muscle insulin sensitivity index
Patients who remained as diabetics
- O-PLS DA:
Orthogonal Projection on the Latent Structure Discriminant Analysis
Oral glucose tolerance test
Principal component analysis
Patients who remitted from T2DM
Receiver operating characteristic
Type 2 diabetes mellitus
DeFronzo RA. Insulin resistance, lipotoxicity, type 2 diabetes and atherosclerosis: the missing links. The claude bernard lecture 2009. Diabetologia. 2010;53(7):1270–87.
WHO. Diabetes, Health Topics. http://www.who.int/diabetes/en/. Accessed 7 Feb 2022.
Amin AP, Spertus JA, Reid KJ, Lan X, Buchanan DM, Decker C, Masoudi FA. The prognostic importance of worsening renal function during an acute myocardial infarction on long-term mortality. Am Heart J. 2010;160(6):1065–71.
Anavekar NS, McMurray JJ, Velazquez EJ, Solomon SD, Kober L, Rouleau JL, White HD, Nordlander R, Maggioni A, Dickstein K, et al. Relation between renal dysfunction and cardiovascular outcomes after myocardial infarction. N Engl J Med. 2004;351(13):1285–95.
Trautner C, Icks A, Haastert B, Plum F, Berger M. Incidence of blindness in relation to diabetes. A population-based study. Diabetes Care. 1997;20(7):1147–53.
Nicholson JK, Lindon JC, Holmes E. ‘Metabonomics’: understanding the metabolic responses of living systems to pathophysiological stimuli via multivariate statistical analysis of biological NMR spectroscopic data. Xenobiotica. 1999;29(11):1181–9.
Pauling L, Robinson AB, Teranishi R, Cary P. Quantitative analysis of urine vapor and breath by gas-liquid partition chromatography. Proc Natl Acad Sci U S A. 1971;68(10):2374–6.
Suhre K. Metabolic profiling in diabetes. J Endocrinol. 2014;221(3):R75-85.
Taylor R, Al-Mrabeh A, Zhyzhneuskaya S, Peters C, Barnes AC, Aribisala BS, Hollingsworth KG, Mathers JC, Sattar N, Lean MEJ. Remission of human type 2 diabetes requires decrease in liver and pancreas fat content but is dependent upon capacity for beta cell recovery. Cell Metab. 2018;28(4):547-556 e543.
Clemente-Postigo M, Roca-Rodriguez Mdel M, Camargo A, Ocana-Wilhelmi L, Cardona F, Tinahones FJ. Lipopolysaccharide and lipopolysaccharide-binding protein levels and their relationship to early metabolic improvement after bariatric surgery. Surg Obes Relat Dis. 2015;11(4):933–9.
Kirwan JP, Munzberg H, Berthoud HR. Mechanisms responsible for metabolic improvements of bariatric surgeries. Diabetes. 2018;67(6):1043–4.
Zhang L, Xue X, Zhai R, Yang X, Li H, Zhao L, Zhang C. Timing of calorie restriction in mice impacts host metabolic phenotype with correlative changes in gut microbiota. mSystems. 2019;4(6):e00348-19.
Delgado-Lista J, Perez-Martinez P, Garcia-Rios A, Alcala-Diaz JF, Perez-Caballero AI, Gomez-Delgado F, Fuentes F, Quintana-Navarro G, Lopez-Segura F, Ortiz-Morales AM, et al. CORonary Diet Intervention with Olive oil and cardiovascular PREVention study (the CORDIOPREV study): rationale, methods, and baseline characteristics: a clinical trial comparing the efficacy of a Mediterranean diet rich in olive oil versus a low-fat diet on cardiovascular disease in coronary patients. Am Heart J. 2016;177:42–50.
Matsuda M, DeFronzo RA. Insulin sensitivity indices obtained from oral glucose tolerance testing: comparison with the euglycemic insulin clamp. Diabetes Care. 1999;22(9):1462–70.
Song Y, Manson JE, Tinker L, Howard BV, Kuller LH, Nathan L, Rifai N, Liu S. Insulin sensitivity and insulin secretion determined by homeostasis model assessment and risk of diabetes in a multiethnic cohort of women: the Women’s Health Initiative Observational Study. Diabetes Care. 2007;30(7):1747–52.
Tang W, Fu Q, Zhang Q, Sun M, Gao Y, Liu X, Qian L, Shan S, Yang T. The association between serum uric acid and residual beta -cell function in type 2 diabetes. J Diabetes Res. 2014;2014:709691.
Abdul-Ghani MA, Matsuda M, Balas B, DeFronzo RA. Muscle and liver insulin resistance indexes derived from the oral glucose tolerance test. Diabetes Care. 2007;30(1):89–94.
Hanson RL, Pratley RE, Bogardus C, Narayan KM, Roumain JM, Imperatore G, Fagot-Campagna A, Pettitt DJ, Bennett PH, Knowler WC. Evaluation of simple indices of insulin sensitivity and insulin secretion for use in epidemiologic studies. Am J Epidemiol. 2000;151(2):190–8.
Ter Horst KW, van Galen KA, Gilijamse PW, Hartstra AV, de Groot PF, van der Valk FM, Ackermans MT, Nieuwdorp M, Romijn JA, Serlie MJ. Methods for quantifying adipose tissue insulin resistance in overweight/obese humans. Int J Obes (Lond). 2017;41(8):1288–94.
Quintana-Navarro GM, Alcala-Diaz JF, Lopez-Moreno J, Perez-Corral I, Leon-Acuna A, Torres-Pena JD, Rangel-Zuniga OA, Arenas de Larriva AP, Corina A, Camargo A, et al. Long-term dietary adherence and changes in dietary intake in coronary patients after intervention with a Mediterranean diet or a low-fat diet: the CORDIOPREV randomized trial. Eur J Nutr. 2020;59(5):2099–110.
Fernandez-Ballart JD, Pinol JL, Zazpe I, Corella D, Carrasco P, Toledo E, Perez-Bauer M, Martinez-Gonzalez MA, Salas-Salvado J, Martin-Moreno JM. Relative validity of a semi-quantitative food-frequency questionnaire in an elderly Mediterranean population of Spain. Br J Nutr. 2010;103(12):1808–16.
Martinez-Gonzalez MA, Fernandez-Jarne E, Serrano-Martinez M, Wright M, Gomez-Gracia E. Development of a short dietary intake questionnaire for the quantitative estimation of adherence to a cardioprotective Mediterranean diet. Eur J Clin Nutr. 2004;58(11):1550–2.
American Diabetes Association Professional Practice C. 2 Classification and diagnosis of diabetes: standards of medical care in diabetes-2022. Diabetes Care. 2022;45(Suppl 1):S17–38.
Bylesjo M, Rantalainen M, Nicholson JK, Holmes E, Trygg J. K-OPLS package: kernel-based orthogonal projections to latent structures for prediction and interpretation in feature space. BMC Bioinformatics. 2008;9:106.
Cloarec O, Dumas ME, Craig A, Barton RH, Trygg J, Hudson J, Blancher C, Gauguier D, Lindon JC, Holmes E, et al. Statistical total correlation spectroscopy: an exploratory approach for latent biomarker identification from metabolic 1H NMR data sets. Anal Chem. 2005;77(5):1282–9.
Lifestyle Management. <em>Standards of Medical Care in Diabetes—2018</em>. Diabetes Care. 2018;41(Supplement 1):S38.
Taylor R, Al-Mrabeh A, Zhyzhneuskaya S, Peters C, Barnes AC, Aribisala BS, Hollingsworth KG, Mathers JC, Sattar N, Lean MEJ. Remission of human type 2 diabetes requires decrease in liver and pancreas fat content but is dependent upon capacity for β cell recovery. Cell Metab. 2018;28(4):547-556.e543.
Martin-Timon I, Sevillano-Collantes C, Segura-Galindo A, Del Canizo-Gomez FJ. Type 2 diabetes and cardiovascular disease: have all risk factors the same strength? World J Diabetes. 2014;5(4):444–70.
Lee WJ, Hur KY, Lakadawala M, Kasama K, Wong SK, Chen SC, Lee YC, Ser KH. Predicting success of metabolic surgery: age, body mass index, C-peptide, and duration score. Surg Obes Relat Dis. 2013;9(3):379–84.
Still CD, Wood GC, Benotti P, Petrick AT, Gabrielsen J, Strodel WE, Ibele A, Seiler J, Irving BA, Celaya MP, et al. Preoperative prediction of type 2 diabetes remission after Roux-en-Y gastric bypass surgery: a retrospective cohort study. Lancet Diabetes Endocrinol. 2014;2(1):38–45.
Rangel-Zuniga OA, Vals-Delgado C, Alcala-Diaz JF, Quintana-Navarro GM, Krylova Y, Leon-Acuna A, Luque RM, Gomez-Delgado F, Delgado-Lista J, Ordovas JM, et al. A set of miRNAs predicts T2DM remission in patients with coronary heart disease: from the CORDIOPREV study. Mol Ther Nucleic Acids. 2021;23:255–63.
Sanchez-Ceinos J, Rangel-Zuniga OA, Clemente-Postigo M, Podadera-Herreros A, Camargo A, Alcala-Diaz JF, Guzman-Ruiz R, Lopez-Miranda J, Malagon MM. miR-223-3p as a potential biomarker and player for adipose tissue dysfunction preceding type 2 diabetes onset. Mol Ther Nucleic Acids. 2021;23:1035–52.
Vals-Delgado C, Alcala-Diaz JF, Molina-Abril H, Roncero-Ramos I, Caspers MPM, Schuren FHJ, Van den Broek TJ, Luque R, Perez-Martinez P, Katsiki N, et al. An altered microbiota pattern precedes type 2 diabetes mellitus development: From the CORDIOPREV study. J Adv Res. 2022;35:99–108.
Taskinen MR. Diabetic dyslipidemia. Atheroscler Suppl. 2002;3(1):47–51.
Han CL, Qu CZ. Cardiovascular risk and safety evaluation of a dual peroxisome proliferator-activated receptor-alpha/gamma agonist, aleglitazar, in patients with type 2 diabetes: a meta-analysis. J Cardiovasc Pharmacol. 2020;75(4):351–7.
Sun Y, Yan D, Cui L, Li G, Sun Y, Hao Z. Sex-disparities in risk factors and atherosclerosis cardiovascular disease in diabetic patients. Postgrad Med. 2021;133(8):860–4.
Ozder A. Lipid profile abnormalities seen in T2DM patients in primary healthcare in Turkey: a cross-sectional study. Lipids Health Dis. 2014;13:183.
Li K, Sinclair AJ, Zhao F, Li D. Uncommon fatty acids and cardiometabolic health. Nutrients. 2018;10(10):1559.
Villasanta-Gonzalez A, Alcala-Diaz JF, Vals-Delgado C, Arenas AP, Cardelo MP, Romero-Cabrera JL, Rodriguez-Cantalejo F, Delgado-Lista J, Malagon MM, Perez-Martinez P, et al. A plasma fatty acid profile associated to type 2 diabetes development: from the CORDIOPREV study. Eur J Nutr. 2022;61(2):843–57.
Wu JH, Micha R, Imamura F, Pan A, Biggs ML, Ajaz O, Djousse L, Hu FB, Mozaffarian D. Omega-3 fatty acids and incident type 2 diabetes: a systematic review and meta-analysis. Br J Nutr. 2012;107(Suppl 2):S214-227.
Wishart DS, Tzur D, Knox C, Eisner R, Guo AC, Young N, Cheng D, Jewell K, Arndt D, Sawhney S, et al. HMDB: the Human Metabolome Database. Nucleic Acids Res. 2007;35(Database issue):D521-526.
Ferderbar S, Pereira EC, Apolinario E, Bertolami MC, Faludi A, Monte O, Calliari LE, Sales JE, Gagliardi AR, Xavier HT, et al. Cholesterol oxides as biomarkers of oxidative stress in type 1 and type 2 diabetes mellitus. Diabetes Metab Res Rev. 2007;23(1):35–42.
Moore LD, Le T, Fan G. DNA methylation and its basic function. Neuropsychopharmacology. 2013;38(1):23–38.
Wang Q, Tang SB, Song XB, Deng TF, Zhang TT, Yin S, Luo SM, Shen W, Zhang CL, Ge ZJ. High-glucose concentrations change DNA methylation levels in human IVM oocytes. Hum Reprod. 2018;33(3):474–81.
Wang Y, Chen L, Pandak WM, Heuman D, Hylemon PB, Ren S. High glucose induces lipid accumulation via 25-hydroxycholesterol DNA-CpG methylation. iScience. 2020;23(5):101102.
Nagao K, Nakamitsu K, Ishida H, Yoshinaga K, Nagai T, Mizobe H, Kojima K, Yanagita T, Beppu F, Gotoh N. Comparison of the lipid-lowering effects of four different n-3 highly unsaturated fatty acids in HepG2 cells. J Oleo Sci. 2014;63(10):979–85.
Elsner M, Gehrmann W, Lenzen S. Peroxisome-generated hydrogen peroxide as important mediator of lipotoxicity in insulin-producing cells. Diabetes. 2011;60(1):200–8.
The CIBEROBN is an initiative of the Instituto de Salud Carlos III, Madrid, Spain. Antonio Camargo is supported by an ISCIII research grant (Programa Miguel-Servet CP14/00114 and CPII19/00007). The funders had no role in the design and conduct of the study; the collection, management, analysis, or interpretation of the data; preparation, review, or approval of the manuscript; or the decision to submit the manuscript for publication. We wish to acknowledge Biobank from the ‘Sistema Sanitario Público de Andalucía’ (Córdoba, Andalusia, Spain) for providing the human biological samples. We also thank EASP (Escuela Andaluza de Salud Publica), in Granada Spain, for performing the randomization process.
The CORDIOPREV study is supported by the Ministerio de Economia y Competitividad, Spain, under the grants AGL2012/39615, PIE14/00005, and PIE14/00031 associated to J.L.-M.; AGL2015-67896-P to J.L.-M. and A.C.; CP14/00114 to A.C.; PI19/00299 to A.C.; DTS19/00007 to A.C.; FIS PI13/00023 to J.D.-L., PI16/01777 to F.P.-J. and P.P.-M.; Antonio Camargo is supported by an ISCIII research contract (Programa Miguel-Servet CPII19/00007); Marina Mora-Ortiz has received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement No 847468; ‘Fundacion Patrimonio Comunal Olivarero’, Junta de Andalucía (Consejería de Salud, Consejeria de Agricultura y Pesca, Consejería de Innovacion, Ciencia y Empresa), ‘Diputaciones de Jaen y Cordoba’, ‘Centro de Excelencia en Investigación sobre Aceite de Oliva y Salud’ and ‘Ministerio de Medio Ambiente, Medio Rural y Marino’, Gobierno de España; ‘Consejeria de Innovación, Ciencia y Empresa, Proyectos de Investigación de Excelencia’, Junta de Andalucía under the grant CVI-7450 obtaiend by J.L.-M.; and we would also like to thank the ‘Fondo Europeo de Desarrollo Regional (FEDER)’.
Ethics approval and consent to participate
CORDIOPREV is registered at ClinicalTrials.gov (number NCT00924937). The trial protocol and subsequent revisions were approved by the Reina Sofia University Hospital Ethics Committee (ref. number 1496/27/03/2009), following the Helsinki Declaration and good clinical practices. All patients signed a written informed consent to participate in the study.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Baseline characteristics of the study population.
Ranking with the 12 most significant metabolites derived from the O-PLS DA model. Figure S2. Adjusted Cox for the analysis of the patient risk scored grouped in ascending tertiles of glycated haemoglobin.
About this article
Cite this article
Mora-Ortiz, M., Alcala-Diaz, J.F., Rangel-Zuñiga, O.A. et al. Metabolomics analysis of type 2 diabetes remission identifies 12 metabolites with predictive capacity: a CORDIOPREV clinical trial study. BMC Med 20, 373 (2022). https://doi.org/10.1186/s12916-022-02566-z
- Insulin resistance
- Prospective human study