- Open Access
Insulin treatment improves liver histopathology and decreases expression of inflammatory and fibrogenic genes in a hyperglycemic, dyslipidemic hamster model of NAFLD
Journal of Translational Medicine volume 19, Article number: 80 (2021)
Non-alcoholic fatty liver disease (NAFLD) and non-alcoholic steatohepatitis (NASH) are highly prevalent comorbidities in patients with Type 2 diabetes. While many of these patients eventually will need treatment with insulin, little is known about the effects of insulin treatment on histopathological parameters and hepatic gene expression in diabetic patients with co-existing NAFLD and NASH. To investigate this further, we evaluated the effects of insulin treatment in NASH diet-fed hamsters with streptozotocin (STZ) -induced hyperglycemia.
Forty male Syrian hamsters were randomized into four groups (n = 10/group) receiving either a NASH-inducing (high fat, fructose and cholesterol) or control diet (CTRL) for four weeks, after which they were treated with STZ or sham-injected and from week five treated with either vehicle (CTRL, NASH, NASH-STZ) or human insulin (NASH-STZ-HI) for four weeks by continuous s.c. infusion via osmotic minipumps.
NASH-STZ hamsters displayed pronounced hyperglycemia, dyslipidemia and more severe liver pathology compared to both CTRL and NASH groups. Insulin treatment attenuated dyslipidemia in NASH-STZ-HI hamsters and liver pathology was considerably improved compared to the NASH-STZ group, with prevention/reversal of hepatic steatosis, hepatic inflammation and stellate cell activation. In addition, expression of inflammatory and fibrotic genes was decreased compared to the NASH-STZ group.
These results suggest that hyperglycemia is important for development of inflammation and profibrotic processes in the liver, and that insulin administration has beneficial effects on liver pathology and expression of genes related to inflammation and fibrosis in a hyperglycemic, dyslipidemic hamster model of NAFLD.
The prevalence of non-alcoholic fatty liver disease (NAFLD) as well as its more progressed form non-alcoholic steatohepatitis (NASH) is increasing world-wide  and is in the general population estimated to be approximately 25% and 1.5–6.5%, respectively . However, in patients with Type 2 diabetes mellitus (T2DM), the prevalence of NAFLD and NASH diagnosed by imaging modalities was recently estimated to be 55.5% and 37.3%, respectively . T2DM is one of the most important risk factors for developing NAFLD and NASH. It is a predictor of progressive liver disease, as well as both overall and liver-related mortality [3,4,5,6]. Many patients with T2DM will eventually need insulin treatment but whether exogenous insulin accelerates or attenuates development of or co-existing NAFLD and NASH is not known. A hypothetical adverse effect of insulin treatment on progression of NAFLD would be concerning, as NAFLD has been associated with an increased incidence of cardiovascular disease [7, 8] (although causality remains controversial [9, 10]), which accounts for the majority of deaths in T2DM .
Previous clinical trials where hepatic effects of insulin therapy in diabetic patients were explored, focused on hepatic fat content and non-specific circulating markers of the disease (e.g., aminotransferases) with no evaluation of liver histology and hence on progressive liver disease [12,13,14,15,16,17,18,19]. Insulin is known to stimulate hepatic de novo lipogenesis and might therefore be expected to increase hepatic lipid content. Nonetheless, the above-mentioned studies collectively indicate that conventional insulin treatment results in a decrease in hepatic fat content, with the largest effect observed in the initial phases of treatment [12,13,14,15,16,17,18,19], likely due to insulin-mediated inhibition of lipolysis in the adipose tissue, resulting in reduced flux of free fatty acids (FFA) to the liver . A reduction in hepatic fat content could potentially exert beneficial effects on e.g. hepatic inflammation and fibrosis, however, some retrospective studies have indicated an association between insulin therapy in diabetes and increased risk of hepatic fibrosis and hepatocellular cancer [21,22,23]. Whether these findings are causally linked to the insulin therapy per se or due to other factors is not known. The effects of insulin treatment on hepatic parameters in diabetic animal models of NAFLD/NASH have been examined in a few preclinical studies and these studies report varying effects of insulin treatment on hepatic inflammation and fibrosis [24,25,26,27]. Therefore, the effects of insulin treatment on the development and progression of NAFLD in a diabetic setting remain unclear. Furthermore, the previous studies were performed in mouse or rat models. It is well described that there are important differences between human lipid metabolism and that of mice and rats [28,29,30], which could complicate the interpretation and compromise the translational value of results from mouse and rat models in which NAFLD/NASH is induced by lipid-enriched diets.
The aim of this study was to investigate hepatic effects of insulin therapy in a NASH diet- and streptozotocin (STZ)-induced hamster model of diabetic NAFLD. The hamster was chosen as a model, because of the similarity in plasma lipid and lipoprotein metabolism compared to humans [30,31,32]. Our hypotheses were that (i) hyperglycemia accelerates liver pathology compared to effects of a high-fat diet alone and that (ii) insulin therapy improves liver pathology by preventing or reversing hepatic steatosis, inflammation and fibrosis. Using RNA-sequencing (RNAseq), we furthermore investigated global hepatic expression patterns in the model and compared it to gene expression patterns in human NAFLD, to explore the translational potential of the hamster model.
Materials and methods
Forty male Syrian hamsters were procured from Janvier Labs (Le Genest-Saint-Isle, France) at the age of six or seven weeks at arrival and weighing ~ 80–100 g. The hamsters were housed four animals/cage and were acclimatized for a week prior to initiation of the study. During the acclimatization period the animals had unrestricted access to standard rodent chow (1324 Altromin, Brogaarden, Denmark), and non-chlorinated, non-acidified tap water. Temperature in the animal rooms was set between 20 and 25 °C, with a light/dark cycle of 12/12 h, a relative humidity of 30–70% and air changes 8–15 times/h. The study was approved by the Danish Animal Experiments Inspectorate in accordance with European Union Directive 2010/63/EU.
After the acclimatization period, animals were block randomized based on body weight (restricted by cages of four) into four groups (n = 10/group), and then randomly assigned to receive either a purified control diet (D16010104, Research Diets, NJ, US, one group, Fig. 1) or a high-fat/high-fructose/high-cholesterol NASH diet (D16010102, Research Diets, NJ, US, three groups, Fig. 1). Ten animals per group were estimated to be sufficient to obtain reliable estimates of mean value and variability of the blood glucose lowering effect of human insulin. After four weeks on diets, two of the groups receiving the NASH diet were randomly assigned to be injected subcutaneously with 40 mg/kg STZ (Sigma-Aldrich, Soeborg, Denmark) suspended in 0.1 M sodium citrate buffer (pH 4.5) in the morning for three consecutive days (Fig. 1). The group receiving the control diet (CTRL) and the remaining NASH diet-fed group (NASH) were injected s.c. with the sodium citrate buffer at the same time points. One week after the STZ/sodium citrate buffer injections, blood glucose and HbA1c was measured by tail vein sampling in all animals. In addition, animals were weighed and scanned using quantitative magnetic resonance (qMR). To ensure comparable blood glucose and body weight levels, all STZ-injected animals were stratified by blood glucose and body weight into two new groups intended to receive either vehicle-treatment (NASH-STZ, n = 10) or treatment with human insulin (NASH-STZ-HI, n = 10). The hamsters were first treated with 5 mg/kg Rimadyl (Bela-Pharm GmbH, Vechta, Germany) and 5 mg/kg Baytril (KVP Pharma + Veterinär Produkte GmbH, Kiel, Germany) by s.c. injection. Then they were anaesthetized in isoflurane (Baxter, IL, USA) and a small transverse incision of approximately 1 cm was made in the subcutis between the shoulder blades. By blunt dissection, a pouch of approximately 1 × 5 cm was made in the subcutis, extending caudally from the incision site. An Alzet osmotic minipump (model: 2ML4, ALZET Osmotic Pumps, CA, US) containing either vehicle (CTRL, NASH, NASH-STZ) or human insulin (NASH-STZ-HI, dose: 40 nmol/kg/day) was then inserted in the pouch. The incision site was subsequently closed using surgical clips (Becton–Dickinson, MD, US). Following the surgical procedure, the animals were treated once daily for two days with 5 mg/kg Rimadyl (5 mg/kg). The s.c. osmotic minipumps remained implanted for four weeks (Fig. 1), during which blood glucose was monitored once a week. After four weeks of treatment a large blood sample was collected from the sublingual vein for evaluation of plasma parameters, and the animals were qMR-scanned. The animals were euthanized by exsanguination while in deep isoflurane anaesthesia. Immediately thereafter, the liver was excised and weighed, and samples were collected from the middle section of the left lateral lobe, the right medial lobe and the caudal lobe and transferred to 10% neutral buffered formalin (Hounisen Laboratorieudstyr A/S, Skanderborg, Denmark) for histopathological analysis. Another sample was collected from the right medial lobe and snap-frozen for cryo-sectioning. Six punch biopsies (diameter: 4 mm/biopsy) were then collected from the left lateral lobe and snap-frozen in liquid nitrogen for later RNA-extraction and biochemical analyses. Lastly, the epididymal fat depots (left and right) were excised and weighed.
From the STZ-injection period (week five) until termination of the study (week nine), three animals from the NASH-STZ group were euthanized due to severe hyperglycemia. Furthermore, one animal from the CTRL group was euthanized before the treatment period due to the presence of an abscess (incidental finding). Finally, one animal from the NASH-STZ-HI group was euthanized after insertion of the osmotic minipumps, due to weight loss, and failure to recover completely from the surgical procedure.
Body weight and body composition
Body weight was monitored once weekly throughout the study period. For quantification of changes in total fat mass and lean mass, animals were qMR-scanned on the day of pump insertion (treatment start) and again at study termination, using an EchoMRI Body Composition Analyser (EchoMRI, Houston, TX, USA) as described previously .
Blood and plasma samples
Blood samples for assessment of blood glucose and HbA1c were collected from the tail vein in heparin-coated capillary tubes (Vitrex Medical, Herlev, Denmark) once weekly throughout the treatment period. Blood glucose was assessed with the glucose oxidase method at a Biosen Apparatus (EKF Diagnostics, Barleben, Germany) according to manufacturer’s instructions. HbA1c was determined using a Cobas 6000 c501 instrument (Roche Diagnostics GmbH, Mannheim, Germany) according to manufacturer’s instructions. In addition, larger blood samples were collected into K3-EDTA microvette tubes (Sarstedt AG & Co, Nümbrecht, Germany) from the sublingual vein in conscious, non-fasted animals on the morning of the day of euthanization. After centrifugation, plasma was isolated and kept at − 20 °C until further analysis. Quantification of very-low-density lipoprotein triglyceride (VLDL-TG), low-density lipoprotein cholesterol (LDL-C), high-density-lipoprotein cholesterol (HDL-C) and total cholesterol in plasma was done using gel-filtration high performance liquid chromatography at LipoSEARCH (Skylight Biotech Inc, Akita, Japan). Plasma levels of free fatty acids (FFA), haptoglobin, alanine aminotransferase (ALT) and aspartate aminotransferase (AST) were measured using a Cobas 6000 c501 instrument (Roche Diagnostics GmbH, 68206 Mannheim, Germany) according to the manufacturer’s instructions. Plasma concentration of endogenous hamster insulin was measured as described previously with an assay designed for detection of rat insulin . The results are therefore presented as “rat insulin immunoactivity equivalents” (RIIE) as the cross reactivity to hamster insulin was not determined. In samples where equivalents were below the lower limit of quantification (LLOQ, 20 pM), these were set to be LLOQ/2 (10 pM) for subsequent statistical analysis, as described previously . Plasma concentration of human insulin was measured in the NASH-STZ-HI group as described previously . The LLOQ for the human insulin assay was 15 pM.
To determine hepatic levels of triglyceride (TG), total cholesterol (TC) and glycogen, liver biopsies from the left lateral lobe were homogenized in a 0.15 M sodium acetate buffer (pH = 4.9) containing 0.75% Triton-X100 (Sigma-Aldrich, Soeborg, Denmark), placed on a heating block (90–100 °C) for two min and subsequently cooled on ice. For each sample the cooled homogenate was then split into two aliquots. Amyloglucosidase (Sigma Aldrich) was added to one of the aliquots, which was then placed on a heating block (50 °C) for 2 h to facilitate the breakdown of glycogen to glucose. Both aliquots were subsequently centrifuged at 9000g for 10 min and the supernatants were analyzed using a Cobas 6000 c501 instrument (Roche Diagnostics GmbH, Mannheim, Germany) according to manufacturer’s instructions.
Sections of liver from the left lateral, right medial and caudal lobe were embedded in paraffin . Sections of 3 µm thickness were then stained with Mayer’s Haematoxylin and Eosin (H&E, Sigma-Aldrich). Cryo-sections were stained with Oil Red O (Sigma-Aldrich) to evaluate and confirm the presence of hepatic steatosis. Presence of hepatocyte ballooning was evaluated on H&E stains. Hepatic inflammation was visualized using immunohistochemistry (IHC) to detect CD68-positive (CD68+), CD11b-positive (CD11b+) and CD45-positive (CD45+) cells. Details on the different IHC protocols are shown in Additional file 1: Table S1. For a small subset of animals, a double-staining of CD68 (IHC) and Oil Red O was performed on cryo-sections from the left lateral liver lobe, to characterize the content of CD68+ cells. Collagen deposition and hepatic stellate cell activation was evaluated using histochemical staining with Picro Sirius Red (Sigma-Aldrich) and IHC for α-SMA (Additional file 1: Table S1), respectively. All sections were scanned using a NanoZoomer 2.0 HT slide scanner (Hamamatsu, Hamamatsu City, Japan) and subsequently evaluated using NanoZoomer Digital Pathology Image Software (Hamamatsu). Image analysis for quantification of inflammation (CD68, CD11b, CD45), collagen deposition (Picro Sirius Red) and hepatic stellate cell activation (α-SMA) was performed with VIS software (Visiopharm, Hoersholm, Denmark). With automated threshold analysis applications, the area with positive staining for each marker of interest was identified in each tissue section and expressed as a percentage of the total area of the tissue section (i.e., the fractional area). For each marker of interest, the same threshold for identification of positive staining was used across all tissue sections in the study.
Extraction of RNA from liver samples and RNA-sequencing
Total RNA was extracted from homogenized liver samples using a RNeasy Mini Kit (Qiagen, Aarhus, Denmark) according to the manufacturer’s instructions. The RNA integrity number (RIN) of each RNA samples was determined with a Bioanalyzer 2100 (Agilent Technologies, CA, US). RNA samples from five animals per group was sent to Qiagen (Aarhus, Denmark) for RNASeq. Average RIN-value of samples sent to Qiagen was 8.3 with SD of 0.4. At Qiagen, library preparation was done using a TruSeq® Stranded mRNA Sample preparation kit (Illumina, CA, US). A starting material of 500 ng of total RNA was mRNA-enriched using the oligodT-bead system and subsequently fragmented by enzymatic treatment. Thereafter, first- and second strand synthesis was performed, and the resulting double stranded cDNA was purified (AMPure XP, Beckman Coulter, IN, US). The cDNA was end repaired, 3′-adenylated and Illumina sequencing adaptors were ligated onto the fragment ends, after which the library was purified (AMPure XP, Beckman Coulter). The libraries size distribution was validated and the subjected to quality control using a 2200 TapeStation system (Agilent Technologies, CA, US). High quality libraries were pooled based in equimolar concentrations based on the Bioanalyzer Smear Analysis tool (Agilent Technologies). The library pool(s) were quantified using quantitative PCR and optimal concentration of the library pool was used to generate clusters on the surface of a flow cell, before sequencing on a NextSeq 500 instrument (75 cycles) according to the manufacturer’s instructions (Illumina). An average of 37 million reads were obtained for each sample. Raw reads obtained from the Illumina pipeline were mapped back to the Golden Hamster reference genome (Ensembl_93/MesAur1.0). The raw read counts were used in R 3.5.1  for data preprocessing. All non-protein coding genes as well as protein-coding genes expressed in less than 25% of the samples were filtered out. In total 11,061 of the 35,020 genes passed this filter and were analysed as described below.
Differential gene expression analysis
Differential gene expression analysis was performed using Qlucore Omics Explorer 3.0 software (Qlucore AB, Lund, Sweden) and Metacore software (Version 19.4) from Clarivate Analytics. Genes and pathways with a Benjamini-Hochberg-adjusted p-value < 0.05 (corresponding to a 5% False Discovery Rate) and fold changes below − 1.5 or above + 1.5 were considered statistically significantly regulated.
Comparison to human NASH data sets
The hamster gene expression dataset was compared to two human gene expression data sets [39, 40]. From the study by Suppli et al.  a list of differentially expressed genes between patients with histologically proven hepatic steatosis and mild inflammation (NAFL) and patients with NASH (characterized by steatosis, inflammation and fibrosis) was extracted and used to generate a heat map of these genes visualizing differences in gene expression between NASH and NASH-STZ hamsters. Furthermore, the authors of this study had defined a panel of 112 so-called “candidate NAFLD genes”, involved in e.g. insulin signaling, monocyte recruitment, inflammation signaling, hepatocellular cell death and stellate cell activation . This panel of genes was probed on the hamster data set, to explore differences in expression of the NAFLD candidate genes between the NASH and the NASH-STZ group. Data from the NAFL and NASH patient groups in the study by Lake et al.  were imported into Qlucore, and analyzed with thresholds similar to the analysis of the hamster RNAseq data (i.e., q-value < 0.05 and fold-change > 1.5-fold). The resulting list of differentially expressed genes between human NASH and NAFL patients was compared to the list of differentially expressed genes between NASH-STZ and NASH hamsters, and genes regulated in both species were plotted on heatmaps, to visualize the gene expression patterns.
Validation of RNA Sequencing results with quantitative real-time PCR
RNA extracted from homogenized liver tissue was used to synthesize cDNA using an iScript cDNA Synthesis Kit (BioRad Laboratories, Denmark) according to the manufacturer’s instructions. Eight genes were selected from the RNA sequencing analysis results for validation by quantitative real-time PCR (qRTPCR), chosen both from lists of top 10 up- and downregulated genes (Gpnmb, Pklr, Spp1, Soat1, Table 1) or from genes represented in significantly enriched pathways (Ccl2, Il1b, Col3a1, Timp1, Additional file 2: Table S2). Expression of each target gene was normalized to the average expression of the two reference genes Hprt1 and Rpl18 using the delta–delta Ct method. Hprt1 and Rpl18 both displayed stable expression across experimental groups. Primer sequences for the eight validation genes and the two reference genes are shown in Table 2 and were designed using the free online software Primer3® (Whitehead Institute for Biomedical Research, Cambridge, MA 02142, United States). The qRT-PCR was performed using the SYBR™ Green method . In brief, 8 µL of Mastermix (containing 1 µL primer mix (5 µM of forward primer and 5 µM of reverse primer), 2 µL RNAse-free water and 5 µL PowerUp™ SYBR™ Green Master Mix) and 2 µL sample cDNA was transferred to a 96-well plate. For negative controls, RNAse-free water was added to wells instead of sample. The qPCR-reaction was carried out on samples (in triplicates) in an Applied Biosystems real-time-PCR thermal cycler (StepOnePlus™ Real-Time PCR System) under the following conditions: 50 °C for 2 min, 95 °C for 5 min and 40 cycles of 95 °C for 10 s, 60 °C for 10 s and 72 °C for 20 s. Melting curves were produced and evaluated for all primer sets. The data were analyzed using StepOne Software vers. 2.3. (Applied Biosystems, CA, US).
For all data (except RNAseq data) statistical analyses and graphs were done using GraphPad Prism version 8.02 (GraphPad Software Inc., La Jolla, CA, US). Data were assessed for normal distribution and variance homology by visual inspection of qq-plots and plots of residuals versus predicted values, respectively. Data which were not normal distributed or with heterogeneous variance were transformed by the natural logarithm. If meeting the above-mentioned assumptions, differences in means between groups were analyzed using one-way ANOVA followed by Tukey’s post hoc tests for multiple pairwise comparisons between groups. Post hoc tests were performed only when the level of F reached statistical significance (p < 0.05). For data not meeting the assumptions required for the use of an ANOVA, a non-parametric Kruskal–Wallis test with subsequent Dunn test for multiple comparisons was used. p-values < 0.05 were considered statistically significant. Data are presented as mean values with SEM in figures and mean values with SD in tables.
Insulin treatment increased body weight, fat mass and lean mass in NASH-STZ hamsters
The influence of diet, hyperglycemia and insulin treatment on body composition is outlined in Fig. 2 and Table 3. The body weight (BW) of NASH-fed hamsters was higher compared to those fed CTRL diet during the diet-induction period (p < 0.05 or less, Fig. 2), but was not significantly different for the remainder of the study. From the start of the STZ-injection period until the start of the insulin treatment period (week five), the NASH-STZ and NASH-STZ-HI groups lost ~ 10% of their pre-STZ BW. During the insulin treatment period (week 5-9), the insulin-treated NASH-STZ-HI group increased in BW compared to NASH-STZ (Fig. 2) and was higher compared to the NASH-STZ group from week seven until study termination (p < 0.05 or less). Fat mass during the treatment period increased only in the insulin-treated group, where the change in fat mass was higher compared to that of the CTRL and the NASH-STZ groups at study termination (p < 0.0001, Table 3). Lean mass increased during the treatment period in all groups, and the change in lean mass in the NASH-STZ-HI group was higher compared to that of the NASH-STZ group at study termination (p < 0.05, Table 3). The epididymal fat depots collected at euthanasia weighed more in NASH and NASH-STZ-HI groups compared to the CTRL and NASH-STZ groups (p < 0.001 or less).
Pronounced hyperglycemia and dyslipidemia was seen in NASH-STZ hamsters and was partly reversed by insulin treatment
Plasma parameters for each group are shown in Table 3. Mean blood glucose was higher in STZ-injected groups compared to the CTRL- and NASH groups at treatment start (p < 0.0001, Fig. 2). Insulin treatment reduced mean blood glucose to normoglycemic levels in the NASH-STZ-HI group within the first week of treatment and remained stable and significantly lower compared to NASH-STZ for the remaining part of the study (p < 0.0001, Fig. 2). The increase in HbA1c during the study period was higher in NASH-STZ hamsters compared to CTRL and NASH hamsters (p < 0.01) and this increase was reversed by insulin treatment (p < 0.0001, Table 3). The levels of endogenous insulin (given as RIIE) were decreased in NASH-STZ hamsters compared to NASH (p < 0.05, Table 3) but was not significantly different from those of the CTRL group. Hyperglycemia resulted in pronounced dyslipidemia in the NASH-STZ group, characterized by increased plasma levels of FFA, VLDL-TG, total cholesterol and LDL-C compared to CTRL and NASH groups (p < 0.05 or less, Table 3) and decreased plasma HDL-C compared to the NASH group (p < 0.0001, Table 3). Insulin treatment significantly improved the dyslipidemia by lowering plasma levels of VLDL-TG, total cholesterol and LDL-C (p < 0.01 or less) and increasing plasma HDL-C (p < 0.0001, Table 3) compared to NASH-STZ. FFA was not significantly decreased by insulin treatment, but a trend towards 50% lower FFA levels in insulin-treated animals was observed. Assessment of ALT and AST in plasma samples from the NASH-STZ group proved difficult due to a high concentration of circulating plasma lipids. Consequently, statistical analysis was not performed on these parameters. However, the average levels of ALT and AST in plasma were numerically higher in the NASH-STZ-group compared to all other groups (Table 3). Plasma levels of 3-β-hydroxybutyrate was higher in NASH-STZ compared to both CTRL and NASH (p < 0.0001) and were decreased with insulin treatment (p < 0.0001, Table 3). Finally, levels of plasma haptoglobin were higher in NASH-STZ compared both CTRL and NASH (p < 0.0001) and were decreased with insulin treatment (p < 0.0001, Table 3).
Liver histopathology in NASH-STZ hamsters resembled NAFLD and insulin treatment had beneficial effects on several histological parameters
Representative macroscopic and histological images of the liver are shown in Figs. 3 and 4. Outputs from image analysis on IHC stains are shown in Fig. 5. Liver weights/100 g BW and hepatic biochemical parameters are shown in Table 3. H&E-stains revealed hepatic steatosis (primarily characterized by pin-point steatosis) in NASH, NASH-STZ and NASH-STZ-HI but not in CTRL (Fig. 3a–d). These findings were confirmed by Oil Red O stains (Fig. 3e–h), and biochemical parameters (hepatic triglyceride and cholesterol, Table 3). In the NASH-STZ animals, the steatosis was distributed throughout the entire liver parenchyma, whereas the steatosis in NASH and NASH-STZ-HI affected primarily zone 1 and 2 in the liver lobules. No ballooning hepatocytes were observed in any of the four groups. Image analysis showed a higher proportion of CD68+ cells in NASH-STZ livers compared to all other groups (p < 0.05 or less, Figs. 5a and 4a, b) and double-staining with CD68 and Oil Red O revealed that CD68 + cells contained lipid species (see smallest insert in Fig. 4a). Area fractions of CD45 + and CD11b+ cells were surprisingly lower in NASH-fed hamsters than CTRL hamsters (p < 0.01, Fig. 5b, c), but in NASH-STZ hamsters the CD45+ and CD11b+ area fractions were increased compared to NASH hamsters (p < 0.001). In insulin-treated animals both area fractions were significantly decreased compared to NASH-STZ hamsters (p < 0.01, Figs. 4c, e and 5b, c). There were no significant differences in hepatic collagen deposition between groups upon image analysis of Picro Sirius Red-stained sections (Figs. 4i, j and 5f), however, a higher fraction of α-SMA-positive (α-SMA+) cells (indicating activation of hepatic stellate cells ) were observed in livers of NASH-STZ animals compared to both CTRL and NASH (p < 0.0001, Fig. 5d) and was decreased after treatment with insulin (p < 0.001, Figs. 4g, h and 5d).
Hepatic expression of inflammatory and fibrogenic genes was increased in NASH-STZ hamsters and decreased by insulin treatment
Venn-diagram visualizing regulation of genes due to the effect of diet (CTRL vs NASH), hyperglycemia (NASH vs NASH-STZ) and insulin treatment (NASH-STZ vs NASH-STZ-HI) are shown in Fig. 6a. Complete lists of genes significantly regulated by diet, hyperglycemia and insulin treatment are available in Additional file 1: Tables S3–S5. Enriched pathways are shown in Fig. 6b–d and pathway-associated regulated genes are shown in Additional file 2: Table S2. Pathways analysis revealed that diet (CTRL vs. NASH) regulated pathways associated with transcriptional control of cholesterol and fatty acid synthesis, protein folding and regulation of immune responses in hepatocytes (Fig. 6b and Additional file 2: Table S2). Hyperglycemia (NASH-STZ vs NASH) impacted regulation of pathways associated with inflammatory and fibrogenic responses in the liver such as “Macrophage and dendritic cell phenotypes”, “Chemokines in inflammation in adipose tissue and liver in obesity in type 2 diabetes and metabolic syndrome X” and “Cell adhesion and ECM remodeling” (Fig. 6c & Additional file 2: Table S2). The majority of the genes in these pathways were upregulated. Insulin treatment (NASH-STZ-HI vs NASH-STZ) impacted regulation of many of the same inflammation- and fibrosis-related pathways affected by hyperglycemia, however, with the majority of pathway-associated genes being downregulated (Fig. 6d). The top 10 significantly regulated genes (up- and down) associated with each comparison (diet, hyperglycemia, insulin treatment) are shown in Table 1.
Gene expression data assessed with RNA sequencing were validated with quantitative RT-PCR
To validate the gene expression data based on RNAseq, eight genes were selected for validation using qRT-PCR (Gpnmb, Soat1, Spp1, Pklr, Ccl2, Il1β, Col3a1 and Timp1). Expression data for these eight genes assessed with qRT-PCR are shown in Table 4 and showed very good agreement between qRT-PCR and RNAseq data.
The liver transcriptome in NASH-STZ hamsters resembled the transcriptome in human NAFLD/NASH
Heat maps visualizing hamster liver expression of genes which in two previous studies were found to be differently expressed between human patients with NAFL and NASH [39, 40] are shown in Fig. 6e and Additional file 3: Figure S1, respectively. 132 genes were found to be differently expressed (q < 0.05) between human NAFL and NASH in the study by Suppli et al. . Of the 132 genes, 89 genes were identified in the hamster gene expression data set, and 37 out of these 89 genes were differently expressed by more than ± 1.5 fold when comparing NASH-STZ to NASH animals (Fig. 6e). Thirty-five of these 37 genes (95%) were similarly regulated in the NASH-STZ vs. NASH group when compared to human NASH vs. NAFL (Fig. 6e). Among others, these were genes involved in ECM remodeling, cell adhesion, immune functions, signal transduction and apoptosis (Fig. 6e). Regulation of the 37 genes in the insulin-treated animals clearly differed from that of the NASH-STZ-group, showing a greater resemblance to the transcriptomic signature of the NASH-fed group (Fig. 6e). The candidate NAFLD gene panel of 112 genes from  was probed on the hamster data set and 96 of these genes were present in the hamster data set. Of these, 38 genes were differently expressed by more than ± 1.5 fold between NASH and NASH-STZ hamsters (Fig. 6f). The 38 candidate genes were involved in monocyte recruitment, hepatocellular death, inflammation signaling, insulin signaling, lipid metabolism or stellate cell activation (Fig. 6f). In the second human data set (), 887 genes were found to be differentially expressed by more than 1.5-fold (q < 0.05) between human patients with NAFL and NASH and of those 635 genes were present in the hamster data set. With a cut-off value of minimum ± 1.5-fold difference in expression, 209 genes were identified. Heat maps of the 209 differentially expressed genes in the human comparison  and in hamsters are seen in Additional file 3: Figure S1. Here, 91% (190/209) genes showed similar expression patterns in human NASH vs NAFL and hamster NASH-STZ vs NASH group. In the human data sets 35 genes were present in both  and  and differently expressed by more than ± 1.5 fold in the comparison of NAFL vs NASH. Sixteen of these genes were also present in the hamster data set and expressed by more than ± 1.5 fold when NASH-STZ and NASH was compared (Additional file 1: Table S6). These genes were involved in functions such as ECM remodeling, immune functions, programmed cell death, signal transduction, cell–cell communication, metabolism of bile acids and post-translational modification.
In this study, we investigated the effects of STZ-induced hyperglycemia on development and progression of NAFLD in hamsters, as well as the effects of insulin treatment on liver pathology and hepatic gene expression in a dyslipidemic, hyperglycemic setting. We found that induction of hyperglycemia exacerbated pathological progression of NAFLD in the hamster when superimposed on a NASH diet. Thus, hyperglycemia per se accelerated dyslipidemia, hepatic and systemic inflammation, activated hepatic stellate cells and increased expression of pro-inflammatory and pro-fibrotic genes. As expected, insulin treatment in this model had beneficial effects on glucose control and dyslipidemia and increased body weight and fat mass, but also significantly improved liver histopathology by decreasing hepatic steatosis, inflammation and stellate cell activation. Furthermore, insulin treatment decreased the hepatic expression of pro-inflammatory and pro-fibrotic genes.
In humans, T2DM is an important risk factor for development of NAFLD and NASH and is also an independent predictor of advanced liver fibrosis and liver-related mortality [3,4,5,6]. In T2DM, hyperglycemia as well as insulin resistance and hyperinsulinemia could all likely contribute to acceleration of liver pathology, however, their individual impact on NAFLD progression in humans is not clear . In the NASH-STZ hamster we do not expect hyperinsulinemia to be a contributing factor, due to the method of diabetes induction (low-dose STZ), however, because of the prolonged hyperglycemia we do expect the NASH-STZ hamster to be insulin resistant as a consequence of glucose toxicity . Insulin resistance was likely also present in the NASH group, where the levels of endogenous insulin were almost twice as high as in the CTRL group (although with the present groups sizes this difference did not reach significance). However, the NASH group did not have progressive liver disease. This suggests that mild insulin resistance did not have a major influence on liver pathology in NASH hamsters and indicates that hyperglycemia, in the absence of hyperinsulinemia, significantly accelerated development of NAFLD in the NASH-STZ hamster.
Hyperglycemia may have accelerated liver pathology in this model by both direct and indirect mechanisms. High levels of circulating glucose could enhance development of hepatic steatosis by increasing substrate availability for hepatic de novo lipogenesis but also by directly inducing expression of carbohydrate response element-binding protein, which is a transcription factor highly capable of activating lipogenic genes in the liver . Furthermore, hypoinsulinemia in the NASH-STZ hamster could likely have increased peripheral lipolysis due to the lack of insulin-mediated lipase inhibition , which in turn results in increased flux of FFA to the liver. In fact, we observed increased levels of plasma FFAs in NASH-STZ hamsters and a significant reduction in mass of the epididymal fat depot, which may indicate increased peripheral lipolysis.
Hyperglycemia has also been shown to possess pro-inflammatory properties in clinical [46,47,48] and preclinical studies [25,26,27]. A 75 g oral glucose load as well as hyperglycemic clamping in healthy human subjects has been shown to induce marked inflammatory changes and oxidative stress with formation of reactive oxygen species (ROS) [46,47,48,49]. These effects occur also when endogenous insulin secretion is suppressed , which to some extent resembles the STZ-induced decrease in the number of functional pancreatic β-cells in our model. ROS has been shown to activate the NF-κB- pathway , one of the key inflammatory pathways in NASH . In accordance with this, the transcription factor NF-κB (denoted as “Rel” in Additional file 2: Table S2) and several of the inflammatory genes it induces (e.g. Il6, Ccl2, Ccr2, Il1b, Vcam) was upregulated in the NASH-STZ-hamster compared to the NASH group. Also, circulating levels of haptoglobin, a suggested biomarker for hepatocyte ballooning in human NASH , were increased and hepatic area fractions of CD68+ , CD11b+ and CD45+ cells were all elevated in hyperglycemic animals. The significant increase in amounts of CD68+ cells and expression levels of Cd68 in the livers of NASH-STZ hamsters indicated marked proliferation of resident macrophages (Kupffer cells) and/or recruitment of mononuclear cells. In human NASH, accumulation of CD68+ cells in the liver increase with disease stage  and macrophages (in particular Kupffer cells) play an important role in the progression towards NASH and present potential pharmacological targets [54, 55]. The STZ-NASH hamster model could therefore potentially be valuable for in evaluation of novel drugs aimed at such targets.
Fibrosis is an important prognostic marker in human NASH , but our hamster model only displayed minimal perisinusoidal fibrosis. However, the area fraction of hepatic α-SMA + cells and expression of genes involved in hepatic stellate cell activation (Acta2 (α-SMA), Col1a1, Col1a2, Col3a1, Col5a3, Col6a1, Mmp9, Timp1 and Timp2) and fibrosis (Col1a1, Col3a1, Timp1, Fap) were increased in the hyperglycemic STZ-NASH hamsters. These findings indicate that fibrotic processes are initiated but not yet fully manifested within the relatively short study period of 9 weeks. In NASH hepatic stellate cells are activated and transformed into myofibroblast-like cells, which then contribute to the fibrogenic processes [57, 58]. Transformation of hepatic stellate cells can be stimulated by high glucose concentrations in vitro [59, 60]. Furthermore, expression of connective tissue growth factor (CTGF) was significantly increased in the NASH-STZ hamster compared to NASH animals (Additional file 1: Table S4). CTGF is involved in liver fibrosis in both humans and animal models [61,62,63], and hyperglycemia has been shown to increase mRNA-levels of CTGF in vivo and in vitro . The absence of advanced fibrosis in our model was therefore likely due to the limited length of the diet- and diabetes-induction periods. More advanced stages of fibrosis may develop if these periods were prolonged. Longer diet- and diabetes induction periods may also result in ballooning of hepatocytes, which is an important hallmark of human NASH .
In this study, insulin treatment significantly decreased hepatic triglyceride levels with ~ 30%. This is in good agreement with results from several clinical studies in T2DM patients [12,13,14,15,16,17,18,19]. The effects of insulin on liver fat have been suggested to be in part due to inhibition of peripheral lipolysis, which decreases the flux of FFAs to the liver . Although insulin is paradoxically also a potent inducer of hepatic de novo lipogenesis , the relative contribution of peripheral lipolysis to liver fat accumulation in insulin-resistant patients with NAFLD has been estimated to exceed that of de novo lipogenesis . We observed a nonsignificant trend towards lower plasma FFA levels in insulin-treated animals, which could explain the decrease in hepatic fat content. The indirect insulin-induced reduction in hepatic lipid content observed in our study and in clinical trials might be theorized to reduce or prevent progression of both inflammation and fibrosis in the liver. However, retrospective studies addressing insulins effect on progression of NAFLD in diabetic patients have reported conflicting results. One study showed insulin treatment to be associated with advanced liver fibrosis , whereas another study associated insulin treatment with improvements in liver fibrosis . In addition, insulin treatment has retrospectively been associated with hepatocellular cancer [22, 23]. However, insulin treatment in the NASH-STZ-HI hamsters clearly prevented or reversed liver pathology to the level seen in NASH hamsters, by decreasing hepatic area fractions of CD68+ , CD45+ and CD11b+ cells and hepatic stellate cell activation. In addition, hepatic expression of a wide range of genes involved in inflammation and fibrosis such as Col3a1 and Timp1 was decreased. In human NASH, expression of Col3a1 has been suggested as a biomarker for discrimination between fibrotic and non-fibrotic liver, and has been shown to be positively correlated with liver inflammation activity . Timp1 is involved in regulation of extracellular matrix remodelling and has been shown to be increased in fibrotic rat and human liver [69, 70]. The top 10 genes significantly upregulated by insulin treatment, were predominantly involved in glucose- (Pppr1rb3, Gck, Pklr) and fat metabolism (Me1, Mogat2, Pnpla3). Pppr1rb3 and Gck are both involved in regulation of glycogen synthesis [71, 72], in good agreement with the increased amount of liver glycogen observed in insulin-treated animals compared to the hyperglycemic NASH-STZ hamsters. Me1, Mogat2, and Pnpla3 are involved in steps relating to fatty acid or triglyceride metabolism and may in part reflect the stimulatory effects of insulin on transcription factor SREBP1-c [73, 74]. The upregulation of these genes was generally mirrored by a downregulation in the NASH-STZ-hamster, likely reflecting restored insulin-signalling. However, in addition to these—perhaps not unexpected—effects of insulin on gene expression, the list of top 10 genes downregulated by insulin treatment revealed decreased expression of Tff3, Gpnmb, Soat1 and Spp1. Expression of these four genes were all increased in the NASH-STZ- hamster and increased expression of these genes has been implicated in progressive liver disease such as fibrosis and hepatocellular carcinoma [75,76,77,78,79,80,81,82,83], which further indicates a beneficial effect of insulin treatment on progressive liver disease.
Previous experimental studies addressing the effects of insulin treatment on hepatic inflammation, oxidative stress and fibrosis in diabetic animal models of NAFLD, have yielded varying outcomes [24,25,26,27]. In C57BL/6 mice fed a high-fat diet and subsequently injected with STZ to induce diabetes, daily s.c. injections with human insulin for 4 weeks decreased hepatic collagen deposition and expression of Timp1 but did not show effect on liver macrophage count . In contrast, treatment with insulin detemir for 6 days in STZ-induced hyperglycemic male dYY mice fed a choline-deficient high-fat diet decreased hepatic Tnfα expression but did not influence Col1a1 expression . In a study with STZ-induced hyperglycemic rats fed a chow diet, insulin was administered via osmotic minipumps for 4 weeks, similar to our study. However, glycemic control in the insulin-treated rats was rather poor, and insulin treatment did not decrease systemic and hepatic inflammation as measured by plasma α-2-macroglobuline and hepatic macrophage infiltration . In another study comparing the effects of pump delivered and injected insulin in STZ-induced hyperglycemic rats treated for 4 weeks , the glycemic control of the pump-treated group was superior to that demonstrated in . In this study hepatic histological analysis showed reduction in both oxidative stress and inflammation. Interestingly, the effects of insulin on hepatic steatosis were either only briefly mentioned or not reported in the experimental studies described above. Discrepancies in outcomes between those and the present study may be attributed to both differences in species, diet, treatment duration, method of insulin administration and/or differences in glycaemic control. In our study, it is a strength that insulin was administered by constant s.c. infusion using osmotic minipumps, mimicking basal insulin treatment in a well-regulated diabetic patient, without large fluctuations in blood glucose. The obtained glycaemic control was markedly better than in previous studies with insulin treatment in mouse and rat models [24,25,26]. Furthermore, the beneficial effects of insulin on liver histopathology observed in our study were confirmed at gene expression level, with decreased expression of genes involved in inflammation and fibrosis, which further strengthens the validity of our findings.
The way insulin treatment improved liver pathology and hepatic gene expression profile in this study is likely tightly related to the overall improvement of glycemic and metabolic status. However, insulin could potentially also exert more direct effects on inflammatory cells and hepatic gene expression. Insulin treatment has been shown to lower serum concentrations of IL-6 and high-sensitive C-reactive protein (hsCRP) in newly diagnosed T2DM patients, independently of its blood glucose lowering properties . Also, a 6-h insulin-infusion in combination with 5% dextran infusion in obese, non-diabetic individuals, showed acute suppressive actions of insulin on intranuclear NF-κB and ROS generation in mononuclear cells, with a concurrent decrease in plasma MCP-l (Ccl2) . Supporting these findings, insulin treatment decreased hepatic expression levels of NF-κB (Rel), Il-6 and Ccl2 in the NASH-STZ-HI hamster.
On a translational level, the expression patterns of genes that were regulated by more than ± 1.5 fold between human NAFL and NASH patients and were present in the hamster data set, overlapped to a large extent with the expression patterns between the NASH and NASH-STZ hamsters. Thus, even though the NASH-STZ hamsters histologically had not developed NASH (only very mild perisinusoidal fibrosis and no ballooning was observed) the transcriptome of NASH-STZ hamsters to some extent resembled the changes in gene expression observed in human NASH. Further supporting a NASH-like transcriptomic phenotype, several NAFLD candidate genes involved in both inflammatory and fibrotic signaling pathways  were significantly upregulated in NASH-STZ hamsters compared to NASH hamsters. However, there were also differences in gene expression between our hamster model and the two human datasets. This is expected when different species and different studies are compared. Furthermore, while we took care in selecting human datasets which we believe were best suited for comparison to the present data from hamsters, it should be remembered that no gold standard human NASH gene expression dataset exists. In the study by Suppli et al., it is for example surprising that expression of several genes associated with monocyte recruitment and inflammation such as Ccr1, Cd14, Cd68, Cd86, Col6a2 and Tlr4 were not upregulated in NASH vs NAFL, although these genes have been implicated in progression of NASH [86,87,88,89,90] and were also upregulated in NASH-STZ hamsters. Furthermore, the transcriptome of NAFL and NASH patients  was overlapping to a large extent, which was confirmed by histological data, with NASH patients displaying rather mild pathological changes, whereas the NAFL patients displayed not only hepatic steatosis but also mild hepatic inflammation. Also, in the study by Suppli et al., NAFL or NASH patients with confirmed diabetes were excluded. Finally, the hamster is considered a non-model organism, with the genome still only being available on a scaffold level (MesAur1.0, Ensembl Genome Browser). As such, inclusion of genes which currently still lacks annotation could potentially increase the overlap of genes present in the hamster data set and expression data from the two human data sets [39, 40].
The present study design did not allow us to determine whether the beneficial effects of insulin on liver pathology were caused by indirect effects of the restoration of glucose and lipid dysregulation or by direct effects at a cellular or molecular level. However, this would be interesting to address in future studies, possibly including other blood-glucose lowering therapies. Also, since this model is based on a relatively short-term induction- and treatment period (4 weeks on NASH-diet, 1 week for STZ-induction of hyperglycemia, followed by additional 4 weeks treatment with insulin or vehicle), it would be interesting to evaluate how the induced phenotype progress during a longer time span, and also if the treatment effects are caused by prevention or reversal of progressive liver disease. Furthermore, the individual contributions from dietary fat, fructose and cholesterol to the observed phenotype would be relevant to explore in future studies. Due to the marked hepatic steatosis accompanied by significant macrophage infiltration and hepatic stellate cell activation and because all features were achieved within a relatively short time span, the NASH-STZ hamster could be of value in investigating effects of drugs targeting these pathological aspects of NASH under hyperglycemic conditions. In addition, the enrichment and upregulation of pathways and genes involved in inflammation and fibrosis could provide insight into drug effects on transcriptional regulation related to these pathological traits.
Our results suggest that hyperglycemia is important for development of inflammation and profibrotic processes in the liver, and we show here that insulin administration has beneficial effects on liver pathology and hepatic gene expression in a hyperglycemic, dyslipidemic hamster model of NAFLD. The effect of insulin treatment could be mediated via a direct effect of insulin on inflammatory cells or through an indirect effect of normalization of blood glucose.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Free fatty acids
Haematoxylin and Eosin
Non-alcoholic fatty liver disease
Rat insulin immunoactivity equivalent
RNA integrity number
Type 2 diabetes mellitus
Quantitative magnetic resonance
Younossi ZM. Non-alcoholic fatty liver disease-a global public health perspective. J Hepatol. 2019;70:531–44.
Younossi ZM, Koenig AB, Abdelatif D, Fazel Y, Henry L, Wymer M. Global epidemiology of nonalcoholic fatty liver disease-meta-analytic assessment of prevalence, incidence, and outcomes. Hepatology. 2016;64:73–84.
Younossi ZM, Golabi P, de Avila L, Paik JM, Srishord M, Fukui N, et al. The global epidemiology of NAFLD and NASH in patients with type 2 diabetes: a systematic review and meta-analysis. J Hepatol. 2019;71:793–801.
Angulo P, Keach JC, Batts KP, Lindor KD. Independent predictors of liver fibrosis in patients with nonalcoholic steatohepatitis. Hepatology. 1999;30:1356–62.
Stepanova M, Rafiq N, Makhlouf H, Agrawal R, Kaur I, Younoszai Z, et al. Predictors of all-cause mortality and liver-related mortality in patients with non-alcoholic fatty liver disease (NAFLD). Dig Dis Sci. 2013;58:3017–23.
Nakahara T, Hyogo H, Yoneda M, Sumida Y, Eguchi Y, Fujii H, et al. Type 2 diabetes mellitus is associated with the fibrosis severity in patients with nonalcoholic fatty liver disease in a large retrospective cohort of Japanese patients. J Gastroenterol. 2014;49:1477–84.
Targher G, Byrne CD, Lonardo A, Zoppini G, Barbui C. Non-alcoholic fatty liver disease and risk of incident cardiovascular disease: a meta-analysis. J Hepatol. 2016;65:589–600.
Wu S, Wu F, Ding Y, Hou J, Bi J, Zhang Z. Association of non-alcoholic fatty liver disease with major adverse cardiovascular events: a systematic review and meta-analysis. Sci Rep. 2016;6:33386.
Stefan N, Haring HU, Cusi K. Non-alcoholic fatty liver disease: causes, diagnosis, cardiometabolic consequences, and treatment strategies. Lancet Diabetes Endocrinol. 2019;7:313–24.
Santos RD, Valenti L, Romeo S. Does nonalcoholic fatty liver disease cause cardiovascular disease? Current knowledge and gaps. Atherosclerosis. 2019;282:110–20.
de Marco R, Locatelli F, Zoppini G, Verlato G, Bonora E, Muggeo M. Cause-specific mortality in type 2 diabetes. The Verona Diabetes Study. Diabetes Care. 1999;22:756–61.
Juurinen L, Tiikkainen M, Hakkinen AM, Hakkarainen A, Yki-Jarvinen H. Effects of insulin therapy on liver fat content and hepatic insulin sensitivity in patients with type 2 diabetes. Am J Physiol Endocrinol Metab. 2007;292:E829–35.
Tang A, Rabasa-Lhoret R, Castel H, Wartelle-Bladou C, Gilbert G, Massicotte-Tisluck K, et al. Effects of insulin glargine and liraglutide therapy on liver fat as measured by magnetic resonance in patients with type 2 diabetes: a randomized trial. Diabetes Care. 2015;38:1339–46.
Bi Y, Zhang B, Xu W, Yang H, Feng W, Li C, et al. Effects of exenatide, insulin, and pioglitazone on liver fat content and body fat distributions in drug-naive subjects with type 2 diabetes. Acta Diabetol. 2014;51:865–73.
Lingvay I, Roe ED, Duong J, Leonard D, Szczepaniak LS. Effect of insulin versus triple oral therapy on the progression of hepatic steatosis in type 2 diabetes. J Investig Med. 2012;60:1059–63.
Lingvay I, Raskin P, Szczepaniak LS. Effect of insulin-metformin combination on hepatic steatosis in patients with type 2 diabetes. J Diabetes Complications. 2007;21:137–42.
Cusi K, Sanyal AJ, Zhang S, Hoogwerf BJ, Chang AM, Jacober SJ, et al. Different effects of basal insulin peglispro and insulin glargine on liver enzymes and liver fat content in patients with type 1 and type 2 diabetes. Diabetes Obes Metab. 2016;18(Suppl 2):50–8.
Liu L, Yan H, Xia M, Zhao L, Lv M, Zhao N, et al. Efficacy of exenatide and insulin glargine on nonalcoholic fatty liver disease in patients with type 2 diabetes. Diabetes Metab Res Rev. 2020;36:e3292.
Shao N, Kuang HY, Hao M, Gao XY, Lin WJ, Zou W. Benefits of exenatide on obesity and non-alcoholic fatty liver disease with elevated liver enzymes in patients with type 2 diabetes. Diabetes Metab Res Rev. 2014;30:521–9.
Frayn KN. Adipose tissue as a buffer for daily lipid flux. Diabetologia. 2002;45:1201–10.
Goh G, Pagadala M, Dasarathy J, Unalp-Arida A, Sargent R. Diabetes mellitus, insulin, sulfonylurea and advanced fibrosis in non-alcoholic fatty liver disease. J Diabetes Metab. 2014;5:2.
Donadon V, Balbi M, Ghersetti M, Grazioli S, Perciaccante A, Della Valentina G, et al. Antidiabetic therapy and increased risk of hepatocellular carcinoma in chronic liver disease. World J Gastroenterol. 2009;15:2506–11.
Zhou YY, Zhu GQ, Liu T, Zheng JN, Cheng Z, Zou TT, et al. Systematic review with network meta-analysis: antidiabetic medication and risk of hepatocellular carcinoma. Sci Rep. 2016;6:33743.
Lo L, McLennan SV, Williams PF, Bonner J, Chowdhury S, McCaughan GW, et al. Diabetes is a progression factor for hepatic fibrosis in a high fat fed mouse obesity model of non-alcoholic steatohepatitis. J Hepatol. 2011;55:435–44.
Harada S, Miyagi K, Obata T, Morimoto Y, Nakamoto K, Kim KI, et al. Influence of hyperglycemia on liver inflammatory conditions in the early phase of non-alcoholic fatty liver disease in mice. J Pharm Pharmacol. 2017;69:698–705.
Dal S, Jeandidier N, Seyfritz E, Bietiger W, Péronet C, Moreau F, et al. Featured article: oxidative stress status and liver tissue defenses in diabetic rats during intensive subcutaneous insulin therapy. Exp Biol Med. 2016;241:184–92.
Schaschkow A, Mura C, Dal S, Langlois A, Seyfritz E, Sookhareea C, et al. Impact of the type of continuous insulin administration on metabolism in a diabetic rat model. J Diabetes Res. 2016;2016:8310516.
Ye P, Cheah IK, Halliwell B. High fat diets and pathology in the guinea pig. Atherosclerosis or liver damage? Biochim Biophys Acta. 2013;1832:355–64.
Guyard-Dangremont V, Desrumaux C, Gambert P, Lallemant C, Lagrost L. Phospholipid and cholesteryl ester transfer activities in plasma from 14 vertebrate species. Relation to atherogenesis susceptibility. Comp Biochem Physiol B: Biochem Mol Biol. 1998;120:517–25.
Zhao Y, Qu H, Wang Y, Xiao W, Zhang Y, Shi D. Small rodent models of atherosclerosis. Biomed Pharmacother. 2020;129:110426.
Liu GL, Fan LM, Redinger RN. The association of hepatic apoprotein and lipid metabolism in hamsters and rats. Comp Biochem Physiol, A: Comp Physiol. 1991;99:223–8.
Arbeeny CM, Meyers DS, Bergquist KE, Gregg RE. Inhibition of fatty acid synthesis decreases very low density lipoprotein secretion in the hamster. J Lipid Res. 1992;33:843–51.
Metzinger MN, Miramontes B, Zhou P, Liu Y, Chapman S, Sun L, et al. Correlation of X-ray computed tomography with quantitative nuclear magnetic resonance methods for pre-clinical measurement of adipose and lean tissues in living mice. Sensors (Basel). 2014;14:18526–42.
Hvid H, Jensen SR, Witgen BM, Fledelius C, Damgaard J, Pyke C, et al. Diabetic phenotype in the small intestine of Zucker diabetic fatty rats. Digestion. 2016;94:199–214.
Beal SL. Ways to fit a PK model with some data below the quantification limit. J Pharmacokinet Pharmacodyn. 2001;28:481–504.
Hvid H, Blouin MJ, Birman E, Damgaard J, Poulsen F, Fels JJ, et al. Treatment with insulin analog X10 and IGF-1 increases growth of colon cancer allografts. PLoS ONE. 2013;8:e79710.
Jensen VS, Tveden-Nyborg P, Zacho-Rasmussen C, Quaade ML, Ipsen DH, Hvid H, et al. Variation in diagnostic NAFLD/NASH read-outs in paired liver samples from rodent models. J Pharmacol Toxicol Methods. 2020;101:106651.
Team RC. R: A language and environment for statistical computing.. R Foundation for Statistical Computing, Vienna, Austria. URL http://www.R-project.org/(Date of access 01/12/2014). Accessed 01/12; 2014.
Suppli MP, Rigbolt KT, Veidal SS, Heebøll S, Eriksen PL, Demant M, et al. Hepatic transcriptome signatures in patients with varying degrees of nonalcoholic fatty liver disease compared with healthy normal-weight individuals. Am J Physiol Gastrointest Liver Physiol. 2019;316:G462–72.
Lake AD, Novak P, Fisher CD, Jackson JP, Hardwick RN, Billheimer DD, et al. Analysis of global and absorption, distribution, metabolism, and elimination gene expression in the progressive stages of human nonalcoholic fatty liver disease. Drug Metab Dispos. 2011;39:1954–60.
Morrison TB, Weis JJ, Wittwer CT. Quantification of low-copy transcripts by continuous SYBR Green I monitoring during amplification. Biotechniques. 1998;24:954–8.
Carpino G, Morini S, Ginanni Corradini S, Franchitto A, Merli M, Siciliano M, et al. Alpha-SMA expression in hepatic stellate cells and quantitative analysis of hepatic fibrosis in cirrhosis and in recurrent chronic hepatitis after liver transplantation. Dig Liver Dis. 2005;37:349–56.
Williams KH, Shackel NA, Gorrell MD, McLennan SV, Twigg SM. Diabetes and nonalcoholic Fatty liver disease: a pathogenic duo. Endocr Rev. 2013;34:84–129.
Kadowaki T, Kasuga M, Akanuma Y, Ezaki O, Takaku F. Decreased autophosphorylation of the insulin receptor-kinase in streptozotocin-diabetic rats. J Biol Chem. 1984;259:14208–16.
Browning JD, Horton JD. Molecular mediators of hepatic steatosis and liver injury. J Clin Invest. 2004;114:147–52.
Mohanty P, Hamouda W, Garg R, Aljada A, Ghanim H, Dandona P. Glucose challenge stimulates reactive oxygen species (ROS) generation by leucocytes. J Clin Endocrinol Metab. 2000;85:2970–3.
Esposito K, Nappo F, Marfella R, Giugliano G, Giugliano F, Ciotola M, et al. Inflammatory cytokine concentrations are acutely increased by hyperglycemia in humans: role of oxidative stress. Circulation. 2002;106:2067–72.
Dhindsa S, Tripathy D, Mohanty P, Ghanim H, Syed T, Aljada A, et al. Differential effects of glucose and alcohol on reactive oxygen species generation and intranuclear nuclear factor-kappaB in mononuclear cells. Metabolism. 2004;53:330–4.
Dandona P, Mohanty P, Chaudhuri A, Garg R, Aljada A. Insulin infusion in acute illness. J Clin Invest. 2005;115:2069–72.
Morgan MJ, Liu ZG. Crosstalk of reactive oxygen species and NF-kappaB signaling. Cell Res. 2011;21:103–15.
Farrell GC, van Rooyen D, Gan L, Chitturi S. NASH is an inflammatory disorder: pathogenic, prognostic and therapeutic implications. Gut Liver. 2012;6:149–71.
Kamada Y, Akita M, Takeda Y, Yamada S, Fujii H, Sawai Y, et al. Serum fucosylated haptoglobin as a novel diagnostic biomarker for predicting hepatocyte ballooning and nonalcoholic steatohepatitis. PLoS ONE. 2013;8:e66328.
Gadd VL, Skoien R, Powell EE, Fagan KJ, Winterford C, Horsfall L, et al. The portal inflammatory infiltrate and ductular reaction in human nonalcoholic fatty liver disease. Hepatology. 2014;59:1393–405.
Kazankov K, Jørgensen SMD, Thomsen KL, Møller HJ, Vilstrup H, George J, et al. The role of macrophages in nonalcoholic fatty liver disease and nonalcoholic steatohepatitis. Nat Rev Gastroenterol Hepatol. 2019;16:145–59.
Schuppan D, Surabattula R, Wang XY. Determinants of fibrosis progression and regression in NASH. J Hepatol. 2018;68:238–50.
Musso G, Gambino R, Cassader M, Pagano G. Meta-analysis: natural history of non-alcoholic fatty liver disease (NAFLD) and diagnostic accuracy of non-invasive tests for liver disease severity. Ann Med. 2011;43:617–49.
Friedman SL. Molecular regulation of hepatic fibrosis, an integrated cellular response to tissue injury. J Biol Chem. 2000;275:2247–50.
Schuppan D, Ruehl M, Somasundaram R, Hahn EG. Matrix as a modulator of hepatic fibrogenesis. Semin Liver Dis. 2001;21:351–72.
Sugimoto R, Enjoji M, Kohjima M, Tsuruta S, Fukushima M, Iwao M, et al. High glucose stimulates hepatic stellate cells to proliferate and to produce collagen through free radical production and activation of mitogen-activated protein kinase. Liver Int. 2005;25:1018–26.
Lin J, Chen A. Curcumin diminishes the impacts of hyperglycemia on the activation of hepatic stellate cells by suppressing membrane translocation and gene expression of glucose transporter-2. Mol Cell Endocrinol. 2011;333:160–71.
Paradis V, Dargere D, Vidaud M, De Gouville AC, Huet S, Martinez V, et al. Expression of connective tissue growth factor in experimental rat and human liver fibrosis. Hepatology. 1999;30:968–76.
Williams EJ, Gaca MD, Brigstock DR, Arthur MJ, Benyon RC. Increased expression of connective tissue growth factor in fibrotic human liver and in activated hepatic stellate cells. J Hepatol. 2000;32:754–61.
Paradis V, Perlemuter G, Bonvoust F, Dargere D, Parfait B, Vidaud M, et al. High glucose and hyperinsulinemia stimulate connective tissue growth factor expression: a potential mechanism involved in progression to fibrosis in nonalcoholic steatohepatitis. Hepatology. 2001;34:738–44.
Yeh MM, Brunt EM. Pathology of nonalcoholic fatty liver disease. Am J Clin Pathol. 2007;128:837–47.
Ferre P, Foufelle F. Hepatic steatosis: a role for de novo lipogenesis and the transcription factor SREBP-1c. Diabetes Obes Metab. 2010;12(Suppl 2):83–92.
Donnelly KL, Smith CI, Schwarzenberg SJ, Jessurun J, Boldt MD, Parks EJ. Sources of fatty acids stored in liver and secreted via lipoproteins in patients with nonalcoholic fatty liver disease. J Clin Invest. 2005;115:1343–51.
Hamaguchi E, Takamura T, Sakurai M, Mizukoshi E, Zen Y, Takeshita Y, et al. Histological course of nonalcoholic fatty liver disease in Japanese patients: tight glycemic control, rather than weight reduction, ameliorates liver fibrosis. Diabetes Care. 2010;33:284–6.
Staten NR, Welsh EA, Sidik K, McDonald SA, Dufield DR, Maqsodi B, et al. Multiplex transcriptional analysis of paraffin-embedded liver needle biopsy from patients with liver fibrosis. Fibrogenesis Tissue Repair. 2012;5:21.
Herbst H, Wege T, Milani S, Pellegrini G, Orzechowski HD, Bechstein WO, et al. Tissue inhibitor of metalloproteinase-1 and -2 RNA expression in rat and human liver fibrosis. Am J Pathol. 1997;150:1647–59.
Benyon RC, Iredale JP, Goddard S, Winwood PJ, Arthur MJ. Expression of tissue inhibitor of metalloproteinases 1 and 2 is increased in fibrotic human liver. Gastroenterology. 1996;110:821–31.
Mehta MB, Shewale SV, Sequeira RN, Millar JS, Hand NJ, Rader DJ. Hepatic protein phosphatase 1 regulatory subunit 3B (Ppp1r3b) promotes hepatic glycogen synthesis and thereby regulates fasting energy homeostasis. J Biol Chem. 2017;292:10444–54.
Petersen MC, Vatner DF, Shulman GI. Regulation of hepatic glucose metabolism in health and disease. Nat Rev Endocrinol. 2017;13:572–87.
Huang Y, He S, Li JZ, Seo YK, Osborne TF, Cohen JC, et al. A feed-forward loop amplifies nutritional regulation of PNPLA3. Proc Natl Acad Sci USA. 2010;107:7892–7.
Shimano H, Yahagi N, Amemiya-Kudo M, Hasty AH, Osuga J, Tamura Y, et al. Sterol regulatory element-binding protein-1 as a key transcription factor for nutritional induction of lipogenic enzyme genes. J Biol Chem. 1999;274:35832–9.
Xiong X, Kuang H, Ansari S, Liu T, Gong J, Wang S, et al. Landscape of intercellular crosstalk in healthy and NASH liver revealed by single-cell secretome gene analysis. Mol Cell. 2019;75(644–60):e5.
Jiang Y, Sun A, Zhao Y, Ying W, Sun H, Yang X, et al. Proteomics identifies new therapeutic targets of early-stage hepatocellular carcinoma. Nature. 2019;567:257–61.
Khoury T, Chadha K, Javle M, Donohue K, Levea C, Iyer R, et al. Expression of intestinal trefoil factor (TFF-3) in hepatocellular carcinoma. Int J Gastrointest Cancer. 2005;35:171–7.
Syn WK, Agboola KM, Swiderska M, Michelotti GA, Liaskou E, Pang H, et al. NKT-associated hedgehog and osteopontin drive fibrogenesis in non-alcoholic fatty liver disease. Gut. 2012;61:1323–9.
Lorena D, Darby IA, Gadeau AP, Leen LL, Rittling S, Porto LC, et al. Osteopontin expression in normal and fibrotic liver. altered liver healing in osteopontin-deficient mice. J Hepatol. 2006;44:383–90.
Gerhard GS, Legendre C, Still CD, Chu X, Petrick A, DiStefano JK. Transcriptomic profiling of obesity-related nonalcoholic steatohepatitis reveals a core set of fibrosis-specific genes. J Endocr Soc. 2018;2:710–26.
Graveel CR, Jatkoe T, Madore SJ, Holt AL, Farnham PJ. Expression profiling and identification of novel genes in hepatocellular carcinomas. Oncogene. 2001;20:2704–12.
Okada H, Kimura MT, Tan D, Fujiwara K, Igarashi J, Makuuchi M, et al. Frequent trefoil factor 3 (TFF3) overexpression and promoter hypomethylation in mouse and human hepatocellular carcinomas. Int J Oncol. 2005;26:369–77.
Katayama A, Nakatsuka A, Eguchi J, Murakami K, Teshigawara S, Kanzaki M, et al. Beneficial impact of Gpnmb and its significance as a biomarker in nonalcoholic steatohepatitis. Sci Rep. 2015;5:16920.
Mao XM, Liu H, Tao XJ, Yin GP, Li Q, Wang SK. Independent anti-inflammatory effect of insulin in newly diagnosed type 2 diabetes. Diabetes Metab Res Rev. 2009;25:435–41.
Dandona P, Aljada A, Mohanty P, Ghanim H, Hamouda W, Assian E, et al. Insulin inhibits intranuclear nuclear factor kappaB and stimulates IkappaB in mononuclear cells in obese subjects: evidence for an anti-inflammatory effect? J Clin Endocrinol Metab. 2001;86:3257–65.
Seki E, Schwabe RF. Hepatic inflammation and fibrosis: functional links and key pathways. Hepatology. 2015;61:1066–79.
Imajo K, Fujita K, Yoneda M, Nozaki Y, Ogawa Y, Shinohara Y, et al. Hyperresponsivity to low-dose endotoxin during progression to nonalcoholic steatohepatitis is regulated by leptin-mediated signaling. Cell Metab. 2012;16:44–54.
Miura K, Yang L, van Rooijen N, Ohnishi H, Seki E. Hepatic recruitment of macrophages promotes nonalcoholic steatohepatitis through CCR2. Am J Physiol Gastrointest Liver Physiol. 2012;302:G1310–21.
Veyel D, Wenger K, Broermann A, Bretschneider T, Luippold AH, Krawczyk B, et al. Biomarker discovery for chronic liver diseases by multi-omics-a preclinical case study. Sci Rep. 2020;10:1314.
Bertola A, Bonnafous S, Anty R, Patouraux S, Saint-Paul MC, Iannelli A, et al. Hepatic expression patterns of inflammatory and immune response genes associated with obesity and NASH in morbidly obese patients. PLoS ONE. 2010;5:e13577.
The authors would like to thank Susanne Juul Rasmussen, Jeanette Bannebjerg Johansen and Karsten Marckstrøm for valuable technical assistance with preparation, handling and staining of histological sections.
The work presented here was supported by Novo Nordisk A/S and the LifePharm Centre of In Vivo Pharmacology at University of Copenhagen.
Ethics approval and consent to participate
The study was approved by the Danish Animal Experiments Inspectorate in accordance with European Union Directive 2010/63/EU.
Consent for publication
The authors disclose the following competing interests: VSJ, CF, CZ, JD, HN, KST, BMV, RKK, EMW and HH are present or former employees of Novo Nordisk A/S and hold shares in the company. VSJ and JL are employed at University of Copenhagen and have collaborated with Novo Nordisk on this project.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Table S1.
Immunohistochemistry protocols. Overview of the immunohistochemistry protocols used in the study. Table S3. Gene expression in NASH vs CTRL hamsters. Fold changes and log2-fold changes in gene expression between NASH and CTRL hamsters. Table S4. Gene expression in NASH-STZ vs NASH hamsters. Fold changes and log2-fold changes in gene expression between NASH-STZ and NASH hamsters. Table S5. Gene expression in NASH-STZ-HI vs NASH-STZ hamsters. Fold changes and log2-fold changes in gene expression between NASH-STZ-HI and NASH-STZ hamsters. Table S6. Overview of genes significantly regulated in all three data sets (two human data sets and one hamster data set). List of genes exhibiting similar (and significant) gene regulation in both the human data sets used as well as the hamster data set.
Additional file 2. Table S2.
Regulation of genes in enriched pathways. Overview of regulated genes in enriched pathways.
Additional file 3: Figure S1.
Regulation of genes in NASH-STZ hamsters and human NASH patients. Heat map of 209 genes differentially regulated between human NAFL (simple steatosis) and NASH patients (Lake et al., 2011, , right panel) compared to the same genes in NASH-fed hamsters and NASH-STZ hamsters (left panel). Of the 209 genes, 190 genes showed similar expression patterns in the two comparisons.
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 http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.
About this article
Cite this article
Jensen, V.S., Fledelius, C., Zachodnik, C. et al. Insulin treatment improves liver histopathology and decreases expression of inflammatory and fibrogenic genes in a hyperglycemic, dyslipidemic hamster model of NAFLD. J Transl Med 19, 80 (2021). https://doi.org/10.1186/s12967-021-02729-1
- Non-alcoholic fatty liver disease
- Non-alcoholic steatohepatitis
- Animal models
- Gene expression
- Insulin treatment
- Insulin therapy