Human Insulin Resistance Is Associated With Increased Plasma Levels of 12α-Hydroxylated Bile Acids
Bile acids (BAs) exert pleiotropic metabolic effects, and physicochemical properties of different BAs affect their function. In rodents, insulin regulates BA composition, in part by regulating the BA 12α-hydroxylase CYP8B1. However, it is unclear whether a similar effect occurs in humans. To address this question, we examined the relationship between clamp-measured insulin sensitivity and plasma BA composition in a cohort of 200 healthy subjects and 35 type 2 diabetic (T2D) patients. In healthy subjects, insulin resistance (IR) was associated with increased 12α-hydroxylated BAs (cholic acid, deoxycholic acid, and their conjugated forms). Furthermore, ratios of 12α-hydroxylated/non–12α-hydroxylated BAs were associated with key features of IR, including higher insulin, proinsulin, glucose, glucagon, and triglyceride (TG) levels and lower HDL cholesterol. In T2D patients, BAs were nearly twofold elevated, and more hydrophobic, compared with healthy subjects, although we did not observe disproportionate increases in 12α-hydroxylated BAs. In multivariate analysis of the whole dataset, controlling for sex, age, BMI, and glucose tolerance status, higher 12α-hydroxy/non–12α-hydroxy BA ratios were associated with lower insulin sensitivity and higher plasma TGs. These findings suggest a role for 12α-hydroxylated BAs in metabolic abnormalities in the natural history of T2D and raise the possibility of developing insulin-sensitizing therapeutics based on manipulations of BA composition.
Bile acids (BAs) regulate multiple aspects of metabolism (1). The oldest known roles for BAs are in cholesterol homeostasis: BAs promote absorption of dietary cholesterol, solubilize biliary cholesterol for excretion, and are themselves the major product of liver cholesterol catabolism. More recently, BAs have been recognized as signaling molecules, repressing their own synthesis as well as that of triglyceride (TG) in the liver and promoting incretin secretion in the gut. Many of these actions occur through activation of BA receptors such as the nuclear receptor FXR and the G protein–coupled receptor TGR5 (1). In addition, BAs can act as mediators of the unfolded protein response (2).
Given the variety of metabolic pathways modulated by BAs, pharmacological manipulation of these molecules has become of interest for obesity and diabetes. Pharmacological use of exogenous BAs is limited by their narrow therapeutic window due to the toxicity of selected secondary BAs and the effect of chenodeoxycholic acid (CDCA) to increase LDL cholesterol (LDL-C) (3). Nonetheless, ursodeoxycholic acid (UDCA) has been gainfully used for dissolution of cholesterol gallstones or treatment of primary biliary cirrhosis (4,5). BA sequestrants are effectively used to reduce LDL-C and improve glycemic control in type 2 diabetes (T2D) (6). FXR and TGR5 agonists have also been developed, although their use has been limited by efficacy (for TGR5 agonists) or side effects (increased LDL-C and pruritus, for FXR agonists) (3,7). There has also been a surge of interest in the effects of bariatric surgery on endogenous BAs, and whether BA signaling can account for some of its insulin-sensitizing effects (8,9).
The primary products of the BA synthetic pathway in humans are CDCA and cholic acid (CA). They in turn are modified by gut microbiota to yield lithocholic acid (LCA) and deoxycholic acid (DCA), two secondary BAs, respectively, as well as other secondary and tertiary BAs (Fig. 1). Conversion of CDCA to CA in the liver requires 12α-hydroxylation by the 12α-hydroxylase, encoded by Cyp8b1. Individual BAs differ in their physicochemical properties, such that individual BAs or subsets of BAs can affect downstream metabolic pathways differentially (1,10). For example, whereas CDCA is a potent agonist of FXR, with a Km in the micromolar range, DCA and LCA have been suggested as partial antagonists of FXR (11).
We have recently shown that in rodents, insulin inhibits the production of 12α-hydroxy BAs by suppressing transcription of Cyp8b1 (12). Conversely, the insulin resistance (IR) characteristic of T2D would be expected to increase levels of 12α-hydroxylated BAs, as has indeed been observed in multiple mouse models of diabetes (13–17). However, in determining the consequences of this shift in BA composition in human diabetes, one should be mindful of important differences in BA metabolism between mammalian species. Unlike in humans, in rodents, CDCA can be 6-hydroxylated, leading to the formation of hydrophilic muricholic acids that are unable to promote dietary cholesterol absorption (18) or agonize FXR (19,20). It is therefore unclear how mechanisms of BA action identified in mice would apply to human diabetes and IR.
The goal of this work was to critically test the hypothesis that there exists a correlation between insulin sensitivity and BA 12α-hydroxylation in humans. To this end, we measured 13 BA species in plasma collected after an overnight fast from 200 healthy individuals with normal glucose tolerance and a range of insulin sensitivities, as assessed by hyperinsulinemic-euglycemic clamps, and from 35 patients with T2D. Secondly, we sought to characterize the relationship between BA subtypes, BA levels, and metabolic parameters.
RESEARCH DESIGN AND METHODS
The study population included 200 nondiabetic subjects (103 women and 97 men) and 35 patients with T2D (18 women and 17 men). The nondiabetic subjects were part of the Relationship between Insulin Sensitivity and Cardiovascular Risk (RISC) study (www.egir.org), which recruited apparently healthy Caucasians in 19 centers in 14 European countries between June 2002 and July 2004 (21). Study participants were between 30 and 60 years of age, and blood pressure (<140/<90 mmHg), plasma cholesterol (<7.8 mmol/L), TGs (<4.6 mmol/L), and fasting and 2-h plasma glucose (<7.0 and 11.1 mmol/L) were within established limits. Exclusion criteria were the presence of overt cardiovascular disease, chronic systemic diseases (hypertension, diabetes, dyslipidemia, inflammatory, and neoplastic diseases), class III obesity, and carotid stenosis >40% or calcified carotid plaques and drug treatment for hypertension, diabetes, dyslipidemia, or obesity or steroid use. The 200 individuals in the current study were selected to cover the entire 8–10-fold range of insulin sensitivity of the original cohort (21). The T2D group was recruited at our Clinical Research Center; their mean diabetes duration was 9 years, their mean HbA1c was 7.9 ± 1.6%, and they were on stable treatment with metformin ± sulfonylureas (nine patients were also receiving insulin). The protocol was approved by the institutional ethics committee of the University of Pisa.
All study subjects received a euglycemic insulin clamp (using an insulin infusion rate of 240 pmol ⋅ min−1 per m2 of body surface area), performed as previously described (22). Insulin sensitivity was calculated as the insulin-mediated whole-body glucose disposal during the final 40 min of the 2-h clamp normalized to the fat-free mass (M) divided by the steady-state plasma insulin concentration during the same time interval (M/I, in units of µmol ⋅ min−1 ⋅ kgffm−1 ⋅ nM−1) (23). In the RISC subjects, insulin secretion was derived from mathematical modeling of the plasma glucose and C-peptide response to an oral glucose tolerance test with the use of C-peptide deconvolution (24). A fatty liver index (FLI) was calculated as previously described (25). Plasma lipid profile was determined by standard methods, plasma hormone concentrations (insulin, proinsulin, glucagon, and adiponectin) were measured by radioimmunoassay, and circulating metabolites were determined by a metabolomics platform (26).
Aliquots (0.1 mL) of plasma collected after an overnight fast were spiked with internal standard (10 μL cholate-d4, 25 μmol/L), mixed with 1 mL ice-cold acetonitrile, vortexed, and centrifuged for 10 min at 11,000g. Supernatant was dried under nitrogen at 45°C and resuspended in 100 μL 55/45 (volume for volume) MeOH/H2O, containing 5 mmol/L ammonium formate. Ten microliters was injected into LC/MS (Waters Quattro Micro with Waters 2795 Alliance HPLC). Quantitative standard curves were used, and deuterated internal standards were used to measure recovery. Subgrouping and hydrophobicity of the measured BAs are listed in Table 1. Hydrophobicity index (HI) was calculated using the published method (27).
Data are given as mean ± SD or median and (interquartile range) for normally or nonnormally distributed variables, respectively. Group differences were analyzed by Kruskal-Wallis or Mann-Whitney test, as appropriate. Principal component analysis was carried out on correlations using a Varimax method for orthogonal factor rotation. Multivariate and partial regression analyses were performed by standard methods. A P value ≤0.05 was considered statistically significant. All statistical analyses were performed using JMP7.0.
The anthropometric and metabolic characteristics of our cohort of nondiabetic subjects are summarized in Table 2. Upon subdividing the group into quartiles of clamp-derived insulin sensitivity (M/I), insulin-resistant subjects had a similar age as insulin-sensitive individuals but were heavier and more hyperinsulinemic. Among BAs, CDCA and its conjugated forms were the most abundant species, accounting for 43 ± 15% of the total BA pool, followed by DCA and its conjugates (29 ± 11% of the pool). Circulating BAs comprised roughly equal parts of primary versus secondary, conjugated versus unconjugated, and 12α-hydroxylated versus non–12α-hydroxylated BAs. The HI averaged 0.47. Across quartiles of M/I, more insulin-resistant individuals tended to have slightly higher total plasma BAs (Table 2).
Based on our findings in experimental animals, we predicted that IR would increase the ratio of 12α-hydroxylated to non–12α-hydroxylated BAs; we therefore analyzed these subcategories. We found that the sum of 12α-hydroxylated species (CA, DCA, and their conjugates) increased significantly with IR, whereas non–12α-hydroxylated species, as a category, were unaltered. The increase of 12α-hydroxylated BAs was not accounted for by any individual species, as five out of six 12α-hydroxylated species showed the same upward trend. Conversely, there were no significant associations of insulin sensitivity with the balance of unconjugated versus conjugated BAs or of secondary versus primary BAs, although the increase of DCA caused trends in these two categories (Table 2). An important factor in determining BA bioactivity is their relative hydrophobicity, which predicts their ability to traverse membranes, bind FXR, and promote lipid and cholesterol absorption in the intestine. Thus, we also examined the hydrophobicity of the BA pool by calculating the estimated HI of the BA pool of each subject (27). We found no significant alteration in hydrophobicity across quartiles of insulin sensitivity (Table 2).
If variations of BA pool composition were a component of the IR syndrome, one would expect to see associations with other metabolic parameters linked to this condition. We used principal components analysis to dissect out relationships between BAs and other metabolic parameters in healthy subjects (Supplementary Table 1). The ratio of 12α-hydroxy/non–12α-hydroxy BAs had its highest loading on a factor characterized, in addition to a lower M/I, by higher BMI, waist-to-hip ratio, and mean blood pressure; increased fasting and postglucose insulin secretion; increased lipolysis; higher glucose, glucagon, and proinsulin; higher TGs; lower HDL cholesterol (HDL-C); and higher γGT, ALAT, and FLI (Fig. 2). Thus, an increase in this ratio was associated with canonical features of IR. In addition, a higher ratio was also associated with changes in blood amino acids (raised glutamate and reduced glycine) and with increased circulating levels of α-hydroxybutyrate, a recently identified amino acid metabolite and marker of metabolic overload and incident diabetes (26).
Total plasma BA levels had their highest (negative) loading on the second independent factor characterized by increased levels of the major fatty acids (oleate, palmitate, and stearate), acylcarnitines, β-hydroxybutyrate, and α-hydroxybutyrate (Fig. 2). Thus, low total BA levels are associated with changes in amino acid and fatty acid utilization and ketogenesis but not closely with parameters of glucose metabolism or insulin sensitivity. The third factor in this analysis, characterized by high levels of acyl-glycerophosphocholines and transaminases, was associated with a high loading of the HI.
Finally, in order to get a quantitative estimate of the impact of a high 12α-hydroxy/non–12α-hydroxy ratio, we compared the upper and lower quartiles of the ratio. Individuals in the upper quartile (a fourfold increased ratio, associated with significant increments in total BA levels and in the ratio of unconjugated/conjugated and secondary/primary species) (Fig. 3) showed no difference in BMI compared with the lowest quartile, but had higher insulin, free fatty acids, TGs, and FLI and lower M value, HDL-C, and adiponectin (Fig. 4). Overall, these data show that increased BA 12α-hydroxylation is associated with the full complement of metabolic features of IR, independently of BMI.
T2D patients in our cohort were older and heavier than healthy controls and showed fasting hyperglycemia and hyperinsulinemia, and severe insulin resistance. These patients were on treatment with antidiabetic medications, a potential confounding factor. In this group, there was a twofold increase in total plasma BAs (Table 2 and Fig. 5), as well as increases of 9 out of the 13 individual BAs measured. However, we found no disproportionate alterations in BA subtypes, such as 12α-hydroxylated versus non–12α-hydroxylated, unconjugated versus conjugated, or secondary versus primary BAs. To correct for potential confounders introduced by differences in BMI, we reanalyzed the data by comparing T2D patients to a subgroup of BMI-matched healthy subjects. However, we still found no disproportionate increase in 12α-hydroxylated BAs (data not shown). In contrast, the HI was significantly higher in the T2D groups than in the healthy control subjects regardless of covariate adjustment (Table 2).
Finally, we examined the relationship between insulin sensitivity and BA 12α-hydroxylation in the whole dataset, including healthy subjects and T2D patients. In multivariate analysis controlling for sex, age, BMI, and glucose tolerance status, we found that higher ratios of 12α-hydroxy/non–12α-hydroxy BAs, but not total BAs, primary versus secondary, or conjugated versus unconjugated BAs, were associated with lower insulin sensitivity (Fig. 6A) and higher plasma TG levels (Fig. 6B). As plasma TGs were twice as high in T2D patients as in controls (1.59 [0.81] vs. 0.85 [0.56] mmol/L, P < 0.0001), we tested whether the HI might be related to the hypertriglyceridemia. In a multivariate regression model adjusting for sex, age, and BMI but also for fasting plasma glucose concentrations and the 12α-hydroxyl/non–12α-hydroxyl BA, a higher HI was positively related to higher plasma TG levels (Fig. 6C).
In this work, we show that in nondiabetic subjects, plasma 12α-hydroxylated BAs are disproportionately increased in association with IR. This finding provides evidence that BA composition is altered within the normal range of insulin sensitivity, and is consistent with insulin-dependent regulation of the 12α-hydroxylase CYP8B1. Our prior work has shown that Cyp8b1 is a target gene of transcription factor FoxO1, whose ability to promote Cyp8b1 transcription is inhibited by insulin via Akt-dependent phosphorylation and nuclear exclusion of FoxO1 (12). These findings dovetail with earlier studies of BA composition, as increased 12α-hydroxylated BAs have been observed in mice lacking hepatic insulin receptors and rodent models of diabetes (13–17). Because FoxO is expected to be more active in IR, owing to decreased Akt-mediated nuclear exclusion, the finding that FoxO1 activates Cyp8b1 reveals a plausible mechanism for upregulation of 12α-hydroxylated BA synthesis in these models.
Surprisingly, a preferential increase in 12α-hydroxylated BAs did not emerge from sampling T2D patients, as we did not find disproportionate increases in these molecules relative to non–12α-hydroxylated BAs despite the marked IR (Table 2). These findings suggest that plasma BA levels and composition are affected by other aspects of T2D in addition to IR. We also note limitations of our study that may confound these results. First, the 35 diabetic patients in our study were taking antidiabetic medications, which may affect BA synthesis and turnover through liver insulin sensitization or other indirect effects. Second, our measurements used fasting rather than postprandial plasma, which may mask some of the differences. A previous study of 16 T2D patients showed increased synthesis of CA and input rate of DCA as compared with BMI-matched controls, causing an increase in postprandial plasma DCA levels and an apparent increase in the 12α-hydroxy/non–12α-hydroxy ratio, but no difference in total plasma BAs (28). In a recent metabolomics screen of 40 subjects with self-reported diabetes (29), DCA, but not CA, was detectable in a higher number of diabetic subjects as compared with a group of 60 control nondiabetic subjects (consistent with the higher levels we observed) (Table 2), but total BAs were not measured; thus a relative preponderance of 12α-hydroxylated species was not reported (29). In addition to medications and postprandial versus fasting plasma, other differences in our study population may account for differences in our findings, such as time since diabetes diagnosis, metabolic control, recruiting period, and geographic location. Moreover, an increased synthesis of 12α-hydroxylated BAs may not be reflected in steady-state plasma levels, because plasma levels are determined by multiple factors, including rate of reabsorption in the small intestine and rate of reuptake in the liver sinusoids. Tracer studies will be required to address these possibilities. On the other hand, evidence from genetically obese T2D ob/ob mice is consistent with the increase in total plasma BAs that we observed; in these mice, the hyperglycemia epigenetically augments expression of Cyp7a1, the rate-limiting enzyme of BA synthesis, whereby BA synthesis and serum BA concentrations are increased (17). FoxO1 has also been implicated in promoting expression of Cyp7a1, although the direction of this effect is debated and may be different between mice and humans (30,31).
Mechanistically, it bears emphasizing that FoxO1 is regulated by both insulin and glucose, and that the lack of a preferential increase in 12α-hydroxylated BAs seen in our study could be explained by hyperglycemia-induced FoxO degradation, and thus inability to fully activate CYP8B1 (32,33).
Of interest is the finding that the HI was significantly higher in T2D patients as compared with nondiabetic subjects, due to the higher circulating levels of DCA and its conjugated forms, relative to CA, CDCA, and their conjugates (Table 2). More hydrophobic BA species are more toxic, especially to the liver (34,35). Hydrophobic BAs are also better able to cross, and are more disruptive to, cell membranes (35,36); are more likely to activate additional BA-sensitive signaling pathways such as protein kinase C, extracellular signal–related kinase 1/2, PXR, and VDR (7,34,37); and are more efficient in absorption of dietary lipids (18). The observed association of the HI with raised transaminases and circulating phospholipid concentrations (Fig. 2) plausibly tracks with this notion. The independent relationship of the HI to TGs is also of interest, although we are unaware of mechanistic data linking BA hydrophobicity to fasting TG levels.
In healthy subjects, we show that an increased ratio of 12α-hydroxy/non–12α-hydroxy BAs is associated with classic features of IR. These include metabolic functions known to be influenced by BAs, including insulin secretion, TG levels, and HDL-C turnover (1,38,39). In comparison, changes in total BA levels were more closely related to changes in fatty acid utilization and ketogenesis. The implications of these findings are that changes in BA types may be an important feature in the pathophysiology of IR, and that the effects of these changes may be dissociable from the effects of altering total BA levels. In addition, the identification of altered 12α-hydroxy BA levels in euglycemic/insulin-resistant individuals may be exploited as a biomarker of hepatic insulin sensitivity and CVD risk.
Mechanistic studies will be required to identify whether and how BA composition affects metabolic abnormalities that arise in IR. However, we speculate that changes in endogenous ligands for BA receptors affect a wide range of BA-regulated pathways. It has been shown that in HepG2 cells, a variety of different BA species can effectively repress Cyp7a1, the rate-limiting step of BA synthesis (40); thus, this pathway is likely to be sensitive to total BA levels. On the other hand, BA species differ considerably in their ability to activate transcription of Bsep, the bile salt export pump, and CDCA is the best activator, through its potent activation of FXR (40). Interestingly, DCA, the most abundant plasma 12α-hydroxylated BA, which we found to be increased in IR, has been shown to interfere with the ability of CDCA to cause recruitment of coactivator SRC-1 to FXR (20). Is it possible that the increases in DCA during IR cause a decline in activation of some FXR pathways? Interventional studies will be required to address this important question.
In addition to alterations in FXR ligands during IR, there may be effects directly on FXR mRNA (Nr1h4), which has been shown to be regulated by glucose and insulin. In three rodent models of hyperglycemia, streptozotocin-induced diabetic rats, diabetic Zucker rats, and liver insulin receptor knockout mice, liver FXR expression is reduced (14,41). Insulin treatment normalizes this defect in streptozotocin-treated rats (41). It is possible that a hyperglycemia-induced reduction in liver FXR mRNA may contribute to metabolic outcomes during diabetes, such as a release of FXR-mediated repression of Cyp7a1 and subsequent increase in BA synthesis (41). Finally, evidence now suggests a role for β-cell FXR in insulin secretion; perhaps changes in FXR ligands or mRNA expression during IR impair glucose-stimulated insulin secretion (38,39,42).
Overall, these findings support a model whereby insulin regulates the production of endogenous 12α-hydroxy BAs that subsequently affect metabolism in multiple ways. The relationship between these metabolites and other phenotypes of IR warrants further investigation. It is conceivable that BA composition could be a target for new therapeutics in IR.
This work was supported in part by Columbia Clinical and Translational Science Award Grant UL1-TR-000040 (National Center for Advancing Translational Sciences/National Institutes of Health [NIH]), DK-63608 (NIH), the Schaefer Scholar Award, a Core Usage Funding Award from Columbia University College of Physicians and Surgeons, and K99-HL-111206 (NIH).
E.F. has served as a speaker and consultant for Boehringer Ingelheim, Merck, Sanofi, Eli Lilly and Company, Johnson & Johnson, Astellas, Daiichi Sankyo, Bristol-Myers-Squibb/AstraZeneca, and Novartis. No other potential conflicts of interest relevant to this article were reported.
R.A.H. conceived of this study, designed and implemented experiments, analyzed data, and wrote the manuscript. B.A. and S.C. performed the in vivo studies, contributed to discussions, and edited the manuscript. D.A. contributed to discussions and edited the manuscript. E.F. is the principal investigator of the RISC study, designed experiments, analyzed data, and wrote the manuscript. R.A.H. and E.F. are the guarantors of this work and, as such, had full access to all the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis.
This work was presented at the 73rd Scientific Sessions of the American Diabetes Association, Chicago, Illinois, 21–25 June 2013.
The authors wish to express their gratitude to the investigators of the RISC study, who originally carried out the phenotyping of the nondiabetic subjects.
This article contains Supplementary Data online at http://diabetes.diabetesjournals.org/lookup/suppl/doi:10.2337/db13-0639/-/DC1.
See accompanying commentary, p. 3987.
- Received April 22, 2013.
- Accepted July 15, 2013.
- © 2013 by the American Diabetes Association.
Readers may use this article as long as the work is properly cited, the use is educational and not for profit, and the work is not altered. See http://creativecommons.org/licenses/by-nc-nd/3.0/ for details.