Dairy Fat Intake, Plasma Pentadecanoic Acid, and Plasma Iso-heptadecanoic Acid Are Inversely Associated With Liver Fat in Children : Journal of Pediatric Gastroenterology and Nutrition

Secondary Logo

Journal Logo

Original Articles: Hepatology

Dairy Fat Intake, Plasma Pentadecanoic Acid, and Plasma Iso-heptadecanoic Acid Are Inversely Associated With Liver Fat in Children

Sawh, Mary Catherine∗,†; Wallace, Martina; Shapiro, Emma∗,§; Goyal, Nidhi P.∗,†; Newton, Kimberly P.∗,†; Yu, Elizabeth L.∗,†; Bross, Craig; Durelle, Janis; Knott, Cynthia||; Gangoiti, Jon A.; Barshop, Bruce A.; Gengatharan, Jivani M.; Meurs, Noah; Schlein, Alexandra#; Middleton, Michael S.#; Sirlin, Claude B.#; Metallo, Christian M.; Schwimmer, Jeffrey B.∗,†

Author Information
Journal of Pediatric Gastroenterology and Nutrition: April 2021 - Volume 72 - Issue 4 - p e90-e96
doi: 10.1097/MPG.0000000000003040


What Is Known/What Is New

What Is Known

  • Pediatric nutrition guidelines recommend that children 2 years and older consume low-fat or nonfat milk to reduce cardiovascular disease and limit excessive weight gain.
  • These guidelines do not directly address nonalcoholic fatty liver disease.

What Is New

  • In children at risk for nonalcoholic fatty liver disease, the average daily intake of dairy fat was significantly inversely correlated with the amount of liver fat.
  • Higher amounts of plasma fatty acids associated with dairy fat, specifically the odd chain fatty acid, pentadecanoic acid, and the branched chain fatty acid, iso-heptadecanoic acid, were associated with lower amounts of liver fat.

Pediatric dietary recommendations regarding dairy fat do not directly address nonalcoholic fatty liver disease (NAFLD). This is a critical gap considering NAFLD is the most common chronic liver disease in children and is thought to be strongly influenced by diet. The American Academy of Pediatrics (AAP) recommends that children at least 2 years consume low-fat or nonfat milk to reduce cardiovascular disease and limit excessive weight gain (1). Longitudinal studies, however, contradict this guideline, demonstrating that in children, consumption of nonfat milk is associated with greater weight gain than is observed with consumption of whole fat milk (2). The role of dietary fat in the development and treatment of NAFLD is controversial. Some pediatric gastroenterologists recommend decreasing dietary fat intake, which typically includes eliminating whole milk and other sources of dairy fat (1); however, dairy fat is not simply a source of cholesterol or calories but rather a complex mixture of bioactive fatty acids. This mixture includes odd-chain fatty acids (OCFA) and monomethyl branched chain fatty acids (BCFA), which may have a beneficial effect on hepatic steatosis and type 2 diabetes (3).

Observational studies in adults demonstrate an association between whole fat dairy consumption and multiple health benefits (4). Several authors suggest that OCFA, such as pentadecanoic acid (C15:0) and heptadecanoic acid (C17:0) mediate the beneficial effects of dairy fat (3). These OCFA are biomarkers of dairy fat intake and in adults are also associated with a lower incidence of type 2 diabetes. Furthermore, authors have speculated that the link between OCFA and decreased diabetes may be through OCFA-mediated reduction in liver fat (5). In addition to OCFA, BCFA, such as iso-heptadecanoic acid (iso-C17:0), are present in dairy fat constituting 2% of total fatty acids in cow's milk. BCFA are saturated fatty acids containing methyl-branches in the iso- or anteiso-positions. The plasma levels of BCFA are lower in adults with obesity (6) but their levels in relation to hepatic steatosis are unknown.

We sought to evaluate the relevance of the pediatric dairy fat recommendation by studying whether dairy fat intake is associated with hepatic steatosis in children at risk for NAFLD. We designed a prospective observational study with the following aims: to evaluate the association between dairy fat intake and the amount of liver fat, and to evaluate the association between plasma levels of dairy-associated OCFA and BCFA with the amount of liver fat. We hypothesized that higher intake of dairy fat is associated with lower liver fat and that plasma OCFAs and BCFAs correlate negatively with liver fat.


Study Cohort

This was an observational, cross-sectional study of children at risk for NAFLD, ages 8 through 17. Participants were recruited from primary care offices and community health centers from July 2015 through December 2017. In order to focus on children at risk for NAFLD, we oversampled male individuals and children with obesity (7–9). The oversampling approach was to enroll a study population that was constituted of 60% to 65% male individuals and that at least 60% of participants had a body mass index (BMI) percentile ≥95th percentile. Exclusion criteria were: inability to complete MRI evaluation (claustrophobia, MRI-contraindicated metal implants, or body circumference greater than the imaging chamber), established diagnosis of chronic liver disease, chronic use of medications known to cause hepatic steatosis or raise liver chemistry, chronic diseases that may have secondary effects on the liver, substance abuse, and pregnancy. The study was approved by the Institutional Review Board of the UC San Diego. Written informed consent from the parent or legal guardian and assent from the participants were obtained.

Clinical and Laboratory Evaluation

Participants made a fasting visit to the Altman Clinical and Translational Research Institute at UC San Diego. Weight and height were measured to the nearest 0.1 kg and 0.1 cm, respectively, with the subjects standing, wearing light clothing without shoes. BMI was calculated as weight (kilograms) divided by height squared (meters). Blood was collected after a 12-hour fast for fatty acid analysis (see below) and the following routine laboratory measures: complete blood count, alanine aminotransferase (ALT), aspartate aminotransferase (AST), gamma-glutamyl transpeptidase (GGT), glucose, insulin, glycosylated hemoglobin, total cholesterol, HDL-cholesterol, and triglycerides.

Dietary Evaluation

All participants completed 3 24-hour dietary recalls. Recalls were performed on nonconsecutive days, capturing 2 weekdays and 1 weekend day within a 2-week period. Recalls were done by a multiple-pass interview approach performed by a registered dietician using Nutrition Data System for Research (NDSR; University of Minnesota Nutrition Coordinating Center, Minneapolis, MN). Assessment of dairy fat intake in grams per day was done using published NDSR methodology (10). Data from the 3 recalls were averaged to yield per day estimates. Total calories per day were calculated and were considered unreliable and excluded if reported intake was less than 50% of the recommended daily caloric intake for age (11,12). We also calculated daily intake (grams per day) of protein, fat, carbohydrates, and sugars.

Fatty Acid Analysis

In addition to dietary recall data, plasma fatty acid profiles provide an objective measure of recent dietary intake. We measured OCFA and BCFA associated with dairy fat including C15:0, C17:0, iso-C17:0, and anteiso-C17:0 (13,14). We also measured 20 other fatty acids selected as the most abundant fatty acids in human plasma (15). A complete list is shown in Table, Supplemental Digital Content, https://links.lww.com/MPG/C162. Plasma total fatty acids were extracted using a Folch-based methanol/chloroform/saline extraction at a ratio of 1 : 2 : 1 with inclusion of d31 C16:0, d3 C14:0, d3 C15:0, d3 C17:0, d3 C18:0, d2 C18:1n9, d3 C22:0, and d4 C24:0 as internal standards (16). Briefly, 250 μL MeOH, 500 μL CHCl3, 250 μL saline, and the fatty acid isotope internal standards were added to 100 μL plasma. This was vortexed for 10 minutes followed by centrifugation at 10,000 g for 5 minutes at 4 °C. The lower chloroform phase was dried under nitrogen and then derivatized to form fatty acid methyl esters (FAMEs) via addition of 500 μL 2% H2SO4 in MeOH and incubation at 50 °C for 2 hours. FAMEs were extracted via addition of 100 μL saturated salt solution and 500 μL hexane and were then analyzed using a Select FAME column (100 m × 0.25 mm i.d.) installed in an Agilent 7890A GC interfaced with an Agilent 5975C MS using the following temperature program: 80 °C initial, increase by 20 °C/minutes to 170 °C, increase by 1 °C/minutes to 204 °C, then 20 °C/minutes to 250 °C and hold for 10 minutes. Fatty acids were quantified from deuterated internal standards or external standard curves.

Magnetic Resonance Imaging Evaluation

MRI examinations were performed on a 3 Tesla scanner (GE Discovery 750, General Electric Healthcare, Wakeshaw, WI) using a previously described advanced magnitude-based confounder-corrected chemical-shift-encoded acquisition and reconstruction technique to estimate proton density fat fraction (PDFF) (17–19). T1 weighting was minimized by using a gradient-recalled echo sequence with a low (10°) flip angle relative to a repetition time (TR) of ≥150 mseconds. Six gradient-recalled echoes were collected at successive nominally out-of-phase and in-phase echo times to allow correction for T2∗ signal decay (20,21). Computer-generated parametric PDFF maps were calculated using least-squares fitting analysis based on a 6 fat-peak spectral model to correct for inter-peak spectral interference (22). For each PDFF parametric map, a 1-cm radius circular region of interest was manually placed in each of the 9 Couinaud liver segments, and a composite mean PDFF value was calculated (8).

Data Analysis

Demographics and other study population characteristics were described. Categorical data were summarized as percentages and counts. One-way analysis of variance was used to calculate mean, median, range, and standard deviation of all continuous data. Correlation analysis was performed using Pearson correlation for parametric data and Spearman rank order correlation for nonparametric data. After bivariate analysis, stepwise multiple linear regression was used to determine the relationship between dairy fat intake (predictor variable) and hepatic MRI-PDFF (outcome variable) while controlling for age, sex, and BMI z-score as potential confounders. This relationship was also evaluated by sensitivity analyses to further control for other dietary factors (total calorie intake, total fat intake, total carbohydrate intake, total sugar intake). To support these analyses, a minimum sample size of 200 participants was projected. In addition, multivariate stepwise linear regression analyses were used to examine fatty acid predictors of MRI-PDFF while controlling for confounders. Variables that were not normally distributed were transformed logarithmically for regression analyses (liver MRI-PDFF and plasma fatty acids). Statistical analysis was performed using IBM SPSS Statistics for Windows, Version 24.0 (IBM Corp, Armonk, NY) and R 3.5.0 (R foundation for statistical computing, Vienna, Austria,). A P value <0.05 was considered significant for all inference testing.


Study Cohort

Flow for screening, exclusion, and inclusion is shown in Figure, Supplemental Digital Content 2, https://links.lww.com/MPG/C163. Characteristics of the study population are shown in Table 1. We studied 237 children with a mean age of 12.9 ± 2.6 years. There were more male than female participants (62% vs 38%). Obesity was present in 71% of participants. Mean ALT (SD) was 35 (49) U/L.

TABLE 1 - Characteristics of study population, n = 237
Characteristic Value
 Age, years, mean (SD) 12.9 (2.6)
Sex, n (%)
 Male 146 (62)
 Female 91 (38)
Race, n (%)
 American Indian 15 (6)
 Black 18 (8)
 White 104 (44)
 Other 100 (42)
Ethnicity, n (%)
 Hispanic 178 (75)
 Non-Hispanic 59 (25)
 Weight, kg, mean (SD) 71.0 (24.7)
 Height, cm, mean (SD) 157.3 (13.8)
 BMI, kg/m2, mean (SD) 28.0 (6.9)
 BMI percentile, mean (SD) 89.6 (20.8)
 BMI z score, mean (SD) 1.7 (1.0)
Nutrition (per day)
 Total calories, mean (SD) 1704 (377)
 Fat, g, mean (SD) 62 (23)
 Protein, g, mean (SD) 72 (21)
 Carbohydrates, g, mean (SD) 214 (53)
 Sugars, g, mean (SD) 85 (42)
 Dairy fat, g, mean (SD) 12.1 (7.9)
 ALT, U/L, mean (SD) 35 (49)
 AST U/L, mean (SD) 31 (27)
 GGT U/L, mean (SD) 25 (28)
 Glucose mg/dL, mean (SD) 87.0 (12)
 HbA1c, %, mean (SD) 5.4 (0.7)
 Insulin, μIU/mL, mean (SD) 27 (23)
 Triglycerides, mg/dL, mean (SD) 114 (69)
 HDL-cholesterol, mg/dL, mean (SD) 45 (11)
 LDL-cholesterol, mg/dL, mean (SD) 88 (25)
 Liver MRI-PDFF, %, mean (SD) 7.9 (8.1)
ALT = alanine aminotransferase; AST = aspartate aminotransferase; BMI = body mass index; GGT = gamma-glutamyl transferase; HDL = high-density lipoprotein; LDL = low-density lipoprotein; MRI-PDFF = magnetic resonance imaging proton density fat fraction; SD = standard deviation; U/L = units/liter.

Dairy Fat Intake

The median dairy fat intake in study participants was 10.6 g/day (range 0.0–44.5 g). The distribution of dairy fat consumed per day is shown in Figure, Supplemental Digital Content 3, https://links.lww.com/MPG/C164. Dairy fat intake was not significantly associated with age (P = 0.148) or sex (P = 0.911) (Figure, Supplemental Digital Content 4, https://links.lww.com/MPG/C165). Dairy fat intake (g/day) was significantly negatively correlated with BMI z-score (r = −0.20, P = 0.002). Dairy fat intake was also significantly negatively correlated with ALT (r = −0.18, P = 0.007), but not HDL (r = 0.07, P = 0.286), LDL (r = −0.07, P = 0.275), or triglycerides (r = −0.09, P = 0.151).

Dairy Fat Intake and Liver Fat

The median MRI liver PDFF was 4.5% (range 0.9%–45.1%). There was a significant negative correlation between grams of dairy fat consumed per day and liver MRI-PDFF (r = −0.162; P = 0.012) (Fig. 1). A multiple linear regression model was calculated to predict liver MRI-PDFF based upon dairy fat intake after controlling for age, sex, ethnicity, and BMI z-score. The model was significant with R2 = 0.297; P < 0.001 (Table 2). The inverse relationship between the amount of dairy fat consumed and liver MRI-PDFF was significant after controlling for the other elements of the model (β −0.152, t −2.675, P = 0.008) (Table 2). In this model, an additional 1 cup of whole milk (7.93 g of dairy fat) would be associated with liver MRI PDFF that is 6.5% points lower. In sensitivity analyses, further controlling for total calories, total fat, total carbohydrate, or total sugar intake attenuated but did not eliminate the relationship between dairy fat intake and liver MRI-PDFF (with total calorie intake, β −0.133, t −2.145, P = 0.033; with total fat intake, β −0.127, t −2.138, P = 0.034; with total carbohydrate intake, β −0.141, t −2.372, P = 0.019; with total sugar intake, β −0.148, t −2.602, P = 0.012) (Fig. 1).

Distribution of daily dairy fat intake by liver parameters. Panel A shows ALT; dairy fat was significantly negatively correlated with ALT (r = −0.18, P = 0.007). Panel B shows liver MRI-PDFF; dairy fat was significantly negatively correlated with liver MRI-PDFF (r = − 0.162; P = 0.012). Panel C shows liver MRI-PDFF separated by children above and below 5% liver MRI-PDFF. Children with MRI liver-PDFF ≥5% had significantly lower dairy fat intake than children with liver MRI-PDFF <5% (mean 9.1 vs 14.2 g/day; P < 0.001). ALT = alanine aminotransferase; MRI-PDFF = magnetic resonance imaging proton density fat fraction.
TABLE 2 - Multivariate regression model of dairy fat intake as a predictor of liver magnetic resonance imaging proton density fat fraction
β t P
Dairy fat, g/day −0.152 −2.675 0.008
Age, years −0.372 −4.006 <0.001
Sex (male) 0.125 2.111 0.036
Ethnicity (Hispanic) 0.401 4.623 <0.001
BMI z-score 0.451 5.200 <0.001
Intercept 2.040 5.695 <0.001
Variables that were not normally distributed were transformed logarithmically for regression analyses (dairy fat and liver MRI-PDFF). Model R2 = 0.297, P < 0.001.BMI = body mass index; MRI-PDFF = magnetic resonance imaging proton density fat fraction.

Plasma Fatty Acids and Liver Fat

The correlation of plasma fatty acids with liver MRI-PDFF was evaluated. Of the 29 plasma fatty acids measured, 20 were significantly correlated with liver MRI-PDFF in univariate analysis (Table 3). Of these, only 3, C15:0, C17:0, and iso-C17:0, had inverse correlations with liver MRI-PDFF. A stepwise multivariable log linear model of plasma fatty acids adjusted for age, sex, ethnicity, and BMI z-score was significantly predictive of liver MRI-PDFF (Table 4, R2 = 0.481; P < 0.001). The covariates sex (male vs female individuals) and BMI z-score were independent predictors of liver MRI-PDFF. Only 5 plasma fatty acids remained independently associated with liver MRI-PDFF. Negative predictors of liver MRI-PDFF were plasma levels of C15:0, iso-C17:0, and C20:4n6. Positive predictors of liver MRI-PDFF were plasma levels of C16:1 and C18:0.

TABLE 3 - Univariate linear regression of plasma fatty acid concentrations as predictors of liver magnetic resonance imaging proton density fat fraction
Fatty acid Name β t P
Significant negative correlations
 Iso-C17:0 Iso-heptadecanoic acid −0.275 −5.838 <0.001
 C15:0 Pentadecanoic acid −0.115 −2.354 0.006
 C17:0 Heptadecanoic acid −0.107 −2.113 0.017
Significant positive correlations
 C16:1 Palmitoleic acid 0.345 7.502 <0.001
 C20:3, n-6 Dihomo-γ-linoleic acid 0.305 6.528 <0.001
 C16:0 Palmitic acid 0.281 5.972 <0.001
 C18:0 Stearic acid 0.266 5.617 <0.001
 C18:1 Oleic acid 0.261 5.505 <0.001
 C20:1n9 11-Eicosenoic 0.253 5.307 <0.001
 C14:0 Myristic acid 0.218 4.530 <0.001
 C20:2 11,14-Eicosadienoic acid 0.181 3.738 <0.001
 C17:1 Heptadecenoate 0.179 3.685 <0.001
 C18:3, n-3 α-Linolenic acid 0.160 3.293 <0.001
 C20:0 Arachidic acid 0.156 3.202 0.001
 C14:1 Myristoleic 0.147 3.012 0.003
 C18:2 Linoleic acid 0.144 2.954 0.003
 C18:3, n-6 γ-Linolenic acid 0.136 2.779 <0.001
 C12:0 Lauric acid 0.116 2.360 0.019
 C20:5, n-3 Eicosapentaenoic acid 0.102 2.313 0.019
No significant correlations
 iso-C16:0 Iso-palmitic acid −0.000 −0.003 0.998
 Anteiso- C17:0 Anteiso-heptadecanoic acid 0.018 0.372 0.512
 Iso-C18:0 Iso-octadecanoate −0.026 −0.527 0.599
 C20:3, n-3 11,14,17-Eicosatrienoic acid 0.068 1.389 0.166
 C20:4, n-6 Arachidonic acid −0.042 −0.849 0.397
 C22:6, n-3 Docosahexaenoic acid 0.059 1.202 0.230
 C22:1n9 Erucic 0.0027 0.544 0.587
 C22:2 13,16-Docosadienoic acid −0.035 −0.714 0.476
 C22:0 Behenic acid 0.017 0.345 0.730
 C24:1n9 Nervonic acid 0.027 0.544 0.587

TABLE 4 - Multivariate model of plasma fatty acids as predictors of liver magnetic resonance imaging proton density fat fraction
Fatty acids β t P
C15:0 −0.247 −2.00 0.041
C16:1 0.322 5.90 <0.001
Iso-C17:0 −0.234 −2.63 0.009
C18:0 0.426 4.68 <0.001
C20:4n6 −0.281 −4.60 <0.001
 Age, y −0.242 −2.834 0.007
 Sex (male) 0.139 2.852 0.005
 Ethnicity (Hispanic) 0.277 3.743 <0.001
 BMI z-score 0.312 5.891 <0.001
Model R2 = 0.481; P < 0.001. BMI = body mass index; C15:0 = pentadecanoic acid; C16:1 = palmitoleic acid; iso-C17:0 = iso-heptadecanoic acid, C18:0 = stearic acid, C20:4n6 = arachidonic acid methyl ester; MRI-PDFF = magnetic resonance imaging proton density fat fraction.


We evaluated a large community-based sample of children at risk for NAFLD and assessed the relationship between dietary dairy fat intake and hepatic steatosis. We also evaluated plasma fatty acids in relationship to hepatic steatosis. We found that the average daily intake of dairy fat was significantly inversely correlated with liver MRI-PDFF. In plasma, the concentrations of the OCFA C15:0 and the BCFA iso-C17:0 were significant negative predictors of liver MRI-PDFF. In addition to C15:0 and iso-C17:0, liver MRI-PDFF was negatively associated with arachidonic acid (C20:4n6) and was positively associated with stearic acid (C18:0) and palmitoleic acid (C16:1).

We found that dairy fat intake was inversely correlated with liver MRI-PDFF in children. We also observed that dairy fat intake was negatively associated with BMI z-score, which is consistent with studies that have shown lower weight gain in children drinking whole milk (2). Importantly, the negative association between dairy fat intake and liver fat was independent of BMI z-score. There have been no prior pediatric studies of the relationship between dairy fat intake and liver fat; however, our data build upon 2 prior adult studies of this relationship. Kratz et al (23) performed an observational study of 17 adults with NAFLD and 15 controls and found a significant negative correlation between dairy fat intake and computed tomography measured liver-spleen ratio, a surrogate measure for liver fat. In the second, an interventional crossover study, Dugan et al (24) evaluated 37 adults who consumed nondairy control foods in 1 phase and 3 servings of dairy per day in the other phase. This study reported significantly lower AST and ALT at the end of the dairy consumption period compared with the control foods phase; however, there were no measurements of hepatic fat performed. In our study, we also found that higher dairy fat intake was associated with lower ALT values.

The OCFA, C15:0 was negatively associated with hepatic steatosis. C15:0 is produced by microbial fatty acid oxidation and de novo lipogenesis in the bovine rumen and the primary source of this for humans is the consumption of dairy fat (25,26). In adults, C15:0 is an established circulating marker of dairy fat intake (27–29), and interventional studies show that circulating levels of C15:0 can be modified by changing the dietary intake of dairy fat (30–32). There are no prior pediatric data on the relationship between the plasma measures of C15:0 and liver fat. Two studies in adults have evaluated C15:0 concentration levels in subjects with NAFLD. Yoo et al evaluated plasma fatty acid profiles in 106 adults with NAFLD and found that plasma C15:0 was significantly lower in patients with higher NAFLD activity scores. Similarly, Kratz et al reported that plasma C15:0 was significantly lower in 17 adults with NAFLD than in 15 controls without NAFLD. In much larger, longitudinal studies, higher plasma levels of C15:0 were associated with lower incidence of both cardiovascular disease and type 2 diabetes (33,34). In a mouse model, supplementation with C15:0 decreased the severity of NAFLD (35). In humans, whether increasing the intake of C15:0 could potentially influence NAFLD is not yet known.

We also demonstrated that lower levels of iso-C17:0 were associated with higher amounts of hepatic steatosis. Prior studies of fatty acid composition in subjects with NAFLD have not evaluated this isomer. Iso-C17:0 is a branched chain fatty acid made by rumen microbiota and found primarily in dairy products (36). In addition to dietary sources, iso-C17:0 can be synthesized de novo from branched chain amino acids (BCAA) (6). In adults, hepatic steatosis was associated with elevated plasma BCAA (6). One reason for simultaneous high BCAA and low BCFA can be a disturbance in the conversion from BCAA to BCFA. Thus, our finding that plasma iso-C17:0 concentrations were significantly lower in children with NAFLD may be because of both lower dietary intake of iso-C17:0 and reduced metabolism of BCAA to BCFA.

In addition to perceptions regarding dairy fat, there are widely held notions regarding the healthfulness of other dietary fats. Dietary advice has progressed from limiting all dietary fat to preferencing “good” fats, such as monounsaturated fats and avoiding “bad” fats, such as saturated fats. Following such dietary recommendations is, however, difficult as we do not eat or drink foods in isolation; rather diets are a complex mixture of many foods. In turn, the fatty acid composition of foods can have a wide range of effects. In addition to C15:0 and iso-C17:0, there were 3 other fatty acids associated with liver MRI-PDFF. Stearic acid, C18:0, was the most strongly positively associated with hepatic steatosis. Stearic acid is a saturated fat that is highly prevalent in many foods, especially beef and pork (37). Stearic acid is also, however, present in dairy fat as well as plant-based oils. Palmitoleic acid, C16:1, was also strongly positively associated with hepatic steatosis. Palmitoleic acid is a monounsaturated fat also found in many foods. The top sources are beef, pork, avocado, butter, salmon, and olive oil (38). Conversely, arachidonic acid was inversely associated with liver MRI-PDDF. Arachadonic acid is an omega-6 polyunsaturated fatty acid found in many foods including eggs, fish, and butter (37). It is difficult to disentangle the source of each of these fatty acids in the individual diet for each participant, and thus it remains challenging to develop fully optimized dietary recommendations to prevent or treat NAFLD.

Strengths of this study include the sample size, detailed phenotyping, and the measurement of liver fat by liver MRI-PDFF. We used a validated approach to assess quantitative dairy fat intake. There are known limitations to the self-report of dietary intake, such as recall bias. Dietary choice may be confounded by many factors and the analyses did not control for these. In addition, dairy fat can be consumed in numerous forms, such as milk, cheese, yogurt, butter, or as part of mixed foods. A larger study would be required to test for the effect of specific foods. To mitigate the limitations of self-reporting, we measured plasma fatty acids known to be associated with dairy fat intake by mass spectroscopy (39). Importantly, the inverse relationship between the dairy-associated fatty acids C15:0 and iso-C17:0 was independent of other fatty acids as well as sex and BMI z-score. This study was cross-sectional and observational, and thus did not evaluate causality. The findings should be considered as hypothesis-generating.


In conclusion, we found that dairy fat intake as well as the plasma levels of the OCFA, C15:0, and the BCFA, iso-heptadecanoic acid (iso-C17:0), were inversely correlated with hepatic steatosis in children. These findings have implications for dietary guidelines and future biomarker studies. Clinical trials of whole fat dairy would be needed to know whether consuming higher amounts of dairy fat decreases liver fat in children. Such studies should consider to what extent changes are mediated by increases in circulating OCFA and BCFA.


1. Gidding SS, Dennison BA, Birch LL, et al. Dietary Recommendations for Children and Adolescents A Guide for Practitioners Consensus Statement From the American Heart Association Endorsed by the American Academy of Pediatrics. Circulation 2005; 112:2061–2075.
2. Beck AL, Heyman M, Chao C, et al. Full fat milk consumption protects against severe childhood obesity in Latinos. Prev Med Rep 2017; 8:1–5.
3. Imamura F, Fretts A, Marklund M, et al. Fatty acid biomarkers of dairy fat consumption and incidence of type 2 diabetes: a pooled analysis of prospective cohort studies. PLoS Med 2018; 15:e1002670.
4. Kummer K, Jensen PN, Kratz M, et al. Full-fat dairy food intake is associated with a lower risk of incident diabetes among American Indians with low total dairy food intake. J Nutr 2019; 149:1238–1244.
5. Yakoob MY, Shi P, Willett WC, et al. Circulating biomarkers of dairy fat and risk of incident diabetes mellitus among men and women in the United States in two large prospective cohorts. Circulation 2016; 133:1645–1654.
6. Wallace M, Green CR, Roberts LS, et al. Enzyme promiscuity drives branched-chain fatty acid synthesis in adipose tissues. Nat Chem Biol 2018; 14:1021–1031.
7. Schwimmer JB, Deutsch R, Kahen T, et al. Prevalence of fatty liver in children and adolescents. Pediatrics 2006; 118:1388–1393.
8. Yu E, Golshan S, Harlow K, et al. Prevalence of nonalcoholic fatty liver disease in children with obesity. J Pediatr 2019; 207:64–70.
9. Spear BA, Barlow SE, Ervin C, et al. Recommendations for treatment of child and adolescent overweight and obesity. Pediatrics 2007; 120:S254–S288.
10. Dugan CE, Barona J, Fernandez ML. Increased dairy consumption differentially improves metabolic syndrome markers in male and female adults. Metab Syndr Relat Disord 2014; 12:62–69.
11. Evans EW, Jacques PF, Dallal GE, et al. The role of eating frequency on total energy intake and diet quality in a low-income, racially diverse sample of schoolchildren. Public Health Nutr 2015; 18:474–481.
12. Phillips KM, Patterson KY, Rasor AS, et al. Quality-control materials in the USDA National Food and Nutrient Analysis Program (NFNAP). Anal Bioanal Chem 2006; 384:1341–1355.
13. Yuzyuk T, Lozier B, Schwarz EL, et al. Intra-individual variability of long-chain fatty acids (C12-C24) in plasma and red blood cells. Prostaglandins Leukot Essent Fatty Acids 2018; 135:30–38.
14. Albani V, Celis-Morales C, O’Donovan CB, et al. Within-person reproducibility and sensitivity to dietary change of C15:0 and C17:0 levels in dried blood spots: data from the European Food4Me Study. Mol Nutr Food Res 2017; 61: 10.1002/mnfr.201700142.
15. Abdelmagid SA, Clarke SE, Nielsen DE, et al. Comprehensive profiling of plasma fatty acid concentrations in young healthy Canadian adults. PLoS One 2015; 10:e0116195(published correction appears in PLoS One. 2015;10(5);e0128167).
16. Metallo CM, Gameiro PA, Bell EL, et al. Reductive glutamine metabolism by IDH1 mediates lipogenesis under hypoxia. Nature 2011; 481:380–384.
17. Schwimmer JB, Middleton MS, Behling C, et al. Magnetic resonance imaging and liver histology as biomarkers of hepatic steatosis in children with nonalcoholic fatty liver disease. Hepatology 2015; 61:1887–1895.
18. Bydder M, Yokoo T, Hamilton G, et al. Relaxation effects in the quantification of fat using gradient echo imaging. Magn Reson Imaging 2008; 26:347–359.
19. Liu CY, McKenzie CA, Yu H, et al. Fat quantification with IDEAL gradient echo imaging: correction of bias from T(1) and noise. Magn Reson Med 2007; 58:354–364.
20. Yu H, Shimakawa A, McKenzie CA, et al. Multiecho water-fat separation and simultaneous R2∗ estimation with multifrequency fat spectrum modeling. Magn Reson Med 2008; 60:1122–1134.
21. Yu H, McKenzie CA, Shimakawa A, et al. Multiecho reconstruction for simultaneous water-fat decomposition and T2∗ estimation. J Magn Reson Imaging 2007; 26:1153–1161.
22. Hamilton G, Yokoo T, Bydder M, et al. In-vivo characterization of the liver fat 1H MR spectrum. NMR Biomed 2011; 784–790.
23. Kratz M, Baars T, Guyenet S. The relationship between high-fat dairy consumption and obesity, cardiovascular, and metabolic disease. Eur J Nutr 2013; 52:1–24.
24. Dugan CE, Aguilar D, Park YK, et al. 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:255–261.
25. Crowe FL, Key TJ, Appleby PN, et al. Dietary fat intake and risk of prostate cancer in the European Prospective Investigation into Cancer and Nutrition. Am J Clin Nutr 2008; 87:1405–1413.
26. Smedman A, Gustafsson IB, Berglund L, et al. Pentadecanoic acid in serum as a marker for intake of milk fat: relations between intake of milk fat and metabolic risk factors. Am J Clin Nutr 1999; 69:22–29.
27. Risérus U, Marklund M. Milk fat biomarkers and cardiometabolic disease. Curr Opin Lipidol 2017; 28:46–51.
28. Wolk A, Vessby B, Ljung H, et al. Evaluation of a biological marker of dairy fat intake. Am J Clin Nutr 1998; 68:291–295.
29. Wolk A, Furuheim M, Vessby B. Fatty acid composition of adipose tissue and serum lipids are valid biological markers of dairy fat intake in men. J Nutr 2001; 131:828–833.
30. Vessby B, Uusitupa M, Hermansen K, et al. KANWU StudySubstituting dietary saturated for monounsaturated fat impairs insulin sensitivity in healthy men and women: the KANWU study. Diabetologica 2001; 44:312–319.
31. Brevik A, Veierød MB, Drevon CA, et al. Evaluation of the odd fatty acids 15:0 and 17:0 in serum and adipose tissue as markers of intake of milk and dairy fat. Eur J Clin Nutr 2005; 59:1417–1422.
32. Wennersberg MH, Smedman A, Turpeinen AM, et al. Dairy products and metabolic effects in overweight men and women: results from a 6-mo intervention study. Am J Clin Nutr 2009; 90:960–996.
33. Otto MC, Lemaitre RN, Song X, et al. Serial measures of circulating biomarkers of dairy fat and total and cause-specific mortality in older adults: the Cardiovascular Health Study. Am J Clin Nutr 2018; 108:476–484.
34. Mozaffarian D, Otto MC, Lemaitre RN, et al. trans-Palmitoleic acid, other dairy fat biomarkers, and incident diabetes: the Multi-Ethnic Study of Atherosclerosis (MESA). Am J Clin Nutr 2013; 97:854–861.
35. Yoo W, Gjuka D, Stevenson HL, et al. Fatty acids in non-alcoholic steatohepatitis: focus on pentadecanoic acid. PLoS One 2017; 12:e0189965.
36. Vlaeminck B, Fievez V, Tamminga S, et al. Milk odd- and branched-chain fatty acids in relation to the rumen fermentation pattern. J Dairy Sci 2006; 89:3954–3964.
37. Wholefoodcatalog.info. 2020. Foods High In Arachidonic Acid | Whole Food Catalog. https://wholefoodcatalog.info/nutrient/arachidonic_acid/foods/high/. Accessed July 18, 2020.
38. Wholefoodcatalog.info. 2020. Palmitoleic Acid Content Of Foods | Whole Food Catalog. https://wholefoodcatalog.info/nutrient/palmitoleic_acid/foods/. Accessed July 18, 2020.
39. Foster E, Bradley J. Methodological considerations and future insights for 24-hour dietary recall assessment in children. Nutr Res 2018; 51:1–11.

milk; nonalcoholic steatohepatitis; nutrition; obesity

Supplemental Digital Content

Copyright © 2020 by European Society for Pediatric Gastroenterology, Hepatology, and Nutrition and North American Society for Pediatric Gastroenterology, Hepatology, and Nutrition