- Open Access
Experimental non-alcoholic steatohepatitis in Göttingen Minipigs: consequences of high fat-fructose-cholesterol diet and diabetes
Journal of Translational Medicine volume 17, Article number: 110 (2019)
Non-alcoholic fatty liver disease (NAFLD) is the most common liver disease in humans, and ranges from steatosis to non-alcoholic steatohepatitis (NASH), the latter with risk of progression to cirrhosis. The Göttingen Minipig has been used in studies of obesity and diabetes, but liver changes have not been described. The aim of this study was to characterize hepatic changes in Göttingen Minipigs with or without diabetes, fed a diet high in fat, fructose, and cholesterol to see if liver alterations resemble features of human NAFLD/NASH.
Fifty-four male castrated minipigs (age 6 to 7 months) were distributed into four groups and diet-fed for 13 months. Groups were: lean controls fed standard diet (SD, n = 8), a group fed high fat/fructose/cholesterol diet (FFC, n = 16), a group fed high fat/fructose/cholesterol diet but changed to standard diet after 7 months (diet normalization, FFC/SD, n = 16), and a streptozotocin-induced diabetic group fed high fat/fructose/cholesterol diet (FFCDIA, n = 14). At termination, blood samples for analyses of circulating biomarkers and liver tissue for histopathological assessment and analyses of lipids and glycogen content were collected.
In comparison with SD and FFC/SD, FFC and FFCDIA pigs developed hepatomegaly with increased content of cholesterol, whereas no difference in triglyceride content was found. FFC and FFCDIA groups had increased values of circulating total cholesterol and triglycerides and the hepatic circulating markers alkaline phosphatase and glutamate dehydrogenase. In the histopathological evaluation, fibrosis (mainly located periportally) and inflammation along with cytoplasmic alterations (characterized by hepatocytes with pale, granulated cytoplasm) were found in FFC and FFCDIA groups compared to SD and FFC/SD. Interestingly, FFC/SD also had fibrosis, a feature not seen in SD. Only two FFC and three FFCDIA pigs had > 5% steatosis, and no hepatocellular ballooning or Mallory–Denk bodies were found in any of the pigs.
Fibrosis, inflammation and cytoplasmic alterations were characteristic features in the livers of FCC and FFCDIA pigs. Overall, diabetes did not exacerbate the hepatic changes compared to FFC. The limited presence of the key human-relevant pathological hepatic findings of steatosis and hepatocellular ballooning and the variation in the model, limits its use in preclinical research without further optimisation.
It has recently been estimated that human non-alcoholic fatty liver disease (NAFLD) has a 25% worldwide prevalence and is expected to become the major reason for liver transplantations in the western world [1, 2].
NAFLD ranges from simple steatosis to non-alcoholic steatohepatitis (NASH), the latter characterized by additional inflammation, hepatocellular hydropic degeneration (also called ballooning) and eventually Mallory–Denk bodies associated with fibrosis which can lead to cirrhosis. Type 2 diabetes is considered the primary risk factor for progression of simple steatosis to advanced stages of NASH . In addition, NAFLD seems to be a central risk factor for the development of cardiovascular disease  and other complications also related to diabetes [5, 6]. It can be difficult to diagnose NAFLD given that clinical signs usually are sparse even at late stages where liver damage is substantial; and changes in circulating liver biomarkers are often non-specific. Histopathological evaluation of samples from liver biopsies is still considered the gold standard for diagnosing the disease even though new non-invasive alternatives are emerging .
The hepatic disease mechanisms are still largely unknown, and effective treatment modalities are lacking. Therefore, various animal models have been investigated to help unravel the pathogenesis and for testing new pharmaceutical drug candidates or lifestyle interventions. Diet- or chemically-induced rodent models of NASH/NAFLD have mainly been used, and several mouse strains exist with spontaneous or transgenic mutations that pinpoint different signalling pathways [8, 9]. Larger animal models have also been examined, as anatomy, physiology or simply size can be beneficial depending on the research question of interest [10, 11].
Histopathological features resembling NASH in humans have been reported in Ossabaw minipigs [10, 12] fed a diet high in fat, fructose and cholesterol. Others have studied Bama minipigs and Microminipigs on high fat/high sucrose diet or high fat/high cholesterol diet with a supplement of cholic acid, respectively, and reported microvesicular steatosis and inflammation, but none or limited hepatocellular ballooning and fibrosis [11, 13]. The Göttingen Minipig is widely used in studies of obesity and metabolic syndrome [14,15,16], but diet-induced histopathological changes in the liver have not been characterized.
In this study, the aim was to investigate the liver changes in Göttingen Minipigs fed a diet with high content of fat, fructose and cholesterol for over a year. Hepatic gross morphology, histopathology and tissue content of lipids and glycogen as well as changes in relevant circulating biomarkers were investigated. Furthermore, it was evaluated if diabetes would exacerbate the alterations and if a diet change to standard diet would limit pathologic changes. The hypothesis was that Göttingen Minipigs fed a diet high in fat, fructose and cholesterol would develop hepatic changes that resemble human NAFLD including NASH.
Materials and methods
Study setup and animals
Castrated male Göttingen Minipigs (Ellegaard Göttingen Minipigs A/S, Dalmose, Denmark) (n = 54 in total) aged 6 to 7 months were weight stratified and distributed into four groups (Fig. 1a) and fed once daily for 13 months. The pigs were part of a larger study (n = 84 in total) looking at different aspects of obesity  and diabetes related complications (two intervention groups were not included in the present study). The included groups were: lean control pigs (SD, n = 8) fed standard diet (Mini-pig, SDS, UK (Table 1)); a group fed high fat/fructose/cholesterol (2%) diet (5B4L) for the first 5 months and subsequently changed to a similar diet with 1% cholesterol (9G4U) for the next 8 months [both diets from Test diet®, Missouri, USA (Table 1)] (high fat/fructose/cholesterol, FFC, n = 16); a group fed the same diets as FFC for the first 7 months but returned to standard diet the last 6 months of the study (diet-normalization, FFC/SD, n = 16) and a diabetic group fed a high fat/fructose/cholesterol (1%) diet throughout the study (9G4U) (FFCDIA, n = 14). Details of diet feeding in the four groups can be seen in Fig. 1b. Straw was used as bedding material, and the animals had access to fresh drinking water at all times.
The diabetic group had type 1-like diabetes induced with streptozotocin administered intravenously once daily for three consecutive days  (modified version of Gerrity et al. optimized for Göttingen Minipigs; 60 mg/kg bodyweight). Diabetic animals were treated once daily with a long acting insulin (insulin glargine, Lantus®, Sanofi-Aventis Deutschland GmbH, Frankfurt am Main, Germany) subcutaneously immediately after the daily diet ration was offered in the morning; targeting fasting morning glucose levels of 14–16 mM.
Three to eight weeks before termination, blood for quantification of circulating biomarkers and intravenous glucose tolerance test (IVGTT) were sampled in plain serum tubes and EDTA coated tubes. Blood samples were kept on ice for maximum 30 min (EDTA) or kept at room temperature (serum) for 1 h before centrifugation for 10 min, 2000×g at 4 °C. After centrifugation EDTA plasma and serum were pipetted into appropriate tubes and stored at − 80 °C. The IVGTT procedure including calculation of intravenous glucose tolerance index (KG) and area under the curve for insulin response (AUCInsulin) were performed as previously described . Homeostasis model assessment of insulin resistance (HOMA-IR) was calculated as fasting insulin (µU/mL) × fasting glucose (mM)/22.5  in all animals except for the diabetic animals since their fasting glucose values were also influenced by the long-acting insulin glargine. Total body fat percentage (BF%) was estimated using dual energy x-ray absorptiometry (DEXA) (Lunar prodigy, GE Healthcare, Brøndby, Denmark) 1 to 6 weeks before termination. All animals were euthanized by exsanguination in general anaesthesia (mixture of zolazepam, tiletamine, ketamine, xylazine and butorphanol as previously described by Pedersen et al. ). After termination, the liver was removed and weighed. A sample (1 × 1 × 1–2 cm) from each of the four main liver lobes [Lobus hepatis sinister medialis (SM) and lateralis (SL), Lobus hepatis dexter medialis (DM) and lateralis (DL)] was collected mid-lobe and fixed in 10% neutral buffered formalin for 24–32 h, processed, and embedded in paraffin. From one liver lobe (SM) samples were obtained, snap frozen in liquid nitrogen and kept at − 80 °C or cryo fixated in optimum cutting temperature compound (OCT) (Tissue-Tek® O.C.T. Compound, Sakura® Finetek, Alphen aan den Rijn, The Netherlands) on dry ice. Animals were fasted overnight before blood sampling, DEXA assessment, and termination.
Analysis of circulating biomarkers
Plasma concentrations of alkaline phosphatase (ALP), alanine transaminase (ALT), aspartate transaminase (AST), glutamate dehydrogenase (GLDH), total cholesterol (TC), triglycerides (TG), glucose (GLU), fructosamine (FRA), and albumin (ALB) were quantified using an autoanalyzer Cobas 6000® (Roche A/S, Hvidovre, Denmark). A modified version of a previously described method  was used for quantification of C-reactive protein in serum (CRP).
Biochemical analysis of liver tissue content of triglycerides, cholesterol and glycogen
Snap frozen SM liver samples were homogenized as previously described  and the content of cholesterol, triglycerides and glycogen were measured by Cobas 6000® in duplicates. Six samples were analyzed six times for the determination of intra-assay variation.
The liver samples were fixed in neutral buffered formalin for 24–32 h. From each lobe one six mm tissue core was punched. The cores were embedded in paraffin in a recognisable pattern [tissue micro array (TMA)]. Three µm thick sections were cut and routinely stained with hematoxylin and eosin (HE), picro-sirius red (PSR), and periodic acid-Schiff (PAS) with and without diastase pre-treatment. Selected elements from the NAFLD activity score (NAS)  were used to evaluate steatosis (HE), fibrosis (PSR) and inflammation (HE). The latter was evaluated by the number of inflammatory foci (defined as at least 5 extravascular inflammatory cells in a cluster). Cytoplasmic alterations (CA) in the hepatocytes (HE) and the glycogen content (PAS) were scored from 0 (none) to 3 (severe). For all findings, the localization or zonal distribution was recorded. In Table 2, an overview of the different elements in the qualitative histopathological assessment is provided. All assessments were performed in a blinded manner by one observer and, moreover, eight samples were scored three times for the determination of intra-observer variation.
Immunohistochemical (IHC) staining using a polyclonal goat anti-ionized calcium binding adapter molecule 1 (Iba1) antibody (Abcam #ab5076) was applied to detect macrophages . Antigen retrieval on deparaffinized slides was done by microwave heating in TEG-buffer pH 9.0 (Ampliqon A/S, Odense, Denmark) for 15 min before endogenous peroxidase was inhibited with hydrogen peroxide (3%) (Merck, Darmstadt, Germany) for 10 min. Both avidin and biotin (Dako Biotin blocking system, Dako A/S, Glostrup, Denmark) was used for 10 min each for blockage of endogenous biotin. Pre-incubation with donkey serum (7%) (Jackson Immunoresearch, West Grove PA, USA) in Tris-buffered saline (TBS) (Ampliqon A/S, Odense, Denmark) + Tween 20 (TBS + T) (Merck, Darmstadt, Germany) and skimmed milk (5%) (BD, Kgs. Lyngby, Denmark) for blockage of unspecific antibody binding was performed for 30 min. Primary goat antibody against Iba1 (1:1500) (Abcam, Cambridge, UK, #5076) was applied for 120 min; thereafter secondary antibody (biotinylated donkey-anti-goat 1:1000) (Jackson Immunoresearch, West Grove PA, USA, #705-065-147) was applied for 30 min. Both antibodies were dissolved in a TBS + T + 0.5% skimmed milk solution. ABC-detection system (VECTASTAIN® Elite® ABC-HRP Kit, Vector Laboratories, Burlingame, CA, USA) was used for 30 min followed by chromogen (DAB + chromogen, Dako A/S, Glostrup, Denmark) for 10 min before counterstaining was performed with Mayer’s hematoxylin (Sigma-Aldrich/Merck, Darmstadt, Germany). Between every step, except blockage of unspecific antibody binding, slides were washed with TBS + T. A negative control slide was used without primary antibody, and two different non-conjugated goat IgG were used for check of unspecific staining [‘ChromPure Goat IgG’ (1:10,000), Jackson Immunoresearch, West Grove PA, USA, #005-000-003 and ‘Normal Goat IgG Control’ (1:1000), R&D Systems, Oxon, UK, #AB-108-C].
Also, 10 µm sections from OCT embedded frozen blocks from lobe SM were cut and stained with Oil red O (ORO) in order to visualize the distribution of lipid and confirm steatosis found on HE sections. Image analysis software (Visiopharm, Hørsholm, Denmark) was used to quantify fibrosis (PSR-VIS), inflammation (Iba1-VIS) and lipid (ORO-VIS) by measuring the stained areas. For all TMAs (PSR-VIS and Iba1-VIS), the measured area was adjusted for the number of lobules in order to normalise for the difference in lobular diameter. The lobular diameter was estimated as the number of all lobules represented (both whole and fractioned) per TMA, since all areas of TMAs were nearly the same; a low number indicating a larger lobular diameter. Total tissue area for frozen OCT samples (ORO-VIS) varied in size which led to correction of the stained area with total tissue area.
Shapiro–Wilk test was used to test whether data were normal distributed in groups and due to the majority of the evaluated parameters not being normally distributed, all data were expressed in tables as medians and interquartiles.
Differences among the four liver lobes in histopathology were investigated by a linear model taking repeated measurements into account.
Group differences in basic characteristics, circulating biomarkers, liver tissue content, and histopathological parametric data were tested using analysis of variance (ANOVA) with Welch adjustment for unequal variance between groups and Tukey’s correction for multiple comparisons. When needed, an appropriate transformation of the response variables was performed to obtain variance homogeneity and normal distribution of residuals in the ANOVA analyses. Fisher’s exact test was used to test categorical data for group differences.
In addition, linear regression was used to test for associations between histopathological data and circulating biomarkers or tissue content. Associations between histopathological scores were evaluated using Fisher’s exact test.
Statistical analyses were performed using SAS 9.4 software (SAS Institute Inc., Cary, NC, USA) and the significance level was set to p < 0.05. Graphs were created with GraphPad Prism (GraphPad Software Inc., La Jolla, CA, USA).
Forty-three pigs completed the study. Eleven pigs were excluded (Fig. 1a). Five pigs were euthanized according to predefined endpoints and six died suddenly: one SD (euthanized due to procedure related complications), two FFC (one euthanized due to procedure related complications and another died suddenly with a tentative diagnosis of circulatory failure with accumulation of serosanguinous fluid or blood in body cavities), two FFC/SD [one due to lameness, another due to spontaneous acute bleeding (diagnosed thrombocytopenia purpura)] and six FFCDIA (two due to severe necrosis/inflammation of pancreas), where one died suddenly, and four died suddenly due to pulmonary embolus (n = 3), or a tentative diagnosis of circulatory failure with accumulation of serosanguinous fluid or blood in body cavities (n = 1). Diagnoses were based on clinical signs and gross pathology.
Higher body weight (BW), liver weight (LW) and total body fat percentages (BF%) were found in FFC and FFCDIA groups compared to the SD group, whereas only BW and BF% were higher in the FFC/SD group (Table 3). By macroscopic inspection, the majority of livers from FFC and FFCDIA pigs were enlarged, pale red or yellow in color and with a bulgy appearance compared to SD. Livers from FFC/SD animals did not differ markedly from SD in size, color and gross appearance (Fig. 2).
Circulating hepatic biomarkers ALP and GLDH were increased in FFC and FFCDIA groups compared to SD, whereas the FFC group had lower values of ALT (Table 3). No statistically significant difference was found for AST or the inflammation markers CRP and ALB. FFC and FFCDIA groups had higher TG and TC than SD and FFC/SD groups with TC further elevated in FFCDIA compared to FFC. As expected, fasting GLU and FRA values were significantly elevated in the FFCDIA group compared to the three other groups, and in addition the FFCDIA group had significantly lower KG and AUCInsulin than the other groups (Table 3). Furthermore, KG in FFC was significantly decreased compared to SD and FFC/SD, suggesting impaired glucose metabolism in the two groups FFC and FFCDIA. The lower KG despite a numerically ~ twofold higher AUCinsulin in the FFC group compared to SD indicates that the FFC diet and following obesity induces insulin resistance, which is also supported by the significantly higher HOMA-IR in this group compared to SD. The additional lowering of KG in the FFCDIA group compared to the FFC group is due to the insulin depletion conferred by the STZ-procedure. The only difference between FFC/SD and SD was a significantly higher HOMA-IR in the FFC/SD, indicating that with this index insulin resistance had not completely normalised yet.
Biochemical analyses of liver tissue content
Liver tissue content of triglyceride content did not differ between the groups (Fig. 3a), but cholesterol was higher in FFC and FFCDIA groups compared to SD and FFC/SD groups (Fig. 3b). Regarding glycogen, FFCDIA had lower liver content than FFC (Fig. 3c).
Only few and minor differences in histopathological changes were observed among the four liver lobes indicating that the alterations were evenly distributed in the liver. Only slightly higher steatosis scores were found in SL (4 out of 43 pigs had score ≥ 1) and SM (5/43) compared to DL and DM (both 0/43) (p < 0.0001), and regarding inflammatory foci, DM (20/43) scored higher than DL (13/43) (p = 0.03). Further detailed evaluation of the liver lobe SM showed that steatosis, predominantly of macrovesicular type, was found in pigs from FFC and FFCDIA groups with a centrilobular distribution, but only five animals (FFC = 2, FFCDIA = 3) had ≥ 5% parenchymal involvement and no significant difference was found among groups (Fig. 4).
An increased content of collagen was found in FFC, FFCDIA and FFC/SD groups in comparison to SD. The majority of these pigs had periportal fibrosis; however, four pigs (all FFC) had both periportal and randomly distributed lobular fibrosis. In addition, five pigs (FFC = 2, FFC/SD = 1, FFCDIA = 2) showed porto-central bridging fibrosis (Fig. 4). None of the cases showed centrilobular pericellular fibrosis.
The number of lobular inflammatory foci was increased in FFC and FFCDIA groups compared to both SD and FFC/SD (Fig. 4).
Hepatocellular cytoplasmic alterations, characterized by hepatocytes with pale granular cytoplasm (Fig. 5), were a characteristic feature in most FFC and FFCDIA pigs, in contrast to both SD and FFC/SD, where no alterations were seen (Fig. 4). These changes were homogenous and widely distributed and appear to be similar to those seen in Ossabaw minipigs using the same diet [10, 12]. The changes have some resemblance to the hepatocellular ballooning in human NASH such as rounded hepatocyte enlargement and presence of pale cytoplasm but lack key characteristics such as cytoplasmatic vacuoles and Mallory–Denk bodies . Further, no difference among groups was found for glycogen staining with PAS with or without pre-treatment with diastasis (Fig. 4); in contrast to the biochemical analysis that showed FFCDIA had less liver tissue glycogen content than FFC.
For all TMAs, total tissue area was nearly the same whereas the number of lobules within that area varied between pigs from seven to fifty-four (lower number = larger lobular diameter). The FFCDIA group had larger lobular diameter than the three other groups (all p < 0.01), while the lobular diameter in the FFC group was also larger than in the FFC/SD group (p = 0.003) (Fig. 6). Quantification of collagen content using PSR revealed excessive collagen for FFC, FFCDIA and FFC/SD compared to SD (Fig. 6). When staining with ORO, larger areas of lipids were found for FFC and FFCDIA groups compared to both SD and FFC/SD groups (Fig. 6); however, a large degree of heterogeneity was seen within the groups (Fig. 7). IHC staining for Iba1 indicated increased number of macrophages in both FFC and FFCDIA groups as compared to SD and FFC/SD groups (Fig. 6). Figure 8 shows typical examples of histopathological findings for each group.
Associations between circulating biomarkers and histopathology
Regarding associations between circulating biomarkers and parametric histopathological variables (Table 4), lipid parameters TC and TG were associated with both ORO-VIS and Iba1-VIS, whereas GLDH was associated with ORO-VIS and ALP to Iba1-VIS.
Comparison between circulating biomarkers and categorical histopathological features (Table 4), revealed association between both circulating lipids TG and TC and the inflammation score and especially the CA score. For circulating hepatic markers, associations were seen between most markers and the CA score, fibrosis score, and inflammation score.
Associations between histopathological findings
For the histopathological features applied, the cytoplasmic alteration score was associated with all other features except the PAS-score (Table 5). In fact, the PAS-score showed no association to any of the other histopathological findings. Scores for both fibrosis and inflammatory foci also had a strong association to the quantification of fibrosis (PSR-VIS) (Table 5).
Associations between biochemical analyses of liver tissue content and histopathological findings
As expected, tissue triglyceride content showed a strong association to the quantification of lipid (ORO-VIS) and glycogen content was strongly associated to the PAS score (Table 6), even though the PAS ± diastasis staining technique can be challenging in terms of accuracy. In addition, a significant association was found between liver tissue content of cholesterol and most histopathological findings, especially CA score.
Intra-assay and intra-observer variations
Lastly the intra-assay coefficients of variance (CV), for all three biochemical tissue content analyses were below 5%. For the histopathological assessment of fibrosis, inflammatory foci, cytoplasmic alterations and PAS positive cells the assigned score differed no more than one when performed by the same observer. No intra-observer difference was found for steatosis score.
This study demonstrated that Göttingen Minipigs fed a high fat diet with fructose and cholesterol became obese and developed hepatomegaly with hepatic fibrosis, inflammation and cytoplasmic alterations when compared to animals fed a normal diet. However, only few animals developed marked steatosis, and hepatocellular hydropic degeneration (ballooning) and Mallory–Denk bodies were not observed in any animal. Streptozotocin-induced diabetes did not exacerbate the changes in circulating biomarkers or hepatic histopathology, compared to non-diabetic animals fed a similar diet. The group changed to standard diet for 6 months had no hepatic changes, except from excess collagen deposition.
The large number of animals, the more than 1-year study period, and the comprehensive investigation of several parameters in the present study strengthen the characterization of the hepatic changes in obese Göttingen Minipigs with and without diabetes.
As expected, the FFC and FFCDIA groups developed obesity with higher BW, BF% and dyslipidemia (defined as increased plasma TG and TC) compared to the two groups fed standard diet. The FFCDIA group had significantly higher levels of fasting GLU and FRA confirming hyperglycemia. Plasma ALP and GLDH were increased, whereas an increase in circulating concentrations of ALT and AST was not found, in fact plasma ALT was lower in FFC pigs compared to SD and FFC/SD. Studies of Ossabaw minipigs on the same diet also found no difference in ALT compared to control animals, whereas AST was increased. However, in pigs, GLHD is considered a more reliable marker of acute liver damage, as compared to ALT, which is non-specific in pigs . Also, in human patients with NAFLD, ALT values have not been found to correlate with the degree of histopathological changes . Table 7 provides an overview of relevant NAFLD/NASH characteristics including selected circulating biomarkers from human and pig studies.
The LW was highly increased in both FFC and FFCDIA groups compared to SD, and despite elevated liver content of cholesterol in both groups, unexpectedly no statistically significant differences were found for triglycerides liver content among the four groups.
Furthermore, the degree of histopathological alterations was less than expected, especially regarding steatosis, despite the presence of metabolic disturbances indicated by the changes in circulating markers. Few animals displayed more than 5% macrovesicular steatosis which is a criterion for the diagnosis of NAFLD in humans [28, 29]. Rodent models on a high fat diet usually report abundant macro- and microvesicular steatosis [30, 31], but previous studies in Ossabaw minipigs on the same diet (5B4L) also reported lack of macrovesicular steatosis, despite extensive liver injury [10, 12]. Others have reported extensive microvesicular steatosis with little or no presence of macrovesicular steatosis in pigs on high fat plus high cholesterol or high sucrose diets [11, 13]. Microvesicular steatosis was rarely seen in our pigs, and the low level of both types of steatosis is consistent both with the biochemical analysis and the quantification of lipid on ORO stained slides. This modest accumulation of lipids in the liver could reflect the fact that the liver in pigs, in contrast to humans and rodents, is not the primary site of de novo lipogenesis . Also, the overnight fasting period before termination might have had some influence on the lipid content in the liver. However, so far it is unclear what is responsible for the significantly increased LW in FFC and FFCDIA groups and it needs to be clarified if hepatocyte hypertrophy or hyperplasia is the reason for the present hepatomegaly.
Hepatocellular ballooning, another criterion for the diagnosis of human NASH, was not seen in our study, but extensive cytoplasmic alterations in hepatocytes were present in FFC and FFCDIA groups. It was speculated that these alterations could be caused by an accumulation of glycogen. Glycogenic hepatopathy is a condition with massive cytoplasmic deposition of glycogen in hepatocytes leading to hepatomegaly, and is mostly seen in patients with poorly controlled type 1-diabetes. The diagnosis is confirmed with liver biopsies and staining for glycogen with PAS . In this study, staining for glycogen with PAS however did not show correlation with CA in either FFCDIA or FFC, and the FFCDIA group also showed decreased glycogen content in their liver tissue in comparison to the FFC group. Instead a strong association between CA and cholesterol content in the liver persisted, and perhaps CA could reflect a functional adaptation to the increased cholesterol load. In humans, it has recently been suggested that dietary cholesterol activates the hepatic stellate cells thereby promoting fibrosis especially if hepatocyte uptake or biliary excretion of cholesterol is inhibited .
Mallory–Denk bodies is a characteristic feature of NASH in humans, but can be difficult to identify and often additional IHC staining has to be performed, e.g., using ubiquitin or cytokeratin 8/18 antibody. Lackner et al. found that CK 8/18 was diminished or absent in ballooned hepatocytes compared to normal hepatocytes. They also found that Mallory–Denk bodies were not always present . Others suggest double immunohistochemical staining of CK8/18 and ubiquitin for the optimal detection of hepatocellular injury in human liver . In some mouse models of NASH, Mallory–Denk bodies have been reported morphologically [30, 36, 37], but until now no previous studies with pigs [10,11,12] or other rodent species [31, 38] have identified Mallory–Denk bodies along with hepatocellular ballooning. Immunohistochemical staining of CK8/18 or ubiquitin in liver tissue of porcine NAFLD/NASH models has not been reported. However, in a mouse model of NAFLD/NASH staining for CK8/18 has been used  and found less or absent staining in hepatocytes with ballooning.
The predominantly periportal fibrosis seen in our minipigs is coherent with findings in Ossabaw minipigs . This is in contrast to humans where perisinusoidal fibrosis around the central vein dominates and periportal fibrosis occurs only in later disease stages. A possible explanation for this peripheral deposition of collagen could be the localization of the porcine hepatic stellate cells. These cells are responsible for collagen production when activated and heterogeneity of stellate cells in the porcine liver was described by Wake et al.  who found that desmin positive stellate cells were most abundant in the peripheral regions of the classical lobules. Interestingly, two subtypes of NAFLD have been described in pediatric patients; type 1 resembles ‘adult’ NAFLD whereas type 2 is characterized by steatosis, portal inflammation and portal fibrosis  and is known to be the most prevalent type in children. The mechanism behind this pediatric portal fibrosis is unknown but perhaps a common pathway for portal fibrosis in children and pigs exist. The hepatic inflammation seen in FFC and FFCDIA groups was not a reflection of a systemic inflammatory state as the circulating inflammatory markers CRP and ALB showed no difference between the four groups.
Interestingly, the FFC/SD group also differed from SD in terms of fibrosis, despite being fed the same standard diet for the last 6 months. This indicates that fibrosis developed in the first 7 months of high fat/high cholesterol-feeding and either did not progress further during the 6 months of healthy dieting, or maybe even regressed.
Diabetes did not exacerbate the hepatic histopathology, although the FFCDIA pigs had substantial dyslipidemia with elevated levels of TC and TG and a liver tissue cholesterol content exceeding that of the FFC group. Previous unpublished studies indicated that 2% cholesterol was needed in order to achieve the desired elevated level of circulating total cholesterol in normal minipigs, whereas diabetic minipigs reached the same level on a 1% cholesterol diet, which was also the case here. Six diabetic pigs were terminated prematurely nearly reducing the group by half, which could have led to lack of statistical power when comparing the FFCDIA group with the FFC group.
Our diabetic animals had a type 1-like diabetes phenotype and were treated with long acting insulin analogue to keep blood glucose around 15 mM. A lower incidence of NAFLD has been found in patients with type 1 diabetes compared to patients with type 2 diabetes, and it is suggested that insulin treatment’s inhibiting effect on lipolysis is responsible for the decreased level of free fatty acids accumulating in the liver. Insulin resistance induced by excessive caloric intake is also known to play a role in hepatic lipid accumulation and the development of NAFLD. Only FFCDIA had elevated fasting blood glucose, but insulin resistance can be present in muscle and liver without hyperglycemia in human patients with NAFLD . Although insulin resistance was not directly assessed in the current study, KG, AUCinsulin and HOMA-IR values indicated that FFC and likely FFCDIA animals were insulin resistant on a whole-body level, but since no profound steatosis was present, hepatic insulin resistance is less likely to be a major factor in this model.
A limitation to this study was, that 11 animals ended the study prematurely and necropsy reports mentioned enlarged pale livers in all FFC (n = 2) and FFCDIA (n = 6) pigs, but no further analyses were performed. It could be speculated that these excluded animals especially from the FFCDIA group may have exhibited more pronounced hepatic changes thus biasing our results as only animals with less pronounced changes were able to complete the full study period. Another limitation is that no liver biopsies were taken at different time points during the study, especially before the FFC/SD group was changed to standard diet, making it difficult to conclude if the differences between FFC and FFC/SD were due to lack of progression or regression of marked changes already present at the intervention time point.
Overall, this diet-induced obese Göttingen Minipig model with or without diabetes poses some challenges in terms of translatability with human NASH, because some of the cardinal characteristics (abundant steatosis, hepatocellular ballooning and Mallory–Denk bodies, and zone 3 fibrosis) were missing. However, to explain some of the rather unexpected elements of this model, further studies are needed. Additional staining with e.g. CK8/18, ubiquitin, p62 or sonic hedgehog markers could be performed in order to validate if Mallory–Denk bodies (or their precursors) were present in this minipig model. Moreover, evaluation of serum CK8/18 may elucidate if depletion at tissue level could give as consequence an increase in the circulation of keratin fragments that are major components of Mallory–Denk bodies . By elucidating the pathological pathways perhaps combined with gene expression and electron microscopy, the mechanisms leading to CA in hepatocytes and the portal/periportal fibrosis could be further clarified.
In conclusion, a diet-induced obese Göttingen Minipig model with and without diabetes has been evaluated as a model of human NAFLD including NASH. Even though fibrosis, inflammation and cytoplasmic alterations occurred, cardinal histopathological characteristics such as hepatocellular ballooning and abundant steatosis were lacking despite the presence of overt obesity and dyslipidaemia. Furthermore, diabetes did not seem to exacerbate the hepatic changes developed by the high fat/fructose/cholesterol diet alone. The limited presence of key human-relevant pathological hepatic findings and variation in the model, limits its use in preclinical research without further optimisation.
- AUCinsulin :
area under the curve for insulin response
total body fat percentage
dual-energy x-ray absorptiometry
lobus hepatis dexter lateralis
lobus hepatis dexter medialis
high fat/high cholesterol diet group
- FFCDIA :
diet normalization group
hematoxylin and eosin stain
homeostasis model assessment of insulin resistance
ionized calcium binding adapter molecule 1
quantification of Iba1 stained inflammation
- KG :
intravenous glucose tolerance index
non-alcoholic fatty live disease
optimum cutting temperature compound
oil red O stain
quantification of ORO stained lipid
periodic acid-Schiff stain
picro-sirius red stain
quantification of PSR stained collagen
lean control group fed standard diet
lobus hepatis sinister lateralis
lobus hepatis sinister medialis
tissue micro arrays
Anstee QM, Targher G, Day CP. Progression of NAFLD to diabetes mellitus, cardiovascular disease or cirrhosis. Nat Rev Gastroenterol Hepatol. 2013;10:330–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.
Hazlehurst JM, Woods C, Marjot T, Cobbold JF, Tomlinson JW. Non-alcoholic fatty liver disease and diabetes. Metabolism. 2016;65:1096–108.
Targher G, Bertolini L, Padovani R, Rodella S, Tessari R, Zenari L, et al. Prevalence of nonalcoholic fatty liver disease and its association with cardiovascular disease among type 2 diabetic patients. Diabetes Care. 2007;30:1212–8.
Yasui K, Sumida Y, Mori Y, Mitsuyoshi H, Minami M, Itoh Y, et al. Nonalcoholic steatohepatitis and increased risk of chronic kidney disease. Metabolism. 2011;60:735–9.
Mantovani A, Mingolla L, Rigolon R, Pichiri I, Cavalieri V, Zoppini G, et al. Nonalcoholic fatty liver disease is independently associated with an increased incidence of cardiovascular disease in adult patients with type 1 diabetes. Int J Cardiol. 2016;225:387–91.
de Ledinghen V, Hiriart JB, Vergniol J, Merrouche W, Bedossa P, Paradis V. Controlled attenuation parameter (CAP) with the XL probe of the fibroscan (R): a comparative study with the M probe and liver biopsy. Dig Dis Sci. 2017;62:2569–77.
Van Herck MA, Vonghia L, Francque SM. Animal models of nonalcoholic fatty liver disease—a starter’s guide. Nutrients. 2017. https://doi.org/10.3390/nu9101072.
Lau JK, Zhang X, Yu J. Animal models of non-alcoholic fatty liver disease: current perspectives and recent advances. J Pathol. 2017;241:36–44.
Liang T, Alloosh M, Bell LN, Fullenkamp A, Saxena R, Van Alstine W, et al. Liver injury and fibrosis induced by dietary challenge in the Ossabaw miniature swine. PLoS ONE. 2015;67:78. https://doi.org/10.1371/journal.pone.0124173.
Yamada S, Kawaguchi H, Yamada T, Guo X, Matsuo K, Hamada T, et al. Cholic acid enhances visceral adiposity, atherosclerosis and nonalcoholic fatty liver disease in microminipigs. J Atheroscler Thromb. 2017;24:1150–66.
Lee L, Alloosh M, Saxena R, Van Alstine W, Watkins BA, Klaunig JE, et al. Nutritional model of steatohepatitis and metabolic syndrome in the Ossabaw miniature swine. Hepatology. 2009;50:56–67.
Yang SL, Xia JH, Zhang YY, Fan JG, Wang H, Yuan J, et al. Hyperinsulinemia shifted energy supply from glucose to ketone bodies in early nonalcoholic steatohepatitis from high-fat high-sucrose diet induced Bama minipigs. Sci Rep. 2015;5:13980.
Mentzel CM, Alkan F, Keinicke H, Jacobsen MJ, Gorodkin J, Fredholm M, et al. Joint profiling of miRNAs and mRNAs reveals miRNA mediated gene regulation in the Gottingen minipig obesity model. PLoS ONE. 2016. https://doi.org/10.1371/journal.pone.0167285.
Johansen T, Hansen HS, Richelsen B, Malmlof R. The obese Gottingen minipig as a model of the metabolic syndrome: dietary effects on obesity, insulin sensitivity, and growth hormone profile. Comp Med. 2001;51:150–5.
Christoffersen BO, Grand N, Golozoubova V, Svendsen O, Raun K. Gender-associated differences in metabolic syndrome-related parameters in Gottingen minipigs. Comp Med. 2007;57:493–504.
Andreasen LJ, Krog S, Ludvigsen TP, Nielsen OL, Moller JE, Christoffersen BO, et al. Dietary normalization from a fat, fructose and cholesterol-rich diet to chow limits the amount of myocardial collagen in a Gottingen Minipig model of obesity. Nutr Metab. 2018;15:64.
Gerrity RG, Natarajan R, Nadler JL, Kimsey T. Diabetes-induced accelerated atherosclerosis in swine. Diabetes. 2001;50:1654–65.
Matthews DR, Hosker JP, Rudenski AS, Naylor BA, Treacher DF, Turner RC. Homeostasis model assessment: insulin resistance and β-cell function from fasting plasma glucose and insulin concentrations in man. Diabetologia. 1985;28:412–9.
Pedersen SF, Ludvigsen TP, Johannesen HH, Lofgren J, Ripa RS, Hansen AE, et al. Feasibility of simultaneous PET/MR in diet-induced atherosclerotic minipig: a pilot study for translational imaging. Am J Nucl Med Mol Imaging. 2014;4:448–58.
Heegaard PM, Pedersen HG, Jensen AL, Boas U. A robust quantitative solid phase immunoassay for the acute phase protein C-reactive protein (CRP) based on cytidine 5′-diphosphocholine coupled dendrimers. J Immunol Methods. 2009;343:112–8.
Tveden-Nyborg P, Birck MM, Ipsen DH, Thiessen T, Feldmann LB, Lindblad MM. Diet-induced dyslipidemia leads to nonalcoholic fatty liver disease and oxidative stress in guinea pigs. Transl Res. 2016;168:146–60.
Kleiner DE, Brunt EM, Van Natta M, Behling C, Contos MJ, Cummings OW, et al. Design and validation of a histological scoring system for nonalcoholic fatty liver disease. Hepatology. 2005;41:1313–21.
Wijesundera KK, Izawa T, Tennakoon AH, Murakami H, Golbar HM, Katou-Ichikawa C, et al. M1- and M2-macrophage polarization in rat liver cirrhosis induced by thioacetamide (TAA), focusing on Iba1 and galectin-3. Exp Mol Pathol. 2014;96:382–92.
Lackner C, Gogg-Kamerer M, Zatloukal K, Stumptner C, Brunt EM, Denk H. Ballooned hepatocytes in steatohepatitis: the value of keratin immunohistochemistry for diagnosis. J Hepatol. 2008;48:821–8.
Jordan HL, Register TC, Tripathi NK, Bolliger AP, Everds N, Zelmanovic D, et al. Nontraditional applications in clinical pathology. Toxicol Pathol. 2014;42:1058–68.
Mofrad P, Contos MJ, Haque M, Sargeant C, Fisher RA, Luketic VA, et al. Clinical and histologic spectrum of nonalcoholic fatty liver disease associated with normal ALT values. Hepatology. 2003;37:1286–92.
Chalasani N, Younossi Z, Lavine JE, Charlton M, Cusi K, Rinella M, et al. The diagnosis and management of nonalcoholic fatty liver disease: practice guidance from the american association for the study of liver diseases. Hepatology. 2018;67:328–57.
Bedossa P. Utility and appropriateness of the fatty liver inhibition of progression (FLIP) algorithm and steatosis, activity, and fibrosis (SAF) score in the evaluation of biopsies of nonalcoholic fatty liver disease. Hepatology. 2014;60:565–75.
Asgharpour A, Cazanave SC, Pacana T, Seneshaw M, Vincent R, Banini BA, et al. A diet-induced animal model of non-alcoholic fatty liver disease and hepatocellular cancer. J Hepatol. 2016;65:579–88.
Ipsen DH, Tveden-Nyborg P, Rolin B, Rakipovski G, Beck M, Mortensen LW, et al. High-fat but not sucrose intake is essential for induction of dyslipidemia and non-alcoholic steatohepatitis in guinea pigs. Nutr Metab. 2016;13:51.
Bergen WG, Mersmann HJ. Comparative aspects of lipid metabolism: impact on contemporary research and use of animal models. J Nutr. 2005;135:2499–502.
Irani NR, Venugopal K, Kontorinis N, Lee M, Sinniah R, Bates TR. Glycogenic hepatopathy is an under-recognised cause of hepatomegaly and elevated liver transaminases in type 1 diabetes mellitus. Intern Med J. 2015;45:777–9.
Schwabe RF, Maher JJ. Lipids in liver disease: looking beyond steatosis. Gastroenterology. 2012;142:8–11.
Guy CD, Suzuki A, Burchette JL, Brunt EM, Abdelmalek MF, Cardona D, et al. Costaining for keratins 8/18 plus ubiquitin improves detection of hepatocyte injury in nonalcoholic fatty liver disease. Hum Pathol. 2012;43:790–800.
Nishida T, Tsuneyama K, Fujimoto M, Nomoto K, Hayashi S, Miwa S, et al. Spontaneous onset of nonalcoholic steatohepatitis and hepatocellular carcinoma in a mouse model of metabolic syndrome. Lab Invest. 2013;93:230–41.
Kucukoglu O, Guldiken N, Chen Y, Usachov V, El-Heliebi A, Haybaeck J, et al. High-fat diet triggers Mallory–Denk body formation through misfolding and crosslinking of excess keratin 8. Hepatology. 2014;60:169–78.
Ichimura M, Masuzumi M, Kawase M, Sakaki M, Tamaru S, Nagata Y, et al. A diet-induced Sprague–Dawley rat model of nonalcoholic steatohepatitis-related cirrhosis. J Nutr Biochem. 2017;40:62–9.
Wake K, Sato T. Intralobular heterogeneity of perisinusoidal stellate cells in porcine liver. Cell Tissue Res. 1993;273:227–37.
Schwimmer JB, Behling C, Newbury R, Deutsch R, Nievergelt C, Schork NJ, et al. Histopathology of pediatric nonalcoholic fatty liver disease. Hepatology. 2005;42:641–9.
Gastaldelli A. Insulin resistance and reduced metabolic flexibility: cause or consequence of NAFLD? Clin Sci. 2017;131:2701–4.
Strnad P, Paschke S, Jang KH, Ku NO. Keratins: markers and modulators of liver disease. Curr Opin Gastroenterol. 2012;28:209–16.
Li SJ, Ding ST, Mersmann HJ, Chu CH, Hsu CD, Chen CY. A nutritional nonalcoholic steatohepatitis minipig model. J Nutr Biochem. 2016;28:51–60.
CSP: planning and performing experiment, data analyses including histopathological assessment and imaging analyses, interpretation, writing manuscript. BØC: study design, planning experiment, interpretation, writing manuscript. RKKI: planning and supervising histopathological assessment and imaging analyses, interpretation, writing manuscript. TPL: study design, planning and performing experiment, data analyses, interpretation, writing manuscript. NEZ: interpretation, writing manuscript. HDP: study design, planning experiment, interpretation, writing manuscript. MV: planning and supervising histopathological assessment, interpretation of histopathological data, writing of manuscript. LHO: study design, planning of experiment, supervising data analyses, interpretation, writing manuscript. All authors read and approved the final manuscript.
Novo Nordisk A/S and the Lifepharm Centre at the University of Copenhagen are acknowledged for funding the study and PhD grant for CSP. Acknowledgement is given to Belinda Bringtoft, Betina Gjedsted Andersen, Christina Tirsdal Kjempff, Jeanette Jensen, University of Copenhagen, Anders Hagedal Uhrenfeldt, Anita Svendsen, Ann-Charlott Kemp, Bettina Brandrup, Camilla Ingvorsen, Chalotte Lasak Krøyer Pedersen, Hanne Toftelund, Helle Wagner, Inge Rubach, Johnny Kaltoft Nisted, Kamilla Larsen, Lotte Schmidt Marcher, Pernille Kristiane Pedersen, Susanne Juul Rasmussen, student helpers, Novo Nordisk A/S for skilled technical assistance performing the experiment or analyses. Torben Martinussen, University of Copenhagen, is acknowledged for helpful guidance on statistical analyses. Acknowledged is also Marianne Kjestine Petersen, University of Copenhagen, for everyday technical support. Fellow Ph.D. students Rikke Lindgaard Thomsen and Simone Krog Christensen are acknowledged for valuable assistance performing the experiment and support during data analyses and manuscript writing.
HDP is full time employed at Ellegaard Göttingen Minipigs A/S.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Consent for publication
Ethics approval and consent to participate
This study using a porcine model was approved by the Danish authority under the Danish Ministry of Environment and food.
Financial support: Novo Nordisk A/S, Måløv, Denmark and LIFEPHARM—Novo Nordisk & LIFE In Vivo Pharmacology Centre, University of Copenhagen, Denmark have funded this study.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.