- Research article
- Open Access
Non-alcoholic fatty liver disease associated with gallstones in females rather than males: a longitudinal cohort study in Chinese urban population
BMC Gastroenterology volume 14, Article number: 213 (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.
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.
A large-scale longitudinal cohort study was set up in 2005, comprising urban Han Chinese from middle to upper socioeconomic strata who attended a routine health check-up at Centers for Health Management of Shandong Provincial Qianfoshan Hospital and Shandong Provincial Hospital. The geographic area served by the two hospitals covered the city of Jinan, the capital of Shandong province with a population of approximately 7 million. The health check-up database contained relatively rich and educated persons in Jinan, representing the middle to upper class population of Shandong province. Based on the large-scale longitudinal cohort, a sub-cohort was selected from those free of gallstones at baseline. A total of 11,200 participants with at least three health checks in the 6-year follow-up were included in the sub-cohort between 2005 and 2010 (Figure 1).
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.
The distributions of potential confounding factors of participants grouped by NAFLD status at baseline are shown in Table 1. The baseline prevalence of NAFLD was 30.55% in males and 13.11% in females. At baseline, statistically significant differences between non-NAFLD and NAFLD groups for the listed variables were observed except for red blood cell distribution width, platelet distribution width, and thrombocytocrit in males, and serum creatinine, mean corpuscular hemoglobin, and red blood cell distribution width in females. Additional file 2: Table S1 shows the distributions of potential confounding factors at the baseline survey, and their distributions between subjects with gallstones and subjects without gallstones at each follow-up interval.
Table 2 shows the frequency and corresponding event rate per 1000 person-years of gallstone events among subjects with and without NAFLD. Of the 4713 individuals who suffered NAFLD, 289 (6.1%) developed gallstones during the follow-up years, leading to an incidence density of 31.75 per 1000 person-years. For the individuals without NAFLD, 209 of 6487 (3.2%) developed gallstones, with the event rate decreasing to 6.96 per 1000 person-years. When subgroup analysis was performed in males and females, the frequency of gallstones in the NAFLD group was over 20 times higher than in the non-NAFLD group in females, whereas in males the frequency of gallstones in both groups was similar.
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).
Table 3 shows the summarized results of GEE analyses regarding the association between NAFLD and gallstones (for detailed information, see Additional file 3: Tables S2, Additional file 4: Table S3, Additional file 5: Table S4, Additional file 6: Table S5, Additional file 7: Table S6 and Additional file 8: Table S7). The following results were obtained. (1) NAFLD was associated with gallstones (relative risk (RR) = 1.427, P < 0.001); after adjusting for age in all participants (Additional file 3: Table S2), this association was still significant (RR = 1.238, 95% confidence interval (CI) =1.003–1.528, P = 0.047), after further adjusting for other potential confounders in the multiple GEE model (Additional file 6: Table S5). (2) When using gender-specific analyses the association became stronger in females, not only in the age-adjusted simple GEE model (RR =1.807, 95% CI = 1.277–2.558, P = 0.001) but also in the multiple GEE model (RR = 1.707, 95% CI = 1.245–2.341, P = 0.001) after further adjusting for other potential confounders (Additional file 4: Tables S3 and Additional file 7: Table S6), whereas no association was detected in males (Additional file 5: Tables S4 and Additional file 8: Table S7). Moreover, the interaction between gender and NAFLD status was statistically significant (P = 0.016). These results indicate that NAFLD was an independent risk factor for gallstones in the female rather than the male population.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Schafer JL: Multiple imputation: a primer. Stat Methods Med Res. 1999, 8: 3-15. 10.1191/096228099671525676.
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.
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.
Stinton LM, Shaffer EA: Epidemiology of gallbladder disease: cholelithiasis and cancer. Gut Liver. 2012, 6: 172-187. 10.5009/gnl.2012.6.2.172.
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.
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.
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.
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.
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.
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.
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.
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.
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.
The authors declare that they have no competing interests.
In our work, YL, FX and CZ designed the study and directed its implementation, including quality assurance and control. FT, HL and DZ did the clinical exam and collected the data. LW, WS helped analyzing the data. JL participated much of the above work and led the writing of the paper. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 3: Table S2.: Single-predictor generalized estimating equation (GEE) models in all participants. (DOC 50 KB)
Additional file 4: Table S3.: Single-predictor generalized estimating equation (GEE) models in female. (DOC 49 KB)
Additional file 5: Table S4.: Single-predictor generalized estimating equation (GEE) models in male. (DOC 48 KB)
Additional file 6: Table S5.: Results of multiple generalized estimating equation (GEE) analysis for non-alcoholic fatty liver disease(NAFLD) and gallstones in all participants after adjusting other potential confounding factors. (DOC 36 KB)
Additional file 7: Table S6.: Results of multiple generalized estimating equation (GEE) analysis for non-alcoholic fatty liver disease (NAFLD) and gallstones in female after adjusting other potential confounding factors. (DOC 34 KB)
Additional file 8: Table S7.: Results of multiple generalized estimating equation (GEE) analysis for non-alcoholic fatty liver disease (NAFLD) and gallstones in male after adjusting other potential confounding factors. (DOC 38 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
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.
About this article
Cite this article
Liu, J., Lin, H., Zhang, C. et al. Non-alcoholic fatty liver disease associated with gallstones in females rather than males: a longitudinal cohort study in Chinese urban population. BMC Gastroenterol 14, 213 (2014). https://doi.org/10.1186/s12876-014-0213-y
- Non-alcoholic fatty liver disease (NAFLD)
- Longitudinal cohort study
- Generalized estimated equation (GEE)