Advertisement
Clinical Research and Public HealthMetabolism Free access | 10.1172/JCI141828
1Center for Human Nutrition and Atkins Center of Excellence in Obesity Medicine, and
2Mallinckrodt Institute of Radiology, Washington University School of Medicine, St. Louis, Missouri, USA.
3Departments of Medicine and Pharmacology, University of California, San Diego, La Jolla, California, USA.
Address correspondence to: Samuel Klein, Center for Human Nutrition, Washington University School of Medicine, 660 South Euclid Avenue, Campus Box 8031, St. Louis, Missouri 63110, USA. Phone: 314.362.8708; Email: sklein@wustl.edu.
Authorship note: VC and SCB contributed equally to this work.
Find articles by Cifarelli, V. in: JCI | PubMed | Google Scholar
1Center for Human Nutrition and Atkins Center of Excellence in Obesity Medicine, and
2Mallinckrodt Institute of Radiology, Washington University School of Medicine, St. Louis, Missouri, USA.
3Departments of Medicine and Pharmacology, University of California, San Diego, La Jolla, California, USA.
Address correspondence to: Samuel Klein, Center for Human Nutrition, Washington University School of Medicine, 660 South Euclid Avenue, Campus Box 8031, St. Louis, Missouri 63110, USA. Phone: 314.362.8708; Email: sklein@wustl.edu.
Authorship note: VC and SCB contributed equally to this work.
Find articles by Beeman, S. in: JCI | PubMed | Google Scholar
1Center for Human Nutrition and Atkins Center of Excellence in Obesity Medicine, and
2Mallinckrodt Institute of Radiology, Washington University School of Medicine, St. Louis, Missouri, USA.
3Departments of Medicine and Pharmacology, University of California, San Diego, La Jolla, California, USA.
Address correspondence to: Samuel Klein, Center for Human Nutrition, Washington University School of Medicine, 660 South Euclid Avenue, Campus Box 8031, St. Louis, Missouri 63110, USA. Phone: 314.362.8708; Email: sklein@wustl.edu.
Authorship note: VC and SCB contributed equally to this work.
Find articles by Smith, G. in: JCI | PubMed | Google Scholar
1Center for Human Nutrition and Atkins Center of Excellence in Obesity Medicine, and
2Mallinckrodt Institute of Radiology, Washington University School of Medicine, St. Louis, Missouri, USA.
3Departments of Medicine and Pharmacology, University of California, San Diego, La Jolla, California, USA.
Address correspondence to: Samuel Klein, Center for Human Nutrition, Washington University School of Medicine, 660 South Euclid Avenue, Campus Box 8031, St. Louis, Missouri 63110, USA. Phone: 314.362.8708; Email: sklein@wustl.edu.
Authorship note: VC and SCB contributed equally to this work.
Find articles by Yoshino, J. in: JCI | PubMed | Google Scholar
1Center for Human Nutrition and Atkins Center of Excellence in Obesity Medicine, and
2Mallinckrodt Institute of Radiology, Washington University School of Medicine, St. Louis, Missouri, USA.
3Departments of Medicine and Pharmacology, University of California, San Diego, La Jolla, California, USA.
Address correspondence to: Samuel Klein, Center for Human Nutrition, Washington University School of Medicine, 660 South Euclid Avenue, Campus Box 8031, St. Louis, Missouri 63110, USA. Phone: 314.362.8708; Email: sklein@wustl.edu.
Authorship note: VC and SCB contributed equally to this work.
Find articles by Morozov, D. in: JCI | PubMed | Google Scholar
1Center for Human Nutrition and Atkins Center of Excellence in Obesity Medicine, and
2Mallinckrodt Institute of Radiology, Washington University School of Medicine, St. Louis, Missouri, USA.
3Departments of Medicine and Pharmacology, University of California, San Diego, La Jolla, California, USA.
Address correspondence to: Samuel Klein, Center for Human Nutrition, Washington University School of Medicine, 660 South Euclid Avenue, Campus Box 8031, St. Louis, Missouri 63110, USA. Phone: 314.362.8708; Email: sklein@wustl.edu.
Authorship note: VC and SCB contributed equally to this work.
Find articles by
Beals, J.
in:
JCI
|
PubMed
|
Google Scholar
|
1Center for Human Nutrition and Atkins Center of Excellence in Obesity Medicine, and
2Mallinckrodt Institute of Radiology, Washington University School of Medicine, St. Louis, Missouri, USA.
3Departments of Medicine and Pharmacology, University of California, San Diego, La Jolla, California, USA.
Address correspondence to: Samuel Klein, Center for Human Nutrition, Washington University School of Medicine, 660 South Euclid Avenue, Campus Box 8031, St. Louis, Missouri 63110, USA. Phone: 314.362.8708; Email: sklein@wustl.edu.
Authorship note: VC and SCB contributed equally to this work.
Find articles by Kayser, B. in: JCI | PubMed | Google Scholar
1Center for Human Nutrition and Atkins Center of Excellence in Obesity Medicine, and
2Mallinckrodt Institute of Radiology, Washington University School of Medicine, St. Louis, Missouri, USA.
3Departments of Medicine and Pharmacology, University of California, San Diego, La Jolla, California, USA.
Address correspondence to: Samuel Klein, Center for Human Nutrition, Washington University School of Medicine, 660 South Euclid Avenue, Campus Box 8031, St. Louis, Missouri 63110, USA. Phone: 314.362.8708; Email: sklein@wustl.edu.
Authorship note: VC and SCB contributed equally to this work.
Find articles by
Watrous, J.
in:
JCI
|
PubMed
|
Google Scholar
|
1Center for Human Nutrition and Atkins Center of Excellence in Obesity Medicine, and
2Mallinckrodt Institute of Radiology, Washington University School of Medicine, St. Louis, Missouri, USA.
3Departments of Medicine and Pharmacology, University of California, San Diego, La Jolla, California, USA.
Address correspondence to: Samuel Klein, Center for Human Nutrition, Washington University School of Medicine, 660 South Euclid Avenue, Campus Box 8031, St. Louis, Missouri 63110, USA. Phone: 314.362.8708; Email: sklein@wustl.edu.
Authorship note: VC and SCB contributed equally to this work.
Find articles by Jain, M. in: JCI | PubMed | Google Scholar
1Center for Human Nutrition and Atkins Center of Excellence in Obesity Medicine, and
2Mallinckrodt Institute of Radiology, Washington University School of Medicine, St. Louis, Missouri, USA.
3Departments of Medicine and Pharmacology, University of California, San Diego, La Jolla, California, USA.
Address correspondence to: Samuel Klein, Center for Human Nutrition, Washington University School of Medicine, 660 South Euclid Avenue, Campus Box 8031, St. Louis, Missouri 63110, USA. Phone: 314.362.8708; Email: sklein@wustl.edu.
Authorship note: VC and SCB contributed equally to this work.
Find articles by
Patterson, B.
in:
JCI
|
PubMed
|
Google Scholar
|
1Center for Human Nutrition and Atkins Center of Excellence in Obesity Medicine, and
2Mallinckrodt Institute of Radiology, Washington University School of Medicine, St. Louis, Missouri, USA.
3Departments of Medicine and Pharmacology, University of California, San Diego, La Jolla, California, USA.
Address correspondence to: Samuel Klein, Center for Human Nutrition, Washington University School of Medicine, 660 South Euclid Avenue, Campus Box 8031, St. Louis, Missouri 63110, USA. Phone: 314.362.8708; Email: sklein@wustl.edu.
Authorship note: VC and SCB contributed equally to this work.
Find articles by
Klein, S.
in:
JCI
|
PubMed
|
Google Scholar
|
Authorship note: VC and SCB contributed equally to this work.
Published November 9, 2020 - More info
BACKGROUND. Data from studies conducted in rodent models have shown that decreased adipose tissue (AT) oxygenation is involved in the pathogenesis of obesity-induced insulin resistance. Here, we evaluated the potential influence of AT oxygenation on AT biology and insulin sensitivity in people.
METHODS. We evaluated subcutaneous AT oxygen partial pressure (pO2); liver and whole-body insulin sensitivity; AT expression of genes and pathways involved in inflammation, fibrosis, and branched-chain amino acid (BCAA) catabolism; systemic markers of inflammation; and plasma BCAA concentrations, in 3 groups of participants that were rigorously stratified by adiposity and insulin sensitivity: metabolically healthy lean (MHL; n = 11), metabolically healthy obese (MHO; n = 15), and metabolically unhealthy obese (MUO; n = 20).
RESULTS. AT pO2 progressively declined from the MHL to the MHO to the MUO group, and was positively associated with hepatic and whole-body insulin sensitivity. AT pO2 was positively associated with the expression of genes involved in BCAA catabolism, in conjunction with an inverse relationship between AT pO2 and plasma BCAA concentrations. AT pO2 was negatively associated with AT gene expression of markers of inflammation and fibrosis. Plasma PAI-1 increased from the MHL to the MHO to the MUO group and was negatively correlated with AT pO2, whereas the plasma concentrations of other cytokines and chemokines were not different among the MHL and MUO groups.
CONCLUSION. These results support the notion that reduced AT oxygenation in individuals with obesity contributes to insulin resistance by increasing plasma PAI-1 concentrations and decreasing AT BCAA catabolism and thereby increasing plasma BCAA concentrations.
TRIAL REGISTRATION. ClinicalTrials.gov NCT02706262.
FUNDING. This study was supported by NIH grants K01DK109119, T32HL130357, K01DK116917, R01ES027595, P42ES010337, DK56341 (Nutrition Obesity Research Center), DK20579 (Diabetes Research Center), DK052574 (Digestive Disease Research Center), and UL1TR002345 (Clinical and Translational Science Award); NIH Shared Instrumentation Grants S10RR0227552, S10OD020025, and S10OD026929; and the Foundation for Barnes-Jewish Hospital.
Insulin-resistant glucose metabolism is the most common metabolic abnormality associated with obesity, and is an important risk factor for type 2 diabetes and cardiovascular disease (1). Although insulin sensitivity generally declines with increasing adiposity, there is considerable heterogeneity in the severity of insulin resistance at any given value for body mass index (BMI) (2); in fact, some people with obesity are even considered to be metabolically healthy with normal insulin sensitivity (3). It is possible that differences in adipose tissue (AT) oxygenation contribute to the variability in insulin resistance. Studies conducted in rodent models have shown that the expansion of AT mass without adequate neovascularization can cause AT hypoxia, leading to AT inflammation and fibrosis, and systemic insulin resistance (4), whereas increasing AT oxygen content decreases AT inflammation and fibrosis and increases whole-body insulin sensitivity in obese mice (5). The link between alterations in AT biology and systemic insulin action suggests a communication between AT and other organs. The systemic signaling mechanism could be related to increased plasma concentrations of proinflammatory cytokines and chemokines secreted by AT immune cells (6), increased plasma branched-chain amino acids (BCAAs) caused by decreased AT BCAA oxidative metabolism (7–11), or both.
The effect of excess adiposity on AT oxygenation in people is unclear because of conflicting data from different studies that have reported that AT partial pressure of oxygen (pO2) was lower (5, 12–15) or not different (16, 17) in people with obesity than in people who are lean, and that AT pO2 was higher (5), lower (17), or not different (13, 14) in people with metabolically healthy obesity (MHO) than in those with metabolically unhealthy obesity (MUO). In addition, the relationship between AT pO2 and whole-body insulin sensitivity in people is unclear because of conflicting data from different studies that reported AT pO2 was inversely (16, 17) or not associated (12, 14) with insulin sensitivity. The reason for the differences between studies could be related to differences in the methods used to assess AT oxygen content, inadequate numbers of participants needed to detect statistically significant effects, and differences in the criteria used to define insulin-resistant and insulin-sensitive cohorts.
The purpose of the present study was to evaluate the potential impact of AT oxygenation on insulin sensitivity; on AT biology (markers of inflammation, fibrosis, and BCAA catabolism); and on systemic markers of inflammation and plasma BCAA concentrations in 3 distinct groups of participants, rigorously characterized as metabolically healthy lean (MHL), MHO, or MUO based on BMI, insulin sensitivity, plasma triglyceride (TG) concentration, and intrahepatic TG (IHTG) content. We hypothesized that AT pO2 would (a) progressively decline from the MHL to the MHO to the MUO group; (b) correlate positively with hepatic and whole-body insulin sensitivity; (c) correlate negatively with AT markers of inflammation and fibrosis and negatively with plasma markers of inflammation; and (d) correlate positively with AT markers of BCAA catabolism and negatively with plasma BCAA concentrations.
Subject flow, body composition, and metabolic characteristics. Among the 58 participants who enrolled in the study, 46 completed the study testing procedures (MHL n = 11, MHO n = 15, and MUO n = 20; Figure 1). The MHO and MUO groups were purposely matched on BMI and percentage body fat. There was no difference in subcutaneous abdominal adipose tissue (SAAT) volume between these 2 groups, but intra-abdominal adipose tissue (IAAT) volume and IHTG content were greater in the MUO than in the MHO group (Table 1). By design, IHTG content was similar in the MHO and MHL groups; however, IAAT volume was greater in the MHO than in the MHL group. Average adipocyte size in the MUO group was not different from that in the MHO group, but adipocyte size was larger in both the MHO and MUO groups than in the MHL group. Fasting plasma TG and glucose concentrations, hemoglobin A1c, and 2-hour oral glucose tolerance test (OGTT) plasma glucose concentration were greater in the MUO than in the MHO and MHL groups, with no difference between the MHO and MHL groups. Plasma insulin and the homeostasis model assessment of insulin resistance (HOMA-IR) progressively increased from the MHL to the MHO to the MUO group, whereas plasma HDL-cholesterol, hepatic insulin sensitivity index (HISI) assessed as the reciprocal of the product of basal endogenous glucose production rate and basal plasma insulin concentration, and whole-body insulin sensitivity assessed as the glucose rate of disposal relative to plasma insulin concentration (glucose Rd/I) during a hyperinsulinemic-euglycemic clamp procedure progressively decreased from the MHL to the MHO to the MUO group (Table 1). Participants in the MUO group were about 8 years older than participants in the MHL group. However, including participant age as a covariate in our statistical analyses did not change the statistical significance of any of our outcomes or alter any conclusions.
AT oxygen tension and HIF1A gene expression are affected by adiposity and positively associated with insulin sensitivity. AT pO2 decreased progressively from the MHL (64.7 ± 3.6 mmHg) to the MHO (56.7 ± 2.8 mmHg) to the MUO group (48.6 ± 1.9 mmHg) (Figure 2A). To assess whether the lower AT pO2 was associated with physiological evidence of hypoxia, we evaluated AT gene expression of hypoxia-inducible factor-1α (HIF-1α), a transcription factor that mediates cellular response to low oxygen tension (18). We found a progressive increase in HIF1A gene expression from the MHL to the MHO to the MUO group with greater gene expression in the MUO group than in the MHL group (Figure 2B) and an inverse association between AT pO2 and HIF1A gene expression (Supplemental Figure 1; supplemental material available online with this article; https://doi.org/10.1172/JCI141828DS1). In addition, AT pO2 among all participants was positively correlated with both hepatic and whole-body insulin sensitivity (Figure 2, C and D).
AT oxygen tension is reduced in MUO and is positively associated with insulin sensitivity. (A and B) AT oxygen partial pressure (AT pO2) (A) and HIF1A gene expression (B) in metabolically healthy lean (MHL; n = 11), metabolically healthy obese (MHO; n = 15), and metabolically unhealthy obese (MUO; n = 20) groups. Data are means ± SEM. One-way ANOVA and Fisher’s least significant difference procedure were used to identify significant mean differences between groups. *Value significantly different from the MHL value, P < 0.001. †Value significantly different from the MHO value, P < 0.05. #Linear trend, P < 0.001. (C and D) Relationship between AT pO2 and hepatic insulin sensitivity index (HISI), assessed as 1000 divided by the product of endogenous glucose rate of appearance (in μmol/kg fat-free mass/min) and plasma insulin concentration (in μU/mL) in the overnight fasting state (C), and between pO2 and whole-body insulin sensitivity, assessed as glucose rate of disposal (Rd; in nmol/kg fat-free mass/min) divided by plasma insulin (I) concentration (in μU/mL) during a hyperinsulinemic-euglycemic clamp procedure (D), in MHL (white circles; n = 11), MHO (gray circles; n = 15), and MUO (black circles; n = 20) participants. Associations between AT pO2 and HISI and between AT pO2 and glucose Rd/I were determined using Pearson’s correlation coefficient.
AT oxygen tension is positively correlated with AT BCAA metabolism and protein synthesis and negatively correlated with AT inflammation, fibrosis, and angiogenesis. An unbiased global assessment of AT gene expression in relation to AT pO2 was evaluated using RNA sequencing. We found that 1982 genes were significantly (P < 0.005) associated with AT pO2. Analysis of the 1207 genes positively associated with AT pO2 revealed a biological pathway involved in BCAA catabolism to be the most highly enriched (Figure 3A). The link between AT pO2 and BCAA metabolism was further supported by (a) higher fasting plasma concentrations of the BCAAs leucine, isoleucine, and valine in the MUO group than in the MHL and MHO groups (Figure 4A); (b) negative associations between plasma BCAA concentrations and AT pO2 (Figure 4B); (c) positive associations between AT pO2 and the expression of genes that encode key enzymes involved in catabolism of BCAAs (Figure 4C); and (d) negative associations between AT HIF1A expression and the expression of genes that encode key enzymes involved in BCAA catabolism (Supplemental Figure 2). The remaining pathways enriched with genes positively associated with AT pO2 involved the regulation of protein synthesis (Figure 3A). Expression of individual genes encoding translation initiation and elongation factors was also positively associated with AT pO2 (Figure 5).
AT biological pathways associated with AT oxygen tension. (A and B) RNA sequencing was conducted in AT samples obtained from MHL (n = 11), MHO (n = 14), and MUO (n = 20) participants. Functional enrichment analyses were performed to identify biological pathways significantly (false discovery rate < 0.01) enriched with genes positively (red) (A) and negatively (blue) (B) associated with AT oxygen tension.
AT oxygen tension is negatively associated with plasma BCAA concentrations and positively associated with AT expression of genes involved in BCAA catabolism. (A) Fasting plasma leucine, isoleucine, and valine concentrations in MHL (n = 11), MHO (n = 15), and MUO (n = 20) groups. Data are means ± SEM. One-way ANOVA and Fisher’s least significant difference procedure were used to identify significant differences between groups. *Value significantly different from the corresponding MHL value, P < 0.05. †Value significantly different from the corresponding MHO value, P < 0.05. #Significant linear trend for the 3 groups, P < 0.05. (B and C) Relationship between AT pO2 and plasma leucine, isoleucine, and valine concentrations (B) and AT gene expression of key enzymes involved in BCAA catabolism (C) in MHL (white circles; n = 11), MHO (gray circles; n = 15 in B and n = 14 in C), and MUO (black circles; n = 20) participants. Associations between AT pO2 and plasma BCAA concentrations and between AT pO2 and AT gene expression were determined using Pearson’s correlation coefficient.
AT oxygen tension is positively associated with AT expression of genes involved in regulating protein synthesis. Relationship between AT pO2 and expression of AT genes involved in regulating protein synthesis in MHL (white circles; n = 11), MHO (gray circles; n = 14), and MUO (black circles; n = 20) participants. Associations between AT pO2 and AT gene expression were determined using Pearson’s correlation coefficient.
Genes that were negatively associated with AT pO2 were highly enriched in multiple biological pathways involved in proinflammatory responses (Figure 3B). Moreover, AT pO2 was inversely correlated with the expression of the proinflammatory markers SERPINE1, CD68, CD74, CXCL16, CCL2, CCL5, TNF, and IL6 (Figure 6). In general, an increase in AT markers of inflammation did not translate to an increase in plasma concentrations of cytokines and chemokines (Table 2). Plasma concentrations of proinflammatory IFN-γ, TNF-α, IL-6, CCL-2, and CCL-5 were not greater in the MUO group than in the MHO and MHL groups, and plasma concentration of the antiinflammatory cytokine IL-10 was not different among the 3 groups. However, plasma concentration of PAI-1, which is the protein product of AT SERPINE1, was greater in the MUO than in the MHO and MHL groups and was inversely correlated with AT pO2 (Table 2).
AT oxygen tension is inversely associated with AT expression of genes related to inflammation. Relationship between AT pO2 and AT expression of genes related to inflammation in MHL (white circles; n = 11), MHO (gray circles; n = 14), and MUO (black circles; n = 20) participants. Associations between AT pO2 and AT gene expression were determined using Pearson’s correlation coefficient.
Genes negatively associated with AT pO2 were also enriched in pathways related to extracellular matrix (ECM) remodeling (Figure 3B), with expression of genes regulating the formation, cross-linking, and stabilization of the ECM inversely correlated with AT pO2 (Figure 7). In addition, genes negatively associated with AT pO2 were enriched in a pathway involved in regulating angiogenesis (Figure 3B), which was driven by genes regulating inflammation and ECM remodeling in this gene set (Supplemental Table 1). In contrast, VEGFA gene expression, which is a major regulator of angiogenesis (19, 20), progressively decreased from the MHL to the MHO to the MUO group and was positively associated with AT pO2 (Figure 8). Our VEGFA results are consistent with data from previous studies that found higher AT HIF1A gene expression but lower VEGFA expression in a rodent model of obesity (21) and lower AT VEGFA expression in people with obesity than in people who are lean (12, 16).
AT oxygen tension is inversely associated with AT expression of genes involved in ECM remodeling. Relationship between AT pO2 and AT expression of genes involved in ECM remodeling in MHL (white circles; n = 11), MHO (gray circles; n = 14), and MUO (black circles; n = 20) participants. Associations between AT pO2 and AT gene expression were determined using Pearson’s correlation coefficient.
AT oxygen tension is positively associated with AT VEGFA expression. (A) AT VEGFA gene expression in MHL (n = 11), MHO (n = 14), and MUO (n = 20) groups. Data are means ± SEM. One-way ANOVA and Fisher’s least significant difference procedure were used to identify significant mean differences between groups. *Value significantly different from the MHL value, P < 0.001. †Value significantly different from the MHO value, P < 0.05. #Linear trend, P < 0.001. (B) Relationship between AT pO2 and AT VEGFA gene expression in MHL (white circles; n = 11), MHO (gray circles; n = 14), and MUO (black circles; n = 20) participants. Associations between AT pO2 and AT gene expression were determined using Pearson’s correlation coefficient.
Studies conducted in rodent models have shown that AT hypoxia is involved in the pathogenesis of obesity-induced insulin resistance due to alterations in AT biological processes. In the present study, we evaluated whether decreased AT pO2 in people with obesity is also associated with alterations in AT biology and a concomitant impairment in insulin action. To this end, we evaluated subcutaneous abdominal AT pO2, hepatic and whole-body insulin sensitivity, AT markers of inflammation, fibrosis, and BCAA catabolism, and related circulating factors proposed to be involved in the pathogenesis of insulin resistance in 3 distinct groups of participants that were stratified by insulin sensitivity and adiposity. We found that AT pO2 progressively declined and AT expression of HIF1A progressively increased from the MHL to the MHO to the MUO group and positively correlated with hepatic and whole-body insulin sensitivity. In addition, AT pO2 was negatively associated with AT gene expression of markers of inflammation and ECM remodeling. Plasma PAI-1 concentration increased progressively from the MHL to the MHO to the MUO group and was negatively associated with AT pO2. In contrast, most plasma concentrations of cytokines and chemokines were not different among the MHL, MHO, and MUO groups and did not correlate with AT pO2. We also found positive associations between AT pO2 and the expression of genes involved in BCAA catabolism, in conjunction with an inverse relationship between AT pO2 and plasma BCAA concentrations. Together, these data demonstrate a link between AT oxygenation and AT inflammation, fibrogenesis, and BCAA catabolism, and suggest that decreased AT pO2 in people contributes to insulin resistance by increasing AT PAI-1 production (encoded by SERPINE1) and decreasing AT BCAA catabolism, thereby increasing plasma PAI-1 and BCAA concentrations.
The adequacy of adipocyte oxygenation is determined by the rate of arterial O2 delivery to AT, the ability of O2 to diffuse from AT capillaries to adipocytes, and the rate of adipocyte O2 consumption. People with obesity are therefore at increased risk of adipocyte hypoxia because (a) systemic arterial O2 content is usually decreased in people with obesity due to pulmonary dysfunction (22, 23); (b) AT blood flow is much lower in people with obesity than in people who are lean (16, 24–26); (c) AT capillary density decreases with excess adiposity (12, 13, 16); and (d) adipocyte hypertrophy increases the O2 diffusion distance from capillaries to adipocytes (27). However, data from studies conducted in vivo suggest that O2 consumption per gram of AT is lower in people with obesity than in those who are lean (16, 22). This decrease in O2 consumption is likely related to decreased AT energy requirements, because people with obesity have fewer, but larger lipid-laden adipocytes per gram of AT (28). We found that AT pO2 was highest in the MHL group and that AT pO2 was lower in people with MUO than MHO even though there was no difference in AT mass or adipocyte size between the groups with obesity. These findings suggest that decreased AT O2 delivery, increased O2 consumption, or both contribute to decreased AT pO2 in people with MUO compared with people with MHO.
We found a positive correlation between AT pO2 and whole-body insulin sensitivity with respect to glucose metabolism, which is inconsistent with results from previous studies that found that AT pO2 was negatively associated (16, 17) or not associated (12, 14) with insulin sensitivity. The discrepancy in findings among studies could be related to differences in the techniques used to measure AT pO2 (i.e., optochemical [refs. 16, 17] versus O2 electrode [refs. 12–15] methods), differences in the criteria used to define insulin-resistant and insulin-sensitive cohorts, and inadequate numbers of subjects to detect statistically significant effects. Nonetheless, our results are consistent with previous studies that found that (a) AT oxygenation was decreased in genetic and diet-induced rodent models of obesity and insulin resistance (29–31); (b) prolonged (16–24 hours) hypoxia causes insulin resistance in cultured adipocytes (30, 32); (c) experimentally induced reductions in AT O2 consumption improved glucose tolerance and insulin sensitivity in rodent models of obesity (5, 33); (d) hyperbaric oxygen therapy increases insulin sensitivity (34, 35); and (e) acute hypoxia causes insulin resistance in people (36).
Our study cannot determine whether the relationship we observed between AT pO2 and insulin sensitivity is causal or simply an association. Reduced AT oxygenation in rodent models of obesity has been proposed to cause systemic insulin resistance, at least in part, by inducing AT inflammation (increased gene expression of proinflammatory markers [refs. 5, 21, 29], macrophage infiltration [ref. 29], and polarization toward proinflammatory M1-like macrophages [refs. 37, 38]) and fibrosis (21, 39). We found an inverse relationship between AT oxygenation and the expression of genes involved in inflammation, consistent with data from studies conducted in cultured adipocytes that found that hypoxia increased the expression of genes involved in inflammation, including TNFA, IL6, and SERPINE1 (40, 41). However, with the exception of plasma PAI-1, we did not detect a corresponding increase in circulating proinflammatory cytokines and chemokines. Plasma PAI-1 concentrations progressively increased from the MHL to the MHO to the MUO group and were negatively associated with AT pO2. Increased plasma PAI-1 is associated with insulin resistance (42) and is an independent risk factor for developing T2D (43, 44). Moreover, whole-body and adipocyte-specific knockouts of PAI-1 improve insulin action (45, 46), whereas adipocyte-specific PAI-1 overexpression causes insulin resistance (47). It is also possible that other mediators associated with AT inflammation and macrophage infiltration, such as AT macrophage-derived exosomes that cause systemic insulin resistance in rodent models (48), also cause insulin resistance in people.
Increased plasma BCAA concentrations are associated with insulin resistance, prediabetes, and type 2 diabetes in people (49–52), and weight loss–induced improvements in insulin sensitivity are associated with a decrease in plasma BCAA concentrations (53, 54). It is likely that multiple mechanisms are responsible for BCAA-induced insulin resistance (55), including mammalian target of rapamycin–mediated (mTOR-mediated) inhibition of insulin-stimulated tyrosine phosphorylation of insulin receptor substrate 1 (IRS-1) (56, 57); production of 3-hydroxy-isobutyrate, a metabolite of valine that stimulates fatty acid uptake in skeletal muscle (58); and reduced hepatic secretion of FGF21 that has an insulin-sensitizing effect on whole-body glucose uptake (59). Plasma BCAAs are regulated in part by BCAA catabolism in AT. Gene and protein expression of enzymes involved in BCAA oxidation and BCAA oxidation rates is downregulated in AT of rodent models of obesity with high plasma BCAA concentrations, and transplantation of normal AT into mice that are globally defective in peripheral BCAA metabolism causes a marked decrease in plasma BCAA concentrations (8, 9). Moreover, results from studies conducted in cultured adipocytes demonstrate that hypoxia suppresses BCAA catabolism and oxidation (7, 11), providing a mechanistic link between AT oxygenation and BCAA metabolism. Consistent with these findings, we found that AT pO2 was positively associated with expression of genes that encode key enzymes involved in BCAA catabolism, and that plasma BCAA concentrations were inversely correlated with AT pO2. In addition, AT expression of pathways and specific genes that regulate transcription and translation was positively associated with AT pO2, which is consistent with data from studies conducted in cell culture that found that hypoxia inhibited the expression of genes involved in translation and protein synthesis (60–62). Together, these data suggest that AT oxygenation is involved in regulating AT amino acid and protein metabolism and that reduced AT pO2 contributes to systemic insulin resistance, in part, by reducing AT catabolism of BCAA and increasing circulating BCAAs.
It is important to consider whether the differences we observed in AT pO2 among our cohorts were adequate to induce a physiological response to hypoxia. We used oxygen-sensitive fiber optic probes, which measure O2 partial pressures in AT interstitial space but do not measure intracellular adipocyte pO2. This technical limitation is common to all studies conducted in people, which report AT extracellular rather than intracellular pO2 (12–17). Studies conducted in cell systems in vitro have shown that intracellular pO2 is lower than extracellular pO2, but the relative difference is maintained across a wide range of oxygen tensions (63). Therefore, our measurement of interstitial AT pO2 likely provides a reliable index of intra-adipocyte pO2. In addition, we assessed AT HIF1A gene expression as an index of physiological AT hypoxia. The cellular response to hypoxia in adipocytes is mediated through an increase in HIF-1α, a transcription factor that acts as a molecular oxygen sensor (64). Although considerable posttranslational modification of HIF-1α influences its activity (65), an increase in HIF1A gene expression is also a robust marker of AT hypoxia (31, 33, 66). We found a progressive increase in AT HIF1A gene expression from the MHL to the MHO to the MUO group. Moreover, AT pO2 was inversely associated with AT gene expression of several HIF-1α targets that cause AT inflammation and fibrosis induced by hypoxia (33, 67, 68), including PAI-1 (69), lysyl oxidase (21), and connective tissue growth factor (70). Accordingly, our data suggest that the reduction in AT pO2 observed in the MHO and MUO groups was sufficient to stimulate a physiological response to hypoxia.
In summary, we found a progressive decrease in AT oxygen tension and a progressive increase in HIF1A gene expression from people who are MHL to MHO to MUO. AT pO2 was negatively correlated with AT expression of genes involved in inflammation and ECM formation and positively correlated with AT expression of genes involved BCAA catabolism. In addition, AT pO2 was negatively correlated with both plasma PAI-1 and BCAA concentrations, and positively correlated with both hepatic and whole-body insulin sensitivity. Together, these data suggest that decreased AT pO2 contributes to insulin resistance in people with obesity by causing a constellation of alterations in AT biology that lead to an increase in plasma PAI-1 and BCAA concentrations (Figure 9). Therefore, differences in the decrease in AT pO2 associated with AT expansion provide a potential mechanism that contributes to the differences in insulin sensitivity between people with metabolically healthy and metabolically unhealthy obesity.
Proposed effects of decreased AT oxygen tension on AT biology and insulin sensitivity in people with obesity. An increase in AT mass associated with obesity decreases AT pO2, which increases AT HIF-1α expression and triggers a cascade of alterations in AT biology, including a decrease in BCAA catabolism and an increase in inflammation and fibrosis. These changes lead to an increase in circulating BCAAs and PAI-1 that impair systemic insulin action. Differences in AT pO2 among people with obesity provide a potential mechanism to help explain the differences in insulin sensitivity between people with metabolically healthy and metabolically unhealthy obesity.
Subjects
Forty-six men and women (age 37.9 ± 1.2 years old) participated in this study, which was conducted in the Clinical Translational Research Unit (CTRU) at Washington University School of Medicine in St. Louis (Missouri, USA) from May 2016 to April 2019. Participants were recruited using the Volunteers for Health database at Washington University School of Medicine and by local postings. Written, informed consent was obtained from all participants before their participation in this study, which was approved by the Human Research Protection Office of Washington University and registered in ClinicalTrials.gov (NCT02706262). All participants completed a screening evaluation that included a medical history and physical examination, standard blood tests, hemoglobin A1c (HbA1c), an oral glucose tolerance test (OGTT), and assessment of IHTG content by MRI. The following inclusion criteria were required for each of the 3 groups: (a) MHL (n = 11, 6 women) had BMI 18.5–24.9 kg/m2, plasma TG concentration <150 mg/dL, fasting plasma glucose concentration <100 mg/dL, 2-hour OGTT plasma glucose concentration <140 mg/dL, HbA1c ≤5.6%, and IHTG content <4%; (b) MHO (n = 15, 13 women) had BMI 30–49.9 kg/m2, plasma TG concentration <150 mg/dL, fasting plasma glucose concentration <100 mg/dL, 2-hour OGTT plasma glucose concentration <140 mg/dL, HbA1c ≤5.6%, and IHTG content <4%; and (c) MUO (n = 20, 15 women) had BMI 30–49.9 kg/m2, prediabetes (fasting plasma glucose concentration ≥100 mg/dL, 2-hour OGTT plasma glucose concentration ≥140 mg/dL, and/or HbA1c ≥5.7%), and IHTG content ≥5%. Potential participants who had a history of diabetes or liver disease other than nonalcoholic fatty liver disease, were taking medications that could affect the study outcome measures, or consumed excessive amounts of alcohol (>21 drinks per week for men and >14 drinks per week for women) were excluded. Three participants (1 person with MHO and 2 people with MUO) had a history of obstructive sleep apnea.
Experimental procedures
Body composition analyses. Body fat mass and fat-free mass (FFM) were determined by dual-energy x-ray absorptiometry (DXA; Lunar iDXA, GE Healthcare Lunar). Subcutaneous abdominal and intra-abdominal AT volumes and IHTG content were determined by MRI (3-T superconducting magnet, Siemens).
Hyperinsulinemic-euglycemic clamp procedure. Subjects were admitted to the inpatient unit of the Washington University CTRU at 1800 hours for about 48 hours. Subjects were given standard meals, each containing one-third of their estimated daily energy requirements (71), at 1900 hours on the day of admission (day 0) and at 0700 hours, 1300 hours, and 1900 hours on day 1. A hyperinsulinemic-euglycemic clamp procedure was performed on day 2, after subjects fasted for 11 hours overnight. At 0600 hours, catheters were inserted into an arm vein for the infusion of stable isotopically labeled glucose, insulin, and dextrose and into a radial artery for blood sampling. At 0700 hours, a primed (8.0 μmol/kg) continuous (0.08 μmol/kg/min) infusion of [U-13C]glucose (Cambridge Isotope Laboratories Inc.) was started and maintained for 210 minutes. At 1030 hours, insulin was infused for 210 minutes at a rate of 50 mU/m2 body surface area (BSA)/min (initiated with a 2-step priming dose of 200 mU/m2 BSA/min for 5 minutes followed by 100 mU/m2 BSA/min for 5 minutes). The infusion of [U-13C]glucose was stopped during insulin infusion because of the expected decrease in hepatic glucose production (72). Euglycemia (~100 mg/dL) was maintained by infusion of 20% dextrose, enriched to approximately 1% with [U-13C]glucose. Blood samples were obtained before the beginning of the tracer infusion to assess background plasma glucose tracer enrichment and plasma cytokine, chemokine, and BCAA concentrations, and every 6–7 minutes during the last 20 minutes of the basal and insulin infusion periods to assess glucose kinetics and plasma insulin concentrations; additional blood samples were obtained every 10 minutes throughout the clamp procedure to monitor blood glucose concentrations.
AT pO2 measurement and biopsies. AT pO2 was measured during the basal period of the clamp procedure. The periumbilical area was cleaned with alcohol and anesthetized with 1% lidocaine. A small (~0.5 cm) skin incision was made with a scalpel and a sterile, 250-μm-diameter, oxygen-sensitive fiber optic probe (OxyLite, Oxford Optronix Ltd.), which was inserted through the incision 3 cm into subcutaneous abdominal AT (SAAT). This probe measures interstitial pO2 at the tip of the probe and provides an accurate assessment of pO2 across the physiological range in AT (73, 74). The fiber optic probe was held in place for 2 minutes, and measurements of pO2 were obtained every 10 seconds; AT pO2 was calculated as the average pO2 recorded during steady-state conditions, defined as less than a 10% variance in pO2 for at least 40 seconds. An AT biopsy was then obtained with a conchotome forceps and placed in 10% formalin before being embedded in paraffin; additional SAAT was obtained by manual aspiration using a 4-mm liposuction cannula (Tulip Medical Products) connected to a 60-cc syringe. Aspirated tissue was immediately rinsed in ice-cold saline and frozen in liquid nitrogen before being stored at –80°C until further processing.
Sample analysis and calculations
Plasma substrate, insulin, and cytokine concentrations. Plasma glucose concentration was determined using the glucose oxidase method (YSI Inc.). Plasma insulin, HbA1c, and lipid profile were measured in the Washington University Core Laboratory for Clinical Studies, as previously described (75). The homeostasis model assessment of insulin resistance (HOMA-IR) was calculated as the product of fasting plasma insulin (in μU/mL) and glucose (in mmol/L) concentrations divided by 22.5 (76). Plasma BCAA concentrations were determined using liquid chromatography–tandem mass spectrometry as previously described (77). Plasma cytokine and chemokine concentrations were determined using commercially available magnetic bead suspension assays (MilliporeSigma and R&D Systems) and a Luminex 200 analyzer (Luminex Corp.).
Glucose kinetics. Plasma glucose tracer-to-tracee ratio (TTR) was determined by gas chromatography–mass spectrometry as previously described (78). Hepatic insulin sensitivity was calculated as the reciprocal of the product of basal plasma insulin concentration (in μU/mL) and basal endogenous glucose production rate (in μmol/kg FFM/min) determined by division of the glucose tracer infusion rate by the average plasma glucose TTR during the last 20 minutes of the basal period of the clamp procedure (79). The total glucose rate of disposal (Rd) during insulin infusion was assumed to be equal to the sum of the endogenous glucose rate of appearance into the bloodstream and the rate of infused glucose during the last 20 minutes of the clamp procedure. Whole-body insulin sensitivity was calculated as glucose Rd per kilogram FFM divided by the average plasma insulin concentration (glucose Rd/I) during the final 20 minutes of the clamp procedure (79).
Adipocyte size. Adipocyte areas were assessed by evaluation of histological sections of AT stained with H&E. Bright-field images were captured with a NanoZoomer whole-slide imaging system (Hamamatsu Photonics). Semiautomated segmentation of the images was performed with MIPAR Image Analysis Software version 3.0.3 (MIPAR Software LLC) with manual correction to obtain the size of 100 adipocytes.
AT RNA sequencing. Total RNA was isolated from frozen SAAT samples using QIAzol lysis reagent and an RNeasy Mini Kit (QIAGEN) in combination with an RNase-free DNase Set (QIAGEN) as we previously described (80). One SAAT sample from the MHO group was excluded from further analysis because of poor RNA quality attributable to RNA degradation. Library preparation on the remaining samples was performed with total RNA and cDNA fragments sequenced on an Illumina NovaSeq 6000 (Illumina) with all samples sequenced on the same instrument and analyzed during the same run to avoid batch effects. The counts per million (CPM) reads were used to select genes that were significantly (P < 0.005) positively and negatively associated with AT pO2 using Pearson’s correlation. Functional enrichment analyses were performed on genes significantly associated with AT pO2 using the Database for Annotation, Visualization and Integrated Discovery (DAVID) Bioinformatics Resources 6.8 (http://david.abcc.ncifcrf.gov/) with Gene Ontology biological process (GO BP) Direct category as we previously described (80). GO terms with a false discovery rate less than 0.01 were considered to be significantly enriched, and redundant GO terms were removed using the REVIGO program (http://revigo.irb.hr/). The RNA sequencing data reported in this paper were deposited in the NCBI’s Gene Expression Omnibus database (GEO GSE152991).
Statistics
One-way ANOVA was used to compare subject characteristics and outcome measures among MHL, MHO, and MUO groups, and Fisher’s least significant difference procedure was used to identify significant mean differences where appropriate. Polynomial contrasts were performed for selected variables to describe the linear trend from the MHL to the MHO to the MUO group. ANCOVA with age as a covariate was also used to assess the potential importance of the statistically significant difference in mean age between the MHL and MUO groups on the outcomes measures. Relationships between AT pO2 and hepatic and whole-body insulin sensitivity, AT gene expression, and plasma BCAA, cytokine, and chemokine concentrations were determined using Pearson’s correlation coefficient. All statistical tests were 2-sided, and a P value less than 0.05 was considered statistically significant. Data are reported as means ± SEM unless otherwise noted. Statistical analyses were performed with SPSS (version 25, IBM). Based on the interindividual variability in interstitial AT pO2 we previously reported (5), we estimated that 15 subjects in each group would be needed to detect between-group differences in AT pO2 of 10 mmHg using a 2-sided test with a power of 0.8 and an α value of 0.05. These computations were performed with G*Power 3.1.9.2 (81).
Study approval
Subjects provided written, informed consent before participating in this study, which was approved by the Human Research Protection Office at Washington University School of Medicine in St. Louis.
SCB, GIS, and JWB conducted the studies. VC, JY, DM, and JDW performed sample analyses. MJ and BWP supervised the sample analyses. VC, SCB, GIS, JY, JWB, BDK, and SK analyzed the data. VC, SCB, GIS, and SK performed the statistical analyses and wrote the manuscript. SK designed and supervised the studies and obtained funding for the work. SK is the guarantor of this work and, as such, had full access to all the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis. All authors critically reviewed and edited the manuscript.
The authors thank Jennifer Shew, Frieda Custodio, and Adewole Okunade for their technical assistance; Janet Winkelmann, Sally Torbitzky, and the nurses of the Clinical Translational Research Unit for their assistance in conducting the studies; and the study subjects for their participation. This study was supported by NIH grants K01DK109119, T32HL130357, K01DK116917, R01ES027595, P42ES010337, DK56341 (Nutrition Obesity Research Center), DK20579 (Diabetes Research Center), DK052574 (Digestive Disease Research Center), and UL1TR002345 (Clinical and Translational Science Award); NIH Shared Instrumentation Grants S10RR0227552, S10OD020025, and S10OD026929; and the Foundation for Barnes-Jewish Hospital.
Address correspondence to: Samuel Klein, Center for Human Nutrition, Washington University School of Medicine, 660 South Euclid Avenue, Campus Box 8031, St. Louis, Missouri 63110, USA. Phone: 314.362.8708; Email: sklein@wustl.edu.
SCB’s present address is: School of Biological and Health Systems Engineering, Arizona State University, Tempe, Arizona, USA.
Conflict of interest: SK has received grants from Janssen Pharmaceuticals, personal fees from Pfizer, and personal fees from Merck Sharp & Dohme Corp.
Copyright: © 2020, American Society for Clinical Investigation.
Reference information: J Clin Invest. 2020;130(12):6688–6699.https://doi.org/10.1172/JCI141828.
Adipose tissue oxygenation and metabolic health in obesity: Time to move from association to causation?
Gijs H. Goossens