Dynamic, sex-differential STAT5 and BCL6 binding to sex-biased, growth hormone-regulated genes in adult mouse liver

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.

Transcriptional programming of translation by BCL6 controls skeletal muscle proteostasis

Skeletal muscle is dynamically controlled by the balance of protein synthesis and degradation. Here we discover an unexpected function for the transcriptional repressor B cell lymphoma 6 (BCL6) in muscle proteostasis and strength in mice. Skeletal muscle-specific Bcl6 ablation in utero or in adult mice results in over 30% decreased muscle mass and force production due to reduced protein synthesis and increased autophagy, while it promotes a shift to a slower myosin heavy chain fibre profile. Ribosome profiling reveals reduced overall translation efficiency in Bcl6-ablated muscles. Mechanistically, tandem chromatin immunoprecipitation, transcriptomic and translational analyses identify direct BCL6 repression of eukaryotic translation initiation factor 4E-binding protein 1 (Eif4ebp1) and activation of insulin-like growth factor 1 (Igf1) and androgen receptor (Ar). Together, these results uncover a bifunctional role for BCL6 in the transcriptional and translational control of muscle proteostasis.

Tmem263 deletion disrupts the GH/IGF-1 axis and causes dwarfism and impairs skeletal acquisition

Genome-wide association studies (GWAS) have identified a large number of candidate genes believed to affect longitudinal bone growth and bone mass. One of these candidate genes, TMEM263, encodes a poorly characterized plasma membrane protein. Single nucleotide polymorphisms in TMEM263 are associated with bone mineral density in humans and mutations are associated with dwarfism in chicken and severe skeletal dysplasia in at least one human fetus. Whether this genotype-phenotype relationship is causal, however, remains unclear. Here, we determine whether and how TMEM263 is required for postnatal growth. Deletion of the Tmem263 gene in mice causes severe postnatal growth failure, proportional dwarfism, and impaired skeletal acquisition. Mice lacking Tmem263 show no differences in body weight within the first 2 weeks of postnatal life. However, by P21 there is a dramatic growth deficit due to a disrupted growth hormone (GH)/insulin-like growth factor 1 (IGF-1) axis, which is critical for longitudinal bone growth. Tmem263-null mice have low circulating IGF-1 levels and pronounced reductions in bone mass and growth plate length. The low serum IGF-1 in Tmem263-null mice is associated with reduced hepatic GH receptor (GHR) expression and GH-induced JAK2/STAT5 signaling. A deficit in GH signaling dramatically alters GH-regulated genes and feminizes the liver transcriptome of Tmem263-null male mice, with their expression profile resembling wild-type female, hypophysectomized male, and Stat5b-null male mice. Collectively, our data validates the causal role for Tmem263 in regulating postnatal growth and raises the possibility that rare mutations or variants of TMEM263 may potentially cause GH insensitivity and impair linear growth.

Peak Scores Significantly Depend on the Relationships between Contextual Signals in ChIP-Seq Peaks

Chromatin immunoprecipitation followed by massively parallel DNA sequencing (ChIP-seq) is a central genome-wide method for in vivo analyses of DNA-protein interactions in various cellular conditions. Numerous studies have demonstrated the complex contextual organization of ChIP-seq peak sequences and the presence of binding sites for transcription factors in them. We assessed the dependence of the ChIP-seq peak score on the presence of different contextual signals in the peak sequences by analyzing these sequences from several ChIP-seq experiments using our fully enumerative GPU-based de novo motif discovery method, Argo_CUDA. Analysis revealed sets of significant IUPAC motifs corresponding to the binding sites of the target and partner transcription factors. For these ChIP-seq experiments, multiple regression models were constructed, demonstrating a significant dependence of the peak scores on the presence in the peak sequences of not only highly significant target motifs but also less significant motifs corresponding to the binding sites of the partner transcription factors. A significant correlation was shown between the presence of the target motifs FOXA2 and the partner motifs HNF4G, which found experimental confirmation in the scientific literature, demonstrating the important contribution of the partner transcription factors to the binding of the target transcription factor to DNA and, consequently, their important contribution to the peak score.

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.

Tmem263 deletion disrupts the GH/IGF-1 axis and causes dwarfism and impairs skeletal acquisition

Genome-wide association studies (GWAS) have identified a large number of candidate genes believed to affect longitudinal bone growth and bone mass. One of these candidate genes, TMEM263, encodes a poorly characterized plasma membrane protein. Single nucleotide polymorphisms in TMEM263 are associated with bone mineral density in humans and mutations are associated with dwarfism in chicken and severe skeletal dysplasia in at least one human fetus. Whether this genotype-phenotype relationship is causal, however, remains unclear. Here, we determine whether and how TMEM263 is required for postnatal growth. Deletion of the Tmem263 gene in mice causes severe postnatal growth failure, proportional dwarfism, and impaired skeletal acquisition. Mice lacking Tmem263 show no differences in body weight within the first two weeks of postnatal life. However, by P21 there is a dramatic growth deficit due to a disrupted GH/IGF-1 axis, which is critical for longitudinal bone growth. Tmem263-null mice have low circulating IGF-1 levels and pronounced reductions in bone mass and growth plate length. The low serum IGF-1 in Tmem263-null mice is associated with reduced hepatic GH receptor (GHR) expression and GH-induced JAK2/STAT5 signaling. A deficit in GH signaling dramatically alters GH-regulated genes and feminizes the liver transcriptome of Tmem263-null male mice, with their expression profile resembling a wild-type female, hypophysectomized male, and Stat5b-null male mice. Collectively, our data validates the causal role for Tmem263 in regulating postnatal growth and raises the possibility that rare mutations or variants of TMEM263 may potentially cause GH insensitivity and impair linear growth.

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.

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.

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.

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.

Mice with humanized livers reveal the role of hepatocyte clocks in rhythmic behavior

The synchronization of circadian clock depends on a central pacemaker located in the suprachiasmatic nuclei. However, the potential feedback of peripheral signals on the central clock remains poorly characterized. To explore whether peripheral organ circadian clocks may affect the central pacemaker, we used a chimeric model in which mouse hepatocytes were replaced by human hepatocytes. Liver humanization led to reprogrammed diurnal gene expression and advanced the phase of the liver circadian clock that extended to muscle and the entire rhythmic physiology. Similar to clock-deficient mice, liver-humanized mice shifted their rhythmic physiology more rapidly to the light phase under day feeding. Our results indicate that hepatocyte clocks can affect the central pacemaker and offer potential perspectives to apprehend pathologies associated with altered circadian physiology.

Androgen-Dependent Expression of CUX2 mRNA in the Pig Liver Is Associated with That of Drug Metabolizing Enzymes and Drug Transporters

We previously identified androgen-dependent sex differences in the mRNA expression of drug metabolizing enzymes (DMEs), including CYPs, sulfotransferases and uridine 5'-diphospho-glucuronosyltransferases, and drug transporters in the pig liver and kidney. To elucidate the mechanism for such sex differences in pigs, we herein focused on the key regulators cut-like homeobox 2 (Cux2), B-cell lymphoma 6 (Bcl6), and signal transducer and activator of transcription 5b (Stat5b), which are reported to be responsible for the sex-biased gene expression of Cyps in the mouse liver. We used real-time RT-PCR to examine androgen-dependent sex differences in the mRNA levels of these regulators in the liver and kidney basically using Meishan and Landrace pigs. Significant sex differences (male > female) in the level of CUX2 mRNA were detected in the liver of both breeds, and levels were significantly decreased in males by castration and increased in castrated males and intact females by administering testosterone propionate. No such clear androgen-dependent sex differences in hepatic BCL6 or STAT5B mRNA expression were observed in either breed. In the kidney, androgen-dependent gene expression of these regulators was not observed. In the liver, CUX2 mRNA expression closely correlated with that of DMEs and drug transporters, which were previously shown to have androgen-dependent expression. Together, these findings demonstrate that hepatic CUX2 mRNA is expressed in an androgen-dependent manner, and strongly suggest that CUX2 plays a key role in the androgen-dependent gene expression of hepatic DMEs and drug transporters.

Liver Injury and Cell Survival in Non-Alcoholic Steatohepatitis Regulated by Sex-Based Difference through B Cell Lymphoma 6

The liver is a crucial organ for maintaining homeostasis in living organisms and is the center of various metabolic functions. Therefore, abnormal metabolic activity, as in metabolic syndrome, leads to pathological conditions, such as abnormal accumulation of lipids in the liver. Inflammation and cell death are induced by several stresses in the fatty liver, namely steatohepatitis. In recent years, an increase in non-alcoholic steatohepatitis (NASH), which is not dependent on excessive alcohol intake, has become an issue as a major cause of liver cirrhosis and liver cancer. There are several recent findings on functional sex-based differences, NASH, and cell stress and death in the liver. In particular, NASH-induced liver injury and tumorigeneses were suppressed by B cell lymphoma 6, the transcriptional factor regulating sex-based liver functional gene expression. In this review, we discuss cell response to stress and lipotoxicity in NASH and its regulatory mechanisms.

NIK/MAP3K14 in hepatocytes orchestrates NASH to hepatocellular carcinoma progression via JAK2/STAT5 inhibition

OBJECTIVE

Nonalcoholic fatty liver disease (NAFLD) ranges from steatosis to nonalcoholic steatohepatitis (NASH), which often progresses to hepatocellular carcinoma (HCC) through a largely undefined mechanism. NASH and HCC depend on inflammatory signaling, whose master regulator is the NFκB transcription factor family, activated by canonical and non-canonical pathways.

METHODS

Here, we investigated non-canonical NFκB-inducing kinase (NIK/MAP3K14) in metabolic NASH, NASH to HCC transition, and DEN-induced HCC. To this end, we performed dietary and chemical interventions in mice that were analyzed via single nucleus sequencing, gene expression and histochemical methods. Ultimately, we verified our mouse results in human patient samples.

RESULTS

We revealed that hepatocyte-specific NIK deficiency (NIKLKO) ameliorated metabolic NASH complications and reduced hepatocarcinogenesis, independent of its role in the NFκB pathway. Instead, hepatic NIK attenuated hepatoprotective JAK2/STAT5 signaling that is a prerequisite for NASH and NASH to HCC progression in mice and humans.

CONCLUSIONS

Our data suggest NIK-mediated inhibitory JAK2 phosphorylation at serine 633 that might be amenable for future therapeutic interventions in patients.

An evolutionary trade-off between host immunity and metabolism drives fatty liver in male mice

Adaptations to infectious and dietary pressures shape mammalian physiology and disease risk. How such adaptations affect sex-biased diseases remains insufficiently studied. In this study, we show that sex-dependent hepatic gene programs confer a robust (~300%) survival advantage for male mice during lethal bacterial infection. The transcription factor B cell lymphoma 6 (BCL6), which masculinizes hepatic gene expression at puberty, is essential for this advantage. However, protection by BCL6 protein comes at a cost during conditions of dietary excess, which result in overt fatty liver and glucose intolerance in males. Deleting hepatic BCL6 reverses these phenotypes but markedly lowers male survival during infection, thus establishing a sex-dependent trade-off between host defense and metabolic systems. Our findings offer strong evidence that some current sex-biased diseases are rooted in ancient evolutionary trade-offs between immunity and metabolism.

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.

Genome-wide association study reveals the genetic basis of growth trait in yellow catfish with sexual size dimorphism

Sexual size dimorphism has been widely observed in a large number of animals including fish species. Genome-wide association study (GWAS) is a powerful tool to dissect the genetic basis of complex traits, whereas the sex-differences in the genomics of animal complex traits have been ignored in the GWAS analysis. Yellow catfish (Pelteobagrus fulvidraco) is an important aquaculture fish in China with significant sexual size dimorphism. In this study, GWAS was conducted to identify candidate SNPs and genes related to body length (BL) and body weight (BW) in 125 female yellow catfish from a breeding population. In total, one BL-related SNP and three BW-related SNPs were identified to be significantly associated with the traits. Besides, one of these SNPs (Chr15:19195072) was shared in both the BW and BL traits in female yellow catfish, which was further validated in 185 male individuals and located on the exon of stat5b gene. Transgenic yellow catfish and zebrafish that expressed yellow catfish stat5b showed increased growth rate and reduction of sexual size dimorphism. These results not only reveal the genetic basis of growth trait and sexual size dimorphism in fish species, but also provide useful information for the marker-assisted breeding in yellow catfish.

Constitutively Active STAT5b Feminizes Mouse Liver Gene Expression

STAT5 is an essential transcriptional regulator of the sex-biased actions of GH in the liver. Delivery of constitutively active STAT5 (STAT5CA) to male mouse liver using an engineered adeno-associated virus with high tropism for the liver is shown to induce widespread feminization of the liver, with extensive induction of female-biased genes and repression of male-biased genes, largely mimicking results obtained when male mice are given GH as a continuous infusion. Many of the STAT5CA-responding genes were associated with nearby (< 50 kb) sites of STAT5 binding to liver chromatin, supporting the proposed direct role of persistently active STAT5 in continuous GH-induced liver feminization. The feminizing effects of STAT5CA were dose-dependent; moreover, at higher levels, STAT5CA overexpression resulted in some histopathology, including hepatocyte hyperplasia, and increased karyomegaly and multinuclear hepatocytes. These findings establish that the persistent activation of STAT5 by GH that characterizes female liver is by itself sufficient to account for the sex-dependent expression of a majority of hepatic sex-biased genes. Moreover, histological changes seen when STAT5CA is overexpressed highlight the importance of carefully evaluating such effects before considering STAT5 derivatives for therapeutic use in treating liver disease.

B-cell lymphoma 6 alleviates nonalcoholic fatty liver disease in mice through suppression of fatty acid transporter CD36

Nonalcoholic fatty liver disease (NAFLD) is an ubiquitous disease that exists across a wide spectrum ranging from steatosis, steatohepatitis, advanced fibrosis, and liver cirrhosis. Hallmarks of NAFLD are lipid accumulation, insulin resistance, and chronic low-grade inflammation. However, there currently are no medications approved for NAFLD. B-cell lymphoma 6 (BCL6) is a transcriptional inhibitor that is vital for germinal center B-cell formation. Our study identified BCL6 as a critical modulator of hepatic lipid metabolism and appears to contribute to the initiation and progression of NAFLD. In our research, we induced hepatic BCL6 overexpression using adeno-associated virus (AAV), as well as conditional liver-specific BCL6 knockout mice (BCL6-CKO). With these models, we noted that BCL6 overexpression improved insulin resistance and hepatic steatosis in mice models maintained on a HFD diet. Conversely, these parameters worsened in the livers of mice with downregulated BCL6 levels. Mechanistically, the translocase fatty acid CD36 was determined to be a transcriptional target of BCL6 that influences its role in hepatic steatosis. BCL6 bound directly to the CD36 promoter region, restraining CD36 transcription under physiological conditions. We conclude that the hepatocyte BCL6 inhibits the NAFLD progression in mice, including deranged lipid accumulation and glucose metabolism, through a CD36-dependent manner. These results indicate that BCL6 may potentially be targeted in NAFLD treatment.

Effects of Overexpression of Fibroblast Growth Factor 15/19 on Hepatic Drug Metabolizing Enzymes

Fibroblast growth factors 15 (FGF15) and 19 (FGF19) are endocrine growth factors that play an important role in maintaining bile acid homeostasis. FGF15/19-based therapies are currently being tested in clinical trials for the treatment of nonalcoholic steatohepatitis and cholestatic liver diseases. To determine the physiologic impact of long-term elevations of FGF15/19, a transgenic mouse model with overexpression of Fgf15 (Fgf15 Tg) was used in the current study. The RNA sequencing (RNA-seq) analysis revealed elevations of the expression of several genes encoding phase I drug metabolizing enzymes (DMEs), including Cyp2b10 and Cyp3a11, in Fgf15 Tg mice. We found that the induction of several Cyp2b isoforms resulted in increased function of CYP2B in microsomal metabolism and pharmacokinetics studies. Because the CYP2B family is known to be induced by constitutive androstane receptor (CAR), to determine the role of CAR in the observed inductions, we crossed Fgf15 Tg mice with CAR knockout mice and found that CAR played a minor role in the observed alterations in DME expression. Interestingly, we found that the overexpression of Fgf15 in male mice resulted in a phenotypical switch from the male hepatic expression pattern of DMEs to that of female mice. Differences in secretion of growth hormone (GH) between male and female mice are known to drive sexually dimorphic, STAT5b-dependent expression patterns of hepatic genes. We found that male Fgf15 Tg mice presented with many features similar to GH deficiency, including lowered body length and weight, Igf-1 and Igfals expression, and STAT5 signaling. SIGNIFICANCE STATEMENT: The overexpression of Fgf15 in mice causes an alteration in DMEs at the mRNA, protein, and functional levels, which is not entirely due to CAR activation but associated with lower GH signaling.

Regulation of Sexually Dimorphic Expression of Major Urinary Proteins

Male house mice excrete large amounts of protein in their urinary scent marks, mainly composed of Major Urinary Proteins (MUPs), and these lipocalins function as pheromones and pheromone carriers. Here, we review studies on sexually dimorphic MUP expression in house mice, including the proximate mechanisms controlling MUP gene expression and their adaptive functions. Males excrete 2 to 8 times more urinary protein than females, though there is enormous variation in gene expression across loci in both sexes. MUP expression is dynamically regulated depending upon a variety of factors. Males regulate MUP expression according to social status, whereas females do not, and males regulate expression depending upon health and condition. Male-biased MUP expression is regulated by pituitary secretion of growth hormone (GH), which binds receptors in the liver, activating the JAK2-STAT5 signaling pathway, chromatin accessibility, and MUP gene transcription. Pulsatile male GH secretion is feminized by several factors, including caloric restriction, microbiota depletion, and aging, which helps explain condition-dependent MUP expression. If MUP production has sex-specific fitness optima, then this should generate sexual antagonism over allelic expression (intra-locus sexual conflict) selectively favoring sexually dimorphic expression. MUPs influence the sexual attractiveness of male urinary odor and increased urinary protein excretion is correlated with the reproductive success of males but not females. This finding could explain the selective maintenance of sexually dimorphic MUP expression. Producing MUPs entails energetic costs, but increased excretion may reduce the net energetic costs and predation risks from male scent marking as well as prolong the release of chemical signals. MUPs may also provide physiological benefits, including regulating metabolic rate and toxin removal, which may have sex-specific effects on survival. A phylogenetic analysis on the origins of male-biased MUP gene expression in Mus musculus suggests that this sexual dimorphism evolved by increasing male MUP expression rather than reducing female expression.

Disruption of the circadian clock component BMAL1 elicits an endocrine adaption impacting on insulin sensitivity and liver disease

SignificanceWhile increasing evidence associates the disruption of circadian rhythms with pathologic conditions, including obesity, type 2 diabetes, and nonalcoholic fatty liver diseases (NAFLD), the involved mechanisms are still poorly described. Here, we show that, in both humans and mice, the pathogenesis of NAFLD is associated with the disruption of the circadian clock combined with perturbations of the growth hormone and sex hormone pathways. However, while this condition protects mice from the development of fibrosis and insulin resistance, it correlates with increased fibrosis in humans. This suggests that the perturbation of the circadian clock and its associated disruption of the growth hormone and sex hormone pathways are critical for the pathogenesis of metabolic and liver diseases.

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.

Towards Understanding the Direct and Indirect Actions of Growth Hormone in Controlling Hepatocyte Carbohydrate and Lipid Metabolism

Growth hormone (GH) is critical for achieving normal structural growth. In addition, GH plays an important role in regulating metabolic function. GH acts through its GH receptor (GHR) to modulate the production and function of insulin-like growth factor 1 (IGF1) and insulin. GH, IGF1, and insulin act on multiple tissues to coordinate metabolic control in a context-specific manner. This review will specifically focus on our current understanding of the direct and indirect actions of GH to control liver (hepatocyte) carbohydrate and lipid metabolism in the context of normal fasting (sleep) and feeding (wake) cycles and in response to prolonged nutrient deprivation and excess. Caveats and challenges related to the model systems used and areas that require further investigation towards a clearer understanding of the role GH plays in metabolic health and disease are discussed.

Beyond the X Factor: Relevance of Sex Hormones in NAFLD Pathophysiology

Non-alcoholic fatty liver disease (NAFLD) is a major health issue worldwide, being frequently associated with obesity, unbalanced dietary regimens, and reduced physical activity. Despite their greater adiposity and reduced physical activity, women show a lower risk of developing NAFLD in comparison to men, likely a consequence of a sex-specific regulation of liver metabolism. In the liver, sex differences in the uptake, synthesis, oxidation, deposition, and mobilization of lipids, as well as in the regulation of inflammation, are associated with differences in NAFLD prevalence and progression between men and women. Given the major role of sex hormones in driving hepatic sexual dimorphism, this review will focus on the role of sex hormones and their signaling in the regulation of hepatic metabolism and in the molecular mechanisms triggering NAFLD development and progression.

Sex Dimorphism of Nonalcoholic Fatty Liver Disease (NAFLD) in Pparg-Null Mice

Men with nonalcoholic fatty liver disease (NAFLD) are more exposed to nonalcoholic steatohepatitis (NASH) and liver fibrosis than women. However, the underlying molecular mechanisms of NALFD sex dimorphism are unclear. We combined gene expression, histological and lipidomic analyses to systematically compare male and female liver steatosis. We characterized hepatosteatosis in three independent mouse models of NAFLD, ob/ob and lipodystrophic fat-specific (PpargFΔ/Δ) and whole-body PPARγ-null (PpargΔ/Δ) mice. We identified a clear sex dimorphism occurring only in PpargΔ/Δ mice, with females showing macro- and microvesicular hepatosteatosis throughout their entire life, while males had fewer lipid droplets starting from 20 weeks. This sex dimorphism in hepatosteatosis was lost in gonadectomized PpargΔ/Δ mice. Lipidomics revealed hepatic accumulation of short and highly saturated TGs in females, while TGs were enriched in long and unsaturated hydrocarbon chains in males. Strikingly, sex-biased genes were particularly perturbed in both sexes, affecting lipid metabolism, drug metabolism, inflammatory and cellular stress response pathways. Most importantly, we found that the expression of key sex-biased genes was severely affected in all the NAFLD models we tested. Thus, hepatosteatosis strongly affects hepatic sex-biased gene expression. With NAFLD increasing in prevalence, this emphasizes the urgent need to specifically address the consequences of this deregulation in humans.

Distinct roles of androgen receptor, estrogen receptor alpha, and BCL6 in the establishment of sex-biased DNA methylation in mouse liver

Sexual dimorphism in gene regulation, including DNA methylation, is the main driver of sexual dimorphism in phenotypes. However, the questions of how and when sex shapes DNA methylation remain unresolved. Recently, using mice with different combinations of genetic and phenotypic sex, we identified sex-associated differentially methylated regions (sDMRs) that depended on the sex phenotype. Focusing on a panel of validated sex-phenotype dependent male- and female-biased sDMRs, we tested the developmental dynamics of sex bias in liver methylation and the impacts of mutations in the androgen receptor, estrogen receptor alpha, or the transcriptional repressor Bcl6 gene. True hermaphrodites that carry both unilateral ovaries and contralateral testes were also tested. Our data show that sex bias in methylation either coincides with or follows sex bias in the expression of sDMR-proximal genes, suggesting that sex bias in gene expression may be required for demethylation at certain sDMRs. Global ablation of AR, ESR1, or a liver-specific loss of BCL6, all alter sDMR methylation, whereas presence of both an ovary and a testis delays the establishment of male-type methylation levels in hermaphrodites. Moreover, the Bcl6-LKO shows dissociation between expression and methylation, suggesting a distinct role of BCL6 in demethylation of intragenic sDMRs.

Transcriptional regulation of N6-methyladenosine orchestrates sex-dimorphic metabolic traits

Males and females exhibit striking differences in the prevalence of metabolic traits including hepatic steatosis, a key driver of cardiometabolic morbidity and mortality. RNA methylation is a widespread regulatory mechanism of transcript turnover. Here, we show that presence of the RNA modification N⁶-methyladenosine (m⁶A) triages lipogenic transcripts for degradation and guards against hepatic triglyceride accumulation. In male but not female mice, this protective checkpoint stalls under lipid-rich conditions. Loss of m⁶A control in male livers increases hepatic triglyceride stores, leading to a more 'feminized' hepatic lipid composition. Crucially, liver-specific deletion of the m⁶A complex protein Mettl14 from male and female mice significantly diminishes sex-specific differences in steatosis. We further surmise that the m⁶A installing machinery is subject to transcriptional control by the sex-responsive BCL6-STAT5 axis in response to dietary conditions. These data show that m⁶A is essential for precise and synchronized control of lipogenic enzyme activity and provide insights into the molecular basis for the existence of sex-specific differences in hepatic lipid traits.

Increased accuracy of genomic predictions for growth under chronic thermal stress in rainbow trout by prioritizing variants from GWAS using imputed sequence data

Through imputation of genotypes, genome‐wide association study (GWAS) and genomic prediction (GP) using whole‐genome sequencing (WGS) data are cost‐efficient and feasible in aquaculture breeding schemes. The objective was to dissect the genetic architecture of growth traits under chronic heat stress in rainbow trout (Oncorhynchus mykiss) and to assess the accuracy of GP based on imputed WGS and different preselected single nucleotide polymorphism (SNP) arrays. A total of 192 and 764 fish challenged to a heat stress experiment for 62 days were genotyped using a customized 1 K and 26 K SNP panels, respectively, and then, genotype imputation was performed from a low‐density chip to WGS using 102 parents (36 males and 66 females) as the reference population. Imputed WGS data were used to perform GWAS and test GP accuracy under different preselected SNP scenarios. Heritability was estimated for body weight (BW), body length (BL) and average daily gain (ADG). Estimates using imputed WGS data ranged from 0.33 ± 0.05 to 0.55 ± 0.05 for growth traits under chronic heat stress. GWAS revealed that the top five cumulatively SNPs explained a maximum of 0.94%, 0.86% and 0.51% of genetic variance for BW, BL and ADG, respectively. Some important functional candidate genes associated with growth‐related traits were found among the most important SNPs, including signal transducer and activator of transcription 5B and 3 (STAT5B and STAT3, respectively) and cytokine‐inducible SH2‐containing protein (CISH). WGS data resulted in a slight increase in prediction accuracy compared with pedigree‐based method, whereas preselected SNPs based on the top GWAS hits improved prediction accuracies, with values ranging from 1.2 to 13.3%. Our results support the evidence of the polygenic nature of growth traits when measured under heat stress. The accuracies of GP can be improved using preselected variants from GWAS, and the use of WGS marginally increases prediction accuracy.

Hepatic transcriptional responses to fasting and feeding

Mammals undergo regular cycles of fasting and feeding that engage dynamic transcriptional responses in metabolic tissues. Here we review advances in our understanding of the gene regulatory networks that contribute to hepatic responses to fasting and feeding. The advent of sequencing and -omics techniques have begun to facilitate a holistic understanding of the transcriptional landscape and its plasticity. We highlight transcription factors, their cofactors, and the pathways that they impact. We also discuss physiological factors that impinge on these responses, including circadian rhythms and sex differences. Finally, we review how dietary modifications modulate hepatic gene expression programs.

Redundant and non-redundant cytokine-activated enhancers control Csn1s2b expression in the lactating mouse mammary gland

Enhancers are transcription factor platforms that synergize with promoters to control gene expression. Here, we investigate enhancers that activate gene expression several hundred-fold exclusively in the lactating mouse mammary gland. Using ChIP-seq for activating histone marks and transcription factors, we identify two candidate enhancers and one super-enhancer in the Csn1s2b locus. Through experimental mouse genetics, we dissect the lactation-specific distal enhancer bound by the mammary-enriched transcription factors STAT5 and NFIB and the glucocorticoid receptor. While deletions of canonical binding motifs for NFIB and STAT5, individually or combined, have a limited biological impact, a non-canonical STAT5 site is essential for enhancer activity during lactation. In contrast, the intronic enhancer contributes to gene expression only in late pregnancy and early lactation, possibly by interacting with the distal enhancer. A downstream super-enhancer, which physically interacts with the distal enhancer, is required for the functional establishment of the Csn1s2b promoter and gene activation. Lastly, NFIB binding in the promoter region fine-tunes Csn1s2b expression. Our study provides comprehensive insight into the anatomy and biology of regulatory elements that employ the JAK/STAT signaling pathway and preferentially activate gene expression during lactation.

Enhancers and promoters work together to actively regulate gene expression affecting several biological processes. Here, the authors provide molecular insights into the regulation of enhancers and super-enhancers in the Csn1s2b locus during lactation.

OXR1A, a Coactivator of PRMT5 Regulating Histone Arginine Methylation

Oxidation resistance gene 1 (OXR1) protects cells against oxidative stress. We find that male mice with brain-specific isoform A knockout (Oxr1A^-/-) develop fatty liver. RNA sequencing of male Oxr1A^-/- liver indicates decreased growth hormone (GH) signaling, which is known to affect liver metabolism. Indeed, Gh expression is reduced in male mice Oxr1A^-/- pituitary gland and in rat Oxr1A^-/- pituitary adenoma cell-line GH3. Oxr1A^-/- male mice show reduced fasting-blood GH levels. Pull-down and proximity ligation assays reveal that OXR1A is associated with arginine methyl transferase PRMT5. OXR1A-depleted GH3 cells show reduced symmetrical dimethylation of histone H3 arginine 2 (H3R2me2s), a product of PRMT5 catalyzed methylation, and chromatin immunoprecipitation (ChIP) of H3R2me2s shows reduced Gh promoter enrichment. Finally, we demonstrate with purified proteins that OXR1A stimulates PRMT5/MEP50-catalyzed H3R2me2s. Our data suggest that OXR1A is a coactivator of PRMT5, regulating histone arginine methylation and thereby GH production within the pituitary gland.

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.

STAT5 Regulation of Sex-Dependent Hepatic CpG Methylation at Distal Regulatory Elements Mapping to Sex-Biased Genes

Growth hormone-activated STAT5b is an essential regulator of sex-differential gene expression in mouse liver; however, its impact on hepatic gene expression and epigenetic responses is poorly understood. Here, we found a substantial, albeit incomplete loss of liver sex bias in hepatocyte-specific STAT5a/STAT5b (collectively, STAT5)-deficient mouse liver. In male liver, many male-biased genes were downregulated in direct association with the loss of STAT5 binding; many female-biased genes, which show low STAT5 binding, were derepressed, indicating an indirect mechanism for repression by STAT5. Extensive changes in CpG methylation were seen in STAT5-deficient liver, where sex differences were abolished at 88% of ∼1,500 sex-differentially methylated regions, largely due to increased DNA methylation upon STAT5 loss. STAT5-dependent CpG hypomethylation was rarely found at proximal promoters of STAT5-dependent genes. Rather, STAT5 primarily regulated the methylation of distal enhancers, where STAT5 deficiency induced widespread hypermethylation at genomic regions enriched for accessible chromatin, enhancer histone marks (histone H3 lysine 4 monomethylation [H3K4me1] and histone H3 lysine 27 acetylation [H3K27ac]), STAT5 binding, and DNA motifs for STAT5 and other transcription factors implicated in liver sex differences. Thus, the sex-dependent binding of STAT5 to liver chromatin is closely linked to the sex-dependent demethylation of distal regulatory elements linked to STAT5-dependent genes important for liver sex bias.

Genetic factors contributing to extensive variability of sex-specific hepatic gene expression in Diversity Outbred mice

Sex-specific transcription characterizes hundreds of genes in mouse liver, many implicated in sex-differential drug and lipid metabolism and disease susceptibility. While the regulation of liver sex differences by growth hormone-activated STAT5 is well established, little is known about autosomal genetic factors regulating the sex-specific liver transcriptome. Here we show, using genotyping and expression data from a large population of Diversity Outbred mice, that genetic factors work in tandem with growth hormone to control the individual variability of hundreds of sex-biased genes, including many long non-coding RNA genes. Significant associations between single nucleotide polymorphisms and sex-specific gene expression were identified as expression quantitative trait loci (eQTLs), many of which showed strong sex-dependent associations. Remarkably, autosomal genetic modifiers of sex-specific genes were found to account for more than 200 instances of gain or loss of sex-specificity across eight Diversity Outbred mouse founder strains. Sex-biased STAT5 binding sites and open chromatin regions with strain-specific variants were significantly enriched at eQTL regions regulating correspondingly sex-specific genes, supporting the proposed functional regulatory nature of the eQTL regions identified. Binding of the male-biased, growth hormone-regulated repressor BCL6 was most highly enriched at trans-eQTL regions controlling female-specific genes. Co-regulated gene clusters defined by overlapping eQTLs included sets of highly correlated genes from different chromosomes, further supporting trans-eQTL action. These findings elucidate how an unexpectedly large number of autosomal factors work in tandem with growth hormone signaling pathways to regulate the individual variability associated with sex differences in liver metabolism and disease.

Growth Hormone Receptor Regulation in Cancer and Chronic Diseases

The GHR signaling pathway plays important roles in growth, metabolism, cell cycle control, immunity, homeostatic processes, and chemoresistance via both the JAK/STAT and the SRC pathways. Dysregulation of GHR signaling is associated with various diseases and chronic conditions such as acromegaly, cancer, aging, metabolic disease, fibroses, inflammation and autoimmunity. Numerous studies entailing the GHR signaling pathway have been conducted for various cancers. Diverse factors mediate the up- or down-regulation of GHR signaling through post-translational modifications. Of the numerous modifications, ubiquitination and deubiquitination are prominent events. Ubiquitination by E3 ligase attaches ubiquitins to target proteins and induces proteasomal degradation or starts the sequence of events that leads to endocytosis and lysosomal degradation. In this review, we discuss the role of first line effectors that act directly on the GHR at the cell surface including ADAM17, JAK2, SRC family member Lyn, Ubc13/CHIP, proteasome, βTrCP, CK2, STAT5b, and SOCS2. Activity of all, except JAK2, Lyn and STAT5b, counteract GHR signaling. Loss of their function increases the GH-induced signaling in favor of aging and certain chronic diseases, exemplified by increased lung cancer risk in case of a mutation in the SOCS2-GHR interaction site. Insight in their roles in GHR signaling can be applied for cancer and other therapeutic strategies.

Highly homologous mouse Cyp2a4 and Cyp2a5 genes are differentially expressed in the liver and both express long non-coding antisense RNAs

The mammalian Cytochrome P450 (Cyp) gene superfamily encodes enzymes involved in numerous metabolic pathways and are frequently expressed in the liver. Despite the remarkably high sequence similarity of Cyp2a4 and Cyp2a5 genes and their surrounding genomic regions, they exhibit differences in expression in the adult mouse liver. For example, Cyp2a4 is highly female-biased whereas Cyp2a5 is only moderately female-biased and Cyp2a4, but not Cyp2a5, is activated in liver cancer. We hypothesized that the limited sequence differences may help us identify the basis for this differential expression. An antisense expressed sequence tag had been uniquely annotated to the Cyp2a4 gene which led us to investigate this transcript as a possible regulator of this gene. We characterized the full-length antisense transcript and also discovered a similar transcript in the Cyp2a5 gene. These transcripts are nuclear long noncoding RNAs that are expressed similarly to their sense mRNA counterparts. This includes the sex-biased and liver tumor differences seen between the Cyp2a4 and Cyp2a5 genes, but we also find that these two genes and their antisense transcripts are expressed within different zones of the liver structure. Interestingly, while the differences in sex-biased expression of the mRNAs are established 1-2 months after birth, the antisense transcripts exhibit these expression differences earlier, at 3-4 weeks after birth. By analyzing published genomic data, we have identified candidate transcription factor binding sites that could account for differences in Cyp2a4/Cyp2a5 expression. Taken together, these studies characterize the first antisense RNAs within the Cyp supergene family and identify potential transcriptional and post-transcriptional mechanisms governing different Cyp2a4 and Cyp2a5 expression patterns in mouse liver.

Non-alcoholic Fatty Liver Disease as a Canonical Example of Metabolic Inflammatory-Based Liver Disease Showing a Sex-Specific Prevalence: Relevance of Estrogen Signaling

There is extensive evidence supporting the interplay between metabolism and immune response, that have evolved in close relationship, sharing regulatory molecules and signaling systems, to support biological functions. Nowadays, the disruption of this interaction in the context of obesity and overnutrition underlies the increasing incidence of many inflammatory-based metabolic diseases, even in a sex-specific fashion. During evolution, the interplay between metabolism and reproduction has reached a degree of complexity particularly high in female mammals, likely to ensure reproduction only under favorable conditions. Several factors may account for differences in the incidence and progression of inflammatory-based metabolic diseases between females and males, thus contributing to age-related disease development and difference in life expectancy between the two sexes. Among these factors, estrogens, acting mainly through Estrogen Receptors (ERs), have been reported to regulate several metabolic pathways and inflammatory processes particularly in the liver, the metabolic organ showing the highest degree of sexual dimorphism. This review aims to investigate on the interaction between metabolism and inflammation in the liver, focusing on the relevance of estrogen signaling in counteracting the development and progression of non-alcoholic fatty liver disease (NAFLD), a canonical example of metabolic inflammatory-based liver disease showing a sex-specific prevalence. Understanding the role of estrogens/ERs in the regulation of hepatic metabolism and inflammation may provide the basis for the development of sex-specific therapeutic strategies for the management of such an inflammatory-based metabolic disease and its cardio-metabolic consequences.

Impact of 3D genome organization, guided by cohesin and CTCF looping, on sex-biased chromatin interactions and gene expression in mouse liver

Several thousand sex-differential distal enhancers have been identified in mouse liver; however, their links to sex-biased genes and the impact of any sex-differences in nuclear organization and chromatin interactions are unknown. To address these issues, we first characterized 1847 mouse liver genomic regions showing significant sex differential occupancy by cohesin and CTCF, two key 3D nuclear organizing factors. These sex-differential binding sites were primarily distal to sex-biased genes but rarely generated sex-differential TAD (topologically associating domain) or intra-TAD loop anchors, and were sometimes found in TADs without sex-biased genes. A substantial subset of sex-biased cohesin-non-CTCF binding sites, but not sex-biased cohesin-and-CTCF binding sites, overlapped sex-biased enhancers. Cohesin depletion reduced the expression of male-biased genes with distal, but not proximal, sex-biased enhancers by >10-fold, implicating cohesin in long-range enhancer interactions regulating sex-biased genes. Using circularized chromosome conformation capture-based sequencing (4C-seq), we showed that sex differences in distal sex-biased enhancer-promoter interactions are common. Intra-TAD loops with sex-independent cohesin-and-CTCF anchors conferred sex specificity to chromatin interactions indirectly, by insulating sex-biased enhancer-promoter contacts and by bringing sex-biased genes into closer proximity to sex-biased enhancers. Furthermore, sex-differential chromatin interactions involving sex-biased gene promoters, enhancers, and lncRNAs were associated with sex-biased binding of cohesin and/or CTCF. These studies elucidate how 3D genome organization impacts sex-biased gene expression in a non-reproductive tissue through both direct and indirect effects of cohesin and CTCF looping on distal enhancer interactions with sex-differentially expressed genes.

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.

Sex-differences in the expression of genes in liver have a direct impact on liver diseases whose incidence and severity is sex-biased, and is controlled by hormones that regulate chemical alterations to histone proteins used to package chromosomal DNA. However, a direct demonstration of the functional importance of such sex differences in histone protein modifications has been elusive. Here, we address this question using a mouse model deficient in two enzymes, Ezh1/Ezh2, which generate the histone repressive mark H3K27me3. Remarkably, although H3K27me3 marks are formed by Ezh1/Ezh2 throughout the genome, loss of liver Ezh1/Ezh2 preferentially disrupts the control of sex-biased genes, with expression increasing in male mouse liver for many female-biased genes and decreasing for many male-biased genes. Sex-biased H3K27me3 repressive marks were abolished, and there was a gain of active histone marks at gene enhancers. We also found increased expression of many genes associated with liver fibrosis and hepatocellular carcinoma, which may help explain the increased sensitivity of Ezh1/Ezh2-deficient livers to hepatotoxic chemicals whose exposure may lead to sex differences in liver disease incidence and susceptibility. Thus, our findings highlight the potential role of sex differences in histone modifications catalyzed by Ezh1/Ezh2 in widespread sex differences in liver diseases.

Regulation of gene expression by growth hormone

Growth hormone (GH) plays a pivotal role in many physiological processes in humans, and in other mammalian and non-mammalian vertebrate species, through actions on somatic growth, tissue development and repair, and intermediary metabolism. This review will focus on mechanisms of GH actions on gene expression, primarily from the perspective of the genes that encode proteins stimulated by GH to regulate somatic growth, especially insulin-like growth factor 1 (IGF-I), but also others that are induced or repressed by GH. Topics to be discussed will include a brief overview of GH-mediated signal transduction pathways and how these cascades alter the functions of responsive transcription factors, with a specific focus on STAT5B, a key member of the signal transducers and activators of transcription family, characterization of essential GH-regulated genes, and elucidation of mechanisms of their regulation from biochemical, genetic, and genomic perspectives.

Growth Hormone Neuroprotection Against Kainate Excitotoxicity in the Retina is Mediated by Notch/PTEN/Akt Signaling

Purpose

In the retina, growth hormone (GH) promotes axonal growth, synaptic restoration, and protective actions against excitotoxicity. Notch signaling pathway is critical for neural development and participates in the retinal neuroregenerative process. We investigated the interaction of GH with Notch signaling pathway during its neuroprotective effect against excitotoxic damage in the chicken retina.

Methods

Kainate (KA) was used as excitotoxic agent and changes in the mRNA expression of several signaling markers were determined by qPCR. Also, changes in phosphorylation and immunoreactivity were determined by Western blotting. Histology and immunohistochemistry were performed for morphometric analysis. Overexpression of GH was performed in the quail neuroretinal-derived immortalized cell line (QNR/D) cell line. Exogenous GH was administered to retinal primary cell cultures to study the activation of signaling pathways.

Results

KA disrupted the retinal cytoarchitecture and induced significant cell loss in several retinal layers, but the coaddition of GH effectively prevented these adverse effects. We showed that GH upregulates the Notch signaling pathway during neuroprotection leading to phosphorylation of the PI3K/Akt signaling pathways through downregulation of PTEN. In contrast, cotreatment of GH with the Notch signaling inhibitor, DAPT, prevented its neuroprotective effect against KA. We identified binding sites in Notch1 and Notch2 genes for STAT5. Also, GH prevented Müller cell transdifferentiation and downregulated Sox2, FGF2, and PCNA after cotreatment with KA. Additionally, GH modified TNF receptors immunoreactivity suggesting anti-inflammatory actions.

Conclusions

Our data indicate that the neuroprotective effects of GH against KA injury in the retina are mediated through the regulation of Notch signaling. Additionally, anti-inflammatory and antiproliferative effects were observed.

Hepatic Tumor Formation in Adult Mice Developmentally Exposed to Organotin

BACKGROUND

Tributyltin (TBT) is a persistent and bioaccumulative environmental toxicant. Developmental exposure to TBT has been shown to cause fatty liver disease (steatosis), as well as increased adiposity in many species, leading to its characterization as an obesogen.

OBJECTIVE

We aimed to determine the long-term effects of developmental TBT exposure on the liver.

METHODS

C57BL/6J mice were exposed to a dose of TBT (0.5 mg / kg body weight per day; 3.07 μ M ) below the current developmental no observed adverse effect level (NOAEL) via drinking water, or drinking water alone, provided to the dam from preconception through lactation. Sires were exposed during breeding and lactation. Pups from two parity cycles were included in this study. Animals were followed longitudinally, and livers of offspring were analyzed by pathological evaluation, immunohistochemistry, immunoblotting, and RNA sequencing.

RESULTS

Developmental exposure to TBT led to increased adiposity and hepatic steatosis at 14 and 20 weeks of age and increased liver adenomas at 45 weeks of age in male offspring. Female offspring displayed increased adiposity as compared with males, but TBT did not lead to an increase in fatty liver or tumor development in female offspring. Liver tumors in male mice were enriched in pathways and gene signatures associated with human and rodent nonalcoholic fatty liver disease (NAFLD) and hepatocellular carcinoma (HCC). This includes down-regulation of growth hormone receptor (GHR) and of STAT5 signaling, which occurred in response to TBT exposure and preceded liver tumor development.

CONCLUSIONS

These data reveal a previously unappreciated ability of TBT to increase risk for liver tumorigenesis in mice in a sex-specific manner. Taken together, these findings provide new insights into how early life environmental exposures contribute to liver disease in adulthood. https://doi.org/10.1289/EHP5414.

Loss of Transcriptional Repression by BCL6 Confers Insulin Sensitivity in the Setting of Obesity

Accumulation of visceral adiposity is directly linked to the morbidity of obesity, while subcutaneous body fat is considered more benign. We have identified an unexpected role for B cell lymphoma 6 (BCL6), a critical regulator of immunity, in the developmental expansion of subcutaneous adipose tissue. In adipocyte-specific knockout mice (Bcl6^AKO), we found that Bcl6 deletion results in strikingly increased inguinal, but not perigonadal, adipocyte size and tissue mass in addition to marked insulin sensitivity. Genome-wide RNA expression and DNA binding analyses revealed that BCL6 controls gene networks involved in cell growth and fatty acid biosynthesis. Using deuterium label incorporation and comprehensive adipokine and lipid profiling, we discovered that ablation of adipocyte Bcl6 enhances subcutaneous adipocyte lipogenesis, increases levels of adiponectin and fatty acid esters of hydroxy fatty acids (FAHFAs), and prevents steatosis. Thus, our studies identify BCL6 as a negative regulator of subcutaneous adipose tissue expansion and metabolic health.

Senagolage et al. identify BCL6 as a key regulator of body fat distribution. BCL6 directly represses fatty acid biosynthetic and growth genes in adipocytes. Mice constitutively lacking adipocyte Bcl6 exhibit expansion of their subcutaneous adipose tissue, enhanced insulin sensitivity, and protection from hepatic steatosis.

Octopus-toolkit: a workflow to automate mining of public epigenomic and transcriptomic next-generation sequencing data

Octopus-toolkit is a stand-alone application for retrieving and processing large sets of next-generation sequencing (NGS) data with a single step. Octopus-toolkit is an automated set-up-and-analysis pipeline utilizing the Aspera, SRA Toolkit, FastQC, Trimmomatic, HISAT2, STAR, Samtools, and HOMER applications. All the applications are installed on the user's computer when the program starts. Upon the installation, it can automatically retrieve original files of various epigenomic and transcriptomic data sets, including ChIP-seq, ATAC-seq, DNase-seq, MeDIP-seq, MNase-seq and RNA-seq, from the gene expression omnibus data repository. The downloaded files can then be sequentially processed to generate BAM and BigWig files, which are used for advanced analyses and visualization. Currently, it can process NGS data from popular model genomes such as, human (Homo sapiens), mouse (Mus musculus), dog (Canis lupus familiaris), plant (Arabidopsis thaliana), zebrafish (Danio rerio), fruit fly (Drosophila melanogaster), worm (Caenorhabditis elegans), and budding yeast (Saccharomyces cerevisiae) genomes. With the processed files from Octopus-toolkit, the meta-analysis of various data sets, motif searches for DNA-binding proteins, and the identification of differentially expressed genes and/or protein-binding sites can be easily conducted with few commands by users. Overall, Octopus-toolkit facilitates the systematic and integrative analysis of available epigenomic and transcriptomic NGS big data.

Widespread epigenetic changes to the enhancer landscape of mouse liver induced by a specific xenobiotic agonist ligand of the nuclear receptor CAR

CAR (Nr1i3), a liver nuclear receptor and xenobiotic sensor, induces drug, steroid and lipid metabolism and dysregulates genes linked to hepatocellular carcinogenesis, but its impact on the liver epigenome is poorly understood. TCPOBOP, a halogenated xenochemical and highly specific CAR agonist ligand, induces localized chromatin opening or closing at several thousand mouse liver genomic regions, discovered as differential DNase-hypersensitive sites (ΔDHS). Active enhancer and promoter histone marks induced by TCPOBOP were enriched at opening DHS and TCPOBOP-inducible genes. Enrichment of CAR binding and CAR motifs was seen at opening DHS and their inducible drug/lipid metabolism gene targets, and at many constitutively open DHS located nearby. TCPOBOP-responsive cell cycle and DNA replication genes co-dependent on MET/EGFR signaling for induction were also enriched for CAR binding. A subset of opening DHS and many closing DHS mapping to TCPOBOP-responsive target genes did not bind CAR, indicating an indirect mechanism for their changes in chromatin accessibility. TCPOBOP-responsive DHS were also enriched for induced binding of RXRA, CEBPA and CEBPB, and for motifs for liver-enriched factors that may contribute to liver-specific transcriptional responses to TCPOBOP exposure. These studies elucidate the enhancer landscape of TCPOBOP-exposed liver and the widespread epigenetic changes that are induced by both direct and indirect mechanisms linked to CAR activation. The global maps of thousands of environmental chemical-induced epigenetic changes described here constitute a rich resource for further research on xenochemical effects on liver chromatin states and the epigenome.

Sex-Biased lncRNAs Inversely Correlate With Sex-Opposite Gene Coexpression Networks in Diversity Outbred Mouse Liver

Sex differences in liver gene expression are determined by pituitary growth hormone secretion patterns, which regulate sex-dependent liver transcription factors and establish sex-specific chromatin states. Hypophysectomy (hypox) identifies two major classes of liver sex-biased genes, defined by their sex-dependent positive or negative responses to pituitary hormone ablation. However, the mechanisms that underlie each hypox-response class are unknown. We sought to discover candidate, regulatory, long noncoding RNAs (lncRNAs) controlling responsiveness to hypox. We characterized gene structures and expression patterns for 15,558 mouse liver-expressed lncRNAs, including many sex-specific lncRNAs regulated during postnatal development or subject to circadian regulation. Using the high natural allelic variance of Diversity Outbred (DO) mice, we discovered tightly coexpressed clusters of sex-specific protein-coding genes (gene modules) in male and female DO liver. Remarkably, many gene modules were strongly enriched for sex-specific genes within a single hypox-response class, indicating that the genetic heterogeneity of DO mice encompasses responsiveness to hypox. Moreover, several distant gene modules were enriched for gene subsets of the same hypox-response class, highlighting the complex regulation of hypox-responsiveness. Finally, we identified eight sex-specific lncRNAs with strong negative regulatory potential, as indicated by their strong negative correlation of expression across DO mouse livers with that of protein-coding gene modules enriched for genes of the opposite sex bias and inverse hypox-response class. These findings reveal an important role for genetic factors in regulating responsiveness to hypox, and present testable hypotheses for the roles of sex-biased liver lncRNAs in controlling the sex-bias of liver gene expression.

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.