Constitutively Active STAT5b Feminizes Mouse Liver Gene Expression

HDI-STARR-seq: Condition-specific enhancer discovery in mouse liver in vivo

Background

STARR-seq and other massively-parallel reporter assays are widely used to discover functional enhancers in transfected cell models, which can be confounded by plasmid vector-induced type-I interferon immune responses and lack the multicellular environment and endogenous chromatin state of complex mammalian tissues.

Results

Here, we describe HDI-STARR-seq, which combines STARR-seq plasmid library delivery to the liver, by hydrodynamic tail vein injection (HDI), with reporter RNA transcriptional initiation driven by a minimal Albumin promoter, which we show is essential for mouse liver STARR-seq enhancer activity assayed 7 days after HDI. Importantly, little or no vector-induced innate type-I interferon responses were observed. Comparisons of HDI-STARR-seq activity between male and female mouse livers and in livers from males treated with an activating ligand of the transcription factor CAR (Nr1i3) identified many condition-dependent enhancers linked to condition-specific gene expression. Further, thousands of active liver enhancers were identified using a high complexity STARR-seq library comprised of ~ 50,000 genomic regions released by DNase-I digestion of mouse liver nuclei. When compared to stringently inactive library sequences, the active enhancer sequences identified were highly enriched for liver open chromatin regions with activating histone marks (H3K27ac, H3K4me1, H3K4me3), were significantly closer to gene transcriptional start sites, and were significantly depleted of repressive (H3K27me3, H3K9me3) and transcribed region histone marks (H3K36me3).

Conclusions

HDI-STARR-seq offers substantial improvements over current methodologies for large scale, functional profiling of enhancers, including condition-dependent enhancers, in liver tissue in vivo, and can be adapted to characterize enhancer activities in a variety of species and tissues by selecting suitable tissue- and species-specific promoter sequences.

HDI-STARR-seq: Condition-specific enhancer discovery in mouse liver in vivo

STARR-seq and other massively-parallel reporter assays are widely used to discover functional enhancers in transfected cell models, which can be confounded by plasmid vector-induced type-I interferon immune responses and lack the multicellular environment and endogenous chromatin state of complex mammalian tissues. Here, we describe HDI-STARR-seq, which combines STARR-seq plasmid library delivery to the liver, by hydrodynamic tail vein injection (HDI), with reporter RNA transcriptional initiation driven by a minimal Albumin promoter, which we show is essential for mouse liver STARR-seq enhancer activity assayed 7 days after HDI. Importantly, little or no vector-induced innate type-I interferon responses were observed. Comparisons of HDI-STARR-seq activity between male and female mouse livers and in livers from males treated with an activating ligand of the transcription factor CAR (Nr1i3) identified many condition-dependent enhancers linked to condition-specific gene expression. Further, thousands of active liver enhancers were identified using a high complexity STARR-seq library comprised of ~50,000 genomic regions released by DNase-I digestion of mouse liver nuclei. When compared to stringently inactive library sequences, the active enhancer sequences identified were highly enriched for liver open chromatin regions with activating histone marks (H3K27ac, H3K4me1, H3K4me3), were significantly closer to gene transcriptional start sites, and were significantly depleted of repressive (H3K27me3, H3K9me3) and transcribed region histone marks (H3K36me3). HDI-STARR-seq offers substantial improvements over current methodologies for large scale, functional profiling of enhancers, including condition-dependent enhancers, in liver tissue in vivo, and can be adapted to characterize enhancer activities in a variety of species and tissues by selecting suitable tissue- and species-specific promoter sequences.

Modifications of the GH Axis Reveal Unique Sexually Dimorphic Liver Signatures for Lcn13, Asns, Hamp2, Hao2, and Pgc1a

Growth hormone (GH) modifies liver gene transcription in a sexually dimorphic manner to meet liver metabolic demands related to sex; thus, GH dysregulation leads to sex-biased hepatic disease. We dissected the steps of the GH regulatory cascade modifying GH-dependent genes involved in metabolism, focusing on the male-predominant genes Lcn13, Asns, and Cyp7b1, and the female-predominant genes Hao2, Pgc1a, Hamp2, Cyp2a4, and Cyp2b9. We explored mRNA expression in 2 settings: (i) intact liver GH receptor (GHR) but altered GH and insulin-like growth factor 1 (IGF1) levels (NeuroDrd2KO, HiGH, aHepIGF1kd, and STAT5bCA mouse lines); and (ii) liver loss of GHR, with or without STAT5b reconstitution (aHepGHRkd, and aHepGHRkd + STAT5bCA). Lcn13 was downregulated in males in most models, while Asns and Cyp7b1 were decreased in males by low GH levels or action, or constant GH levels, but unexpectedly upregulated in both sexes by the loss of liver Igf1 or constitutive Stat5b expression. Hao, Cyp2a4, and Cyp2b9 were generally decreased in female mice with low GH levels or action (NeuroDrd2KO and/or aHepGHRkd mice) and increased in HiGH females, while in contrast, Pgc1a was increased in female NeuroDrd2KO but decreased in STAT5bCA and aHepIGF1kd females. Bioinformatic analysis of RNAseq from aHepGHRkd livers stressed the greater impact of GHR loss on wide gene expression in males and highlighted that GH modifies almost completely different gene signatures in each sex. Concordantly, we show that altering different steps of the GH cascade in the liver modified liver expression of Lcn13, Asns, Cyp7b1, Hao2, Hamp2, Pgc1a, Cyp2a4, and Cyp2b9 in a sex- and gene-specific manner.

Plasma growth hormone pulses induce male-biased pulsatile chromatin opening and epigenetic regulation in adult mouse liver

Sex differences in plasma growth hormone (GH) profiles, pulsatile in males and persistent in females, regulate sex differences in hepatic STAT5 activation linked to sex differences in gene expression and liver disease susceptibility, but little is understood about the fundamental underlying, GH pattern-dependent regulatory mechanisms. Here, DNase-I hypersensitivity site (DHS) analysis of liver chromatin accessibility in a cohort of 18 individual male mice established that the endogenous male rhythm of plasma GH pulse-stimulated liver STAT5 activation induces dynamic, repeated cycles of chromatin opening and closing at several thousand liver DHS and comprises a novel mechanism conferring male bias to liver chromatin accessibility. Strikingly, a single physiological replacement dose of GH given to hypophysectomized male mice restored, within 30 min, liver STAT5 activity and chromatin accessibility at 83% of the dynamic, pituitary hormone-dependent male-biased DHS. Sex-dependent transcription factor binding patterns and chromatin state analysis identified key genomic and epigenetic features distinguishing this dynamic, STAT5-driven mechanism of male-biased chromatin opening from a second GH-dependent mechanism operative at static male-biased DHS, which are constitutively open in male liver. Dynamic but not static male-biased DHS adopt a bivalent-like epigenetic state in female liver, as do static female-biased DHS in male liver, albeit using distinct repressive histone marks in each sex, namely, H3K9me3 at male-biased DHS in female liver and H3K27me3 at female-biased DHS in male liver. Moreover, sex-biased H3K36me3 marks are uniquely enriched at static sex-biased DHS, which may serve to keep these sex-dependent hepatocyte enhancers free of H3K27me3 repressive marks and thus constitutively open. Pulsatile chromatin opening stimulated by endogenous, physiological hormone pulses is thus one of two distinct GH-determined mechanisms for establishing widespread sex differences in hepatic chromatin accessibility and epigenetic regulation, both closely linked to sex-biased gene transcription and the sexual dimorphism of liver function.

Plasma Growth Hormone Pulses Induce Male-biased Pulsatile Chromatin Opening and Epigenetic Regulation in Adult Mouse Liver

Sex-differences in plasma growth hormone (GH) profiles, pulsatile in males and persistent in females, regulate sex differences in hepatic STAT5 activation linked to sex differences in gene expression and liver disease susceptibility, but little is understood about the fundamental underlying, GH pattern-dependent regulatory mechanisms. Here, DNase hypersensitivity site (DHS) analysis of liver chromatin accessibility in a cohort of 18 individual male mice established that the endogenous male rhythm of plasma GH pulse-stimulated liver STAT5 activation induces dynamic, repeated cycles of chromatin opening and closing at several thousand liver DHS and comprises a novel mechanism conferring male bias to liver chromatin accessibility. Strikingly, a single physiological replacement dose of GH given to hypophysectomized male mice restored, within 30 min, liver STAT5 activity and chromatin accessibility at 83% of the pituitary hormone-dependent dynamic male-biased DHS. Sex-dependent transcription factor binding patterns and chromatin state analysis identified key genomic and epigenetic features distinguishing this dynamic, STAT5-driven mechanism of male-biased chromatin opening from a second GH-dependent mechanism operative at static male-biased DHS, which are constitutively open in male liver. Dynamic but not static male-biased DHS adopt a bivalent-like epigenetic state in female liver, as do static female-biased DHS in male liver, albeit using distinct repressive histone marks in each sex, namely, H3K27me3 at female-biased DHS in male liver, and H3K9me3 at male-biased DHS in female liver. Moreover, sex-biased H3K36me3 marks are uniquely enriched at static sex-biased DHS, which may serve to keep these sex-dependent hepatocyte enhancers free of H3K27me3 repressive marks and thus constitutively open. Pulsatile chromatin opening stimulated by endogenous, physiological hormone pulses is thus one of two distinct GH-determined mechanisms for establishing widespread sex differences in hepatic chromatin accessibility and epigenetic regulation, both closely linked to sex-biased gene transcription and the sexual dimorphism of liver function.

Direct androgen receptor control of sexually dimorphic gene expression in the mammalian kidney

Mammalian organs exhibit distinct physiology, disease susceptibility, and injury responses between the sexes. In the mouse kidney, sexually dimorphic gene activity maps predominantly to proximal tubule (PT) segments. Bulk RNA sequencing (RNA-seq) data demonstrated that sex differences were established from 4 and 8 weeks after birth under gonadal control. Hormone injection studies and genetic removal of androgen and estrogen receptors demonstrated androgen receptor (AR)-mediated regulation of gene activity in PT cells as the regulatory mechanism. Interestingly, caloric restriction feminizes the male kidney. Single-nuclear multiomic analysis identified putative cis-regulatory regions and cooperating factors mediating PT responses to AR activity in the mouse kidney. In the human kidney, a limited set of genes showed conserved sex-linked regulation, whereas analysis of the mouse liver underscored organ-specific differences in the regulation of sexually dimorphic gene expression. These findings raise interesting questions on the evolution, physiological significance, disease, and metabolic linkage of sexually dimorphic gene activity.

RISING STARS: Sex differences in toxicant-associated fatty liver disease

Based on biological sex, the consequential health outcomes from exposures to environmental chemicals or toxicants can differ in disease pathophysiology, progression, and severity. Due to basal differences in cellular and molecular processes resulting from sexual dimorphism of organs including the liver and additional factors influencing 'gene-environment' interactions, males and females can exhibit different responses to toxicant exposures. Associations between environmental/occupational chemical exposures and fatty liver disease (FLD) have been well-acknowledged in human epidemiologic studies and their causal relationships demonstrated in experimental models. However, studies related to sex differences in liver toxicology are still limited to draw any inferences on sex-dependent chemical toxicity. The purpose of this review is to highlight the present state of knowledge on the existence of sex differences in toxicant-associated FLD (TAFLD), discuss potential underlying mechanisms driving these differences, implications of said differences on disease susceptibility, and emerging concepts. Chemicals of interest include various categories of pollutants that have been investigated in TAFLD, namely persistent organic pollutants, volatile organic compounds, and metals. Insight into research areas requiring further development is also discussed, with the objective of narrowing the knowledge gap on sex differences in environmental liver diseases. Major conclusions from this review exercise are that biological sex influences TAFLD risks, in part due to (i) toxicant disruption of growth hormone and estrogen receptor signaling, (ii) basal sex differences in energy mobilization and storage, and (iii) differences in chemical metabolism and subsequent body burden. Finally, further sex-dependent toxicological assessments are warranted for the development of sex-specific intervention strategies.

Direct and systemic actions of growth hormone receptor (GHR)-signaling on hepatic glycolysis, de novo lipogenesis and insulin sensitivity, associated with steatosis

BACKGROUND

Evidence is accumulating that growth hormone (GH) protects against the development of steatosis and progression of non-alcoholic fatty liver disease (NAFLD). GH may control steatosis indirectly by altering systemic insulin sensitivity and substrate delivery to the liver and/or by the direct actions of GH on hepatocyte function.

APPROACH

To better define the hepatocyte-specific role of GH receptor (GHR) signaling on regulating steatosis, we used a mouse model with adult-onset, hepatocyte-specific GHR knockdown (aHepGHRkd). To prevent the reduction in circulating insulin-like growth factor 1 (IGF1) and the subsequent increase in GH observed after aHepGHRkd, subsets of aHepGHRkd mice were treated with adeno-associated viral vectors (AAV) driving hepatocyte-specific expression of IGF1 or a constitutively active form of STAT5b (STAT5bCA). The impact of hepatocyte-specific modulation of GHR, IGF1 and STAT5b on carbohydrate and lipid metabolism was studied across multiple nutritional states and in the context of hyperinsulinemic:euglycemic clamps.

RESULTS

Chow-fed male aHepGHRkd mice developed steatosis associated with an increase in hepatic glucokinase (GCK) and ketohexokinase (KHK) expression and de novo lipogenesis (DNL) rate, in the post-absorptive state and in response to refeeding after an overnight fast. The aHepGHRkd-associated increase in hepatic KHK, but not GCK and steatosis, was dependent on hepatocyte expression of carbohydrate response element binding protein (ChREBP), in re-fed mice. Interestingly, under clamp conditions, aHepGHRkd also increased the rate of DNL and expression of GCK and KHK, but impaired insulin-mediated suppression of hepatic glucose production, without altering plasma NEFA levels. These effects were normalized with AAV-mediated hepatocyte expression of IGF1 or STAT5bCA. Comparison of the impact of AAV-mediated hepatocyte IGF1 versus STAT5bCA in aHepGHRkd mice across multiple nutritional states, indicated the restorative actions of IGF1 are indirect, by improving systemic insulin sensitivity, independent of changes in the liver transcriptome. In contrast, the actions of STAT5b are due to the combined effects of raising IGF1 and direct alterations in the hepatocyte gene program that may involve suppression of BCL6 and FOXO1 activity. However, the direct and IGF1-dependent actions of STAT5b cannot fully account for enhanced GCK activity and lipogenic gene expression observed after aHepGHRkd, suggesting other GHR-mediated signals are involved.

CONCLUSION

These studies demonstrate hepatocyte GHR-signaling controls hepatic glycolysis, DNL, steatosis and hepatic insulin sensitivity indirectly (via IGF1) and directly (via STAT5b). The relative contribution of these indirect and direct actions of GH on hepatocytes is modified by insulin and nutrient availability. These results improve our understanding of the physiologic actions of GH on regulating adult metabolism to protect against NAFLD progression.

Dysregulation of murine long noncoding single-cell transcriptome in nonalcoholic steatohepatitis and liver fibrosis

LncRNAs comprise a heterogeneous class of RNA-encoding genes typified by low expression, nuclear enrichment, high tissue-specificity, and functional diversity, but the vast majority remain uncharacterized. Here, we assembled the mouse liver noncoding transcriptome from >2000 bulk RNA-seq samples and discovered 48,261 liver-expressed lncRNAs, a majority novel. Using these lncRNAs as a single-cell transcriptomic reference set, we elucidated lncRNA dysregulation in mouse models of high fat diet-induced nonalcoholic steatohepatitis and carbon tetrachloride-induced liver fibrosis. Trajectory inference analysis revealed lncRNA zonation patterns across the liver lobule in each major liver cell population. Perturbations in lncRNA expression and zonation were common in several disease-associated liver cell types, including nonalcoholic steatohepatitis-associated macrophages, a hallmark of fatty liver disease progression, and collagen-producing myofibroblasts, a central feature of liver fibrosis. Single-cell-based gene regulatory network analysis using bigSCale2 linked individual lncRNAs to specific biological pathways, and network-essential regulatory lncRNAs with disease-associated functions were identified by their high network centrality metrics. For a subset of these lncRNAs, promoter sequences of the network-defined lncRNA target genes were significantly enriched for lncRNA triplex formation, providing independent mechanistic support for the lncRNA-target gene linkages predicted by the gene regulatory networks. These findings elucidate liver lncRNA cell-type specificities, spatial zonation patterns, associated regulatory networks, and temporal patterns of dysregulation during hepatic disease progression. A subset of the liver disease-associated regulatory lncRNAs identified have human orthologs and are promising candidates for biomarkers and therapeutic targets.

Direct androgen receptor regulation of sexually dimorphic gene expression in the mammalian kidney

Mammalian organs exhibit distinct physiology, disease susceptibility and injury responses between the sexes. In the mouse kidney, sexually dimorphic gene activity maps predominantly to proximal tubule (PT) segments. Bulk RNA-seq data demonstrated sex differences were established from 4 and 8 weeks after birth under gonadal control. Hormone injection studies and genetic removal of androgen and estrogen receptors demonstrated androgen receptor (AR) mediated regulation of gene activity in PT cells as the regulatory mechanism. Interestingly, caloric restriction feminizes the male kidney. Single-nuclear multiomic analysis identified putative cis-regulatory regions and cooperating factors mediating PT responses to AR activity in the mouse kidney. In the human kidney, a limited set of genes showed conserved sex-linked regulation while analysis of the mouse liver underscored organ-specific differences in the regulation of sexually dimorphic gene expression. These findings raise interesting questions on the evolution, physiological significance, and disease and metabolic linkage, of sexually dimorphic gene activity.

Growth Hormone Accelerates Recovery From Acetaminophen-Induced Murine Liver Injury

Background and Aims

Acetaminophen (APAP) overdose is the leading cause of acute liver failure, with one available treatment, N-acetyl cysteine (NAC). Yet, NAC effectiveness diminishes about ten hours after APAP overdose, urging for therapeutic alternatives. This study addresses this need by deciphering a mechanism of sexual dimorphism in APAP-induced liver injury, and leveraging it to accelerate liver recovery via growth hormone (GH) treatment. GH secretory patterns, pulsatile in males and near-continuous in females, determine the sex bias in many liver metabolic functions. Here, we aim to establish GH as a novel therapy to treat APAP hepatotoxicity.

Approach and Results

Our results demonstrate sex-dependent APAP toxicity, with females showing reduced liver cell death and faster recovery than males. Single-cell RNA sequencing analyses reveal that female hepatocytes have significantly greater levels of GH receptor expression and GH pathway activation compared to males. In harnessing this female-specific advantage, we demonstrate that a single injection of recombinant human GH protein accelerates liver recovery, promotes survival in males following sub-lethal dose of APAP, and is superior to standard-of-care NAC. Alternatively, slow-release delivery of human GH via the safe nonintegrative lipid nanoparticle-encapsulated nucleoside-modified mRNA (mRNA-LNP), a technology validated by widely used COVID-19 vaccines, rescues males from APAP-induced death that otherwise occurred in control mRNA-LNP-treated mice.

Conclusions

Our study demonstrates a sexually dimorphic liver repair advantage in females following APAP overdose, leveraged by establishing GH as an alternative treatment, delivered either as recombinant protein or mRNA-LNP, to potentially prevent liver failure and liver transplant in APAP-overdosed patients.

Sex matters in liver fat regulation

Growth hormone orchestrates a complex, sex-dependent balancing act.

Interplay Between GH-regulated, Sex-biased Liver Transcriptome and Hepatic Zonation Revealed by Single-Nucleus RNA Sequencing

The zonation of liver metabolic processes is well-characterized; however, little is known about the cell type-specificity and zonation of sexually dimorphic gene expression or its growth hormone (GH)-dependent transcriptional regulators. We address these issues using single-nucleus RNA-sequencing of 32 000 nuclei representing 9 major liver cell types. Nuclei were extracted from livers from adult male and female mice; from males infused with GH continuously, mimicking the female plasma GH pattern; and from mice exposed to TCPOBOP, a xenobiotic agonist ligand of the nuclear receptor CAR that perturbs sex-biased gene expression. Analysis of these rich transcriptomic datasets revealed the following: 1) expression of sex-biased genes and their GH-dependent transcriptional regulators is primarily restricted to hepatocytes and is not a feature of liver nonparenchymal cells; 2) many sex-biased transcripts show sex-dependent zonation within the liver lobule; 3) gene expression is substantially feminized both in periportal and pericentral hepatocytes when male mice are infused with GH continuously; 4) sequencing nuclei increases the sensitivity for detecting thousands of nuclear-enriched long-noncoding RNAs (lncRNAs) and enables determination of their liver cell type-specificity, sex-bias and hepatocyte zonation profiles; 5) the periportal to pericentral hepatocyte cell ratio is significantly higher in male than female liver; and 6) TCPOBOP exposure disrupts both sex-specific gene expression and hepatocyte zonation within the liver lobule. These findings highlight the complex interconnections between hepatic sexual dimorphism and zonation at the single-cell level and reveal how endogenous hormones and foreign chemical exposure can alter these interactions across the liver lobule with large effects both on protein-coding genes and lncRNAs.