Holistic Characterization of Single Hepatocyte Transcriptome Responses to High Fat Diet

During nutritional overload and obesity, hepatocyte function is grossly altered, and a subset of hepatocytes begins to accumulate fat droplets, leading to non-alcoholic fatty liver disease (NAFLD). Recent single cell studies revealed how non-parenchymal cells, such as macrophages, hepatic stellate cells, and endothelial cells, heterogeneously respond to NAFLD. However, it remains to be characterized how hepatocytes, the major constituents of the liver, respond to nutritional overload in NAFLD. Here, using droplet-based single cell RNA-sequencing (Drop-seq), we characterized how the transcriptomic landscape of individual hepatocytes is altered in response to high-fat diet (HFD) and NAFLD. We showed that the entire hepatocytes population undergoes substantial transcriptome changes upon HFD, although the patterns of alteration were highly heterogeneous with zonation-dependent and –independent effects. Periportal (zone 1) hepatocytes downregulated many zone 1-specific marker genes, while a small number of genes mediating gluconeogenesis were upregulated. Pericentral (zone 3) hepatocytes also downregulated many zone 3-specific genes; however, they upregulated several genes that promote HFD-induced fat droplet formation, consistent with findings that zone 3 hepatocytes accumulate more lipid droplets. Zone 3 hepatocytes also upregulated ketogenic pathways as an adaptive mechanism to HFD. Interestingly, many of the top HFD-induced genes, which encode proteins regulating lipid metabolism, were strongly co-expressed with each other in a subset of hepatocytes, producing a variegated pattern of spatial co-localization that is independent of metabolic zonation. In conclusion, our dataset provides a useful resource for understanding hepatocellular alteration during NAFLD at single cell level.

[1]  T. Baumert,et al.  Single-cell genomics and spatial transcriptomics: Discovery of novel cell states and cellular interactions in liver physiology and disease biology , 2020, Journal of hepatology.

[2]  B. Blagoev,et al.  Multi-omics characterization of a diet-induced obese model of non-alcoholic steatohepatitis , 2020, Scientific Reports.

[3]  C. Ponting,et al.  Single-Cell Transcriptomics Uncovers Zonation of Function in the Mesenchyme during Liver Fibrosis , 2019, Cell reports.

[4]  C. Ponting,et al.  Resolving the fibrotic niche of human liver cirrhosis at single cell level , 2019, Nature.

[5]  Jiandie D. Lin,et al.  Landscape of Intercellular Crosstalk in Healthy and NASH Liver Revealed by Single-Cell Secretome Gene Analysis. , 2019, Molecular cell.

[6]  Paul J. Hoffman,et al.  Comprehensive Integration of Single-Cell Data , 2018, Cell.

[7]  F. Tacke,et al.  Single Cell RNA Sequencing Identifies Subsets of Hepatic Stellate Cells and Myofibroblasts in Liver Fibrosis , 2019, Cells.

[8]  Lai Guan Ng,et al.  Dimensionality reduction for visualizing single-cell data using UMAP , 2018, Nature Biotechnology.

[9]  Christoph Hafemeister,et al.  Comprehensive integration of single cell data , 2018, bioRxiv.

[10]  Gary D Bader,et al.  Single cell RNA sequencing of human liver reveals distinct intrahepatic macrophage populations , 2018, Nature Communications.

[11]  A. Saltiel,et al.  Lipotoxicity induces hepatic protein inclusions through TANK binding kinase 1–mediated p62/sequestosome 1 phosphorylation , 2018, Hepatology.

[12]  Kevin R. Moon,et al.  Recovering Gene Interactions from Single-Cell Data Using Data Diffusion , 2018, Cell.

[13]  Nancy R. Zhang,et al.  SAVER: Gene expression recovery for single-cell RNA sequencing , 2018, Nature Methods.

[14]  Dongmin Jung,et al.  ShinyGO: a graphical enrichment tool for ani-mals and plants , 2018, bioRxiv.

[15]  D. Lauffenburger,et al.  Hepatic Dysfunction Caused by Consumption of a High-Fat Diet. , 2017, Cell reports.

[16]  L. Campbell,et al.  G0S2: A small giant controller of lipolysis and adipose-liver fatty acid flux. , 2017, Biochimica et biophysica acta. Molecular and cell biology of lipids.

[17]  J. C. Love,et al.  Seq-Well: portable, low-cost RNA sequencing of single cells at high throughput. , 2017, Nature methods.

[18]  B. Wang,et al.  Ablation of cytochrome P450 omega-hydroxylase 4A14 gene attenuates hepatic steatosis and fibrosis , 2017, Proceedings of the National Academy of Sciences.

[19]  S. Davies,et al.  Lipid zonation and phospholipid remodeling in nonalcoholic fatty liver disease , 2017, Hepatology.

[20]  I. Amit,et al.  Single-cell spatial reconstruction reveals global division of labor in the mammalian liver , 2016, Nature.

[21]  L. Grøntved,et al.  High fat diet-induced changes of mouse hepatic transcription and enhancer activity can be reversed by subsequent weight loss , 2016, Scientific Reports.

[22]  J. McManaman,et al.  Perilipin-2 Deletion Impairs Hepatic Lipid Accumulation by Interfering with Sterol Regulatory Element-binding Protein (SREBP) Activation and Altering the Hepatic Lipidome* , 2016, The Journal of Biological Chemistry.

[23]  W. Birchmeier.  Orchestrating Wnt signalling for metabolic liver zonation , 2016, Nature Cell Biology.

[24]  A. D. Jones,et al.  Liver-specific loss of Perilipin 2 alleviates diet-induced hepatic steatosis, inflammation, and fibrosis. , 2016, American journal of physiology. Gastrointestinal and liver physiology.

[25]  D. Erion,et al.  Hepatocyte-Specific Disruption of CD36 Attenuates Fatty Liver and Improves Insulin Sensitivity in HFD-Fed Mice. , 2016, Endocrinology.

[26]  B. Stanger.  Probing hepatocyte heterogeneity , 2015, Cell Research.

[27]  K. Deisseroth,et al.  Hybrid Periportal Hepatocytes Regenerate the Injured Liver without Giving Rise to Cancer , 2015, Cell.

[28]  R. Nusse,et al.  Self-renewing diploid Axin2+ cells fuel homeostatic renewal of the liver , 2015, Nature.

[29]  Evan Z. Macosko,et al.  Highly Parallel Genome-wide Expression Profiling of Individual Cells Using Nanoliter Droplets , 2015, Cell.

[30]  S. Schuster,et al.  Zonation of hepatic fatty acid metabolism - The diversity of its regulation and the benefit of modeling. , 2015, Biochimica et biophysica acta.

[31]  S. Enosawa,et al.  Improvement of Infusion Process in Cell Transplantation: Effect of Shear Stress on Hepatocyte Viability Under Horizontal and Vertical Syringe Orientation. , 2015, Cell medicine.

[32]  X. Hou,et al.  Adipogenic changes of hepatocytes in a high-fat diet-induced fatty liver mice model and non-alcoholic fatty liver disease patients , 2015, Endocrine.

[33]  Yoshihisa Takahashi,et al.  Histopathology of nonalcoholic fatty liver disease/nonalcoholic steatohepatitis. , 2014, World journal of gastroenterology.

[34]  T. Auguet,et al.  Molecular pathways in non-alcoholic fatty liver disease , 2014, Clinical and experimental gastroenterology.

[35]  S. Kersten.  Integrated physiology and systems biology of PPARα , 2014, Molecular metabolism.

[36]  Melvin E Andersen,et al.  A map of the PPARα transcription regulatory network for primary human hepatocytes. , 2014, Chemico-biological interactions.

[37]  T. Czyzyk,et al.  Targeted Disruption of G0/G1 Switch Gene 2 Enhances Adipose Lipolysis, Alters Hepatic Energy Balance, and Alleviates High-Fat Diet–Induced Liver Steatosis , 2014, Diabetes.

[38]  Brenda S. Hijmans,et al.  Zonation of glucose and fatty acid metabolism in the liver: mechanism and metabolic consequences. , 2014, Biochimie.

[39]  A. Diehl,et al.  NAFLD, NASH and liver cancer , 2013, Nature Reviews Gastroenterology &Hepatology.

[40]  S. Wölfl,et al.  Zonation of Nitrogen and Glucose Metabolism Gene Expression upon Acute Liver Damage in Mouse , 2013, PloS one.

[41]  Eliseo Guallar,et al.  Prevalence of nonalcoholic fatty liver disease in the United States: the Third National Health and Nutrition Examination Survey, 1988-1994. , 2013, American journal of epidemiology.

[42]  Thomas R. Gingeras,et al.  STAR: ultrafast universal RNA-seq aligner , 2013, Bioinform..

[43]  Guangchuang Yu,et al.  clusterProfiler: an R package for comparing biological themes among gene clusters. , 2012, Omics : a journal of integrative biology.

[44]  S. Colnot,et al.  Transcription dynamics in a physiological process: β-catenin signaling directs liver metabolic zonation. , 2011, The international journal of biochemistry & cell biology.

[45]  O. Cummings,et al.  Relationship of steatosis grade and zonal location to histological features of steatohepatitis in adult patients with non-alcoholic fatty liver disease. , 2008, Journal of hepatology.

[46]  Geoffrey E. Hinton,et al.  Visualizing Data using t-SNE , 2008 .

[47]  J. Dyck,et al.  Increased Hepatic CD36 Expression Contributes to Dyslipidemia Associated With Diet-Induced Obesity , 2007, Diabetes.

[48]  B. Jayaprakasam,et al.  Regulation of hepatic fatty acid elongase and desaturase expression in diabetes and obesity Published, JLR Papers in Press, June 21, 2006. , 2006, Journal of Lipid Research.

[49]  S. Kersten,et al.  Peroxisome proliferator-activated receptor alpha mediates the effects of high-fat diet on hepatic gene expression. , 2006, Endocrinology.

[50]  B. Ahrén,et al.  The high-fat diet-fed mouse: a model for studying mechanisms and treatment of impaired glucose tolerance and type 2 diabetes. , 2004, Diabetes.

[51]  Insuk Sohn,et al.  Hepatic gene expression profiles in a long-term high-fat diet-induced obesity mouse model. , 2004, Gene.

[52]  Ling-jia Wang,et al.  Loss of cell anchorage triggers apoptosis (anoikis) in primary mouse hepatocytes. , 2002, Molecular genetics and metabolism.

[53]  N. Katz,et al.  Metabolic heterogeneity of hepatocytes across the liver acinus. , 1992, The Journal of nutrition.

[54]  K. Jungermann,et al.  Functional specialization of different hepatocyte populations. , 1989, Physiological reviews.

[55]  S. Thorgeirsson,et al.  In situ hybridization studies on expression of albumin and alpha-fetoprotein during the early stage of neoplastic transformation in rat liver. , 1987, Cancer research.

[56]  K. Jungermann,et al.  Functional Hepatocellular Heterogeneity , 2007, Hepatology.