- Research article
- Open Access
Associations of coffee consumption with markers of liver injury in the insulin resistance atherosclerosis study
BMC Gastroenterology volume 15, Article number: 88 (2015)
Coffee consumption has been associated with reduced risk of developing type 2 diabetes mellitus (T2DM) however, the mechanism for this association has yet to be elucidated. Non-alcoholic fatty liver disease (NAFLD) characterizes and predicts T2DM yet the relationship of coffee with this disorder remains unclear. Our aim was to investigate the associations of coffee with markers of liver injury in 1005 multi-ethnic, non-diabetic adults in the Insulin Resistance Atherosclerosis Study.
Dietary intake was assessed using a validated 114-item food frequency questionnaire. Alanine aminotransferase (ALT), aspartate aminotransferase (AST) and fetuin-A were determined in fasting blood samples and the validated NAFLD liver fat score was calculated. Multivariate linear regression assessed the contribution of coffee to variation in markers of liver injury.
Caffeinated coffee showed significant inverse associations with ALT (β = −0.08, p = 0.0111), AST (β = −0.05, p = 0.0155) and NAFLD liver fat score (β = −0.05, p = 0.0293) but not with fetuin-A (β = 0.04, p = 0.17). When the highest alcohol consumers were excluded, these associations remained (ALT β = −0.11, p = 0.0037; AST β = −0.05, p = 0.0330; NAFLD liver fat score β = −0.06, p = 0.0298). With additional adjustment for insulin sensitivity, the relationship with ALT remained significant (ALT β = −0.08, p = 0.0400; AST β = −0.03, p = 0.20; NAFLD liver fat score β = −0.03, p = 0.27). There were no significant associations of decaffeinated coffee with liver markers.
These analyses indicate a beneficial impact of caffeinated coffee on liver morphology and/or function, and suggest that this relationship may mediate the well-established inverse association of coffee with risk of T2DM.
Coffee consumption has consistently been associated with a reduction in risk of developing type 2 diabetes mellitus (T2DM), however few studies have attempted to identify the pathophysiological pathway underlying this relationship [1, 2]. Exploring the impact of coffee consumption on major risk factors for T2DM will be important in uncovering potential mediating relationships underlying the observed beneficial association of coffee with this chronic disease. Non-alcoholic fatty liver disease (NAFLD) has been identified as a predictor of diabetes risk according to a meta-analysis of prospective population-based studies using transaminase concentrations . In addition, three studies reported a significant association of ultrasound-identified NAFLD with incident T2DM [4–6]. However, the relationship of coffee with NAFLD remains unclear . Evidence suggests a beneficial impact of coffee on the liver in individuals with alcoholic and viral liver diseases [8–11] and hepatocellular carcinoma [8, 9, 12], although the relationship of coffee with liver markers is largely unstudied in the context of liver disease of metabolic origin. Existing epidemiological studies support a beneficial impact of coffee consumption on the liver in apparently healthy populations, however the majority of these data were from studies of male populations in Asia [13–18]. A recent study of participants in the National Health and Nutrition Examination Survey (NHANES) extended these investigations to the United States and reported an inverse association of coffee consumption with liver enzymes, including alanine aminotransferase (ALT), aspartate aminotransferase (AST), gamma-glutamyl transferase (GGT) and alkaline phosphatase (ALP) . Examining these relationships in a population with more extensive metabolic characterization, including insulin resistance and pre-diabetes, would be informative.
The existing literature on coffee and liver injury markers is limited in several aspects. First, little is known about the individual effects of caffeinated vs decaffeinated coffee. Xiao et al. (2014) were the first to differentiate between caffeinated and decaffeinated coffee in their examination of the independent relationships of recent intake, assessed by 24 h recall, of these coffee types with markers of liver injury . Additional studies that distinguish between caffeinated and decaffeinated coffee consumption using measures that capture usual, rather than recent, intake of these coffee types are needed in order to assess the independent relationships of habitual intake of coffee with markers of liver injury. Second, previous investigations focused primarily on associations of coffee with routinely measured liver enzymes; examination of coffee’s relationship with other biomarkers of liver function and biomarker-based scores predictive of NAFLD would be instructive. Fetuin-A, a biomarker for inflammation and liver function, has been positively associated with insulin resistance and an increased risk of T2DM [20–22]. Investigations of the relationship of coffee consumption with fetuin-A have been limited to two studies which suggest a beneficial impact of coffee, but these findings need to be corroborated [20, 23]. Biomarker-based predictive scores derived from simple, non-invasive, clinical data are emerging as a practical, yet robust, method for predicting NAFLD [24, 25]. The validated NAFLD liver fat score was recently identified as the best non-invasive prediction score for identifying NAFLD when compared to similar tools . To our knowledge, associations of coffee with this predictive score have yet to be examined. In light of these important knowledge gaps, the objective of the current study was to investigate the associations of usual caffeinated and decaffeinated coffee consumption with markers of NAFLD, including ALT, AST, fetuin-A, and the validated NAFLD liver fat score , in a multi-ethnic cohort.
Study population and design
The Insulin Resistance Atherosclerosis Study (IRAS) is an epidemiological investigation of the relationship of insulin resistance (IR) with cardiovascular risk factors in a large multi-ethnic population . The design of IRAS has been published in detail . Briefly, data collection was initiated in October 1992 with participants recruited from four geographic areas in the United States (San Antonio TX, San Luis Valley CO, Oakland and Los Angeles CA). The recruitment target was to achieve equal participation across categories of sex, ethnicity (non-Hispanic white, Hispanic, African American), age and glucose tolerance status (normal, impaired glucose tolerance, diabetes). In total, the baseline IRAS population included 1625 participants, from which written informed consent was obtained. This study was conducted in accordance with the Helsinki Declaration and the institutional review boards at Kaiser Permanente Division of Research (Oakland, CA), University of California (Los Angeles, CA), University of Texas Health Science Center at San Antonio (San Antonio, TX), the University of Colorado (Denver, CO), and the Bowman Gray School of Medicine (Winston-Salem, NC) approved the project.
For the current cross-sectional study, the population was restricted to individuals without T2DM, as defined by the 1999 World Health Organization criteria . Further exclusions were made on the basis of missing liver enzyme data or extremes of energy intake, resulting in a final study population of n = 1005. Fetuin-A data were not available for all individuals however participant characteristics did not differ substantially between those for which fetuin-A data was available (n = 650) and those for which it was not (n = 355) (data not shown).
Participant assessments were conducted across two 4 h visits approximately 1 week apart and participants were asked to fast for 12 h prior to each visit . Diabetes status was assessed on the first visit using a 75 g oral glucose tolerance test (OGTT). Blood was collected for fasting and 2-h glucose samples. IR was assessed at the second visit using the insulin-modified frequently sampled intravenous glucose tolerance test (FSIGTT), and insulin sensitivity (SI), was calculated using mathematical modeling methods (MINMOD version 3.0, 1994, Los Angeles, CA), as has been described in previous publications [26, 28–30]. Socio-demographic information was collected by participant report and questionnaires assessed energy expenditure, alcohol consumption and smoking behaviours .
Dietary intake over the past year was assessed using a validated interviewer-administered 114-item food frequency questionnaire (FFQ) modified for the IRAS from the National Cancer Institute’s Health Habits and History Questionnaire to include a number of regional and ethnic foods [26, 32]. This measure was validated in a sub-sample of the population through the administration of eight 24-h recalls administered on randomly selected days over a 1 year period. Results demonstrated desirable validity and confirmed the appropriateness of this instrument for use in this population . Although the validation did not specifically address coffee, it is established that coffee drinking behaviour is accurately captured in FFQs (American women r = 0.8, men r = 0.9, respectively, vs dietary records) [33, 34].
Caffeinated and decaffeinated coffee consumption were assessed individually and participants reported on both the frequency and quantity of consumption . Frequency was reported using a nine category scale with responses ranging from “never or less than 1 per month” to “six or more times per day” . Individuals indicated whether their portion sizes were “small, medium or large compared with other men or women about your age” . Frequency and portion size data were integrated in order to provide a single value reflective of intake by weighting the intake frequency by a factor of 0.5, 1 or 1.5 for the reported portion size of small, medium or large, respectively. Thus, one serving of an item reflects one participant-identified medium serving. Intakes of caffeinated and decaffeinated coffee were summed for a measure of total coffee. Categories of vegetable, fruit and whole grain intakes were calculated according to groupings established in previous work in IRAS . Nutrient analysis of FFQ data was conducted using the HHHQ-DIETSYS analysis software (version 3.0, 1993; National Cancer Institute, Bethesda, MD, USA).
Fasting blood samples were used to determine markers of liver injury. Liver enzymes ALT and AST were measured using standard methods at the central IRAS laboratory with a Paramax PLA instrument (Baxter) . Fetuin-A was determined using a sandwich-format immunoassay developed by Tethys (Tethys Bioscience, Emeryville, CA) with antibodies from R&D systems (Minneapolis, USA). The assay had a detection limit of 0.5 ng/ml and an inter-assay coefficient of variation of 11 %. NAFLD liver fat score was calculated according to the equation developed and validated by Kotronen et al. which includes fasting serum insulin, ALT, AST and presence of the metabolic syndrome (International Diabetes Federation harmonized definition) [24, 38].
Statistical analyses were conducted using SAS version 9.2 (SAS Institute, Cary, NC). Participant characteristics were tested for differences across categories of total coffee intake. Chi-square tests were used to test differences in frequencies of categorical variables. For continuous variables, differences in means were tested using ANOVA unless the variable was highly skewed. In this case, the Kruskal-Wallis test for non-parametric variables was used.
Linear regression analyses were conducted to assess the independent relationships of caffeinated and decaffeinated coffee consumption, modeled as continuous variables, with markers of liver injury. The distributions of outcome variables were evaluated and natural log transformations were used for all outcome variables as they resulted in more normal distributions. Due to the presence of negative values, a constant of 6.4 was added to all NAFLD liver fat score data before the transformation. The independent associations of caffeinated and decaffeinated coffee consumption with outcome variables were analyzed using staged multivariate models. A base model adjusted for age, sex, and ethnicity. A second model additionally adjusted for energy intake and expenditure, BMI, smoking, alcohol consumption and education. Finally, a third model additionally adjusted for a number of dietary factors including whole grain consumption, vegetable and fruit intake, percent of energy intake from saturated and polyunsaturated fat as well as intake of the alternate type of coffee consumed (caffeinated or decaffeinated), and sugar-sweetened beverages (regular soda and lemonade/sweetened mineral water). To test for potential modifiers of the relationship between coffee consumption and markers of liver injury, interaction terms were tested for age, sex, ethnicity, BMI and SI, and were judged to be statistically significant if the p value was <0.05.
A sensitivity analysis was performed in which multivariate linear regression analyses were conducted with total coffee consumption as the exposure variable. Because of the well-known association between high alcohol consumption and elevations in liver transaminases, sensitivity analyses were performed in which regression models were repeated when individuals in the two highest categories of alcohol consumption were excluded (≥1 alcoholic drink/day). In addition, a fourth mechanistic regression model which additionally adjusted for insulin sensitivity was added in order to assess the contribution of SI to the observed association of coffee and the outcome variables.
Participant characteristics by categories of total coffee consumption are presented in Table 1. In this population, 79 % of individuals were consuming at least one type of coffee. Amounts of coffee consumed however were relatively modest with a median total coffee intake of 0.50 servings/day [interquartile range (IQR), 0.04-0.79]. Approximately 66 % of the population consumed caffeinated coffee while decaffeinated coffee intake was less common with only 28 % of individuals consuming this type. Ethnicity, level of education and smoking and alcohol consumption behaviours differed significantly across categories of total coffee consumption (Table 1). The highest category of total coffee consumption had the lowest proportion of participants with impaired glucose tolerance, however SI, BMI, waist circumference, family history of diabetes and energy expended did not differ significantly across categories of total coffee consumption. Energy intake, percent of energy from saturated fat and regular soda consumption differed significantly across categories where the highest energy intakes and saturated fat consumption were found in the highest categories of coffee consumption. Caffeinated and decaffeinated coffee consumption were inversely correlated (Spearman r = −0.20, p <0.0001). ALT and fetuin-A were inversely correlated with SI (ALT Spearman r = −0.27, p <0.0001; fetuin-A Spearman r = −0.16, p <0.0001), and positively correlated with each other (Spearman r = 0.11, p = 0.0057).
Caffeinated coffee consumption was inversely associated with ALT, AST and NAFLD liver fat score, and, in the case of ALT, there was a stronger magnitude of association across staged multivariate models (Table 2, model 3: ALT β = −0.08, p = 0.0111; AST β = −0.05, p = 0.0155; NAFLD liver fat score β = −0.05, p = 0.0293). Caffeinated coffee consumption however, showed no significant association with fetuin-A (Table 2, model 3: β = 0.04, p = 0.17). There were no significant associations of decaffeinated coffee consumption with markers of liver injury.
Interaction and sensitivity analyses
Associations of coffee with markers of liver injury were not modified by age, sex, ethnicity, BMI or SI (all p > 0.05).
Multivariate linear regression analyses were also conducted with total coffee as the exposure variable. For the liver enzymes, results of regression analyses with total coffee consumption were consistent with those of caffeinated coffee consumption, with significant associations of total coffee intake with both ALT and AST (Additional file 1: Table S1, model 3: ALT β = −0.07, p = 0.0177; AST β = −0.05, p = 0.0131). Total coffee consumption was not however, associated with NAFLD liver fat score (Additional file 1: Table S1, model 3: β = −0.04, p = 0.11) or fetuin-A, (Additional file 1: Table S1, model 3: β = 0.02, p = 0.47).
When the heaviest alcohol consumers were excluded, the significant inverse association of caffeinated coffee with ALT, AST and NAFLD liver fat score not only remained, but, in the case of ALT, there was a stronger magnitude of association compared to the initial analysis (Table 3, model 3: ALT β = −0.11, p = 0.0037; AST β = −0.05, p = 0.0330; NAFLD liver fat score β = −0.06; p = 0.0298). With additional adjustment for SI in a fourth (mechanistic) model, the significant inverse association of caffeinated coffee with ALT remained, however the association with AST and NAFLD liver fat score was attenuated to non-significance (Table 3, model 4: ALT β = −0.08, p = 0.0400; AST β = −0.03, p = 0.20; NAFLD liver fat score β = −0.03, p = 0.27). Associations of caffeinated coffee with fetuin-A (Table 3, model 4: β = 0.03, p = 0.44), and decaffeinated coffee with all liver markers remained unchanged from the original analysis (data not shown).
NAFLD has emerged as an important predictor of T2DM [37, 39] and thus it is of interest to investigate the relationship of this disorder with coffee consumption in the context of attempting to identify a potential mediating role on the inverse relationship of coffee with T2DM. Our findings regarding ALT, AST and NAFLD liver fat score shed light on a potential mediating relationship that may underlie this inverse association [1, 2].
Results of the present cross-sectional examination of a multi-ethnic, non-diabetic cohort reveal significant inverse associations of caffeinated coffee consumption with ALT and AST, and in the case of ALT, there was a stronger magnitude of association across multivariate models and with the exclusion of individuals whose transaminase levels may be affected by alcohol consumption. Analyses also revealed an inverse relationship of caffeinated coffee with NAFLD liver fat score, a non-invasive validated tool for identifying NAFLD using routinely available clinical and laboratory data . These findings are consistent with accumulating evidence demonstrating a beneficial impact of coffee consumption on the liver. Specifically, a beneficial effect of coffee on the risk of NAFLD as defined by elevated liver transaminases  or abdominal imaging [41, 42], a reduction in risk of fibrosis associated with coffee consumption in patients with NAFLD  and NASH , and inverse associations of coffee with ALT and GGT [11, 13–19]. Previous studies of associations of coffee with liver enzymes however have been conducted largely in Japanese male populations or among individuals with alcoholic and viral liver diseases. Investigations of associations of coffee with liver markers in liver diseases of metabolic origin (ie. NAFLD) are scarce. Our findings are consistent with the only other study of healthy individuals conducted in North America  and extend these investigations to a large multi-ethnic population with more complete metabolic characterization, including directly-measured insulin resistance and diabetes status assessed by OGTT. Our assessment of usual, rather than recent, intake of caffeinated and decaffeinated coffee consumption allowed us to better characterize coffee drinking behavior in this population and assess the independent relationships of habitual intake of coffee types with markers of liver injury.
Our findings suggest that caffeinated coffee’s beneficial impact on NAFLD may potentially mediate its inverse association with T2DM. The relationship of NAFLD with IR and oxidative stress is well established . It is possible that the antioxidant properties of coffee may help to offset reactive oxygen species, thereby contributing to a reduction in IR and subsequent fatty acid deposition in the liver [45, 46]. Results of our fourth mechanistic model support this notion as additional adjustment for SI attenuated associations of caffeinated coffee with markers of liver injury. Interestingly, for ALT, the more sensitive and liver-specific transaminase, the association with caffeinated coffee consumption remained significant in this model suggesting that this relationship is not entirely explained by SI.
Investigations of the relationship of coffee consumption with fetuin-A are extremely limited [20, 23]. To our knowledge, only one study has examined the specific relationship of caffeinated coffee with fetuin-A which similarly reported null associations .
Our analyses did not identify significant associations of decaffeinated coffee with markers of liver injury, which is notable considering the inverse relationship of decaffeinated coffee with T2DM was shown to be stronger than that of caffeinated coffee in a meta-analysis . It is likely that the lack of significant associations in the current analysis may be related to the generally modest consumption of decaffeinated coffee and narrow range of intakes in this population. The low level of decaffeinated coffee consumption in the current study may explain why our findings of a null association of decaffeinated coffee with fetuin-A are discordant with the only other study investigating this relationship . Wedick et al. (2011) reported an inverse association of decaffeinated coffee with fetuin-A in their randomized control trial in which individuals consumed 5 6 oz cups of decaffeinated coffee/day . It is possible that higher levels of consumption than reported in the current study are necessary to document this relationship. Alternatively, it is possible that null associations of decaffeinated coffee with liver markers may be explained by differing mechanistic pathways that are being impacted by these coffee types. It has been suggested for example that decaffeinated coffee may specifically benefit beta cell function rather than IR, which appears to be more strongly linked with caffeinated coffee in this cohort .
The current study has a number of strengths including its large sample size, multi-ethnic population and detailed measures used to assess participant characteristics including diabetes status, assessed using the OGTT, and SI determined by FSIGTT. Further, a validated FFQ was used to measure usual intake of coffee and nutritional covariates. Unlike the majority of previous studies, our FFQ differentiated between caffeinated and decaffeinated coffee intake offering insight into the relationship of usual intake of each coffee type with outcome variables. The examination of multiple markers of liver injury, including liver transaminases, fetuin-A, and NAFLD liver fat score, represents a further strength of this study.
Potential limitations of the current study include the cross-sectional design as this precludes conclusions being drawn regarding the causality or the temporal nature of associations. As an observational study, there is potential for residual confounding. However extensive data were collected on potential covariates and these variables were carefully considered in the development of regression models. The reliance on participant self-report and the semi-quantitative design of the FFQ may introduce a degree of misclassification error in estimates of coffee consumption due to imprecise definitions of portion sizes. Such error, however, would likely be non-differential thus resulting in conservative estimates of associations. Additionally, NAFLD was not assessed using the gold standard liver biopsy, as this technique is not suitable for use in epidemiologic studies. Liver transaminase levels have been shown to be well correlated with directly-measured liver fat by magnetic resonance spectroscopy, and the NAFLD liver fat score is a validated, non-invasive method found to be more robust than similar predictive scores for identifying NAFLD, making these techniques appropriate alternatives to liver biopsy [24, 25, 47].
Finally, it should be emphasized that the magnitude of associations of caffeinated coffee with liver markers in the present study was modest, a finding that may be due in part to misclassification of consumption (which would dilute effect estimates, as described above), as well as the relatively modest consumption of coffee in the population, particularly in comparison to other studies investigating associations of coffee with markers of liver injury; the majority of which reported average consumption of 1 or more cups of coffee/day. In the present study, the median (IQR) of caffeinated and decaffeinated coffee intake were 0.4 (0–0.8) servings/day and 0 (0–0.03) servings/day, respectively. Conducting this study in a population in which coffee types were being consumed at higher levels may have improved our ability to detect associations of coffee with outcome variables particularly if coffee’s beneficial effects are exerted only at higher doses. This may be particularly important for decaffeinated coffee since consumption of this type was quite low in this cohort.
Results of this multi-ethnic cross-sectional epidemiologic study identified significant inverse associations of caffeinated coffee consumption with liver transaminases ALT and AST, and the NAFLD liver fat score. There were no significant associations of decaffeinated coffee intake with markers of liver injury. These findings contribute to a greater understanding of the relationship of coffee consumption with liver injury and in turn, suggest that this relationship may mediate the well-established inverse association of caffeinated coffee and the risk of T2DM.
Food frequency questionnaire
Frequently sampled intravenous glucose tolerance test
Insulin resistance atherosclerosis study
Non-alcoholic fatty liver disease
Oral glucose tolerance test
Type 2 diabetes mellitus
Huxley R, Lee CM, Barzi F, Timmermeister L, Czernichow S, Perkovic V, et al. Coffee, decaffeinated coffee, and tea consumption in relation to incident type 2 diabetes mellitus: a systematic review with meta-analysis. Arch Intern Med. 2009;169(22):2053–63.
Jiang X, Zhang D, Jiang W. Coffee and caffeine intake and incidence of type 2 diabetes mellitus: a meta-analysis of prospective studies. Eur J Nutr. 2014;53(1):25–38.
Fraser A, Harris R, Sattar N, Ebrahim S, Davey Smith G, Lawlor DA. Alanine aminotransferase, gamma-glutamyltransferase, and incident diabetes: the British Women’s Heart and Health Study and meta-analysis. Diabetes Care. 2009;32(4):741–50.
Shibata M, Kihara Y, Taguchi M, Tashiro M, Otsuki M. Nonalcoholic fatty liver disease is a risk factor for type 2 diabetes in middle-aged Japanese men. Diabetes Care. 2007;30(11):2940–4.
Kim CH, Park JY, Lee KU, Kim JH, Kim HK. Fatty liver is an independent risk factor for the development of Type 2 diabetes in Korean adults. Diabet Med. 2008;25(4):476–81.
Okamoto M, Takeda Y, Yoda Y, Kobayashi K, Fujino MA, Yamagata Z. The association of fatty liver and diabetes risk. J Epidemiol. 2003;13(1):15–21.
Saab S, Mallam D. Cox GA,2nd, Tong MJ. Impact of coffee on liver diseases: a systematic review. Liver Int. 2014;34(4):495–504.
Masterton GS, Hayes PC. Coffee and the liver: a potential treatment for liver disease? Eur J Gastroenterol Hepatol. 2010;22(11):1277–83.
Muriel P, Arauz J. Coffee and liver diseases. Fitoterapia. 2010;81(5):297–305.
Klatsky AL, Morton C, Udaltsova N, Friedman GD. Coffee, cirrhosis, and transaminase enzymes. Arch Intern Med. 2006;166(11):1190–5.
Ruhl CE, Everhart JE. Coffee and caffeine consumption reduce the risk of elevated serum alanine aminotransferase activity in the United States. Gastroenterology. 2005;128(1):24–32.
Bravi F, Bosetti C, Tavani A, Gallus S, La Vecchia C. Coffee reduces risk for hepatocellular carcinoma: an updated meta-analysis. Clin Gastroenterol Hepatol. 2013;11(11):1413–1421.e1.
Kono S, Shinchi K, Imanishi K, Todoroki I, Hatsuse K. Coffee and serum gamma-glutamyltransferase: a study of self-defense officials in Japan. Am J Epidemiol. 1994;139(7):723–7.
Honjo S, Kono S, Coleman MP, Shinchi K, Sakurai Y, Todoroki I, et al. Coffee consumption and serum aminotransferases in middle-aged Japanese men. J Clin Epidemiol. 2001;54(8):823–9.
Honjo S, Kono S, Coleman MP, Shinchi K, Sakurai Y, Todoroki I, et al. Coffee drinking and serum gamma-glutamyltransferase: an extended study of Self-Defense Officials of Japan. Ann Epidemiol. 1999;9(5):325–31.
Nakanishi N, Nakamura K, Nakajima K, Suzuki K, Tatara K. Coffee consumption and decreased serum gamma-glutamyltransferase: a study of middle-aged Japanese men. Eur J Epidemiol. 2000;16(5):419–23.
Tanaka K, Tokunaga S, Kono S, Tokudome S, Akamatsu T, Moriyama T, et al. Coffee consumption and decreased serum gamma-glutamyltransferase and aminotransferase activities among male alcohol drinkers. Int J Epidemiol. 1998;27(3):438–43.
Yamashita K, Yatsuya H, Muramatsu T, Toyoshima H, Murohara T, Tamakoshi K. Association of coffee consumption with serum adiponectin, leptin, inflammation and metabolic markers in Japanese workers: a cross-sectional study. Nutrition and Diabetes. 2012;2(4):e33.
Xiao Q, Sinha R, Graubard BI, Freedman ND. Inverse associations of total and decaffeinated coffee with liver enzyme levels in National Health and Nutrition Examination Survey 1999–2010. Hepatology. 2014;60:2091–8.
Wedick NM, Brennan AM, Sun Q, Hu FB, Mantzoros CS, van Dam RM. Effects of caffeinated and decaffeinated coffee on biological risk factors for type 2 diabetes: a randomized controlled trial. Nutr J. 2011;10:93.
Ix JH, Biggs ML, Mukamal KJ, Kizer JR, Zieman SJ, Siscovick DS, et al. Association of fetuin-a with incident diabetes mellitus in community-living older adults: the cardiovascular health study. Circulation. 2012;125(19):2316–22.
Mori K, Emoto M, Yokoyama H, Araki T, Teramura M, Koyama H, et al. Association of serum fetuin-A with insulin resistance in type 2 diabetic and nondiabetic subjects. Diabetes Care. 2006;29(2):468.
Jacobs S, Kroger J, Floegel A, Boeing H, Drogan D, Pischon T, et al. Evaluation of various biomarkers as potential mediators of the association between coffee consumption and incident type 2 diabetes in the EPIC-Potsdam Study. Am J Clin Nutr. 2014;100(3):891–900.
Kotronen A, Peltonen M, Hakkarainen A, Sevastianova K, Bergholm R, Johansson LM, et al. Prediction of non-alcoholic fatty liver disease and liver fat using metabolic and genetic factors. Gastroenterology. 2009;137(3):865–72.
Cheung CL, Lam K, Wong I, Cheung B. Non-invasive score identifies ultrasonography-diagnosed non-alcoholic fatty liver disease and predicts mortality in the USA. BMC Med. 2014;12(1):154.
Wagenknecht LE, Mayer EJ, Rewers M, Haffner S, Selby J, Borok GM, et al. The insulin resistance atherosclerosis study (IRAS) objectives, design, and recruitment results. Ann Epidemiol. 1995;5(6):464–72.
World Health Organization. Definition, Diagnosis and Classification of Diabetes Mellitus and its Complications. Part 1: Diagnosis and Classification of Diabetes Mellitus. 1999.
Pacini G, Bergman RN. MINMOD: a computer program to calculate insulin sensitivity and pancreatic responsivity from the frequently sampled intravenous glucose tolerance test. Comput Methods Programs Biomed. 1986;23(2):113–22.
Saad MF, Anderson RL, Laws A, Watanabe RM, Kades WW, Chen YD, et al. A comparison between the minimal model and the glucose clamp in the assessment of insulin sensitivity across the spectrum of glucose tolerance. Insulin Resistance Atherosclerosis Study. Diabetes. 1994;43(9):1114–21.
Steil GM, Volund A, Kahn SE, Bergman RN. Reduced sample number for calculation of insulin sensitivity and glucose effectiveness from the minimal model. Suitability for use in population studies. Diabetes. 1993;42(2):250–6.
Mayer-Davis EJ, D’Agostino Jr R, Karter AJ, Haffner SM, Rewers MJ, Saad M, et al. Intensity and amount of physical activity in relation to insulin sensitivity: the Insulin Resistance Atherosclerosis Study. JAMA. 1998;279(9):669–74.
Mayer-Davis EJ, Vitolins MZ, Carmichael SL, Hemphill S, Tsaroucha G, Rushing J, et al. Validity and reproducibility of a food frequency interview in a Multi-Cultural Epidemiology Study. Ann Epidemiol. 1999;9(5):314–24.
Salvini S, Hunter DJ, Sampson L, Stampfer MJ, Colditz GA, Rosner B, et al. Food-based validation of a dietary questionnaire: the effects of week-to-week variation in food consumption. Int J Epidemiol. 1989;18(4):858–67.
Feskanich D, Rimm EB, Giovannucci EL, Colditz GA, Stampfer MJ, Litin LB, et al. Reproducibility and validity of food intake measurements from a semiquantitative food frequency questionnaire. J Am Diet Assoc. 1993;93(7):790–6.
Loopstra-Masters RC, Liese AD, Haffner SM, Wagenknecht LE, Hanley AJ. Associations between the intake of caffeinated and decaffeinated coffee and measures of insulin sensitivity and beta cell function. Diabetologia. 2011;54(2):320–8.
Liese AD, Schulz M, Moore CG, Mayer-Davis EJ. Dietary patterns, insulin sensitivity and adiposity in the multi-ethnic Insulin Resistance Atherosclerosis Study population. Br J Nutr. 2004;92(6):973–84.
Hanley AJ, Williams K, Festa A, Wagenknecht LE, D’Agostino Jr RB, Kempf J, et al. Elevations in markers of liver injury and risk of type 2 diabetes: the insulin resistance atherosclerosis study. Diabetes. 2004;53(10):2623–32.
Alberti KG, Eckel RH, Grundy SM, Zimmet PZ, Cleeman JI, Donato KA, et al. Harmonizing the metabolic syndrome: a joint interim statement of the International Diabetes Federation Task Force on Epidemiology and Prevention; National Heart, Lung, and Blood Institute; American Heart Association; World Heart Federation; International Atherosclerosis Society; and International Association for the Study of Obesity. Circulation. 2009;120(16):1640–5.
Sattar N, Scherbakova O, Ford I, O’Reilly DS, Stanley A, Forrest E, et al. Elevated alanine aminotransferase predicts new-onset type 2 diabetes independently of classical risk factors, metabolic syndrome, and C-reactive protein in the west of Scotland coronary prevention study. Diabetes. 2004;53(11):2855–60.
Birerdinc A, Stepanova M, Pawloski L, Younossi ZM. Caffeine is protective in patients with non-alcoholic fatty liver disease. Aliment Pharmacol Ther. 2012;35(1):76–82.
Catalano D, Martines GF, Tonzuso A, Pirri C, Trovato FM, Trovato GM. Protective role of coffee in non-alcoholic fatty liver disease (NAFLD). Dig Dis Sci. 2010;55(11):3200–6.
Gutierrez-Grobe Y, Chavez-Tapia N, Sanchez-Valle V, Gavilanes-Espinar JG, Ponciano-Rodriguez G, Uribe M, et al. High coffee intake is associated with lower grade nonalcoholic fatty liver disease: the role of peripheral antioxidant activity. Ann Hepatol. 2012;11(3):350–5.
Anty R, Marjoux S, Iannelli A, Patouraux S, Schneck AS, Bonnafous S, et al. Regular coffee but not espresso drinking is protective against fibrosis in a cohort mainly composed of morbidly obese European women with NAFLD undergoing bariatric surgery. J Hepatol. 2012;57(5):1090–6.
Molloy JW, Calcagno CJ, Williams CD, Jones FJ, Torres DM, Harrison SA. Association of coffee and caffeine consumption with fatty liver disease, nonalcoholic steatohepatitis, and degree of hepatic fibrosis. Hepatology. 2012;55(2):429–36.
Sakaguchi S, Takahashi S, Sasaki T, Kumagai T, Nagata K. Progression of alcoholic and non-alcoholic steatohepatitis: common metabolic aspects of innate immune system and oxidative stress. Drug Metab Pharmacokinet. 2011;26(1):30–46.
Yki-Jarvinen H. Liver fat in the pathogenesis of insulin resistance and type 2 diabetes. Dig Dis. 2010;28(1):203–9.
Westerbacka J, Corner A, Tiikkainen M, Tamminen M, Vehkavaara S, Hakkinen AM, et al. Women and men have similar amounts of liver and intra-abdominal fat, despite more subcutaneous fat in women: implications for sex differences in markers of cardiovascular risk. Diabetologia. 2004;47(8):1360–9.
JCD was supported by the Ministry of Training, Colleges and Universities Ontario Graduate Scholarship, the Canadian Institutes of Health Research Frederick Banting and Charles Best Canada Graduate Scholarships Master’s award, and the Banting and Best Diabetes Centre Novo Nordisk Studentship. AJH holds a Tier II Canada Research Chair in Diabetes Epidemiology. IRAS was supported by grants U01-HL47887, U01-HL47889, U01-HL47890, U01-HL47892, U01-HL47902, DK-29867 and R01-58329 from the National Heart, Lung, and Blood Institute, and by grant M01-RR-43 from the National Institutes of Health. None of the funders had any role in the design, collection, analysis, interpretation, or presentation of the data.
Guarantor of the article: Anthony J. Hanley.
The authors declare that they have no competing interests.
JCD – conducted literature review and analysis, wrote and edited manuscript. ADL – contributed to the discussion and edited manuscript. CL – designed study, contributed to the discussion and edited manuscript. SMH – designed study, contributed to the discussion and edited manuscript. SMW – involved in data analysis, contributed to the discussion and edited manuscript. SJH – involved in data analysis, contributed to the discussion and edited manuscript. JKS – involved in data analysis, contributed to the discussion and edited manuscript. LEW – designed study, contributed to the discussion and edited manuscript. AJH – contributed to the discussion and edited manuscript. All authors read and approved the final version of the manuscript.
Additional file 1: Table S1.
Multiple linear regression analysis of the association of total coffee consumption with markers of liver injury. (DOC 31 kb)
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.
The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.
To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/.
The Creative Commons Public Domain Dedication waiver (https://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Dickson, J.C., Liese, A.D., Lorenzo, C. et al. Associations of coffee consumption with markers of liver injury in the insulin resistance atherosclerosis study. BMC Gastroenterol 15, 88 (2015). https://doi.org/10.1186/s12876-015-0321-3
- Type 2 diabetes
- NAFLD liver fat score
- Liver enzymes
- Insulin resistance atherosclerosis study