Path analysis model to identify the effect of poor diet quality on NAFLD among Iranian adults from Amol Cohort Study

Nonalcoholic fatty liver disease (NAFLD) is expanding as a global health problem with approximately 25% of the world's population affected by it. Dietary modification is one of the most important strategies for preventing NAFLD. The association between nutrient density and the Healthy Eating Index 2015 (HEI2015) with NAFLD demonstrates that nutrient density is an independent predictor of NAFLD in Iranian adults [fully adjusted model: OR (95% CI)tertile3vs.1: 0.68 (0.54–0.85), P for trend = 0.001]. However, a favorable association between NAFDL and diet quality (HEI 2015) is more pronounced in participants with abdominal obesity [fully adjusted model: OR (95% CI)tertile3vs.1: 0.63 (0.41–0.98), P for trend = 0.03]. Based on the gender-stratified path analysis, diet quality indirectly through Waist-to-Height Ratio (WHtR), C-reactive protein (CRP), and metabolic syndrome in women, and men through WHtR, hemoglobin A1c (HBA1c), CRP, and metabolic syndrome affects NAFLD. Nutrient density directly and indirectly in women through WHtR, CRP, and metabolic syndrome, and in men indirectly through WHtR, hemoglobin A1c, and metabolic syndrome negatively affect NAFLD. Hence, in these subjects; we can provide early dietary intervention and education to prevent progression to NAFLD.


Baseline characteristics of the study participants
Of 2956 (47.09 ± 14.46) recruited in the cohort, 1332 (45.1%) were women.The prevalence of NAFLD was 46% (n = 1360).Details of socio-demographic, anthropometric, biochemical parameters, physical activity and dietary assessment according to gender and disease status are presented in Tables 1 and 2. NAFLD participants had a higher BMI, WC, and WHtR compared to healthy adults.Compared to those without NAFLD, diabetes and metabolic syndrome comorbidities were much more highly prevalent in patients with NAFLD (all P < 0.001).In women with NAFLD, the medication history approached significance.This was significant for men who used glucose-lowering agents (p = 0.001).In both genders, there were high levels of TG, total cholesterol, and LDL in patients with NAFLD, whereas HDL levels were significantly lower in patients with NAFLD compared to patients without NAFLD.There was no difference between the groups in other variables.
The scores of nutrient adequacy (p < 0.001) and healthy eating index (p = 0.01) were significantly lower in women with NAFLD than in healthy women.Despite the low score of the healthy eating index and nutritional adequacy in men with NAFLD compared with healthy subjects, these differences were not significant (p > 0.05).
Additional baseline characteristics and dietary assessments across tertile categories of each dietary index (HEI 2015 and NRF 9.3 ) stratified based on gender and NAFLD status are included in Supplementary Tables 1 & 2. A comparison of anthropometric and biochemical parameters across categories of HEI 2015 showed that in female participants, a higher HEI 2015 score was associated with lower rates of abdominal obesity (p < 0.001), lower HBA 1c (p = 0.04), and a lower intake of dietary energy density (DED) (p < 0.001).
Higher HEI 2015 and NRF 9.3 scores were associated with lower diastolic blood pressure among men compared to women (p = 0.02 and p = 0.01, respectively).A less energy-dense diet was also associated with higher HEI 2015 and NRF 9.3 scores (p < 0.001) in both genders.

Association among healthy eating index and nutrient density with NAFLD
The multiple-adjusted odds ratio (95% confidence interval) in Tables 3 and 4 demonstrated an inverse association of nutrient density with NAFLD for the highest (vs.lowest) tertile of NRF 9.3 [0.68 (0.54-0.85),P for trend = 0.001].When stratified by gender and abdominal obesity, greater nutrient adequacy was associated with a lower risk of NAFLD in participants with abdominal obesity [0.62 (0.40-0.95), P for trend = 0.03] compared to those without abdominal obesity [0.69 (0.52-0.90),P for trend = 0.007].
The results for nutrient quality and lower odds of NAFLD prevalence in both genders were also similar [in men: 0.68 (0.50-0.93),P for trend = 0.01; in women, 0.64 (0.46-0.89),P for trend = 0.01].The findings of the full multiple-adjusted model of HEI 2015 , stratified by gender and abdominal obesity, showed that the favorable association was only presented on subjects with abdominal obesity [0.63, (0.41-0.98),P for trend = 0.03].www.nature.com/scientificreports/
The indirect relationship between HEI 2015 and NAFLD was explained through WHtR, and metabolic syndrome, as well as WHtR, CRP, and metabolic syndrome in women.In men, these effects were partially mediated through WHtR, and metabolic syndrome, as well as WHtR, HBA 1c, and metabolic syndrome.
The standardized path coefficients (β), standardized total effects, along with direct and indirect effects concerning HEI 2015 and NRF 9.3 among the participants of the study were documented in Supplementary Table 3.

Association between nutrient density and diet quality with NAFLD
In the current cross-sectional study, nutrient density (high intake of vitamin A, vitamin C, calcium, potassium, and iron, and the low consumption of sodium and saturated fatty acid) was an independent predictor of reduced prevalence for NAFLD in adult Iranians.The association between diet quality and NAFDL was pronounced in participants with abdominal obesity, where better compliance with the healthy eating index was markedly linked to a lower risk of NAFLD.
Although previous studies have explored the role of individual nutrients in the development, progression, and treatment of NAFLD, nutrient pattern as a total has not been investigated 27,28 .
Total NRF scores (nutrient adequacy) significantly reduced the NAFLD risk by almost 50% in both genders, although there were differences for individual items between patients with NAFLD vs. those without NAFLD.In the current study, low overall intake of nutrients, including vitamin A, vitamin C, potassium, calcium, and iron, and the high consumption of sodium and saturated fatty acid in NAFLD participants, suggest lower the micronutrient density, particularly in women.This pattern was similar only to Iron and SFA among men (Supplementary Table 2).

Table 2.
Components and total scores of HEI-2015 and NRF9.3 according to sex with and without NAFLD among adults participants of Amol Cohort Study (n 2956), Iran, 2016-2017 1 .NAFLD non-alcoholic fatty liver disease, HEI Healthy Eating Index, NRF9.3 Nutrient-Rich Food Index 9.3.Significant at P < .05for Independent t-test for continuous variables. 1Means ± SDs (all such values). 2A higher score indicates a higher diet quality, except for added sugars, saturated fats, sodium components., and refined grains. 3Maximum score for each component of HEI indicated in parentheses. 4Calculated as the sum of all component scores. 5Calculated as the sum of scores for nine nutrients to encourage (i.e., protein, dietary fiber, vitamins A, C, and E, calcium, iron, potassium, and magnesium) minus the sum of scores for three nutrients to limit (i.e., added sugar, saturated fats, and sodium).Adequacy components (score) www.nature.com/scientificreports/These findings have been confirmed in other studies 29 .Panera et al. argued that the composition of the macro and micronutrients is more critical in the etiology and management of NAFLD than total calorie intake 30 .Imamura et al., in a systematic assessment of males and females in 187 countries, concluded that females had better dietary patterns than males 31 .However, there is a dearth of evidence indicating gender differences in the association between diet quality and adequacy and NAFLD risk, and the results remain controversial.According to the HEI components, whole fruit intake was significantly lower in the women with NAFLD than the healthy women, which may suggest that fruit intake reduced NAFLD risk by functioning as a source of vitamin C, vitamin A, calcium, potassium, and dietary fiber 32 .NAFLD is associated with oxidative stress and low-grade inflammation 33 .Micronutrient adequacy protects hepatic cells from lipotoxicity-induced oxidative stress.This status can promote trigger inflammation known to contribute to metabolic dysfunction and disturbing vitamin E trafficking through the gut axis 34 .A higher intake of vitamin C, vitamin A, calcium, and potassium  is a major preventive nutrient for metabolic syndrome and NAFLD 35 .Van Tien et al., in a Multi-institutional Collaborative Cohort of the Japanese population using 1588 subjects, proposed that a nutrient diet rich in vitamins, fiber, iron, and potassium was linked to a lower prevalence of NAFLD 36 .A cohort study assessing dietary intake in NAFLD patients found that recommended intakes of calcium, vitamin A, iron, vitamin B1, vitamin B2, zinc, and magnesium were not met 37 .Aktary et al. 38 , in a cross-sectional investigation of the dietary intake and health profile of a sample of Canadian adults with NAFLD (n = 42), demonstrated that NAFLD patients had poor micronutrients such as magnesium, calcium, vitamin D, and vitamin E. The increased consumption of saturated fatty acid has been shown to cause mitochondrial dysfunction, increased oxidative stress, and lowgrade inflammation 39 .Moreover, a recent meta-analysis documented a positive association between high sodium consumption and a 60% greater risk of NAFLD 40 .
In both genders, NAFLD patients significantly had a lower intake of iron and higher consumption of saturated fatty acids than healthy participants did.The association between select micronutrients and NAFDL warrants further investigation.
The association between diet quality and NAFLD aligned with those in a US population-based study (n = 10,858, mean age = 42.9 years, 47.1% men) 40 .Women (healthy vs. NAFDL participants) had more consumption of seafood, plant proteins, and whole fruits, while healthy men consumed better fatty acid composition (greater ratio of polyunsaturated (PUFAs) and monounsaturated fatty acids (MUFAs) to SFA) relative to NAFLD participants.
The reason the relationship was significant only among those with abdominal obesity might be that these participants were more prone to have NAFLD 41 , and weight loss as primary therapy for most NAFLD patients has been documented 42 .Improving diet quality, for example, incorporating more vegetables, whole fruits, whole grains, and limiting the content of the SFA sources, may have important benefits in preventing weight gain or promoting weight loss in adults of both genders and even being predisposed to obesity 43,44 .Thus, implementing healthy dietary patterns in the clinical setting may be a viable alternative to encourage weight loss and reduce the risk of developing NAFLD; however, more research is required to identify short-term and long-term beneficial effects of these dietary interventions in the context of NAFLD.
Considering these findings, primary prevention of NAFLD by improving micronutrient adequacy may be a more effective and beneficial objective for the patient at risk of developing NAFLD.

Influencing factors of NAFLD
Based on the gender-stratified path analysis, in women, diet quality indirectly through WHtR, CRP, and metabolic syndrome, and in men through WHtR, Hemoglobin A 1c , and metabolic syndrome affected NAFLD.
While nutrient density directly and indirectly (through HBA 1c and metabolic syndrome) reduced NAFLD risk in women, it only had a direct protective effect on NAFLD in men.
The current findings apply to understanding how habitual diet shapes anthropometric indices, metabolic risk profiles, and health outcomes in the Iranian surveyed population.
The findings indicate age is a significant predictor of HEI, physical activity, cardiometabolic parameters (e.g., WHtR, HBA1c), and NAFLD risk.Previous studies revealed differences in lifestyle factors (e.g., eating behaviors and physical activity) based on age and gender 45 .Age directly influences diet quality, and is a significant predictor of physical activity 46 .Research findings affirmed that increased physical activity was correlated with better diet quality measured by HEI 2015 47,48 .
NAFLD risk has been more prevalent in older people 49 .Asian studies also reported that under the age of 50 years, NAFLD was more prevalent in men, but in populations over 50 years, it was higher in women 50,51 .Processes associated with aging are considered possible contributing mechanisms in the pathogenesis of NAFLD and cardiometabolic disorders 52 .
In line with our results, growing evidence demonstrated obesity and metabolic syndrome were independently linked with NAFLD irrespective of other cardiometabolic risk factors [53][54][55] .The NAFLD's pathophysiology concerning obesity involves excess fat deposition in the liver and insulin resistance development, which are pivotal in the progression of NAFLD 56 .In a recent study, the correlation between cardiometabolic disorders and inflammation and the incidence of NAFLD was verified 57 .
Evaluating overall diet quality, rather than specific nutrients or food components, is more effective in identifying diet-disease associations 5 .HEI is an indicator of determining the nutritional balance and predicting health risks.However, there are discrepancies in the risk prediction of disease in diet quality because of unmeasured interactions with various effects of modifiers or mediators 41 .Furthermore, the relationship and interrelationship between diet quality and health risk cannot be precisely calculated using common statistical methods.Path analysis may help evaluate this connection within a conceptual framework by concurrently investigating all relevant regression pathways, including direct and indirect 7 .Applying Path analysis makes it easier to assess the mediating role of diet quality and adequacy, anthropometric and metabolic parameters, and NAFLD risk.Additionally, this method permits a thorough understanding of such a relationship and allows a more precise interpretation of results.
A higher HEI score signifies a more balanced nutritional intake, leading to improved insulin sensitivity, reduced inflammatory markers, and a decreased likelihood of metabolic syndrome 58 .The NRF index, unlike the HEI, evaluates individual foods and simultaneously provides a precise measure of overall diet nutrient density.Individuals can meet their nutrient requirements without excessive energy intake and gaining weight by choosing nutrient-dense foods 59 .Moreover, the NRF index prioritizes nutrients critical for preventing metabolic disorders and NAFLD 60,61 .
Although the rationale for WHtR, CRP, and HBA 1c as strong predictors for NAFLD has been justified by several previous studies 40 , to our knowledge there is no research on metabolic dysfunction being a mediator of NAFLD as the present study is the only Path analysis modeling study in this regard.A meta-analysis reported the superiority of centralized obesity measures, particularly, WHtR, for NAFLD risk detection 62 .Evidence has shown that visceral adiposity is the main adipose depot responsible for fatty liver and is associated with it in a dose-dependent manner 37 .In the present study, high WHtR values were associated with CRP and HBA 1c .Several epidemiologic studies have proved the causal link between obesity and increasing liver disease in individuals 63 .According to pathophysiology and clinical studies, the progression of NAFLD is caused by an imbalance between lipid intake and disposal, which leads to oxidative stress and hepatocyte injury 63 .This finding is important to broaden the discussion about the high level of early inflammatory markers in obese adults and clarify this relationship.CRP acts as a regulator of nitric oxide production in the endothelium and coordinates the production and secretion of various cytokines, increasing the pro-inflammatory activity of different adipokines.The measurement of CRP and HBA 1c were independent predictors of metabolic syndrome in other cohort studies (OR 1.22, 95% CI; 1.12 to 1.35; OR 1.57, 95% CI; 1.35 to 1. 82, respectively) 6,8 .MetS was associated with oxidative stress and chronic low-grade inflammation 63 , and NAFLD is one criterion of MetS 30 .
Additionally, our results revealed a gender inequalities association of diet quality and adequacy with anthropometric and metabolic parameters and NAFLD.The beneficial effects of nutrient adequacy of NAFLD risk in women, directly and indirectly through changes in HBA1c and metabolic syndrome, and in men directly were potentially exerted.
These results reveal the importance of gender-specific interventions to control NAFLD and also the pivotal role of diet adequacy and adequacy in obesity and glycemic control and preventing low-grade systemic inflammation and metabolic syndrome among high-risk individuals.
Key strengths of the study include the relatively large community-based study sample recruited from rural and urban areas of Amol city that afforded us sufficient power to probe small effects, incorporating multiple potential anthropometric and biochemical mediators, and assessing their mediation role simultaneously in the relationship between diet quality and adequacy and the risk of NAFLD for the first time, using a reliable and validated semi-quantitative FFQ 64 developed for the Iranian population, which results in a better representation of the participants' dietary habits.However, some potential limitations of this study need to be acknowledged.First, because of the cross-sectional design of the study, drawing any causal inference from the association would be incorrect.Second, although liver biopsy is a gold standard for diagnosing NAFLD, we used sonography for evaluating NAFLD due to the risks associated with liver biopsy and the impossibility of applying it in populationbased studies.Furthermore, the sensitivity of the ultrasound for the detection of moderate to severe fatty liver is approximately 85%, which keeps it a preferred and practical modality for diagnosing NAFLD in epidemiological settings.Third, since dietary intake and other socio-demographic parameters in Amol may differ from those in other parts of the country, our results cannot be extended to all Iranians.Fourth, other effective factors, including meal and snack patterns and cooking methods, were not investigated in the current study, so the observed associations are not entirely explained.Finally, we could not completely rule out residual confounding due to unknown or unmeasured confounders in this study.

Conclusion
Nutrient density was an independent predictor of NAFLD prevalence in Iranian adults.The association between diet quality (assessed by the HEI 2015 ) and NAFDL was more pronounced in participants with abdominal obesity.The beneficial indirect effects of diet quality and nutrient density on NAFLD prevention were mediated by changing WHtR, HBA1c, CRP, and metabolic syndrome.Therefore, for subjects with MetS, high WHtR, high HBA1c, and CRP, we can provide early dietary intervention and proper education to prevent progression to NAFLD.Future research assessing the longitudinal relationship using prospective study designs is needed to better understand these relationships and confirm the findings in the present study.

Study design, setting, and participants
This cross-sectional study was conducted within the framework of the Amol Cohort Study (AmolCS), a prospective study conducted on rural and urban residents of Amol City in the North of Iran, which evaluated obesity-related metabolic disorders and CVD.The AmolCS was set up in two phases.In the first phase started in 2009, 7104 participants aged 10-90 years through sixteen strata with ten-year intervals (10-19, 20-29, 30-39,  40-49, 50-59, 60-69, 70-79, and 80-89 years) were randomly selected across rural and urban health centers of Amol city.The second phase of the AmolCS, including 5147 adult participants ≥ 18 years of age, was launched in 2017, and the data from the second phase of the cohort study were used in the present analysis.
The exclusion criteria for the participants were pregnancy/lactation, following a specific dietary or physical activity regimen, history of disease including Wilson's disease, autoimmune liver disease, hemochromatosis, virus infection, alcoholic fatty liver, malignancy, thyroid disorder, and autoimmune diseases, as well as participants with significant alcohol consumption (> 30 g/d for men and > 20 g/d for women).Written informed consent was obtained from all participants before the study.Further details on the project are available in the previous studies 65,66 .
The study design and selection flowchart are outlined in Fig. 2. Approval for this study was received from the Iran University of Medical Sciences (IUMS) ethics committee (NO: IR.IUMS.REC.1399.1393).www.nature.com/scientificreports/

Data collection
Written and verbal informed consent was obtained from participants.The documentation of participants' information in the second phase of the cohort study included demographic and lifestyle characteristics, clinical testing results, dietary assessment, and NAFLD diagnosis.

Dietary assessment
A validated semi-quantitative food-frequency questionnaire (FFQ) was used to evaluate the habitual intake of 168 food items 64 .For each food item on the list, participants were asked about the usual frequency of consumption in a commonly used unit or portion size (daily, weekly, and monthly) over the previous year.The consumption intake of each food item was calculated as grams/day by household measures 67 .Nutrient and energy intake was calculated using the food composition table (FCT) of the United States Department of Agriculture (USDA) 68 and the Iranian FCT for traditional Iranian foodstuffs 69 .

Healthy eating index and nutrient density
In order to evaluate the quality of the diet, the Healthy Eating Index 2015 (HEI 2015 ) was calculated using the method explained by the National Cancer Institute and the US Department of Agriculture (USDA) center 70,71 .In this index, nutritional intakes were compared with the US dietary guidelines.Scores could range between 0 and 100, with a higher score suggesting a healthier diet.The Nutrient Rich Food 9.3 score (NRF 9.3 ) was calculated for the whole diet to measure nutrient density.Drewnowski et al. 72 described the details of NRF  72 were used in the present study.

Diagnosis and assessment of NAFLD
All the study subjects underwent ultrasonography of the abdomen to assess the hepatic parenchyma and biliary tree, performed by a single expert radiologist blinded to the clinical, and laboratory data of the participants using an ultrasound system (Esaote SpA, Genova, Italy) with transducer (frequency bandwidth 3-5 MHZ).

Laboratory testing
After an overnight fast of at least 12 h, intravenous blood samples of each participant were collected, with one tube for Ethylene Diamine Tetraacetic Acid (EDTA) anticoagulation and one tube for separation gel coagulation, and then centrifuged at 3000 rpm for 10 min at 4 °C; the aliquots were stored at -80 °C until use.Fasting blood sugar (FBS) was measured by the hexokinase method, and lipid profile containing total cholesterol (TC), highdensity lipoprotein cholesterol (HDLc), and low-density lipoprotein cholesterol (LDLc) were measured by the enzymatic method using an Auto-analyzer BS200 (Mindray, Shenzhen, China) and diagnostic kits (Pars Azmoon Co., Tehran, Iran).Alanine transaminase (ALT), aspartate transaminase (AST), γ-glutamyl transaminase (GGT), and CRP were measured with a rating method.Hepatitis B surface antibody, hepatitis B surface antigen, hepatitis C virus antibody, and hepatitis B core antibody were assessed by Enzyme-linked immunosorbent assay (ELIZA) kits (Pishtaz Teb Co., Tehran, Iran).Ten percent of the blood samples were re-evaluated by the Iranian National Reference Laboratory.The coefficients of variations ranged from 1.7% to 3.8% for all laboratory values.

Anthropometric variables and covariates
For each participant, covariates of demographic and lifestyle characteristics, smoking and alcohol drinking status, and physical activity were collected with the questionnaire-based interview.Data on physical activity was completed using the validated international physical activity questionnaire (IPAQ), as metabolic equivalent minutes per minute per week (MET-min/week) 73 .
Trained assistants employed medical equipment to measure height (cm), weight (kg), waist circumference (WC, cm), systolic blood pressure (SBP), and diastolic blood pressure (DBP).Height and weight were measured with the subjects wearing light clothing and no shoes.Height was recorded at the nearest 0.1 cm and weight was to the nearest 0.1 kg.Body mass index (BMI) was computed as weight (kg) divided by height squared (m 2 ).The girth of the midpoint between the lowest point of the rib and the upper edge of the iliac crest was calculated as waist circumference (WC).The WC measurement was taken to the nearest 0.1 cm.WHtR was respectively calculated as WC divided by height 74 .
Blood pressure was measured 2 times after at least 5 min of rest using the standardized desktop sphygmomanometer.The average blood pressure derived from two measurement readings was used 75 .All variables were collected according to standard interview guidelines and standard protocols 75,76 .

Descriptive analysis
Descriptive statistics included the frequency count (percent) for categorical variables and mean and standard deviation (SD) for continuous variables.The normality of continuous variables was evaluated by the Shapiro-Wilk statistical test.Baseline characteristics and dietary intake of the participants across tertile categories of each dietary index (HEI 2015 and NRF 9. www.nature.com/scientificreports/ a Bonferroni post-hoc analysis to make multiple comparisons for continuous variables and chi-square test for categorical variables.The associations were adjusted for energy intake.

Logistic regression
To examine the association between nutrient density and healthy eating indices with NAFLD, multiple logistic regression was used in several models for all participants.The obtained findings were adjusted for confounding factors, including age, sex, WC, BMI, energy intake, physical activity, and smoking.Further adjustment for lowering serum lipid drugs, lowering hypertension (HPTN) drugs, lowering serum glucose drugs, residual areas, the presence of heart disease, and diabetes was applied in the last model.Stratified analyses by gender, as well as waist circumference status, were also conducted.We used tertile categories as an ordinal variable to assess the trend of odds ratios across increasing tertiles of dietary indxes scores.For potential confounding factors, a univariate analysis was applied, and those with a P-value for entry (Pe) lower than 0.20 were included in the final multiple models.

Path analysis model
A path analysis model was utilized to assess the hypothesized model.Path analysis is based on the maximumlikelihood estimation of the entire system of the hypothesized model and assesses the degree to which the data fits the specified model 77 .In the present analysis, we performed a two-step strategy outlined by Anderson and Gerbing 78 .The initial hypothesized model was evaluated by Path analysis to measure the fit and path coefficients.We computed the standardized regression weights, standardized total effects, as well as direct and indirect effects.
In this study, the significance level was set at 5%, and all reported P-values are based on two-sided tests and the corresponding 95% confidence interval (CI).All the statistical analyses were done via SPSS version 24 (Statistical Package for Social Science, SPSS Inc, Chicago, IL, USA) software.A path analysis model was employed using AMOS 23.0 to build a measurement model and verify the structural relationship between nutrient density and healthy eating index with NAFLD.

Figure 1 .
Figure 1.The research conceptual framework of the relationship between Dietary Indices (HEI 2015 & NRF 9.3 ) and NAFLD considering the effect of mediators and other variables.

Figure 2 .
Figure 2. Flowchart of the study design.

Figure 3 .
Figure 3.The final path analysis model for the relationship between HEI2015 and NAFLD.The numbers on the paths represent standardized regression coefficients (standardized effects).The bold-face coefficients represent values for women, and the coefficients below them represent values for men.The significance level of the comparison of each effect between men and women is depicted by asterisks (*p < 0.05, **p < 0.01, ***p < 0.001).Pink arrows refer to females whereas blue arrows refer to males.HEI Healthy Eating Index, HBA1c Hemoglobin A1c, MetS Metabolic syndrome, CRP c-reactive protein, NAFLD non-alcoholic fatty liver disease.

Table 1 .
Baseline characteristics of the study adult participants by gender with and without NAFLD (n = 2956) Iran, 2016-2017.NAFLD non-alcoholic fatty liver disease, PA physical activity, MET: metabolic equivalent of task, BMI body mass index, ALT Alanine transaminase, AST Aspartate transaminase, GGT Gammaglutamyl transferase, ALKP Alkaline phosphatase, CRP c-reactive protein, DED dietary energy density.Data are Mean ± SD (all such values) unless indicated.Significant at P < .05for independent t-test for continuous variables and chi-square test for dichotomous variables.

Table 3 .
Multiple-adjusted odds ratio and 95% confidence intervals (95% CI) for non-alcoholic fatty liver disease (NAFLD) according to tertiles of dietary index in all adult participants and stratified by sex (n 2956).HEI Healthy Eating Index; NRF9.3, Nutrient-Rich Food Index 9.3 Model 1: adjusted for age and sex (except for sex-stratified analysis).Model 2: Additional adjustment for WC, BMI, energy intake, physical activity, and smoking.Model 3: additional adjustment for lowering serum lipid drugs, lowering HPTN drugs, lowering serum glucose drugs, residual areas, heart disease, diabetes.

Table 4 .
Multiple-adjusted odds ratio and 95% confidence intervals (95% CI) for non-alcoholic fatty liver disease (NAFLD) according to tertiles of dietary index in all adult participants stratified by waist circumference (n = 2956).HEI, Healthy Eating Index; NRF9.3, Nutrient-Rich Food Index 9.3 Model 1: adjusted for age and sex.Model 2: additional adjustment for BMI, energy intake, physical activity, and smoking.Model 3: additional adjustment for lowering serum lipid drugs, lowering HPTN drugs, lowering serum glucose drugs, residual areas, heart disease, diabetes.† Abdominal obesity: waist circumference > 102 cm for men and > 88 cm for women.
9.3 calculation.In brief, the calculation of NRF9.3 is based on nine qualifying nutrients, including protein, fiber, vitamins A, C, and D, calcium, magnesium, potassium, and iron; and three disqualifying nutrients, including saturated fat, added sugar, and sodium.NRF 9.3 was calculated as the sum of the percentage of the reference daily values (RDVs) for qualifying nutrients (NR 9 ) minus the sum of the percentage of maximum recommended value (MRVs) for disqualifying nutrients (Lim 3 ).All daily values calculated per 2000 kcal and the RDVs and MRVs suggested by Drewnowski et al. al. (based on several sources, i.e., WHO and FDA)