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

BuDDI: BulkDeconvolution withDomainInvariance to predict cell-type-specific perturbations from bulk

While single-cell experiments provide deep cellular resolution within a single sample, some single-cell experiments are inherently more challenging than bulk experiments due to dissociation difficulties, cost, or limited tissue availability. This creates a situation where we have deep cellular profiles of one sample or condition, and bulk profiles across multiple samples and conditions. To bridge this gap, we propose BuDDI (BUlk Deconvolution with Domain Invariance). BuDDI utilizes domain adaptation techniques to effectively integrate available corpora of case-control bulk and reference scRNA-seq observations to infer cell-type-specific perturbation effects. BuDDI achieves this by learning independent latent spaces within a single variational autoencoder (VAE) encompassing at least four sources of variability: 1) cell type proportion, 2) perturbation effect, 3) structured experimental variability, and 4) remaining variability. Since each latent space is encouraged to be independent, we simulate perturbation responses by independently composing each latent space to simulate cell-type-specific perturbation responses. We evaluated BuDDI's performance on simulated and real data with experimental designs of increasing complexity. We first validated that BuDDI could learn domain invariant latent spaces on data with matched samples across each source of variability. Then we validated that BuDDI could accurately predict cell-type-specific perturbation response when no single-cell perturbed profiles were used during training; instead, only bulk samples had both perturbed and non-perturbed observations. Finally, we validated BuDDI on predicting sex-specific differences, an experimental design where it is not possible to have matched samples. In each experiment, BuDDI outperformed all other comparative methods and baselines. As more reference atlases are completed, BuDDI provides a path to combine these resources with bulk-profiled treatment or disease signatures to study perturbations, sex differences, or other factors at single-cell resolution.

Single-cell spatial multi-omics and deep learning dissect enhancer-driven gene regulatory networks in liver zonation

In the mammalian liver, hepatocytes exhibit diverse metabolic and functional profiles based on their location within the liver lobule. However, it is unclear whether this spatial variation, called zonation, is governed by a well-defined gene regulatory code. Here, using a combination of single-cell multiomics, spatial omics, massively parallel reporter assays and deep learning, we mapped enhancer-gene regulatory networks across mouse liver cell types. We found that zonation affects gene expression and chromatin accessibility in hepatocytes, among other cell types. These states are driven by the repressors TCF7L1 and TBX3, alongside other core hepatocyte transcription factors, such as HNF4A, CEBPA, FOXA1 and ONECUT1. To examine the architecture of the enhancers driving these cell states, we trained a hierarchical deep learning model called DeepLiver. Our study provides a multimodal understanding of the regulatory code underlying hepatocyte identity and their zonation state that can be used to engineer enhancers with specific activity levels and zonation patterns.

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.

How single-cell transcriptomics provides insight on hepatic responses to TCDD

The prototypical aryl hydrocarbon receptor (AHR) ligand, 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), has been a valuable model for investigating toxicant-associated fatty liver disease (TAFLD). TCDD induces dose-dependent hepatic lipid accumulation, followed by the development of inflammatory foci and eventual progression to fibrosis in mice. Previously, bulk approaches and in vitro examination of different cell types were relied upon to study the mechanisms underlying TCDD-induced liver pathologies. However, the advent of single-cell transcriptomic technologies, such as single-nuclei RNA sequencing (snRNAseq) and spatial transcriptomics (STx), has provided new insights into the responses of hepatic cell types to TCDD exposure. This review explores the application of these single-cell transcriptomic technologies and highlights their contributions towards unraveling the cell-specific mechanisms mediating the hepatic responses to TCDD.

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.

Single-cell dynamics of liver development in postnatal pigs

The postnatal development of the liver, an essential organ for metabolism and immunity, remains poorly characterized at the single-cell resolution. Here, we generated single-nucleus and single-cell transcriptomes of 84,824 pig liver cells at four postnatal time points: day 30, 42, 150, and 730. We uncovered 23 cell types, including three rare cell types: plasmacytoid dendritic cells, CAVIN3+IGF2+ endothelial cells, and EBF1+ fibroblasts. The latter two were verified by multiplex immunohistochemistry. Trajectory and gene regulatory analyses revealed 33 genes that encode transcription factors associated with hepatocyte development and function, including NFIL3 involved in regulating hepatic metabolism. We characterized the spatiotemporal heterogeneity of liver endothelial cells, identified and validated leucine zipper protein 2 (LUZP2) as a novel adult liver sinusoidal endothelial cell-specific transcription factor. Lymphoid cells (NK and T cells) governed the immune system of the pig liver since day 30. Furthermore, we identified a cluster of tissue-resident NK cells, which displayed virus defense functions, maintained proliferative features at day 730, and manifested a higher conservative transcription factor expression pattern in humans than in mouse liver. Our study presents the most comprehensive postnatal liver development single-cell atlas and demonstrates the metabolic and immune changes across the four age stages.

Sex-biased gene expression across mammalian organ development and evolution

Sexually dimorphic traits are common among mammals and are specified during development through the deployment of sex-specific genetic programs. Because little is known about these programs, we investigated them using a resource of gene expression profiles in males and females throughout the development of five organs in five mammals (human, mouse, rat, rabbit, and opossum) and a bird (chicken). We found that sex-biased gene expression varied considerably across organs and species and was often cell-type specific. Sex differences increased abruptly around sexual maturity instead of increasing gradually during organ development. Finally, sex-biased gene expression evolved rapidly at the gene level, with differences between organs in the evolutionary mechanisms used, but more slowly at the cellular level, with the same cell types being sexually dimorphic across species.

GHR disruption in mature adult mice alters xenobiotic metabolism gene expression in the liver

BACKGROUND

Lifelong reduction of growth hormone (GH) action extends lifespan and improves healthspan in mice. Moreover, congenital inactivating mutations of GH receptor (GHR) in mice and humans impart resistance to age-associated cancer, diabetes, and cognitive decline. To investigate the consequences of GHR disruption at an adult age, we recently ablated the GHR at 6-months of age in mature adult (6mGHRKO) mice. We found that both, male and female 6mGHRKO mice have reduced oxidative damage, with males 6mGHRKO showing improved insulin sensitivity and cancer resistance. Importantly, 6mGHRKO females have an extended lifespan compared to controls.

OBJECTIVE AND METHODS

To investigate the possible mechanisms leading to health improvements, we performed RNA sequencing using livers from male and female 6mGHRKO mice and controls.

RESULTS

We found that disrupting GH action at an adult age reduced the gap in liver gene expression between males and females, making gene expression between sexes more similar. However, there was still a 6-fold increase in the number of differentially expressed genes when comparing male 6mGHRKO mice vs controls than in 6mGHRKO female vs controls, suggesting that GHR ablation affects liver gene expression more in males than in females. Finally, we found that lipid metabolism and xenobiotic metabolism pathways are activated in the liver of 6mGHRKO mice.

CONCLUSION

The present study shows for the first time the specific hepatic gene expression profile, cellular pathways, biological processes and molecular mechanisms that are driven by ablating GH action at a mature adult age in males and females. Importantly, these results and future studies on xenobiotic metabolism may help explain the lifespan extension seen in 6mGHRKO mice.

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.

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.

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.

Hepatic Energy Metabolism under the Local Control of the Thyroid Hormone System

The energy homeostasis of the organism is orchestrated by a complex interplay of energy substrate shuttling, breakdown, storage, and distribution. Many of these processes are interconnected via the liver. Thyroid hormones (TH) are well known to provide signals for the regulation of energy homeostasis through direct gene regulation via their nuclear receptors acting as transcription factors. In this comprehensive review, we summarize the effects of nutritional intervention like fasting and diets on the TH system. In parallel, we detail direct effects of TH in liver metabolic pathways with regards to glucose, lipid, and cholesterol metabolism. This overview on hepatic effects of TH provides the basis for understanding the complex regulatory network and its translational potential with regards to currently discussed treatment options of non-alcoholic fatty liver disease (NAFLD) and non-alcoholic steatohepatitis (NASH) involving TH mimetics.

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.

Chromosomal and gonadal sex drive sex differences in lipids and hepatic gene expression in response to hypercholesterolemia and statin treatment

BACKGROUND

Biological sex impacts susceptibility and presentation of cardiovascular disease, which remains the leading cause of death for both sexes. To reduce cardiovascular disease risk, statin drugs are commonly prescribed to reduce circulating cholesterol levels through inhibition of cholesterol synthesis. The effectiveness of statin therapy differs between individuals with a sex bias in the frequency of adverse effects. Limited information is available regarding the mechanisms driving sex-specific responses to hypercholesterolemia or statin treatment.

METHODS

Four Core Genotypes mice (XX and XY mice with ovaries and XX and XY mice with testes) on a hypercholesteremic Apoe-/- background were fed a chow diet without or with simvastatin for 8 weeks. Plasma lipid levels were quantified and hepatic differential gene expression was evaluated with RNA-sequencing to identify the independent effects of gonadal and chromosomal sex.

RESULTS

In a hypercholesterolemic state, gonadal sex influenced the expression levels of more than 3000 genes, and chromosomal sex impacted expression of nearly 1400 genes, which were distributed across all autosomes as well as the sex chromosomes. Gonadal sex uniquely influenced the expression of ER stress response genes, whereas chromosomal and gonadal sex influenced fatty acid metabolism gene expression in hypercholesterolemic mice. Sex-specific effects on gene regulation in response to statin treatment included a compensatory upregulation of cholesterol biosynthetic gene expression in mice with XY chromosome complement, regardless of presence of ovaries or testes.

CONCLUSION

Gonadal and chromosomal sex have independent effects on the hepatic transcriptome to influence different cellular pathways in a hypercholesterolemic environment. Furthermore, chromosomal sex in particular impacted the cellular response to statin treatment. An improved understanding of how gonadal and chromosomal sex influence cellular response to disease conditions and in response to drug treatment is critical to optimize disease management for all individuals.