Prevalence and risk factors of nonalcoholic fatty liver disease in patients with inflammatory bowel diseases: A cross-sectional and longitudinal analysis

Abstract

BACKGROUND

Nonalcoholic fatty liver disease (NAFLD) is common in the German population, with an even higher prevalence in inflammatory bowel disease patients.

AIM

To investigate the risk factors for NAFLD in inflammatory bowel disease patients.

METHODS

This monocentric retrospective study with a cross-sectional and a longitudinal part included 694 patients. Inclusion criteria were diagnosed inflammatory bowel disease, age ≥ 18 years, availability of at least one abdominal ultrasound. Patients with infectious or suspected alcoholic fatty liver disease were excluded. NAFLD was defined by increased echogenicity at liver ultrasound. Demographic characteristics, disease activity and medications were analyzed as potential risk factors. Parameters influencing the course of NAFLD were identified by a generalized linear mixed model.

RESULTS

Forty-eight percent of Crohn’s disease (CD) patients and 44% of ulcerative colitis patients suffered from NAFLD. Its occurrence was associated with greater age, hypertension and body mass index (BMI) in both groups, and with higher disease activity and dyslipidemia in CD. 2467 ultrasound results were included in the longitudinal analysis. Risk factors for NAFLD were age, BMI, higher disease activity, bowel resection(s), endoscopic activity and azathioprine use in CD; and BMI and endoscopic activity in ulcerative colitis.

CONCLUSION

NAFLD was highly prevalent in this cohort of German inflammatory bowel disease patients. Its risk increased mainly with rising age and BMI. This analysis provides a rationale for non-invasive liver screening in inflammatory bowel disease patients.

Key Words: Crohn’s disease; Ulcerative colitis; Inflammatory bowel disease; Nonalcoholic fatty liver disease; Nonalcoholic steatohepatitis; Ultrasound

Core Tip: Nonalcoholic fatty liver disease (NAFLD) has a high prevalence. Aim is to investigate risk factors for NAFLD in inflammatory bowel disease. 48% of Crohn’s disease (CD) and 44% of ulcerative colitis patients suffered from NAFLD. Its occurrence was associated with greater age and body mass index (BMI) in both groups, and with higher disease activity in CD. Risk factors for NAFLD in the longitudinal analysis were age, BMI, higher disease activity, bowel resection(s), endoscopic activity and azathioprine use in CD; and BMI and endoscopic activity in ulcerative colitis.

INTRODUCTION

Nonalcoholic fatty liver disease (NAFLD) is a common disorder with a prevalence of 25% worldwide[1] and 22.9% in the German population. Its prevalence in Germany is expected to reach 26.4% by the year 2030[2]. In patients with inflammatory bowel diseases (IBDs) – consisting of Crohn’s disease (CD) and ulcerative colitis (UC) – an even higher NAFLD prevalence of up to 40% has been described[3-5].

The main risk factor for NAFLD is the metabolic syndrome[6]. Ultrasound of the liver is the most common diagnostic tool to detect NAFLD with a sensitivity of 85% and specificity of 94%[7]. Hepatic steatosis may progress to nonalcoholic steatohepatitis (NASH), which can result in the development of cirrhosis and hepatocellular carcinoma. NAFLD is associated with an increased overall mortality and independently predicts the risk for future cardiovascular disease events[8,9].

IBD patients develop NAFLD based on fewer metabolic risk factors in comparison with non-IBD NAFLD patients[10]. In IBD patients, hypertension, obesity, small bowel surgery, and use of corticosteroids at the time of imaging were associated with NAFLD, while among patients receiving anti-TNFα therapy, NAFLD has been shown to be less prevalent[10]. In contrast, another study found a high prevalence of NAFLD among IBD patients under anti-TNFα treatment[11].

Thus it may be assumed that the risk factors for NAFLD in IBD patients differ from those in the general population. Even underweight IBD patients showed an increased prevalence of NAFLD[12]. A history of bowel resection(s) and parenteral nutrition have also been described as risk factors for NAFLD[13,14]. Treatment with high doses of corticosteroids–as still frequently applied in IBD patients–has been related to hepatic steatosis as well[15].

The key goal of this study was to investigate the prevalence and risk factors of NAFLD in patients suffering from IBD. The superordinate aim of this project is to optimize preventive and therapeutic measures of NAFLD in the identified population at risk, and thus to reduce the incidence of NASH or hepatocellular carcinoma in IBD patients.

MATERIALS AND METHODS

Study design and data extraction

This is an uncontrolled, monocentric retrospective study including patients diagnosed with IBD at a South German tertiary referral center for the care of IBD. The Ethics Committee of the University of Heidelberg (protocol number: S-299/2018) approved the study. Inclusion criteria were as follows: (1) Age ≥ 18 years; (2) Diagnosis of IBD according to European Crohn’s and Colitis Organisation criteria[16]; (3) Therapy at the IBD outpatient clinic of the University Hospital Heidelberg; (4) At minimum one visit to the IBD outpatient department between June 1, 2014 and May 31, 2018; and (5) Availability of at least one documented abdominal ultrasound examination including profound assessment of the liver performed prior to May 31, 2018, as this date was defined as the cut-off time point for data acquisition. Exclusion criteria were: (1) Diagnosis of any infectious liver disease; (2) Suspected alcoholic fatty liver disease; or (3) Documented abuse of alcohol (daily intake of alcohol of > 10 g/d in women and > 20 g/d in men). All data were collected from computerized medical course records. The available daily consultation lists of the IBD outpatient clinic between June 1, 2014 and May 31, 2018 were screened, to identify eligible patients.

All data concerning NALFD diagnosis and grading of NAFLD severity were extracted from the liver ultrasound reports of the patients. NAFLD severity was categorized into four groups: (1) No steatosis; (2) Mild steatosis; (3) Moderate steatosis; and (4) Severe steatosis, as specified in the ultrasound report.

Demographic and disease-related parameters of all included patients were entered into a Microsoft Excel spreadsheet. As standard at the IBD outpatient department, all IBD patients who were under immunosuppressive treatment underwent a routine abdominal ultrasound every year. Patients who visited the clinic for the sole purpose of a second opinion usually did not receive an abdominal ultrasound examination. In cases of diagnosed NALFD, the patients were routinely advised about lifestyle modifications, including nutrition and physical training.

Definitions

The Montreal classification for CD and UC[17] was used to categorize disease extent. Disease activity was graded using the Harvey Bradshaw Index (HBI)[18] for CD patients and the Simple Clinical Colitis Activity Index (SCCAI)[19] for UC patients. A HBI and SCCAI of > 4, respectively, were considered to represent active inflammation. HBI and SCCAI scores have been routinely used at our outpatient clinic since 2017 to monitor disease activity and evaluate therapeutic effects. For data collected prior to 2017, the indices were retrospectively calculated on the basis of all written information available in the patients’ records.

Studied variables included age at diagnosis of IBD, age at index abdominal ultrasound, body mass index (BMI), sex, presence of diabetes mellitus, dyslipidemia or hypertension, smoking status, CD location and phenotype, extent of UC, laboratory results [plasma concentrations of aspartate aminotransferase, alanine aminotransferase, gamma-glutamyltransferase (γGT), C-reactive protein, white blood cell count, fecal calprotectin concentration], HBI, SCCAI, inflammation detected by intestinal ultrasound or colonoscopy, medications, history of intestinal surgery, parenteral nutrition, concomitant diseases, and infections.

A sonographically measured bowel wall thickness of > 3 mm in the small bowel and of > 4 mm in the colon were considered to represent inflammation. Endoscopically active disease was considered whenever the mucosa was not described as normal.

Liver ultrasound

In the cross-sectional component of this study, ultrasound results from May 31, 2017 through May 31, 2018 were included. In the longitudinal part of the study, all results of any ultrasound examination ever performed at our center in the included patients were considered.

Study endpoints

The primary study end points were (1) To determine the prevalence of NAFLD in the cross-sectional part of the study; (2) To reveal factors associated with changes of the presence and severity of NAFLD during the further disease course in the longitudinal part of the study.

The secondary aim of the study was to identify potential risk factors for NAFLD in the cross-sectional component of the study.

Statistical analyses

All statistical analyses were performed using “R” Version 3.5.2 (2018-12-20). Normality of data distribution was first assessed using the Shapiro-Wilk test. None of the data were normally distributed, so the results were calculated as medians and interquartile ranges.

To identify potential risk factors for NAFLD, the Mann-Whitney-U-test and Kruskal-Wallis-test were used for ordinal and continuous variables, and X²-tests for categorical variables. For the multivariable assessment of risk factors for NAFLD, a logistic regression analysis was performed. A backwards elimination was used with the following parameters: BMI, age, sex, diabetes mellitus, hypertension, current smoking, dyslipidemia, activity index showing inflammation (HBI and SCCAI), elevated C-reactive protein level, inflammation in colonoscopy, previous bowel resection(s), as well as therapy with steroids, azathioprine, infliximab, adalimumab, vedolizumab, ustekinumab, mesalazine or sulfasalazine. A linear mixed model for repeated measurements was applied for the identification of potential risk factors which influence the changes in NAFLD: The change of NAFLD status was the dependent variable, while the previously identified and known confounders of BMI, age, sex, diabetes mellitus, hypertension, dyslipidemia, current smoking, and use of steroids were fixed factors. Parameter estimates were determined together with 95% confidence intervals and associated P values.

In the linear mixed model for repeated measurements, all ultrasound examinations were included. The model considered the different numbers of ultrasound examinations per patient as well as changes in clinical, laboratory and treatment parameters and their impact on NAFLD. Because of the exploratory character of the study, there was no adjustment for multiple testing. The interpretation of P values is descriptive. P values < 0.05 were considered as statistically significant.

RESULTS

Characteristics of Crohn’s disease patients

651 different CD patients attended the IBD outpatient clinic between June 1, 2014 and May 31, 2018. Among these, 184 patients did not undergo any ultrasound examination, which was mainly due to their receiving a second opinion for therapy of CD. Exclusion of 22 patients was due to the diagnosis of other liver diseases, and five patients due to alcohol abuse. Thus, 440 CD patients who underwent at least one abdominal ultrasound examination were included. Among these 440 patients, 138 patients were excluded from the cross-sectional part of the study because they did not undergo any ultrasound examination between May 31, 2017 and May 31, 2018, leaving 302 patients included in the cross-sectional part of the study. In the longitudinal part of the study, 1620 abdominal ultrasound results from 440 patients were included. Figure 1A shows a flow diagram of all patient inclusions and exclusions.

Figure 1 Flow diagram of inclusion and exclusion. A: Crohn’s disease patients; B: Ulcerative colitis patients.

Forty-nine percent of the included CD patients in the cross-sectional study were male, the median age at the time of the ultrasound examination was 42 years, and the median BMI was 24.5 kg/m².

Results of the cross-sectional study in CD patients

In the cross-sectional study, fifty-two percent of the CD patients did not suffer from NAFLD. Mild NAFLD was revealed in 18%, moderate NAFLD in 16%, and severe NAFLD in 14% of the patients.

Comparison of demographic and disease-related parameters between the subgroups with and without a diagnosis of NAFLD determined by the ultrasound results is shown in Table 1.

Table 1 Comparison of baseline characteristics between the subgroups of patients with nonalcoholic fatty liver disease versus those without nonalcoholic fatty liver disease in Crohn’s disease.

Variable	NAFLD (n = 145)	No NAFLD (n = 157)	P value
Male, n (%)	75 (52)	74 (47)	0.495
Age at abdominal ultra sound (yr), median (IQR)	49.0 (20.0)	35.0 (19.0)	< 0.001c
BMI (kg/m²), median (IQR)	28.0 (7.6)	23.0 (4.4)	< 0.001c
Diabetes mellitus, n (%)	7 (5)	1 (1)	0.057
Dyslipidemia, n (%)	12 (8)	3 (2)	0.0131
Hypertension, n (%)	50 (34)	20 (13)	< 0.001c
Current smoking, n (%)	35 (24)	29 (18)	0.288
Smoking in the past, n (%)	28 (19)	31 (20)	1.000
Smoking in general, n (%)	63 (44)	60 (38)	0.420
HBI (IQR), n = 217	3.0 (5.0)	1.0 (3.0)	< 0.001c
HBI > 4, n (%), n = 217	46 (40)	15 (15)	< 0.001c
Inflammation in colonoscopy, n (%), n = 83	27 (61)	20 (51)	0.330
Inflammation in bowel ultrasound, n (%), n = 101	54 (82)	47 (85)	0.772
CRP, mg/L (IQR), n = 293	2.5 (7.0)	3.0 (7.6)	0.846
White blood cells, /nl (IQR), n = 293	7.8 (3.0)	7.8 (2.8)	0.737
Fecal calprotectin, µg/g (IQR), n = 172	111 (230)	112 (249)	0.638
γGT, U/l (IQR)	24.0 (25.0)	16 .0(16.5)	< 0.001c
Elevated transaminases, n (%), n = 260	18 (16)	9 (6)	0.069
Disease duration at abdominal ultrasound (yr), median (IQR)	16.0 (21.0)	11.0 (12.0)	0.002b
History of bowel resection(s), n (%)	72 (50)	54 (45)	0.010a
Biological treatment at the time of ultrasound, n (%)	104 (72)	103 (66)	0.308
Treatment with adalimumab at the time of ultrasound, n (%)	53 (37)	57 (36)	0.767
Treatment with infliximab at the time of ultrasound, n (%)	29 (20)	21 (13)	0.164
Treatment with vedolizumab at the time of ultrasound, n (%)	11 (8)	10 (6)	0.850
Treatment with ustekinumab at the time of ultrasound, n (%)	14 (10)	18 (11)	0.746
Treatment with azathioprine at the time of ultrasound, n (%)	19 (13)	18 (11)	0.796
Treatment with budesonide at the time of ultrasound, n (%)	7 (5)	12 (8)	0.576
Treatment with systemic steroids < 7.5 mg at the time of ultrasound, n (%)	6 (4)	8 (5)	0.903
Treatment with systemic steroids > 7.5 mg at the time of ultrasound, n (%)	10 (7)	15 (10)	0.530
Treatment with systemic steroids at the time of ultrasound, n (%)	14 (10)	21 (13)	0.407
Treatment with mesalazine at the time of ultrasound, n (%)	25 (17)	19 (12)	0.271
Treatment with sulfasalazine at the time of ultrasound, n (%)	13 (9)	3 (2)	0.013a

^aP < 0.05.

^bP < 0.01.

^cP < 0.001. BMI: Body mass index; CRP: C-reactive protein; γGT: Gamma-glutamyl transferase; HBI: Harvey bradshaw index; IQR: Interquartile range; NAFLD: Nonalcoholic fatty liver disease.

The BMI differed significantly between the subgroups: 23.0 kg/m² was the median BMI of patients without NAFLD as compared to 28.0 kg/m² in patients with NAFLD (P < 0.001). Age at abdominal ultrasound was 35 years in patients without NAFLD as compared to 49 years in patients with NAFLD (P < 0.001). An increased HBI > 4 was found in 15% of the patients without NAFLD, compared to 40% with NAFLD (P < 0.001). Other significant differences were revealed in dyslipidemia, hypertension, γGT, disease duration, history of bowel resection(s), and treatment with sulfasalazine. Laboratory results and endoscopic or ultrasound results suggesting active inflammation were not associated with the presence of NAFLD, while elevated plasma liver enzyme concentrations tended to be more frequent in NAFLD patients (P = 0.069).

Using a multivariable logistic regression model, age at abdominal ultrasound, increased HBI, and increased BMI were identified as independent risk factors for NAFLD in CD patients (Table 2).

Table 2 Results from multivariable logistic regression model (backwards elimination) to identify risk factors for nonalcoholic fatty liver disease in Crohn’s disease patients.

Parameter	Odds ratio	95%CI	P value
Intercept	< 0.01	< 0.01; < 0.01	< 0.001c
Elevated HBI	3.80	1.75-8.70	0.001b
Age at abdominal ultrasound (yr)	1.05	1.02-1.08	< 0.001c
BMI (kg/m²)	1.25	1.15-1.37	< 0.001c
History of bowel resection(s)	1.47	0.74-2.92	0.274
Hypertension	0.81	0.31-2.09	0.668
Current smoking	1.17	0.51-2.69	0.707
Dyslipidemia	3.48	0.52-39.13	0.245

^bP < 0.01.

^cP < 0.001. BMI: Body mass index; CRP: C-reactive protein; HBI: Harvey Bradshaw Index; OR: Odds ratio; CI: Confidence interval.

Using a multivariable logistic regression model considering BMI and disease duration, disease duration was also significantly associated with NAFLD (P < 0.001). These results were evaluated in a separate model, as age and disease duration are too similar to be included in the same model.

Time-dependent course of NAFLD and relation to CD therapy

In the 440 included CD patients, 1620 abdominal ultrasound examinations were performed (Figure 1A). The minimum number of examinations per patient was one, while the maximum number of examinations per patient was 12 (Table 3). Of 1180 time intervals, a NAFLD improvement was found in 141 time intervals (11.9%), while NAFLD deteriorated in 159 intervals (13.5%), 880 time intervals showed no change (74.6%). The mean length of a time interval was three years and seven months. Using a linear mixed model for repeated measurements, the parameters at the time intervals between the ultrasound examinations were inspected. In Table 4, the results of the linear mixed model for repeated measurements analyzing the NAFLD status are shown.

Table 3 Number of abdominal ultrasound examinations per patient.

Ultrasound examinations (n)	CD patients	UC patients
1	119	83
2	62	50
3	57	31
4	53	21
5	48	14
6	37	9
7	24	14
8	17	17
9	14	7
10	4	13
11	2	3
12	2	1

CD: Crohn’s disease; UC: Ulcerative colitis.

Table 4 Results of mixed model analysis to identify factors associated with nonalcoholic fatty liver disease in Crohn’s disease patients.

Confounder model
Variable	Estimate	95%CI	P value
Female sex	0.70	0.09-1.31	0.293
Age at abdominal ultrasound (in years)	1.07	1.04-1.09	< 0.001c
BMI at abdominal ultrasound (in years)	1.32	1.30-1.34	< 0.001c
Diabetes mellitus	2.36	0.11-4.61	0.509
Hypertension	0.78	-0.07-1.62	0.586
Dyslipidemia	10.82	8.85-12.80	0.049a
Current smoking	1.91	1.28-2.53	0.109
Steroid use	0.74	0.27-1.21	0.287
Parameters included in model with age and BMI
Elevated HBI	2.99	2.41-3.58	< 0.001c
Inflammation in colonoscopy	16.40	14.59-18.11	0.002b
Inflammation in bowel ultrasound	1.46	1.30-1.61	0.344
Elevated CRP concentration	1.17	1.01-1.33	0.502
Elevated white blood cell count	0.96	0.10-1.12	0.295
Elevated fecal calprotectin concentration	1.34	1.19-1.50	0.455
History of bowel resection(s)	2.03	1.42-2.64	0.023a
Biologic treatment	1.01	0.85-1.17	0.959
Adalimumab treatment	1.09	0.93-1.25	0.709
Infliximab treatment	1.64	1.48-1.79	0.135
Vedolizumab treatment	1.40	1.24-1.56	0.523
Ustekinumab treatment	1.13	0.01-2.25	0.833
Golimumab treatment	1.34	1.19-1.50	0.874
Azathioprine therapy	1.98	1.36-2.61	0.031a
Methotrexate treatment	1.15	0.99-1.30	0.843
6-Mercaptopurine treatment	1.08	0.92-1.24	0.940
Sulfasalazine treatment	0.53	0.38-0.69	0.310
Mesalazine treatment	1.14	0.98-1.29	0.688

^aP < 0.05.

^bP < 0.01.

^cP < 0.001. BMI: Body mass index; CRP: C-reactive protein; HBI: Harvey bradshaw index; CI: Confidence interval.

First, potential confounders as BMI, age, sex, diabetes mellitus, hypertension, dyslipidemia, current smoking, and steroid use were included in the model. BMI [Odds Ratio (OR): 1.32, P < 0.001] age (OR: 1.07, P < 0.001) and dyslipidemia (OR: 10.82, P = 0.049) were statistically associated with progression of NAFLD. In a second step, an additional linear mixed model including all parameters plus BMI and age was calculated. Increased HBI, inflammation in colonoscopy, history of bowel resection(s), and azathioprine use were significantly associated with NAFLD progression.

Characteristics of ulcerative colitis patients

400 different UC patients attended the IBD outpatient department between June 1, 2014 and May 31, 2018. Of these, 128 patients suffering from UC did not undergo any ultrasound examination at our clinic mainly due to the fact that they came to receive a second opinion. Sixteen patients were excluded due to other liver diseases, and two patients due to alcohol abuse. Thus, 254 patients who underwent more than one abdominal ultrasound examination at the IBD outpatient clinic were included in the study. Of these 254 patients, 101 patients were excluded from the cross-sectional part of the study because they did not undergo any ultrasound examination between May 31, 2017 and May 31, 2018. Thus, 153 patients were included in the cross-sectional part of the study. In the longitudinal part of the study, 847 abdominal ultrasound results from 254 patients were considered. Figure 1B shows a flow diagram of inclusions and exclusions for all UC patients.

Fifty percent of the UC patients included in the cross-sectional part of the study were male, the median age at the time point of ultrasound was 40 years, while the median BMI was 25.4 kg/m².

Results of the cross-sectional study in UC patients

In the cross-sectional study, fifty-six percent of the included UC patients showed no signs of NAFLD, while mild NAFLD was diagnosed in 20%, moderate NAFLD in 16%, and severe NAFLD in 8%.

Demographic and clinical parameters were compared between the patients who were diagnosed with NAFLD and those who were not. Table 5 demonstrates the differences of demographic and clinical parameters of patients with NAFLD compared to patients without NAFLD.

Table 5 Comparison of baseline characteristics between the subgroups of patients with nonalcoholic fatty liver disease versus those without nonalcoholic fatty liver disease in ulcerative colitis.

Variable	NAFLD (n = 68)	No NAFLD (n = 85)	P value
Male, n (%)	40 (59)	37 (44)	0.086
Age at abdominal ultrasound (yr), median (IQR)	48.5 (18.5)	37.0 (21.0)	0.005b
BMI (kg/m²), median (IQR)	28.0 (5.7)	22.8 (5.3)	< 0.001c
Diabetes mellitus, n (%)	6 (9)	1 (1)	0.057
Dyslipidemia, n (%)	2 (3)	2 (2)	1.000
Hypertension, n (%)	22 (32)	6 (7)	< 0.001c
Current smoking, n (%)	4 (6)	1 (1)	0.242
Smoking in the past, n (%)	20 (30)	11 (13)	0.021a
Smoking in general, n (%)	24 (35)	12 (14)	0.004b
SCCAI (IQR), n = 124	2.0 (4.0)	1.0 (3.0)	0.253
SCCAI > 4, n (%), n = 124	14 (24)	9 (14)	0.156
Inflammation in colonoscopy, n (%), n = 45	15 (60)	9 (45)	0.123
Inflammation in bowel ultrasound, n (%), n = 28	10 (91)	15 (88)	0.812
CRP, mg/l (IQR), n = 146	2.9 (8.0)	1.0 (5.1)	0.114
White blood cells, /nl (IQR), n = 147	7.5 (2.9)	7.0 (3.2)	0.300
Fecal calprotectin, µg/g (IQR), n = 87	103 (390)	147 (359)	0.612
γGT, U/l (IQR)	23.0 (40.0)	18.0 (25.0)	0.007b
Elevated transaminases, n (%), n = 147	8 (13)	7 (8)	0.387
Disease duration at abdominal ultrasound (yr), median (IQR)	9.0 (12.5)	10.0 (30.5)	0.408
History of bowel resection(s), n (%)	6 (9)	12 (14)	0.449
Biological treatment at the time of ultrasound, n (%)	33 (49)	52 (61)	0.161
Treatment with adalimumab at the time of ultrasound, n (%)	8 (12)	13 (15)	0.694
Treatment with infliximab at the time of ultrasound, n (%)	7 (10)	20 (24)	0.055
Treatment with vedolizumab at the time of ultrasound, n (%)	19 (28)	19 (22)	0.544
Treatment with azathioprine at the time of ultrasound, n (%)	13 (19)	14 (16)	0.831
Treatment with budesonide at the time of ultrasound, n (%)	3 (4)	3 (4)	1.000
Treatment with systemic steroids < 7.5 mg at the time of ultrasound, n (%)	7 (10)	5 (6)	0.480
Treatment with systemic steroids > 7.5 mg at the time of ultrasound, n (%)	10 (15)	12 (14)	1.000
Treatment with systemic steroids at the time of ultrasound, n (%)	13 (19)	16 (19)	1.000
Treatment with mesalazine at the time of ultrasound, n (%)	42 (62)	42 (49)	0.173
Treatment with sulfasalazine at the time of ultrasound, n (%)	2 (3)	5 (6)	0.634

^aP < 0.05.

^bP < 0.01.

^cP < 0.001. BMI: Body mass index; CRP: C-reactive protein; γGT: Gamma-glutamyl transferase; IQR: Interquartile range; NAFLD: Nonalcoholic fatty liver disease; SCCAI: Simple Clinical Colitis Activity Index.

The BMI diverged significantly between patients with NAFLD and those without: patients without NAFLD had a median BMI of 22.8 kg/m² as compared to 28.0 kg/m² in patients with NAFLD (P < 0.001). Age at abdominal ultrasound was 37.0 years in patients without NAFLD as compared to 48.5 years in patients with NAFLD (P = 0.005). In patients with NAFLD, hypertension and smoking in genereal were more common compared to patients without NAFLD (P < 0.001). An SCCAI > 4 was found in 14% of the patients without NAFLD as compared to 24% in those with NAFLD (P = 0.156). Another significant difference was revealed in plasma γGT concentrations, which were 23 U/L in NAFLD patients and 18 U/L in patients without NALFD (P = 0.007). There was no association of laboratory results as well as endoscopic or ultrasound results proving inflammation with the occurrence of NAFLD.

Making use of a multivariable logistic regression model, increased BMI was identified as the only independent risk factor for NAFLD in UC patients (Table 6).

Table 6 Results from multivariable logistic regression model (backwards elimination) to identify risk factors for nonalcoholic fatty liver disease in ulcerative colitis patients.

Parameter	Odds ratio	95%CI	P value
Intercept	< 0.01	< 0.01; < 0.01	< 0.001c
Age at abdominal ultrasound (yr)	1.01	0.98-1.04	0.689
BMI (kg/m²)	1.25	1.14-1.38	< 0.001c
Hypertension	3.01	0.99-10.04	0.058
Adalimumab therapy	2.10	0.68-7.20	0.213

^cP < 0.001. BMI: Body mass index; OR: Odds ratio; CI: Confidence interval.

A multivariable logistic regression model considering BMI and disease duration did not identify disease duration as a risk factor of NAFLD in UC patients (P = 0.066).

Time-dependent course of NAFLD and relation to UC therapy

In the cohort of the 253 UC patents, 847 ultrasound examinations were performed in total (Figure 1B). The minimum number of examinations per patient was one, while the maximum number of examinations per patient was 12 (Table 3). In 593 time intervals, NAFLD improved over 77 time intervals (13.0%), while it worsened over 72 intervals (12.1%), and 444 time intervals showed no change (74.9%). The mean length of a time interval was three years and nine months. A linear mixed model for repeated measurements analyzing NAFLD change over time showed the results summarized in Table 7.

Table 7 Results from multivariable logistic regression model (backwards elimination) to identify risk factors for nonalcoholic fatty liver disease in ulcerative colitis patients.

Confounder model
Variable	Estimate	95%CI	P value
Female sex	0.62	0.50-0.65	0.159
Age at ultrasound (in years)	0.99	0.98-1.03	0.637
BMI at ultrasound (in years)	1.28	0.71-1.85	< 0.001c
Diabetes mellitus	2.14	0.95-4.79	0.441
Hypertension	2.89	2.51-4.10	0.037a
Dyslipidemia	1.46	0.98-2.84	0.695
Current smoking	1.77	1.26-2.40	0.525
Steroid use	0.88	0.20-1.54	0.652
Parameters included in model with age and BMI
Elevated SCCAI	1.58	1.42-1.74	0.145
Inflammation in colonoscopy	3.60	2.69-4.50	0.005b
Inflammation in bowel ultrasound	2.31	2.16-2.47	0.162
Elevated CRP concentration	0.99	0.83-1.14	0.959
Elevated white blood cell count	1.01	0.85-1.17	0.782
Elevated fecal calprotectin concentration	1.06	0.90-1.22	0.923
History of bowel resection(s)	1.66	1.50-1.82	0.332
Biologic treatment	0.69	0.53-0.85	0.167
Adalimumab treatment	0.74	0.58-0.90	0.426
Infliximab treatment	0.77	0.62-0.93	0.484
Vedolizumab treatment	1.14	0.98-1.30	0.739
Ustekinumab treatment	0.19	0.03-0.35	0.407
Golimumab treatment	1.14	0.98-1.30	0.892
Azathioprine therapy	1.19	1.03-1.35	0.597
Methotrexate treatment	1.50	1.34-1.66	0.734
6-Mercaptopurine treatment	0.94	0.79-1.10	0.947
Sulfasalazine treatment	1.43	1.27-1.58	0.491
Mesalazine treatment	0.88	-0.25-2.00	0.635

^aP < 0.05.

^bP < 0.01.

^cP < 0.001. BMI: Body mass index; CRP: C-reactive protein; SCCAI: Simple Clinical Colitis Activity Index; CI: Confidence interval.

In a first step, the potential confounders BMI, age, sex, diabetes mellitus, hypertension, dyslipidemia, current smoking, and steroid use were included in the model. BMI and hypertension were statistically associated with progression of NAFLD (OR: 1.28, P < 0.001; OR: 2.89, P = 0.037), while age was not (OR: 1.01, P = 0.667). In a second step, a linear mixed model including all the different parameters plus BMI and age was calculated. Significant associations were revealed between progression of NAFLD and endoscopically detected inflammation.

DISCUSSION

Our study is one of the first large retrospective studies also including longitudinal analyses to explore changes of NAFLD severity during the course of IBD in up to 18 years of follow-up. As key findings, we confirmed that higher BMI and age are risk factors for the development of NAFLD in IBD patients, while clinical disease activity as expressed by higher HBI scores proved to be another independent risk factor in CD patients.

In patients with steatosis of > 33% of hepatocytes, liver ultrasound is a very reliable tool to detect steatosis. However, the method is not suitable to exclude mild fatty infiltration of the liver[20]. Therefore, we are aware that our results may underestimate the prevalence of NAFLD.

In the cross-sectional part of the study, 157 (52%) of the CD patients had no NAFLD, while mild NAFLD was diagnosed in 54 (18%), moderate NAFLD in 50 (16%), and severe NAFLD in 41 (14%) of the patients. In UC patients, corresponding results were 85 (56%), 30 (20%), 25 (16%), and 13 (8%).

In the longitudinal part of this study, 2467 hepatic ultrasound results were included. Age, dyslipidemia and BMI in CD patients and BMI and hypertension in UC patients were identified as risk factors in the mixed logistic regression model. Also, in CD patients, higher HBI, endoscopic disease activity, history of bowel resection(s), and use of azathioprine were risk factors for NAFLD. In UC patients, endoscopic disease activity was significantly associated with NAFLD in the mixed logistic regression model.

In a recent smaller cohort of Greek IBD patients, NAFLD was revealed in 14.2% of the included CD patients and in 26.2% of the included UC patients[21]. In comparison, our data identify a higher prevalence of NAFLD in both patient groups. The same is the case if we compare our results with those of Bargiggia et al[3] from 2003. Concerning the Greek study, the discrepancy might be explained by diverging dietary habits. Looking at the study by Bargiggia et al[3], it should be noted that nutrition has relevantly changed over the last 15 years. One more important factor to bear in mind is the improvement of IBD therapies and thus the better management of malnutrition in IBD patients over time.

A recently published meta-analysis included one study from Japan and 18 studies from Western countries: It found an overall prevalence of NAFLD in 27.5% of IBD patients and identified type 2 diabetes, hypertension, obesity, insulin resistance, metabolic syndrome, chronic kidney disease, methotrexate use, surgery for IBD, and longer duration of IBD as risk factors for its development[22].

We also identified greater age in CD patients as a risk factor for NAFLD. This finding may be due to a higher prevalence of metabolic risk factors with rising age.

In Germany, the mean BMI in the general population in the year 2013 was 26.5 kg/m² in females and 27.2 kg/m² in males[23]. Median BMI in our study as compared to the mean BMI in Germany was lower: 25.5 kg/m² in CD patients and 25.7 kg/m² in UC patients. That is why higher BMI alone cannot explain the high prevalence of NAFLD in our IBD patients. Furthermore, metabolic syndrome was only diagnosed in 23% of IBD patients in the United States[24].

In contrast to the results of the above-cited meta-analysis, we found in our cohort no association between methotrexate use and NAFLD, but between azathioprine use and an increase of NAFLD severity (OR: 1.98). In our cohort, azathioprine was prescribed far more frequently than methotrexate. All other drug therapies of IBD including corticosteroids were not associated with the course of NAFLD in our patients. Interestingly, in contrast, a study from 2013 identified corticosteroids as potential risk factor for NAFLD[10], while a more recent study showed no difference of steroid use in NAFLD and non-NAFLD IBD patients[25].

The positive correlation of clinical and endoscopic disease activity with NAFLD deterioration in the longitudinal part of the study leads us to the assumption that the systemic, extraintestinal inflammation caused by IBD may promote the development of NAFLD. Also, a history of bowel resection(s) resulted in a progression of NAFLD over time. This observation supports the hypothesis that bowel inflammation may aggravate NAFLD, as surgery is usually performed in patients with aggressive disease courses.

Adipokines also play a role in intestinal inflammation[26]. There may be a link between visceral fatty tissue, NALFD and intestinal inflammation, all factors influencing one another in a detrimental way. For example, apelin is upregulated in the intestinal epithelial cells of IBD patients and highly expressed in the mesenteric adipose tissue of Crohn’s disease patients[27,28]. Another study group showed that apelin is promoting hepatic fibrosis[29].

There are several limitations of the study. Due to its retrospective, uncontrolled nature, risk of bias is given. We focused on patients visiting a university IBD outpatient clinic. 184 CD and 128 UC patients did not undergo any abdominal ultrasound examination due to only getting a second opinion, which may produce a selection bias. Patients who did not receive immunosuppressive therapy did not undergo routine ultrasound examinations as well. Another limitation is the examiner dependency of the ultrasound examination results. Further limitations are that vibration controlled transient elastography or MR spectroscopy were not routinely performed to grade steatosis and fibrosis; nor was liver histology examined. Differentiation between NAFLD and NASH is only possible by histology. Also, nothing is known about the nutrition or the exact metabolic status of the patients. However, elevated plasma γGT concentrations and the tendency of elevated transaminases in CD patients in our study support the correct discrimination of NAFLD by liver ultrasound.

One of the strengths of our study is the multitude of liver ultrasound examinations performed over a long duration. Thus, our study implies the possibility to reveal the influence of parallel therapies for IBD and the course of NAFLD under real-world conditions. Our data is relatively homogeneous, due to restriction to patients treated at a university center for IBD. In the cross-sectional component of our study, the data was retrieved exclusively from May 2017 until 2018, when data documentation at our center had been optimized by new standard operation procedures for the outpatient clinic.

CONCLUSION

We conclude from our results that NAFLD is very common in IBD patients, especially in older patients and patients with higher BMI. Disease activity and azathioprine use in CD may increase the prevalence of NAFLD. As visceral fat seems to perpetuate inflammation, IBD patients with NAFLD may be at risk of more aggressive IBD course, so they may require more intense anti-inflammatory treatment strategies. Our study provides a rationale for regular ultrasound screening of the liver in IBD patients. Further prospective studies are necessary to detect the role of disease activity and the influence of different medication in NAFLD in IBD patients.

ARTICLE HIGHLIGHTS

Research background

Nonalcoholic fatty liver disease (NAFLD) is common in the German population, with an even higher prevalence in inflammatory bowel disease patients.

Research motivation

Finding reasons for the high prevalence of NAFLD in inflammatory bowel disease patients.

Research objectives

Identifying risk factors for NAFLD in inflammatory bowel disease patients.

Research methods

Analysis of abdominal ultrasound examinations of inflammatory bowel disease patients and electronically available data was performed. Demographic characteristics, disease activity and medications were analyzed as potential risk factors in a cross-sectional design. Parameters influencing the course of NAFLD were identified by a generalized linear mixed model.

Research results

Forty-eight percent of Crohn’s disease (CD) patients and 44% of ulcerative colitis patients suffered from NAFLD. Its occurrence was associated with greater age, hypertension and body mass index (BMI) in both groups, and with higher disease activity and dyslipidemia in CD. 2467 ultrasound results were included in the longitudinal analysis. Risk factors for NAFLD were age, BMI, higher disease activity, bowel resection(s), endoscopic activity and azathioprine use in CD; and BMI and endoscopic activity in ulcerative colitis.

Research conclusions

This study shows that NAFLD is highly prevalent in a cohort of German inflammatory bowel disease patients. Its risk increased mainly with rising age and BMI. The analysis provides a rationale for non-invasive liver screening in inflammatory bowel disease patients.

Research perspectives

A randomized controlled prospective study is required to further verify the discovered risk factors.