The association between cardiorespiratory fitness, liver fat and insulin resistance in adults with or without type 2 diabetes: a cross-sectional analysis

Background Exercise-induced improvements in cardiorespiratory fitness (CRF) often coincide with improvements in insulin sensitivity and reductions in liver fat content. However, there are limited data concerning the relationship between CRF and liver fat content in adults with varying degrees of metabolic dysfunction. Methods The aim of this study was to examine the association between CRF, liver fat content, and insulin resistance in inactive adults with obesity and with or without type 2 diabetes (T2D), via cross-sectional analysis. CRF was determined via a graded exercise test. Liver fat content was assessed via proton magnetic resonance spectroscopy and insulin resistance was assessed via homeostatic model of insulin resistance (HOMA-IR). A partial correlation analysis, controlling for age and gender, was performed to determine the association between CRF, demographic, cardiometabolic, and anthropometric variables. Independent t tests were performed to compare cardiometabolic outcomes between participants with T2D and participants without T2D. Results Seventy-two adults (46% male) with a mean age of 49.28 ± 10.8 years, BMI of 34.69 ± 4.87 kg/m2, liver fat content of 8.37 ± 6.90%, HOMA-IR of 3.07 ± 2.33 and CRF of 21.52 ± 3.77 mL/kg/min participated in this study. CRF was inversely associated with liver fat content (r = − 0.28, p = 0.019) and HOMA-IR (r = − 0.40, p < 0.001). Participants with T2D had significantly higher liver fat content (+ 3.66%, p = 0.024) and HOMA-IR (+ 2.44, p < 0.001) than participants without T2D. Participants with T2D tended to have lower CRF than participants without T2D (− 1.5 ml/kg/min, p = 0.094). Conclusion CRF was inversely associated with liver fat content and insulin resistance. Participants with T2D had lower CRF than those without T2D, however, the difference was not statistically significant. Further longitudinal studies are required to elucidate the relationship between CRF and the progression of obesity-related diseases such as T2D. Registration: ACTRN12614001220651 (retrospectively registered on the 19th November 2014) and ACTRN12614000723684 (prospectively registered on the 8th July 2014). Supplementary Information The online version contains supplementary material available at 10.1186/s13102-021-00261-9.


Background
Obesity plays a significant role in the development of many chronic diseases such as hypertension, coronary heart disease, numerous cancers, and type 2 diabetes (T2D) [48]. While the relationship between obesity and increased cardiometabolic risk is well established [1], the location of adipose tissue, particularly in and/or around the liver, heart, muscles, and pancreas, also known as ectopic fat, is significantly more predictive of adverse health outcomes such as hypertension and hyperglycaemia [36]. Metabolic dysfunction-associated fatty liver disease (MAFLD), characterised by excessive liver fat (LF) content, is strongly associated with insulin resistance and is highly prevalent in individuals with T2D [41]. As T2D and MAFLD share a similar underlying pathological process, they often present alongside other disorders such as obesity, dyslipidaemia, and hypertension [49]. A known strategy for managing these risk factors is increased physical activity [23,43], which importantly leads to improved cardiorespiratory fitness (CRF) [8]. However, individuals with T2D and/or MAFLD reportedly experience disease-related aerobic impairments, which contribute to, or are further affected by, mitochondrial dysfunction, cardiac dysfunction, insulin resistance, and diastolic dysfunction -which may manifest as low CRF [6,44].
Low CRF is a well-established risk factor for all-cause morbidity and mortality [18] and increases an individual's risk for developing T2D [40]. Multiple studies have shown that individuals who engage in structured exercise can improve CRF, insulin sensitivity, and reduce ectopic fat [13,[29][30][31], however, it is unclear whether these improvements are mediated through changes in CRF or through simply undertaking regular exercise. Recent observational data have shown that while physical activity levels were not associated with LF [16], low CRF was strongly and independently associated with MAFLD prevalence [9]. However, these findings are yet to be confirmed in participants with varying T2D-status using goldstandard LF quantification techniques and maximal exercise tests to determine CRF.
Decoupling CRF from physical activity-related improvements has proven difficult due to the interrelated nature of the two measures [24]. It could therefore be assumed that individuals who have higher physical activity and CRF levels also partake in other healthy behaviours, which provide further protection from a variety of obesity-related complications [34]. Furthermore, it is becoming increasingly accepted that both modifiable ,such as physical activity levels, and non-modifiable factors, such as gene-specific variations, contribute to variations in CRF [5,33], with the relative genetic contribution to CRF reported to be~50% [3,33]. Consequently, it is important to explore the association between CRF and cardiometabolic risk in adults with similar levels of selfreported physical activity in order to decouple physical activity from CRF and better assess CRF-related benefits independently. There are limited studies which have investigated the relationship between CRF and insulin resistance in adults with and without T2D, and even fewer studies assessing the relationship between CRF and LF, as quantified via gold-standard proton magnetic spectroscopy ( 1 H-MRS). Therefore, the aim of this crosssectional study, was to determine the association between CRF, LF, insulin resistance, and other cardiometabolic outcomes in inactive adults with obesity, and with or without T2D. A secondary aim was to compare cardiometabolic and CRF differences between adults with T2D to those without T2D. It was hypothesised that CRF would be inversely associated with LF content and insulin resistance. It was also hypothesised that inactive adults with obesity and T2D would have significantly lower CRF than inactive adults with obesity but without T2D.

Participants
The participants of this study were recruited via electronic bulletins, clinical databases, and media advertisements between June 2011 and February 2019. Eligible volunteers were between the ages of 18-65 years, had a BMI ≥ 30.0 kg/m 2 , self-reported being physically inactive (exercising < 3 days/week) and/or not currently meeting physical activity guidelines [4]. All participants were screened by a medical practitioner prior to enrolment and were excluded if there was evidence of an unstable cardiac condition, uncontrolled hypertension, or uncontrolled blood glucose. The analysis included 72 volunteers from two larger intervention trials (ACTRN12614001220651 and ACTRN12614000723684) for which the results have been published elsewhere [14,31,46]. Eligible participants were screened via telephone interview and those who met the inclusion criteria and provided written informed consent, were enrolled in the study and were assessed at the University of Sydney (NSW, Australia). The study conformed to the ethical guidelines of the 1975 Declaration of Helsinki and the procedures were approved by the University of Sydney Human Research Ethics Committee.

Anthropometry and blood pressure assessment
Relevant anthropometric data were collected in accordance with international standards [38]. The height of the participants was measured via stadiometer (SECA model 220 Telescopic Height Rod, Hamburg, Germany). Participants body weight was measured using Tanita BC-418 Body Composition Analyzer (Tanita Corporation, Tokyo, Japan) to the nearest 0.1 kg and BMI (kg/m 2 ) was calculated. Waist circumference (WC) was measured (SECA Model 201, Hamburg, Germany) thrice horizontally, between the inferior margin of the ribs and the superior border of the iliac crest after expiration but before inspiration. After 10 to 15 mins of quiet sitting, systolic (SBP) and diastolic (DBP) blood pressures were measured manually on each arm with a sphygmomanometer (Welch Allyn ® 767 Series Aneroid; New York, USA) and the average of three measures was recorded from the arm which produced the highest SBP and DBP readings.

Biochemical parameters
Venous blood was collected after an overnight fast (> 10 h) for the purpose of determining fasting blood glucose (FBG), insulin, high-sensitivity C-reactive protein (CRP), total cholesterol (TC), triglycerides (TG), high-density lipoprotein cholesterol (HDL), lowdensity lipoprotein cholesterol (LDL), alanine aminotransferase (ALT), and aspartate aminotransferase (AST). Concentration of plasma free fatty acids (FFA) was measured using stored plasma. All biochemical and lipid assessments were completed by a private accredited laboratory. Insulin resistance was assessed via the homeostatic model of insulin resistance (HOMA-IR) [45].

Proton magnetic resonance spectroscopy ( 1 H-MRS)
LF% was measured via 1 H-MRS using a Phillips Intera 1.5 Telsa Achieva MRI system (Philips Medical Systems, Best, Netherlands). Spectral data were post-processed by an assessor (NAJ), who was blinded to participant details, using a magnetic resonance user interface software (jMRUI, version 5.2; www.jmrui.eu) [25,37]. LF% ≥ 5.5% was considered consistent with the presence of metabolic dysfunction-associated fatty liver disease (MAFLD) [39]. The LF% quantification methodology employed in this study is comprehensively detailed elsewhere [31].

Cardiorespiratory fitness
CRF was assessed via a graded maximal exercise test on an electronically-braked cycle ergometer (Lode Corival, Netherlands) under the supervision of an Accredited Exercise Physiologist [35]. All tests incorporated a threemin warm up at 35 W and 65 W for women and men, respectively, and workloads were incrementally adjusted by 25 W every 150 s until volitional fatigue as described elsewhere [15]. Heart rate, blood pressure, and ratings of perceived exertion (RPE) were obtained at each stage of exercise, with RPE measured using the Borg scale [2]. The test was terminated when the pedalling rate fell below 50 revolutions per min despite encouragement, or volitional fatigue. Peak work capacity (W peak ) was measured [20] and peak oxygen consumption (VO 2Peak ) estimated as described previously [10]. CRF was assessed within one-week of LF assessment.

Statistical analysis
Data were analysed using Statistical Package for the Social Sciences (SPSS version 24.0; IBM Corp., Armonk, NY, USA). All data are reported as the mean values ± standard deviation (SD) unless otherwise stated. Independent t tests were performed to compare differences in CRF and other biochemical and anthropometric measures between individuals with and without T2D. Homogeneity of variances was assessed via Levene's test for equality of variance and χ 2 for categorical data. Abnormally distributed primary outcome measures were transformed via natural logarithm prior to between-group comparisons. Partial correlation coefficients (r), controlling for age and sex, were used for correlations between continuous variables, and where one of the variables was dichotomous categorical data and the other was continuous. Associations between two categorical variables were assessed by χ 2 test for independence, whereby the φ coefficient determined the magnitude of the correlation. The magnitude of correlations were qualitatively assessed as: trivial (r < 0.1), small (r > 0.1 to 0.3), moderate (r > 0.3 to 0.4), strong (r > 0.5 to 0.7), very strong (r > 0.7 to 0.9), nearly perfect (r > 0.9), and perfect (r = 1.0) [11]. Male and female participants were divided into CRF quartiles, respectively, and analysis of variance was employed to determine significant differences in cardiometabolic outcomes between CRF quartiles (lowest fitness, IQR, highest fitness). Least significant difference post-hoc comparisons were used to identify and compare significant differences between CRF quartiles. Effect sizes were calculated as standardised difference in the means and expressed as Cohen's d. Statistical significance was set at p < 0.05. A two-tailed sensitivity analysis using the effect size of the difference in LF content between participants with T2D versus participants without T2D revealed that the study achieved 78% power (G-Power software; University of Trier, Trier, Germany).

Participants
Participant characteristics are summarised in Table 1. A total of 72 (33 male and 39 female) volunteers participated in this study. The mean age was 49.28 ± 10.18 years, BMI 34.69 ± 4.87 kg/m 2 , WC 108.55 ± 14.43 cm, LF was 8.37 ± 6.90%, and CRF was 21.52 ± 3.77 ml/kg/min. The prevalence of MAFLD was 43% in participants without T2D , 65% in participants with T2D, and 54% amongst all participants. Participants with T2D had significantly higher WC, LF%, HOMA-IR, FBG, ALT, TC, HDL, LDL, and FFA than those without T2D (p < 0.05 for all). There were no differences in other variables between participants with or without T2D. Male participants had significantly higher WC, TC, HDL, LDL, CRP, HOMA-IR and SBP than female participants (Supplementary Table 1).

Partial correlations
The associations between CRF, demographic, anthropometric, and cardiometabolic variables are summarised in Table 2

CRF and cardiometabolic outcomes
Participant characteristics stratified by CRF are summarised in Table 3. Participants with the lowest CRF had significantly higher LF% and HOMA-IR than those with the highest CRF (+ 5.31%, p = 0.021, + 2.50, p = 0.001, respectively) ( Fig. 1). When compared to participants with the highest CRF, those with the lowest CRF also had significantly higher BMI (+ 7.54 kg/m 2 , p < 0.

Discussion
This study is one of the first to show significant associations between key cardiometabolic risk factors, such as LF content, quantified via gold-standard 1 H-MRS, and CRF in individuals with and without T2D. The analyses showed that the prevalence of MAFLD in participants with T2D was Data presented as mean (SD). ES effect size, CI confidence interval, M male, F female, MAFLD metabolic dysfunction-associated fatty liver disease, Y yes, N no, BMI body mass index, LF% liver fat percentage, SBP systolic blood pressure, DBP diastolic blood pressure, AST aspartate aminotransferase, ALT alanine aminotransferase, CRP high-sensitivity C-reactive protein, FBG fasting blood glucose, HOMA-IR homeostatic model assessment of insulin resistance, HDL highdensity lipoprotein cholesterol, LDL low-density lipoprotein cholesterol, FFA free fatty acids, VO 2peak Peak oxygen consumption. a HOMA-IR and insulin measures reported for 37 participants with normal glucose tolerance and 29 participants with T2D as 6 participants with T2D were undergoing exogenous insulin therapy Table 2 Partial correlations between graded exercise test measured cardiorespiratory fitness, demographic, cardiometabolic and anthropometric variables  65%, which is consistent with previous reports [50]. When compared to participants without T2D, those with T2D had significantly higher cardiometabolic abnormalities. Furthermore, adults with T2D also had lower CRF than those without T2D, however, the difference was not statistically significant. Further analyses showed that even relatively small variations in CRF were associated with increased LF content and insulin resistance, which were significantly higher in individuals with relatively low CRF than those with relatively high CRF. Similarly, individuals with the highest CRF had lower cardiometabolic and inflammatory abnormalities than those with the lowest CRF.

The association between cardiorespiratory fitness and liver fat
In accordance with the original hypothesis, CRF was inversely associated with LF content. This finding supports previously reported data, where Kantartzis and colleagues showed that LF content was inversely associated with CRF (r = − 0.22, p = 0.005) in a mixed sample of adults with MAFLD or adults at risk of metabolic disease [12]. Similarly, the findings from The Young Finns study showed that for every 1 ml/kg/ min increase in CRF, the risk of MAFLD significantly declined (risk ratio = 0.90, 95% confidence interval 0.88 to 0.93; p < 0.001), however, LF content was measured up to three years after initial assessment of CRF and done so using ultrasound, which is less accurate than 1 H-MRS [27]. Another study reported low CRF was inversely associated with increasing MAFLD activity and steatohepatitis severity measured via liver biopsy and sampled up to four months after CRF assessment [19]. The results of the study reported herein add to previous findings and suggest that even relatively small variations in CRF are associated with increased LF in inactive adults with or at risk of T2D. Importantly, the Data presented as mean (SD). IQR interquartile range, ES effect size, CI confidence interval, M male, F female, MAFLD metabolic dysfunction-associated fatty liver disease, Y yes, N no, BMI body mass index, LF% liver fat percentage, SBP systolic blood pressure, DBP diastolic blood pressure, AST aspartate aminotransferase, ALT alanine aminotransferase, CRP high-sensitivity C-reactive protein, FBG fasting blood glucose, HOMA-IR homeostatic model assessment of insulin resistance, HDL high-density lipoprotein cholesterol, LDL low-density lipoprotein cholesterol, FFA free fatty acids, VO 2 Peak oxygen consumption. a HOMA-IR and insulin measures reported for 37 participants with normal glucose tolerance and 29 participants with T2D as 6 participants with T2D were undergoing exogenous insulin therapy. b significant difference between lowest and highest fitness; c significant difference between lowest fitness and IQR assessments of CRF and LF content were undertaken within a narrow timeframe (< 1 week) and LF% was quantified using gold-standard 1 H-MRS, thus highlighting the novelty and methodological rigour of the present study.

The association between cardiorespiratory fitness and insulin resistance and other cardiometabolic outcomes
In accordance with the original hypothesis, CRF was inversely associated with insulin resistance. Further analyses revealed that, on average, individuals with relatively poor fitness had abnormally high levels of insulin resistance, blood pressure, inflammation and LF content, whereas individuals with relatively high fitness, only had slightly abnormal LF content (11.40% vs. 6.08% for low vs. high CRF, respectively). These data, which are supported by previous findings [12,13,22,[29][30][31]47], suggest that CRF may play an important role in the context of metabolic disease. A recent meta-analysis reported that CRF was inversely associated with T2D prevalence in a dosedependent manner [28]. The results of the current study showed that while participants with T2D had lower CRF than those without T2D, the difference was not statistically significant. Furthermore, as physical activity has been shown to be inversely associated with LF content independent of BMI [17], only inactive participants were included in this study in an attempt to control for higher levels of physical activity -which incur cardiometabolic benefits. Because of this, the mean level of CRF of participants was quite low at 21.5 mL/kg/min. Reports show that CRF < 29.1 mL/kg/min increases the likelihood of developing metabolic syndrome six-fold [21], consequently a greater number of participants with higher levels of CRF are required to provide more robust results.

Mechanisms
While the mechanistic interplay between low CRF, MAFLD, and T2D remains unclear, it is purported that low CRF and ensuing mitochondrial defects contribute to the incomplete oxidation of fatty acids, which contribute to the accumulation of fatty acid by-products, such as ceramides and diacylglycerol in skeletal muscle and liver cells. The intracellular accumulation of these by-products impair insulin signalling pathways and contribute to insulin resistance [7,26]. Importantly, T2Drelated exercise intolerance appears to be reversed by regular exercise, which is made evident by the amelioration of skeletal muscle mitochondrial impairments, as well as improved insulin sensitivity, and CRF [29,32,42].

Strengths
This study adds to existing literature by highlighting the importance of CRF for metabolic health in adults with or at risk of metabolic disease. While previous studies have shown similar findings [12,27], this is the first study to do so involving a mixed sample of inactive adults with obesity and with or without T2D. Furthermore, CRF and LF were measured within a week of each other, whereas previous studies had measured CRF up to three years after LF assessment [27]. Finally, this study quantified LF% using 1 H-MRS, which is currently considered the gold-standard non-invasive technique for LF quantification.

Limitations
This study has limitations that should be considered when interpreting the results. Firstly, the results of this study, by nature, incorporated measures of CRF and LF content at a time-specific point and did not track the progression of any outcome to determine their relative importance in the development of MAFLD or T2D disease progression. Secondly, this study was completed using baseline data from previous interventional studies which prevented an a priori power analysis. However, the two-tailed sensitivity analysis revealed that the study achieved 78% power, which is just shy of the commonly accepted 80% value. Thirdly, this study assessed the amount of LF% per se and the methodology employed cannot determine the amount of fibrosis or classification of more severe liver diseases such as non-alcoholic steatohepatitis and/or their association with CRF. Additionally, while 1 H-MRS is currently considered the goldstandard non-invasive measurement technique for LF%, HOMA-IR is comparatively more limited and cannot provide inference into tissue-specific impairments in insulin sensitivity. Fourthly, although CRF was assessed using a validated graded exercise test model [10], the gold-standard of aerobic capacity testing involves direct measures of gas analysis, and where possible, this method should be implemented. Finally, while an attempt was made to control for high levels of physical activity by only recruiting individuals who reported to be inactive (exercising < 3 days/week), inter-participant variations in physical activity levels likely contributed to the associations between CRF and LF.

Conclusions
The results of this study showed that CRF was inversely associated with 1 H-MRS-quantified LF content in inactive adults with obesity. CRF was also inversely associated with insulin resistance and other key cardiometabolic risk factors. Adults with T2D had lower CRF than adults without T2D, however, the difference did not achieve statistical significance. These findings suggest that CRF may play a key role in metabolic dysfunction, however, further longitudinal studies are required to elucidate the relationship between CRF and the progression of obesity-related diseases such as MAFLD and T2D.