Association between total body muscle percentage and prevalence of non-alcoholic fatty liver disease in Korean adults findings from an 18-year follow-up: a prospective cohort study
Article information
Abstract
Background
This study aimed to elucidate the association between total lean muscle mass and the incidence of non-alcoholic fatty liver disease (NAFLD) in the adult Korean population.
Methods
Utilizing data derived from the 18-year prospective cohort of the Korean Genome and Epidemiology Study, NAFLD was diagnosed via the hepatic steatosis index with an established cutoff value of 36. Lean muscle mass was assessed via bioelectrical impedance analysis and subsequently divided into tertiles. A generalized mixed model with a logit link was employed for repeated measures data analysis, accounting for potential confounders.
Results
Analysis encompassed 7,794 participants yielding 49,177 measurements. The findings revealed a markedly increased incidence of NAFLD in the lower tertiles of muscle mass, specifically, tertile 1 (odds ratio [OR], 20.65; 95% confidence interval [CI], 9.66–44.11) and tertile 2 (OR, 4.57; 95% CI, 2.11–9.91), in comparison to tertile 3. Age-dependent decreases in the OR were observed within the tertile 1 group, with ORs of 10.12 at age of 40 years and 4.96 at age of 80 years. Moreover, each 1%-point increment in total muscle mass corresponded with an estimated OR of 0.87 (95% CI, 0.82–0.93) for NAFLD resolution.
Conclusions
The study demonstrates a significant association between total muscle mass and NAFLD prevalence among Korean adults. Given the potential endocrine role of muscle mass in NAFLD pathogenesis, interventions aimed at enhancing muscle mass might serve as an effective public health strategy for mitigating NAFLD prevalence.
Introduction
Non-alcoholic fatty liver disease (NAFLD) comprises a spectrum of liver conditions occurring in the absence of significant alcohol consumption. These conditions range from non-alcoholic fatty liver to non-alcoholic steatohepatitis (NASH) and NASH cirrhosis, all typified by excessive hepatic fat accumulation [1]. The global prevalence rate of NAFLD, the most common chronic liver disease, is a striking 32% among adults. In the Republic of Korea, the prevalence has risen sharply, largely attributable to westernized dietary habits, diminished physical activity, and the ensuing increase in obesity and diabetes rates. Further, NAFLD has been identified as an independent risk factor for numerous systemic diseases, including cardiovascular disease, diabetes mellitus, metabolic syndrome, chronic kidney disease, and various malignancies [2-4], thereby elevating its public health and socioeconomic impact.
The role of skeletal muscle extends beyond locomotion, serving as a secretory organ that produces myokines, substances influencing other organs in an autocrine, paracrine, or endocrine hormone-like manner [5,6]. Specific myokines, such as interleukin 6 (IL-6), secreted from skeletal muscle, modulate the development of NAFLD by acting as endocrine organs that regulate various metabolic processes, including glucose metabolism in the liver or pancreas and lipolysis in adipose tissue [7].
Previous research has established a correlation between muscle mass and NAFLD, with a majority of such studies being cross-sectional [8]. Nevertheless, there is a growing call for longitudinal research to map this association over time, thereby yielding a more holistic understanding of the muscle mass-NAFLD interplay. Longitudinal studies are crucial for elucidating the dynamic nature of this relationship, which holds significant implications for NAFLD prevention and management. Therefore, our study aimed to examine the association between muscle mass and NAFLD using a Korean population cohort dataset.
Methods
Ethical statements: Prior to the survey, all participants were informed that they had been randomly chosen to participate in the Korean Genome and Epidemiology Study (KoGES) survey with the right to refuse to be involved in further analyses, and signed informed consents were obtained. We received the data in fully anonymized form. This study was conducted after obtaining an approval from the Institutional Review Board (IRB) of Yeungnam University Hospital, (IRB No: 2023-05-057). The study was carried out in accordance with the ethical standards of the Helsinki Declaration.
1. Study participants
This study utilized a subset of data derived from the Korean Genome and Epidemiology Study (KoGES), an ongoing initiative by the National Institute of Health (NIH) and the Korea Disease Control and Prevention Agency since 2001 [9]. Permission was secured to access and scrutinize laboratory results and participant questionnaires from the KoGES cohort. The dataset under analysis encompassed a subset of the KoGES cohort consisting of adults aged 40 to 70 years at baseline, residing in the South Korean cities of Ansan and Anseong. Biennial data collection commenced in 2001, with this study utilizing data from the ninth follow-up survey conducted during 2017 to 2018. It should be noted that participant attendance varied per examination phase, with an absence of data explaining non-participation in each phase or exclusion from specific assessments.
Initial participant selection for this study entailed excluding individuals who consumed over 20 g of alcohol per day, given that the NAFLD diagnosis was confined to non-drinkers or those consuming less than 20 g of alcohol daily [10]. Those with a prior cancer diagnosis were likewise excluded. Further exclusions comprised cases with missing data among the primary confounding variables considered in this study, including the variables instrumental to the operational definition of NAFLD and those absent from the bioelectrical impedance analysis (BIA) of body composition at each measurement time point.
2. Sample collection and laboratory analysis
Laboratory analysis results were procured and utilized from the NIH. Blood samples from the KoGES cohort were drawn from participants and collected in a serum separator tube and two ethylenediaminetetraacetic acid tubes. To optimize long-term storage, both serum and plasma fractions were prepared and divided into six to 10 vials, each containing 300 to 500 μL. All samples were then transported to the National Biobank of Korea for meticulous storage, ensuring their availability for future research applications. Venous blood samples were obtained from study participants after a fasting period of ≥12 hours. Fasting glucose levels were assessed using a glucose oxidase-based assay. Aspartate aminotransferase (AST) and alanine aminotransferase (ALT) activities were quantified using enzymatic methods (Advia 1650 system; Siemens, Tarrytown, NY, USA). For samples procured prior to September 2002, analyses of glucose, AST, and ALT were performed with a HITACHI 7600 (Hitachi, Tokyo, Japan), with data adjusted to account for systematic errors due to equipment transitions.
3. Assessment of body composition
Total muscle mass was assessed using multifrequency BIA (MF-BIA, InBody 3.0; Biospace, Seoul, Korea). Unlike traditional BIA methods, which often provide partial measurements and rely on estimations for whole body composition, MF-BIA posits that the human body comprises five interconnected cylinders, enabling direct impedance measurements from these bodily compartments. An eight-point tetrapolar tactile electrode system was employed to measure impedance at four specific frequencies (5, 50, 250, and 500 kHz) across five segments (right arm, left arm, trunk, right leg, and left leg). These readings were utilized to calculate segmental lean body mass by assessing the intracellular and extracellular water components of total body water. The total muscle weight (in kg) was then divided by the total body weight (in kg) to yield the total body muscle percentage (TBM%).
4. Survey and anthropometric assessment
In this study, information regarding smoking habits, medical history (including diabetes mellitus and cancer), and alcohol consumption was collected through structured questionnaires. The survey incorporated a categorization of smoking status into “current nonsmoker,” “former smoker,” and “current smoker.” For alcohol intake, the questionnaire inquired about the frequency of alcohol use and daily consumption. This latter measure was determined considering the volume and type of alcoholic beverages consumed, their respective alcohol content, and the specific gravity of ethanol. Anthropometric measurements, such as height and weight, were recorded by trained personnel, which enabled the computation of body mass index (BMI; kg/m2) using the formula: weight (kg) divided by the square of height (m).
5. Non-alcoholic fatty liver disease categorization based on laboratory data
NAFLD was diagnosed using the hepatic steatosis index (HSI), a diagnostic tool developed and validated specifically for the Korean population [11]. A standard cutoff value of 36 was applied in the HSI formula:
In this formula, diabetes was defined either by a self-reported diagnosis in the survey or by fasting glucose levels equal to or exceeding 126 mg/dL in the corresponding year. Participants in survey phases 2 and 3 were questioned, “Have you been diagnosed with diabetes since the previous examination?” Thus, the classification of individuals with diabetes necessitated consideration of responses from each phase and previous examinations. Participants consuming over 20 g of alcohol per day were excluded at the participant selection stage; hence those exceeding the cutoff value were classified as having NAFLD.
6. Statistical analysis
Initially, participant characteristics at each time point were delineated based on the presence or absence of NAFLD, focusing specifically on relevant and confounding variables. Preliminary analysis revealed a significant discrepancy in body weight between NAFLD and non-NAFLD groups, necessitating the use of a weight-adjusted indicator for data analysis. TBM% was categorized by computing the tertile cutoff for each sex, derived from baseline measurements. Consequently, muscle mass at every time point was categorized into tertiles using these cutoff values.
Given the longitudinal nature of the study, repeated measures, and potential missing data, along with the dynamic variability of NAFLD prevalence and TBM% over time, a generalized linear mixed model was applied for multivariable analysis. This model selection ensured the accommodation of potential non-independence of measurements within individuals. With the outcome variable being binary, representing NAFLD presence or absence, a logit link function was employed. Random effects were incorporated for individuals to account for individual-specific repeated measurements, including random interactions. An adaptive Gauss-Hermite quadrature approximation was not utilized for model estimation. Additionally, an interaction term was included between TBM% and age to explore the variance in the relationship between TBM% and NAFLD according to age. Consideration of potential confounding variables [12]—age, sex, smoking status (nonsmoker, former smoker, current smoker), diabetes mellitus (either self-reported or fasting glucose, ≥126 mg/dL), and obesity (BMI, ≥25 kg/m2)—was applied, as these could possibly influence both the independent and dependent variables and are not in the causal pathway between them.
To further explore the potential impact of muscle development on the resolution of NAFLD, given the longitudinal nature of the study, we sought to determine the odds ratio (OR) for NAFLD in subsequent evaluations among individuals initially diagnosed with NAFLD, based on the change in TBM% from baseline. A model was constructed using the discrepancy between baseline and follow-up TBM% as an independent variable, excluding the initial data. Generally, a declining trend in the number of participants was observed in tertile 3 over each assessment period, potentially indicative of aging effects seen in longitudinally tracked subjects over an extended timeframe, using the baseline as the reference point. To adjust for this effect, age at each evaluation, rather than baseline age, was accounted for as a confounding variable in the mixed-effects model. While not the primary focus of the investigation, the results for total fat mass proportion (total fat/weight, %), total muscle adjusted by height (total muscle/height2), total muscle mass (kg), and waist circumference (cm) were also reported, utilizing the same statistical methodology. The results were presented as ORs, with statistical significance established at p<0.05.
All statistical analysis was performed using R software version 4.3.0 (R Foundation for Statistical Computing, Vienna, Austria; https://r-project.org), incorporating the lme4 package for the generalized linear mixed model analysis.
Results
1. General characteristics of participants
From an initial pool of 10,030 participants in the KoGES cohort at baseline, the study excluded 1,567 individuals identified as high-risk drinkers (consuming >20 g/day of alcohol), 237 individuals with a history of cancer, and 432 individuals with incomplete data relevant to the study. Consequently, the study encompassed a total of 7,794 individuals. Due to attrition, the number of measured cases varied over time: 6,607 in phase 2, 5,796 in phase 3, 5,224 in phase 4, 5,229 in phase 5, 4,899 in phase 6, 4,608 in phase 7, 4,627 in phase 8, and 4,393 in phase 9 (Fig. 1).
Table 1 provides the cohort’s baseline characteristics. At baseline, the NAFLD and non-NAFLD groups exhibited no significant age difference, with mean ages of 52.7±9.1 and 52.5±8.5 years, respectively. A higher female representation was observed in the NAFLD group, comprising 1,349 individuals (67.3%), compared to 3,433 individuals (59.3%) in the non-NAFLD group. Notably, the prevalence of obesity significantly varied across groups, with 1,550 individuals (26.8%) in the non-NAFLD group versus 1,821 (90.9%) in the NAFLD group.
For the purposes of categorizing TBM% into three groups, we established cutoff points based on the baseline (2001–2002) tertiles of TBM% in the study’s participants, separated by sex. The tertile cutoff values for males were 72.1% and 76.4%, while for females, they were 62.2% and 66.0%. Among the NAFLD group, the respective proportions for each tertile were 1,373 (68.5%), 500 (25.0%), and 130 (6.5%), which were significantly different from those of the non-NAFLD group, with values of 1,223 (21.1%), 2,092 (36.1%), and 2,476 (42.8%).
The NAFLD group demonstrated higher AST and ALT levels, with averages of 31.1±22.4 U/L and 39.0±45.8 U/L, respectively, compared to 27.8±15.8 U/L and 22.7±14.8 U/L respectively in the non-NAFLD group. Detailed results and statistical metrics for variables in relation to NAFLD status at each time point, including TBM% changes (% point), are provided in Supplementary Table 1.
2. Association between total body muscle percentage and non-alcoholic fatty liver disease
The results from the generalized linear mixed model analysis revealed significant correlations between TBM% and NAFLD presence when categorized based on tertile cutoffs. Relative to individuals in tertile 3, those in tertile 1 demonstrated an OR for NAFLD of 20.65 (95% CI, 9.66–44.11), suggesting a substantial risk elevation. Similarly, the OR was 4.57 (95% CI, 2.11–9.91) for individuals in tertile 2 compared to tertile 3, implying another significant risk increase. Concerning age, NAFLD saw a significant annual decrease by 0.97 (95% CI, 0.95–0.98). The interaction between TBM% in tertile 1 and age yielded a statistically significant result of 0.98 (95% CI, 0.97–1.00; p<0.01). However, the interaction between age and tertile 2 was not statistically significant (p=0.06) (Table 2). Accounting for these interaction terms, individuals in tertile 1 had an OR for NAFLD of 10.12 (95% CI, 7.85–13.05) at age of 40 years, 8.47 (95% CI, 7.21–9.95) at age of 50 years, 7.09 (95% CI, 6.72–8.92) at age of 60 years, 5.93 (95% CI, 4.68–7.51) at age of 70 years, and 4.96 (95% CI, 3.49–7.06) at age of 80 years, relative to those in tertile 3 (Fig. 2).
In the model examining the influence of TBM% alterations on NAFLD resolution among initial NAFLD cases, a 1%-point increase in TBM% corresponded to an OR of 0.87 (95% CI, 0.82–0.93) for NAFLD. Yet, the interaction between TBM% at different time points was not statistically significant (Table 3).
Models incorporating other body composition indices are provided exclusively Supplementary Table 2, presenting the adjusted ORs for each index. Total muscle mass, total muscle adjusted by height, and fat proportion exhibited positive associations with NAFLD prevalence.
Discussion
In this study, we scrutinized 18 years’ worth of cohort data to discern the relationship between TBM% and NAFLD prevalence. Individuals grouped within the highest tertile (tertile 3) of TBM% presented with notably reduced NAFLD prevalence when contrasted with those in tertiles 1 and 2. The effect size, represented by the discrepancy in NAFLD prevalence between tertiles 1 and 3, exhibited a diminishing trend as age increased. Yet, even at the advanced age of 80 years, the effect size remained substantial and statistically significant with an OR of 4.96. Conversely, the NAFLD prevalence difference between tertile 2 and tertile 3 did not demonstrate notable age-related variations. With an OR of 4.57, the effect size for the second tertile was generally lower relative to tertile 3, thereby suggesting a dose-dependent reduction in NAFLD risk with increased muscle mass, notwithstanding the influence of aging. Furthermore, the enhancement of muscle mass among individuals with NAFLD was correlated with NAFLD resolution, irrespective of TBM%. However, when assessing factors such as total muscle mass and muscle mass adjusted by height—factors that do not account for body weight—conflicting findings emerged. These inconsistencies might be attributed to the substantial weight differential, with the NAFLD group displaying approximately 10 kg more weight compared to their non-NAFLD counterparts. Therefore, it appears appropriate to employ the weight-adjusted indicator TBM% in such analyses.
Skeletal muscle, often perceived as an endocrine and paracrine organ, secretes myokines which exert effects on various organs and systems, including the liver, adipose tissue, pancreas, bones, and cardiovascular system. Moreover, skeletal muscle, serving as a primary organ responsible for insulin-mediated glucose uptake, has a close association with insulin resistance [13]. The condition of reduced skeletal muscle mass, known as sarcopenia, has been linked to metabolic syndrome [14], cardiovascular disease [15], diabetes mellitus [16], and even increased overall mortality [17], as reported in several systematic review studies.
Sarcopenia shares several risk factors with NAFLD and other metabolic disorders, such as obesity and metabolic syndrome, albeit via distinct mechanisms. Pathways involving systemic inflammation and insulin resistance are implicated in both disorders, with a decrease in skeletal muscle mass contributing to decreased insulin sensitivity and energy expenditure, thereby promoting increased hepatic gluconeogenesis, elevated hepatic free fatty acid uptake, and reduced fatty acid oxidation [18]. Moreover, the suppression of muscle peroxisome proliferator-activated receptor γ coactivator-1 α [19], IL-6, irisin [20], along with other molecules secreted by skeletal muscle, is associated with heightened risk of metabolic disorders and NAFLD [21]. However, potential pathways through which NAFLD impacts skeletal muscle have been proposed [22], and in epidemiological studies seeking to ascertain causality direction, a clear consensus remains elusive [23]. Some schools of thought suggest an association rooted in a common etiological basis [24].
Epidemiological research has thoroughly investigated the correlation between muscle mass and NAFLD or NASH. Past studies have examined the association between NAFLD or NASH and reduced muscle mass within the context of sarcopenia as well as reported the relationship between skeletal mass index and NAFLD prevalence. A 2020 systematic review, incorporating 19 studies for meta-analysis, demonstrated significantly lower skeletal mass index in NAFLD patients and a significant relationship between sarcopenia and increased prevalence of NAFLD, NASH, and NAFLD-related fibrosis [8]. Although the majority of studies examining the NAFLD and sarcopenia relationship, including those in the systematic review, were of cross-sectional design, some longitudinal studies have also been reported. Notably, Kim et al. [25] analyzed a cohort dataset of 12,624 individuals recruited from a medical institution and found the highest tertile of skeletal muscle mass to significantly decrease future NAFLD incidence (adjusted hazard ratio [aHR], 0.44; 95% CI, 0.38–0.51) in the no-NAFLD baseline group, while increasing the future NAFLD resolution rate in the NAFLD baseline group (aHR, 2.09; 95% CI, 1.02–4.28). In a separate study, Lee et al. [26] analyzed a cohort of individuals undergoing health screenings at a different hospital, identified the presence or absence of NAFLD via ultrasound, and discovered that individuals without NAFLD had an increased risk of developing NAFLD over a 10-year period based on muscle and adipose tissue mass levels. The study revealed that higher appendicular skeletal muscle mass in non-obese individuals was associated with a lower prevalence of NAFLD (OR, 1.38; 95% CI, 1.34–2.45).
Our study, in contrast to prior research, holds the advantage of sampling a broad population rather than a patient group from medical facilities, thus mitigating selection bias. Moreover, we leveraged an expansive dataset collected over an 18-year duration. The use of mixed models enabled the inclusion of individuals diagnosed with NAFLD at baseline while accommodating the dynamic fluctuations in NAFLD and muscle mass throughout the follow-up period. The study also considered age-associated changes in the relationship. As age increases, the diminishing OR of NAFLD in relation to TBM% could be ascribed to the increased influence of age-correlated NAFLD risk factors such as dyslipidemia [27] and altered thyroid function [28]. These conditions may exert a heightened impact on NAFLD development among the elderly, thereby weakening the correlation between TBM% and NAFLD. Nevertheless, it is significant that a noteworthy OR was observed even among older subjects. This finding reinforces that the relationship between TBM% and NAFLD remains robust within the aging population, suggesting that muscle mass continues to play a pivotal role in NAFLD pathogenesis, notwithstanding the amplified effect of age-related risk factors.
Additionally, our findings underscore that temporal alterations in muscle mass could modulate NAFLD resolution. This constitutes a strength of our study, given its longitudinal structure and the model’s capacity to account for changes in both TBM% and NAFLD status. Regardless of TBM% per se, shifts in muscle mass were linked to a decrease in NAFLD prevalence. In light of previously understood mechanisms, it can be inferred that an individual-level increase in muscle mass may lessen the NAFLD risk, thereby suggesting the potential efficacy of public health initiatives aimed at augmenting muscle mass for NAFLD reduction.
A limitation of this study lies in the diagnostic modality employed for NAFLD categorization. While we utilized the HSI index, a validated tool for NAFLD diagnosis among Koreans, we did not use direct diagnostic tools such as ultrasound. Additionally, the usage of a secondary dataset not primarily intended for NAFLD research could introduce potential misclassification, as it may not definitively exclude hepatic conditions resulting from medication use, viral infection, autoimmune disease, among others. Nevertheless, given the availability of data on alcohol consumption, we adopted a more stringent criterion to minimize misclassification. Diverse criteria have been applied to distinguish NAFLD from alcoholic fatty liver disease based on alcohol consumption, with some studies implementing thresholds of 30 g/day [29] or 40 g/day [30]. In our study, we opted for a relatively strict threshold of 20 g/day [10], aiming to enhance the specificity of the operational definition of NAFLD at the cost of potentially decreased sensitivity. Another limitation pertains to our usage of TBM% as a surrogate measure, instead of skeletal muscle, which bears clinical significance for the diagnosis of sarcopenia. Furthermore, muscle mass measurements were not obtained through techniques such as dual-energy X-ray absorptiometry, computed tomography, or magnetic resonance imaging, and the lack of control over factors like hydration in muscle measurement using BIA presents another limitation. However, it is noteworthy that TBM% estimated using BIA is known to correlate well with dual-energy X-ray absorptiometry results [31]. In order to mitigate potential overestimation issues with BIA, our study analysis was based on tertiles of muscle mass rather than muscle weight or proportion, thereby reducing the potential misclassification of participants due to differences in muscle measurement methods. Therefore, the impact of misclassification based on muscle measurement methods is likely minimal in our analysis. Another limitation to consider is that while the study elucidates the association between NAFLD and TBM%, it does not establish a causal relationship in which muscle mass precedes NAFLD.
Notwithstanding its limitations, this study underscores that TBM% is an important independent factor associated with NAFLD among the Korean adult demographic. This highlights the integral role muscle tissue plays in managing NAFLD within this population. Considering the presumptive causality suggested by muscle in light of previous mechanistic investigations, alongside the association of muscle mass alterations with NAFLD resolution, interventions directed toward muscle mass enhancement could substantially impact NAFLD reduction. Owing to the relative ease of quantifying muscle mass via BIA devices, this factor can be leveraged in the realm of public health to manage NAFLD, through venues such as nationwide health examinations or health promotion initiatives.
Notes
Conflicts of interest
No potential conflict of interest relevant to this article was reported.
Acknowledgments
Data in this study were from the Korean Genome and Epidemiology Study (KoGES; 4851-302), National Institute of Health, Korea Disease Control and Prevention Agency, Republic of Korea.
Funding
None.
Author contributions
Conceptualization: BCA, JHH, KOB; Data curation: CYP; Formal analysis: JHH, KOB; Methodology: CYP, KOB; Project administration, Supervision: KOB; Visualization, Resources: JHH; Investigation: BCA; Writing-original draft: BCA; Writing-review & editing: CYP, JHH, KOB.
Supplementary materials
Supplementary Tables 1 and 2 can be found via https://doi.org/10.12701/jyms.2023.00605.