Dietary fatty acids and risk of non-alcoholic steatohepatitis: A national study in the United States

Background: Non-alcoholic steatohepatitis (NASH), the early invertible stage of non-alcoholic fatty liver disease, has become a public health challenge due to the great burden and lack of effective treatment. Dietary nutrients are one of the modifiable factors to prevent and slow down disease progression. However, evidence linking dietary fatty acids intake and risk of NASH is lacking.

Objectives: This study aimed to examine the association between dietary total saturated fatty acids (SFAs), monounsaturated fatty acids (MUFAs), polyunsaturated fatty acids (PUFAs), their subtypes, the ratio of unsaturated (UFAs) to SFAs, and the risk of NASH among a nationwide population in the United States.

Methods: This cross-sectional study was conducted among 4,161 adults in the national health and nutrition examination survey in 2017–2018 cycle. Moreover, NASH was defined by transient elastography. Dietary fatty acids were assessed using a validated 24-h food recall method. Logistic regression models were used to estimate odds ratios (ORs) and 95% confidence intervals (95% CIs).

Results: A total of 2,089 (50.2%) participants with NASH were identified. Compared with participants in the bottom tercile of dietary intakes of total PUFAs, those in the highest tercile had lower risk of NASH, with an adjusted OR of 0.67 (95% CI: 0.46–0.97). Similar associations were found between the subtype of PUFA 18:3 and NASH, while the fully adjusted OR in the highest tercile was 0.67 (95% CI: 0.47–0.96). Interactions of dietary PUFAs and body mass index (BMI) could be found influencing NASH risk. Stronger associations of dietary total PUFAs intakes with NASH risk were found in obese participants (OR, 95% CI: 0.41, 0.22–0.75) than in the non-obese participants (OR, 95% CI: 1.00, 0.70–1.43; p-interaction = 0.006). Similar effects on risk of NASH were also observed between BMI and dietary intakes of PUFA 18:3. However, no significant associations were observed between NASH risk and dietary total SFAs, MUFAs, their subtypes as well as the ratio of UFAs to SFAs.

Conclusion: Dietary intakes of total PUFAs, as well as its subtype of PUFA 18:3, were inversely associated with risk of NASH. The further large prospective studies need to be conducted to confirm the findings of this study.

Introduction

Non-alcoholic fatty liver disease (NAFLD) has emerged as one of the most common causes of chronic liver disease, as well as the most quickly growing promoters to liver morbidity and mortality in the United States (1). The prevalence of NAFLD has increased rapidly over the past three decades (1, 2) and the number of people suffering from NAFLD are expected to be more than 100 million by 2030 (3). Usually, parallels to the prevalence of obesity, type 2 diabetes and cardiovascular disease, NAFLD decreases life expectancy and increases requirements of liver transplantation (4, 5). Moreover, NAFLD progresses from steatosis, non-alcoholic steatohepatitis (NASH), to liver fibrosis and even ultimately to hepatocellular carcinoma (4), while approximately 20–30% of individuals may progress to NASH (6). It is estimated that NASH prevalence will increase up to 56% by 2030 in the United States at a greater rate than the other European and Asian countries (3). Due to the persistent cellular damage and excess fat deposition, patients with NASH are more likely to progress to the irreversible advanced fibrosis (4). Given the great burden as well as the lack of effective treatment, it is important to prevent and slow down disease progression at the early stage. The dietary nutrients are considered to be one of the effective and modifiable factors (7).

Fatty acids are composed of saturated fatty acids (SFAs) and unsaturated fatty acids (UFAs), the latter includes monounsaturated fatty acids (MUFAs) and polyunsaturated fatty acids (PUFAs). They can also be subdivided into a number of subtypes according to the position and numbers of carbon atoms and double bonds. Fatty acids play an important and indispensable role in human diet and body metabolism, which not only are sources of essential fatty acids and precursors of bioactive substances but also involve in bio-membrane structure formation, signal transmission, and lipid transportation. Liver is the main metabolic organ of fatty acids. Improper dietary fatty acids consumption or de novo lipogenesis beyond liver capacity has the consequences of abnormal cellular lipid composition and toxic lipid accumulation, leading to organelle dysfunction, cellular damage, inflammation, and occurrence of diseases (8). However, the results from epidemiological studies linking the relationships between dietary fatty acids and NAFLD are sparse and contradictory. For example, a systematic review reported that patients with NAFLD had disturbed fatty acids metabolisms compared to healthy controls (6). Evidences from the previous studies demonstrated that people with NASH possessed elevated total SFAs, while total MUFAs and PUFAs may be protective (9). However, a non-linear association between total PUFAs intakes and NAFLD was found among Chinese Han adults, the total PUFAs was positively associated with risk of NAFLD at certain dose of total PUFAs consumption (10). The reasons for the discrepancies could be ascribed to differences of sample sizes, population, and different detect methods of NAFLD.

Liver biopsy, the gold standard of NAFLD assessments, is unsuitable for population screening due to its invasion. Other methods of detection such as hepatic steatosis index (11), fatty liver index (12), NAFLD liver fat score (13), SteatoTest (14), which are calculated by individual conditions and blood biomarkers, have several shortages such as limited sensitivity, specificity, popularization or expensive price. A non-invasive and effective approach to distinguish NAFLD as well as its stages is urgently needed. Compared with liver biopsy, transient elastography is a non-invasive, convenient and fast method to assess NAFLD from NASH to liver fibrosis, but few studies apply it to diagnosis of NAFLD in a large population.

Furthermore, different fatty acids subtypes have different effects on liver health. For example, as for hepatic fatty acids compositions, a decrease in the ratio of SFA 18:0 to SFA 16:0 was associated with the steatosis score and insulin resistance, while higher ratio of MUFA 16:1 to SFA 16:0 was associated with lobular inflammation and hepatocellular ballooning in patients with NASH (15). However, none of the existing studies have investigated the association between different subtypes of dietary fatty acids and the risk of NAFLD. Therefore, in this study, we aimed to examine the association between dietary total SFAs, MUFAs, PUFAs, their subtypes, the ratio of UFAs to SFAs, and the risk of NASH (the early invertible stage of NAFLD), which was determined by liver ultrasound transient elastography among a nationwide population in the United States.

Materials and methods

Study population

The analysis of this study was conducted based on data from national health and nutrition examination survey (NHANES), which was carried out by the Center for Disease Control and Prevention of the United States. Detailed information could be found from the official website¹. In this study, we included participants in the 2017–2018 NHANES cycle (n = 9,254). Participants were excluded if they were underage (<18 years, n = 3,398), had unavailable data on transient elastography (n = 737) or dietary assessment (n = 445), were hepatitis C virus (n = 44) or hepatitis B virus (n = 18) infected, or had heavy alcohol consumption (>30 g/day for men and >20 g/day for women, n = 451). Finally, 4,161 participants were included (Figure 1). Due to the discrepancies between adults and minors in many aspects, we only compare the characteristics of adults (≥18 years) in the included and the excluded participants. The included participants did not differ by most of the basic characteristics from adults excluded from this analysis (n = 1,695, Supplementary Table 1). The institutional review board approval of the National Center for Health Statistics and informed consent was obtained from all participants before data collection.

FIGURE 1

Figure 1. Flowchart of selection of participants from national health and nutrition examination survey (NHANES). Abbreviations: NASH, non-alcoholic steatohepatitis.

Assessment and definition of non-alcoholic steatohepatitis

The liver ultrasound transient elastography was first used to provide objective measures for hepatic steatosis in the NHANES Mobile Examination Center (MEC) in the 2017–2018 cycle (16). The participants who were aged ≥12 years, were able to lie down, were not pregnant, had no implanted electronic medical device, and had no lesions where measurements would be taken were eligible to test. Using the FibroScan model 502 V2 Touch equipped with a medium or extra-large probe, controlled attenuation parameter (CAP) was measured and recorded as the indicator for hepatic steatosis according to the fatness in the liver. The elastography exam was performed by well-trained NHANES health technicians according to the manufacture guidelines. The inter-rater reliability between health technician and reference examiners was 0.94 (mean differences 4.5 ± 19.8 dB/m), and variances within and between the device machines and the probes over time were under control (intra-machine coefficient of variation was 1.2–3.2%; inter-machine intra class correlation was 2–22%) for CAP. High accuracy of the CAP measurements for the detection of steatosis compared to biopsy has been reported in the previous studies (17–19).

In this study, those CAP scores equal or greater than 263 dB/m (had a 96% positive predictive value) were defined as cases of NASH (20), others as controls (non-NASH). The cut-off point was determined before statistical analysis.

Assessment of dietary fatty acids

Assessment of dietary intakes was undertaken by trained interviewers using a validated 24-h food recall method in the MEC, and was repeated by telephone 3–10 days later. Values of dietary intakes in these two discrete days were averaged to represent their dietary status, while values of the first time were used for participants with a lack of dietary data of the second time (12.2% in the current analysis). The intakes of total energy, SFAs, MUFAs, PUFAs, and their subtypes were estimated according to the United States Department of Agriculture’s Food and Nutrient Database for Dietary Studies. The healthy eating index-2015 (HEI-2015) was calculated according to the MyPyramid Equivalents Database 2.0 for United States Department of Agriculture Survey Foods to reflect the quality of diet comprehensively (Supplementary Table 2). A higher score between 0–100 indicated a better dietary quality. The ratio of UFAs to SFAs was calculated as (PUFAs + MUFAs)/SFAs, and therefore divided into three levels according to the scoring standards of HEI-2015.

Covariates collection

Sociodemographic information including age, sex (men and women), ethnicity (non-Hispanic white, non-Hispanic black, Mexican American, and others), marital status (married/living with partner, separated/divorced/widowed, and never married), education levels (less than high school, high school or equivalent, and college or above), and family income-to-poverty ratio (<1.3, 1.3–3.5, and >3.5) were collected. Anthropometric measurements including height, weight, and waist circumference were measured. Body mass index (BMI) was calculated as weight (kg)/height squared (m²). Never smokers were participants who smoked less than 100 cigarettes in their lifetime. Former smokers were those who had given up smoking before the interview, and current smokers were those who smoked more than 100 cigarettes in their whole life and kept the habit of smoking at the time of interview. The participants who had regular exercise were defined as those who reported that they had moderate- or vigorous-intensity physical activities at least 10 min in a typical week, in succession with small or large increases in heart rate or breathing. Those who were exposed to oral corticosteroid medication for more than 180 days were defined as the presence of use of oral corticosteroid. Hypertension was defined as elevated blood pressure (systolic/diastolic blood pressure equal or higher than 140/90 mm Hg), self-reported hypertension diagnosis by clinician or taking anti-hypertensive drugs. Diabetes mellitus was defined as fasting plasma glucose concentration ≥ 7.0 mmol/L, glycosylated hemoglobin level ≥ 6.5%, self-reported diabetes diagnosis, or use of diabetic pills (including insulin). Dyslipidemia was defined if any of the following status was matched: (1) Total cholesterol ≥ 200 mg/dl, (2) triglyceride ≥ 150 mg/dl, (3) low-density lipoprotein cholesterol ≥ 130 mg/dl, (4) high-density lipoprotein cholesterol < 40 mg/dl or <50 mg/dl for men and women, respectively, (5) self-reported taking prescribed lipid-modifying medication. Cardiovascular disease was defined as self-reported diagnosis of angina, congestive heart failure, coronary heart disease, heart attack, or stroke. Cancer was defined as self-reported diagnosis of any kind of cancer by clinician during the whole lifetime.

Three consecutive blood pressure measurements were obtained after resting quietly in a seated position for 5 min in the MEC according to the physician examination procedures manual (21). A fourth determination would be taken if a blood pressure measurement was interrupted or incomplete. The values of these three or four readings were averaged to represent their blood pressure status. Blood collection took place in the MEC under standardized conditions at each survey location including collecting, processing, storing, and shipping. Laboratory parameters including alanine aminotransferase (ALT), alkaline phosphatase (ALP), aspartate aminotransferase (AST), and gamma-glutamyl transferase (GGT), glucose, glycosylated hemoglobin, total cholesterol, triglyceride, low-density lipoprotein cholesterol, and high-density lipoprotein cholesterol were measured using corresponding methods described in the official website (22).

Statistical analysis

Given to the complex sampling design, appropriate sample weight was conducted according to NHANES analytic guidelines in the current analysis. Non-normally distributed data were natural logarithm transformed and estimates of fatty acids intakes were adjusted for energy intakes using the residuals method (23) before further analysis. The basic characteristics and dietary intakes of fatty acids of participants (overall, non-NASH, and NASH) were described by the weighted mean and standard error (SE) for continuous variables, as well as counts and weighted frequencies for categorical variables. The differences of basic characteristics between participants with and without NASH were compared by general linear models and Chi-squared test as appropriate. Analyses of covariance (ANCOVA) controlling for sex, age, and BMI were used for comparison of the mean differences in dietary intakes of energy and fatty acids.

Participants were divided into three groups according to terciles in the non-NASH group of dietary intakes of total SFAs, MUFAs, PUFAs and their subtypes. Logistic regression models were performed to examine the association between the dietary intakes of fatty acids and the risk of NASH. According to the previous studies (24) and the specialized knowledge, several covariates were selected for adjustments to minimize the residual confounding. Minimally adjusted models included age (continuous), sex (categorized), and BMI (continuous). Other potential risk factors, including ethnicity (categorized), marital status (categorized), education levels (categorized), family income-to-poverty ratio (categorized), waist circumference (continuous), smoking status (categorized), regular exercise (categorized), use of oral corticosteroid (categorized), HEI-2015 (continuous), ALT (continuous), ALP (continuous), AST (continuous), GGT (continuous), prevalence of hypertension, diabetes mellitus, dyslipidemia, cardiovascular disease, and cancer (categorized) were additionally adjusted in the fully adjusted models. Odds ratios (ORs) and 95% confidence intervals (CIs) were calculated with the lowest terciles as the reference. Since only a few covariates were missing with a small portion, observations with missing data were automatically excluded from the corresponding adjusted models. Sensitivity analysis was conducted among participants without use of oral corticosteroid (n = 4,118). Stratified analysis was performed to examine whether the association between terciles of dietary fatty acids intakes and risk of NASH was different in participants with various characteristics. Interactions were estimated by including the multiplicative interaction terms.

The data were analyzed from January 2022 to May 2022. Statistical analyses were performed using the R software 4.1.0 (the “survey” package) and the SPSS, version 25.0 (IBM Corp., Armonk, NY, United States). All tests were two-sided and p < 0.05 was considered statistically significant.

Results

Basic characteristics and dietary intakes of fatty acids of study participants

Weighted distributions of sociodemographic information, lifestyle, laboratory parameters, and prevalence of several chronic diseases for overall population, participants with and without NASH were shown in Table 1. Of the 4,161 study participants, 48.6% were men and mean (SE) age was 47.5 (0.8) years. The participants in the current study tended to be obese with a mean BMI of 29.8 (0.3) kg/m², and a total of 2,089 (50.2%) participants with NASH were identified. Compared with non-NASH participants, those with NASH were more likely to be older, had a higher BMI, higher waist circumference, higher levels of ALT, ALP, AST, and GGT, higher prevalence of hypertension, diabetes mellitus, dyslipidemia, cardiovascular disease, and cancer, had lower education levels, and lower frequencies of regular exercise (p < 0.05). A significantly greater proportion of participants with NASH were men, were Mexican American, were married, and were former or current smokers (p < 0.05). No significant differences were observed in family income-to-poverty ratio and use of oral corticosteroid between participants with and without NASH (p > 0.05).

TABLE 1

Table 1. Basic characteristics of participants of this study.

With regard to dietary intakes of overall participants, energy intakes controlling for sex, age, and BMI were 1,980.53 (12.21) kcal/day. Mean (SE) of dietary total SFAs, MUFAs and PUFAs adjusted for energy intakes, sex, age, and BMI were 23.21 (0.11) g/day, 24.97 (0.11) g/day, and 17.21 (0.09) g/day (Table 2). Compared with non-NASH participants, those with NASH had lower dietary intakes of total MUFAs, total PUFAs, MUFA 18:1, PUFA 18:2, PUFA 20:4, and lower score in HEI-2015 (p < 0.05). Other dietary intakes including energy intake, total SFAs, the ratio of UFAs and SFAs, and other 16 subtypes of fatty acids did not show significant difference between participants with and without NASH.

TABLE 2

Table 2. Dietary fatty acids intakes of participants of this study.

Dietary intakes of fatty acids and risk of non-alcoholic steatohepatitis

Associations between terciles of dietary total SFAs, MUFAs, PUFAs, their subtypes, the ratio of UFAs to SFAs and NASH were presented in Table 3. An inverse association between dietary total PUFAs and NASH risk was found, with an OR of 0.67 (95% CI: 0.46–0.97) at the highest tercile in comparison with the bottom tercile after adjustments for potential risk factors. Similar associations were found between the subtype of PUFAs 18:3 and NASH, while the fully adjusted OR in the highest tercile was 0.67 (95% CI: 0.47–0.96). Participants who were in the second tercile of total MUFAs and MUFA 18:1, rather than in the highest tercile, had a lower risk of NASH compared with the first tercile after full adjustments. However, dietary intakes of total SFAs, other subtypes, as well as the ratio of UFAs to SFAs did not show significant associations with NASH risk after adjustment for potential risk factors in any tercile. Moreover, results remained largely unchanged in the sensitivity analysis when restricting participants to those who did not use oral corticosteroid (Supplementary Table 3).

TABLE 3

Table 3. Odds ratios (ORs) and 95% confidence intervals (CIs) of non-alcoholic steatohepatitis by terciles of dietary intakes of fatty acids among controls (n = 4,161).

Interactions and stratified analyses

Influences of dietary total PUFAs intakes on risk of NASH stratified by selected potential risk factors were shown in Table 4. Multiplicative interactions were statistically significant between terciles of dietary intakes of total PUFAs and BMI (<30 kg/m² or ≥30 kg/m²) on associations with risk of NASH (p-interaction = 0.004). Stronger associations of dietary total PUFAs intakes with NASH risk were found in obese participants with a BMI ≥ 30 kg/m² (OR, 95% CI: 0.41, 0.22–0.75), while significant associations could not be observed in non-obese participants with a BMI < 30 kg/m² (OR, 95% CI: 1.00, 0.70–1.43). Null significant multiplicative interactions between other potential risk factors and terciles of dietary total PUFAs were identified on risk of NASH. Similar effects on risk of NASH were also observed between BMI and dietary intakes of PUFA 18:3 (p-interaction = 0.015, Figure 2). Moreover, inverse associations of dietary intakes of PUFA 18:3 and NASH risk were stronger among participants without presence of cardiovascular disease (OR, 95% CI: 0.66, 0.47–0.94), compared to their counterparts (OR, 95% CI: 0.64, 0.25–1.63, p-interaction = 0.042).

TABLE 4

Table 4. Odds ratios (ORs) and 95% confidence intervals (CIs) of non-alcoholic steatohepatitis by terciles of dietary total polyunsaturated fatty acids among controls stratified by covariates.

FIGURE 2

Figure 2. Odds ratios (ORs) and 95% confidence intervals (CIs) of non-alcoholic steatohepatitis by terciles of dietary PUFA 18:3 among controls stratified by covariates. Abbreviations: PUFA, polyunsaturated fatty acid NASH, non-alcoholic steatohepatitis; HEI-2015, healthy eating index-2015; and BMI, body mass index. Cut-off points of waist circumference were 102 cm for men and 88 cm for women. Odds ratios were adjusted for age, sex, BMI, ethnicity, marital status, education levels, family income-to-poverty ratio, waist circumference, smoking status, regular activities, use of oral corticosteroid, energy intakes, HEI-2015, ALT, ALP, AST, GGT, hypertension, diabetes mellitus, dyslipidemia, cardiovascular disease, and cancer. Stratified factors were not included in the corresponding models.

Discussion

In this nationally representative study with 4,161 adults in the United States, we explored the associations between risk of NASH and dietary intakes of fatty acids, including total SFAs, MUFAs, PUFAs, their common subtypes, and the ratio of UFAs to SFAs. The dietary intakes of total PUFAs, as well as its subtype of PUFA 18:3, were inversely associated with risk of NASH, independent of several potential covariates, especially in those with obesity (BMI ≥ 30 kg/m²). No significant associations were observed between NASH risk and dietary total SFAs, MUFAs, their subtypes as well as the ratio of UFAs to SFAs.

A few studies have explored the relationships between fatty acids and risk of liver diseases or liver-related indices; however, the results were still inconsistent and inconclusive. For example, high abundance of hepatic total SFAs and MUFAs was observed in humans with NAFLD and mice with NASH, which implied that higher total SFAs and MUFAs may be associated with hepatic lipotoxicity and inflammation (25). However, intervention study using diets full with MUFAs showed decreased cholesterol, triglycerides, and increased HDL-cholesterol levels in participants with NAFLD (26). In this study, no significant associations were found between dietary total SFAs, MUFAs, and NASH risk. In addition, an animal study found that moderate intakes of fatty acids with the high ratio of UFAs to SFAs could inhibit liver lipogenesis and steatosis (27), although no significant association was found between the ratio of UFAs to SFAs and NASH risk in our study. Further studies need to be conducted to confirm these findings.

In regard to the relationship between total PUFAs and risk of NAFLD, a systematic review and meta-analysis of 13 studies, consisting of 668 patients with NAFLD, found that total PUFAs or fish oil supplementation may affect serum ALT levels and improve liver function (28). Similar associations could also be observed in other studies of dietary or supplementation with n-3 fatty acids (29–32). However, contradictory results could also be found. A case–control study conducted in 971 Chinese Han adults found that total PUFAs intakes were positively associated with the risk of NAFLD (10). Another cross-sectional study of 233 American children found that dietary long-chain n-3 fatty acids were inversely associated with portal and lobular inflammation, although no significant effects could be found on NASH, which was assessed using serum ALT and histological parameters (33).

The reasons including sample sizes, population, methods of dietary assessment, detection methods of NAFLD may account for diversities among the findings of the aforementioned studies. Particularly, none of those studies have explored the associations of their subtypes with risk of NASH, which was necessary since different fatty acids subtypes exerted different or even opposite effects on liver health (15), due to discrepancies of the length of carbon chain, straight or branched chain, position and numbers of double bonds (34). To the best of our knowledge, this is the first study to investigate associations between dietary total SFAs, MUFAs, PUFAs, their subtypes, the ratio of UFAs to SFAs, and risk of NASH. No significant associations were observed between NASH risk and dietary total SFAs, MUFAs, their subtypes as well as the ratio of UFAs to SFAs. However, dietary intakes of total PUFAs, as well as its subtype of PUFA 18:3, were inversely associated with risk of NASH, independent of several potential covariates, especially in those with obesity (BMI ≥ 30 kg/m²). Several biologic mechanisms could explain the favorable associations between dietary intakes of total PUFAs, its subtype of PUFA 18:3, and the risk of NASH. NASH was characterized as lipid deposition and hepatic inflammation. Oxidative stress, insulin resistance, lipid peroxidation, abnormal secretion of proinflammatory cytokines and adipokines, and intestinal dysbiosis played important roles in the occurrence and development of NAFLD (8, 35–38). In line with our findings on the protective impact of total PUFAs on NASH, total PUFAs have been shown to exert anti-inflammatory effects in both in vitro and in vivo studies. Supplementation with n-3 PUFAs inhibited lipogenesis, attenuated hepatic oxidative stress, decreased inflammation and increased insulin sensitivity, further to preserve hepatic architecture and prevent hepatic steatosis (39). The n-3 PUFAs could downregulate sterol regulatory element-binding protein-1c (SREBP-1c) and upregulate peroxisome proliferator-activated receptor-alpha (PPAR-α) function, and therefore decreased de novo lipogenesis and increased free fatty acids oxidation, improved the biochemical and ultrasonographic manifestations of patients with NAFLD (40). Rats fed with a high-fat, high-calorie solid diet were observed with an increased expression of hepatic adiponectin and PPAR-α, a reduction of tumor necrosis factor-alpha (TNF-α), and therefore an improvement of fatty liver and the degree of liver injury after supplementation of n-3 PUFAs (41).

Polyunsaturated fatty acid 18:3, so called α-linolenic acid, was a kind of important n-3 PUFAs with anti-inflammatory and antioxidant effects for the human body. In a 6-month, randomized, placebo-controlled, double-blind trial, supplementation of n-3 PUFAs impacted on plasma lipid profile in patients with NASH, which was specific in reduction of triglycerides, and therefore improved plasma lipidomic and hepatic proteomic markers of lipogenesis, mitochondrial functions and endoplasmic reticulum stress (42, 43). Further mechanisms about effects of PUFA 18:3 on inflammation need to be certified in future studies.

There may be several kinds of common risk factors and links between overweight, obesity, other related metabolic diseases and NAFLD (44). In our study, an inverse association of dietary intakes of PUFA 18:3 with NASH was observed only among participants without cardiovascular disease. Increasing dietary intakes of PUFA 18:3 may reduce the risk of NASH in participants without presence of cardiovascular disease. The BMI was another risk factor for NASH, and decreases in hepatic fat content were partially attributed to favorable changes in BMI (45). Interactions could be found between dietary PUFAs and obesity on NASH risk in our study. Stronger associations of dietary PUFAs intakes with NASH risk were found in participants with obesity rather than those without obesity, suggesting that obese people might benefit more from increasing dietary intakes of PUFAs.

Our study has notable strengths. First, our observational study was based on a national, representative population with large sample size in the United States. Complex sampling design as well as appropriate sample weight method we conducted increased the reliability and generalizability of our findings. In addition, NASH was diagnosed by transient elastography, which was a non-invasive and fast method with high sensitivity and specificity by comparison with liver biopsy (17–19).

Several limitations need to be considered. First, although we adjusted for multiple potential confounders, including demographic information, lifestyle, medication use and history of chronic diseases, residual confounding cannot be eliminated fully. Furthermore, stratified analysis performed in our study may result in potential statistical power loss, it was desirable to conduct studies with larger sample sizes to validate our finding in the future. In addition, since the progression of NASH was long, dietary assessment in a long term was more appropriate to explore the relationships between dietary intakes of fatty acids and risks of NASH. However, in the analysis of this study, dietary intakes of fatty acids were assessed with a 24-h food recall method in two inconsecutive days, future studies with repeated dietary assessment in a long term are warranted. Moreover, restricted by the observational study design, we cannot definitively exclude the possibility that our findings may be affected by residual confounding and reverse causality. Further prospective or interventional studies need to be carried out to verify our finding about associations between dietary fatty acids and NASH.

Conclusion

In conclusion, inverse associations were observed between dietary intakes of total PUFAs, as well as its subtype of PUFA 18:3, and risk of NASH after adjusting potential confounders. Further large prospective studies need to be conducted to confirm our findings.

Data availability statement

The original contributions presented in this study are included in the article/Supplementary material, further inquiries can be directed to the corresponding authors.

Ethics statement

The studies involving human participants were reviewed and approved by National Center for Health Statistics. The patients/participants provided their written informed consent to participate in this study.

Author contributions

Z-YL and H-LZ created the analytic design. Z-YL performed data extraction. X-TL, Z-YL, Y-DW, and T-TZ analyzed the data. X-TL drafted the manuscript. X-TL, Z-YL, and H-LZ critically revised the manuscript. All authors read and approved the final manuscript.

Funding

This work was funded by the Basic and Applied Basic Research Foundation of Guangdong Province, China (Grant No. 2020A1515110682).

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.

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.

Acknowledgments

We thank all the participants, the investigators, and the staff of the National Health and Nutrition Examination Survey for their valuable contribution.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnut.2022.952451/full#supplementary-material

Abbreviations

NAFLD, non-alcoholic fatty liver disease; NASH, non-alcoholic steatohepatitis; SFA, saturated fatty acid; UFA, unsaturated fatty acid; MUFA, monounsaturated fatty acid; PUFA, polyunsaturated fatty acid; NHANES, national health and nutrition examination survey; MEC, mobile examination center; CAP, controlled attenuation parameter; HEI-2015, healthy eating index-2015; BMI, body mass index; ALT, alanine aminotransferase; ALP, alkaline phosphatase; AST, aspartate aminotransferase; GGT, gamma-glutamyl transferase; SE, standard error; OR, odds ratio; CI, confidence interval; SREBP-1c, sterol regulatory element-binding protein-1c; PPAR- α, peroxisome proliferator-activated receptor-alpha; TNF- α, tumor necrosis factor-alpha; NLRP3, NOD-like receptor protein 3.

Footnotes

1. ^ https://www.cdc.gov/nchs/nhanes/

References

1. Paik JM, Golabi P, Younossi Y, Mishra A, Younossi ZM. Changes in the global burden of chronic liver diseases from 2012 to 2017: the growing impact of Nafld. Hepatology. (2020) 72:1605–16. doi: 10.1002/hep.31173

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Younossi ZM, Stepanova M, Afendy M, Fang Y, Younossi Y, Mir H, et al. Changes in the prevalence of the most common causes of chronic liver diseases in the United States from 1988 to 2008. Clin Gastroenterol Hepatol. (2011) 9:524–30e1;quize60. doi: 10.1016/j.cgh.2011.03.020

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Estes C, Anstee QM, Arias-Loste MT, Bantel H, Bellentani S, Caballeria J, et al. Modeling Nafld disease Burden in China, France, Germany, Italy, Japan, Spain, United Kingdom, and United States for the Period 2016–2030. J Hepatol. (2018) 69:896–904. doi: 10.1016/j.jhep.2018.05.036

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Komolafe O, Buzzetti E, Linden A, Best LM, Madden AM, Roberts D, et al. Nutritional supplementation for nonalcohol-related fatty liver disease: a network meta-analysis. Cochrane Database Syst Rev. (2021) 7:CD013157. doi: 10.1002/14651858.CD013157.pub2

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Kardashian A, Dodge JL, Terrault NA. Food insecurity is associated with mortality among U.S. Adults with nonalcoholic fatty liver disease and advanced fibrosis. Clin Gastroenterol Hepatol. (2021) S1542–3565:1267–1262. doi: 10.1016/j.cgh.2021.11.029

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Piras C, Noto A, Ibba L, Deidda M, Fanos V, Muntoni S, et al. Contribution of metabolomics to the understanding of nafld and nash syndromes: a systematic review. Metabolites. (2021) 11:694. doi: 10.3390/metabo11100694

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Younossi Z, Anstee QM, Marietti M, Hardy T, Henry L, Eslam M, et al. Global burden of nafld and nash: trends, predictions, risk factors and prevention. Nat Rev Gastroenterol Hepatol. (2018) 15:11–20. doi: 10.1038/nrgastro.2017.109

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Musso G, Cassader M, Paschetta E, Gambino R. Bioactive lipid species and metabolic pathways in progression and resolution of nonalcoholic steatohepatitis. Gastroenterology. (2018) 155:282–302e8. doi: 10.1053/j.gastro.2018.06.031

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Willis SA, Bawden SJ, Malaikah S, Sargeant JA, Stensel DJ, Aithal GP, et al. The role of hepatic lipid composition in obesity-related metabolic disease. Liver Int. (2021) 41:2819–35. doi: 10.1111/liv.15059

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Xie Y, Tian H, Xiang B, Li D, Liu J, Cai Z, et al. Total polyunsaturated fatty acid intake and the risk of non-alcoholic fatty liver disease in chinese han adults: a secondary analysis based on a case-control study. BMC Gastroenterol. (2021) 21:451. doi: 10.1186/s12876-021-02039-2

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Lee JH, Kim D, Kim HJ, Lee CH, Yang JI, Kim W, et al. Hepatic steatosis index: a simple screening tool reflecting nonalcoholic fatty liver disease. Dig Liver Dis. (2010) 42:503–8. doi: 10.1016/j.dld.2009.08.002

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Bedogni G, Bellentani S, Miglioli L, Masutti F, Passalacqua M, Castiglione A, et al. The fatty liver index: a simple and accurate predictor of hepatic steatosis in the general population. BMC Gastroenterol. (2006) 6:33. doi: 10.1186/1471-230X-6-33

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Kotronen A, Peltonen M, Hakkarainen A, Sevastianova K, Bergholm R, Johansson LM, et al. Prediction of non-alcoholic fatty liver disease and liver fat using metabolic and genetic factors. Gastroenterology. (2009) 137:865–72. doi: 10.1053/j.gastro.2009.06.005

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Poynard T, Ratziu V, Naveau S, Thabut D, Charlotte F, Messous D, et al. The diagnostic value of biomarkers (Steatotest) for the prediction of liver steatosis. Comp Hepatol. (2005) 4:10. doi: 10.1186/1476-5926-4-10

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Yamada K, Mizukoshi E, Sunagozaka H, Arai K, Yamashita T, Takeshita Y, et al. Characteristics of hepatic fatty acid compositions in patients with nonalcoholic steatohepatitis. Liver Int. (2015) 35:582–90. doi: 10.1111/liv.12685

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Centers for Disease Control and Prevention.National Health and NutritionExamination Survey (NHANES) 2017–2018 Data Documentation, Codebook, and Frequencies. (2020). Available online at: https://wwwn.cdc.gov/Nchs/Nhanes/2017-2018/LUX_J.htm (accessed June 13, 2022).

Google Scholar

17. Myers RP, Pollett A, Kirsch R, Pomier-Layrargues G, Beaton M, Levstik M, et al. Controlled attenuation parameter (Cap): a noninvasive method for the detection of hepatic steatosis based on transient elastography. Liver Int. (2012) 32:902–10. doi: 10.1111/j.1478-3231.2012.02781.x

PubMed Abstract | CrossRef Full Text | Google Scholar

18. de Ledinghen V, Wong GL, Vergniol J, Chan HL, Hiriart JB, Chan AW, et al. Controlled attenuation parameter for the diagnosis of steatosis in non-alcoholic fatty liver disease. J Gastroenterol Hepatol. (2016) 31:848–55. doi: 10.1111/jgh.13219

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Sasso M, Audiere S, Kemgang A, Gaouar F, Corpechot C, Chazouilleres O, et al. Liver steatosis assessed by controlled attenuation parameter (Cap) measured with the Xl probe of the fibroscan: a pilot study assessing diagnostic accuracy. Ultrasound Med Biol. (2016) 42:92–103. doi: 10.1016/j.ultrasmedbio.2015.08.008

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Siddiqui MS, Vuppalanchi R, Van Natta ML, Hallinan E, Kowdley KV, Abdelmalek M, et al. Vibration-controlled transient elastography to assess fibrosis and steatosis in patients with nonalcoholic fatty liver disease. Clin Gastroenterol Hepatol. (2019) 17:156.e–63.e. doi: 10.1016/j.cgh.2018.04.043

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Centers for Disease Control and Prevention.National Health and Nutrition Examination Survey (NHANES) 2017-2018 Procedure Manuals. (2020). Available online at: https://wwwn.cdc.gov/nchs/nhanes/continuousnhanes/manuals.aspx?BeginYear=2017 (accessed June 13, 2022).

Google Scholar

22. Centers for Disease Control and Prevention.National Health and NutritionExamination Survey (NHANES) 2017-2018 Laboratory Methods. (2020). Available online at: https://wwwn.cdc.gov/Nchs/Nhanes/continuousnhanes/labmethods.aspx?BeginYear=2017 (accessed June 13, 2022).

Google Scholar

23. Willett W, Stampfer MJ. Total energy intake: implications for epidemiologic analyses. Am J Epidemiol. (1986) 124:17–27. doi: 10.1093/oxfordjournals.aje.a114366

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Kim D, Konyn P, Cholankeril G, Ahmed A. Physical activity is associated with nonalcoholic fatty liver disease and significant fibrosis measured by fibroscan. Clin Gastroenterol Hepatol. (2022) 20: e1438–e1455. doi: 10.1016/j.cgh.2021.06.029

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Spooner MH, Jump DB. Omega-3 Fatty acids and nonalcoholic fatty liver disease in adults and children: where do we stand? Curr Opin Clin Nutr Metab Care. (2019) 22:103–10. doi: 10.1097/MCO.0000000000000539

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Barrera F, George J. The role of diet and nutritional intervention for the management of patients with nafld. Clin Liver Dis. (2014) 18:91–112. doi: 10.1016/j.cld.2013.09.009

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Siddiqui RA, Xu Z, Harvey KA, Pavlina TM, Becker MJ, Zaloga GP. Comparative study of the modulation of fructose/sucrose-induced hepatic steatosis by mixed lipid formulations varying in unsaturated fatty acid content. Nutr Metab (Lond). (2015) 12:41. doi: 10.1186/s12986-015-0038-x

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Yu L, Yuan M, Wang L. The effect of Omega-3 unsaturated fatty acids on non-alcoholic fatty liver disease: a systematic review and meta-analysis of rcts. Pak J Med Sci. (2017) 33:1022–8. doi: 10.12669/pjms.334.12315

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Musa-Veloso K, Venditti C, Lee HY, Darch M, Floyd S, West S, et al. Systematic review and meta-analysis of controlled intervention studies on the effectiveness of long-chain omega-3 fatty acids in patients with nonalcoholic fatty liver disease. Nutr Rev. (2018) 76:581–602. doi: 10.1093/nutrit/nuy022

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Mitrovic M, Sistilli G, Horakova O, Rossmeisl M. Omega-3 phospholipids and obesity-associated nafld: potential mechanisms and therapeutic perspectives. Eur J Clin Invest. (2022) 52:e13650. doi: 10.1111/eci.13650

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Lu W, Li S, Li J, Wang J, Zhang R, Zhou Y, et al. Effects of omega-3 fatty acid in nonalcoholic fatty liver disease: a meta-analysis. Gastroenterol Res Pract. (2016) 2016:1459790. doi: 10.1155/2016/1459790

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Corte CD, Iasevoli S, Strologo AD, Sanseviero M, Nobili V. Omega-3 fatty acids and fatty liver disease in children. Adv Food Nutr Res. (2018) 85:59–77. doi: 10.1016/bs.afnr.2018.03.001

PubMed Abstract | CrossRef Full Text | Google Scholar

33. St-Jules DE, Watters CA, Brunt EM, Wilkens LR, Novotny R, Belt P, et al. Estimation of fish and omega-3 fatty acid intake in pediatric nonalcoholic fatty liver disease. J Pediatr Gastroenterol Nutr. (2013) 57:627–33. doi: 10.1097/MPG.0b013e3182a1df77

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Kazaz S, Miray R, Lepiniec L, Baud S. Plant monounsaturated fatty acids: diversity, biosynthesis, functions and uses. Prog Lipid Res. (2022) 85:101138. doi: 10.1016/j.plipres.2021.101138

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Videla LA, Valenzuela R. Perspectives in liver redox imbalance: toxicological and pharmacological aspects underlying iron overloading, nonalcoholic fatty liver disease, and thyroid hormone action. Biofactors. (2022) 48:400–15. doi: 10.1002/biof.1797

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Pafili K, Roden M. Nonalcoholic fatty liver disease (Nafld) from pathogenesis to treatment concepts in humans. Mol Metab. (2021) 50:101122. doi: 10.1016/j.molmet.2020.101122

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Liu W, Baker RD, Bhatia T, Zhu L, Baker SS. Pathogenesis of nonalcoholic steatohepatitis. Cell Mol Life Sci. (2016) 73:1969–87. doi: 10.1007/s00018-016-2161-x

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Byrne CD. Fatty liver: role of inflammation and fatty acid nutrition. Prostaglandins Leukot Essent Fatty Acids. (2010) 82:265–71. doi: 10.1016/j.plefa.2010.02.012

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Jeyapal S, Kona SR, Mullapudi SV, Putcha UK, Gurumurthy P, Ibrahim A. Substitution of linoleic acid with alpha-linolenic acid or long chain N-3 polyunsaturated fatty acid prevents western diet induced nonalcoholic steatohepatitis. Sci Rep. (2018) 8:10953. doi: 10.1038/s41598-018-29222-y

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Bouzianas DG, Bouziana SD, Hatzitolios AI. Potential treatment of human nonalcoholic fatty liver disease with long-chain omega-3 polyunsaturated fatty acids. Nutr Rev. (2013) 71:753–71. doi: 10.1111/nure.12073

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Svegliati-Baroni G, Candelaresi C, Saccomanno S, Ferretti G, Bachetti T, Marzioni M, et al. A model of insulin resistance and nonalcoholic steatohepatitis in rats: role of peroxisome proliferator-activated receptor-alpha and n-3 polyunsaturated fatty acid treatment on liver injury. Am J Pathol. (2006) 169:846–60. doi: 10.2353/ajpath.2006.050953

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Nogueira MA, Oliveira CP, Ferreira Alves VA, Stefano JT, Rodrigues LS, Torrinhas RS, et al. Omega-3 polyunsaturated fatty acids in treating non-alcoholic steatohepatitis: a randomized, double-blind, placebo-controlled trial. Clin Nutr. (2016) 35:578–86. doi: 10.1016/j.clnu.2015.05.001

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Okada L, Oliveira CP, Stefano JT, Nogueira MA, Silva I, Cordeiro FB, et al. Omega-3 pufa modulate lipogenesis, er stress, and mitochondrial dysfunction markers in nash - proteomic and lipidomic insight. Clin Nutr. (2018) 37:1474–84. doi: 10.1016/j.clnu.2017.08.031

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Allen AM, Therneau TM, Larson JJ, Coward A, Somers VK, Kamath PS. Nonalcoholic fatty liver disease incidence and impact on metabolic burden and death: A 20 year-community study. Hepatology. (2018) 67:1726–36. doi: 10.1002/hep.29546

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Lee CH, Fu Y, Yang SJ, Chi CC. Effects of omega-3 polyunsaturated fatty acid supplementation on non-alcoholic fatty liver: a systematic review and meta-analysis. Nutrients. (2020) 12:2769. doi: 10.3390/nu12092769

PubMed Abstract | CrossRef Full Text | Google Scholar