Sex-dependent pituitary growth hormone (GH) secretory patterns determine the sex-biased expression of >1,000 genes in mouse and rat liver, affecting lipid and drug metabolism, inflammation, and disease. A fundamental biological question is how robust differential expression can be achieved for hundreds of sex-biased genes simply based on the GH input signal pattern: pulsatile GH stimulation in males versus near-continuous GH exposure in females. STAT5 is an essential transcriptional mediator of the sex-dependent effects of GH in the liver, but the mechanisms that underlie its sex-dependent actions are obscure. Here we elucidate the dynamic, sex-dependent binding of STAT5 and the GH/STAT5-regulated repressor BCL6 to mouse liver chromatin genome wide, revealing a counteractive interplay between these two regulators of sex differences in liver gene expression. Our findings establish a close correlation between sex-dependent STAT5 binding and sex-biased target gene expression. Moreover, sex-dependent STAT5 binding correlated positively with sex-biased DNase hypersensitivity and H3-K4me1 and H3-K4me3 (activating) marks, correlated negatively with sex-biased H3-K27me3 (repressive) marks, and was associated with sexdifferentially enriched motifs for HNF6/CDP factors. Importantly, BCL6 binding was preferentially associated with repression of female-biased STAT5 targets in male liver. Furthermore, BCL6 and STAT5 common targets but not BCL6 unique targets showed strong enrichment for lipid and drug metabolism. These findings provide a comprehensive, genome-wide view of the mechanisms whereby these two GH-regulated transcription factors establish and maintain sex differences affecting liver physiology and disease. The approaches used here to characterize sex-dependent STAT5 and BCL6 binding can be applied to other condition-specific regulatory factors and binding sites and their interplay with cooperative chromatin binding factors.

Targeted disruption of the signal transducer and activator of transcription 5b gene (STAT5b) leads to decreased expression in male mouse liver of a male-predominant cytochrome (Cyp) 2d protein, whereas female-predominant Cyp2b proteins are increased. Presently, we characterize the effects of STAT5b deficiency on 15 specific, individual Cyp RNAs and other sexually dimorphic liver gene products. All seven male-specific RNAs investigated were decreased to normal female levels in STAT5b-deficient male liver, whereas five of eight female-specific RNAs, designated class I female genes, were increased in expression up to 200-fold or more. STAT5b deficiency had a much more modest effect on the expression of these genes in females. Hypophysectomy and GH replacement studies demonstrated positive GH pulse regulation of all seven male RNAs and negative GH pulse regulation of class I, but not class II, female RNAs in wild-type, but not in STAT5b-deficient, male mice. A majority of the sex-specific genes responded in parallel to the loss of STAT5b and the loss of hepatocyte nuclear factor 4alpha, indicating that both transcription factors are essential and suggesting they may coregulate sexually dimorphic liver gene expression. Continuous GH treatment of intact male mice, which overrides the endogenous male, pulsatile plasma GH pattern, down-regulated all seven male RNAs and induced expression of the five class I female RNAs within 4-7 d; however, induction of class II female RNAs was delayed until d 7-14. Given the slow responses of all 15 genes to changes in plasma GH status, GH regulation of sex-specific Cyp expression is proposed to be indirect and mediated by STAT5b- and hepatocyte nuclear factor 4alpha-dependent factors that may include repressors of female-specific Cyps and other targets of GH action.

The transcriptional repressor Bcl6 is a male-specific rat liver gene product and one of 24 early GH-response genes encoding DNA-binding proteins. Presently, the sex specificity of Bcl6 was shown to emerge at puberty, when hepatic Bcl6 mRNA was induced in males and repressed in females by the female plasma GH profile. Hepatic Bcl6 mRNA was increased to near-normal male levels in hypophysectomized females and was extinguished in intact males given a continuous GH infusion (female-like GH pattern). Bcl6 was also repressed in adult male somatostatin-deficient mice, where plasma GH profiles are female like. Hepatic Bcl6 RNA was rapidly down-regulated by GH pulse treatment, both in hypophysectomized male rats and in primary rat hepatocytes. Bcl6 was substantially induced in female mice deficient in hepatic signal transducer and activator of transcription (STAT)5a/STAT5b, suggesting that these STAT transcriptional mediators of GH signaling repress Bcl6. Indeed, STAT5 was bound to Bcl6 STAT5-binding region-B, previously associated with Bcl6 repression, in both male and female liver chromatin. STAT5 also bound to Bcl6 region-A in male chromatin but only during a plasma GH pulse. Analysis of primary transcripts (heterogeneous nuclear RNA) across the Bcl6 gene revealed a novel mechanism of GH-dependent sex specificity, with two apparent blocks in Bcl6 transcription elongation seen in female liver and in continuous GH-treated male liver, one early in intron 4 and one in exon 5, which together reduced transcription beyond exon 5 more than 300-fold. Finally, Bcl6 was bound to a subset of STAT5-binding sites in male liver chromatin, including a Socs2 STAT5-binding site where Bcl6 binding increased substantially between plasma GH pulses, i.e. when STAT5 binding was low. Bcl6 and STAT5 binding are thus inversely coordinated by the endogenous pulses of pituitary GH release, suggesting this male-specific transcriptional repressor modulates hepatic GH signaling to select STAT5 target genes.

Hepatocyte-specific, albumin-Cre recombinase-mediated deletion of the entire mouse Stat5a-Stat5b locus was carried out to evaluate the role of signal transducer and activator of transcription 5a and 5b (STAT5ab) in the sex-dependent transcriptional actions of GH in the liver. The resultant hepatocyte STAT5ab-deficient mice were fertile, and unlike global STAT5b-deficient male mice, postnatal body weight gain was normal, despite a 50% decrease in serum IGF-I. Whole-liver STAT5ab RNA decreased by approximately 65-85%, and residual STAT5 immunostaining was observed in a minority of the hepatocytes, indicating incomplete excision by Cre-recombinase. Quantitative PCR analysis of 20 sexually dimorphic, liver-expressed genes revealed significant down-regulation of 10 of 11 male-specific genes in livers of male hepatocyte STAT5ab-deficient mice. Class I female-specific liver genes were markedly up-regulated (de-repressed), whereas the expression of class II female genes, belonging to the Cyp3a subfamily, was unaffected by the loss of hepatocyte STAT5ab. STAT5ab is thus required in the liver for positive regulation of male-specific genes and for negative regulation of a subset of female-specific genes. Continuous GH infusion strongly induced (>500-fold) the class II female gene Cyp3a16 in both wild-type and hepatocyte STAT5ab-deficient male mice, indicating sex-specific transcriptional regulation by GH that is STAT5ab independent. In contrast, hepatocyte STAT5ab deficiency abolished the strong suppression of the male-specific Cyp2d9 by continuous GH seen in control mouse liver. Analysis of global STAT5a-deficient mice indicated no essential requirement of STAT5a for expression of these sex-specific liver Cyp genes. Thus, the major loss of liver sexual dimorphism in hepatocyte STAT5ab-deficient mice can primarily be attributed to the loss of STAT5b.

Plasma GH profiles regulate the sexually dimorphic expression of cytochromes P450 and many other genes in rat and mouse liver; however, the proximal transcriptional regulators of these genes are unknown. Presently, we characterize three liver transcription factors that are expressed in adult female rat and mouse liver at levels up to 16-fold [thymus high-mobility group box protein (Tox)], 73-fold [tripartite motif-containing 24 (Trim24)/transcription initiation factor-1alpha (TIF1alpha)], and 125-fold [cut-like 2 (Cutl2)/cut homeobox 2 (Cux2)] higher than in adult males, depending on the strain and species, with Tox expression only detected in mice. In rats, these sex differences first emerged at puberty, when the high prepubertal expression of Cutl2 and Trim24 was extinguished in males but was further increased in females. Rat hepatic expression of Cutl2 and Trim24 was abolished by hypophysectomy and, in the case of Cutl2, was restored to near-female levels by continuous GH replacement. Cutl2 and Trim24 were increased to female-like levels in livers of intact male rats and mice treated with GH continuously (female GH pattern), whereas Tox expression reached only about 40% of adult female levels. Expression of all three genes was also elevated to normal female levels or higher in male mice whose plasma GH profile was feminized secondary to somatostatin gene disruption. Cutl2 and Trim24 both responded to GH infusion in mice within 10-24 h and Tox within 4 d, as compared with at least 4-7 d required for the induced expression of several continuous GH-regulated cytochromes P450 and other female-specific hepatic genes. Cutl2, Trim24, and Tox were substantially up-regulated in livers of male mice deficient in either of two transcription factors implicated in GH regulation of liver sex specificity, namely, signal transducer and activator of transcription 5b (STAT5b) and hepatocyte nuclear factor 4alpha (HNF4alpha), with sex-specific expression being substantially reduced or lost in mice deficient in either nuclear factor. Cutl2 and Trim24 both display transcriptional repressor activity and could thus contribute to the loss of GH-regulated, male-specific liver gene expression seen in male mice deficient in STAT5b or HNF4alpha. Binding sites for Cutl1, whose DNA-binding specificity is close to that of Cutl2, were statistically overrepresented in STAT5b-dependent male-specific mouse genes, lending support to this hypothesis.

Phylogenetic footprinting was used to predict functional transcription factor binding sites (TFBS) for signal transducer and activator of transcription (STAT) 5, a GH-activated transcription factor, in the GH-responsive genes IGF-I, SOCS2, and HNF6. Each gene, including upstream (100 kb) and downstream regions (25 kb), was aligned across four species and searched for conserved STAT5-binding sites using TFBS matrices. Predicted sites were classified as paired or single and whether or not they matched the STAT5 consensus sequence TTCN(3)GAA. Fifty-seven of the predicted genomic regions were assayed by chromatin immunoprecipitation from male rat liver with high STAT5 activity. STAT5 binding was enriched (up to 24-fold) at eight genomic regions of IGF-I, including three novel regions in the second intron, and at four regions of SOCS2, including three novel upstream sites. STAT5 binding to HNF6 was modestly enriched (up to 3-fold) at one consensus site and two novel, nonconsensus sites. Overall, 14 of 17 identified sites were paired STAT5 sites. STAT5 binding to these sites was dynamic in male rat liver, cycling on and off in response to each plasma GH pulse. Moreover, sex-specific STAT5 binding was apparent; in female rat liver, where nuclear STAT5 activity is generally low, STAT5 binding to IGF-I and SOCS2 was limited to high-affinity sites. Analysis of the verified STAT5 binding sites indicated that STAT5 TFBS matrix 459 in combination with a STAT5 consensus sequence was the best predictor of STAT5 binding to these three genes. Using these criteria, multiple novel STAT5 binding sites were identified and then verified in several other GH-inducible genes, including MUP genes, where male-specific gene expression was associated with male-specific STAT5 binding to multiple low-affinity STAT5 sites.

Many genes are expressed in mammalian liver in a sexually dimorphic manner. DNA microarray analysis has shown that growth hormone (GH) and its sex-dependent pattern of pituitary secretion play a major role in establishing the sexually dimorphic patterns of liver gene expression. However, GH may exert effects on protein post-translational modification and nuclear localization that are not reflected at the mRNA level. To investigate these potential effects of GH, we used two-dimensional gel electrophoresis followed by LC-MS/MS to: 1) identify rat liver nuclear proteins whose abundance or state of post-translational modification displays sex-dependent differences; and 2) determine the role of the plasma GH profile in establishing these differences. Nuclear extracts prepared from livers of individual male (n ‫؍‬ 9) and female (n ‫؍‬ 5) adult rats, and from males given GH by continuous infusion for 7 days to feminize liver gene expression (n ‫؍‬ 5 rats), were resolved by two-dimensional electrophoresis. Image analysis of SYPRO Ruby-stained gels revealed 165 sexually dimorphic protein spots that differ in normalized volume between male and female groups by >1.5-fold at p < 0.05. Sixty of these proteins exhibited female-like changes in spot abundance following continuous GH treatment. Comparison of male and GH-treated male groups revealed 130 proteins that displayed >1.5-fold differences in abundance, with 60 of these GH-responsive spots being sexually dimorphic. Thus, GH plays an important role in establishing the sex-dependent differences in liver nuclear protein content. Twenty-eight of the sexually dimorphic and/or GH-regulated protein spots were identified by LC-MS/MS. Proteins identified include regucalcin, nuclear factor 45, and heterogeneous nuclear ribonucleoproteins A3, D-like, and K, in addition to proteins such as GST, normally associated with cytosolic extracts but also reported to be localized in the nucleus.

No abstract