Prevalence and determinants of non-alcoholic fatty liver disease in lifelines: A large Dutch population cohort

Abstract

Background & aims

Non-alcoholic fatty liver disease is an increasing health issue that develops rather unnoticed with obesity, type 2 diabetes mellitus and metabolic syndrome. We investigated prevalence, determinants and associated metabolic abnormalities of non-alcoholic fatty liver disease in the largest population-based cohort to date.

Methods

Biochemical characteristics, type 2 diabetes mellitus and metabolic syndrome were determined in the Lifelines Cohort Study (N = 167,729), a population-based cohort in the North of the Netherlands. Non-alcoholic fatty liver disease was defined as Fatty Liver Index (FLI)≥60. Exclusion criteria were age <18 years, immigrants, missing data to assess FLI and metabolic syndrome, excessive alcohol use, previous-diagnosed hepatitis or cirrhosis and non-fasting blood sampling.

Results

Out of 37,496 included participants (median age 44 years, 62.1% female), 8,259 (22.0%) had a FLI≥60. Individuals with a FLI≥60 were more often male, older, obese, had higher levels of hemoglobinA1c, fasting glucose, liver enzymes, total cholesterol, low-density lipoprotein cholesterol, triglycerides, c-reactive protein and leucocytes and lower high-density lipoprotein cholesterol (all P<0.0001). Participants with a FLI≥60 showed higher prevalence of type 2 diabetes mellitus (9.3% vs. 1.4%), metabolic syndrome (54.2% vs. 6.2%), impaired renal function (20.1% vs. 8.7%) and cardiovascular disease (4.6% vs. 1.6%) (all P<0.0001). Multivariable logistic analysis showed that smoking, hemoglobin, leucocytes, c-reactive protein, platelets, alanine aminotransferase, alkaline phosphatase, albumin, impaired renal function (OR 1.27, 95%CI 1.15–1.41), metabolic syndrome (OR 11.89, 95%CI 11.03–12.82) and its individual components hyperglycemia (OR 2.53, 95%CI 2.34–2.72), hypertension (OR 1.89, 95%CI 1.77–2.01) and reduced high-density lipoprotein cholesterol (OR 3.44, 95%CI 3.22–3.68) were independently associated with suspected non-alcoholic fatty liver disease (all P<0.0001).

Conclusion

Twenty-two percent (22.0%) of the population in the North of the Netherlands is suspected to suffer from non-alcoholic fatty liver disease, coinciding with a significant increased risk of type 2 diabetes mellitus, metabolic syndrome, cardiovascular disease and impaired renal function.

About the Authors:

Eline H. van den Berg

Contributed equally to this work with: Eline H. van den Berg, Marzyeh Amini

* E-mail: [email protected]

Affiliation: Department of Gastroenterology and Hepatology, University Medical Center Groningen, University of Groningen, Groningen, The Netherlands

Marzyeh Amini

Contributed equally to this work with: Eline H. van den Berg, Marzyeh Amini

Affiliation: Department of Epidemiology, University Medical Center Groningen, University of Groningen, Groningen, The Netherlands

Tim C. M. A. Schreuder

Affiliation: Department of Gastroenterology and Hepatology, University Medical Center Groningen, University of Groningen, Groningen, The Netherlands

Robin P. F. Dullaart

Affiliation: Department of Endocrinology, University Medical Center Groningen, University of Groningen, Groningen, The Netherlands

Klaas Nico Faber

Affiliation: Department of Gastroenterology and Hepatology, University Medical Center Groningen, University of Groningen, Groningen, The Netherlands

Behrooz Z. Alizadeh

Affiliations Department of Gastroenterology and Hepatology, University Medical Center Groningen, University of Groningen, Groningen, The Netherlands, Department of Epidemiology, University Medical Center Groningen, University of Groningen, Groningen, The Netherlands

Hans Blokzijl

Affiliation: Department of Gastroenterology and Hepatology, University Medical Center Groningen, University of Groningen, Groningen, The Netherlands

Introduction

Non-alcoholic fatty liver disease (NAFLD) is characterized by hepatic steatosis in the absence of excessive alcohol consumption. The spectrum of NAFLD ranges from simple steatosis to non-alcoholic steatohepatitis (NASH), fibrosis and ultimately cirrhosis with its known complications, such as decompensation and hepatocellular carcinoma (HCC)[1]. In patients with NASH, progression to fibrosis occurs in 40.8% with a liver specific mortality hazard ratio of 2.6[2]. As a result of the global obesity epidemic, NAFLD is an increasing relevant public health issue and emerging as the most common cause of chronic liver disease in Western countries. It is expected to become the most important indication for liver transplantation in the near future[3]. Although most patients with NAFLD are not at risk of dying from liver disease, they have a substantial increased risk of early morbidity and mortality[1,4,5]. NAFLD frequently co-exists with metabolic disorders and the association with the metabolic syndrome (MetS) is strong[6]. Another condition associated with NAFLD is cardiovascular disease, with increased intima-media thickness and carotid plaques representing progressive atherosclerosis[7].

In European countries the prevalence of NAFLD has been reported to range widely from 3.1–41.2% (S1 Table), and is likely to increase over the coming years[8–24]. Nonetheless, there are only few, mostly small-sized European epidemiological studies analyzing the prevalence of NAFLD in the general population[8,9,13,20], with the largest survey being performed in only 4,222 participants[20]. Other studies were performed in selected categories adapted from a general population[10–12,14–19,21–24].

Given considerable variation in reported prevalence numbers derived from rather small cohorts, and the increasing incidence of NAFLD with its serious consecutive complications and comorbidity, the present study was initiated to establish a comprehensive sufficiently powered analysis on the prevalence of NAFLD. Here, we aimed to investigate the prevalence, determinants and comorbid conditions of NAFLD in a large population-based cohort from the North of the Netherlands.

Methods

Study design

This cross-sectional study was conducted within the framework of the Lifelines Cohort Study[25–27]. The Lifelines Cohort Study is a multi-disciplinary prospective population-based cohort study of 167,729 persons living in the North of the Netherlands. It employs a broad range of investigative procedures in assessing the biomedical, socio-demographic, behavioral, physical and psychological factors which contribute to the health and disease of the general population, with a special focus on multi-morbidity and complex genetics. Participants were recruited via general practitioners, subsequently family members were invited to participate and finally, adults could self-register to participate. All participants provided written informed consent. The medical ethics committee of the University of Groningen, the Netherlands, approved the study[25–27].

Study participants

Subjects of Western-European origin were included. All study participants were aged between 18–91 years at time of enrollment. Exclusion criteria were participants <18 years, those with missing data required to calculate the Fatty Liver Index (FLI)[28] (described below) and to determine MetS components, non-fasting participants at time of blood collection, immigrants, participants with self-reported excessive alcohol use and those previous diagnosed with hepatitis or cirrhosis. Information about nationality, fasting state, smoking, medication use, alcohol consumption, hepatitis B virus infection and cirrhosis was extracted from the self-administered questionnaires. Participants were assumed to be of Western-European origin if his/her birth country and that of both parents was the Netherlands, which is in accordance with the definition of Statistics of the Netherlands[27]. Participants were considered normal drinkers when daily alcoholic intake was ≤1 drink in females and ≤2 drinks in males[29]. Current smokers consisted of participants with active smoking or smoking in the past month.

Data collection and measurements

Data was collected in the Lifelines Cohort Study between 2006–2013. Questionnaires were collected, anthropometry and blood pressure were measured and biomaterial (blood) was collected at the Lifelines research sites. A standardized protocol was used to obtain blood pressure and anthropometric measurements (height, weight and waist circumference). Systolic and diastolic blood pressures were measured 10 times during a period of 10 minutes, using an automated Dinamap Monitor (GE Healthcare, Freiburg, Germany). The size of the cuff was chosen according to the arm circumference. The average of the final three readings was used for each blood pressure parameter. Anthropometric measurements were measured without shoes. Body weight was measured to the nearest 0.5 kg. Height and waist circumference were measured to the nearest 0.5 cm. Height was measured with a stadiometer placing their heels against the rod and the head in Frankfort Plane position. Waist circumference was measured in standing position with a tape measure all around the body at the level midway between the lower rib margin and the iliac crest[25,26].

Venous blood samples were collected between 8.00–10.00 a.m. into heparin-containing tubes, centrifuged at 1,885xg and the plasma aliquots were processed for laboratory measurements at the same day and stored at -80°C. Hemoglobin, total leucocytes and platelets were measured using routine procedures on a XE2100-system (Sysmex, Japan). High-sensitivity c-reactive protein (CRP) was measured with CardioPhase hs CRP (Siemens, BNII, Germany) and from 2012 with CRPL3 on a Roche Modular P chemistry analyzer. Total cholesterol, low-density lipoprotein (LDL) cholesterol, high-density lipoprotein (HDL) cholesterol and triglycerides (TG) were measured using routine procedures on a Roche Modular P chemistry analyzer. Glucose was assayed with the UV-test hexokinase method on a Roche Modular P chemistry analyzer and hemoglobin A1c (HbA1c) was measured with high performance liquid chromatography (HPLC) (Roche). Gamma-glutamyltransferase (GGT), alkaline phosphatase (ALP), alanine aminotransferase (ALT) and aspartate aminotransferase (AST) were quantified according to the recommendation of the International Federation of Clinical Chemistry on a Roche Modular Platform. ALT and AST were measured with pyridoxal phosphate activation. Albumin was measured with a BCG albumin assay kit for colorimetric testing on a Roche Modular P chemistry analyzer. All laboratory measurements were performed with standardized laboratory measurements and quality assessment control at the Department of Laboratory Medicine of the University Medical Center Groningen, the Netherlands[25,26].

Definition of NAFLD

For the diagnosis of NAFLD the algorithm of the Fatty Liver Index (FLI) was used. The FLI was calculated according to the formula published by Bedogni[28]. FLI = (e^0.953*^{loge (triglycerides +0.139}*^BMI+0.718*^{loge (GGT)+0.053}*^{waist circumference–15.745})/(1+e^0.953*^{loge (triglycerides)+0.139}*^{BM +0.718}*^{loge(GGT +0.053}*^{waist circumference–15.745})*100, where GGT is gamma-glutamyltransferase. The optimal cut-off value for the FLI has been documented to be 60 with an accuracy of 0.84, a sensitivity of 61% and a specificity of 86% for detecting NAFLD as determined by ultrasonography[28]. A FLI≥60 was thus used as a proxy of NAFLD. The 2016 EASL-EASD-EASO NAFLD guideline recommends that for larger scale screening studies, serum biomarkers are the preferred diagnostic tool with the FLI currently considered to be one of the best validated steatosis scores[29].

Definition of comorbid diseases

Computational models for the determination of comorbid diseases were used. For the definition of obesity the body mass index (BMI) was used, calculated as weight (kg) divided by height squared (m²). The diagnosis of type 2 diabetes mellitus (T2DM) was confirmed when a subject had either self-reported on T2DM, used glucose lowering medication, had a fasting glucose (FG) ≥7.0 mmol/L or a HbA1c ≥47.5 mmol/mol. MetS was defined by the revised diagnostic criteria from the American Heart Association by the National Cholesterol Education Program Adult Treatment Panel III[30] and consist of five criteria: (1) enlarged waist circumference (males ≥102 cm and females ≥88 cm), (2) elevated TG (≥1.7 mmol/L) and/or medication use for elevated TG, (3) reduced HDL cholesterol (males <1.0 mmol/L and females <1.3 mmol/L) and/or medication use for reduced HDL cholesterol, (4) elevated blood pressure (systolic blood pressure ≥130 mmHg or diastolic blood pressure ≥85 mmHg) and/or medication use for hypertension, (5) elevated fasting glucose (≥5.6 mmol/L) and/or medication use for elevated glucose. Participants were diagnosed with MetS when at least three out of five criteria were present[30]. The presence of a self-reported history of myocardial infarction, percutaneous coronary intervention, coronary artery bypass surgery, stroke or the diagnosis of narrowing of one or both carotid arteries was defined as atherosclerotic cardiovascular disease. Chronic impaired renal function was defined by calculating the estimated glomerular filtration rate (eGFR) < 60 ml/min/1.73m², using the Modification of Diet in Renal Disease (MDRD) Study Equation [31,32].

Data analyses and statistical modeling

Statistical analyses was performed with SPSS (version 22.0, SPSS Inc., Chicago, IL, USA). Data are expressed in means with standard deviations (SD), medians with interquartile ranges (IQR) and in numbers with percentages. Normality of distribution was assessed and checked for skewness. Variables were compared between FLI≥60 and FLI<60 groups using Student T-test, Mann-Whitney U test Chi-square test. To preclude interactions with the dependent factor FLI≥60, all variables in the equation defining the FLI (i.e. BMI, waist circumference, TG and GGT) were excluded in multivariable analyses. Due to correlations ≥0.5; AST (correlation with ALT), glucose (correlation with HbA1c) and total cholesterol (correlation with LDL cholesterol) were excluded from multivariable analyses and residual variables were made to exclude remaining interactions. For continuous variables a Z-score was calculated and used in multivariable analyses. Stepwise binary logistic regression analyses was performed to disclose the independent association of a FLI≥60. Results are presented by odds ratio (OR) with 95% confidence intervals (CI). To account for the number of independent tests, we applied a Bonferroni correction. Two-sided P-values of <0.001 (0.05/60) were considered statistically significant, given the use of 60 independent tests embedded in 4 multivariable models.

Results

From the 167,729 participants of the Lifelines Cohort Study, 152,180 participants were older than 18 years and 50,704 participants were eligible for our study with necessary available biomedical data concerning the calculation for the FLI and MetS. After applying exclusion criteria, the final study group consisted of 37,496 participants (Fig 1). The median age of the study group was 44 years, with a median BMI of 25.5 kg/m² and was predominantly female (62.1%). Population characteristics are presented in Table 1.

[Figure omitted. See PDF.]

Fig 1. Flow chart of the study population.

https://doi.org/10.1371/journal.pone.0171502.g001

[Figure omitted. See PDF.]

Table 1. Population characteristics.

https://doi.org/10.1371/journal.pone.0171502.t001

Suspected NAFLD was defined by FLI≥60. Suspected NAFLD was observed in 22.0% (8,259 participants) of the study group. Table 2 shows the clinical and laboratory characteristics in subjects with and without suspected NAFLD (FLI<60). Those with suspected NAFLD were older (median age 47 years) and more likely to be male; corresponding prevalence numbers were 32.7% in all males and 15.7% in all females, respectively. As expected, in the group with suspected NAFLD, more obese participants were detected (median BMI of 30.8 kg/m²) compared to those with a FLI<60 (median BMI of 24.4 kg/m²). T2DM (9.3% vs. 1.4%, P<0.0001) and MetS (54.2% vs. 6.2%, P<0.0001) were more prevalent in subjects with a FLI≥60. Significant differences for each individual MetS component were also present (all P<0.0001). Cardiovascular disease (4.6% vs. 1.6%, P<0.0001) and impaired renal function (20.1% vs. 8.7%, P<0.0001) were also more prevalent in subjects with a FLI≥60. In subjects with a FLI≥60, hemoglobin, total leucocytes, CRP, platelets, ALT, AST, GGT, ALP, HbA1c, FG, LDL cholesterol, TG and total cholesterol values were significantly higher and HDL cholesterol and albumin were lower. After adjusting for age and sex, these differences remained significant (all P<0.001).

[Figure omitted. See PDF.]

Table 2. Clinical and laboratory characteristics in subjects with and without non-alcoholic fatty liver disease estimated by the Fatty Liver Index (FLI ≥ 60).

https://doi.org/10.1371/journal.pone.0171502.t002

In order to disclose the independent associations of a FLI≥60 with clinical and biochemical characteristics subsequent stepwise multivariable logistic regression was performed (Tables 3 and 4). In age- and sex- adjusted analysis, impaired renal function, current smoking, hemoglobin, total leucocytes, CRP, platelets, ALT, ALP, albumin, HDL cholesterol and LDL cholesterol were all independent factors associated with a FLI≥60 (all P<0.01) (Table 3). Of note, HbA1c (OR 1.10, 95%CI 1.07–1.14, P<0.0001) and T2DM (OR 2.31, 95%CI 1.97–2.70, P<0.0001) were both independently associated with a FLI≥60 (Table 3; model 1 vs. model 2). In consecutive analysis, the presence of MetS and its individual components were included (Table 4). Waist circumference, HDL cholesterol, TG, HbA1c and T2DM were excluded to preclude interactions of variables accounted for the MetS components and concurrent presence in the equation of the FLI. After inclusion of MetS in the model, independent associations of a FLI≥60 were found with impaired renal function, smoking, hemoglobin, total leucocytes, CRP, platelets, ALT, ALP and albumin (all P<0.0001) (Table 4, model 1 and 2). LDL cholesterol and a cardiovascular disease history were only significantly associated with FLI when the individual components of MetS were added (P<0.01) (Table 4, model 2). Besides a very strong association of a FLI≥60 with the presence of MetS (OR 11.89, 95%CI 11.03–12.82, P<0.0001) (Table 4; model 1), a FLI≥60 was also independently associated with all of the remaining individual MetS components; hyperglycemia (OR 2.53, 95%CI 2.34–2.72, P<0.0001), hypertension (OR 1.89, 95%CI 1.77–2.01, P<0.0001) and low HDL cholesterol (OR 3.44, 95%CI 3.22–3.68, P<0.0001) (Table 4; model 2).

[Figure omitted. See PDF.]

Table 3. Multivariable logistic regression analyses demonstrating independent associations of non-alcoholic fatty liver disease estimated by the Fatty Liver Index (FLI ≥ 60) with current smoking, HbA1c and type 2 diabetes mellitus.

https://doi.org/10.1371/journal.pone.0171502.t003

[Figure omitted. See PDF.]

Table 4. Multivariable logistic regression analyses demonstrating independent associations of non-alcoholic fatty liver disease estimated by the Fatty Liver Index (FLI ≥ 60) with the presence of current smoking, the metabolic syndrome and its individual components.

https://doi.org/10.1371/journal.pone.0171502.t004

Discussion

In this large population based cross-sectional study among almost 40,000 subjects from the Northern part of the Netherlands, the prevalence of NAFLD and its associated metabolic derangements were studied demonstrating that 22% of this adult Western-European population is suspected to suffer from NAFLD. These individuals were more likely to be men, older and suffering from hypertension, T2DM, MetS, cardiovascular disease and impaired renal function. Laboratory tests revealed significant increased glucose, ALT and ALP levels and decreased HDL cholesterol. Further, current smoking, higher levels of hemoglobin, CRP and total leucocytes count were also independently associated with suspected NAFLD.

Previous European studies that have investigated the prevalence of NAFLD in the general populations reported outcomes ranging from 17.9–29.9% (S1 Table)[8,9,13,20]. Gastaldeli et al. found a prevalence of 17.9% in 1,307 participants from 14 different European countries by the use of the FLI (FLI>60)[8]. Three other single country cohorts of the general populations from Germany[20], Spain[13] and Italy[9] demonstrated a NAFLD prevalence by ultrasonography of 29.9%[20], 25.8%[13] and 22.6%[9], respectively. However, these studies represented only 4,222[20], 766[13] and 598[9] participants. Other small European prevalence studies used specific categories of general populations introducing potential bias (e.g. hospitalized patients, heavy drinkers, obese subjects and deceased patients) (S1 Table)[10–12,14–19,21–24]. A recent meta-analysis of different European prevalence studies (including those with selected subgroups) found an overall NAFLD prevalence of 23.7% in 16,735 included subjects[2], corresponding with the prevalence estimate of 22% in this study. NAFLD is less prevalent in Western-Europe when compared to other regions, which show incremental prevalence in North America (24.1%), Asia (27.4%), South America (30.5%) and the Middle East (31.8%)[2]. To date, the Lifelines cohort study with nearly 40,000 participants is the largest study investigating the prevalence of NAFLD in a Western-European cohort. By additionally excluding immigrants, subjects with excessive alcohol use, as well as previously diagnosed hepatitis or cirrhosis, the presently studied cohort is representative in demonstrating a most accurate prevalence figure and coinciding abnormalities in subjects with suspected NAFLD.

When compared to the European prevalence of T2DM and MetS in subjects with suspected NAFLD, prevalence of T2DM was less (9.3% vs. 17.7%) and MetS was more prevalent (54.2% vs. 38.3%) in the Northern region of the Netherlands[2]. This difference could be explained by other studies including subgroup populations, resulting in a combination of different ethnicities and heterogeneity in diagnostic procedures (radiological imaging, ICD codes, self-reported diagnosis and biomarkers) for establishing NAFLD, T2DM and MetS.

All liver enzymes appeared to be increased in the suspected NAFLD group. For ALT, this could be explained by its association with visceral fat, steatosis, inflammation and fibrosis[5]. Remarkably, within the suspected NAFLD group the medians with IQR and means with standard deviations (Table 2) of these liver enzymes were all within the normal reference range used in daily clinical practice. When using the upper limit of normal ALT, 80.3% of subjects in the suspected NAFLD group had normal ALT levels. Others have confirmed these findings. In 79% of subjects with hepatic steatosis[33] and in up to 59% of those with NASH and advanced fibrosis, normal ALT levels were found[34]. This clearly demonstrates the limitations in using ALT levels as a surrogate marker for diagnosing NAFLD and discriminating simple steatosis from steatohepatitis.

A strong association between current smoking, hemoglobin, inflammatory markers (e.g. CRP and total leucocyte count) and suspected NAFLD was found. An association of smoking with NAFLD has not been uniformly reported[5], but may be a confounding environmental stressor. Previous studies have demonstrated that smokers have a higher BMI, increased insulin resistance and that smoking is associated with central fat accumulation, dyslipidemia and concomitant T2DM and MetS, predisposing to comorbidities and risk factors for NAFLD[5,35]. Smoking has been linked to increased hepatic lipid accumulation by modulating the activity of AMPK and SREBP-1, which represent pathways involved in lipid synthesis[36]. The association of a higher hemoglobin level and NAFLD has also been previously demonstrated[37,38] and has been related to progression of NAFLD to NASH and fibrosis[39]. Suggested mechanisms resulting in increased hemoglobin levels are hepatic hypoxia, oxidative stress, formation of reactive oxygen species and lipid peroxidation[37,39]. The association of (subclinical) elevated inflammatory markers and the presence of NAFLD has also been reported in other studies[40]. This may be explained by increased visceral adipose tissue conferring a pro-inflammatory state[41,42]. Also, hepatic free fatty acid oxidation generates oxygen radicals with subsequent lipid peroxidation, cytokine induction and mitochondrial dysfunction, which all conceivably promote inflammation and cause hepatocyte apoptosis and cellular injury. Finally, genetic and gut-derived bacterial factors (in combination with increased intestinal permeability) have an impact on systemic low-grade inflammation[41,43].

The FLI score was used to discriminate between suspected NAFLD and non-NAFLD in this study. The FLI is a well-accepted diagnostic tool for NAFLD, but it is clear that the FLI score is not an absolute measure of hepatic fat accumulation. While histological assessment of liver tissue is still the golden standard for diagnosing NAFLD, liver biopsies have well-known limitations with respect to invasiveness and sampling variability[44] and cannot be performed in very large-scale studies. Alternative, non-invasive strategies for the evaluation of NAFLD are serum biomarkers or the use of imaging techniques. However, imaging techniques are time consuming, expensive and also not feasible in large observational studies. Given these considerations, the recent EASL-EASD-EASO NAFLD guidelines have adopted that serum biomarkers are the preferred diagnostic tool for larger scale screening studies[29]. For the identification of participants with NAFLD in this study, the FLI was used, which was developed from data of the Dionysos Nutrition & Liver Study in Northern Italy[28]. The FLI is one of the three best-validated steatosis biomarkers in the new international accepted guideline[29], has a good steatosis predicting value (AUROC 0.83)[45], and is accurate in detecting NAFLD (accuracy of 0.84 and specificity of 86% for a FLI≥60)[28].

This study is unique in its cohort size of nearly 40,000 participants, which enabled careful calculations on effect sizes, sufficiently powered subgroup analysis and sufficient statistical power to investigate associations. All participants included in the Lifelines Cohort Study have been well examined, with extensive validated questionnaires, standardized anthropometric and laboratory measurements performed in serum samples in one certified laboratory with ditto equipment and quality assessment control for all samples[26]. In addition, included participants in this study had similar distributions of sex, age, BMI, T2DM and MetS compared with the whole Lifelines cohort, so results can be reflected to the total study population. Furthermore, the Lifelines study population has been previously validated, the risk of selection bias is low, is representative and can be generalized to the population of the North of the Netherlands[27].

Several methodological aspects and limitations also need to be addressed. First, this is a cross-sectional study. Thus cause-effect relationships cannot be established with certainty. Second, although the FLI score is an accepted diagnostic tool for NAFLD, it is not an absolute measure of hepatic fat accumulation and thus over- and underestimation of NAFLD could have occurred. Moreover, since the formula of the FLI contains the variables GGT, triglycerides, waist circumference and BMI, the associations of these variables with suspected NAFLD cannot be appropriately ascertained. Finally, since ancestry, alcohol intake, medication use and medical history were based on self-administered questionnaires, misreporting by individuals cannot be excluded. However, considering the large number of subjects, this limitation does not materially affect the interpretation of the presented results.

Conclusions

In this large study cohort of almost 40,000 subjects performed in the Northern part of the Netherlands, NAFLD is a major suspected health problem. NAFLD is suspected in 22% of a general European population and this group has an increased risk of having T2DM, MetS and a history of cardiovascular disease and impaired renal function. Future analysis of these subjects regarding the development of fibrosis and other population-based studies are mandatory to better understand the natural history of NAFLD and prevent and treat its complications.

Supporting Information

[Figure omitted. See PDF.]

S1 Table. Overview of studies on non-alcoholic fatty liver disease prevalence in Europe.

https://doi.org/10.1371/journal.pone.0171502.s001

(DOCX)

Acknowledgments

The authors wish to acknowledge the services of the Lifelines Cohort Study and Biobank, the contributing research centers delivering data, the participating general practitioners and pharmacists and all the study participants.

The manuscript is based on data from the Lifelines Cohort Study. Lifelines adheres to standards for open data availability. The data catalogue of Lifelines is publicly accessible at www.lifelines.net. All international researchers can apply for data at the Lifelines research office. The Lifelines system allows access for reproducibility of the study results.

Author Contributions

1. Conceptualization: EB MA TS RD KNF BA HB.

2. Data curation: EB MA BA.

3. Formal analysis: EB MA.

4. Investigation: EB.

5. Methodology: EB MA BA HB.

6. Supervision: HB.

7. Validation: MA BA.

8. Visualization: EB.

9. Writing – original draft: EB.

10. Writing – review & editing: EB MA TS RD KNF BA HB.

Citation: van den Berg EH, Amini M, Schreuder TCMA, Dullaart RPF, Faber KN, Alizadeh BZ, et al. (2017) Prevalence and determinants of non-alcoholic fatty liver disease in lifelines: A large Dutch population cohort. PLoS ONE 12(2): e0171502. https://doi.org/10.1371/journal.pone.0171502

References

1. Loomba R, Sanyal AJ. The global NAFLD epidemic. Nat Rev Gastroenterol Hepatol. Nature Publishing Group; 2013;10: 686–690.

2. Younossi ZM, Koenig AB, Abdelatif D, Fazel Y, Henry L, Wymer M. Global Epidemiology of Non-Alcoholic Fatty Liver Disease-Meta-Analytic Assessment of Prevalence, Incidence and Outcomes. Hepatology. 2015.

3. Wong RJ, Aguilar M, Cheung R, Perumpail RB, Harrison SA, Younossi ZM, et al. Nonalcoholic steatohepatitis is the second leading etiology of liver disease among adults awaiting liver transplantation in the United States. Gastroenterology. 2015;148: 547–555. pmid:25461851

4. Rafiq N, Bai C, Fang Y, Srishord M, McCullough A, Gramlich T, et al. Long-term follow-up of patients with nonalcoholic fatty liver. Clin Gastroenterol Hepatol. 2009;7: 234–238. pmid:19049831

5. Satapathy SK, Sanyal AJ. Epidemiology and Natural History of Nonalcoholic Fatty Liver Disease. Semin Liver Dis. 2015;35: 221–235. pmid:26378640

6. Marchesini G, Bugianesi E, Forlani G, Cerrelli F, Lenzi M, Manini R, et al. Nonalcoholic fatty liver, steatohepatitis, and the metabolic syndrome. Hepatology. 2003;37: 917–923. pmid:12668987

7. Fracanzani AL, Burdick L, Raselli S, Pedotti P, Grigore L, Santorelli G, et al. Carotid Artery Intima-media Thickness in Nonalcoholic Fatty Liver Disease. The American Journal of Medicine. 2008;121: 72–78. pmid:18187076

8. Gastaldelli A, Kozakova M, Højlund K, Flyvbjerg A, Favuzzi A, Mitrakou A, et al. Fatty liver is associated with insulin resistance, risk of coronary heart disease, and early atherosclerosis in a large European population. Hepatology. 2009;49: 1537–1544. pmid:19291789

9. Bedogni G, Miglioli L, Masutti F, Tiribelli C, Marchesini G, Bellentani S. Prevalence of and risk factors for nonalcoholic fatty liver disease: the Dionysos nutrition and liver study. Hepatology. 2005;42: 44–52. pmid:15895401

10. Bellentani S, Saccoccio G, Masutti F, Crocè LS, Brandi G, Sasso F, et al. Prevalence of and risk factors for hepatic steatosis in Northern Italy. Ann Intern Med. 2000;132: 112–117. pmid:10644271

11. Zois C-D, Baltayiannis G-H, Bekiari A, Goussia A, Karayiannis P, Doukas M, et al. Steatosis and steatohepatitis in postmortem material from Northwestern Greece. World J Gastroenterol. 2010;16: 3944–3949. pmid:20712056

12. Armstrong MJ, Houlihan DD, Bentham L, Shaw JC, Cramb R, Olliff S, et al. Presence and severity of non-alcoholic fatty liver disease in a large prospective primary care cohort. J Hepatol. 2012;56: 234–240. pmid:21703178

13. Caballería L, Pera G, Auladell MA, Torán P, Muñoz L, Miranda D, et al. Prevalence and factors associated with the presence of nonalcoholic fatty liver disease in an adult population in Spain. Eur J Gastroenterol Hepatol. 2010;22: 24–32. pmid:19730384

14. Kanerva N, Sandboge S, Kaartinen NE, Männistö S, Eriksson JG. Higher fructose intake is inversely associated with risk of nonalcoholic fatty liver disease in older Finnish adults. Am J Clin Nutr. 2014;100: 1133–1138. pmid:25099548

15. Loguercio C, De Girolamo V, de Sio I, Tuccillo C, Ascione A, Baldi F, et al. Non-alcoholic fatty liver disease in an area of southern Italy: main clinical, histological, and pathophysiological aspects. J Hepatol. 2001;35: 568–574. pmid:11690701

16. Radu C, Grigorescu M, Crisan D, Lupsor M, Constantin D, Dina L. Prevalence and associated risk factors of non-alcoholic fatty liver disease in hospitalized patients. J Gastrointestin Liver Dis. 2008;17: 255–260. pmid:18836616

17. Suomela E, Oikonen M, Virtanen J, Parkkola R, Jokinen E, Laitinen T, et al. Prevalence and determinants of fatty liver in normal-weight and overweight young adults. The Cardiovascular Risk in Young Finns Study. Ann Med. 2015;47: 40–46. pmid:25333756

18. Tarnoki AD, Tarnoki DL, Bata P, Littvay L, Osztovits J, Jermendy G, et al. Heritability of non-alcoholic fatty liver disease and association with abnormal vascular parameters: a twin study. Liver Int. 2012;32: 1287–1293. pmid:22651705

19. van der Voort EAM, Koehler EM, Dowlatshahi EA, Hofman A, Stricker BH, Janssen HLA, et al. Psoriasis is independently associated with nonalcoholic fatty liver disease in patients 55 years old or older: Results from a population-based study. J Am Acad Dermatol. 2014;70: 517–524. pmid:24373781

20. Völzke H, Robinson D-M, Kleine V, Deutscher R, Hoffmann W, Ludemann J, et al. Hepatic steatosis is associated with an increased risk of carotid atherosclerosis. World J Gastroenterol. 2005;11: 1848–1853. pmid:15793879

21. Ludwig U, Holzner D, Denzer C, Greinert A, Haenle MM, Oeztuerk S, et al. Subclinical and clinical hypothyroidism and non-alcoholic fatty liver disease: a cross-sectional study of a random population sample aged 18 to 65 years. BMC Endocr Disord. 2015;15: 41. pmid:26276551

22. Kotronen A, Yki-Järvinen H, Männistö S, Saarikoski L, Korpi-Hyövälti E, Oksa H, et al. Non-alcoholic and alcoholic fatty liver disease—two diseases of affluence associated with the metabolic syndrome and type 2 diabetes: the FIN-D2D survey. BMC Public Health. 2010;10: 237. pmid:20459722

23. Papatheodoridis GV, Goulis J, Christodoulou D, Manolakopoulos S, Raptopoulou M, Andrioti E, et al. High prevalence of elevated liver enzymes in blood donors: associations with male gender and central adiposity. Eur J Gastroenterol Hepatol. 2007;19: 281–287. pmid:17353691

24. Pendino GM, Mariano A, Surace P, Caserta CA, Fiorillo MT, Amante A, et al. Prevalence and etiology of altered liver tests: a population-based survey in a Mediterranean town. Hepatology. 2005;41: 1151–1159. pmid:15841464

25. Stolk RP, Rosmalen JGM, Postma DS, de Boer RA, Navis G, Slaets JPJ, et al. Universal risk factors for multifactorial diseases: LifeLines: a three-generation population-based study. Eur J Epidemiol. 2008;23: 67–74. pmid:18075776

26. Scholtens S, Smidt N, Swertz MA, Bakker SJL, Dotinga A, Vonk JM, et al. Cohort Profile: LifeLines, a three-generation cohort study and biobank. Int J Epidemiol. 2015;44: 1172–1180. pmid:25502107

27. Klijs B, Scholtens S, Mandemakers JJ, Snieder H, Stolk RP, Smidt N. Representativeness of the LifeLines Cohort Study. PLoS ONE. 2015;10: e0137203. pmid:26333164

28. 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. pmid:17081293

29. European Association for the Study of the Liver (EASL). Electronic address: [email protected], European Association for the Study of Diabetes (EASD), European Association for the Study of Obesity (EASO). EASL-EASD-EASO Clinical Practice Guidelines for the management of non-alcoholic fatty liver disease. J Hepatol. 2016.

30. Grundy SM, Cleeman JI, Daniels SR, Donato KA, Eckel RH, Franklin BA, et al. Diagnosis and management of the metabolic syndrome: an American Heart Association/National Heart, Lung, and Blood Institute Scientific Statement. Circulation. 2005;112: 2735–2752. pmid:16157765

31. Levey AS, Coresh J, Greene T, Stevens LA, Zhang YL, Hendriksen S, et al. Using standardized serum creatinine values in the modification of diet in renal disease study equation for estimating glomerular filtration rate. Ann Intern Med. 2006;145: 247–254. pmid:16908915

32. Levey AS, Coresh J, Greene T, Marsh J, Stevens LA, Kusek JW, et al. Expressing the Modification of Diet in Renal Disease Study equation for estimating glomerular filtration rate with standardized serum creatinine values. Clin Chem. 2007;53: 766–772. pmid:17332152

33. Browning JD, Szczepaniak LS, Dobbins R, Nuremberg P, Horton JD, Cohen JC, et al. Prevalence of hepatic steatosis in an urban population in the United States: impact of ethnicity. Hepatology. 2004;40: 1387–1395. pmid:15565570

34. Verma S, Jensen D, Hart J, Mohanty SR. Predictive value of ALT levels for non-alcoholic steatohepatitis (NASH) and advanced fibrosis in non-alcoholic fatty liver disease (NAFLD). Liver Int. 2013;33: 1398–1405. pmid:23763360

35. Chiolero A, Faeh D, Paccaud F, Cornuz J. Consequences of smoking for body weight, body fat distribution, and insulin resistance. Am J Clin Nutr. 2008;87: 801–809. pmid:18400700

36. Yuan H, Shyy JY-J, Martins-Green M. Second-hand smoke stimulates lipid accumulation in the liver by modulating AMPK and SREBP-1. J Hepatol. 2009;51: 535–547. pmid:19556020

37. Trak-Smayra V, Dargere D, Noun R, Albuquerque M, Yaghi C, Gannagé-Yared M-H, et al. Serum proteomic profiling of obese patients: correlation with liver pathology and evolution after bariatric surgery. Gut. 2009;58: 825–832. pmid:18403495

38. Yu C, Xu C, Xu L, Yu J, Miao M, Li Y. Serum proteomic analysis revealed diagnostic value of hemoglobin for nonalcoholic fatty liver disease. J Hepatol. 2012;56: 241–247. pmid:21756851

39. Yilmaz Y, Senates E, Ayyildiz T, Colak Y, Tuncer I, Ovunc AOK, et al. Characterization of nonalcoholic fatty liver disease unrelated to the metabolic syndrome. Eur J Clin Invest. 2012;42: 411–418. pmid:21913918

40. Foroughi M, Maghsoudi Z, Khayyatzadeh S, Ghiasvand R, Askari G, Iraj B. Relationship between non-alcoholic fatty liver disease and inflammation in patients with non-alcoholic fatty liver. Adv Biomed Res. 2016;5: 28. pmid:27014655

41. Lyon CJ, Law RE, Hsueh WA. Minireview: adiposity, inflammation, and atherogenesis. Endocrinology. 2003;144: 2195–2200. pmid:12746274

42. Dullaart RPF, Gruppen EG, Connelly MA, Lefrandt JD. A pro-inflammatory glycoprotein biomarker is associated with lower bilirubin in metabolic syndrome. Clin Biochem. 2015;48: 1045–1047. pmid:26129880

43. Day CP. Non-alcoholic fatty liver disease: current concepts and management strategies. Clin Med (Lond). 2006;6: 19–25.

44. Rockey DC, Caldwell SH, Goodman ZD, Nelson RC, Smith AD, American Association for the Study of Liver Diseases. Liver biopsy. Hepatology (Baltimore, Md.). 2009. pp. 1017–1044.

45. Fedchuk L, Nascimbeni F, Pais R, Charlotte F, Housset C, Ratziu V, et al. Performance and limitations of steatosis biomarkers in patients with nonalcoholic fatty liver disease. Aliment Pharmacol Ther. 2014;40: 1209–1222. pmid:25267215