ORIGINAL RESEARCH article

Front. Nutr., 09 September 2022

Sec. Nutrition and Metabolism

Volume 9 - 2022 | https://doi.org/10.3389/fnut.2022.962834

Non-alcoholic fatty liver disease and dairy products consumption: Results from FASA Persian cohort study

  • 1. Noncommunicable Diseases Research Center, Fasa University of Medical Sciences, Fasa, Iran

  • 2. Department of Nutrition, School of Medicine, Zanjan University of Medical Sciences, Zanjan, Iran

  • 3. National Nutrition and Food Technology Research Institute, Faculty of Nutrition Sciences and Food Technology, Shahid Beheshti University of Medical Sciences, Tehran, Iran

  • 4. Department of Public Health, School of Health, Larestan University of Medical Sciences, Larestan, Iran

  • 5. Student Research Committee, Shiraz University of Medical Sciences, Shiraz, Iran

  • 6. Health Policy Research Center, Institute of Health, Shiraz University of Medical Sciences, Shiraz, Iran

Article metrics

View details

37

Citations

56,4k

Views

1,6k

Downloads

Abstract

Background/objectives:

There are limited data on the association between dairy products consumption and nonalcoholic fatty liver disease (NAFLD). This study was conducted to evaluate the association between total intake of different dairy products and fatty liver index (FLI), a marker of subclinical fatty liver.

Methods:

A total of 7,540 adults were included in this population-based cohort study. Dairy products consumption was evaluated by a validated interview questionnaire for food intake frequency. The FLI was calculated using the standard formula. Liver enzyme levels, lipid profiles, glycemic profiles and demographic characteristics were recorded for all participants. Univariate and multiple logistic regression models were used to respectively assess the mean percentage difference of mean FLI and odds ratios (ORs) for subclinical NAFLD across quantiles of dairy consumption.

Results:

The mean age of all participants was 48.81 ± 9.631 years. FLI measurements for men and women were 26.71 ± 23.39 and 39.99 ± 26.64 respectively, which was significantly higher in women (P < 0.05). Multiple logistic regression analysis demonstrated that the amount of milk consumption was an independent preventive predictor of FLI (OR = 0.96; 95% CI: 0.94–0.99), conversely, it did not predict higher levels of liver enzymes. In term of cheese intake, participants in the third tertile of cheese intake had significantly lower FLI than lower tertiles (P = 0.01). However, there wasn't any significant association between cheese intake and the odds of FLI in the multivariate model (P > 0.05). We didn't find any significant association between yogurt consumption and NAFLD indicators (P > 0.05).

Conclusion:

Higher milk consumption was inversely associated with FLI. However, there wasn't any significant association between other types of dairy products and NAFLD indicators.

Introduction

Nonalcoholic Fatty Liver Disease (NAFLD) is a clinical histopathological condition in which triglyceride accumulates in the liver cells of patients with little or no history of alcohol consumption (1). This disease usually happens when fat accumulation reaches 5–10% of the liver's weight. NAFLD has a wide histological spectrum ranging from simple steatosis to non-alcoholic steatohepatitis (NASH) and more severe complications like cirrhosis or liver fibrosis (2, 3).

The total Prevalence of NAFLD based on a Meta-analytic study from 22 countries, has been reported as 25.24 % (4). In Asia, the number of patients with NAFLD is gradually increasing (5) due to lifestyle changes (high-fat diet, low physical activity, central obesity, and type 2 diabetes mellitus), which is comparable to Western countries (6).

NAFLD can lead to cirrhosis, cardiovascular disease, type 2 diabetes mellitus, and chronic kidney disease if left untreated (7).

Today, various studies have shown that improving the microbial profile of the intestine plays an important role in the prevention and treatment of NAFLD by reducing oxidative stress and inflammation, as well as decreasing triglycerides (TG) accumulation in the liver (812). As a result, the chance of hepatic steatosis will decrease. Dairy products are widely included in daily diet due to their high nutrition contents of protein, fat, minerals, and vitamins (13). Currently, it is recommended to limit the consumption of dairy products due to high levels of saturated fat and cholesterol, in order for reducing the risk of cardiovascular disease. However, a cross-sectional analysis of the Oslo Health Study has shown that increasing the frequency of consumption of some dairy products, such as cheese, is positively associated with improved serum concentrations of high-density lipoprotein (HDL) and is inversely related to serum triglyceride levels (14). Moreover, the consumption of some dairy products, especially those fortified with probiotics, appears to reduce the level of low-density lipoprotein (LDL) cholesterol in human intervention studies (15).

This study aims to focus on the relevance of dairy products consumption with NAFLD and fatty liver index in a large-population cohort in Fasa, Iran.

Patients and methods

Study population

In a population-based cohort, at least 10,000 people within the age range of 35–70 years old from Sheshdeh, the suburb of Fasa city and its 24 satellite villages were recruited. A detailed demographic, socioeconomic, anthropometric, nutrition, and medical history was obtained for each individual besides limited physical examinations and determination of physical activity and sleep patterns supplemented by body composition and electrocardiographic records (16). Participants older than 35 years with one of the following conditions were included in the study: alanine aminotransferase level >30 U/L, Gamma-glutamyltransferase >35 U/L.

Individuals with evidence of alcohol consumption of more than 20 grams per day, hepatitis, drug-induced liver disease, Wilson's disease, hereditary hemochromatosis, encephalopathy, or variceal bleeding were excluded from the study. Routine laboratory assessments were done, and a comprehensive biobank was compiled for future biological investigations. All data were stored online using a dedicated software.

Measurement of clinical parameters and biochemical analysis

Fatty liver index (FLI) is a standard predictor of hepatic steatosis severity. The following formula was designed and validated by Bogoni et al. for calculating FLI using laboratory and anthropometric data, including gamma glutamyl transferase concentration, triglyceride, Body mass index and waist circumference (17).

FLI = (e 0.953 × loge (TG) + 0.139 × BMI + 0.718 × loge (GGT) + 0.053 × WC – 15.745) / (1 + e 0.953 × loge (TG) + 0.139 × BMI + 0.718 × loge (GGT) + 0.053 × WC – 15.745) × 100

To calculate the amount of dairy intake, a validated food frequency questionnaire (FFQ) was used. Different dairy products consumption was recorded separately according to the amount of daily, weekly, monthly and annual intake units. The questionnaire included 100 items with a specified serving size, 7 response frequency categories from “almost never eat” to “twice or more per day” for foods, and 8 response frequency categories from “almost never drink” to “Four or more times a day” for beverages. Total daily energy and nutrient intake are calculated based on the latest available table of Iranian food composition. Frequency of consumption of dairy products (milk, yogurt, cheese, doogh, curd and cream and any other dairy products) in the last month was assessed using 7 categories of answers as follows: Almost never, less than once a week, once a week, 2–3 times a week, 4–6 times a week, once a day and twice a day. Detailed information about the design, foods included and the validity of this questionnaire has been reported elsewhere (18). Previous validation studies of this FFQ designed specifically for Iranian adults revealed good correlation between dietary intakes assessed by a similar FFQ and multiple days of 24-h dietary recalls (18, 19).

Statistical analysis

The characteristics of the participants across this study were showed in quartiles of FLI. Continuous variables were presented as mean ± SD and categorical variables as frequency and percentages (%). Data were analyzed using SPSS (version 18; SPSS Inc., Chicago, IL, USA) by running the student's t-test and one way ANOVA, and chi-squared tests. Odds ratios (ORs) and corresponding 95% confidence intervals (CIs) were estimated from univariate and multiple logistic regression for being classified with an undesirable level of FLI, ALT, AST and GGT and increasing quartile groups (Q1–Q4) for intake of dairy products. Multiple logistic regression models using backward elimination method were applied, adjusted for age, energy, physical activity, BMI, vitamin D levels, selenium levels, smoking status, hypertension, hypercholesterolemia, and diabetes mellitus, and for each relevant food items. Statistical tests were two-sided and P < 0.05 was considered statistically significant.

Results

Characteristics of study participants by sex and age

Characteristics of the study participants are presented in the Table 1. Data from 7,540 participants (59% female) were included in the final analysis. The mean age of all participants in this study was 48.81 ± 9.631 years. BMI was significantly higher in females as compared with males (26.63 ± 4.79 kg/m2 vs. 23.27 ± 4.17 kg/m2; P < 0.001). FLI was higher in females than males (39.991 ± 26.64 vs. 26.708 ± 23.39; P < 0.05). Also, serum levels of HDL, TC, LDL, FBS, SBP and DBP were significantly higher in females than males (P < 0.001). Dietary intakes of study participants across quartile of FLI score are presented in Table 2. A greater FLI score was significantly associated with the higher intake of carbohydrates, fruits, vegetables, meats, beans, dairy, vitamin D and selenium (P < 0.05).

Table 1

VariablesTotal (n = 7,540)Men (n = 3,048)Women (n = 4,492)P-value
Age48.81 ± 9.63149.48 ± 9.66648.36 ± 9.581<0.001
BMI (kg/m2)25.27 ± 4.8423.27 ± 4.1726.63 ± 4.79<0.001
Waist circumference (cm)92.274 ± 11.873487.483 ± 10.793095.526 ± 11.4660<0.001
Current smoker (%)1,985 (26.3 %)1,770 (58.1 %)215 (4.8 %)<0.001
Physical activities (Met-h/day)41.379 ± 11.07845.622 ± 14.34638.499 ± 6.772<0.001
Total cholesterol (mg/dl)182.70 ± 37.808174.03 ± 34.775188.58 ± 38.645<0.001
Triglycerides (mg/dl)122.89 ± 71.818121.63 ± 72.313123.74 ± 71.4760.210
HDL-cholesterol (mg/dl)51.66 ± 16.02647.64 ± 14.52554.39 ± 16.419<0.001
LDL-cholesterol (mg/dl)106.41 ± 31.882101.99 ± 29.820109.41 ± 32.874<0.001
Fasting glucose (mg/dl)90.90 ± 25.81588.66 ± 21.22092.42 ± 28.417<0.001
AST (U/l)20.20 ± 4.58421.07 ± 4.63319.60 ± 4.454<0.001
ALT (U/l)18.02 ± 5.38518.65 ± 5.67717.59 ± 5.133<0.001
GGT (U/l)17.08 ± 7.88218.58 ± 8.31116.07 ± 7.408<0.001
ALP (U/I)203.27 ± 62.184205.58 ± 58.175201.71 ± 64.7230.008
SBP (mmHg)110.69 ± 18.655109.19 ± 17.692111.70 ± 19.216<0.001
DBP (mmHg)74.08 ± 11.91673.15 ± 11.65474.72 ± 12.050<0.001
FLI34.621 ± 26.20226.708 ± 23.39839.991 ± 26.641<0.001

Baseline characteristics of study participants.

BMI, body mass index; HDL, high-density lipoprotein; LDL, low-density lipoprotein; AST, aspartate aminotransferase; ALT, alanine aminotransferase; GGT, gamma-glutamyl transferase; ALP, alkaline phosphatase; SBP, systolic blood pressure; DBP, diastolic blood pressure; FLI; fatty liver index.

Data are mean ± standard deviation or as n (%), The data analyzed using Student's t-test or chi-square test.

Table 2

Total FLI, quartiles
Quartile 1 (n = 1,885)Quartile 2 (n = 1,884)Quartile 3 (n = 1886)Quartile 4 (n = 1885)P-value
FLI6.395 ± 2.85219.195 ± 4.79940.314 ± 7.26372.569 ± 12.454<0.001
BMI, kg/m219.999 ± 2.21123.630 ± 2.13526.435 ± 2.22531.042 ± 3.783<0.001
Age47.49 ± 9.46948.94 ± 9.94849.16 ± 9.81249.65 ± 9.147<0.001
Smoking status, %779 (41.4%)490 (26.0%)407 (21.6%)309 (16.4%)<0.001
Diabetes, %81 (4.3%)173 (9.2%)237 (12.6%)339 (18.0%)<0.001
Hypertension, %147 (7.8%)301 (16.0%)410 (21.7%)651 (34.5%)<0.001
Hypercholesterolemia, %153 (8.1%)249 (13.2%)327 (17.3%)455 (24.1%)<0.001
physical activity (Met-h/day)44.523 ± 12.92242.266 ± 11.87340.320 ± 9.87538.407 ± 8.059<0.001
Energy, kcal2,949.230 ± 1,159.7842,976.097 ± 1,114.8902,886.361 ± 1,121.0602,909.666 ± 1,145.8910.071
Protein, g94.558 ± 40.84292.953 ± 37.82792.357 ± 38.29592.013 ± 38.6330.191
Carbohydrate, g501.27 ± 204.83490.19 ± 197.97491.97 ± 206.10470.85 ± 188.48<0.001
Fruit, g/d353.096 ± 316.174379.641 ± 302.302399.193 ± 318.931423.674 ± 335.956<0.001
Vegetables, g/d597.414 ± 338.962636.174 ± 331.090664.631 ± 350.200717.518 ± 403.156<0.001
Grains, g/d372.769 ± 204.430381.838 ± 212.351379.374 ± 211.520382.936 ± 211.9810.447
Meats and beans g/d150.558 ± 91.025145.492 ± 83.309140.905 ± 80.312140.087 ± 83.842<0.001
Dairy, g/d201.501 ± 181.891200.343 ± 162.967209.012 ± 177.231223.612 ± 184.515<0.001
Milk, g/d41.894 ± 64.49741.918 ± 63.98942.078 ± 79.58041.478 ± 69.0280.994
Cheese, g/d18.448 ± 19.93316.865 ± 16.24417.243 ± 17.72119.509 ± 21.458<0.001
Yogurt, g/d88.754 ± 92.84889.789 ± 87.99795.806 ± 94.337103.108 ± 100.118<0.001
Calcium, mg1,417.222 ± 712.9581,393.980 ± 657.0891,400.386 ± 671.8941,398.778 ± 663.3750.737
Riboflavin, mg2.453 ± 1.0752.399 ± 1.0022.378 ± 0.9922.393 ± 1.0050.120
Vitamin D, IUs1.435 ± 1.1781.412 ± 1.0301.348 ± 1.0161.345 ± 0.9750.014
Vitamin E, mg8.388 ± 4.1248.533 ± 4.0968.299 ± 3.7818.484 ± 4.0190.282
Selenium152.895 ± 77.523147.591 ± 73.298144.941 ± 73.542142.119 ± 72.438<0.001

Participant characteristic and dietary intake according to FLI quartile.

FLI; fatty liver index; BMI, body mass index.

Data are mean ± standard deviation or as n (%), The data analyzed using Student's t-test or chi-square test.

Association of fatty liver index with dairy products

Crude and multiple-adjusted models for mean intake of dairy intake across quartiles of FLI are shown in Table 3. After controlling for potential confounders, participants in the highest quartile of FLI didn't have higher odds for dairy products intake [(OR = 1.00; 95% CI: 0.99, 1.00; P = 0.084) for males and (OR = 1.00; 95% CI: 0.99, 1.01; P = 0.306) for females]. Also, both females and males in the highest quartile of FLI didn't have higher odds for consumption of milk, yogurt, cheese, doogh, curd and cream (P > 0.05). We also examined the association between ALT and AST quartiles and odds of dairy consumption. The results showed that individuals in the highest ALT quartile did not have odds of more dairy products consumption. Participants in the highest quartile of AST had higher odds only for cream intake (OR= 0.97; 95% CI: 0.95, 0.99; Ptrend = 0.011). There wasn't any significant association between other dairy products and ALT and AST quartiles (P > 0.05).

Table 3

SexCrude modelAdjusted modela
FLIQ1 (M = 1,112, W = 773)Q2 (M = 818, W = 1,066)Q3 (M = 649, W = 1,237)Q4 (M = 469, W = 1,416)P-trendQ1 (M = 1,112, W = 773)Q2 (M = 818, W = 1,066)Q3 (M = 649, W = 1,237)Q4 (M = 469, W = 1,416)P-trend
Dairy productsM1 (Ref)1.00 (0.99–1.00) P = 0.6111.00 (1.00–1.00) P = 0.5011.00 (1.00–1.00) P = 0.0840.2011 (Ref)1.00 (0.99–1.00) P = 0.3501.00 (0.99–1.00) P = 0.6861.00 (0.99–1.00) P = 0.7010.811
W1 (Ref)1.00 (1.00–1.00) P = 0.0361.00 (1.00–1.00) P = 0.0001.00 (1.00–1.00) P = 0.000<0.0011 (Ref)1.00 (1.00–1.00) P = 0.2911.00 (0.99–1.00) P = 0.4471.00 (0.99–1.01) P = 0.3060.671
MilkM1 (Ref)1.00 (0.99–1.00) P = 0.8551.00 (0.99–1.00) P = 0.8960.99 (0.99–1.00) P = 0.3160.7771 (Ref)1.00 (0.99–1.00) P = 0.9531.00 (0.99–1.00) P = 0.8130.99 (0.99–1.00) P = 0.1880.404
W1 (Ref)1.00 (0.99–1.00) P = 0.3131.00 (1.00–1.00) P = 0.1921.00 (1.00–1.00) P = 0.1130.4211 (Ref)1.00 (1.00–1.01) P = 0.1191.00 (0.99–1.00) P = 0.2631.00 (0.99–1.01) P = 0.3460.452
CheeseM1 (Ref)0.99 (0.99–0.99) P = 0.0300.99 (0.98–0.99) P = 0.0281.00 (0.99–1.01) P = 0.3180.0061 (Ref)0.99 (0.99–1.01) P = 0.2910.99 (0.98–1.01) P = 0.2151.01 (0.99–1.02) P = 0.2740.224
W1 (Ref)1.00 (0.99–1.01) P = 0.9911.00 (0.99–1.01) P = 0.1661.01 (1.00–1.02) P = 0.000<0.0011 (Ref)0.99 (0.98–1.01) P = 0.5681.00 (0.99–1.01) P = 0.9721.01 (0.99–1.02) P = 0.3100.126
YogurtM1 (Ref)1.00 (0.99–1.00) P = 0.6031.00 (1.00–1.00) P = 0.0571.00 (1.00–1.00) P = 0.0120.0471 (Ref)1.00 (0.99–1.00) P = 0.5931.00 (0.99–1.00) P = 0.9911.00 (0.99–1.00) P = 0.8610.890
W1 (Ref)1.00 (1.00–1.00) P = 0.1821.00 (1.00–1.00) P = 0.0011.00 (1.00–1.01) P = 0.000<0.0011 (Ref)0.99 (0.99–1.00) P = 0.5420.99 (0.99–1.00) P = 0.4981.00 (0.99–1.00) P = 0.8380.798
DoughM1 (Ref)1.00 (0.99–1.00) P = 0.4371.00 (0.99–1.00) P = 0.6361.00 (1.00–1.00) P = 0.1390.2291 (Ref)0.99 (0.99–1.00) P = 0.3211.00 (0.99–1.00) P = 0.9371.00 (0.99–1.00) P = 0.9720.667
W1 (Ref)1.00 (1.00–1.00) P = 0.1821.00 (1.00–1.00) P = 0.0011.00 (1.00–1.01) P = 0.000<0.0011 (Ref)1.00 (1.00–1.01) P = 0.0881.00 (0.99–1.00) P = 0.1441.00 (0.99–1.00) P = 0.1420.384
CurdM1 (Ref)0.98 (0.96–1.00) P = 0.0890.98 (0.96–1.01) P = 0.1291.01 (0.98–1.03) P = 0.4040.0581 (Ref)0.97 (0.95–1.01) P = 0.1460.98 (0.94–1.02) P = 0.3890.99 (0.95–1.04) P = 0.8780.391
W1 (Ref)0.99 (0.97–1.03) P = 0.4470.99 (0.97–1.01) P = 0.4611.00 (0.98–1.02) P = 0.9850.7101 (Ref)0.98 (0.95–1.01) P = 0.2950.99 (0.95–1.03) P = 0.4920.95 (0.89–1.01) P = 0.1010.262
CreamM1 (Ref)0.97 (0.95–0.99) P = 0.0420.98 (0.96–1.01) P = 0.1250.96 (0.93–0.99) P = 0.0220.0291 (Ref)1.01 (0.96–1.06) P = 0.6931.02 (0.97–1.07) P = 0.4841.02 (0.97–1.08) P = 0.4230.791
W1 (Ref)1.00 (0.98–1.03) P = 0.7541.01 (0.98–1.03) P = 0.6820.99 (0.97–1.02) P = 0.9320.9251 (Ref)0.99 (0.96–1.02) P = 0.5920.98 (0.94–1.03) P = 0.4600.99 (0.95–1.04) P = 0.9620.622
Dairy productsM1 (Ref)1.00 (1.00–1.00) P = 0.1711.00 (1.00–1.00) P = 0.6321.00 (1.00–1.00) P = 0.0020.0061 (Ref)1.00 (0.99–1.00) P = 0.8731.00 (0.99–1.00) P = 0.1531.00 (1.00–1.00) P = 0.3350.063
W1 (Ref)1.00 (1.00–1.00) P = 0.4361.00 (1.00–1.00) P = 0.2161.00 (1.00–1.00) P = 0.1780.5141 (Ref)1.00 (0.99–1.00) P = 0.8471.00 (0.99–1.00) P = 0.6651.00 (0.99–1.00) P = 0.7290.975
MilkM1 (Ref)1.00 (0.99–1.00) P = 0.6390.99 (0.99–1.00) P = 0.4071.00 (0.99–1.00) P = 0.3970.3181 (Ref)0.99 (0.99–1.00) P = 0.1340.99 (0.99–1.00) P = 0.0541.00 (0.99–1.00) P = 0.7030.107
W1 (Ref)1.00 (0.99–1.00) P = 0.7841.00 (0.99–1.00) P = 0.9341.00 (0.99–1.00) P = 0.2930.6691 (Ref)1.00 (0.99–1.00) P = 0.9711.00 (0.99–1.00) P = 0.5511.00 (0.99–1.00) P = 0.4880.641
CheeseM1 (Ref)1.00 (0.99–1.01) P = 0.0990.99 (0.99–1.01) P = 0.7751.01 (1.00–1.01) P = 0.0340.0301 (Ref)1.00 (0.99–1.01) P = 0.1830.99 (0.99–1.00) P = 0.4171.00 (0.99–1.01) P = 0.1640.053
W1 (Ref)0.99 (0.99–1.00) P = 0.5761.00 (0.99–1.00) P = 0.6971.01 (1.00–1.01) P = 0.0350.0521 (Ref)0.99 (0.99–1.00) P = 0.5541.00 (0.99–1.01) P = 0.9481.00 (0.99–1.01) P = 0.0860.122
YogurtM1 (Ref)1.00 (0.99–1.00) P = 0.5911.00 (0.99–1.00) P = 0.7471.00 (1.00–1.00) P = 0.0680.2361 (Ref)1.00 (0.99–1.01) P = 0.5180.99 (0.99–1.00) P = 0.2341.00 (0.99–1.00) P = 0.9820.489
W1 (Ref)1.00 (1.00–1.00) P = 0.2361.00 (1.00–1.00) P = 0.1731.00 (0.99–1.00) P = 0.8470.4381 (Ref)1.00 (0.99–1.00) P = 0.8151.00 (0.99–1.00) P = 0.8110.99 (0.99–1.00) P = 0.1390.324
DoughM1 (Ref)1.00 (1.00–1.01) P = 0.0051.00 (1.00–1.01) P = 0.0351.00 (1.00–1.01) P = 0.0000.0011 (Ref)1.00 (1.00–1.01) P = 0.1351.00 (0.99–1.00) P = 0.5291.00 (1.00–1.01) P = 0.0980.077
W1 (Ref)1.00 (0.99–1.00) P = 0.7891.00 (1.00–1.00) P = 0.2401.00 (1.00–1.00) P = 0.2410.5221 (Ref)1.00 (0.99–1.00) P = 0.6571.00 (0.99–1.00) P = 0.9811.00 (0.99–1.00) P = 0.9160.968
CurdM1 (Ref)0.99 (0.97–1.01) P = 0.6940.99 (0.97–1.02) P = 0.5381.01 (0.98–1.02) P = 0.4750.5081 (Ref)0.99 (0.97–1.02) P = 0.4800.98 (0.96–1.01) P = 0.3521.00 (0.98–1.02) P = 0.7690.481
W1 (Ref)0.99 (0.97–1.01) P = 0.3340.99 (0.98–1.01) P = 0.6131.00 (0.98–1.02) P = 0.9520.7231 (Ref)0.98 (0.97–1.01) P = 0.2200.99 (0.97–1.01) P = 0.3310.99 (0.97–1.02) P = 0.7700.579
CreamM1 (Ref)0.98 (0.96–1.01) P = 0.2540.98 (0.96–1.01) P = 0.1900.99 (0.96–1.01) P = 0.2670.5221 (Ref)0.98 (0.96–1.01) P = 0.3220.98 (0.96–1.01) P = 0.2520.99 (0.97–1.01) P = 0.4360.657
W1 (Ref)1.01 (0.98–1.02) P = 0.5700.99 (0.97–1.02) P = 0.9451.00 (0.98–1.02) P = 0.9060.9121 (Ref)1.01 (0.98–1.03) P = 0.4641.00 (0.97–1.02) P = 0.8721.01 (0.98–1.03) P = 0.5800.846
Dairy productsM1 (Ref)1.00 (0.99–1.00) P = 0.4391.00 (0.99–1.00) P = 0.6331.00 (1.00–1.00) P = 0.3180.1911 (Ref)1.00 (0.99–1.00) P = 0.2731.00 (0.99–1.00) P = 0.1631.00 (0.99–1.00) P = 0.9850.280
W1 (Ref)1.00 (0.99–1.00) P = 0.1591.00 (0.99–1.00) P = 0.5201.00 (0.99–1.00) P = 0.1320.3761 (Ref)1.00 (0.99–1.00) P = 0.2151.00 (0.99–1.00) P = 0.2891.00 (0.99–1.00) P = 0.0970.366
MilkM1 (Ref)0.99 (0.99–1.00) P = 0.1821.00 (0.99–1.00) P = 0.6311.00 (0.99–1.00) P = 0.5870.1991 (Ref)0.99 (0.99–1.00) P = 0.1370.99 (0.99–1.00) P = 0.2871.00 (0.99–1.00) P = 0.8660.323
W1 (Ref)1.00 (0.99–1.00) P = 0.6581.00 (0.99–1.00) P = 0.5581.00 (0.99–1.00) P = 0.9370.8031 (Ref)1.00 (0.99–1.00) P = 0.7081.00 (0.99–1.00) P = 0.7121.00 (0.99–1.00) P = 0.8310.897
CheeseM1 (Ref)1.00 (0.99–1.01) P = 0.6210.99 (0.99–1.00) P = 0.7611.00 (0.99–1.01) P = 0.2390.3621 (Ref)1.00 (0.99–1.01) P = 0.7470.99 (0.99–1.00) P = 0.5111.00 (0.99–1.01) P = 0.4840.450
W1 (Ref)0.99 (0.99–1.00) P = 0.0990.99 (0.99–1.00) P = 0.6870.99 (0.99–1.00) P = 0.1040.2251 (Ref)0.99 (0.99–1.00) P = 0.2151.00 (0.99–1.00) P = 0.8770.99 (0.99–1.01) P = 0.4760.496
YogurtM1 (Ref)1.00 (0.99–1.00) P = 0.8861.00 (0.99–1.00) P = 0.4171.00 (0.99–1.00) P = 0.5510.4541 (Ref)1.00 (0.99–1.00) P = 0.9640.99 (0.99–1.00) P = 0.1411.00 (0.99–1.00) P = 0.9170.297
W1 (Ref)1.00 (0.99–1.00) P = 0.5711.00 (0.99–1.00) P = 0.7720.99 (0.99–1.00) P = 0.0570.1801 (Ref)1.00 (0.99–1.00) P = 0.7131.00 (0.99–1.00) P = 0.5840.99 (0.99–1.00) P = 0.0500.199
DoughM1 (Ref)0.99 (0.99–1.00) P = 0.3131.00 (0.99–1.00) P = 0.4981.00 (0.99–1.00) P = 0.2750.1091 (Ref)0.99 (0.98–1.00) P = 0.2041.00 (0.99–1.00) P = 0.9631.00 (0.99–1.00) P = 0.7080.296
W1 (Ref)0.99 (0.99–1.00) P = 0.1180.99 (0.99–1.00) P = 0.1281.00 (0.99–1.00) P = 0.6550.2941 (Ref)0.99 (0.99–1.00) P = 0.1230.99 (0.99–1.00) P = 0.0551.00 (0.99–1.00) P = 0.3780.188
CurdM1 (Ref)1.00 (0.97–1.02) P = 0.9880.99 (0.96–1.01) P = 0.4220.99 (0.97–1.01) P = 0.5560.7841 (Ref)0.99 (0.97–1.02) P = 0.8760.98 (0.96–1.01) P = 0.2770.99 (0.96–1.01) P = 0.3750.612
W1 (Ref)0.99 (0.98–1.01) P = 0.7560.98 (0.96–1.00) P = 0.1591.00 (0.98–1.02) P = 0.9840.4801 (Ref)0.99 (0.98–1.01) P = 0.7440.98 (0.95–1.00) P = 0.0570.99 (0.97–1.01) P = 0.5760.236
CreamM1 (Ref)0.98 (0.97–1.01) P = 0.3290.96 (0.93–0.98) P = 0.0050.97 (0.95–0.99) P = 0.0310.0171 (Ref)0.98 (0.96–1.01) P = 0.2930.96 (0.94–0.98) P = 0.0040.97 (0.95–0.99) P = 0.0200.011
W1 (Ref)1.01 (0.98–1.03) P = 0.3721.01 (0.98–1.02) P = 0.9601.01 (0.98–1.03) P = 0.6020.8131 (Ref)1.01 (0.98–1.04) P = 0.3121.01 (0.98–1.03) P = 0.5511.01 (0.98–1.03) P = 0.4720.726

Odds ratio (95% CI) from multiple logistic regression by dairy food type.

W, women; M, men.

Odds ratio (95% CI) are from multiple logistic regression for being classified with an undesirable level of FLI, ALT, AST and GGT and increasing quintile groups (Q1–Q4) for intake of dairy products. Q1, lowest intake, is the reference category. Statistically significant p-values are in supers.

a

Multivariable-adjusted model adjusted for age (years), energy (kcals/day), physical activity (minutes/week), BMI (kg/m2), Vitamin D levels, Selenium levels, smoking status (yes, no), hypertension (yes, no), hypercholesterolemia (yes, no), and diabetes (yes, no), and for each relevant food items.

The results of univariate and multiple analyses for determining independent association between fatty liver index and dairy consumption are shown in the Table 4. The mean FLI in the tertile 3 milk intake was 34.08 (95% CI: 32.99, 35.17) that wasn't significantly different from tertile 1 (P = 0.171). Mean liver enzymes did not differ significantly in different milk consumption tertiles (P = 0.08 for AST, P = 0.174 for ALT, and P = 0.491 for GGT).

Table 4

Milk
Tertile 1 (n = 2,952)Tertile 2 (n = 2,368)Tertile 3 (n = 2,220)
Mean (95% CI)P–value ANOVA
FLI35.32 (34.35–36.29)34.25 (33.23–35.27)34.08 (32.99–35.17)0.171
AST20.07 (19.91–20.24)20.20 (20.02–20.39)20.35 (20.16–20.55)0.088
ALT17.96 (17.77–18.15)17.92 (17.71–18.14)18.20 (17.97–18.42)0.174
GGT17.21 (16.92–17.51)16.97 (16.68–17.25)17.03 (16.68–17.37)0.491
OR (95% CI)a
FLI1 (Ref)0.99 (0.99–1.00)0.2300.96 (0.94–0.99)0.025
AST1 (Ref)1.00 (0.99–1.02)0.6081.00 (0.99–1.02)0.223
ALT1 (Ref)0.99 (0.98–1.01)0.5291.00 (0.99–1.01)0.826
GGT1 (Ref)0.99 (0.98–1.00)0.2600.99 (0.99–1.00)0.439
Cheese
Tertile 1 (n =3,075)Tertile 2 (n =1,996)Tertile 3 (n =2,469)
Mean (95% CI)P–value ANOVA
FLI35.55 (32.65–34.46)35.73 (34.57–36.89)35.04 (33.98–36.10)0.010
AST20.14 (19.98–20.30)20.22 (20.03–20.42)20.24 (20.06–20.43)0.678
ALT17.75 (17.56–17.94)18.16 (17.93–18.40)18.24 (18.02–18.45)0.001
GGT16.84 (16.59–17.08)17.32 (16.96–17.67)17.19 (16.85–17.54)0.071
OR (95% CI)a
FLI1 (Ref)1.00 (0.99–1.01)0.1371.00 (0.99–1.01)0.936
AST1 (Ref)1.01 (0.99–1.02)0.3261.00 (0.98–1.01)0.850
ALT1 (Ref)1.01 (0.99–1.02)0.1191.00 (0.99–1.01)0.373
GGT1 (Ref)1.01 (0.99–1.01)0.1101.00 (0.99–1.01)0.598
Yogurt
Tertile 1 (n =2,514)Tertile 2 (n =2,938)Tertile 3 (n =2,088)
Mean (95% CI)P–value ANOVA
FLI32.73 (31.71–33.74)34.75 (33.81–35.69)36.70 (35.56–37.85)0.000
AST20.13 (19.96–20.31)20.37 (20.20–20.54)20.03 (19.84–20.22)0.026
ALT17.84 (17.63–18.06)18.10 (17.90–18.30)18.11 (17.88–18.34)0.143
GGT16.86 (16.58–17.13)17.05 (16.75–17.35)17.39 (17.03–17.74)0.073
OR (95% CI)a
FLI1 (Ref)1.00 (0.99–1.00)0.5771.00 (0.99–1.01)0.879
AST1 (Ref)1.01 (0.99–1.02)0.2010.98 (0.97–1.00)0.073
ALT1 (Ref)0.99 (0.98–1.01)0.8150.99 (0.97–1.00)0.099
GGT1 (Ref)0.96 (1.00–0.99)0.9641.00 (0.99–1.01)0.523

Univariate and multiple analyses to determine the independent association between fatty liver index and dairy consumption.

a

Multivariable-adjusted model adjusted for sex, age (years), energy (kcals/day), physical activity (minutes/week), BMI (kg/m2), Vitamin D levels, Selenium levels, smoking status (yes, no), hypertension (yes, no), hypercholesterolemia (yes, no), and diabetes (yes, no), and for each relevant food items.

The association between the amount of dairy consumption and severity of fatty liver disease is shown in Table 4. When controlled for confounding factors, multiple logistic regression analysis demonstrated that the amount of milk consumption was the independent preventive predictor of FLI (OR = 0.96; 95% CI: 0.94, 0.99; P = 0.025), but was not the predictors of higher levels of AST, ALT and GGT (P > 0.05). For the cheese intake, participants in the tertile 3 had significantly lower FLI than the lowest tertile 35.04 (95% CI: 33.98, 36.10; P = 0.01). Individuals in the highest tertile of cheese intake had higher levels of ALT than lowest tertile 18.24 mg/dl (95% CI: 18.02, 18.45; P = 0.001). However, there wasn't any significant association between cheese intake and odds of fatty liver index (P > 0.05).

Participants in the highest tertile of yogurt intake had significantly higher FLI scores 36.70 (95% CI: 35.56, 37.85; P < 0.001) and AST concentration 20.03 mg/dl (95% CI: 19.84, 20.22; P = 0.026). We did not find a significant association between yogurt consumption and fatty liver index or biomarkers (P > 0.05).

Discussion

The results of this study clearly indicate that participants with higher FLI score consume more dairy, cheese and yogurt compared to the subjects with lower FLI. Moreover, there is no strong correlation between dairy consumption and fatty liver indicators and only milk consumption is a preventive indicator of FLI. Also, in our study, the consumption of dairy products in women had more protective effects against NAFLD than in men.

The impacts of dairy products on the amounts of liver-derived lipoproteins have been widely studied and discussed previously (20, 21). Several studies have looked at dietary factors, including dairy intake as the key interventional variable on NAFLD-related parameters in NAFLD patients. In some observational studies, researchers have reported an inverse association between low-fat dairy intake and the severity of NAFLD. Ferolla et al. in a cross-sectional study of 96 patients with NAFLD observed an inverse association between low fat dairy intake and NAFLD severity (22). Similar findings were found in another study (23). Conversely, it has been shown that some dietary regiments which contain adequate amounts of dairy products such as DASH (Dietary Approaches to Stop Hypertension) can have a positive effect on reducing the symptoms of NAFLD. Razavi Zade et al. conducted a clinical trial with 3 daily servings with low dairy products to evaluate the effects of DASH diet on the NAFLD severity and showed adherence to DASH eating pattern leads to a significant reduction in serum levels of ALT, AST, and insulin resistance (24). However, part of the positive effects of DASH diet on fatty liver can be attributed to high amounts of fibers and phytochemicals in the DASH diet (25).

In the present study, except for milk and FLI correlation, we did not observe any correlations between dairy products and fatty liver biomarkers. Contrary to our findings, in a 6-week randomized controlled trial, Dugan et al. found a significant reduction in ALT, AST, hepatic steatosis index, and mRNA expression of IL-6 and IL-1β in peripheral blood mononuclear cells (PBMCs), following the consumption of 3 daily servings of low-fat dairy (296 mL 1% milk, 170 g non-fat yogurt, 56.7 g 2% cheese) as compared with isocaloric carbohydrate-based control foods (26). A recent meta-analysis was done on the association between food groups and the risk of NAFLD which included three cross sectional studies and one case control study. The cross-sectional studies didn't find a significant correlation between dairy consumption and the likelihood of NAFLD. However, the case–control study showed a positive association between dairy product consumption and the possibility of NAFLD (27).

The diversity in the type of dairy products and their fat percentages may partly account for the difference observed in the results. Most previous studies that have reported significant changes on fatty liver related indicators following dairy consumption have examined the effect of low-fat dairy products, while in our study participants consumed both high-fat and low-fat dairy products. In fact, the diet of Iranian people especially those in rural areas is mostly based on high consumption of high-fat dairy products (28). Esmaillzadeh et al. studied 486 healthy Iranian women aged 40–60 years and found an inverse association between low fat dairy intake and C-reactive protein, IL-6 and soluble vascular cell adhesion molecule (29). Moreover, they reported an increase in serum amyloid A and soluble vascular cell adhesion molecule-1 among women with high fat dairy consumption (30). High-fat dairy products contain high levels of saturated fatty acids which play a key role in dyslipidemia, insulin resistance and inflammation, that altogether are major risk factors for NAFLD (31).

In contrast with our findings, Zhang et al. in a cross-sectional study of 24,389 adults found that participants with yogurt consumption more than 4 times per weeks had lower likelihood for NAFLD (32). Other similar studies have shown a preventative role for yogurt consumption in the development of other chronic diseases such as cardiovascular disease and type 2 diabetes (29, 33). Most studies that have shown the positive effect of yogurt consumption on liver enzymes have used a variety of probiotic yogurts (9, 34).

In our study, participants in the higher FLI quartile consumed significantly more yogurt than subjects in the lower quartile of FLI. The present study was a cross-sectional study and due to the nature of cross-sectional studies it is not possible to identify the causal relationship between two variables. However, it seems that this difference might be attributed to the fact that people who are at risk of NAFLD increase their yogurt consumption to get more probiotics because of the positive role of bacteria in yogurt against NAFLD, especially probiotic yogurt. The positive effect of yogurt on fatty liver has been attributed to the presence of various probiotic compounds in yogurt. Animal studies have shown that some bacterial strains, such as Lactobacillus, have the ability to reduce the inflammation caused by high levels of lipopolysaccharide (LPS) and subsequent hepatic toll-like receptor 4 (TLR4) activation in NAFLD patients (35). Additionally, part of the beneficial effects of dairy products, especially low-fat dairy products, might be due to high amounts of nutrients, especially magnesium, potassium, protein and calcium in dairy products, which in turn can increase the whole-body fat oxidation (36, 37).

In the present study, we conducted sex-specific analyses which showed that the consumption of dairy products in women had stronger protective effects against NAFLD than in men. While the exact mechanisms for this difference are unclear, several studies have reported a higher risk of NALFD in men and post-menopausal women than in pre-menopausal women (38, 39). Interestingly, experimental results suggest a protective role for estradiol in hepatic injuries by suppressing lipid accumulation and liver fibrosis (40). The presence of estrone (E1) and 17β-estradiol (E2) in raw whole cow's milk has been demonstrated (41). Dairy products have been estimated to account for up to 60% of estrogens in a German diet (42). Some studies have shown that the estrogen in dairy products may increase the concentration of estrogen in the serum (43).

The results of some studies contradicted with the findings of our study. Lee et al. in a study of 5171 adults showed higher dairy protein intake was significantly and inversely associated with the risk of incident NAFLD in men and women aged ≥50 years (44). They suggested that replacing macronutrients equivalent to carbohydrate intake with dairy protein intake may have contributed to lowering the risk of developing NAFLD (32, 44). It has also been suggested that some of the beneficial effects of dairy products in preventing NAFLD are related to the whey protein found in dairy products. Several human and animal studies have shown that whey protein or dairy products reduce weight and fat mass (4547). Other mechanisms have been suggested for the beneficial effects of dairy products against NAFLD, some of which are related to insulin resistance. The population-based prospective Coronary Artery Risk Development in Young Adults study found an inverse association between frequency of dairy intake and development of insulin resistance syndrome (48). Also, a meta-analysis study showed that the consumption of dairy products, especially low-fat dairy products had a beneficial effect on Homeostatic Model Assessment for Insulin Resistance (HOMA-IR), waist circumference, and body weight (49).

This study has several important advantages including the large sample size, and extensive information on lifestyle and dietary factors, which allowed us to control for many potential confounders. In this study, we adjusted the results for several variables, especially BMI which is a crude measure of body fat and is associated with NAFLD and may negatively affect the accuracy of the results (32). On the other hand, the present study was not free of limitations. Due to the structure of the questionnaires used, we were not able to examine low-fat and high-fat dairy products separately. Also, because of the nature of questionnaires filled in interviews, there is always possibility of bias which can affect the accuracy of the results. Another limitation of the current study was the presence of some potential confounding factors such as diabetes, hyperlipidemia, and high blood pressure and their variability between FLI quartile. However, in the adjusted models, we adjusted the effects of these confounding factors. Finally, because of the difficulties in budget management and the high sample size, we were not able to assess the liver through imaging studies. Evaluation of hepatic steatosis and fibrosis could increase the accuracy of the results in future studies.

Conclusion

In summary, this population-based cohort study didn't show any strong association between dairy products consumption and NAFLD indicators. However, a modest inverse correlation was observed between milk consumption and FLI. Therefore, based on the results of this study, consumption of appropriate amounts of milk, especially low-fat types (at least one unit of milk more than 5–6 times a week) can play a positive role in preventing NAFLD. The results of the present study seem to be in line with the existing recommendations for consumption of a healthy dietary pattern to prevent NAFLD.

Funding

This study was supported by the Deputy of Research and Technology of Fasa University of Medical Sciences, Fasa, Iran (No. 400067).

Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Statements

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The study protocol was approved by the Research Council and the Ethics Committee of Fasa University of Medical Sciences (IR.FUMS.REC.1400.085). The patients/participants provided their written informed consent to participate in this study.

Author contributions

ZK, MF, RH, and RT make substantial contributions to conception, design, acquisition of data, analysis and interpretation of data, and drafting of this study. ZK, MR, AA, and MK participate in drafting the article or revising it. MF and RH gives final approval of the version to be submitted and any revised version. MF and RT are the guarantors of this work. All authors have read and approved the manuscript.

Acknowledgments

The authors thank all the study participants for their valuable contributions. This paper was derived from a thesis for fulfillment of obtaining a MD degree.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

    Abbreviations

  • NAFLD

    Nonalcoholic fatty liver disease

  • FLI

    Fatty liver index

  • OR

    Odds ratio

  • FFQ

    Food frequency questionnaire

  • BMI

    body mass index

  • HDL

    high-density lipoprotein

  • LDL

    low-density lipoprotein

  • AST

    aspartate aminotransferase

  • ALT

    alanine aminotransferase

  • GGT

    gamma-glutamyl transferase

  • ALP

    alkaline phosphatase

  • SBP

    systolic blood pressure

  • DBP

    diastolic blood pressure

  • ANOVA

    One-way Analysis of Variance.

References

  • 1.

    MatteoniCAYounossiZMGramlichTBoparaiNLiuYCMcculloughAJ. Nonalcoholic fatty liver disease: a spectrum of clinical and pathological severity. Gastroenterology. (1999) 116:14139. 10.1016/S0016-5085(99)70506-8

  • 2.

    Organization WH. Global Status Report on Noncommunicable Diseases 2014. Geneva: World Health Organization (2014).

  • 3.

    RezendeACSouzaLGJardimTVPerilloNBAraújoYCLDe SouzaSGet al. Is waist-to-height ratio the best predictive indicator of hypertension incidence? A cohort study. BMC Public Health. (2018) 18:111. 10.1186/s12889-018-5177-3

  • 4.

    YounossiZMKoenigABAbdelatifDFazelYHenryLWymerM. Global epidemiology of nonalcoholic fatty liver disease—meta-analytic assessment of prevalence, incidence, and outcomes. Hepatology. (2016) 64:7384. 10.1002/hep.28431

  • 5.

    WongS-WChanW-K. Epidemiology of non-alcoholic fatty liver disease in Asia. Indian J Gastroenterol. (2020) 39:18. 10.1007/s12664-020-01018-x

  • 6.

    AshtariSPourhoseingholiMAZaliMR. Non-alcohol fatty liver disease in Asia: Prevention and planning. World J Hepatol. (2015) 7:1788. 10.4254/wjh.v7.i13.1788

  • 7.

    MantovaniAScorlettiEMoscaAAlisiAByrneCDTargherG. Complications, morbidity and mortality of nonalcoholic fatty liver disease. Metabolism. (2020) 111:154170. 10.1016/j.metabol.2020.154170

  • 8.

    KimD-HKimHJeongDKangI-BChonJ-WKimH-Set al. Kefir alleviates obesity and hepatic steatosis in high-fat diet-fed mice by modulation of gut microbiota and mycobiota: targeted and untargeted community analysis with correlation of biomarkers. J Nutr Biochem. (2017) 44:3543. 10.1016/j.jnutbio.2017.02.014

  • 9.

    BakhshimoghaddamFShateriKSinaMHashemianMAlizadehM. Daily consumption of synbiotic yogurt decreases liver steatosis in patients with nonalcoholic fatty liver disease: a randomized controlled clinical trial. J Nutr. (2018) 148:127684. 10.1093/jn/nxy088

  • 10.

    ChenYFengRYangXDaiJHuangMJiXet al. Yogurt improves insulin resistance and liver fat in obese women with nonalcoholic fatty liver disease and metabolic syndrome: a randomized controlled trial. Am J Clin Nutr. (2019) 109:16119. 10.1093/ajcn/nqy358

  • 11.

    KongCYLiZMHanBZhangZYChenHLZhangSLet al. Diet consisting of balanced yogurt, fruit, and vegetables modifies the gut microbiota and protects mice against nonalcoholic fatty liver disease. Mol Nutr Food Res. (2019) 63:1900249. 10.1002/mnfr.201900249

  • 12.

    EgresiADrexlerDHagymásiKBlázovicsAJakabZKocsisIet al. The potential role of organic and conventional yoghurt consumption in the treatment of non-alcoholic fatty liver disease. Orv Hetil. (2020) 161:146674. 10.1556/650.2020.31839

  • 13.

    Colín-RamírezERivera-MancíaSInfante-VázquezOCartas-RosadoRVargas-BarrónJMaderoMet al. Protocol for a prospective longitudinal study of risk factors for hypertension incidence in a Mexico City population: the Tlalpan 2020 cohort. BMJ Open. (2017) 7:e016773. 10.1136/bmjopen-2017-016773

  • 14.

    HøstmarkATTomtenSE. The Oslo Health Study: cheese intake was negatively associated with the metabolic syndrome. J Am Coll Nutr. (2011) 30:18290. 10.1080/07315724.2011.10719959

  • 15.

    GuoZLiuXZhangQShenZTianFZhangHet al. Influence of consumption of probiotics on the plasma lipid profile: a meta-analysis of randomised controlled trials. Nutr Metabol Cardiovasc Dis. (2011) 21:84450. 10.1016/j.numecd.2011.04.008

  • 16.

    FarjamMBahramiHBahramaliEJamshidiJAskariAZakeriHet al. A cohort study protocol to analyze the predisposing factors to common chronic non-communicable diseases in rural areas: Fasa Cohort Study. BMC Public Health. (2016) 16:18. 10.1186/s12889-016-3760-z

  • 17.

    KraussRMWinstonMFletcherBJGrundySM. Obesity: impact on cardiovascular disease. Circulation. (1998) 98:14726. 10.1161/01.CIR.98.14.1472

  • 18.

    KeshteliAHEsmaillzadehARajaieSAskariGFeinle-BissetCAdibiP. A dish-based semi-quantitative food frequency questionnaire for assessment of dietary intakes in epidemiologic studies in Iran: design and development. Int J Prev Med. (2014) 5:29. Available online at: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3915470/

  • 19.

    MalekshahAKimiagarMSaadatian-ElahiMPourshamsANouraieMGoglaniGet al. Validity and reliability of a new food frequency questionnaire compared to 24 h recalls and biochemical measurements: pilot phase of Golestan cohort study of esophageal cancer. Eur J Clin Nutr. (2006) 60:9717. 10.1038/sj.ejcn.1602407

  • 20.

    HuthPJParkKM. Influence of dairy product and milk fat consumption on cardiovascular disease risk: a review of the evidence. Adv Nutr. (2012) 3:26685. 10.3945/an.112.002030

  • 21.

    AstrupA. Yogurt and dairy product consumption to prevent cardiometabolic diseases: epidemiologic and experimental studies. Am J Clin Nutr 99 1235s-1242s. (2014). 10.3945/ajcn.113.073015

  • 22.

    FerollaSMFerrariTCLimaMLReisTOTavaresWCJrCoutoOFet al. Dietary patterns in Brazilian patients with nonalcoholic fatty liver disease: a cross-sectional study. Clinics. (2013) 68:117. 10.6061/clinics/2013(01)OA03

  • 23.

    LeiSLiuZWYunLCaiGZhangHSongLJet al. The prevalence of nonalcoholic fatty liver disease and its association with lifestyle/dietary habits among university faculty and staff in Chengdu. Biomed Environ Sci. (2012) 25:38391. 10.3967/0895-3988.2012.04.002

  • 24.

    Razavi ZadeMTelkabadiMHBahmaniFSalehiBFarshbafSAsemiZ. The effects of DASH diet on weight loss and metabolic status in adults with non-alcoholic fatty liver disease: a randomized clinical trial. Liver Int. (2016) 36:56371. 10.1111/liv.12990

  • 25.

    DreherML. Dietary patterns, foods, nutrients and phytochemicals in non-alcoholic fatty liver disease. In: Dietary Patterns and Whole Plant Foods in Aging and Disease. New York, NY: Springer (2018). p. 291311.

  • 26.

    DuganCEAguilarDParkYKLeeJYFernandezML. Dairy consumption lowers systemic inflammation and liver enzymes in typically low-dairy consumers with clinical characteristics of metabolic syndrome. J Am Coll Nutr. (2016) 35:25561. 10.1080/07315724.2015.1022637

  • 27.

    HeKLiYGuoXZhongLTangS. Food groups and the likelihood of non-alcoholic fatty liver disease: a systematic review and meta-analysis. Br J Nutr. (2020) 124:113. 10.1017/S0007114520000914

  • 28.

    Salehi-SahlabadiASadatSBeigrezaeiSPourmasomiMFeiziAGhiasvandRet al. Dietary patterns and risk of non-alcoholic fatty liver disease. BMC Gastroenterol. (2021) 21:112. 10.1186/s12876-021-01612-z

  • 29.

    ChenMSunQGiovannucciEMozaffarianDMansonJEWillettWCet al. Dairy consumption and risk of type 2 diabetes: 3 cohorts of US adults and an updated meta-analysis. BMC Med. (2014) 12:215. 10.1186/s12916-014-0215-1

  • 30.

    EsmaillzadehAAzadbakhtL. Dairy consumption and circulating levels of inflammatory markers among Iranian women. Public Health Nutr. (2010) 13:1395402. 10.1017/S1368980009992126

  • 31.

    FanJGCaoHX. Role of diet and nutritional management in non-alcoholic fatty liver disease. J Gastroenterol Hepatol. (2013) 28:817. 10.1111/jgh.12244

  • 32.

    ZhangSFuJZhangQLiuLLuMMengGet al. Association between habitual yogurt consumption and newly diagnosed non-alcoholic fatty liver disease. Eur J Clin Nutr. (2020) 74:4919. 10.1038/s41430-019-0497-7

  • 33.

    BuendiaJRLiYHuFBCabralHJBradleeMLQuatromoniPAet al. Regular yogurt intake and risk of cardiovascular disease among hypertensive adults. Am J Hypertens. (2018) 31:55765. 10.1093/ajh/hpx220

  • 34.

    NabaviSRafrafMSomiMHHomayouni-RadAAsghari-JafarabadiM. Effects of probiotic yogurt consumption on metabolic factors in individuals with nonalcoholic fatty liver disease. J Dairy Sci. (2014) 97:738693. 10.3168/jds.2014-8500

  • 35.

    XueLHeJGaoNLuXLiMWuXet al. Probiotics may delay the progression of nonalcoholic fatty liver disease by restoring the gut microbiota structure and improving intestinal endotoxemia. Sci Rep. (2017) 7:113. 10.1038/srep45176

  • 36.

    MelansonESharpTSchneiderJDonahooWGrunwaldGHillJ. Relation between calcium intake and fat oxidation in adult humans. Int J Obes. (2003) 27:196203. 10.1038/sj.ijo.802202

  • 37.

    ZemelMB. Proposed role of calcium and dairy food components in weight management and metabolic health. Phys Sportsmed. (2009) 37:2939. 10.3810/psm.2009.06.1707

  • 38.

    YangJDAbdelmalekMFPangHGuyCDSmithADDiehlAMet al. Gender and menopause impact severity of fibrosis among patients with nonalcoholic steatohepatitis. Hepatology. (2014) 59:140614. 10.1002/hep.26761

  • 39.

    LonardoANascimbeniFBallestriSFairweatherDWinSThanTAet al. Sex differences in nonalcoholic fatty liver disease: state of the art and identification of research gaps. Hepatology. (2019) 70:145769. 10.1002/hep.30626

  • 40.

    LeeCKimJJungY. Potential therapeutic application of estrogen in gender disparity of nonalcoholic fatty liver disease/nonalcoholic steatohepatitis. Cells. (2019) 8:1259. 10.3390/cells8101259

  • 41.

    AbshiriniMSiassiFKoohdaniFQorbaniMGolpour-HamedaniSKhosraviSet al. Association between dairy consumption and menopausal symptoms: A cross-sectional study among Iranian postmenopausal women. Int Dairy J. (2020) 105:104688. 10.1016/j.idairyj.2020.104688

  • 42.

    HartmannSLacornMSteinhartH. Natural occurrence of steroid hormones in food. Food Chem. (1998) 62:720. 10.1016/S0308-8146(97)00150-7

  • 43.

    Pape-ZambitoDRobertsRKensingerR. Estrone and 17β-estradiol concentrations in pasteurized-homogenized milk and commercial dairy products. J Dairy Sci. (2010) 93:253340. 10.3168/jds.2009-2947

  • 44.

    LeeJHLeeHSAhnSBKwonYJ. Dairy protein intake is inversely related to development of non-alcoholic fatty liver disease. Clin Nutr. (2021) 40:525260. 10.1016/j.clnu.2021.08.012

  • 45.

    FrestedtJLZenkJLKuskowskiMAWardLSBastianED. A whey-protein supplement increases fat loss and spares lean muscle in obese subjects: a randomized human clinical study. Nutr Metab. (2008) 5:17. 10.1186/1743-7075-5-8

  • 46.

    EllerLKReimerRA. Dairy protein attenuates weight gain in obese rats better than whey or casein alone. Obesity. (2010) 18:70411. 10.1038/oby.2009.300

  • 47.

    TrottierSKMacphersonREKnuthCMTownsendLKPepplerWTMikhaeilJSet al. Dairy attenuates weight gain to a similar extent as exercise in rats fed a high-fat, high-sugar diet. Obesity. (2017) 25:170715. 10.1002/oby.21941

  • 48.

    PereiraMAJacobsDRJrVan HornLSlatteryMLKartashovAILudwigDS. Dairy consumption, obesity, and the insulin resistance syndrome in young adults: the CARDIA study. JAMA. (2002) 287:20819. 10.1001/jama.287.16.2081

  • 49.

    SocholKMJohnsTSButtarRSRandhawaLSanchezEGalMet al. The effects of dairy intake on insulin resistance: a systematic review and meta-analysis of randomized clinical trials. Nutrients. (2019) 11:2237. 10.3390/nu11092237

Summary

Keywords

nonalcoholic fatty liver diseases, dairy, yogurt, cheese, cohort

Citation

Keshavarz Z, Rahimlou M, Farjam M, Homayounfar R, Khodadost M, Abdollahi A and Tabrizi R (2022) Non-alcoholic fatty liver disease and dairy products consumption: Results from FASA Persian cohort study. Front. Nutr. 9:962834. doi: 10.3389/fnut.2022.962834

Received

06 June 2022

Accepted

22 August 2022

Published

09 September 2022

Volume

9 - 2022

Edited by

Zachary Clayton, University of Colorado Boulder, United States

Reviewed by

Mee Young Hong, San Diego State University, United States; Devin Wahl, Colorado State University, United States

Updates

Copyright

*Correspondence: Reza Tabrizi Reza Homayounfar

†These authors have contributed equally to this work

This article was submitted to Nutrition and Metabolism, a section of the journal Frontiers in Nutrition

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics