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

Genome-wide identification of STATs and analysis of their role in sex determination in Pacific oysters (Crassostrea gigas)

STAT (signal transducer and activator of the transcription) proteins, are a group of highly conserved transcription factors and fundamental components of the JAK-STAT signaling pathway. They play crucial roles in a variety of biological processes, such as immunity, proliferation, differentiation, and growth. However, little information is known regarding their role in gonad development and sex determination in mollusks. In this study, we identified 3 STAT genes in Pacific Oyster Crassostrea gigas. Phylogenetic analysis showed that STATs from mollusks were highly conserved, and most of them had four identical motif regions, except for the STAT1 and STAT3 predicted sequences from Crassostrea hongkongensis. Tissue expression analysis indicated CgSTAT1 had a high expression level in most tissues, while CgSTAT3 had a low expression level in most tissues. Expression analysis of early developmental stages showed CgSTAT1 had a higher expression level from egg to D shaped larva and a lower expression level in subsequent stages. In contrast CgSTAT1, CgSTAT2 had a reverse expression pattern. Expression analysis of different developmental stages of diploid gonads indicated that CgSTAT1 had a higher expression level at the S1 and S3 stages relative to the S2 stage in females, while in males the S3 stage had a higher expression than than the S2 stage. The expression level of CgSTAT1 between diploids and triploids in females differed significantly, but there were no significant differences in males. Expression of CgSTAT2 differed significantly between diploid and triploid males. These data suggest an important role for STATs in sex differentiation in diploid and triploid oysters. Our study is the first to explore the role of STATs in sex differentiation and gonadal development in oysters, and will help us better understand the molecular mechanisms of sex differentiation in shellfish.

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

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

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

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

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

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

Assessing the Differential Methylation Analysis Quality for Microarray and NGS Platforms

Differential methylation (DM) is actively recruited in different types of fundamental and translational studies. Currently, microarray- and NGS-based approaches for methylation analysis are the most widely used with multiple statistical models designed to extract differential methylation signatures. The benchmarking of DM models is challenging due to the absence of gold standard data. In this study, we analyze an extensive number of publicly available NGS and microarray datasets with divergent and widely utilized statistical models and apply the recently suggested and validated rank-statistic-based approach Hobotnica to evaluate the quality of their results. Overall, microarray-based methods demonstrate more robust and convergent results, while NGS-based models are highly dissimilar. Tests on the simulated NGS data tend to overestimate the quality of the DM methods and therefore are recommended for use with caution. Evaluation of the top 10 DMC and top 100 DMC in addition to the not-subset signature also shows more stable results for microarray data. Summing up, given the observed heterogeneity in NGS methylation data, the evaluation of newly generated methylation signatures is a crucial step in DM analysis. The Hobotnica metric is coordinated with previously developed quality metrics and provides a robust, sensitive, and informative estimation of methods' performance and DM signatures' quality in the absence of gold standard data solving a long-existing problem in DM analysis.

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.

Genetic variation in the Y chromosome and sex-biased DNA methylation in somatic cells in the mouse

Several lines of evidence suggest that the presence of the Y chromosome influences DNA methylation of autosomal loci. To better understand the impact of the Y chromosome on autosomal DNA methylation patterns and its contribution to sex bias in methylation, we identified Y chromosome dependent differentially methylated regions (yDMRs) using whole-genome bisulfite sequencing methylation data from livers of mice with different combinations of sex-chromosome complement and gonadal sex. Nearly 90% of the autosomal yDMRs mapped to transposable elements (TEs) and most of them had lower methylation in XY compared to XX or XO mice. Follow-up analyses of four reporter autosomal yDMRs showed that Y-dependent methylation levels were consistent across most somatic tissues but varied in strains with different origins of the Y chromosome, suggesting that genetic variation in the Y chromosome influenced methylation levels of autosomal regions. Mice lacking the q-arm of the Y chromosome (B6.NPYq-2) as well as mice with a loss-of-function mutation in Kdm5d showed no differences in methylation levels compared to wild type mice. In conclusion, the Y-linked modifier of TE methylation is likely to reside on the short arm of Y chromosome and further studies are required to identify this gene.

Identification of mutations in porcine STAT5A that contributes to the transcription of CISH

Identification of causative genes or genetic variants associated with phenotype traits benefits the genetic improvement of animals. CISH plays a role in immunity and growth, however, the upstream transcriptional factors of porcine CISH and the genetic variations in these factors remain unclear. In this study, we firstly identified the minimal core promoter of porcine CISH and confirmed the existence of STATx binding sites. Overexpression and RT-qPCR demonstrated STAT5A increased CISH transcriptional activity (P < 0.01) and mRNA expression (P < 0.01), while GATA1 inhibited CISH transcriptional activity (P < 0.01) and the following mRNA expression (P < 0.05 or P < 0.01). Then, the putative functional genetic variations of porcine STAT5A were screened and a PCR-SSCP was established for genotype g.508A>C and g.566C>T. Population genetic analysis showed the A allele frequency of g.508A>C and C allele frequency of g.566C>T was 0.61 and 0.94 in Min pigs, respectively, while these two alleles were fixed in the Landrace population. Statistical analysis showed that Min piglets with CC genotype at g.566C>T or Hap1: AC had higher 28-day body weight, 35-day body weight, and ADG than TC or Hap3: CT animals (P < 0.05, P < 0.05). Further luciferase activity assay demonstrated that the activity of g.508A>C in the C allele was lower than the A allele (P < 0.05). Collectively, the present study demonstrated that STAT5A positively regulated porcine CISH transcription, and SNP g.566C>T in the STAT5A was associated with the Min piglet growth trait.

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.

Growth Hormone and Insulin-Like Growth Factor 1 Regulation of Nonalcoholic Fatty Liver Disease

Patients with obesity have a high prevalence of nonalcoholic fatty liver disease (NAFLD), representing a spectrum of simple steatosis to nonalcoholic steatohepatitis (NASH), without and with fibrosis. Understanding the etiology of NAFLD is clinically relevant since NAFLD is an independent risk factor for diabetes and cardiovascular disease. In addition, NASH predisposes patients to the development of cirrhosis and hepatocellular carcinoma, and NASH cirrhosis represents the fastest growing indication for liver transplantation in the United States. It is appreciated that multiple factors are involved in the development and progression of NAFLD. Growth hormone (GH) and insulin-like growth factor 1 (IGF1) regulate metabolic, immune, and hepatic stellate cell function, and alterations in the production and function of GH is associated with obesity and NAFLD/NASH. Therefore, this review will focus on the potential role of GH and IGF1 in the regulation of hepatic steatosis, inflammation, and fibrosis.

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.

Sex differences in the intergenerational inheritance of metabolic traits

Strong evidence suggests that early-life exposures to suboptimal environmental factors, including those in utero, influence our long-term metabolic health. This has been termed developmental programming. Mounting evidence suggests that the growth and metabolism of male and female fetuses differ. Therefore, sexual dimorphism in response to pre-conception or early-life exposures could contribute to known sex differences in susceptibility to poor metabolic health in adulthood. However, until recently, many studies, especially those in animal models, focused on a single sex, or, often in the case of studies performed during intrauterine development, did not report the sex of the animal at all. In this review, we (a) summarize the evidence that male and females respond differently to a suboptimal pre-conceptional or in utero environment, (b) explore the potential biological mechanisms that underlie these differences and (c) review the consequences of these differences for long-term metabolic health, including that of subsequent generations.

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.

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.

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.

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.