Hippokratia 2018, 22(1):10-16
Eremić-Kojić N1,5, Đerić M 1,5, Govorčin ML2,5, Balać D 3,5, Kresoja M 4, Kojić-Damjanov S1,5
1Centre for Laboratory Medicine, 2Centre for Radiology, Clinical Centre of Vojvodina, 3Centre for Hygiene and Human Ecology, Institute of Public Health of Vojvodina, 4Department of Mathematics and Informatics, Faculty of Sciences, University of Novi Sad, 5Faculty of Medicine, University of Novi Sad, Novi Sad, Serbia
Background: In order to optimize the identification of persons with non-alcoholic fatty liver disease (NAFLD), several algorithms for hepatic steatosis were developed. These available algorithms, as well as an algorithm, derived using biochemical and anthropometric data of our participants, are compared in a cross-sectional pilot study.
Material and methods: We included 77 participants with abdominal obesity: 43 with NAFLD and 33 without NAFLD. Body mass index (BMI), waist circumference (WC) and hip circumference (HC), systolic and diastolic blood pressure were assessed. Fibrinogen, high sensitive C-reactive protein (hsCRP), aspartate aminotransferase (AST), alanine transaminase (ALT), gamma-glutamyl transferase (GGT), uric acid, ferritin, glucose, insulin, homocysteine, lipid status parameters, apolipoprotein A-I, apolipoprotein B and Lp(a)-lipoprotein were measured. Fatty liver was assessed by ultrasound with the presence or absence of hepatic steatosis. Discovering the most significant factor in the presence of NAFLD is assessed through logistic regression modeling. The predictor variables were chosen according to an algorithm derived from conducted factor analysis and other available algorithms for hepatic steatosis.
Results: Participants with NAFLD had significantly higher BMI (34.38 ± 9.73 vs 28.05 ± 4.79 kg/m2, p =0.001), WC (108.05 ± 11.47 vs 96.15 ± 14.27 cm, p =0.001), HC (114.93 ± 11.01 vs 108.21 ± 9.82 cm, p =0.050), systolic (128.98 ± 8.67 vs 122.42 ± 10.62 mmHg, p =0.010) and diastolic blood pressure (83.64 ± 5.94 vs 78.33 ± 7.57 mmHg, p =0.001), AST (23.93 ± 6.91 vs 21.70 ± 5.21 U/L, p =0.014), ALT (30.50 ± 13.70 vs 23.00 ± 11.75 U/L, p =0.007), hsCRP (4.34 ± 5.56 vs 2.98 ± 2.34mg/l, p =0.004) and uric acid (358.02 ± 83.29 vs 296.78 ± 84.54µmol/l, p =0.001), in comparison non NAFLD. Logistic regression model with algorithm derived from factor analysis showed the best performance. From other available algorithms, only fatty liver index (FLI) and hepatic steatosis index (HSI) had statistically significant discriminatory power.
Conclusions: Elevation of WC, HC, BMI, DBP, SBP, Fbg, hsCRP, glucose, and uric acid, incorporated in our hepatic steatosis prediction model, had the best predictive power among all assessed algorithms. HIPPOKRATIA 2018, 22(1): 10-16.
Key words: Non-alcoholic fatty liver disease, NAFLD, algorithms, body mass index, abdominal obesity
Corresponding author: Eremić-Kojić Nevena, MD, Department of Pathophysiology and Laboratory Medicine, Faculty of Medicine, University of Novi Sad, 3 Hajduk Veljkova Street, 21000 Novi Sad, Serbia, tel/fax: +38121525289, +381641333433, e-mail: email@example.com
Non-alcoholic fatty liver disease (NAFLD) has become popular in the recent years because of its high prevalence rate of more than 30 % in the general population and up to 70 % in the high-risk groups, such as morbidly obese individuals and patients with type 2 diabetes mellitus (T2DM)1.
NAFLD represents a broad spectrum of disease from isolated steatosis (NAFL: non-alcoholic fatty liver) up to a more severe form, known as non-alcoholic steatohepatitis (NASH) which is associated with lobular inflammation, ballooning, and progression to fibrosis and cirrhosis2. According to the American Association for the Study of Liver Disease (AASLD) NAFLD is evidence of hepatic steatosis in more than 5 % of hepatocytes; either by imaging or histology in the absence of other causes of secondary hepatic fat accumulation, such as significant alcohol consumption, steatogenic medication, viral hepatitis or hereditary disorders3.
NAFLD is strongly associated with progressive hepatic disease, but also with cardiometabolic disorders. It may accelerate the development of T2DM and cardiovascular disease4. Therefore, it is vital to identify patients with hepatic steatosis for early risk intervention. Although liver biopsy is considered to be the gold standard for diagnosis of NAFLD, it is invasive and therefore not applicable as a screening method in the routine clinical practice. Proton magnetic resonance spectroscopy (1H-MRS) is considered the most accurate non-invasive method for measuring liver fat content5. However, it is only available in highly specialized institutions due to its high cost. Routine ultrasound is also being used to diagnose NAFLD, but it has its limitations. Ultrasound is subjective, and operator dependent, which has a low sensitivity for the detection of mild steatosis and there may be technical issues in morbidly obese patients. Therefore, there has been considerable interest in blood markers alone or in combination with clinical parameters which might help in identifying patients with NAFLD2. Several steatosis algorithms have been reported in the literature since 2005, and their diagnostic performances for prediction of hepatic steatosis are reported in Table 16-10.
The current study aimed to compare the different existing steatosis algorithms as well as an algorithm that was made according to biochemical and anthropometric data of participants from our study group.
Material and methods
The cross-sectional pilot study was conducted from July 2016 until July 2017. We evaluated 77 subjects of both sexes (42 female and 35 male). Inclusion criterion for the participants was abdominal obesity [waist circumference (WC) ≥102 cm in males and ≥88 cm in females according to the National Cholesterol Educational Program Adult Treatment Panel III (NCEP ATP III) criteria for abdominal obesity]11. Participants were both genders, between 18 and 55 years old, with females in their reproductive period. Patients with cardiovascular diseases other than hypertension, renal diseases, infective, hepatic, malignant, and autoimmune diseases were excluded from this study. We also excluded smokers, subjects in acute stress situations or with infections [high sensitive C-reactive protein (hsCRP) >10 mg/l], as well as those with recent (i.e. within three months) weight and physical activity changes. Subjects which used drugs that could affect biomarker levels of endothelial dysfunction, lipid metabolism, glucose metabolism, and menstrual cycle were also excluded. Patients that consumed alcohol more than 20g/day were excluded3. Participants with positive surface antigen of the hepatitis B (HB) virus (HBsAg), HB antibodies (anti-HBs) or hepatitis C antibodies (anti-HCV) were also excluded.
Fatty liver was assessed as presence or absence of hepatic steatosis and the grading of hepatic steatosis obtained by ultrasound, identified by a professional operator using the National Health and Nutrition Examination Survey (NHANES) III: Hepatic Steatosis, Ultrasound Images Assessment, Procedures Manual, November 201012. According to that protocol, participants with hepatic steatosis were graded into three degrees of steatosis: mild, moderate, and severe. The study was conducted according to the Declaration of Helsinki and approved by the Ethical Committee of the Clinical Centre of Vojvodina (00-55/336-21/05/2012). Written consent was obtained from every participant.
All participants attended the Department for Nutrition and Food Safety, Center of Hygiene and Human Ecology Institute of Public Health of Vojvodina for anthropometric measurements and clinical examination. Data was collected for each patient regarding the average daily alcohol consumption via a seven-day alcohol consumption survey. Twenty four hours before blood sampling, participants were asked to refrain from strenuous physical activity and consumption of alcoholic and caffeinated beverages. Venous blood was drawn from the antecubital vein after a 12h overnight fast. Analyses were performed immediately after sampling. Ultrasound examination was done the same day using GE healthcare LOGIQ7 (GE Yokogawa Medical Systems Ltd, Tokyo, Japan) and convex probe of 4 Hz.
Anthropometric and clinical measurements
Body height (BH) in centimeters was determined using Martins anthropometer with 0.1 cm precision. Body weight (BW) in kilograms was counted using a decimal scale with movable weights and precision 0.1 kg. Body mass index (BMI) was calculated as BW/BH2 (kg/m2). WC was measured using a non-elastic tape band with 0.1 cm precision. WC was measured in standing position in the horizontal plane midway between the lowest rib and the iliac crest. Hip circumference (HC) was measured using a non-elastic tape band with 0.1 cm precision in standing position around the widest portion of the buttocks, with the tape parallel to the floor13. The blood pressure was measured according to standard protocol (Riva Rocci method, mercury sphygmomanometer), in sitting position after ten minutes of rest. For the calculation of the average daily alcohol consumption, >20 g/day of ethanol, the following formula was used:
D (g) = F x Vol (ml) x 0.8 g/ml
where D is grams of ethanol, F is the fraction of ethanol (%v/v), Vol is the volume of drink in milliliters, and 0.8 g/ml is the specific gravity if ethanol.
Standard biochemical methods were used to assess blood glucose levels, alanine aminotransferase (ALT), aspartate aminotransferase (AST), gamma-glutamyl transpeptidase (GGT), uric acid, immunoturbidimetry for determination of hsCRP, and ferritin on ADVIA® 1800 Chemistry (Siemens Healthcare Diagnostics, Marburg, Germany). Parameters of lipid status: total cholesterol (TC), triglyceride (TG), and high-density cholesterol (HDL-C) were determined using standard biochemical methods while apolipoprotein A-I (apo A-I), apolipoprotein B (apo B), and Lp(a)-lipoprotein [Lp(a)] using an immunoturbidimetry method (ARCHITECT ci4100, Abbott Laboratories, USA). Low-density lipoprotein cholesterol (LDL-C) was calculated using the Friedewald equation, and in cases TG ≥4.5 mmol/l was determined using the direct method on ARCHITECT ci4100. Coagulametric Clauss method was used to assess fibrinogen on BCxp (Siemens Healthcare Diagnostics, Marburg, Germany). Insulin was measured using an electrochemiluminescent immunoassay (ECLIA) on Elecsys 2010 (COBAS – Roche Diagnostics, USA). HBsAg and anti-HCV antibodies were obtained using enzyme-linked fluorescence assay (EFLA) using mini Vidas (bioMerieux, Marcy-l’Etoile, France). Homeostatic model assessment for insulin resistance (HOMA-IR) was calculated as fasting insulin (μIU/ml) x fasting glucose (mmol/l) / 22.5 and homeostatic model assessment for beta cell function (HOMA-%B) as 20 × fasting insulin (μIU/ml)/fasting glucose (mmol/ml) – 3.514,15.
Fatty liver index (FLI) was calculated as: [(e0.953 × loge (TG) + 0.139 × BMI + 0.718 × loge (GGT) + 0.053 × (WC) – 15.745) / (1 + e0.953 × loge (TG) + 0.139 × BMI + 0.718 × loge (GGT) + 0.053 × (WC) – 15.745)] × 1006. Hepatic steatosis index (HSI) was calculated as 8 × ALT / AST + BMI (+2 if T2DM was present; +2 if female sex)7. Lipid accumulation product (LAP) was calculated for males as: (WC [cm] – 65) × (TG [mmol/L]), and for females as: (WC [cm] – 58) × (TG [mmol/L])8. Triglyceride and glucose index (TyG) was calculated as: Ln (TG (mg/dl) x glucose (mg/dl)/2)9 and visceral adiposity index (VAI) for males as: (WC/ 39.68 + (1.88 × BMI) x (TG/1.03) x (1.31/HDL-C), and for females as: (WC/ 36.58 + (1.89 × BMI) x (TG/0.81) x (1.52/HDL-C)16.
Summary statistics are presented as mean and standard deviation for metric variables and as frequencies for categorical variables. Kolmogorov-Smirnov normality test was used to test if a variable follows a normal distribution. Independent T-test and Mann-Whitney U-test were used to compare differences between two independent groups, depending on the distribution of the dependent variables. Spearman correlation coefficient was used to test the association between two ranked variables. Factor analysis was used to extract four-factor structure from biochemical and anthropometrical parameters obtained in the sample. Discovering factors that have a significant influence on the presence of NAFLD is assessed through logistic regression modeling. A total of six logistic regression models were built. In the first model, the independent variable was chosen according to the selected extracted factor. In all other models, independent variables were available diagnostic algorithms for hepatic steatosis. The regression analyses results are presented as coefficients (B), standard errors (S.E.), tests for significance (Wald, df, p), and odds ratios. The odds ratio above 1 indicates a positive association between independent and dependent variables, i.e. increase of the likelihood of dependent variable. The ratio below 1 standard, describes negative associations, i.e., that predictor decreases the likelihood of dependent variable. A significant level of 5 % is considered in each test (IBM SPSS Statisticsfor Windows, version 20.0, IBM Corp., Armonk, NY, USA).
In this study, we included 77 participants (35 male and 42 female), 44 with NAFLD (23 male and 21 female), and 33 (12 male and 21 female) without NAFLD. Anthropometric and biochemical measurements are shown in Table 2. Subjects with NAFLD had significantly higher BMI, WC, HC, and systolic (SBP) and diastolic blood pressure (DBP). Levels of AST, ALT, uric acid, and hsCRP were significantly higher in participants with NAFLD. There were significantly higher mean values of algorithms for hepatic steatosis FLI, HSI, and LAP in participants with NAFLD.
Correlation analysis between algorithms for hepatic steatosis and degree of NAFLD and HOMA-IR is shown in Table 3. There is a positive statistically significant correlation between the degree of NAFLD determined by ultrasound and FLI, TyG, LAP, and HSI. There were no significant correlations between the degree of NAFLD and VAI. Finally, given the close pathogenic link between steatosis and insulin resistance, we investigated if steatosis scores independently correlated with HOMA-IR. There is a statistically positive correlation between HOMA-IR and all the algorithms of hepatic steatosis except HSI.
According to biochemical and anthropometrical parameters determined in our subjects using factor analysis, four-factor structure was extracted explaining 56.713 % total variance. All factor loadings show values above 0.3 on the main factor; the KMO values of 0.6 indicate distinct and reliable factors. The Barlett test was significant (χ2 =1145.382, df =210, p <0.001) indicating that correlations between items are significantly different from zero. Table 4 shows pattern matrix and indicates the following factor structure: A – obesity (WC, TC, BMI, DBP, SBP, Fbg, hsCRP, glucose, uric acid); B – atherogenic factor (LDL-C, TC, APOB/A-I, Ne/Ly); C – liver damage (ALT, GGT, ferritin, AST), and D – insulin resistance (TG, HDL-C, HOMA-%B, HOMA-IR).
In order to determine the best predictor of the presence of NAFLD, six logistic regression models have been developed. Each model has a different predictor variable. The first model predicts the presence of NAFLD using variable derived from FACTOR A as the predictor variable. The other models use hepatic steatosis algorithms: FLI, HSI, TyG, LAP, and VAI as independent variables, respectively. The results are given in Table 5. The estimated coefficient for variable derived from FACTOR A, in the first model, is 1.455 and the variable has a statistically significant influence on the presence of NAFLD. Increasing the value of FACTOR A increases the chance for presence of NAFLD for more than four times. This model has an overall classification rate of 70.1 %. The second model predicts the presence of NAFLD using FLI as an independent variable with an overall classification rate of 66.2 % and the estimated regression coefficient of 0.031. The results show that increasing the values of FLI will significantly increase the probability of having NAFLD. The third model uses HSI as an independent variable for prediction of NAFLD presence with an overall classification rate of 55.8 % and a regression coefficient of 0.107. In all other built models, the predictors did not show statistically significant influence on the prediction of NAFLD presence.
Considering the hepatic and metabolic consequences of fat accumulation in the liver1-3, there is an essential medical need for a simple, accurate and cost-effective biomarker of liver fat content. There have been several proposed diagnostic scores as surrogate markers to diagnose hepatic steatosis: FLI, HSI, LAP, VAI, and TyG index6-10.
We evaluated five available algorithms for estimation of hepatic steatosis (FLI, HSI, LAP, VAI, and TyG). According to biochemical and anthropometrical parameters determined in our subject group using factor analysis, four factors have been extracted, only algorithm derived from FACTOR A – obesity (WC, HC, BMI, DBP, SBP, Fbg, hsCRP, glucose, and uric acid) is significantly positively associated with NAFLD with the best overall classification rate 70.1 %. In our study group, all parameters incorporated in factor A were significantly higher in subjects with NAFLD. Uric acid was especially interesting factor because it was significantly higher in participants with NAFLD and grouped itself in FACTOR A. The association of high uric acid levels and NAFLD is still not clear, but there are several mechanisms that may contribute to NAFLD pathogenesis mainly through interaction with insulin resistance, induction of intracellular and mitochondrial stress and activation of the NLRP3 inflammasome17. Oxidative stress inhibits enzyme of Krebs cycle aconitase and leads to accumulation of citrate and stimulation of ATP citrate lyase resulting in increased fat synthesis and inhibition of enoyl CoA hydratase which leads to impaired beta-oxidation of fatty acids and inhibition of AMPK-activated protein kinase18.
The second best predictor of NAFLD in our study group was FLI algorithm. These results were in line with the findings of previous studies19. Four parameters were incorporated into FLI: two anthropometric (BMI and WC) and two biochemical (GGT and TG)17 that are part of the metabolic syndrome phenotype and its prevalence increases in obese individuals20. In some studies, GGT was considered as an independent predictor of NAFLD considering that its prevalence increases in NAFLD. On the other hand among the markers of dyslipidemia TG is strongly associated with NAFLD21.
In our study, there were no significant differences between levels of GGT and TG in participants with and without NAFLD. This may explain why FLI did not have better discriminatory power in comparison with our algorithm. HSI was the third best predictor of NAFLD in our study group. Overall classification rate was moderate 55.8 %. HSI had relatively lower predictive power according to other studies, too7-10.
Other models (VAI, LAP, TyG) did not have statistically predictive power but they all significantly correlated with insulin resistance as assessed by HOMA-IR. Adipose tissue dysfunction (increased lipolysis, imbalanced adipocytokine production) which is indirectly expressed by VAI could directly participate in both liver steatosis and induction of inflammation and an inverse correlation with insulin sensitivity. It was anticipated that VAI would be associated with liver fat and that demonstrated in histologically defined steatosis in genotype 1 infected patients with HCV22 but not in patients with NAFLD23. TyG was effective in identifying the risk of insulin resistance assessed by HOMA-IR in Chinese population24 and glucose clamp studies in Mexican population25. It is proposed that TyG is not a good measure of peripheral insulin resistance but rather hepatic insulin resistance since it correlates well with hepatic fat content26. In our study, there was a significant correlation with TyG index and degree of NAFLD determined by ultrasound as well as in Chinese study9. Only three parameters were incorporated in LAP: gender, WC, and TG. WC is a surrogate measure of visceral fat which is thought to play a key role in insulin resistance and lipotoxicity. TG are commonly elevated in the presence of insulin resistance and strongly associated with hepatic triacylglycerol content. LAP may be valuable in recognizing patients with insulin resistance along with ectopic lipid deposition as well as liver fat accumulation. In our study, there was a significant correlation with LAP index and degree of NAFLD determined by ultrasound as well as in Italian study8.
This study has several limitations. It is a cross-sectional pilot study with a small number of patients. Presence and degree of NAFLD were determined using ultrasound which has limited sensitivity, specificity, and is operator-dependent.
In conclusion, the elevation of all anthropometrical and biochemical parameters which are incorporated in the first model for prediction of hepatic steatosis (WC, HC, BMI, DBP, SBP, Fbg, hsCRP, glucose, and uric acid) showed the best discriminatory power in the prediction of hepatic steatosis determined by ultrasound.
From the algorithms derived, the results of international studies, FLI and HSI had the best discriminatory power in predicting hepatic steatosis. Other algorithms did not have statistically significant discriminatory power but correlated with the degree of insulin resistance assessed by HOMA-IR which is thought to be a predisposition for the development of NAFLD. Further trials with a larger number of participants are needed to construct an adequate algorithm for prediction of hepatic steatosis.
Conflict of interest
The authors state that they have no financial interests and conflict of interests to disclose regarding any of the products or methods used in this study. This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.
1. Williams CD, Stengel J, Asike MI, Torres DM, Shaw J, Contreras M, et al. Prevalence of nonalcoholic fatty liver disease and nonalcohol-ic steatohepatitis among a largely middle-aged population utilizing ultrasound and liver biopsy: a prospective study. Gastroenterology. 2011; 140: 124-131.
2. Stren C, Castera L. Non-invasive diagnosis of hepatic steatosis. Hepatol Int. 2017; 11: 70-78.
3. Chalasani N, Younossi Z, Lavine JE, Diehl AM, Brunt EM, Cusi K, et al. The diagnosis and management of non-alcoholic fatty liver disease: Practice guideline by the American Association for the Study of Liver Diseases, American College of Gastroenterology, and the American Gastroenterological Association. Am J Gastroenterol. 2012; 107: 811-826.
4. Targher G, Day CP, Bonora E. Risk of cardiovascular disease in patients with nonalcoholic fatty liver disease. New Engl J Med. 2010; 363: 1341-1350.
5. Szczepaniak LS, Nurenberg P, Leonard D, Browning JD, Reingold JS, Grundy S, et al. Magnetic resonance spectroscopy to measure hepatic triglyceride content: prevalence of hepatic steatosis in the general population. Am J Physiol Endocrinol Metab. 2005; 288: E462-E468.
6. Bedogni G, Bellentani S, Miglioli L, Masutti F, Passalacqua M, Castiglione A, et al. The Fatty Liver Index: a simple and accurate predic-tor of hepatic steatosis in the general population. BMC Gastroenterol. 2006; 6: 33.
7. Lee JH, Kim D, Kim HJ, Lee CH, Yang JI, Kim W, et al. Hepatic steatosis index: a simple screening tool reflecting nonalcoholic fatty liver disease. Dig Liver Dis. 2010; 42: 503-508.
8. Bedogni G, Kahn HS, Bellentani S, Tiribelli C. A simple index of lipid overaccumulation is a good marker of liver steatosis. BMC Gas-troenterol. 2010; 10: 98.
9. Zhang S, Du T, Zhang J, Lu H, Lin X, Xie J, et al. The triglyceride and glucose index (TyG) is an effective biomarker to identify nonal-coholic fatty liver disease. Lipids Health Dis. 2017; 16: 15.
10. Fedchuk L, Nascimbeni F, Pais R, Charlotte F, Housset C, Ratziu V. Performance and limitations of steatosis biomarkers in patients with nonalcoholic fatty liver disease. Aliment Pharmacol Ther. 2014; 40: 1209-1222.
11. Expert Panel on Detection, Evaluation, and Treatment of High Blood Cholesterol in Adults. Executive Summary of The Third Report of The National Cholesterol Education Program (NCEP) Expert Panel on Detection, Evaluation, And Treatment of High Blood Cholesterol In Adults (Adult Treatment Panel III). JAMA. 2001; 285: 2486-2497.
12. National Health and Nutrition Examination Survey (NHANES) III. Hepatic Steatosis Ultrasound Images Assessment Procedures Manual. 2010. Available at: http://www.cdc.gov/nchs/data/nhanes/nhanes3/Hepatic_Steatosis_Ultrasound_Procedures_Manual.pdf, last accessed: 28/07/2017.
13. Waist Circumference and Waist-Hip Ratio. Report of a WHO Expert Consultation. 2008. Available at: https://apps.who.int/iris/bitstream/handle/10665/44583/9789241501491_eng.pdf?sequence=1, last accessed: 28/07/2017.
14. Matthews DR, Hosker JP, Rudenski AS, Naylor BA, Treacher DF, Turner RC. Homeostasis model assessment: insulin resistance and beta-cell function from fasting plasma glucose and insulin concentrations in man. Diabetologia. 1985; 28: 412-419.
15. Garg MK, Dutta MK, Mahalle N. Study of beta-cell function (by HOMA model) in metabolic syndrome. Indian J Endocrinol Metab. 2011; 15: S44-S49.
16. Amato MC, Pizzolanti G, Torregrossa V, Misiano G, Milano S, Giordani C. Visceral adiposity index (VAI) is predictive of an altered adipokine profile in patients with type 2 diabetes. PLoS One. 2014; 9: e91969.
17. Lombardi R, Pisano G, Fargion S. Role of Serum Uric Acid and Ferritin in the Development and Progression of NAFLD. Int J Mol Sci. 2016; 17: 548.
18. Kanbay M, Jensen T, Solak Y, Le M, Roncal-Jimenez C, Rivard C, et al. Uric acid in metabolic syndrome: From an innocent bystander to a central player. Eur J Intern Med. 2016; 29: 3-8.
19. Motamed N, Sohrabi M, Ajdarkosh H, Hemmasi G, Maadi M, Sayeedian FS, et al. Fatty liver index vs waist circumference for predicting non-alcoholic fatty liver disease. World J Gastroenterol. 2016; 22: 3023-3030.
20. Milić S, Lulić D, Štimac D. Non-alcoholic fatty liver disease and obesity: biochemical, metabolic and clinical presentations. World J Gastroenterol. 2014; 20: 9330-9337.
21. Bedogni G, Miglioli L, Masutti F, Tiribellli C, Marchesini G, Bellentani S. Prevalence of and risk factors for nonalcoholic fatty liver disease: the Dionysos nutrition and liver study. Hepatology. 2005; 42: 44-52.
22. Petta S, Di Marco V, Di Stefano R, Cabibi D, Cammà C, Marchesini G, et al. TyG index, HOMA score and viral load in patients with chronic hepatitis C due to genotype 1. J Viral Hepat. 2011; 18: e372-e380.
23. Musso G, Cassader M, Gambino R. Diagnostic accuracy of adipose insulin resistance index and visceral adiposity index for progressive liver histology and cardiovascular risk in nonalcoholic fatty liver disease. Hepatology. 2012; 56: 788-789.
24. Du T, Yuan G, Zhang M, Zhou X, Sun X, Yu X. Clinical usefulness of lipid ratios, visceral adiposity indicators, and the triglycerides and glucose index as risk markers of insulin resistance. Cardiovasc Diabetol. 2014; 13: 146.
25. Guerrero-Romero F, Simental-Mendía LE, González-Ortiz M, Martínez-Abundis E, Ramos-Zavala MG, Hernández-González SO, et al. The product of triglycerides and glucose, a simple measure of insulin sensitivity. Comparison with the euglycemic hyperinsulinemic clamp. J Clin Endocrinol Metab. 2010; 95: 3347-3351.
26. Simental-Mendía LE, Simental-Mendía E, Hernández-Rodrígez H, Rodrígez-Moran M, Guerrero-Romero F. The product of triglycerides and glucose as biomarker for screening simple steatosis and NASH in asymptomatic women. Ann Hepatol. 2016; 15: 715-720.