Metabolic signature of obesity-associated insulin resistance and type 2 diabetes

1Obesity has become a global health care problem due to associated comorbidities including type 2 diabetes mellitus (T2DM), coronary artery disease (CAD), non-alcoholic fatty liver disease (NAFLD) and cancer. However, a subset of obese individuals exhibit fewer comorbidities than their equally obese counterparts including maintaining their insulin sensitivity as well as having a healthier lipid profile. The underlying protective mechanisms of the metabolically healthy obesity, also known as insulin sensitive obesity, remain unknown.

2Previous studies have suggested that lower levels of inflammatory mediators play a role in the protective phenotype of obese insulin sensitive (OIS) individuals compared to their pathologically obese counterparts, also known as obese insulin resistant (OIR) individuals. Other reports have suggested that OIS individuals show fewer markers of oxidative stress. These two mediators (inflammation and oxidative stress) could potentially be influenced by various genetic and environmental factors. Although evidence of the genetic component remains limited, the environmental effect of certain pollutants and various medications has been previously established.

3Advancement in metabolomic tools including mass spectrometry (MS) technologies has allowed the identification of novel metabolic mediators of disease progression, including obesity associated insulin resistance and T2DM. Recent evidence showed that adipose tissue from OIS, OIR and T2DM individuals exhibit a unique lipidomic signature associated with an increased risk of obesity-associated insulin resistance. Furthermore, metabolomics studies in individuals with T2DM have revealed several diabetes-associated metabolites, including 1,5-anhydroglucitol (1,5-AG), mannose and glucose. Additionally, lipidomics analysis of plasma samples from young adults has revealed that waist circumference was associated with levels of several sphingomyelins, diacylphosphatidylcholines and lysophosphatidylcholines, whereas HOMA-IR was associated with specific diacylphosphatidylcholines, lysophosphatidylcholines and diacylphosphatidylcholines. However, no metabolomics studies have compared the metabolic differences in blood between lean healthy controls, OIS, OIR and T2DM. Such an approach can provide a deeper understanding of the underlying protective mechanisms in those lower risk individuals, and help in the design of novel diagnostic and therapeutic strategies targeting those at higher risk of disease.

4The aim of this study was to employ untargeted metabolomics analysis of blood samples from lean, OIS, OIR and obese-T2DM individuals in order to investigate the metabolic pathways underlying obesity-associated insulin resistance and T2DM.

5Interleukin 6 (IL-6) and leptin ELISAs were from R&D systems (Abingdon, UK). Insulin ELISA was from Mercodia Diagnostics (Uppsala, Sweden). Other chemicals and reagents were from Sigma (Munich, Germany).

6One hundred and seven individuals (75 obese and 32 lean) were recruited at Al Emadi hospital and Hamad Medical Corporation. Lean participants were healthy females visiting the clinic for acne concerns. Obese participants were amongst patients undergoing weight reduction surgery. Subject inclusion criteria included males and females aged over 18 years and under 65 years of age. Subject exclusion criteria included malignancy or other terminal illness, poorly compliant patients, from whatever cause, inability to give informed consent, or involvement in other research projects. All individuals gave their written informed consent. Protocols were approved by Institutional Review Boards of the Anti-Doping Laboratory Qatar (X2017000224) and Weill Cornell Medicine-Qatar (15-00007). Measurements of body mass index (BMI), systolic blood pressure (SBP), diastolic blood pressure (DBP) and mean arterial blood pressure (MAP) were recorded. Fasting blood samples were obtained from all participants. Plasma cholesterol (total, HDL, LDL and triacylglycerol), fasting blood glucose (FBG) and liver function enzymes (total protein, ALP, AST, ALT and bilirubin) were measured by COBAS INTEGRA (Roche Diagnostics, Basil). IL-6, leptin and insulin were determined using commercially available ELISA. Insulin resistance was computed by homeostatic model assessment (HOMA-IR, https://www.dtu.ox.ac.uk/homacalculator/) using 30th percentile (HOMA-IR = 2.4) as a threshold point. Accordingly, obese subjects (BMI > 30) were dichotomized into IS (HOMA-IR < 2.4, n = 20, 6 males and 14 females), IR (HOMA-IR > 2.4, n = 41, 15 males and 26 females) and 14 clinically diagnosed T2DM patients (9 males and 5 females) according to the definition of the American Diabetes Association (ADA) “Standards of Medical Care in Diabetes”.

7Metabolomics profiling was performed using established protocols at Metabolon, Durham, NC, USA. All methods employed a Waters ACQUITY ultra-performance liquid chromatography (UPLC) and a Thermo Scientific Q-Exactive high resolution/accurate mass spectrometer interfaced with a heated electrospray ionization (HESI-II) source and Orbitrap mass analyzer operated at 35,000 mass resolution. The detailed description of the liquid chromatography-mass spectrometry (LC–MS) methodology was previously described. Briefly, serum samples from the 107 participants were methanol extracted to remove the protein fraction. The resulting extract was divided into five fractions: two for analysis by two separate reverse phase (RP)/UPLC-MS/MS methods with positive ion mode electrospray ionization (ESI), one for analysis by RP/UPLC-MS/MS with negative ion mode ESI, one for analysis by hydrophilic interaction chromatography (HILIC)/UPLC-MS/MS with negative ion mode ESI, and one sample was reserved for backup. Raw data was extracted, peak-identified, and quality control-processed using Metabolon’s hardware and software. Compounds were identified by comparison to library entries of purified standards or recurrent unknown entities with more than 3300 commercially available purified standard compounds. Library matches for each compound were checked for each sample and corrected if necessary.

8Statistical analyses were carried out using IBM SPSS version 25, R version 3.2.1 and SIMCA 13.0.1 software (Umetrics, Sweden). Variables with skewed distributions were log transformed or taken the square root of as appropriate to ensure normality. Comparisons were performed with t-test, Wilcoxon–Mann–Whitney and 1-way ANOVA as appropriate. Significance was defined as p ≤ 0.05. Non-parametric tests were used for comparing ordinal or non-normal variables. Metabolomics data were log-transformed to ensure normality. Batch correction was performed by Metabolon by rescaling each metabolite’s median to 1. Principle component analysis (PCA) was performed using version 2.14, http://www.r-project.org/. PCA revealed two main components (PC1 and PC2) that together captured 27% of the variance in the data. Linear regression was performed to identify significant metabolites differentiating study groups (OIS vs OIR and T2DM) and (lean = 0, OIS = 1, OIR = 2, T2DM = 3, denoting disease progression) using the R statistical package (version 2.14, http://www.r-project.org/) after correcting for age, gender, BMI and principle components (PC1 and PC2). PCs represent common signals by the metabolites that contribute to the overall variance in the data and uncover fingerprints of confounders allowing their incorporation into the model by assigning them quantifiable measures. In the first model, the variable study group is categorical whereas the variable disease progression in the second group is continuous. Pathway enrichment analyses were carried out using Chi square tests to identify pathways with metabolites enriched at the top of the list of metabolites ranked by p-value from the linear model since Bonferroni level of significance was not observed. Orthogonal partial least square discriminant analysis (OPLS-DA) was used to compare lean, OIS, OIR and T2DM groups using SIMCA 14 with percentage of missing metabolite values across the samples of 50%. A partial correlation analysis was used to determine metabolic traits of disease (age, BMI, blood pressure, lipids, glucose/insulin/HOMA-IR and liver function enzymes) that exhibit best association with metabolites showing significantly differing levels between disease groups using IBM SPSS version 25, R version 3.2.1.

9Thirty-two lean (BMI = 22.7 ± 2.5 kg/m2, all females) and seventy-five obese and morbidly obese (BMI = 45 ± 6.7 kg/m2, 45 females and 30 males) individuals were recruited at Hamad Medical Corporation and Al Emadi hospital, respectively. Lean individuals were younger and had significantly lower levels of SBP, MAP, triglycerides, triglycerides/HDL ratio, FBG, ALP, ALT and AST than obese individuals. Among obese participants, OIR individuals showed higher FBG than expected, suggesting a high prevalence of undiagnosed T2DM within this group. Therefore, subsequent analyses considered OIR and T2DM groups as one group (all IR) as both groups share obesity and insulin resistance. OIS subjects showed significantly lower MAP and levels of triacylglycerols, FBG, insulin and HOMA-IR than their equally obese all IR (OIR + T2DM) counterparts (Table 1).

10Non-targeted metabolomics of serum samples from the 107 participants was applied to identify metabolites that differentiate OIS vs OIR and OIS vs OIR + T2DM individuals to reveal a metabolic signature of obesity-associated insulin resistance and T2DM. Initial analysis revealed no significant differences in levels of metabolites between OIS and OIR due to their small group sizes (data not shown); however, when combining OIR + T2DM, the linear model revealed 27 metabolites exhibiting significant differences between OIS and OIR + T2DM groups (Table 2). These included metabolites associated with glycolysis, gluconeogenesis and pyruvate metabolism (glucose and 1,5 AG), histidine metabolism (1-methylhistamine, 1-ribosyl-imidazoleacetate and formiminoglutamate) and phospholipid metabolism (choline, glycerophosphoethanolamine and glycerophosphorylcholine). Since the Bonferroni level of significance was not achieved for any of the identified associations, pathway enrichment analysis was performed based on identifying pathways reported by nominally significant metabolites more frequently than can be attributed to random chance. Among the significantly altered metabolic pathways, the phospholipids metabolic pathway was significantly over-represented based on enrichment analysis of the nominally significant metabolites from the group comparisons (p = 3.9E−7). The corresponding metabolites associated with the phospholipids metabolic pathway differentiating OIS from OIR + T2DM included choline, glycerophosphoethanolamine and glycerophosphorylcholine (GPC) (highlighted in Table 2). Figure 1 illustrates levels of significant metabolites that belong to enriched pathways in different study groups. Figure 1 demonstrates higher levels of choline, glycophosphoethanolamine and GPC in OIS compared OIR + T2DM and lean groups. Levels of these metabolites in individual groups are also shown in Additional file 1: Fig. S1.

11An additional a linear model was used to assess the significance of metabolites associated with increased risk of obesity-associated insulin resistance and T2DM (as defined in the method section). Sixty-six metabolites exhibited significant differences with disease progression. The list of metabolites and their associated pathways are shown in Additional file 2: Table S1. These included metabolites associated with glycolysis (glucose), mannose metabolism (mannose), monohydroxy fatty acid (3-hydroxylaurate, 3-hydroxyoctanoate, 3-hydroxydecanoate and 3-hydroxymyristate), medium chain fatty acids (laurate) and urea cycle; arginine and proline metabolism (ADMA + SDMA) among others. Enriched metabolic pathways included glycolysis, gluconeogenesis and the pyruvate metabolic pathway (p = 0.02), fatty acid monohydroxy metabolic pathway, urea cycle metabolic pathway (p = 0.04) and arginine and proline metabolic pathway (p = 0.05). Subsequently, metabolites that showed significant differences with disease progression (Additional file 2: Table S1) within these enriched pathways were identified (Table 3). Figure 2 demonstrates patterns of increased (3-hydroxylaurate, 3-hydrocyoctanoate, 3-hydroxydecanoate, 3-hydroxymyristate, and glucose) or decreased (1,5-AG, ADMA + SDMA, homoarginine, ornithine, 2-oxoarginine) metabolites with disease progression.

12An orthogonal partial least square discriminate analysis (OPLS-DA) comparing subjects from lean, OIS, OIR and T2DM was used for ease of visualization. The model revealed three class-discriminatory components accounting for 48% of the variation in the data due to participant groups (Fig. 3). The score plot in Fig. 3a indicates an x-axis separating the lean group from OIS, OIR and T2DM; the latter group being rather separated along the y-axis. The corresponding loading plot, shown in Fig. 3b, indicates enriched pathways’ associated metabolites significantly differentiating OIS and OIR + T2DM and those associated with disease progression as per linear models. Specifically, higher glucose, choline, GPC, 3-hydroxymyristate and 3-hydroxylaurate and lower 1,5-AG, dimethylarginine (ADMA + SDMA), homoarginine, ornithine and 2-oxoarginine are indicated.

13A partial correlation analysis was used to determine traits of disease best associated with metabolites showing significantly differing levels between disease groups. In essence, the correlation between each of such metabolites and each trait was evaluated after correcting for the effect of all other remaining traits. The correlations that remained significant after such correction are listed in Table 4. The trait of liver function enzymes (ALP and ALT), BMI, TAGs, leptin, insulin and HOMA-IR showed the most significant correlations with levels of metabolites differentiating OIS and OIR + T2DM and those associated with disease progression.

14Obesity triggers a cascade of biochemical changes that increase the risk of various comorbidities including insulin resistance and T2DM. However, some obese individuals seem to be protected against obesity-associated comorbidities. Understanding the underlying mechanisms of this apparently protective phenotype could provide a therapeutic strategy to mitigate the comorbidities associated with pathological obesity. Various studies have investigated the potential mechanisms underlying differences among lean, obese-IS, obese-IR and obese-T2DM individuals, however no study has compared the differences in the metabolic signature among these groups as a means to identify potential diagnostic and therapeutic targets. In this study, untargeted metabolomics analysis of serum samples from lean, OIS, OIR and obese-T2DM individuals was utilized to investigate the metabolic pathways underlying progression of insulin resistance and T2DM. Our novel data indicate that the phospholipid metabolites (choline, glycerophosphoethanolamine and glycerophosphorylcholine) were significantly altered when comparing OIS and OIR + T2DM. Additionally, our data confirmed metabolic changes in several metabolic pathways with obesity-associated insulin resistance and T2DM, including fatty acid and arginine metabolism as well as metabolic markers of liver disease, vascular disease, and diabetes. Therefore, the novel metabolites reported here differentiate the metabolically healthy obese group (OIS) from the pathological obese group (OIR + T2DM) and confirm known biomarkers of obesity-associated insulin resistance with potential diagnostic and therapeutic applications. The causative nature of the identified correlations between metabolites and insulin resistance cannot be ruled out particularly as it is recognized that free fatty acids, for instance, increase insulin resistance. Therefore, future in vitro and in vivo functional studies are warranted where the effects of these metabolites on inducing insulin resistance could confirm their functional relevance.

15Since there were no difference between OIS and OIR likely due to their small group sizes, the analysis was repeated by comparing OIS and combined OIR and T2DM groups as the latter two groups were matched for obesity and insulin resistance. Three phospholipids were found to differentiate between OIS and OIR + T2DM. These included increased levels of choline and GPC in OIS compared to OIR + T2DM and lean groups, suggesting a protective role in obesity-associated insulin resistance. GPC is a natural precursor of phospholipids and a metabolite derived from phosphatidylcholine. It contributes the most to circulating choline levels; therefore, GPC serves as a precursor for acetylcholine. The latter is an important neurotransmitter and a vasodilator that shows a different microvascular reactivity between IR and IS nondiabetic women. Previous studies have reported that dietary choline levels can also lower the risk of fatty liver disease and liver damage. Glycerophosphoethanolamine was another metabolite that differentiated OIS from OIR + T2DM. Glycerophosphoethanolamine represents a membrane degradation product that has been linked to chronic liver disease. The novel associations between higher levels of these phospholipid metabolites and obesity-associated insulin sensitivity could therefore reflect decreased risk of microvascular disease, small vessel disease, lipotoxic cardiac diseases and non-alcoholic liver disease in the OIS group compared to OIR + T2DM group of participants.

16When comparing the metabolic profiles of lean, obese-IS, IR and T2DM individuals, several metabolites significantly changed with disease progression. These included metabolites that were previously reported in association with insulin resistance and T2DM such as glucose and 1,5-AG. Other identified metabolites were reported in association with comorbidities of insulin resistance and T2DM including fatty acid metabolic disorders (such as 3-hydroxylaurate), impairment of liver function and diabetic status (such as 3-hydroxymyristate and homoarginine) and vascular disease (such as dimethylarginine). Other novel arginine metabolites were also found to be significantly changed with disease progression including ornithine (a precursor of arginine, also a medication for hepatic encephalopathy) and 2-oxoarginine (a metabolite of arginine catabolism and a marker of argininemia). Novel metabolites in association with disease progression were also identified including medium chain fatty acids 3-hydroxyoctanoate and 3-hydroxydecanoate that have been reported to be involved in beta-oxidation of longer-chain fatty acids. Previous reports have associated increased plasma levels of 3-hydroxyoctanoate in patients with an inherited deficiency of long-chain 3-hydroxyacyl-CoA dehydrogenase, as a marker of various clinical cases such as recurrent myoglobinuria, hypoketotic hypoglycemic encephalopathy, hypertrophic/dilatative cardiomyopathy, sudden infant death, and fulminant hepatic failure.

17When considering correlations between the identified metabolites and classical mediators of metabolic disease such as age, BMI, lipids, FBG, insulin, HOMA-IR and liver function enzymes, a partial correlation analysis revealed several significant associations. Choline, previously shown to be lower in hepatic damage, was found to positively correlate with leptin and ALT. Despite its positive correlation with ALT, choline was found to be higher in OIS compared to OIR + T2DM, indicating a relationship between this metabolite and the protective phenotype of OIS individuals that requires further investigation. On the other hand, glycerophosphoethanolamine was found to be associated with BMI, suggesting increased levels of this membrane degradation product with obesity.

18As expected, glucose and 1,5-AG, previously shown to be associated with T2DM, were found to correlate significantly with levels of insulin and circulating triacylglycerol levels. When considering metabolites that were significantly associated with obesity-related comorbidities, a significant correlation between levels of 3-hydroxylaurate and ALP, was revealed. This suggests that 3-hydroxylaurate, a medium chain fatty acid that is associated with intolerance to prolonged fasting and recurrent episodes of hypoglycemic coma, may constitute a novel marker of fatty liver disease. Similarly, 3-hydroxyocanoate was also found to be associated with ALP and BMI, suggesting that it may also be a novel marker of obesity-associated fatty liver disease. 3-hydroxydecanoate was also found to be associated with BMI, suggesting increased levels of another medium-chain fatty acid with a role in the beta-oxidation and obesity. Ornithine, previously shown to be associated with hepatic damage, was found to be associated with leptin and ALT, providing a further evidence of its association with obesity associated non-alcoholic fatty liver disease.

19This has a number of limitations including the relatively low number of participants per group and the cross-sectional nature of the study limited the interpretation of the findings from a pathophysiological point of view. The observational nature of the findings requires functional validation before suggesting any causalities, especially as some findings were based on weak to moderate associations. Furthermore, since blood samples were collected at multiple sites, a batch effect may have occurred, but this was mitigated by standardized protocols for sample collection, processing and storage. It is possible that other unmeasured factors may have impacted our data including dietary habits, medication/supplements and other unknown environmental factors; however, inclusion of principle components in the regression model may have captured part of these potential confounding factors. Finally, controls were not matched for age and gender compared to the study groups, adding an additional variable; however, both age and gender were corrected for in the analysis, but their influence over metabolic differences cannot be ruled out.

20In the comparison between equally obese insulin sensitive and insulin resistance individuals, phospholipid metabolites including choline, glycerophosphoethanolamine and glycerophosphorylcholine (GPC) were significantly altered. In addition, several metabolites were identified and were confirmatory for insulin resistance and T2DM (such as glucose and 1,5-AG) or their comorbidities (such as 3-hydroxylaurate, 3-hydroxymyristate, homoarginine and dimethylarginine). This pilot study also identified novel metabolic markers such as the medium chain fatty acids 3-hydroxyoctanoate and 3-hydroxydecanoate and highlighted their potential link to non-alcoholic fatty liver disease, a hallmark of increased risk of obesity-associated insulin resistance. Further studies are needed to confirm these associations in prospective cohorts and to investigate their functional relevance.