Research | Open | Published:
CCL20 is up-regulated in non-alcoholic fatty liver disease fibrosis and is produced by hepatic stellate cells in response to fatty acid loading
Journal of Translational Medicinevolume 16, Article number: 108 (2018)
Nonalcoholic fatty liver disease (NAFLD) is a prevalent complication of extreme obesity. Loading of the liver with fat can progress to inflammation and fibrosis including cirrhosis. The molecular factors involved in the progression from simple steatosis to fibrosis remain poorly understood.
Gene expression profiling using microarray, PCR array, and RNA sequencing was performed on RNA from liver biopsy tissue from patients with extreme obesity. Patients were grouped based on histological findings including normal liver histology with no steatosis, lobular inflammation, or fibrosis, and grades 1, 2, 3, and 4 fibrosis with coexistent steatosis and lobular inflammation. Validation of expression was conducted using quantitative PCR. Serum analysis was performed using ELISA. Expression analysis of hepatocytes and hepatic stellate cells in response to lipid loading were conducted in vitro using quantitative PCR and ELISA.
Three orthogonal methods to profile human liver biopsy RNA each identified the chemokine CCL20 (CC chemokine ligand 20 or MIP-3 alpha) gene as one of the most up-regulated transcripts in NAFLD fibrosis relative to normal histology, validated in a replication group. CCL20 protein levels in serum measured in 224 NAFLD patients were increased in severe fibrosis (p < 0.001), with moderate correlation of hepatic transcript levels and serum levels. Expression of CCL20, but not its cognate receptor CC chemokine receptor 6, was significantly (p < 0.001) increased in response to fatty acid loading in LX-2 hepatic stellate cells, with relative increases greater than those in HepG2 hepatocyte cells.
These results suggest that expression of CCL20, an important inflammatory mediator, is increased in NAFLD fibrosis. CCL20 serves as a chemoattractant molecule for immature dendritic cells, which have been shown to produce many of the inflammatory molecules that mediate liver fibrosis. These data also point to hepatic stellate cells as a key cell type that may respond to lipid loading of the liver.
Obesity is a major risk factor for the development of several metabolic disorders including nonalcoholic fatty liver disease (NAFLD), a group of conditions defined histologically that include lipid accumulation in the liver (steatosis), some degree of inflammation (steatohepatitis), and scarring (fibrosis). The overconsumption of calories leads to storage of triglycerides in adipocytes, as well as in hepatocytes in susceptible individuals. Progression of NAFLD from simple steatosis to hepatic inflammation and liver fibrosis is a complex multi-cellular process involving the interaction of a number of cell types including hepatocytes, Kupffer cells, hepatic stellate cells, sinusoidal endothelial cells as well as infiltrating immune cells . The signals that regulate this complex multicellular process have not yet been well delineated  but likely involve a number of “hits” , including oxidative stress such as lipid peroxidation and mitochondrial dysfunction , production of pro-inflammatory mediators, such as cytokines and chemokines , and the intestinal microbiome producing hepatotoxic substances such as lipopolysaccharide [6, 7].
Acknowledging the molecular and cellular complexity of the pathophysiological processes that are involved in progressing from simple steatosis to more advanced forms of NAFLD , various “omics” approaches have been used including genomic [9, 10], transcriptomic , proteomic , lipidomic , and metabolomics [14, 15]. We applied several orthogonal genome-wide transcriptional profiling methods to RNA isolated from liver biopsies obtained from patients with extreme obesity undergoing bariatric surgery. We used microarray technology, a focused PCR-based array, and RNA sequencing to identify the C-C motif chemokine ligand 20 (CCL20) as a highly up-regulated transcript in NAFLD-associated fibrosis relative to normal liver histology. Circulating CCL20 protein levels in serum were found to be significantly increased in patients with liver fibrosis. To investigate a potential cellular mechanism, CCL20 gene expression was found to be markedly up-regulated in cultured human hepatic stellate cells relative to hepatocyte cells in response to fatty acid loading. CCL20 is the primary ligand for the CC chemokine receptor 6 (CCR6) and serves as a chemoattractant molecule for immature dendritic cells, which have been shown to produce many of the inflammatory molecules that mediate liver fibrosis . CCL20 produced by stellate cells in response to lipid loading may therefore be a key mechanism in the fibrotic progression of NAFLD in response to the increased caloric intake in extreme obesity. These data also suggest that lipid sensing by hepatic stellate cells  may be a novel cellular mechanism in the progression of NAFLD. CCL20 may also be a common pathogenetic mediator in both non-alcoholic and alcoholic liver disease .
Patients and samples
Liver biopsies were obtained from Caucasian individuals enrolled in the Bariatric Surgery Program at the Geisinger Clinic Center for Nutrition and Weight Management . Liver wedge biopsies were obtained intraoperatively as previously described [20,21,22] with a portion fixed in neutral buffered formalin, stained with hematoxylin and eosin, and histologically evaluated as part of clinical standard of care using NASH CRN criteria , as previously described . Patients with serologic, histologic, or other evidence for any chronic liver diseases were excluded from this study. Comprehensive (1–2 h) in-person interviews and psychological evaluations were conducted by certified psychology professionals that included assessment of drug, alcohol, and smoking behaviors to determine clinically significant alcohol intake. Patients determined to have either definite or possible evidence of any addictive behavior were excluded as surgical candidates. In addition, available clinical data, including medical and medication history and diagnostic ICD-9 diagnostic codes, were searched for evidence of drug or alcohol abuse. Samples from patients with evidence of abuse were also excluded from the study. Clinical data were obtained from an obesity database as described previously , and included demographics, clinical measures, ICD-9 codes, medical history, medication use, and common lab results. All study participants provided written informed consent for research, which was conducted according to The Code of Ethics of the World Medical Association (Declaration of Helsinki). The Institutional Review Boards (IRBs) of Geisinger Health System and the Translational Genomics Research Institute approved the research.
RNA isolation and cDNA synthesis
Total RNA was isolated using the RNeasy total RNA isolation kit (Qiagen, Valencia, CA,) according to the manufacturer’s protocol and quantified using the Nanodrop ND-1000 spectrophotometer (Thermo Scientific, Wilmington, DE).
Affymetrix expression profiling
Total RNA was analyzed using the GeneChip human genome U133 plus 2.0 array (Affymetrix, Santa Clara, CA), covering > 47,000 transcripts and variants from the human genome. Due to the built-in redundancy in the microarray design, multiple probe sets mapping to the same transcript yielded a total of 39,000 interrogated loci. All reagents used to make the final hybridization cocktail were from the GeneChip Expression 3′ Amplification One-Cycle Target Labeling and Control Reagents Kit (Affymetrix). Biotin-labeled RNA fragments were produced from a pool of liver RNA with equal amounts from each of nine patients per group pooled based on liver pathology via double-stranded cDNA synthesis followed by in vitro transcription and fragmentation. Hybridization cocktails containing fragmented cRNA, array controls (Affymetrix), bovine serum albumin (BSA), and herring sperm DNA were prepared and hybridized to the array at 45 °C for 16 h. The hybridized arrays were washed, and bound biotin-labeled cRNA detected with a streptavidin–phycoerythrin conjugate. Washing and staining procedures were automated using the Affymetrix Fluidics Station (model 450). Each array was scanned once using the Affymetrix Gene Array Scanner 3000 and the Affymetrix GeneChip Operating Software v1.2 was used to adjust for artifacts, noise, and background. The global method of scaling/normalization was used to for each chip. Fluorescent signal values from the Affymetrix platform were exported to the Spotfire Decision Suite 8.1 software application (https://www.spotfire.com) for analysis. Normalization for chip-to-chip variation and probe-to-probe variability was performed. A Z-score and ANOVA were performed to identify differentially expressed genes.
Total RNA samples were converted to cDNA with the RT2 First Strand Kit (SuperArray Bioscience Corp., Frederick, MD) and examined for cytokine expression level using the Human Inflammatory Cytokines and Receptors PCR Array (Cat. # PAHS-011, SuperArray Bioscience Corp., Frederick, MD) with the RT2 SYBR Green/ROX qPCR Master Mix (SuperArray Bioscience Corp., Frederick, MD) on a 7300 Real-Time PCR System (Applied Biosystems, Foster City, CA) according to the product protocol. Fold change was calculated by determining the ratio of mRNA levels to control values using the Δ threshold cycle (Ct) method (2−ΔΔCt). All data were normalized to GAPDH expression.
Quantitative PCR (qPCR)
Total RNA derived from human tissue samples was reverse-transcribed to cDNA using QuantiTect Rev. Transcription Kit (Qiagen). Quantitative RT-PCR assays were performed in duplicate using designed TaqMan assays for human CCL20 (Hs00171125_m1), CCR6 (Hs00171121_m1), and GAPDH (Hs99999905_m1), all from Life Technologies, in conjunction with the ABI 7500 Fast Real-Time PCR system (Life Technologies, Carlsbad, CA, USA) according to the manufacturer’s protocol. Results were analyzed using the 7500 Software v2.0.6 and DataAssist v3.0 (Life Technologies).
For cell culture analysis, total RNA was extracted from cells with TRIzol reagent (ThermoFisher), and gene expression was measured using commercial TaqMan assays- CCL20 (Hs00355476_m1), CCR6 (Hs01890706_S1) and GAPDH (Hs02758991_g1) (ThermoFisher) and the TaqMan RNA-to-CT 1-step kit and the Step One Plus real–time PCR system. Gene expression levels were normalized to GAPDH and fold-change was determined using the comparative threshold method (ΔΔCT). Two or more replicate experiments were performed.
Serum CCL20 protein levels were measured using commercially available Human CCL20/MIP-3α Colorimetric Sandwich ELISA kits (Cat. # DM3A00, R & D Systems, Minneapolis, MN) according to manufacturer’s protocol (standard curve shown in Additional file 1: Figure S4) using triplicate measurements. Insufficient serum was available for replicate experiments.
HepG2 (ATCC HB-8065) cells were culture in Eagle’s minimum essential medium (ATCC) with 10% FBS. LX-2 cells (SCC064, Millipore) were cultured in DMEM with 10% FBS. Cells were maintained at 37 °C with humidified air and CO2 (5%).
Free fatty acid treatment
Sodium palmitate (2 mmol/l) (P9767, Sigma) was conjugated to fatty acid-free BSA (A7030, Sigma) in a 6:1 molar ratio as described previously . For free fatty acid stimulation, cells were seeded on six-well plates and cultured in complete growth medium to what approximately 80% confluency. Cells were washed twice with pre-warmed phosphate-buffered saline, and treated with various concentrations of free fatty acid diluted into serum–free media, and incubated at 37 °C for 24 h. Experiments were performed in duplicate using triplicate measurements for each experimental condition.
Oil red O staining
Confluent cell monolayers were fixed in phosphate-buffered paraformaldehyde (10%) for 10 min, rinsed with water, followed by 60% isopropyl alcohol, and stained with Oil Red O solution (six parts saturated Oil Red O dye in isopropyl alcohol to four parts water) for 15 min. Stained cells were washed with water, and then incubated for 5 min with 1.0 ml of isopropyl alcohol to dissolve stained oil droplets. The absorbance of the dye-triglyceride complex was measured at 500 nm.
Determination of statistical differences between or among groups was performed using analysis tools within GraphPad Prism 7 (GraphPad Software, Inc., La Jolla, CA). All data are shown as mean ± SD unless otherwise indicated. Statistical significance with p < 0.05 was considered significant with multiple comparison correction performed as indicated. Measurements were performed in triplicate unless otherwise indicated with at least 2 independent experiments.
CCL20 RNA expression is increased in liver RNA from patients with NAFLD fibrosis
A genome-wide strategy was used to identify candidate genes whose expression was dysregulated in NAFLD fibrosis. Hepatic gene expression profiles were generated to identify genes differentially expressed between two histological NAFLD groups. Three pools of three liver RNA samples (nine total samples) from patients with extreme obesity who underwent intra-operative liver biopsy during bariatric surgery who had normal histology and three pools of three liver RNA samples (nine total samples) from patients whose liver histology was classified as grades 1–3 fibrosis (5/9 with bridging fibrosis, Additional file 1: Table S1) were analyzed using the human genome U133 microarray (Affymetrix). This pooling approach provides limited power , thus statistical analysis was not rigorously applied to the data but instead was used to select candidate genes for subsequent validation. Genes whose expression was the most up-regulated on the microarray with signal log ratios of > 3 (Table 1) included two immune related genes, C-C motif chemokine ligand 20 (CCL20) and defensin, alpha 1 (DEFA1), as well as a collagen, type XV, alpha 1 (COL15A1). CCL20 was the most up-regulated gene with a signal log ratio of 5.9.
Based on these data, we then focused analysis on soluble immune-related genes using a PCR-based microarray assay to measure the expression of a panel of 86 human inflammatory cytokine and receptor genes using hepatic RNA samples from a separate set of patients (Additional file 1: Table S2) with severe fibrosis (N = 6) and normal histology (N = 6). CCL20 (Additional file 1: Table S3) was again the most up-regulated gene (23.5-fold). IL8 was the second most up-regulated gene, whose circulating levels have been reported to be increased in patients with extreme obesity and NAFLD fibrosis . CCL20 was also identified as on of the top 20 most up-regulated genes by RNA-sequencing of 142 liver samples from another independent group of participants from the same cohort (manuscript under review).
We then analyzed expression of CCL20 mRNA in an independent group of liver RNA samples with advanced fibrosis using quantitative PCR (QPCR). CCL20 mRNA expression from 32 samples with fibrosis (11 with grade 1, 16 with grade 3 bridging fibrosis, and 5 with grade 4 cirrhosis) was significantly increased (Mann–Whitney two-tailed p < 0.0001) relative to 35 liver samples with normal histology (Fig. 1).
Circulating CCL20 protein levels are increased in patients with NAFLD fibrosis
In order to determine whether liver RNA expression was associated with CCL20 protein expression in vivo, circulating CCL20 levels were determined. CCL20 protein levels in the serum were measured using ELISA in 183 patients (Additional file 1: Table S4) including 106 with normal liver histology, 18 with grade 1 fibrosis, 18 with grade 2 fibrosis, 28 with grade 3 (bridging) fibrosis, and 13 with grade 4 fibrosis (cirrhosis). As shown in Fig. 2, median levels were significantly increased in patients with fibrosis (one-way analysis-of-variance-by-ranks Kruskal–Wallis test p < 0.0001). The CCL20 level of each fibrosis group was also statistically higher (one-way analysis-of-variance-by-ranks Kruskal–Wallis test with Dunn’s adjusted p value p < 0.0001) than the normal group (normal 2.2 pg/ml, grade 1 fibrosis 19.3 pg/ml, grade 2 fibrosis 17.1 pg/ml, grade 3 fibrosis 24.7 pg/ml, grade 4 fibrosis 32.9 pg/ml).
Data were available for 60 patients in whom both CCL20 liver RNA QPCR data and serum protein data were measured. Since correlations can be significantly influenced by outlier values, a robust regression and outlier removal procedure  was performed, which excluding outliers from either the liver RNA expression or CCL20 serum datasets. Correlation was then determined for 44 pairs of samples (Fig. 3). The r2 was 0.104 with a slope that was statistically different from zero (p = 0.032).
CCL20 expression increases in response to fatty acid loading in LX-2 stellate cells
To investigate a potential cellular mechanism for the observed changes in vivo, we selected HepG2 cells as an in vitro model for human hepatocytes and LX-2 cells as an in vitro model for human hepatic stellate cells. Both cell lines have been well characterized in vitro systems for the study of NAFLD . Cells were incubated with physiological levels  of palmitic acid (PA), oleic acid (OA), or a mixture of the two (PA/OA) for 24 h to induce lipid accumulation (Additional file 1: Figures S1 and S2). The highest increases in intracellular lipid staining in both HepG2 and LX-2 cells were observed in response to PA, followed by the PA/OA mixture then OA.
Transcript levels of CCL20 and its cognate receptor CC chemokine receptor 6 (CCR6) were measured in lipid-loaded HepG2 and LX-2 cells using QPCR. Lipopolysaccharide (LPS), known to induce CCL20 expression in primary human hepatic stellate cells  was used as a positive control, which had no effect on lipid loading in either cell type (Additional file 1: Figures S1 and S2). In HepG2 cells (Fig. 4a), CCL20 transcript levels were found to be significantly increased (one-way analysis-of-variance-by-ranks Kruskal–Wallis test p < 0.0001) across treatments. In pair-wise comparisons, CCL20 transcript levels were significantly increased ~fivefold in response to 1 mM oleic acid (Mann–Whitney two-tailed test p = 0.0029), ~twofold in response to 250 μM palmitic acid (Mann–Whitney two-tailed test p < 0.0001), and ~threefold in response to 500 μM palmitic acid (Mann–Whitney two-tailed test p < 0.0001) and the positive control LPS (Mann–Whitney two-tailed test p < 0.0001). In contrast, in LX-2 cells (Fig. 4b), CCL20 levels increased ~ 13-fold with 250 μM (Mann–Whitney two-tailed test p < 0.0001) and the positive control LPS (Mann–Whitney two-tailed test p < 0.0001), ~ 15-fold with 500 μM of palmitic acid (Mann–Whitney two-tailed test p < 0.0001), and ~twofold (Mann–Whitney two-tailed test p < 0.0001) and ~eightfold (Mann–Whitney two-tailed test p < 0.0001) in response to oleic acid or 1 mM palmitic acid/oleic acid (50/50) mixture, respectively. No statistically significant differences were observed in CCR6 transcript levels across the different treatments in HepG2 cells (Fig. 4c) after adjusting for multiple comparisons. In LX-2 cells (Fig. 4d), the only pairwise comparison to achieve statistical significance was ~ 50% decrease in CCR6 transcript levels in response to 500 μM oleic acid (Mann–Whitney two-tailed test p = 0.041).
To determine if the increases in transcription of CCL20 resulted in increased CCL20 protein, ELISA was used to measure CCL20 in conditioned culture medium following 24 and 48 h of fatty acid or LPS treatment of LX-2 cells (Additional file 1: Figure S3). Approximately 50 pmol/ml and 22 pmol/nl of secreted CCL20 was observed in response to 500 μM palmitic acid and 250 μM palmitic acid, respectively, with almost 75 pmol/ml measured after LPS treatment. In response to treatment with a 50/50 palmitic acid/oleic acid mixture, a 1 mM concentration yielded approximately 16 pmol/ml protein, while 500 μM mixture produced < 2 pmol/ml. Treatment with oleic acid alone had no effect on CCL20 protein levels, similar to the results observed for the transcript measurements.
The initial hepatic steatosis that characterizes NAFLD results primarily from the storage of free fatty acids as triglycerides derived from plasma via the diet, de novo lipogenesis, and adipose tissue lipolysis, thus patients with extreme obesity are at particularly high risk . Increased expression of the chemokine CCL20, which belongs to a class of molecules that serve as chemoattractants for the infiltration of immune cells to the injured liver , was associated with NAFLD fibrosis. CCL20 is a key molecule involved in the homing of dendritic cells, which have been shown to play a key role in liver fibrosis in animal models . Our results extend initial observations of an approximate tenfold increase in expression of hepatic CCL20 transcript levels in eight samples from patients with NASH  to a large cohort of patients with extreme obesity. CCL20 also appears to play a key role in alcoholic liver disease , as well an emerging role in viral hepatitis [35, 36]. These data suggest that CCL20 may be a central regulatory molecule in a pathway common to the development of liver injury and fibrosis in alcoholic and non-alcoholic, as well as viral-mediated, liver disease.
CCL20, also known as macrophage inflammatory protein (MIP-3α) , serves as a ligand for its only receptor, CC chemokine receptor 6 (CCR6) . In studies of alcoholic liver disease closely paralleling those reported here, CCL20 was found to be the most up-regulated chemokine gene by microarray gene expression profiling of human liver biopsy tissue . This was followed by the findings that CCL20 hepatic RNA expression and serum levels were increased in patients with alcoholic hepatitis and were associated grade of hepatic fibrosis . In that study, hepatic stellate cells were found to be the primary cell type producing CCL20, which was shown to induce a profibrogenic and proinflammatory expression profile in primary human hepatic stellate cells exposed to CCL20 in vitro. This implies a potential dual role for CCL20, one as a chemoattractant for immune cells, the other as an autocrine factor driving the fibrogenic process.
These data also point to HSCs as key cell type orchestrating liver damage. Induction of CCL20 in response to fatty acid loading suggests that HSCs may serve as a key lipid-sensing function in NAFLD, in addition to the presumed central role of the hepatocyte, in order to recruit dendritic and other immune cells as part of an inflammatory response to lipotoxicity. HSCs store triglycerides, as well as retinol, and express adipocyte-like features . Further implicating stellate cells in lipid sensing is the association of a variant in the PNPLA3 gene, discovered through genome-wide association studies of NAFLD and associated with risk of NAFLD fibrosis , with increased lipid droplet content in primary HSCs and LX-2 cells, as well as with the increased production of several cytokines [42, 43]. Production of CCL20 by HSCs may therefore play a role in the mechanism of PNPLA3-mediated NAFLD.
An important aspect of our results, similar to those for the role of CCL20 in alcoholic liver disease, is the use of ex vivo human samples and cultured human cell lines to define the transcriptional regulation of the CCL20 gene in HSCs, which has not yet been well characterized in vivo . Studies on humans are confounded by an uncontrolled environment, behavioral and genetic heterogeneity, and other variables. However, often used mouse models for either alcoholic or non-alcoholic liver disease lack physiological fidelity to humans. For example, available animal models do not re-capitulate the pathophysiology of severe alcoholic hepatitis that is found in humans . Hepatic gene expression in nine NAFLD mouse models were found to have little overlap with transcriptional changes in human liver biopsies . Even a dietary obese mouse model purporting a high level of fidelity to human NASH fibrosis was characterized by increased serum HDL levels and decreased triglyceride levels , the opposite of that found in the metabolic syndrome in humans , raising concerns about the potential contribution of these starkly different metabolic effects on physiological phenotype. Ambient temperature also appears to play a significant role in the development of NAFLD in mice . Future studies on CCL20 using animal models of NAFLD may therefore be problematic due to these species-specific differences, thus other approaches, such as single cell RNA sequencing and multi-cellular in vitro models that based on human samples may be needed.
The results we obtained from a cohort of primarily females with extreme obesity who are seeking bariatric surgery may not be generalizable to individuals with lower levels of obesity where the pathophysiology of NAFLD may be different. However, the population with this degree of obesity constitutes over 6% of the US population and has been increasing . The setting of bariatric surgery provides a unique opportunity to obtain liver biopsy tissue in the absence of suspected liver disease, which is largely the case for outpatient-based studies of liver tissue for which a clinical indication for biopsy must exist. Obtaining liver biopsies from individuals without suspicion of clinical disease, i.e., healthy subjects, is also more difficult due to the potential morbidity and mortality associated with the procedure. It is also possible that NAFLD in extreme obesity, which appears to be largely a clinically silent disease , may have a different pathophysiology than NAFLD in cohorts with lesser degrees of obesity. We also acknowledge that the relative homogeneity of the cohort, primarily of Caucasian ethnicity and from a single institution, and the retrospective design, may limit the applicability of the findings to other populations.
CCL20 represents a new candidate molecule that may be involved in the pathogenesis of NAFLD fibrosis in the population with obesity that was investigated. Future studies on more heterogeneous populations will be needed in order to determine whether our findings are generalizable. As a soluble inflammatory mediator, CCL20 may also represent a potential diagnostic and/or therapeutic target.
nonalcoholic fatty liver disease
C-C motif chemokine ligand 20
CC chemokine receptor 6
Wynn TA. Cellular and molecular mechanisms of fibrosis. J Pathol. 2008;214(2):199–210.
Gluchowski NL, Becuwe M, Walther TC, Farese RV Jr. Lipid droplets and liver disease: from basic biology to clinical implications. Nat Rev Gastroenterol Hepatol. 2017;14(6):343–55.
Basaranoglu M, Basaranoglu G, Senturk H. From fatty liver to fibrosis: a tale of “second hit”. World J Gastroenterol. 2013;19(8):1158–65.
Tariq Z, Green CJ, Hodson L. Are oxidative stress mechanisms the common denominator in the progression from hepatic steatosis towards non-alcoholic steatohepatitis (NASH)? Liver Int. 2014;34(7):e180–90.
Marra F, Tacke F. Roles for chemokines in liver disease. Gastroenterology. 2014;147(3):577–94.
Wieland A, Frank DN, Harnke B, Bambha K. Systematic review: microbial dysbiosis and nonalcoholic fatty liver disease. Aliment Pharmacol Ther. 2015;42(9):1051–63.
Arslan N. Obesity, fatty liver disease and intestinal microbiota. World J Gastroenterol. 2014;20(44):16452–63.
Haas JT, Francque S, Staels B. Pathophysiology and mechanisms of nonalcoholic fatty liver disease. Annu Rev Physiol. 2016;78:181–205.
Anstee QM, Seth D, Day CP. Genetic factors that affect risk of alcoholic and non-alcoholic fatty liver disease. Gastroenterology. 2016;150(8):1728–44.
Kahali B, Halligan B, Speliotes EK. Insights from genome-wide association analyses of nonalcoholic fatty liver disease. Semin Liver Dis. 2015;35(4):375–91.
Zhu R, Baker SS, Moylan CA, Abdelmalek MF, Guy CD, Zamboni F, Wu D, Lin W, Liu W, Baker RD, Govindarajan S, Cao Z, Farci P, Diehl AM, Zhu L. Systematic transcriptome analysis reveals elevated expression of alcohol-metabolizing genes in NAFLD livers. J Pathol. 2016;238(4):531–42.
Miller MH, Walsh SV, Atrih A, Huang JT, Ferguson MA, Dillon JF. Serum proteome of nonalcoholic fatty liver disease: a multimodal approach to discovery of biomarkers of nonalcoholic steatohepatitis. J Gastroenterol Hepatol. 2014;29(10):1839–47.
Stiuso P, Scognamiglio I, Murolo M, Ferranti P, De Simone C, Rizzo MR, Tuccillo C, Caraglia M, Loguercio C, Federico A. Serum oxidative stress markers and lipidomic profile to detect NASH patients responsive to an antioxidant treatment: a pilot study. Oxid Med Cell Longev. 2014;2014:169216.
Rodriguez-Gallego E, Guirro M, Riera-Borrull M, Hernandez-Aguilera A, Marine-Casado R, Fernandez-Arroyo S, Beltran-Debon R, Sabench F, Hernandez M, del Castillo D, Menendez JA, Camps J, Ras R, Arola L, Joven J. Mapping of the circulating metabolome reveals alpha-ketoglutarate as a predictor of morbid obesity-associated non-alcoholic fatty liver disease. Int J Obes (Lond). 2015;39(2):279–87.
von Schonfels W, Patsenker E, Fahrner R, Itzel T, Hinrichsen H, Brosch M, Erhart W, Gruodyte A, Vollnberg B, Richter K, Landrock A, Schreiber S, Bruckner S, Beldi G, Sipos B, Becker T, Rocken C, Teufel A, Stickel F, Schafmayer C, Hampe J. Metabolomic tissue signature in human non-alcoholic fatty liver disease identifies protective candidate metabolites. Liver Int. 2015;35(1):207–14.
Connolly MK, Bedrosian AS, Clair JMS, Mitchell AP, Ibrahim J, Stroud A, Pachter HL, Bar-Sagi D, Frey AB, Miller G. In liver fibrosis, dendritic cells govern hepatic inflammation in mice via TNF-alpha. J Clin Investig. 2009;119(11):3213–25.
Hetherington AM, Sawyez CG, Zilberman E, Stoianov AM, Robson DL, Borradaile NM. Differential lipotoxic effects of palmitate and oleate in activated human hepatic stellate cells and epithelial hepatoma cells. Cell Physiol Biochem. 2016;39(4):1648–62.
Gao B, Xu M. Chemokines and alcoholic hepatitis: are chemokines good therapeutic targets? Gut. 2014;63(11):1683–4.
Wood GC, Chu X, Manney C, Strodel W, Petrick A, Gabrielsen J, Seiler J, Carey D, Argyropoulos G, Benotti P, Still CD, Gerhard GS. An electronic health record-enabled obesity database. BMC Med Inform Decis Mak. 2012;12(1):45.
Leti F, Malenica I, Doshi M, Courtright A, Van Keuren-Jensen K, Legendre C, Still CD, Gerhard GS, DiStefano JK. High-throughput sequencing reveals altered expression of hepatic microRNAs in nonalcoholic fatty liver disease-related fibrosis. Transl Res. 2015;166(3):304–14.
DiStefano JK, Kingsley C, Wood GC, Chu X, Argyropoulos G, Still CD, Doné SC, Legendre C, Tembe W, Gerhard GS. Genome-wide analysis of hepatic lipid content in extreme obesity. Acta Diabetol. 2015;52(2):373–82.
Gerhard GS, Benotti P, Wood GC, Chu X, Argyropoulos G, Petrick A, Strodel WE, Gabrielsen JD, Ibele A, Still CD, Kingsley C, DiStefano J. Identification of novel clinical factors associated with hepatic fat accumulation in extreme obesity. J Obes. 2014;2014:368210.
Kleiner DE, Brunt EM, Van Natta M, Behling C, Contos MJ, Cummings OW, Ferrell LD, Liu YC, Torbenson MS, Unalp-Arida A, Yeh M, McCullough AJ, Sanyal AJ. Design and validation of a histological scoring system for nonalcoholic fatty liver disease. Hepatology. 2005;41(6):1313–21.
DiStefano JK, Kingsley C, Wood GC, Chu X, Argyropoulos G, Still CD, Doné SC, Legendre C, Tembe W, Gerhard GS. Genome-wide analysis of hepatic lipid content in extreme obesity. Acta Diabetol. 2014;52(2):373–82.
Wood GC, Chu X, Manney C, Strodel W, Petrick A, Gabrielsen J, Seiler J, Carey D, Argyropoulos G, Benotti P, Still CD, Gerhard GS. An electronic health record-enabled obesity database. BMC Med Inform Decis Mak. 2012;12:45.
Readnower RD, Brainard RE, Hill BG, Jones SP. Standardized bioenergetic profiling of adult mouse cardiomyocytes. Physiol Genom. 2012;44(24):1208–13.
Peng X, Wood CL, Blalock EM, Chen KC, Landfield PW, Stromberg AJ. Statistical implications of pooling RNA samples for microarray experiments. BMC Bioinform. 2003;4:26.
du Plessis J, Korf H, van Pelt J, Windmolders P, Vander Elst I, Verrijken A, Hubens G, Van Gaal L, Cassiman D, Nevens F, Francque S, van der Merwe S. Pro-inflammatory cytokines but not endotoxin-related parameters associate with disease severity in patients with NAFLD. PLoS ONE. 2016;11(12):e0166048.
Motulsky HJ, Brown RE. Detecting outliers when fitting data with nonlinear regression—a new method based on robust nonlinear regression and the false discovery rate. BMC Bioinform. 2006;7:123.
Kanuri G, Bergheim I. In vitro and in vivo models of non-alcoholic fatty liver disease (NAFLD). Int J Mol Sci. 2013;14(6):11963–80.
Affo S, Morales-Ibanez O, Rodrigo-Torres D, Altamirano J, Blaya D, Dapito DH, Millan C, Coll M, Caviglia JM, Arroyo V, Caballeria J, Schwabe RF, Gines P, Bataller R, Sancho-Bru P. CCL20 mediates lipopolysaccharide induced liver injury and is a potential driver of inflammation and fibrosis in alcoholic hepatitis. Gut. 2014;63(11):1782–92.
Petrick A, Benotti P, Wood GC, Still CD, Strodel WE, Gabrielsen J, Rolston D, Chu X, Argyropoulos G, Ibele A, Gerhard GS. Utility of ultrasound, transaminases, and visual inspection to assess nonalcoholic fatty liver disease in bariatric surgery patients. Obes Surg. 2015;25(12):2368–75.
Wasmuth HE, Tacke F, Trautwein C. Chemokines in liver inflammation and fibrosis. Semin Liver Dis. 2010;30(3):215–25.
Lukacs-Kornek V, Schuppan D. Dendritic cells in liver injury and fibrosis: shortcomings and promises. J Hepatol. 2013;59(5):1124–6.
Xin J, Ding W, Hao S, Chen X, Zhang J, Jiang L, Zhou Q, Shi D, Zhang L, Xu X, Cao H, Li L, Li J. Serum macrophage inflammatory protein 3alpha levels predict the severity of HBV-related acute-on-chronic liver failure. Gut. 2016;65(2):355–7.
Liu Y, Li L, Liu J, She WM, Shi JM, Li J, Wang JY, Jiang W. Activated hepatic stellate cells directly induce pathogenic Th17 cells in chronic hepatitis B virus infection. Exp Cell Res. 2017;359(1):129–37.
Hieshima K, Imai T, Opdenakker G, Van Damme J, Kusuda J, Tei H, Sakaki Y, Takatsuki K, Miura R, Yoshie O, Nomiyama H. Molecular cloning of a novel human CC chemokine liver and activation-regulated chemokine (LARC) expressed in liver. Chemotactic activity for lymphocytes and gene localization on chromosome 2. J Biol Chem. 1997;272(9):5846–53.
Schutyser E, Struyf S, Van Damme J. The CC chemokine CCL20 and its receptor CCR6. Cytokine Growth Factor Rev. 2003;14(5):409–26.
Affo S, Dominguez M, Lozano JJ, Sancho-Bru P, Rodrigo-Torres D, Morales-Ibanez O, Moreno M, Millan C, Loaeza-del-Castillo A, Altamirano J, Garcia-Pagan JC, Arroyo V, Gines P, Caballeria J, Schwabe RF, Bataller R. Transcriptome analysis identifies TNF superfamily receptors as potential therapeutic targets in alcoholic hepatitis. Gut. 2013;62(3):452–60.
Sauvant P, Cansell M, Atgie C. Vitamin A and lipid metabolism: relationship between hepatic stellate cells (HSCs) and adipocytes. J Physiol Biochem. 2011;67(3):487–96.
Rotman Y, Koh C, Zmuda JM, Kleiner DE, Liang TJ, Nash CRN. The association of genetic variability in patatin-like phospholipase domain-containing protein 3 (PNPLA3) with histological severity of nonalcoholic fatty liver disease. Hepatology. 2010;52(3):894–903.
Pirazzi C, Valenti L, Motta BM, Pingitore P, Hedfalk K, Mancina RM, Burza MA, Indiveri C, Ferro Y, Montalcini T, Maglio C, Dongiovanni P, Fargion S, Rametta R, Pujia A, Andersson L, Ghosal S, Levin M, Wiklund O, Iacovino M, Boren J, Romeo S. PNPLA3 has retinyl-palmitate lipase activity in human hepatic stellate cells. Hum Mol Genet. 2014;23(15):4077–85.
Bruschi FV, Claudel T, Tardelli M, Caligiuri A, Stulnig TM, Marra F, Trauner M. The PNPLA3 I148M variant modulates the fibrogenic phenotype of human hepatic stellate cells. Hepatology. 2017;65(6):1875–90.
Zhao L, Xia J, Wang X, Xu F. Transcriptional regulation of CCL20 expression. Microbes Infect. 2014;16(10):864–70.
Teufel A, Itzel T, Erhart W, Brosch M, Wang XY, Kim YO, von Schonfels W, Herrmann A, Bruckner S, Stickel F, Dufour JF, Chavakis T, Hellerbrand C, Spang R, Maass T, Becker T, Schreiber S, Schafmayer C, Schuppan D, Hampe J. Comparison of gene expression patterns between mouse models of nonalcoholic fatty liver disease and liver tissues from patients. Gastroenterology. 2016;151(3):513–25.
Krishnan A, Abdullah TS, Mounajjed T, Hartono S, McConico A, White T, LeBrasseur N, Lanza I, Nair S, Gores G, Charlton M. A longitudinal study of whole body, tissue, and cellular physiology in a mouse model of fibrosing NASH with high fidelity to the human condition. Am J Physiol Gastrointest Liver Physiol. 2017;312(6):G666–80.
Kaur J. A comprehensive review on metabolic syndrome. Cardiol Res Pract. 2014;2014:943162.
Giles DA, Moreno-Fernandez ME, Stankiewicz TE, Graspeuntner S, Cappelletti M, Wu D, Mukherjee R, Chan CC, Lawson MJ, Klarquist J, Sunderhauf A, Softic S, Kahn CR, Stemmer K, Iwakura Y, Aronow BJ, Karns R, Steinbrecher KA, Karp CL, Sheridan R, Shanmukhappa SK, Reynaud D, Haslam DB, Sina C, Rupp J, Hogan SP, Divanovic S. Thermoneutral housing exacerbates nonalcoholic fatty liver disease in mice and allows for sex-independent disease modeling. Nat Med. 2017;23(7):829–38.
Sturm R, Hattori A. Morbid obesity rates continue to rise rapidly in the United States. Int J Obes (Lond). 2013;37(6):889–91.
Study concept and design: GSG; acquisition of data: XC, QJ, HC, ATP, WS, JG, and FL; analysis and interpretation of data: CDS, DC, TM, JKD, and GSG; drafting of the manuscript: JKD and GSG; critical revision of the manuscript for important intellectual content: XC, QJ, GCW, PB, ATP, WS, JG, CDS, DC, TM, JKD, and GSG; Statistical analysis: GCW and GSG; obtained funding: GSG, CDS, and JKD. All authors read and approved the final manuscript.
Dr. C. Still receives grant and consulting support from Ethicon Endo-Surgery; Dr. Benotti receives consulting support from Ethicon Endo-Surgery. No other authors have any potential competing interests to disclose.
Availability of data and materials
All data generated or analyzed were included in this manuscript and its Additional file. For any additional information, please contact the corresponding author.
Consent for publication
The manuscript does not contain any personal health identifiers. The data provided dataset are anonymous.
Ethics approval and consent to participate
All study participants provided written informed consent for research, which was conducted according to The Code of Ethics of the World Medical Association (Declaration of Helsinki). The Institutional Review Boards (IRBs) of Geisinger Health System and the Translational Genomics Research Institute approved the research.
This work was supported by funding from the National Institutes of Health grant DK091601 (GSG and JKD), the Joseph & Rebecca Goodfriend Endowed Chair In Genetics at Temple University, and the Geisinger Obesity Research Institute.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.