1. Introduction
Non-alcoholic fatty liver disease (NAFLD) is the most prevalent hepatic disease worldwide, and it is highly associated to metabolic disorders including obesity, insulin resistance, diabetes, and metabolic syndrome [
1,
2,
3]. NAFLD comprises a spectrum of histopathological changes characterized by abnormally high accumulation of lipids in the liver found in simple steatosis (SS) that progresses to steatohepatitis (NASH) where inflammation and ballooning are observed [
1]. Fibrosis might be present or absent during both simple steatosis and NASH [
4], however factors determining which patients develop fibrosis are not clear, although metabolic, endocrine, genetic and aging have been suggested among them [
1]. Cirrhosis, hepatocellular carcinoma, and liver failure are the hepatic endpoints of this disease, but the presence of severe fibrosis is considered the most important predictor of NAFLD outcome and death risk [
5,
6,
7].
Insulin-like growth factor (IGF) binding proteins (IGFBP) is a family of secreted proteins whose primary function is binding IGFs in the bloodstream, regulating its bioavailability and half-life [
8,
9]. IGF-independent functions including proliferation, apoptosis and cellular senescence have been described [
8,
9]. These proteins are expressed in most tissues, but the liver is the main source for most of them [
9]. A role during chronic liver disease has been suggested for some IGFBP. IGFBP-1 expression is increased in the liver during post-injury regeneration [
10] and is considered a hepatoprotective factor that prevents apoptosis in hepatocytes [
11]. In humans, IGFBP-3 expression is reduced in hepatocellular carcinoma, compared with cirrhotic tissue. IGFBP-3 is regulated by transforming growth factor (TGF)-β promoting hepatic stellate cells (HSC) migration in vitro and increasing portal hypertension in the bile duct ligation model [
12], and might be a key protein during alcoholic liver disease by inducing lipid droplet and triglyceride accumulation in vitro [
13]. IGFBP-7, a low IGF-affinity protein, is mainly expressed in HSC, being upregulated during transdifferentiation [
14], and activation [
15,
16]. Expression of IGFBP-7 is significantly elevated in liver biopsies from patients with fibrosis and cirrhosis [
17], whereas inhibition or silencing of this protein prevents accumulation of extracellular matrix in experimental models of liver fibrosis [
16,
17]. IGFBP-7 has shown to induce cellular senescence [
18]; in fact, IGFBP-7 expression in HCC tissue is lower compared with healthy tissue [
19]. Regarding metabolic abnormalities including NAFLD; low levels of IGFBP-2 have been reported in obesity, type 2 diabetes and metabolic syndrome [
20], whereas IGFBP-5 is increased in NASH [
21]. Serum IGFBP-3 is decreased in patients with NAFLD [
22] whereas IGF-1/IGFBP-3 ratio has been associated with lower likelihood of NAFLD and lower grade steatosis [
23] and, histopathologic features of the liver biopsy including ballooning and inflammation [
24]. Interestingly, increased expression of IGFBP-1 and IGFBP-7 in the liver might contribute to hepatic insulin resistance [
25] and further fat accumulation.
Accordingly, IGFBPs might have a role during NAFLD, as well as on the development of fibrosis, however this is not clear. The aim of the study was to assess IGFBPs in both serum and liver of mice during the progression of NAFLD and the onset of fibrosis in this disease, we hypothesized that IGFBP, through their IGF-independent actions might have an actual role in the development of the disease.
2. Materials and Methods
2.1. NAFLD induction
Male C57BL/6 mice of 16 weeks of age and weighing 23 ± 3 g were obtained from the Animal Care Facilities at the Experimental Medicine Unit from School of Medicine, UNAM at Hospital General de México, and maintained under controlled conditions. Both food and water were allowed ad libitum. All animals received humane care; all procedures were approved by the Institutional Committee of Care and Use of Laboratory Animals (FM/DI/005/2022 approved on September 6, 2022) and agree with the national guidelines and the ARRIVE guidelines for animal use and care in research. In order to develop different stages of NAFLD, mice were randomly assigned to be fed either a Methionine-Choline Deficient (MCD) or a Methionine-Choline Complete (Control) Diet (MP Biomedicals, CA, USA), during 2 (MCD 2w, n=19; Control 2w, n=7), 8 (MCD 8w, n=17; Control 8w, n=9) or 12 weeks (MCD 12w, n=20; Control 12w, n=8).
2.2. Sample collection
After appropriate time, mice were anesthetized with Xylazine-Ketamine. Blood was collected and allowed to clot at 4°C; serum was obtained by centrifugation at 750 xg for 10 minutes and stored at -80°C until assayed. Liver samples from the left and medial lobes were collected, snap frozen in liquid nitrogen and stored at -80°C until assayed. Samples from the left lobe were divided to either be fixed in 3.7% formaldehyde-Phosphate buffered saline (PBS) and embedded in paraffin or, to be embedded in Tissue-tek OCT (Sakura Finetek, CA, USA) and stored at -20°C until assayed.
2.3. Histological assessment
Histological evaluation of fatty liver disease was performed in Hematoxylin-Eosin-stained sections according to the Non-alcoholic fatty liver disease Activity Score (NAS) [
26]. Fibrosis was assessed in Masson’s Trichrome stained sections and classified according to their fibrosis degree [
26]. Liver fat contents were assessed in frozen sections stained with Oil-Red O (Abcam, MA, USA) and quantified by a morphometric analysis using Image J v1.53k software (NIH, MD, USA) according to [
27]. Briefly: 10 optic fields were captured, and the percentage of red-stained area was calculated as a percentage of the complete area of the optic field.
2.4. Senescence analysis
Cellular senescence was assessed in frozen sections embedded in Tissue-Tek OCT using the Senescence Detection Kit (Abcam, MA, USA). A blue stain was indicative of the activity of the senescence-associated β-galactosidase (SA-β-gal). A morphometric analysis was performed by using Image J v1.53k software.
2.5. Liver protein isolation
Total protein was obtained from each sample using the PBS-Protease Inhibitor Cocktail Set III (Calbiochem, Darmstadt, Germany). Protein integrity was tested by SDS-PAGE prior to the assays.
2.6. IGFBP assessment
IGFBP-1, -2, -3, -5, -6, and -7 from the livers and serum were quantified by multiple suspension arrays (Milliplex MAP, MIGFBPMAG43K, Merck-Millipore, Billerica, MA, USA) according to manufacturer instructions. Bead regions were assigned as follows: IGFBP-1: 27; IGFBP-2: 39; IGFBP-3: 42; IGFBP-5: 55; IGFBP-6: 61; IGFBP-2: 72.
2.7. Statistics
Data were analyzed by SSPS v22 and presented as Mean ± standard deviation (SD). Subjects receiving control diet did not exhibited differences attributable to time of exposure to the diet or aging and were analyzed as one single control group. Subjects receiving MCD diet were categorized according to NAS or to fibrosis stages. One-way ANOVA followed by Tukey’s post-hoc test was performed. Independent Linear regression models for the stage of NAFLD and fibrosis were used to explore the associations between time of exposure to MCD diet, the degree of fibrosis, and the levels of IGFBP (both crude and adjusted by hepatic lipid content). P<0.05 was considered significant.
4. Discussion
IGFBP superfamily has been related to insulin resistance [
25,
34], fat deposition and other cellular processes including senescence [
9]. All these factors have a clear role in the development of NAFLD. Here we assessed a group of IGFBP in both liver and serum according to the progression of NAFLD and fibrosis in the MCD mouse model. Our results show significantly increased serum levels of IGFBP-1, and -2 in mice with NAFLD and fibrosis. IGFBP-7, a low affinity IGFBP, was increased in serum as well as in the liver tissue; interestingly, it was elevated in serum in bNASH and NASH as well as during fibrosis in F1C and F2, suggesting a role for IGFBP-7 in the progression of NAFLD and the onset of fibrosis in this model. Linear regression models confirmed the predictive value of IGFBP-1, -2 and -7 for both, NAFLD stages and fibrosis degree.
As expected, we observed higher fat contents in the liver of all NAFLD subjects compared with controls, although no differences in the amounts of fat were observed among the stages of the disease (
Figure 2b), or among the different degrees of fibrosis observed: F0 to F2 (
Figure 2c).
For IGFBP-1, we observed differences in serum, related to both NAFLD and fibrosis (
Figure 3b, 4b). This protein is known to possess a hepatoprotective role [
11]; IGFBP-1 is among the first genes over expressed after a partial hepatectomy [
10], it also increases during liver disease [
35,
36]. IGFBP-1 is transcriptionally regulated by insulin [
37], and has been linked, as well as IGFBP-2, with insulin sensitivity [
38]. On the other hand, it has been shown that hedgehog, a well-known fibrogenic pathway, also regulates IGFBP-1 the expression [
39], suggesting this peptide might also be a fibrogenic mediator in the liver. Serum levels of IGFBP-1 were a significant crude and adjusted predictor for SS, bNASH and NASH (
Table 1). Regarding fibrosis, hepatic IGFBP-1 was an adjusted predictor for F0 and F2, whereas serum levels of this protein were significant predictors for F0, F1C and F2. We consider these associations to be related to both, hepatic lipid contents and the onset of liver fibrosis, according to Hagstrom et. al, IGFBP-1 levels might increase even more as fibrosis progresses [
29].
IGFBP-2 has been extensively studied in obesity, type 2 diabetes, insulin resistance and metabolic syndrome [
34,
40,
41,
42,
43]. In fact, this peptide has been suggested as a possible protector against obesity and insulin resistance [
44]. However, not much data is available during NAFLD. Here, we observed an increase in IGFBP-2 during NAFLD and fibrosis, in both liver and serum. Interestingly, during fibrosis IGFBP-2 was significantly elevated in the liver during F0 compared with F1C and F2, however both stages exhibited higher expression compared with controls. Hepatic expression of IGFBP-2 was a significant crude and adjusted predictor of SS, bNASH and NASH, whereas its serum levels were crude predictor as well. The increased levels of IGFBP-2 observed in the MCD-induced NAFLD might be explained by the fact that this model is characterized by low bodyweight and lacks insulin and leptin resistance. However, it is noteworthy the effect observed in the liver during fibrosis, where this peptide was a significant predictor of fibrosis progression in both, crude and adjusted by lipid models. Although IGFBP-2 has been suggested to have a role in idiopathic pulmonary fibrosis [
30,
31], its role in hepatic fibrosis is not clear.
Regarding IGFBP-7, this peptide is a well-known inductor of cellular senescence [
45,
46,
47], and it is implicated in HSC activation, which in turn produce ECM accumulation and fibrosis [
15,
16]. Cellular senescence has also been related to HSC activation and NAFLD as well [
33]. IGFBP-7 was significantly increased in the serum of subjects with bNASH and NASH, compared with controls (
Figure 3f). IGFBP-7 progressively increased in both liver and serum according to fibrosis; in the liver its expression was significantly higher in F2 (
Figure 4e) whereas its serum levels rose from F1C (
Figure 4f). When adjusted by lipid content in the liver, serum IGFBP-7 was a predictor of NASH (
Table 1). Regarding fibrosis, IGFBP-7 was a predictor (crude and adjusted by lipids) of F2 (
Table 2).
IGFBP-7 levels were in accordance with cellular senescence, assessed as the percentage of activity of SA-β-gal. During NAFLD progression, higher percentages of senescent cells were detected in bNASH and NASH compared with controls and SS. However, SS also exhibited increased levels of senescent cells compared with controls (
Figure 5a). In a similar manner, SA-β-gal activity was augmented as fibrosis progressed from mild to moderate, but high levels of senescent cells were also observed in F0 (
Figure 5b). Fat depots have also been implicated in the induction of cellular senescence in the affected hepatocytes as a consequence of lipotoxicity, increased oxidative stress as well as DNA damage and telomere erosion [
33]. Our data agrees with the association of increased fat contents in hepatocytes and exhibiting higher percentage of senescence [
33]. However, as seen in SS and F0, hepatic fat depots are not enough to induce the progression of NAFLD, nor fibrosis, instead the increase in IGFBP7 followed by cellular senescence might be a trigger to progress to NASH as well as the onset of ECM accumulation, by activating HSC. Senescence is known to increase in aging organs and tissues [
48], however we do not consider aging a factor affecting our results, since we did not observed differences in SA-β-gal activity in the livers of mice fed control diet during the 2, 8 or 12 weeks, even though they were 6-10 weeks older at sample collection, in fact they were reported as a single control group. Cellular senescence is involved in a range of chronic liver disorders including viral hepatitis B and C [
49], alcoholic liver disease [
32], genetic haemochromatosis [
50] and NAFLD [
33], however the mechanisms inducing senescence in chronic liver disease remain unclear. NAFLD and its derived fibrosis occur simultaneously; however, not all subjects progress at the same rate and fibrosis is not strictly associated to NASH [
4]. Serum levels of IGFBP-7 were significantly increased in the subjects with fibrosis (F1C and F2), but its hepatic expression was only increased in the F2 livers. Regarding senescence and fibrosis, we observed higher percentages of activity of SA-β-gal in the liver from all MCD-fed subjects that increased with the fibrosis progression.
Several factors determine which subjects develop fibrosis during NAFLD, the well documented oxidative stress [
51], endoplasmic reticulum stress, and, as suggested by our findings, IGFBP-7 and its resulting cellular senescence. Other features of NAFLD, including inflammation and ballooning, might be implicated in the increased senescence observed in bNASH and NASH. In our results increased IGFBP-1, -2, -7 and senescence were the most important determinants for the progression to NASH and fibrosis in the MCD-induced NAFLD mouse model.
One limitation of our study is that we observed only mild to moderate fibrosis; further studies in severer stages of fibrosis, cirrhosis and hepatocellular carcinoma derived from NAFLD are needed in order to establish a complete association between fibrosis, IGFBP-7 and senescence during NAFLD. Another limitation is in the MCD model, where NAFLD occurs in the absence of metabolic affection, however we consider our data as valuable since the presence of such metabolic affections might also influence IGFBP expression as well as senescence being a confounding factor when relating them to hepatocellular lipid content and fibrosis.
Author Contributions
“Conceptualization, C.G. and G.G-R.; methodology, C.G., A.C-E, M.G.B-U, I.I.R-B., G.G-R; validation, C.G., D.S-V.; formal analysis, C.G., A.C-E, M.G.B-U, I.I.R-B., D.S-V., G. G-R; investigation, C.G., G.G-R; resources, C.G.; data curation, C.G., A.C-E, M.G.B-U, I.I.R-B., D.S-V., G.G-R; writing—original draft preparation, C.G., D.S-V.; writing—review and editing, C.G., A.C-E., I.I.R-B., D.S-V., G.G-R.; project administration, C.G., G.G-R; funding acquisition, C.G. All authors have read and agreed to the published version of the manuscript.