Loss of Glp2r signaling activates hepatic stellate cells and exacerbates diet-induced steatohepatitis in mice

findings demonstrate that expression of the hepatic Glp2r is induced following high-fat A glucagon-like peptide-2 (GLP-2) analog is used in individuals with intestinal failure who are at risk for liver disease, yet the hepatic actions of GLP-2 are not understood. Treatment of high-fat diet–fed (HFD-fed) mice with GLP-2 did not modify the development of hepatosteatosis or hepatic inflammation. In contrast, Glp2r –/– mice exhibited increased hepatic lipid accumulation, deterioration in glucose tolerance, and upregulation of biomarkers of hepatic inflammation. Both mouse and human liver expressed the canonical GLP-2 receptor (GLP-2R), and hepatic Glp2r expression was upregulated in mice with hepatosteatosis. Cell fractionation localized the Glp2r to hepatic stellate cells (HSCs), and markers of HSC activation and fibrosis were increased in livers of Glp2r –/– mice. Moreover, GLP-2 directly modulated gene expression in isolated HSCs ex vivo. Taken together, these findings define an essential role for the GLP-2R in hepatic adaptation to nutrient excess and unveil a gut hormone-HSC axis, linking GLP-2R signaling to control of HSC activation. Furthermore, we localize a functional GLP-2R to hepatic stellate cells (HSCs) and demonstrate increased markers of stellate cell activation in Glp2r –/– mice. Collectively, these findings reveal a prandial signal trans-duction pathway linking gut-derived GLP-2 to the control of HSC activity, with implications for understanding GLP-2 action in the context of nutrition-associated hepatocellular injury. D (7-AAD) (BD Pharmingen, 559925), and purified by flow cytometry performed by the flow cytometry core at The Centre for Phenogenomics using MoFlo Astrios EQ (Beckman Coulter) equipped with 405 nm, 488 nm, 561 nm, and 640 nm lasers and running Summit v.6.2 software. Live stellate cells were identified by the following fluorescent signal combination: low 7-AAD (detector parameter 488–620/29 nm), high 405 nm violet laser autofluorescence (detector parameter 405–448/59 nm) and low 488 nm blue laser autofluorescence (detector parameter 488–513/26 nm). FACS gates were based on stellate cell autofluorescence lipid droplets. Additional gating to remove FITC range enhanced + – reproduc-ibly above 95%.


Introduction
Glucagon-like peptide-2 (GLP-2) is a 33-amino acid proglucagon-derived peptide that is secreted, together with GLP-1, from the intestinal enteroendocrine L cell (1). GLP-2 is continuously secreted at low basal levels in the fasting or interprandial state, and circulating levels rise briskly within minutes of food ingestion (2). GLP-2 exerts proliferative and cytoprotective actions in the small bowel, resulting in the expansion of the mucosal surface area and enhanced nutrient absorption (3,4). Sustained augmentation of GLP-2R signaling with a degradation-resistant GLP-2 analog, teduglutide, reduces the requirement for parenteral nutrition (PN) in human subjects with short-bowel syndrome (SBS) (5), and teduglutide is approved for the therapy of adults and children with PN-dependent SBS (6).
Although the actions of GLP-2 are predominantly localized to the gut (7), the extraintestinal expression of the GLP-2R raises important questions surrounding the safety of therapy with GLP-2R agonists, and potential actions of teduglutide beyond the gastrointestinal tract. For example, we recently detected robust expression of the GLP-2R in gallbladder and demonstrated that administration of teduglutide markedly inhibits gallbladder emptying in mice (8). These findings, subsequently extended to human studies, are consistent with reports of increased numbers of gallbladder-related adverse events associated with teduglutide therapy (9).
A full-length GLP-2 receptor (GLP-2R) mRNA transcript has also been identified in murine liver, albeit at levels much lower than in the jejunum (10). The GLP-2R was localized to hepatocytes (11); however, GLP-2R mRNA transcripts were not identified within hepatocytes profiled by single-cell RNA sequencing (12). Moreover, the biological actions of GLP-2 on the liver are conflicting, with some preclinical studies demonstrating a detrimental effect on hepatic steatosis (13), whereas other reports highlight a hepatoprotective role for GLP-2, characterized by a reduction in markers of systemic and hepatic inflammation (14). The frequent development of hepatic steatosis in PN-dependent subjects with SBS, associated in some subjects with progression to hepatic failure and liver transplantation, highlights the importance of understanding the liver as a target of GLP-2 action (15,16).
To resolve the hepatic actions of GLP-2, we have now studied GLP-2 action in high-fat, high-fructose, high-cholesterol-fed (HFHC-fed) WT mice and simultaneously examined the hepatic phenotype of Glp2r -/mice. Our findings demonstrate that expression of the hepatic Glp2r is induced following high-fat A glucagon-like peptide-2 (GLP-2) analog is used in individuals with intestinal failure who are at risk for liver disease, yet the hepatic actions of GLP-2 are not understood. Treatment of highfat diet-fed (HFD-fed) mice with GLP-2 did not modify the development of hepatosteatosis or hepatic inflammation. In contrast, Glp2r -/mice exhibited increased hepatic lipid accumulation, deterioration in glucose tolerance, and upregulation of biomarkers of hepatic inflammation. Both mouse and human liver expressed the canonical GLP-2 receptor (GLP-2R), and hepatic Glp2r expression was upregulated in mice with hepatosteatosis. Cell fractionation localized the Glp2r to hepatic stellate cells (HSCs), and markers of HSC activation and fibrosis were increased in livers of Glp2r -/mice. Moreover, GLP-2 directly modulated gene expression in isolated HSCs ex vivo. Taken together, these findings define an essential role for the GLP-2R in hepatic adaptation to nutrient excess and unveil a gut hormone-HSC axis, linking GLP-2R signaling to control of HSC activation. diet (HFD) feeding; however, sustained GLP-2R agonism has little effect on the development of hepatic steatosis. In contrast, loss of the GLP-2R increases steatosis and biomarkers of hepatic inflammation. Furthermore, we localize a functional GLP-2R to hepatic stellate cells (HSCs) and demonstrate increased markers of stellate cell activation in Glp2r -/mice. Collectively, these findings reveal a prandial signal transduction pathway linking gut-derived GLP-2 to the control of HSC activity, with implications for understanding GLP-2 action in the context of nutrition-associated hepatocellular injury.

Results
GLP-2 does not alter the metabolic and inflammatory sequelae of HFHC feeding. Previous studies of metabolic and hepatic health in HFD-fed mice and rats treated with GLP-2 have yielded conflicting results (13,17).
To study the effect of GLP-2 in the setting of energy excess, we administered GLP-2 to mice on (a) a HFHC diet or (b) a 10% fat control diet (CD) matched for other macronutrient sources. On week 15, mice were randomly allocated to treatment with either the long-acting GLP-2 analog h[Gly2]GLP-2, which is identical to teduglutide (hereafter referred to as GLP-2), or saline, for 11 days ( Figure 1A). Body and liver weight were increased by the HFHC diet but were not different in GLP-2-treated mice ( Figure 1B and Supplemental Figure 1A; supplemental material available online with this article; https://doi.org/10.1172/ jci.insight.136907DS1). GLP-2 increased small bowel length, total intestinal weight, and weight per cm as well as gallbladder volume (Supplemental Figure 1, B and C). Glucose tolerance (i.p.) and plasma lipid profiles (LDL, HDL, total cholesterol, and triglyceride [TG]) were unaffected by GLP-2 treatment ( Figure  1, C-E). Hepatic Oil Red O staining ( Figure 1F) and hepatic cholesterol and TG content ( Figure 1G) were increased after exposure to HFHC diet but not different in GLP-2-treated mice. Levels of hepatic mRNA transcripts encoding proteins important for lipid metabolic pathways were frequently modulated by diet but not altered by GLP-2 treatment ( Figure 1H).
Plasma alanine transferase (ALT) levels were higher ( Figure 2A) and hepatic and plasma levels of KC/ GRO, TNF-α, and IL-10 were elevated in HFHC-fed mice; however, none of the measured hepatic or plasma cytokines, including IL-1β, IFN-γ, IL-6, and IL-5, were different in GLP-2-treated mice ( Figure 2B and Supplemental Figure 1D). Similarly, hepatic mRNA transcripts corresponding to markers of inflammation and fibrosis were not different after GLP-2 administration ( Figure 2C). Thus, in HFHC diet-induced steatohepatitis, GLP-2 treatment did not impair indices of glucose and lipid metabolism or hepatic and systemic inflammation.
Glp2r -/mice exhibit impaired glucose tolerance and enhanced hepatic fat deposition. Previous studies have demonstrated that both gain and loss of GLP-2 action promoted hepatic steatosis in HFD-fed mice (13,17). To resolve this discrepancy, we assessed indices of glucose and lipid metabolism in male Glp2r -/and Glp2r +/+ littermate control mice fed a HFHC or CD for 17 weeks ( Figure 3A). HFHC-fed mice exhibited increased weight gain and liver mass; however, these parameters were not different in Glp2r -/mice ( Figure  3B and Supplemental Figure 2A). Glp2r -/mice demonstrated elevated AUC glucose during the second hour of i.p. glucose tolerance ( Figure 3C). Fasting plasma lipid profiles, including LDL, HDL, total cholesterol, and TG, were not different in Glp2r -/versus Glp2r +/+ mice ( Figure 3, D and E).
CD-fed Glp2r -/mice exhibited greater hepatic Oil Red O staining ( Figure 3F) and hepatic cholesterol and TG content were elevated in HFHC diet-fed mice, with no difference between genotypes. However, hepatic cholesterol and TG levels were increased in CD-fed Glp2r -/versus Glp2r +/+ mice, despite similar body weights in both genotypes (Figure 3, B and G). We next assessed levels of mRNA transcripts corresponding to genes important for TG and cholesterol uptake and degradation and de novo lipogenesis. Although the majority of hepatic mRNA transcript levels were not different between genotypes, levels of Hmgcr, Cpt1a, Hadhb, and Pparg were increased in Glp2r -/mice ( Figure 3H).
To assess extrahepatic mechanisms that could modify the extent of hepatic lipid deposition in Glp2r -/mice, we studied oral lipid absorption as well as dipeptidyl peptidase-4 (DPP-4) activity and prandial GLP-1 levels in a separate cohort of mice. Plasma DPP-4 activity was not different (Supplemental Figure 2C), whereas fasting GLP-1 levels trended higher in Glp2r -/mice (P = 0.07) (Supplemental Figure 2D). Plasma TGs were not different in the fasting state or after olive oil challenge (Supplemental Figure 2, E and F); however, plasma GLP-1 levels were higher in Glp2r -/mice after enteral olive oil administration (Supplemental Figure 2G). Accordingly, putative differences in DPP-4 activity and increased GLP-1 levels, which might also reduce postprandial lipemia (18,19), do not likely contribute to elevated lipid deposition and dysglycemia in CD-fed Glp2r -/mice. Hence, Glp2r -/mice develop dysglycemia and increased hepatic lipid deposition, which could not be attributed to known extrahepatic aspects of GLP-2R signaling.
Glp2r is expressed in HSCs. To understand how loss of the GLP-2R might predispose to increased hepatic inflammation and lipid deposition, we assessed hepatic Glp2r expression. The liver expresses low levels of mRNA transcripts corresponding to a full-length GLP-2R (8, 10), previously suggested to be localized to mouse and human hepatocytes (11). We detected expression of full-length GLP2R mRNA transcripts in 3 commercially available human liver total RNA samples ( Figure 5A). The hepatocyte fraction isolated from mouse liver, confirmed by cell-specific marker expression (Supplemental Figure 3A), did not contain Glp2r mRNA transcripts, whereas the nonparenchymal cell (NPC) fraction contained detectable Glp2r mRNA transcripts ( Figure 5B). Since HSCs are not efficiently recovered by standard NPC isolation protocols, we analyzed highly purified HSC fractions (20) ( Figure 5, C-E). The purity of HSC-sorted fractions was confirmed through fluorescent microscopy ( Figure 5E) and mRNA abundance of markers for HSC and nonstellate cell fractions ( Figure 5D). A substantial enrichment of Glp2r expression was observed in highly purified HSCs relative to levels in total liver RNA ( Figure 5F), with expression levels in HSCs approximately 50% of those detected in the jejunum. Conventional RT-PCR followed by Southern blotting identified a full-length Glp2r mRNA transcript in RNA from mouse liver and stellate cells as well as in jejunum ( Figure 5G). There was no detectable expression of full-length Glp2r in RNA isolated from hepatocytes ( Figure 5G and Supplemental Figure 3C).
To probe whether hepatic NPCs other than stellate cells express Glp2r, we used flow cytometry (Supplemental Figure 3B) to sort NPC subpopulations using CD26, CD31, and CD45 surface markers. CD31 + CD45 -(endothelial cells), CD45 + CD31 -(subdivided by green autofluorescence properties to Kupffer cells and other hematopoietic cells), and CD31 -CD45 -CD26 + (comprising cholangiocytes and other nonendothelial, nonhematopoietic NPCs, with the exception of HSCs that are CD26 -) were further identified by (C) Hepatic mRNA abundance, relative to Ppia, of genes relevant to inflammatory pathways. Data are presented as the mean ± SD. *P < 0.05, **P < 0.01, ***P < 0.001, § P < 0.05, significant effect of diet on gene expression variance, using 2-way ANOVA with Tukey correction for multiple comparisons. Hepatic mRNA abundance, relative to Ppia, of genes affecting lipid metabolism pathways. Graphs on the right illustrate relative mRNA levels of Hmgcr, Pparg, Cpt1a, and Hadhb genes, for which genotype had a significant effect on gene expression variability. Data are presented as the mean ± SD. *P < 0.05, **P < 0.01, ***P < 0.001, § P < 0.05 significant effect of diet, # P < 0.05, significant effect of genotype on variance of gene expression, using 2-way ANOVA with Tukey correction for multiple comparisons. The HSC GLP-2R is directly activated by GLP-2 in vitro. Activation of GLP-2R signaling rapidly induces immediate-early gene (IEG) expression in the intestine as well as in heterologous cells expressing the GLP-2R (21). Exposure of isolated murine HSCs to GLP-2 increased levels of Egr1, Nr4a, Rasd1, and Phlda1, consistent with a functional HSC GLP-2R ( Figure 5H and Supplemental Figure 4A).
Regulation of hepatic Glp2r expression. We next analyzed hepatic Glp2r mRNA expression from different groups of mice exposed to various diets. Mice fed a HFHC diet for 17 weeks exhibited an approximately 2-fold increase in hepatic Glp2r expression (Supplemental Figure 4B). A similar increase in hepatic Glp2r mRNA transcripts was detected in mice after 12 weeks of 45% HFD feeding (Supplemental Figure 4B) as well as in mice on a methionine-choline-deficient (MCD) diet. Notably, mRNA levels of the stellate cell marker Pdgfrb, a reflection of HSC abundance, were not different (Supplemental Figure 4B).
Genotype-dependent differences in HSC mRNA transcripts. We next analyzed levels of mRNA transcripts for 11 genes reported to be upregulated in the context of HSC activation (22)(23)(24)(25) in livers from mice fed a MCD diet for 2 weeks. Three of these markers, Timp1, Mmp13, and Adora2b, were upregulated in whole Hepatic mRNA abundance, relative to Ppia, of genes relevant to inflammatory pathways. Graphs on the right illustrate relative mRNA levels of Oas2, Pkr, Cxcr2, and Crp genes, for which genotype had a significant effect on gene expression variability. Data are presented as the mean ± SD. *P < 0.05, **P < 0.01, ***P < 0.001, ****P < 0.001, # P < 0.05 significant effect of genotype on variance of gene expression, using 2-way ANOVA with Tukey correction for multiple comparisons. (G) Full-length mouse Glp2r expression in jejunum, whole liver, hepatocytes, and isolated HSCs. Glp2r mRNA transcripts were amplified by reverse transcription PCR with primers that amplify a 1.65-kb product encompassing the entire Glp2r open reading frame. PCR products were analyzed by agarose gel electrophoresis followed by SYBR green staining (top and bottom). Southern blotting with an internal Glp2r 32 -P-labeled oligonucleotide (middle) was used to detect the PCR product. (H) mRNA abundance, relative to Ppia, of immediate-early genes used as markers of GLP-2 action in isolated HSCs. Hepatic stellate cells were purified from 18 C57BL/6 WT mice at 4-6 months old. Cells, combined from 1-4 mice per data point, were incubated for 20 minutes with either 10% FBS, GLP-2 (50 nM), or saline (Veh). Data are presented as the mean ± SEM. *P < 0.05, **P < 0.01, using ratio-paired t test.
liver (Supplemental Figure 4C). In a repository of data from hepatocyte single-cell RNA sequencing (12), these 3 markers were expressed at negligible levels in hepatocytes. Thus, we infer that stellate cell activation markers can be informative from analysis of whole liver RNA. Consistent with these findings, we detected increased transcript abundance for stellate cell markers Pdgfrb, Mmp13, Timp1, and Col8a1 after HFHC feeding, with higher levels for Timp1 and Col8a1 detected in livers from Glp2r -/mice compared with those from Glp2r +/+ mice ( Figure 6A). However, liver fibrosis, quantified by Sirius red staining ( Figure 6B), was increased in HFHC-fed mice but not different in Glp2r +/+ mice compared with Glp2r -/mice, and mRNA levels of the fibrogenic markers Col1a1 and Acta2 were comparable in Glp2r +/+ mice and Glp2r -/mice (Figure 6B). Furthermore, levels of mRNA transcripts for Col1a1 and Acta2, encoding markers of fibrosis, were also increased in hepatic RNA from 50-week-old Glp2r -/mice ( Figure 6C).
Activation or loss of GLP-2R signaling does not effect the hepatic response to CCl 4 . To assess the importance of the GLP-2R in an independent model of hepatic injury, we used the profibrotic agent CCl 4 . WT male C57/BL6 mice were treated with i.p. CCl 4 every 3 days, for a total of 4 doses (26). Starting 2 days prior to initiation of CCl 4 treatment and until termination 3 days after the last CCl 4 injection, mice were treated daily with GLP-2 or vehicle ( Figure 7A). CCl 4 -treated mice had elevated levels of plasma ALT, increased liver fibrosis (as approximated by Sirius red staining), and elevated levels of mRNA transcripts for several HSC activation markers analyzed in whole liver samples and in purified HSC preparations (Figure 7, B-E). Nevertheless, GLP-2 treatment had no effect on any of these endpoints or on body weight (Figure 7, B-E, and Supplemental Figure 4D). Unexpectedly, HSC Glp2r expression was approximately 8-fold lower in CCl 4 -treated mice ( Figure 7E). We next studied the consequences of CCl 4 administration in Glp2r +/+ and Glp2r -/male mice ( Figure 8A). Plasma ALT, the extent of Sirius red staining, and mRNA levels of HSC activation markers in the whole liver were generally similar between groups (Figure 8, B-E). Interestingly, mRNA levels of Lox, frequently induced in early stages of HSC activation (27), were upregulated in HSCs from Glp2r -/mice ( Figure 8E).

Discussion
The current studies were prompted in part by questions surrounding the use of the GLP-2 analog teduglutide in human subjects with intestinal failure at risk for development of PN-associated hepatosteatosis (28) and the conflicting data surrounding the putative actions of GLP-2 in the liver (13,17,29). GLP-2 rapidly increases postprandial levels of TGs in mice, hamsters, and humans (30)(31)(32), suggesting that increased substrate availability may contribute to hepatosteatosis in the setting of sustained GLP-2R signaling. Indeed, Taher and colleagues demonstrated increased hepatic lipogenesis in HFD-fed hamsters treated with GLP-2 as well as in rodent chow-fed (RC-fed) mice treated with GLP-2 (13). In contrast, treatment of HFDfed mice with the GLP-2R antagonist GLP-2(3-33) exacerbated hepatic steatosis in HFD-fed mice (17), whereas GLP-2 administration attenuated development of hepatosteatosis in PN-fed rats after small bowel resection (29). Our data, using administration of a degradation-resistant GLP-2R agonist, revealed no compelling effect of GLP-2R activation on hepatic steatosis in WT HFHC-fed mice. Collectively, our data are consistent with findings that teduglutide administration produced no change or reduced levels of transaminases and bilirubin in humans with SBS (33).
In contrast, we demonstrated that genetic loss of the Glp2r increased the extent of hepatic steatosis and upregulated hepatic cytokine expression in mice. These findings are consistent with data from Baldassano (17) and Cani (14), who demonstrated a protective role for GLP-2R signaling in HFD-fed mice and ob/ ob mice treated with or without prebiotics and GLP-2 or a GLP-2R antagonist, respectively. Intriguingly, recent evidence links HSC activity to de novo lipogenesis and the development of experimental steatosis (34). Hence, our current findings of hepatosteatosis and increased hepatic inflammation in Glp2r -/mice are consistent with a putative mechanism arising from loss of the HSC GLP-2R or potential consequences of absent GLP-2R signaling affecting control of gut barrier function, inflammation, and hepatic insulin sensitivity (14). The precise mechanisms linking loss of GLP-2R activity to increased hepatic fat accumulation remain uncertain and require additional investigation.
Despite multiple studies examining the roles of GLP-2R in liver biology, the precise cellular localization of the hepatic GLP-2R has remained uncertain. We localized the hepatic GLP-2R to stellate cells and not, as previously reported (11), in hepatocytes. Intriguingly, GLP-2 infusion rapidly enhances the appearance of circulating TG-rich lipoprotein-associated retinyl palmitate (35), establishing a putative link among GLP-2 action, vitamin A bioavailability, and subsequent storage in HSCs. Indeed, the HSC GLP-2R was functional, as GLP-2 directly induced a subset of IEGs in stellate cells ex vivo, findings consistent with previous studies of GLP-2R signaling ex vivo (36). Notably, relative Glp2r expression in purified HSCs approximates levels detected in RNA isolated from the jejunum, a major site of GLP-2R expression (37). Our data are also consistent with the detection of HSC Glp2r expression in studies analyzing single-cell RNA-sequencing data from murine stellate cells (38).
Interestingly, HSC Glp2r transcripts were upregulated after exposure to various diets used for experimental induction of hepatic injury, including a 45% HFD, a HFHC diet, and a MCD diet. These 3 different diets promote hepatic steatosis; however, the MCD diet induces a milder degree of steatosis, no weight gain, and increased hepatic inflammation and fibrosis (39). Conversely, mice with CCl 4 -induced liver injury exhibited a marked reduction in Glp2r expression within stellate cells. Collectively, these data show for the first time to our knowledge that Glp2r expression is dynamically regulated in the injured liver and dietary versus chemical hepatic insults produce divergent regulation of HSC Glp2r expression. These observations reveal the importance of interpreting hepatic GLP-2 action in the context of potential changes in expression of the endogenous GLP-2R.
Current concepts linking the gut to stellate cell activation highlight roles for microbial metabolites, including bile acids (40), as well as bacterial cell wall products, exemplified by lipopolysaccharide (41), which engage HSCs via distinct receptors and signaling pathways. Notably, bacterial metabolites (42) and cell membrane constituents such as LPS (43,44) as well as bile acids (45) are also potent stimulators of GLP-1 and GLP-2 secretion from enteroendocrine L cells (46). Our current findings demonstrate that loss of basal GLP-2R signaling is associated with HSC activation in mice with diet-induced hepatic injury. Collectively, these observations expand the concept of a gut-HSC axis through identification of a functional GLP-2R signaling system within HSCs, diversifying the mechanisms through which nutritional stimuli and microbiota may indirectly modulate HSC activity.
Our studies have several limitations. First, we used whole body Glp2r -/mice with germline loss of the Glp2r that may exhibit developmental or functional compensation, such as elevated levels of GLP-1 as shown here. Tissue-specific and conditional inactivation of the Glp2r selectively within HSCs  Figure 3A (n = 9-10 per group). (B) Representative Sirius red staining of liver sections with quantification of Sirius red-positive staining (left) and liver mRNA abundance, relative to Ppia, of fibrosis markers Col1a1and Acta2 (right) (n = 9-10 per group). Scale bar: 100 μm. (C) mRNA levels, relative to Ppia, of stellate cell activation markers in livers from Glp2r +/+ and Glp2r -/male mice age 50 ± 3.5 weeks old (n = 4 each). Data are presented as the mean ± SD. *P < 0.05, **P < 0.01, ***P < 0.001, ****P < 0.0001, using 2-way ANOVA with Tukey correction for multiple comparisons or Student's t test where appropriate.
would further establish the role of the HSC GLP-2R in the phenotypes observed in our studies. Moreover, we studied predominantly younger mice, and stellate cell activation is known to become more prominent with increasing age (47). Additionally, manipulation of central nervous system GLP-2R signaling has been shown to modulate hepatic insulin sensitivity (48), a pathway that was not directly examined in the current studies. In addition, we were unable to detect GLP-2R protein within HSCs due to limitations in the sensitivity and specificity of available GLP-2R antisera (7). Nevertheless, our experiments establish the presence of a full-length Glp2r mRNA and functional GLP-2R in HSCs, findings with potential importance for understanding how gut-derived signals effect HSC biology in a broader context.
Given the use of teduglutide in human children and adults with SBS and hepatosteatosis, taken together with the emerging importance of HSC activation in the pathophysiology of nonalcoholic fatty liver disease, our findings have translational relevance for understanding how the gut, microbiota, and microbial metabolites modify stellate cell adaptive responses to nutritional stress, inflammation, and hepatic injury.

Animal care and genotyping
Mice were housed (2-5 mice per cage), under a 12-hour-light/dark cycle in the Toronto Centre for Phenogenomics (TCP). Experiments were started in mice at 7-9 weeks of age or older if involving stellate cell isolation, as indicated in the figure legends. Age-matched mice were randomly assigned to experimental groups. All experiments were done in male mice.
WT C57Bl/6 male mice were purchased from The Jackson Laboratory, and Glp2r -/mice, generated in-house (49) on a C57Bl/6 background, were bred at the TCP animal facility. Genotyping was performed by PCR on tail DNA as previously described (49).
Rationale for choice of diet. Fructose-containing diets promote hepatic de novo lipogenesis and fat oxidation and modify the effects of a HFD on the development of insulin resistance (50). They are also associated with higher degrees of hepatic inflammation and fibrosis (51).
Additional cohorts. For data shown in Supplemental Figure 4, B and C, 12-week-old WT C57BL/6 male mice were maintained on a MCD or methionine-choline-sufficient diet for 14 days (Cedarlane, A02082002Bi or A02082003Bi, respectively). For data shown in Supplemental Figure 4B, 8-week-old WT mice were fed RC or 45% HFD (Research Diets, D12451i) for 12 weeks. For data shown in Figure 6C, age-matched Glp2r +/+ and Glp2r -/mice were bred and housed in the TCP animal facility and fed a regular chow diet until age 1 year.
Excluded mice. One mouse from the experiment described in Figures 1 and 2 on CD was excluded prior to allocation to GLP-2 treatment due to unexplained extremely poor weight gain. One mouse from the experiment described in Figures 3 and 4 from the Glp2r +/+ group fed HFHC was excluded from analyses after repeat validation genotyping after sacrificing revealed Glp2r heterozygosity.

GLP-2 treatment
Mice were injected subcutaneously with 0.2 mg/Kg/d of the long-acting GLP-2 analog h[Gly2]GLP-2. Custom synthesized human [Gly2]GLP-2, used for the experiment described in Figures 1 and 2, was from Pepceuticals Ltd. Human [Gly2]GLP-2 (Teduglutide), used for the experiment described in Figure 7 and Supplemental Figure 4D, was a gift from Shire- Takeda Inc. GLP-2 was administered twice daily at 9 am and 5 pm for 11 days (for the experiment described in Figures 1 and 2) or once daily at 9-10 am for 14 days (for the experiment described in Figure 7). Control mice were treated with saline in the same volume. GLP-2 doses and length of treatment were based on previous studies describing the intestinotrophic effects.

Blood and tissue collection
Blood was taken via the tail vein during metabolic studies. For isolation of liver fractions, mice were anesthetized with Avertin. For other terminal studies, mice were sacrificed by CO 2 inhalation, and blood was collected via cardiac puncture into 1-ml syringes washed with heparin (10 units/ml) and centrifuged at 14,000 g for 10 minutes at 4°C to collect the plasma. Plasma samples were stored at -80°C until further analyses.
Intestines were flushed with PBS and omental fat was removed prior to weight and length measurements. Sections were obtained from the proximal jejunum for RNA and protein analyses. Gallbladder volume measurement was performed as described previously (8). Briefly, the gallbladder was ligated to prevent emptying and resected prior to liver manipulation. A photograph of the gallbladder was then taken next to a ruler, and volume was calculated as ellipsoid volume, assuming both short diameters to be equal. The liver was resected and weighed after removal of gallbladder.

Glucose tolerance test
After a 5-to 6-hour fast, glucose (1.5 mg/g body wt; MilliporeSigma) was administered through i.p. injection. A blood sample was drawn from a tail vein at 0, 10, 20, 30, 60, 90, and 120 minutes after glucose administration, and blood glucose levels were measured using a Glucometer Elite blood glucometer (Bayer).

Lipid tolerance test
RC diet-fed Glp2r +/+ and Glp2r -/littermate mice (8-12 weeks old) were gavaged orally with 200 μl olive oil (MilliporeSigma) following a 5-to 6-hour fast. Tail vein blood samples were collected into heparin-coated capillary tubes (a) before gavage only for DPP-4 activity; (b) before and 1, 2, and 3 hours after gavage for measurement of TG; and (c) before and 10 minutes after gavage for measurement of active GLP-1. Plasma was prepared from the collected blood samples as described above.
All other analytes in plasma were performed by the pathology core at The Centre for Phenogenomics. The Beckman Coulter AU480 clinical chemistry analyzer, which measures analytes in plasma or serum samples using photometry testing, was used in combination with appropriate reagents (ALT, AST, TG, cholesterol, LDL, and HDL), calibrators (Beckman Coulter Lyophilized Chemistry Calibrator Levels 1 and 2), and quality control materials (Bio-Rad Liquid Assayed Multiqual 1 and 3).

Measurement of plasma and hepatic cytokine levels
Cytokines were measured in 25 μl mouse plasma or in 40-70 mg liver tissue homogenized in 500 μl ice-cold lysis buffer (50 mM Tris HCl, pH 8, 1 mM EDTA, 10% glycerol [wt/vol], 0.02% Brij-35 [wt/vol]). Lysis buffer was supplemented with 1 mM dithiothreitol and protease inhibitor for cytokine assay only. Tissue and plasma (containing 10% TED) levels of cytokines were measured using the MSD V-PLEX Mouse Proinflammatory Panel assay kit (Meso Scale). Cytokine concentrations in tissue were normalized to total protein amount measured by Bradford Assay.

Liver histology
Liver samples were fixed in 10% buffered formalin for 48 hours and then transferred into 70% ethanol until processing. For Oil Red O staining, liver tissue was flash frozen in OCT compound. Samples were paraffin embedded, sectioned at 4-μm thickness, slide mounted, and stained with H&E, Picrosirius red stain, or Oil Red O using standard protocols. Sections were scanned using the ScanScope CS System (Aperio Technologies), and digital images were captured at ×10 magnification using ImageScope Software (Aperio Technologies) or an OlyVIA 2.9 system (OLYMPUS). Immunostaining was performed by the pathology core at The Centre for Phenogenomics.
Tissue fibrosis was analyzed in Sirius red-stained samples by red pixel count corresponding to stained collagen using the Aperio Imagescope Positive pixel count program that was calibrated once and applied uniformly to 5 ×10 images from each sample. The percentage of red pixels in total pixel area was determined. Tissue lipids were quantified in a similar manner using Oil Red O-stained sections. Three to four images per sample were analyzed using ImageJ software (NIH) reporting red and total pixel counts.

RNA isolation, cDNA, and qPCR
Tissue samples were dissected, weighed, frozen in liquid nitrogen, and stored at −80°C. Frozen tissues were homogenized in Tri Reagent (Molecular Research Center) using a TissueLyser II system (Qiagen), and total RNA was extracted using the manufacturer's protocol. Reverse transcription was performed with RNA treated with DNase I (Thermo Fisher Scientific, EN0521) using random hexamers (Thermo Fisher Scientific, 58875) and SuperScript III (Thermo Fisher Scientific, 18080-044). The resultant cDNA was used to assess mRNA expression by real-time quantitative PCR (QuantStudio 5 System, Thermo Fisher Scientific) with TaqMan Fast Advanced Master Mix (Thermo Fisher Scientific, 4444557) and TaqMan Gene Expression Assays (Thermo Fisher Scientific). The specific gene expression assays used are listed in Supplemental Table 1. Quantification of transcript levels was performed by the 2 -ΔCtT method, and expression levels for each gene were normalized to Ppia or Actb (β-actin). RNA extracted from purified HSCs was normalized to Ppia or Pdgfrb.

Glp2r full-length mRNA studies
For human studies, the expression of full-length GLP2R and GAPDH transcripts was assessed by RT-PCR in total RNA from the human colon cancer cell line DLD-1, as previously confirmed (54), to stably express the human GLP-2R and in human liver total RNA purchased from the specified commercial sources (Lifetech, Agilent, Clontech). PCR products were analyzed by agarose gel electrophoresis followed by SYBR green staining and by Southern blot analysis using an internal GLP2R-specific 32 P-labeled oligonucleotide probe. The presence of genomic DNA contamination in samples was monitored, in parallel, with control reactions in which the reverse transcriptase (RT) enzyme was omitted from the reaction mixture (RT-). Expression of mouse Glp2r and Gapdh transcripts was assessed following a similar protocol in mouse jejunum (positive control), liver, isolated hepatocytes, and purified HSCs. The human GLP2R PCR primer set is as follows: forward 5′-TTGTGAAGGTGCACGAGGAA-3′ and reverse 5′-CACCCTAGATCT-CACTCTCT-3′. Primers were amplified for 35 cycles at an annealing temperature of 60°C. The product size of 1,687 bp spanned the entire GLP-2R open reading frame. The Gapdh PCR primer set is as follows: forward 5′-GACCACAGTCCATGACATCACT-3′ and reverse 5′-TCCACCACCCTGTTGCTG-TAG-3′. Primers were amplified for 25 cycles at an annealing temperature of 60°C. The product size was 450 bp. This primer set is identical in human and mouse. The mouse Glp2r primer set is as follows: forward 5′-GCCCAGTAGATGCAGAGAGG-3′ and reverse 5′-AGTTGCCAAGCTGTGGTGAT-3′. Amplification was set for 37 cycles at an annealing temperature of 60°C. The product size of 1,661 bp spanned the entire GLP-2R open reading frame.

Hepatocyte and hepatic NPC isolation
Primary hepatocyte isolation was as previously described (55). For purification of the hepatic NPC fraction, the hepatocytes were depleted by 2-minute centrifugation at 50 g, and the supernatant containing NPCs was isolated and centrifuged for 5 minutes at 300 g. The pellet was resuspended in hepatocyte wash medium (Thermo Fisher Scientific, 17704024). After 4 cycles of hepatocyte depletion, the NPC pellet was enriched in cells expressing macrophage and endothelial cell markers and depleted of albumin-expressing hepatocytes (Supplemental Figure 3A) and was collected for further analysis.

Stellate cell isolation and FACS purification
Isolation of stellate cells from mouse livers was performed following the protocol described in ref. 20. Briefly, the inferior vena cava was cannulated, the portal vein was resected, and the superior vena cava was clamped. Where appropriate, a segment of the left liver lobe was clamped prior to cannulation and resected for histology and mRNA/protein expression analyses. Following cannulation, the liver was perfused with a series of buffers containing EGTA, pronase, and collagenase; removed; and mechanically disrupted in a pronase plus collagenase mix at 37°C. Stellate cells were then isolated by Nycodenz-based (Alere Technologies, 1002424) gradient isolation at 1380 g at 4°C.
The stellate enriched fraction was removed, stained with the vital dye 7-aminoactinomycin D (7-AAD) (BD Pharmingen, 559925), and purified by flow cytometry performed by the flow cytometry core at The Centre for Phenogenomics using MoFlo Astrios EQ (Beckman Coulter) equipped with 405 nm, 488 nm, 561 nm, and 640 nm lasers and running Summit v.6.2 software. Live stellate cells were identified by the following fluorescent signal combination: low 7-AAD (detector parameter 488-620/29 nm), high 405 nm violet laser autofluorescence (detector parameter 405-448/59 nm) and low 488 nm blue laser autofluorescence (detector parameter 488-513/26 nm). FACS gates were based on stellate cell autofluorescence in the violet range owing to cell-type specific vitamin A laden lipid droplets. Additional gating to remove cells based on autofluorescence in the FITC range enhanced purification from Kupffer cells and other minor contaminants. The postsorting violet + FITCgate was confirmed to be highly enriched for HSC markers and depleted for markers of other hepatic cell populations ( Figure 5D). Postsorting HSC viability was sampled frequently by 7-AAD and was reproducibly above 95%.

Flow cytometry for fractionation of hepatic NPC populations
After collection of NPC-enriched fractions, cells were passed through a 40-μm filter, collected by centrifugation at 300 g, and resuspended in hepatocyte wash medium (Invitrogen-Gibco). Following 4 cycles of hepatocyte depletion, the NPC pellet was resuspended in FACS-sorting buffer (0.5 M EDTA, 25 mM HEPES, and 1% FCS in PBS); incubated for 30 minutes at 4°C in antibody mixes based on CD-26-PE (Biolegend, 137804), CD-45-APC (Biolegend, 109814), and CD-31-FITC (Biolegend, 102406); and then washed and sorted. Flow cytometry was performed by the flow cytometry core at The Centre for Phenogenomics using equipment and software as detailed above.
In vitro transcriptional response of IEGs in stellate cells IEG expression is activated by the stellate cell purification process (56) as well as by FBS that is typically added during cell sorting and recovery. Hence, we purified HSCs using standard protocols (20), and, following purification, pelleted the cells at 580 g and resuspended them in 1.5 ml DMEM plus charcoal Dextran-treated (MilliporeSigma, C6241) bovine serum albumin (MilliporeSigma, A7030) BSA 1% for 50 minutes at 37°C. Following incubation, samples were divided into 3 equal parts and 10% FBS (positive control), GLP-2 (50 nM), or vehicle (PBS, in the same volume as GLP-2) was added. The dose of GLP-2 was based on previously established doses for direct induction of GLP-2R signaling in cells ex vivo (36).
To validate that incubation time sufficiently recovered IEG expression to levels amenable to activation, we tested whether the positive control cells, treated with 10% FBS demonstrated expression at least 5%