This article has Open Peer Review reports available.
Non-alcoholic fatty liver disease associated with gallstones in females rather than males: a longitudinal cohort study in Chinese urban population
© Liu et al.; licensee BioMed Central. 2014
Received: 2 November 2013
Accepted: 5 December 2014
Published: 13 December 2014
Whether non-alcoholic fatty liver disease (NAFLD) is a risk factor for gallstones remains uncertain. Few longitudinal or cohort studies have been used to identify this relationship. The aim of this study was to confirm the association between NAFLD and gallstones in a longitudinal cohort of urban dwellers in China.
To elucidate the association between NAFLD and gallstones, we fitted a generalized estimating equation (GEE) model in a large-scale longitudinal cohort over 6 years, which included 11,200 participants with at least three regular health check-ups.
A total of 498 cases of gallstones occurred during the 6-year follow-up, which resulted in a total incidence density of 12.73 per 1000 person-years (498/39, 135.5 person-years). The GEE analyses confirmed and clarified the association between NAFLD and gallstones (relative risk (RR) = 1.2381, 95% confidence interval (CI) = 1.003–1.528, P = 0.047) after adjusting for other potential confounding factors, especially in females (RR = 1.707, 95% CI = 1.245–2.341, P = 0.001).
NAFLD is associated with gallstones in an urban Chinese population from the middle to upper socioeconomic strata. Moreover, this association is more strongly apparent in females than in males. Further cohort studies must be conducted to confirm this association in the general population.
KeywordsGallstones Non-alcoholic fatty liver disease (NAFLD) Longitudinal cohort study Generalized estimated equation (GEE)
Gallstones are common digestive disorders that constitute a significant health problem in Western countries . In recent years, their presence has also been increasing in China. For persons older than 20 years, the prevalence of gallstone disease is rapidly increasing to as high as 3.8% in southern China and 6.1% in northern China . Some factors including ethnicity ,, genetic susceptibility , age, gender, obesity , insulin resistance , and type 2 diabetes mellitus , have been reported to be risk factors for gallstone development. However, the pathogenic mechanisms of gallstone formation have not been fully ascertained. Further explorations of risk factors for gallstone development are needed to aid in the early identification and prediction of gallstones.
Non-alcoholic fatty liver disease (NAFLD) describes a wide extent of liver conditions ranging from simple steatosis to non-alcoholic steatohepatitis and cirrhosis. The prevalence of NAFLD is reported as approximately 15% in Asians . It is not yet clear as to what role NAFLD may play in gallstone formation. Various cross-sectional studies have indicated that NAFLD is associated with gallstones -. The prevalence of gallstone disease was reported to be higher in NAFLD patients than in the general adult population in Italy ,. NAFLD was found to be an independent risk predictor for gallstone disease in Slovakian patients with metabolic risk factors . The association between fatty liver and gallstones was also found in the Chinese  and Japanese  populations. However, conclusions from cross-sectional studies remain insufficient to confirm whether NAFLD is an independent risk factor for gallstones. There are three potential hypotheses to explain why gallstones are always accompanied by NAFLD. (1) Given that both gallstones and NAFLD have a high prevalence in the general population, it is likely to be a chance co-occurrence. (2) Gallstones and NAFLD share common risk factors -, including obesity, insulin resistance, and type 2 diabetes mellitus. (3) NAFLD is an independent risk factor for gallstones.
Sustained efforts are presently directed at clarifying the relationship between NAFLD and gallstones, with the goal of identifying high-risk groups and providing clues for exploring the pathogenic mechanisms of gallstone formation. Results from a large population-based case–control study in China have shown that a history of liver conditions, especially liver cirrhosis, is a risk factor for gallstone development . However, conclusions from such a case–control study are still not sufficient to confirm whether NAFLD is an independent risk factor for gallstones. Therefore, we conducted a large-scale longitudinal cohort study, based on health check-ups, in an urban Han Chinese population from middle to upper socioeconomic strata to further confirm the association between NAFLD and gallstones. A simple generalized estimating equation (GEE) model was first used to select potential confounding factors of gallstone development, and the multiple GEE model was further adopted to detect the association between NAFLD and gallstones. The longitudinal design allowed us to use repeated observations of the same set of variables during the follow-up, and the GEE model was adjusted to explore the inherent correlations among the observations.
Patients’ characteristics and laboratory tests
All participants underwent a general health questionnaire, physical examination, and fasting laboratory assay. The general health questionnaire covered smoking, alcohol intake, sleeping quality, and physical activity. The physical examination included height, weight, systolic blood pressure, and diastolic blood pressure. The body mass index (BMI) was calculated as weight (kg) divided by height (m) squared. After 5 minutes of seated rest, blood pressure was measured from the right arm of the seated participant using an automated sphygmomanometer. Fasting laboratory assays included total glucose, total cholesterol, low-density lipoprotein, high-density lipoprotein, triglycerides, serum total protein, serum albumin, serum globulins, blood urea nitrogen, serum creatinine, hemoglobin, mean corpuscular hemoglobin, red blood cell distribution width, white blood cell count, platelet distribution width, mean platelet volume, and thrombocytocrit. This study was approved by the Ethics Committee of the School of Public Health, Shandong University, and written informed consent for participation in the study was obtained from all participants.
Definition of gallstones
Ultrasonography was used to diagnose gallstones based on the sonographic evidence of gallstones (one or several echogenic, distally shadowing, possibly movable structures in the gallbladder). Cholecystectomy was defined as the absence of a gallbladder on abdominal ultrasonography. As the individuals in our study were healthy, there were few cases of cholecystectomy after applying exclusion criteria such as fewer than three routine health check-ups. Therefore, cholecystectomy was excluded from our study.
Definition of non-alcoholic fatty liver disease
According to the revised definition for NAFLD by the Chinese Hepatology Association dating from February 2006 , NAFLD was diagnosed by abdominal ultrasonography as brightness in the liver and a diffusely echogenic change in the liver parenchyma, with other known etiologic factors of chronic liver disease (alcohol intake >20 g/day, hepatitis B antigen or hepatitis C virus antibody positive, and autoimmune hepatitis) excluded.
Missing data imputation
Multiple imputations were used to account for missing values. Since our imputation method was selected depending on the pattern of missing data and the type of the imputed variables, the Markov chain Monte Carlo method was chosen according to the MI Procedure of SAS 9.2 . Data were imputed only for covariates. Before imputation, the extent of missing data in all variables was less than 20%; in particular, less than 17% for drinking, smoking, quality of sleep, physical activity, triglycerides, high-density lipoprotein, and low-density lipoprotein, and 10% for other variables. After imputation, all variables had less than 10% missing observations; in particular, less than 6% for drinking, smoking, quality of sleep, physical activity, triglycerides, high-density lipoprotein, and low-density lipoprotein, and 1% for other variables.
Statistical descriptions are expressed as mean ± standard deviation for continuous variables and as frequencies for categorical variables. Summary statistics were obtained for variables of interest at both baseline and follow-up. Mean values were compared by t test for unequal variances. Frequencies were compared by the chi-squared test.
As the age and gender had been confirmed as the chief risk factors for gallstones, the association between gallstones and other factors were examined after adjusting for baseline age; moreover, males and females were analyzed both together and separately. Simple GEE models were first used to select factors associated with gallstones, and multiple GEE models were further adopted to detect the association between NAFLD and gallstones . The GEE models used ‘Logit’ as the link function. Furthermore, the GEE model was fitted to test the regression coefficient of interaction terms between gender and NAFLD status. All statistical analyses were performed using SAS 9.2 (SAS Institute, Cary, NC, USA). A two-sided P value of less than 0.05 was considered to be statistically significant.
Additional file 1: Figure S1 shows the structure of the cohort during the 6 years of follow-up. Of the 11,200 subjects fulfilling enrollment criteria (5227 males and 5973 females), a total of 498 gallstone events occurred over the 6-year follow-up, leading to a total incidence density of 12.73 per 1000 person-years (498/39, 135.5 person-years), with 7.1, 14.5, 16.6, 18.7, and 2.6 per 1000 person-years of follow-up in the 1st, 2nd, 3rd, 4th, and 5th year, respectively.
Distribution of potential confounding factors of participants grouped by NAFLD status at baseline
43.14 ± 14.73
51.11 ± 13.89
40.65 ± 12.61
54.11 ± 13.29
23.47 ± 3.08
27.38 ± 2.82
22.83 ± 3.06
27.35 ± 3.23
120.81 ± 19.06
135.42 ± 19.49
115.94 ± 17.67
136.09 ± 21.31
27.44 ± 4.00
28.19 ± 4.30
27.81 ± 3.93
29.54 ± 4.41
46.35 ± 2.46
46.59 ± 2.40
46.02 ± 2.40
45.78 ± 2.56
4.86 ± 1.26
5.32 ± 1.22
4.51 ± 1.16
4.92 ± 1.22
77.56 ± 14.39
83.81 ± 14.09
70.28 ± 10.12
70.98 ± 10.59
5.01 ± 0.93
5.56 ± 1.30
4.90 ± 0.81
5.73 ± 1.47
4.89 ± 0.94
5.32 ± 0.95
4.88 ± 0.97
5.55 ± 1.04
1.13 ± 0.87
2.06 ± 1.51
0.99 ± 0.79
1.91 ± 1.39
1.38 ± 0.33
1.23 ± 0.31
1.47 ± 0.32
1.35 ± 0.32
2.76 ± 0.73
3.12 ± 0.71
2.69 ± 0.74
3.22 ± 0.77
142.43 ± 14.85
151.64 ± 13.22
134.20 ± 11.51
138.62 ± 10.96
29.85 ± 1.93
29.98 ± 1.70
29.51 ± 2.10
29.41 ± 2.06
41.41 ± 2.58
41.49 ± 2.46
41.35 ± 2.52
41.28 ± 2.30
6.30 ± 1.50
7.08 ± 1.57
6.14 ± 1.47
6.95 ± 1.62
12.30 ± 1.72
12.30 ± 1.70
12.31 ± 1.72
12.16 ± 1.57
10.41 ± 0.82
10.35 ± 0.80
10.46 ± 0.81
10.32 ± 0.76
0.25 ± 0.08
0.25 ± 0.07
0.26 ± 0.09
0.27 ± 0.06
Gallstone events by non-alcoholic fatty liver disease (NAFLD) status
Gallstone events (%)
Event rate (per 1000 person-years)
Additional file 3: Tables S2, Additional file 4: Table S3, Additional file 5: Table S4 show the results of age-adjusted simple GEE models in all participants, and in males and females, respectively. In all participants, BMI, systolic blood pressure, serum albumin, serum globulins, glucose, and triglycerides were statistically significant (Additional file 3: Table S2). In gender-specific analyses, BMI, systolic blood pressure, serum albumin, serum globulins, glucose, triglycerides, and white blood cell count were statistically significant for males (Additional file 4: Table S3). Serum albumin, serum globulins, and glucose were statistically significant for females (Additional file 5: Table S4).
Results of generalized estimating equation (GEE) analysis for non-alcoholic fatty liver disease (NAFLD) and gallstones with their risk ratio (RR) and 95% confidence intervals (CI)
In this longitudinal cohort study, our findings have mainly shown that NAFLD is an independent risk factor for gallstone development, especially in females. As the pathogenic mechanisms of gallstone formation have not been fully ascertained, accurate assessment of a high-risk population becomes crucial to decreasing the prevalence of gallstones. Approximately 10% of females with NAFLD developed gallstones during the 6-year follow-up, making this a high-risk population for gallstones.
The relevant findings of this study are as follows. (1) An elevated incidence density of gallstones in the NAFLD group was observed in comparison with the non-NAFLD group. This finding is in accordance with the results of cross-sectional studies in various populations from Italy, Japan, and China -,. (2) The time sequence of advancement from NAFLD to the presence of gallstones has been confirmed, providing solid evidence that NAFLD is an independent risk factor for gallstones in females.
Specifically, females with NAFLD are about two times as likely as females without NAFLD to have gallstones, after adjustments for potential confounding factors (see Table 3), but such a trend is not observed in males. Various studies have shown that gallstones are a female-predominant disease , whereas NAFLD is predominant in males . However, whether males or females are more susceptible to gallstones in a NAFLD population remains unknown. Although a cross-sectional study showed that the association between NAFLD and gallstones was stronger in males than in females , the time sequence between NAFLD and gallstones was not ascertained, owing to the disadvantageous study design. Moreover, this discrepancy might also be due to differences in ethnicity. Our present longitudinal cohort study confirms that NAFLD is associated with gallstones in females rather than males (see Table 3), implying a continuing such trend in the future. Therefore, more attention should be paid to the female population with NAFLD in compiling clinical guidelines. Based on various studies ,-, the authors propose two hypotheses to explain the gender difference regarding the association between NAFLD and gallstones. (1) Males tend to develop an antioxidant reaction through overload of uric acid while females react more easily through glucose metabolism, which is associated with biliary cholesterol secretion and subsequent promotion of cholesterol gallstone formation. (2) Estrogen aggravates lipoapoptosis in females with NAFLD, which results in supersaturated cholesterol of bile, further promoting the formation of crystalline precipitate and ultimately leading to gallstone formation. However, these hypotheses need further research on their molecular mechanisms before they can be verified.
The advantage of the present study is that it is based on a large-scale longitudinal cohort design, which ascertains the time sequence from NAFLD to gallstones to confirm that NAFLD is an independent risk factor for gallstone development. Although the causal relationship between NAFLD and gallstones could not be examined by this longitudinal study, it has nonetheless provided solid evidence that NAFLD is an independent risk factor for gallstones in female subjects. However, some disadvantages do exist in our study. Owing to the inherent weakness of our routine health check-up database, we were unable to obtain information on family history and the use of medication. Considering that the study cohort contained only the urban population of middle to upper socioeconomic strata in Shandong province, and that no information on alanine aminotransferase and aspartate aminotransferase were included, further larger validation cohort studies are needed to confirm our results.
NAFLD is associated with gallstones in an urban Chinese population of middle to upper socioeconomic strata. In gender-specific analyses, this association was strongly confirmed in females but not in males. Further cohort studies are required to confirm this association in the general population.
This work was supported by grants from the National Nature Science Fund (No. 81273177 and No. 81273082) of China. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. We thank Hugh McGonigle who provided language editing services on behalf of Edanz.
- Everhart JE, Khare M, Hill M, Maurer KR: Prevalence and ethnic differences in gallbladder disease in the United States. Gastroenterology. 1999, 117: 632-639. 10.1016/S0016-5085(99)70456-7.View ArticlePubMedGoogle Scholar
- Zeng Q, He Y, Qiang DC, Wu LX: Prevalence and epidemiological pattern of gallstones in urban residents in China. Eur J Gastroenterol Hepatol. 2012, 24: 1459-1460. 10.1097/MEG.0b013e3283583d13.View ArticlePubMedGoogle Scholar
- Shaffer EA: Gallstone disease: epidemiology of gallbladder stone disease. Best Pract Res Clin Gastroenterol. 2006, 20: 981-996. 10.1016/j.bpg.2006.05.004.View ArticlePubMedGoogle Scholar
- Kosters A, Jirsa M, Groen AK: Genetic background of cholesterol gallstone disease. Biochim Biophys Acta. 2003, 1637: 1-19. 10.1016/S0925-4439(02)00173-4.View ArticlePubMedGoogle Scholar
- Everhart JE: Contributions of obesity and weight loss to gallstone disease. Ann Intern Med. 1993, 119: 1029-1035. 10.7326/0003-4819-119-10-199311150-00010.View ArticlePubMedGoogle Scholar
- Tsai CJ, Leitzmann MF, Willett WC, Giovannucci EL: Macronutrients and insulin resistance in cholesterol gallstone disease. Am J Gastroenterol. 2008, 103: 2932-2939. 10.1111/j.1572-0241.2008.02189.x.View ArticlePubMedGoogle Scholar
- Pagliarulo M, Fornari F, Fraquelli M, Zoli M, Giangregorio F, Grigolon A, Peracchi M, Conte D: Gallstone disease and related risk factors in a large cohort of diabetic patients. Dig Liver Dis. 2004, 36: 130-134. 10.1016/j.dld.2003.10.007.View ArticlePubMedGoogle Scholar
- De Alwis NM, Day CP: Non-alcoholic fatty liver disease: the mist gradually clears. J Hepatol. 2008, 48 (Suppl 1): S104-S112. 10.1016/j.jhep.2008.01.009.View ArticlePubMedGoogle Scholar
- Fracanzani AL, Valenti L, Russello M, Miele L, Bertelli C, Bellia A, Masetti C, Cefalo C, Grieco A, Marchesini G, Fargion S: Gallstone disease is associated with more severe liver damage in patients with non-alcoholic fatty liver disease. Plos One. 2012, 7: e41183-10.1371/journal.pone.0041183.View ArticlePubMedPubMed CentralGoogle Scholar
- Loria P, Lonardo A, Lombardini S, Carulli L, Verrone A, Ganazzi D, Rudilosso A, D'Amico R, Bertolotti M, Carulli N: Gallstone disease in non-alcoholic fatty liver: prevalence and associated factors. J Gastroenterol Hepatol. 2005, 20: 1176-1184. 10.1111/j.1440-1746.2005.03924.x.View ArticlePubMedGoogle Scholar
- Koller T, Kollerova J, Hlavaty T, Huorka M, Payer J: Cholelithiasis and markers of nonalcoholic fatty liver disease in patients with metabolic risk factors. Scand J Gastroenterol. 2012, 47: 197-203. 10.3109/00365521.2011.643481.View ArticlePubMedGoogle Scholar
- Chen CH, Huang MH, Yang JC, Nien CK, Etheredge GD, Yang CC, Yeh YH, Wu HS, Chou DA, Yueh SK: Prevalence and risk factors of gallstone disease in an adult population of Taiwan: an epidemiological survey. J Gastroenterol Hepatol. 2006, 21: 1737-1743. 10.1111/j.1440-1746.2006.04381.x.View ArticlePubMedGoogle Scholar
- Nomura H, Kashiwagi S, Hayashi J, Kajiyama W, Ikematsu H, Noguchi A, Tani S, Goto M: Prevalence of gallstone disease in a general population of Okinawa, Japan. Am J Epidemiol. 1988, 128: 598-605.PubMedGoogle Scholar
- Hsing AW, Gao YT, McGlynn KA, Niwa S, Zhang M, Han TQ, Wang BS, Chen J, Sakoda LC, Shen MC, Zhang BH, Deng J, Rashid A: Biliary tract cancer and stones in relation to chronic liver conditions: a population-based study in Shanghai, China. Int J Cancer. 2007, 120: 1981-1985. 10.1002/ijc.22375.View ArticlePubMedGoogle Scholar
- Zeng MD, Fan JG, Lu LG, Li YM, Chen CW, Wang BY, Mao YM: Guidelines for the diagnosis and treatment of nonalcoholic fatty liver diseases. J Dig Dis. 2008, 9: 108-112. 10.1111/j.1751-2980.2008.00331.x.View ArticlePubMedGoogle Scholar
- Schafer JL: Multiple imputation: a primer. Stat Methods Med Res. 1999, 8: 3-15. 10.1191/096228099671525676.View ArticlePubMedGoogle Scholar
- Hubbard AE, Ahern J, Fleischer NL, Van der Laan M, Lippman SA, Jewell N, Bruckner T, Satariano WA: To GEE or not to GEE: comparing population average and mixed models for estimating the associations between neighborhood risk factors and health. Epidemiology. 2010, 21: 467-474. 10.1097/EDE.0b013e3181caeb90.View ArticlePubMedGoogle Scholar
- Ruhl CE, Everhart JE: Relationship of non-alcoholic fatty liver disease with cholecystectomy in the US Population. Am J Gastroenterol. 2013, 108: 952-958. 10.1038/ajg.2013.70.View ArticlePubMedGoogle Scholar
- Stinton LM, Shaffer EA: Epidemiology of gallbladder disease: cholelithiasis and cancer. Gut Liver. 2012, 6: 172-187. 10.5009/gnl.2012.6.2.172.View ArticlePubMedPubMed CentralGoogle Scholar
- Hashimoto E, Tokushige K: Prevalence, gender, ethnic variations, and prognosis of NASH. J Gastroenterol. 2011, 46 (Suppl 1): 63-69. 10.1007/s00535-010-0311-8.View ArticlePubMedGoogle Scholar
- Lonardo A, Trande P: Are there any sex differences in fatty liver? A study of glucose metabolism and body fat distribution. J Gastroenterol Hepatol. 2000, 15: 775-782. 10.1046/j.1440-1746.2000.02226.x.View ArticlePubMedGoogle Scholar
- Xie Y, Wang M, Zhang Y, Zhang S, Tan A, Gao Y, Liang Z, Shi D, Huang Z, Zhang H, Yang X, Lu Z, Wu C, Liao M, Sun Y, Qin X, Hu Y, Li L, Peng T, Li Z, Mo Z: Serum uric acid and non-alcoholic fatty liver disease in non-diabetic Chinese men. Plos One. 2013, 8: e67152-10.1371/journal.pone.0067152.View ArticlePubMedPubMed CentralGoogle Scholar
- Targher G, Byrne CD: Clinical review: nonalcoholic fatty liver disease: a novel cardiometabolic risk factor for type 2 diabetes and its complications. J Clin Endocrinol Metab. 2013, 98: 483-495. 10.1210/jc.2012-3093.View ArticlePubMedGoogle Scholar
- Cirillo DJ, Wallace RB, Rodabough RJ, Greenland P, LaCroix AZ, Limacher MC, Larson JC: Effect of estrogen therapy on gallbladder disease. JAMA. 2005, 293: 330-339. 10.1001/jama.293.3.330.View ArticlePubMedGoogle Scholar
- Henriksson P, Einarsson K, Eriksson A, Kelter U, Angelin B: Estrogen-induced gallstone formation in males. Relation to changes in serum and biliary lipids during hormonal treatment of prostatic carcinoma. J Clin Invest. 1989, 84: 811-816. 10.1172/JCI114240.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang HH, Liu M, Clegg DJ, Portincasa P, Wang DQ: New insights into the molecular mechanisms underlying effects of estrogen on cholesterol gallstone formation. Biochim Biophys Acta. 2009, 1791: 1037-1047. 10.1016/j.bbalip.2009.06.006.View ArticlePubMedPubMed CentralGoogle Scholar
- Puri P, Baillie RA, Wiest MM, Mirshahi F, Choudhury J, Cheung O, Sargeant C, Contos MJ, Sanyal AJ: A lipidomic analysis of nonalcoholic fatty liver disease. Hepatology. 2007, 46: 1081-1090. 10.1002/hep.21763.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.