Abstract

Introduction. The mechanisms underlying the association between insulin resistance and intrahepatic lipid (IHL) accumulation are not completely understood. We sought to determine whether this association was explained by differences in fasting non-esterified fatty acid (NEFA) levels and/or NEFA suppression after oral glucose loading. Materials and Methods. We performed a cross-sectional analysis of 70 healthy participants in the Hertfordshire Physical Activity Trial (39 males, age 71.3 ± 2.4 years) who underwent oral glucose tolerance testing with glucose, insulin, and NEFA levels measured over two hours. IHL was quantified with magnetic resonance spectroscopy. Insulin sensitivity was measured with the oral glucose insulin sensitivity (OGIS) model, the leptin: adiponectin ratio (LAR), and the homeostasis model assessment (HOMA). Results. Measures of insulin sensitivity were not associated with fasting NEFA levels, but OGIS was strongly associated with NEFA suppression at 30 minutes and strongly inversely associated with IHL. Moreover, LAR was strongly inversely associated with NEFA suppression and strongly associated with IHL. This latter association (beta = 1.11 [1.01, 1.21], ) was explained by reduced NEFA suppression ( after adjustment). Conclusions. Impaired postprandial NEFA suppression, but not fasting NEFA, contributes to the strong and well-established association between whole body insulin resistance and liver fat accumulation.

1. Introduction

Excess intrahepatic lipid (IHL) accumulation (nonalcoholic fatty liver disease, NAFLD) is an important component of the spectrum of metabolic derangements associated with central obesity and insulin resistance [1]. Studies elucidating the mechanistic basis for these associations have suggested that increased circulating nonesterified fatty acids (NEFAs) lead to elevations in IHL [2] and that these NEFAs are generated predominantly from lipolysis in white adipose tissue [3]. Inhibition of glucose oxidation by fatty acids is known to be an important feature of the glucose-fatty acid (or Randle) cycle [4], such that increases in plasma NEFA levels during starvation or in the diabetic state result in greater fat oxidation at the expense of glucose oxidation. Although insulin-mediated suppression of circulating NEFA levels is a robust marker of adipocyte insulin sensitivity [5, 6], elevated fasting NEFA levels are generally considered to be the cause of hepatic fat accumulation in obese insulin-resistant states [7]. Furthermore, a recent study of 42 nonobese adults which measured postabsorptive fatty acid disposal with labeled palmitate found no association between this and whole body or peripheral insulin sensitivity [8]. A formal comparison of the strengths of the associations between fasting and suppressed NEFA levels is therefore warranted. Our first objective was to compare the strengths of the associations, if any, between whole body insulin sensitivity (as the exposure) and fasting NEFA levels versus NEFA suppression after oral glucose loading (at 30 and 60 minutes, as the “outcomes”). For this, we conducted a post hoc, cross-sectional analysis of metabolic and anthropometric data from a cohort of healthy older adults who participated in the Hertfordshire Physical Activity Trial (HPAT) [9].

NEFA suppression is impaired in type 2 diabetes [10] and NAFLD [11], but whether the well-established associations between insulin resistance, impaired NEFA suppression, and IHL persist in apparently healthy older individuals without prevalent NAFLD or diabetes is not known. Also, the extent to which impaired NEFA suppression modulates the accumulation of IHL has not previously been determined. Chronic inflammation within adipose tissue has been implicated in the pathogenesis of peripheral insulin resistance [12], and levels of fat-derived hormones (specifically the leptin : adiponectin ratio, LAR) have recently been shown to correlate well with clamp measures of whole body insulin sensitivity [13]. Therefore, our second objective was to estimate the associations between model-based measures of insulin resistance and impaired NEFA suppression as exposures and IHL as the outcome and to determine the extent to which the associations between insulin resistance and IHL were explained by variations in NEFA suppression.

2. Methods

The rationale and design for the Hertfordshire Physical Activity Trial (ISRCTN 60986572) have been described previously [9]. Data reported here relate to post hoc, cross-sectional analyses of volunteers’ anthropometric and metabolic characteristics at the time of their entry into the study. Each participant provided written informed consent. The original study protocol was approved by the Hertfordshire Research Ethics Committee (LREC ref. 05/Q0201/23).

Trial participants were recruited from the Hertfordshire Cohort Study, consisting of men and women born in Hertfordshire, UK, between 1931 and 39 and still residing there [14]. Specifically, those who were deemed to be potentially suitable by their general practitioner for inclusion in a supervised aerobic exercise programme and who lived within ten miles of the exercise facility were invited to participate, as described previously [9]. Those with known diabetes, untreated or unstable ischaemic heart disease, or any medical condition that would preclude participation in an exercise programme were excluded from the trial. However, participants with incident diabetes (diagnosed at the time of entry to this study) were included in these analyses. Recruits attended the clinical research facility after an overnight fast. Of 106 who attended the screening visit, six were deemed to be unsuitable for the study because of poor mobility, preexisting diabetes, symptoms or signs suggestive of untreated ischaemic heart disease, or a combination of these factors were excluded. Of the remaining 100, MR imaging and spectroscopy were not performed on 30 individuals who had claustrophobia, cardiac pacemakers, or metal implants. Thus, 70 participants who enrolled in the study had baseline liver spectroscopy measures performed and constitute the cohort described herein.

All measurements were undertaken by trained staff adhering to standard operating procedures. Weight was measured on a Tanita (Tokyo, Japan) scale and height with a Seca (Hamburg, Germany) wall-mounted stadiometer. Waist circumference was measured using a D-loop nonstretch fibreglass tape measure and defined as the midpoint between the lower costal margin and the level of the superior iliac crests. Blood pressure was measured with an oscillometric device (Omron, Kyoto, Japan) using the right arm, after participants were seated quietly for five minutes. A dual energy X-ray absorptiometry (DEXA) scan (Lunar Prodigy Advance, GE Healthcare, Bedford, UK) was used to measure lean mass and body fat percentage [15]. Magnetic resonance measures of intrahepatic lipid (IHL) and visceral adipose tissue (VAT) were conducted on a whole body Siemens 3T Tim Trio scanner (Erlangen, Germany), as described previously [9]. A questionnaire was used to quantify alcohol consumption in units per week.

A standard 75 g oral glucose tolerance test (OGTT) was performed. Fasting samples were taken for glucose, insulin, C-peptide, lipid profile, NEFA, leptin, and adiponectin. Glucose was measured using a hexokinase assay (Siemens, Frimley, UK). Insulin and C-peptide were measured using a fluorometric autoDELFIA immunoassay (PerkinElmer Life Sciences, Turku, Finland). NEFA levels were measured on a colorimetric assay (Roche Diagnostics, Burgess Hill, UK). Leptin and adiponectin were both measured using a DELFIA assay (R&D Systems Europe, Abingdon, UK). After ingestion of glucose, further samples were taken every 30 minutes over two hours. Samples for glucose and lipid profiles were processed immediately, while those for insulin, C-peptide, NEFA, leptin, and adiponectin were spun and frozen for subsequent batch analysis. All samples were processed in the same laboratory. The oral glucose insulin sensitivity (OGIS) model [16] was used to determine peripheral insulin sensitivity based on dynamic insulin and glucose responses during the OGTT, primarily mediated through insulin effects on muscle. Additionally, we used the leptin : adiponectin ratio [13] as an index of whole body insulin sensitivity and the homeostasis model assessment (HOMA) as an index of hepatic insulin sensitivity [17]. The degree of insulin-mediated NEFA suppression was determined by calculating the percentage reduction in NEFA levels 30 and 60 minutes after glucose loading. The area under the concentration curve (AUC) for NEFA during the OGTT was calculated with the trapezium rule.

The anthropometric and metabolic characteristics of the study participants were summarised using means and standard deviations. To estimate the association between each of these characteristics and IHL, linear regression was used with log(IHL) as the outcome, and each characteristic standardised to have mean 0 and variance 1. The models also included age, gender, alcohol consumption (units per week, self-reported), and, where relevant, MRI-derived visceral fat area. For each exposure, the beta coefficient and 95% confidence limits were exponentiated, giving a ratio of geometric mean IHL per standard deviation increase in the exposure. The associations between standardised measures of insulin sensitivity (OGIS, LAR, and HOMA) and IHL were estimated using the same method, with adjustment for age, gender, MRI-derived visceral fat area, and alcohol consumption. The potential confounding effect of NEFA suppression at 30 and 60 minutes was also explored.

3. Results

Of the 70 HPAT participants included in these analyses, 39 were men. Mean ± SD age was 71.3 ± 2.4 years. Systolic and diastolic blood pressures were 137 ± 18 and 75 ± 9 mmHg, respectively. The median IHL content was 3.6% (range 0.2–34.5%), while 42% of participants had IHL >5.5%, thus exceeding the arbitrary diagnostic threshold for NAFLD [18]. Three participants were found to have incident, asymptomatic type 2 diabetes based on OGTT results and were included in all analyses. Data relating to other anthropometric and metabolic characteristics of the cohort are summarised in Table 1. The associations between these characteristics and IHL are also shown in Table 1. These associations have been standardised in order to allow a comparison of their relative strengths. So, for example, each standard deviation rise in LAR was associated with a 49% increase in IHL, while each standard deviation rise in adiponectin was associated with a 48% reduction in IHL. Different measures of adiposity were positively associated with increased IHL, as expected. There were significant associations between each measure of insulin sensitivity and IHL, in the anticipated directions. NEFA suppression after oral glucose loading was inversely associated with IHL as shown in Table 1 and Figure 1. Results were similar for unadjusted analyses (data not shown).

To assess the relative strengths of the associations between insulin sensitivity and fasting as opposed to suppressed NEFA levels, we standardised OGIS, LAR, and HOMA. There were no significant associations between any indices of insulin sensitivity and fasting NEFA levels (Table 2), nor did HOMA correlate with measures of NEFA suppression. The associations between both OGIS and LAR and NEFA suppression were stronger for the 30-minute than the 60-minute values (Table 2) and were not significant for any of the 120-minute values (data not shown). OGIS and LAR had equivalent (though opposing) strengths of association with NEFA suppression at 30 minutes although the significance of the inverse associations at 60 minutes for LAR was strengthened after adjusting for age, gender, MR VAT, and alcohol consumption.

In order to determine the extent to which impaired NEFA suppression contributes to the observed inverse associations between insulin sensitivity (as exposure) and liver fat (as outcome), we compared these associations before and after adjusting for NEFA suppression (Table 3). After adjusting for NEFA suppression at 30 minutes, the positive association between LAR and IHL was attenuated and lost statistical significance. However, when OGIS or HOMA was treated as the insulin sensitivity measure, there was no attenuation of the association with IHL after adjusting for NEFA suppression. These results were similar after adjusting for NEFA suppression at 60 minutes, as shown. In order to avoid any confounding effects of medications known to influence lipid metabolism, we conducted subgroup analyses excluding those taking statin medications ( , 20%) and those taking beta-blockers ( , 17%), but these did not change our findings (data not shown).

4. Discussion

We found strong and consistent associations between OGIS and LAR (but not HOMA) and NEFA suppression. The absence of any associations with fasting NEFA levels is consistent with other recent observations [19] and suggests that fasting NEFA levels may be influenced by other factors such as catecholamine or growth hormone levels, whereas postprandial NEFA suppression is predominantly determined by insulin. NEFA levels are very promptly suppressed by insulin, so the observation that NEFA suppression at 30 rather than 60 or 120 minutes is a stronger marker of insulin sensitivity is not unexpected.

The similar strengths of associations for LAR and OGIS with NEFA suppression are potentially interesting. (They occurred in opposite directions because the former measures insulin resistance while the latter measures insulin sensitivity.) Given that NEFA suppression is mediated by adipose tissue insulin sensitivity [6], we anticipated that LAR might be more strongly associated with it than OGIS, the latter reflecting insulin-mediated glucose disposal primarily in skeletal muscle, but this was not the case. However, LAR has only previously been shown to correlate with whole body insulin sensitivity rather than fat tissue sensitivity specifically [13], and while leptin and adiponectin are derived exclusively from adipocytes rather than merely acting as markers of adipocyte inflammation, they may actively modulate insulin action in other tissues. The absence of an association between impaired NEFA suppression and HOMA (which is generally regarded as an index of hepatic insulin resistance) suggests that even though the liver is capable of disposal of free fatty acids [20], variations in hepatic insulin sensitivity have a relatively small impact on wholebody fatty acid disposal. Lastly, the finding that the inverse association between NEFAAUC and IHL was twice as strong as that for NEFA suppression at 30 or 60 minutes probably reflects the bidimensional nature of the AUC measure (see Figure 1), where fasting as well as postprandial NEFA levels have a multiplicative effect on this variable.

Our data indicate that excess body fat (particularly visceral fat), insulin resistance, and impaired suppression of NEFA levels after oral glucose loading are all directly associated with increased IHL in older individuals. In order to explore the mechanistic basis for the association between peripheral insulin resistance and NAFLD, we adjusted for NEFA suppression and found that this association was significantly attenuated for LAR but not for OGIS or HOMA. For these analyses, we also adjusted for central adiposity and chose the MRI-derived cross-sectional visceral fat area over other measures of fatness (BMI, body fat %, and waist circumference) because it was most strongly associated with liver fat content (Table 1). However, results were similar when these other measures of fatness were used (data not shown). So reduced suppression of fatty acids explains the association between LAR, but not OGIS or HOMA, and liver fat content. This suggests that while OGIS and LAR are both indices of whole body insulin sensitivity, adipocyte insulin resistance is reflected to a greater extent with LAR than with OGIS or HOMA. These results also suggest that there are other mechanisms apart from impaired NEFA suppression linking insulin resistance and NAFLD, such as dysregulated de novo lipogenesis and hepatocyte endoplasmic reticulum stress [21].

This study has a number of important strengths. Very detailed anthropometric and metabolic characterisation was conducted in each participant, and results were consistent across different measures of body fatness and insulin sensitivity. MR spectroscopy is the most robust noninvasive method for quantifying IHL. Rather than comparing categories of steatosis, body fatness, or glycaemic status, as many studies do, all our measures are continuous and represent the distributions in a relatively healthy cohort of older white participants. The study also has some limitations. It is a post-hoc analysis of data from a subgroup of individuals who were willing to participate in an exercise trial. All participants in the trial were White. Thus, our results may not be generalisable to all older people or those who would be less amenable to an exercise intervention, for whatever reason. Also, only 70% of people who enrolled in the trial had MR imaging at baseline, while others were unwilling or were too large for the scanner, which may have introduced bias.

Nonetheless, the 42% prevalence of NAFLD was higher than we anticipated, and there was a substantial level of metabolic disturbance in this cohort, particularly in relation to the number of individuals with abnormal glucose metabolism during the OGTT. It is important to note that these abnormalities were only detected through participation in the study and were not diagnosed prior to it (and so were not “prevalent” as such). All of these individuals volunteered to participate in a 12-week exercise intervention. Diabetes was one of several exclusion criteria. Nonetheless, three of the 70 individuals (4.3%) had newly diagnosed type 2 diabetes at entry to the study, while a further 27 (38.6%) had abnormal glucose metabolism. None of these had symptoms of hyperglycaemia, nor were they on treatment for it at the time of testing. Therefore, we felt it is appropriate to include them in the analysis. These participants were “apparently healthy,” and within a cohort of this age, a certain level of undiagnosed metabolic disease, be it diabetes or liver steatosis, is to be expected, so to be able to quantify this so precisely in the paper contributes to the novelty of our findings. So our observation too that impaired NEFA suppression does mediate the association between whole body insulin resistance, measured with LAR and liver steatosis.

5. Conclusions

In our experience, there is a widely held perception amongst scientists and clinicians in the field of metabolism that fasting NEFA levels are strongly positively associated with insulin resistance. However, in conducting a formal comparison of the relative strengths of the associations between fasting versus suppressed NEFA levels, we have confirmed that the degree of NEFA suppression is far more strongly associated with indices of insulin resistance, namely, OGIS and LAR, and is thus important from a clinical and pathophysiological point of view. We have been careful to take account of confounding factors such as age [22], sex [8], and other factors likely to influence these associations such as alcohol consumption. We believe that the relatively good health of the participants in this study makes the findings described above, particularly in relation to the high prevalence of NAFLD, even more novel and compelling.

Abbreviations

DEXA: Dual energy X-ray absorptiometry
HOMA: Homeostasis model assessment
HPAT: Hertfordshire Physical Activity Trial
IHL: Intrahepatic lipid
LAR: Leptin : adiponectin Ratio
NAFLD: Nonalcoholic fatty liver disease
OGIS: Oral glucose insulin sensitivity
VAT: Visceral adipose tissue.

Conflict of Interests

The authors declare that there is no conflict of interests associated with this paper.

Authors’ Contribution

Francis M. Finucane contributed to study design, data acquisition, analysis and interpretation, and writing the paper. Stephen J. Sharp contributed to study design, data analysis and interpretation, and writing the paper. Mensud Hatunic contributed to data interpretation and writing the paper. Alison Sleigh contributed to study design, data acquisition and analysis, and writing the paper. Ema De Lucia Rolfe contributed to study design, data acquisition and analysis, and writing the paper. Avan Aihie Sayer contributed to study design, data acquisition, and writing the paper. Cyrus Cooper contributed to study funding, design, data acquisition, and writing the paper. Simon J. Griffin contributed to study funding, design, data acquisition, analysis and interpretation, and writing the paper. David B. Savage contributed to study design, data acquisition, analysis and interpretation, and writing the paper. Nicholas J. Wareham contributed to study funding, design, data acquisition, analysis and interpretation, and writing the paper.

Funding

This work was funded by the Medical Research Council.

Acknowledgments

The authors are grateful to the staff of the Wolfson Brain Imaging Centre and the Wellcome Trust Clinical Research Facility in Addenbrooke’s Hospital and the Field Epidemiology and Data Management Teams from the MRC Epidemiology Unit for assistance with study measures and data management, respectively. They thank the Cambridge University Hospitals NHS Foundation Trust, Department of Clinical Biochemistry, the MRC Centre for Obesity and Related Disorders (CORD), and the NIHR Cambridge Biomedical Research Centre, Core Biochemical Assay Laboratory (CBAL) for carrying out the biochemical assays. They are also very grateful to the study participants for making this work possible.