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.

MAnorm2 for quantitatively comparing groups of ChIP-seq samples.

Eukaryotic gene transcription is regulated by a large cohort of chromatin-associated proteins, and inferring their differential binding sites between cellular contexts requires a rigorous comparison of the corresponding ChIP-seq data. We present MAnorm2, a new computational tool for quantitatively comparing groups of ChIP-seq samples. MAnorm2 uses a hierarchical strategy for normalization of ChIP-seq data and assesses within-group variability of ChIP-seq signals based on an empirical Bayes framework. In this framework, MAnorm2 allows for abundant differential ChIP-seq signals between groups of samples as well as very different global within-group variability between groups. Using a number of real ChIP-seq data sets, we observed that MAnorm2 clearly outperformed existing tools for differential ChIP-seq analysis, especially when the groups of samples being compared had distinct global within-group variability.

Harnessing natural variation to identify cis regulators of sex-biased gene expression in a multi-strain mouse liver model.

Sex differences in gene expression are widespread in the liver, where many autosomal factors act in tandem with growth hormone signaling to regulate individual variability of sex differences in liver metabolism and disease. Here, we compare hepatic transcriptomic and epigenetic profiles of mouse strains C57BL/6J and CAST/EiJ, representing two subspecies separated by 0.5-1 million years of evolution, to elucidate the actions of genetic factors regulating liver sex differences. We identify 144 protein coding genes and 78 lncRNAs showing strain-conserved sex bias; many have gene ontologies relevant to liver function, are more highly liver-specific and show greater sex bias, and are more proximally regulated than genes whose sex bias is strain-dependent. The strain-conserved genes include key growth hormone-dependent transcriptional regulators of liver sex bias; however, three other transcription factors, Trim24, Tox, and Zfp809, lose their sex-biased expression in CAST/EiJ mouse liver. To elucidate the observed strain specificities in expression, we characterized the strain-dependence of sex-biased chromatin opening and enhancer marks at cis regulatory elements (CREs) within expression quantitative trait loci (eQTL) regulating liver sex-biased genes. Strikingly, 208 of 286 eQTLs with strain-specific, sex-differential effects on expression were associated with a complete gain, loss, or reversal of the sex differences in expression between strains. Moreover, 166 of the 286 eQTLs were linked to the strain-dependent gain or loss of localized sex-biased CREs. Remarkably, a subset of these CREs apparently lacked strain-specific genetic variants yet showed coordinated, strain-dependent sex-biased epigenetic regulation. Thus, we directly link hundreds of strain-specific genetic variants to the high variability in CRE activity and expression of sex-biased genes and uncover underlying genetically-determined epigenetic states controlling liver sex bias in genetically diverse mouse populations.

TCDD dysregulation of lncRNA expression, liver zonation and intercellular communication across the liver lobule.

The persistent environmental aryl hydrocarbon receptor agonist and hepatotoxin TCDD (2,3,7,8-tetrachlorodibenzo-p-dioxin) induces hepatic lipid accumulation (steatosis), inflammation (steatohepatitis) and fibrosis. Thousands of liver-expressed, nuclear-localized lncRNAs with regulatory potential have been identified; however, their roles in TCDD-induced hepatoxicity and liver disease are unknown. We analyzed single nucleus (sn)RNA-seq data from control and subchronic (4 wk) TCDD-exposed mouse liver to determine liver cell-type specificity, zonation and differential expression profiles for thousands of lncRNAs. TCDD dysregulated >4000 of these lncRNAs in one or more liver cell types, including 684 lncRNAs specifically dysregulated in liver non-parenchymal cells. Trajectory inference analysis revealed major disruption by TCDD of hepatocyte zonation, affecting >800 genes, including 121 lncRNAs, with strong enrichment for lipid metabolism genes. TCDD also dysregulated expression of >200 transcription factors, including 19 Nuclear Receptors, most notably in hepatocytes and Kupffer cells. TCDD-induced changes in cell-cell communication patterns included marked decreases in EGF signaling from hepatocytes to non-parenchymal cells and increases in extracellular matrix-receptor interactions central to liver fibrosis. Gene regulatory networks constructed from the snRNA-seq data identified TCDD-exposed liver network-essential lncRNA regulators linked to functions such as fatty acid metabolic process, peroxisome and xenobiotic metabolism. Networks were validated by the striking enrichments that predicted regulatory lncRNAs showed for specific biological pathways. These findings highlight the power of snRNA-seq to discover functional roles for many xenobiotic-responsive lncRNAs in both hepatocytes and liver non-parenchymal cells and to elucidate novel aspects of foreign chemical-induced hepatotoxicity and liver disease, including dysregulation of intercellular communication within the liver lobule.

Sex-biased genetic programs in liver metabolism and liver fibrosis are controlled by EZH1 and EZH2.

Sex differences in the incidence and progression of many liver diseases, including liver fibrosis and hepatocellular carcinoma, are associated with sex-biased hepatic expression of hundreds of genes. This sexual dimorphism is largely determined by the sex-specific pattern of pituitary growth hormone secretion, which controls a transcriptional regulatory network operative in the context of sex-biased and growth hormone-regulated chromatin states. Histone H3K27-trimethylation yields a major sex-biased repressive chromatin mark deposited at many strongly female-biased genes in male mouse liver, but not at male-biased genes in female liver, and is catalyzed by polycomb repressive complex-2 through its homologous catalytic subunits, Ezh1 and Ezh2. Here, we used Ezh1-knockout mice with a hepatocyte-specific knockout of Ezh2 to investigate the sex bias of liver H3K27-trimethylation and its functional role in regulating sex-differences in the liver. Combined hepatic Ezh1/Ezh2 deficiency led to a significant loss of sex-biased gene expression, particularly in male liver, where many female-biased genes were increased in expression while male-biased genes showed decreased expression. The associated loss of H3K27me3 marks, and increases in the active enhancer marks H3K27ac and H3K4me1, were also more pronounced in male liver. Further, Ezh1/Ezh2 deficiency in male liver, and to a lesser extent in female liver, led to up regulation of many genes linked to liver fibrosis and liver cancer, which may contribute to the observed liver pathologies and the increased sensitivity of these mice to hepatotoxin exposure. Thus, Ezh1/Ezh2-catalyzed H3K27-trimethyation regulates sex-dependent genetic programs in liver metabolism and liver fibrosis through its sex-dependent effects on the epigenome, and may thereby determine the sex-bias in liver disease susceptibility.

Global analysis of expression, maturation and subcellular localization of mouse liver transcriptome identifies novel sex-biased and TCPOBOP-responsive long non-coding RNAs.

BACKGROUND: While nuclear transcription and RNA processing and localization are well established for protein coding genes (PCGs), these processes are poorly understood for long non-coding (lnc)RNAs. Here, we characterize global patterns of transcript expression, maturation and localization for mouse liver RNA, including more than 15,000 lncRNAs. PolyA-selected liver RNA was isolated and sequenced from four subcellular fractions (chromatin, nucleoplasm, total nucleus, and cytoplasm), and from the chromatin-bound fraction without polyA selection. RESULTS: Transcript processing, determined from normalized intronic to exonic sequence read density ratios, progressively increased for PCG transcripts in going from the chromatin-bound fraction to the nucleoplasm and then on to the cytoplasm. Transcript maturation was similar for lncRNAs in the chromatin fraction, but was significantly lower in the nucleoplasm and cytoplasm. LncRNA transcripts were 11-fold more likely to be significantly enriched in the nucleus than cytoplasm, and 100-fold more likely to be significantly chromatin-bound than nucleoplasmic. Sequencing chromatin-bound RNA greatly increased the sensitivity for detecting lowly expressed lncRNAs and enabled us to discover and localize hundreds of novel regulated liver lncRNAs, including lncRNAs showing sex-biased expression or responsiveness to TCPOBOP a xenobiotic agonist ligand of constitutive androstane receptor (Nr1i3). CONCLUSIONS: Integration of our findings with prior studies and lncRNA annotations identified candidate regulatory lncRNAs for a variety of hepatic functions based on gene co-localization within topologically associating domains or transcription divergent or antisense to PCGs associated with pathways linked to hepatic physiology and disease.

Loss of growth hormone-mediated signal transducer and activator of transcription 5 (STAT5) signaling in mice results in insulin sensitivity with obesity.

Growth hormone (GH) has an important function as an insulin antagonist with elevated insulin sensitivity evident in humans and mice lacking a functional GH receptor (GHR). We sought the molecular basis for this sensitivity by utilizing a panel of mice possessing specific deletions of GHR signaling pathways. Metabolic clamps and glucose homeostasis tests were undertaken in these obese adult C57BL/6 male mice, which indicated impaired hepatic gluconeogenesis. Insulin sensitivity and glucose disappearance rate were enhanced in muscle and adipose of mice lacking the ability to activate the signal transducer and activator of transcription (STAT)5 via the GHR (Ghr-391-/-) as for GHR-null (GHR-/-) mice. These changes were associated with a striking inhibition of hepatic glucose output associated with altered glycogen metabolism and elevated hepatic glycogen content during unfed state. The enhanced hepatic insulin sensitivity was associated with increased insulin receptor ß and insulin receptor substrate 1 activation along with activated downstream protein kinase B signaling cascades. Although phosphoenolpyruvate carboxykinase (Pck)-1 expression was unchanged, its inhibitory acetylation was elevated because of decreased sirtuin-2 expression, thereby promoting loss of PCK1. Loss of STAT5 signaling to defined chromatin immunoprecipitation targets would further increase lipogenesis, supporting hepatosteatosis while lowering glucose output. Finally, up-regulation of IL-15 expression in muscle, with increased secretion of adiponectin and fibroblast growth factor 1 from adipose tissue, is expected to promote insulin sensitivity.-Chhabra, Y., Nelson, C. N., Plescher, M., Barclay, J. L., Smith, A. G., Andrikopoulos, S., Mangiafico, S., Waxman, D. J., Brooks, A. J., Waters, M. J. Loss of growth hormone-mediated signal transducer and activator of transcription 5 (STAT5) signaling in mice results in insulin sensitivity with obesity.

Metronomic cyclophosphamide activation of anti-tumor immunity: tumor model, mouse host, and drug schedule dependence of gene responses and their upstream regulators.

BACKGROUND: Cyclophosphamide (CPA) can activate immunogenic tumor cell death, which induces immune-based tumor ablation and long-term anti-tumor immunity in a syngeneic C57BL/6 (B6) mouse GL261 glioma model when CPA is given on a 6-day repeating metronomic schedule (CPA/6d). In contrast, we find that two other syngeneic B6 mouse tumors, LLC lung carcinoma and B16F10 melanoma, do not exhibit these drug-induced immune responses despite their intrinsic sensitivity to CPA cytotoxicity. METHODS: To elucidate underlying mechanisms, we investigated gene expression and molecular pathway changes associated with the disparate immune responsiveness of these tumors to CPA/6d treatment. RESULTS: Global transcriptome analysis indicated substantial elevation of basal GL261 immune infiltration and strong CPA/6d activation of GL261 immune stimulatory pathways and their upstream regulators, but without preferential depletion of negative immune regulators compared to LLC and B16F10 tumors. In LLC tumors, where CPA/6d treatment is shown to be anti-angiogenic, CPA/6d suppressed VEGFA target genes and down regulated cell adhesion and leukocyte transendothelial migration genes. In GL261 tumors implanted in adaptive immune-deficient scid mice, where CPA/6d-induced GL261 regression is incomplete and late tumor growth rebound can occur, T cell receptor signaling and certain cytokine-cytokine receptor responses seen in B6 mice were deficient. Extending the CPA treatment interval from 6 to 9 days (CPA/9d) - which results in a strong but transient natural killer cell response followed by early tumor growth rebound - induced fewer cytokines and increased expression of drug metabolism genes. CONCLUSIONS: These findings elucidate molecular response pathways activated by intermittent metronomic CPA treatment and identify deficiencies that characterize immune-unresponsive tumor models and drug schedules.

Transcriptional profiling provides insights into metronomic cyclophosphamide-activated, innate immune-dependent regression of brain tumor xenografts.

BACKGROUND: Cyclophosphamide treatment on a six-day repeating metronomic schedule induces a dramatic, innate immune cell-dependent regression of implanted gliomas. However, little is known about the underlying mechanisms whereby metronomic cyclophosphamide induces innate immune cell mobilization and recruitment, or about the role of DNA damage and cell stress response pathways in eliciting the immune responses linked to tumor regression. METHODS: Untreated and metronomic cyclophosphamide-treated human U251 glioblastoma xenografts were analyzed on human microarrays at two treatment time points to identify responsive tumor cell-specific factors and their upstream regulators. Mouse microarray analysis across two glioma models (human U251, rat 9L) was used to identify host factors and gene networks that contribute to the observed immune and tumor regression responses. RESULTS: Metronomic cyclophosphamide increased expression of tumor cell-derived DNA damage, cell stress, and cell death genes, which may facilitate innate immune activation. Increased expression of many host (mouse) immune networks was also seen in both tumor models, including complement components, toll-like receptors, interferons, and cytolysis pathways. Key upstream regulators activated by metronomic cyclophosphamide include members of the interferon, toll-like receptor, inflammatory response, and PPAR signaling pathways, whose activation may contribute to anti-tumor immunity. Many upstream regulators inhibited by metronomic cyclophosphamide, including hypoxia-inducible factors and MAP kinases, have glioma-promoting activity; their inhibition may contribute to the therapeutic effectiveness of the six-day repeating metronomic cyclophosphamide schedule. CONCLUSIONS: Large numbers of responsive cytokines, chemokines and immune regulatory genes linked to innate immune cell recruitment and tumor regression were identified, as were several immunosuppressive factors that may contribute to the observed escape of some tumors from metronomic CPA-induced, immune-based regression. These factors may include useful biomarkers that facilitate discovery of clinically effective immunogenic metronomic drugs and treatment schedules, and the selection of patients most likely to be responsive to immunogenic drug scheduling.

Anti-tumor innate immunity activated by intermittent metronomic cyclophosphamide treatment of 9L brain tumor xenografts is preserved by anti-angiogenic drugs that spare VEGF receptor 2.

BACKGROUND: Metronomic cyclophosphamide given on an intermittent, 6-day repeating schedule, but not on an exposure dose-equivalent daily schedule, activates an anti-tumor innate immune response that leads to major regression of large implanted gliomas, without anti-angiogenesis. METHODS AND APPROACH: Mice bearing implanted 9L gliomas were used to investigate the effects of this 6-day repeating, immunogenic cyclophosphamide schedule on myeloid-derived suppressor cells, which are pro-angiogenic and can inhibit anti-tumor immunity, and to elucidate the mechanism whereby the innate immune cell-dependent tumor regression response to metronomic cyclophosphamide treatment is blocked by several anti-angiogenic receptor tyrosine kinase inhibitors. RESULTS: Intermittent metronomic cyclophosphamide scheduling strongly increased glioma-associated CD11b+ immune cells but not CD11b+Gr1+ myeloid-derived suppressor cells, while bone marrow and spleen reservoirs of the suppressor cells were decreased. The inhibition of immune cell recruitment and tumor regression by anti-angiogenic receptor tyrosine kinase inhibitors, previously observed in several brain tumor models, was recapitulated in the 9L tumor model with the VEGFR2-specific inhibitory monoclonal antibody DC101 (p < 0.01), implicating VEGFR2 signaling as an essential step in metronomic cyclophosphamide-stimulated immune cell recruitment. In contrast, sorafenib, a multi-receptor tyrosine kinase inhibitor with comparatively weak VEGF receptor phosphorylation inhibitory activity, was strongly anti-angiogenic but did not block metronomic cyclophosphamide-induced innate immunity or tumor regression (p > 0.05). CONCLUSIONS: The interference by receptor tyrosine kinase inhibitors in the immunogenic actions of intermittent metronomic chemotherapy is not a consequence of anti-angiogenesis per se, as demonstrated in an implanted 9L tumor model. Furthermore, this undesirable interaction with tyrosine kinase inhibitors can be avoided by using anti-angiogenic drugs that spare the VEGFR2 pathway.

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.

Steatotic liver disease induced by TCPOBOP-activated hepatic constitutive androstane receptor: primary and secondary gene responses with links to disease progression.

Constitutive Androstane Receptor (CAR, Nr1i3), a liver nuclear receptor and xenobiotic sensor, induces drug, steroid and lipid metabolizing enzymes, stimulates liver hypertrophy and hyperplasia, and ultimately, hepatocellular carcinogenesis. The mechanisms linking early CAR responses to later disease development are poorly understood. Here we show that exposure of CD-1 mice to TCPOBOP, a halogenated xenochemical and selective CAR agonist ligand, induces pericentral steatosis marked by hepatic accumulation of cholesterol and neutral lipid, and elevated circulating alanine aminotransferase, indicating hepatocyte damage. TCPOBOP-induced steatosis was weaker in the pericentral region but stronger in the periportal region in females compared to males. Early (1-day) TCPOBOP transcriptional responses were enriched for CAR-bound primary response genes, and for lipogenesis and xenobiotic metabolism and oxidative stress protection pathways; late (2-wk) TCPOBOP responses included many CAR binding-independent secondary response genes, with enrichment for macrophage activation, immune response and cytokine and reactive oxygen species production. Late upstream regulators specific to TCPOBOP-exposed male liver were linked to pro-inflammatory responses and hepatocellular carcinoma progression. TCPOBOP administered weekly to male mice using a high corn oil vehicle activated carbohydrate-responsive transcription factor (MLXIPL)-regulated target genes, dysregulated mitochondrial respiratory and translation regulatory pathways, and induced more advanced liver pathology. Overall, TCPOBOP exposure recapitulates histological and gene expression changes characteristic of emerging steatotic liver disease, including secondary gene responses in liver non-parenchymal cells indicative of transition to a more advanced disease state. Upstream regulators of both the early and late TCPOBOP response genes include novel biomarkers for foreign chemical-induced metabolic dysfunction-associated steatotic liver disease.

Impact of tumor blood flow modulation on tumor sensitivity to the bioreductive drug banoxantrone.

We investigated the hypoxia-dependent cytotoxicity of AQ4N (banoxantrone) using a panel of 13 cancer cell lines and studied its relationship to the expression of the quinone reductase DT-diaphorase (NQO1), which is widely found in cancer cells. We also investigated pharmacologic treatments that increase tumor hypoxia in vivo and their impact on AQ4N chemosensitivity in a solid tumor xenograft model. AQ4N showed ≥ 8-fold higher cytotoxicity under hypoxia than normoxia in cultures of 9L rat gliosarcoma and H460 human non-small-cell lung carcinoma cells but not for 11 other human cancer cell lines. DT-diaphorase protein levels and AQ4N chemosensitivity were poorly correlated across the cancer cell line panel, and AQ4N chemosensitivity was not affected by DT-diaphorase inhibitors. The vasodilator hydralazine decreased tumor perfusion and increased tumor hypoxia in 9L tumor xenografts, and to a lesser extent in H460 tumor xenografts. However, hydralazine did not increase AQ4N-dependent antitumor activity. Combination of AQ4N with the angiogenesis inhibitor axitinib, which increases 9L tumor hypoxia, transiently increased antitumor activity but with an increase in host toxicity. These findings indicate that the capacity to bioactivate AQ4N is not dependent on DT-diaphorase and is not widespread in cultured cancer cell lines. Moreover, the activation of AQ4N cytotoxicity in vivo requires tumor hypoxia that is more extensive or prolonged than can readily be achieved by vasodilation or by antiangiogenic drug treatment.

Changes in mouse uterine transcriptome in estrus and proestrus.

Changes in the CD-1 mouse uterine transcriptome during proestrus and estrus were investigated to help elucidate mechanisms of uterine tissue remodeling during the estrus cycle and their regulation by estrogen and progesterone in preparation of the uterus for pregnancy. Mice were staged beginning at 6 weeks of age, and uterine horns were harvested after monitoring two estrus cycles. Microarray analysis of whole uterine horn RNA identified 2428 genes differentially expressed in estrus compared to proestrus, indicating there is extensive remodeling of mouse uterus during the estrus cycle, affecting ~10% of all protein-encoding genes. Many (~50%) of these genes showed the same differential expression in independent analyses of isolated uterine lumenal epithelial cells. Changes in gene expression associated with structural alterations of the uterus included remodeling of the extracellular matrix, changes in cell keratins and adhesion molecules, activation of mitosis and changes in major histocompatibility complex class II (MHCII) presentation, complement and coagulation cascades, and cytochrome P450 expression. Signaling pathways regulated during the estrus cycle, involving ligand-gated channels, Wnt and hedgehog signaling, and transcription factors with poorly understood roles in reproductive tissues, included several genes and gene networks that have been implicated in pathological states. Many of the molecular pathways and biological functions represented by the genes differentially expressed from proestrus to estrus are also altered during the human menstrual cycle, although not necessarily at the corresponding phases of the cycle. These findings establish a baseline for further studies in the mouse model to dissect mechanisms involved in uterine tissue response to endocrine disruptors and the development of reproductive tract diseases.

TCDD dysregulation of lncRNA expression, liver zonation and intercellular communication across the liver lobule.

The persistent environmental aryl hydrocarbon receptor agonist and hepatotoxin TCDD (2,3,7,8-tetrachlorodibenzo- p -dioxin) induces hepatic lipid accumulation (steatosis), inflammation (steatohepatitis) and fibrosis. Thousands of liver-expressed, nuclear-localized lncRNAs with regulatory potential have been identified; however, their roles in TCDD-induced hepatoxicity and liver disease are unknown. We analyzed single nucleus (sn)RNA-seq data from control and chronic TCDD-exposed mouse liver to determine liver cell-type specificity, zonation and differential expression profiles for thousands of IncRNAs. TCDD dysregulated >4,000 of these lncRNAs in one or more liver cell types, including 684 lncRNAs specifically dysregulated in liver non-parenchymal cells. Trajectory inference analysis revealed major disruption by TCDD of hepatocyte zonation, affecting >800 genes, including 121 IncRNAs, with strong enrichment for lipid metabolism genes. TCDD also dysregulated expression of >200 transcription factors, including 19 Nuclear Receptors, most notably in hepatocytes and Kupffer cells. TCDD-induced changes in cellâ€"cell communication patterns included marked decreases in EGF signaling from hepatocytes to non-parenchymal cells and increases in extracellular matrix-receptor interactions central to liver fibrosis. Gene regulatory networks constructed from the snRNA-seq data identified TCDD-exposed liver network-essential lncRNA regulators linked to functions such as fatty acid metabolic process, peroxisome and xenobiotic metabolic. Networks were validated by the striking enrichments that predicted regulatory IncRNAs showed for specific biological pathways. These findings highlight the power of snRNA-seq to discover functional roles for many xenobiotic-responsive lncRNAs in both hepatocytes and liver non-parenchymal cells and to elucidate novel aspects of foreign chemical-induced hepatotoxicity and liver disease, including dysregulation of intercellular communication within the liver lobule.

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 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.

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.