Cross Talk Between GH-Regulated Transcription Factors HNF6 and CUX2 in Adult Mouse Liver

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.

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.

CUT homeobox genes: transcriptional regulation of neuronal specification and beyond

CUT homeobox genes represent a captivating gene class fulfilling critical functions in the development and maintenance of multiple cell types across a wide range of organisms. They belong to the larger group of homeobox genes, which encode transcription factors responsible for regulating gene expression patterns during development. CUT homeobox genes exhibit two distinct and conserved DNA binding domains, a homeodomain accompanied by one or more CUT domains. Numerous studies have shown the involvement of CUT homeobox genes in diverse developmental processes such as body axis formation, organogenesis, tissue patterning and neuronal specification. They govern these processes by exerting control over gene expression through their transcriptional regulatory activities, which they accomplish by a combination of classic and unconventional interactions with the DNA. Intriguingly, apart from their roles as transcriptional regulators, they also serve as accessory factors in DNA repair pathways through protein-protein interactions. They are highly conserved across species, highlighting their fundamental importance in developmental biology. Remarkably, evolutionary analysis has revealed that CUT homeobox genes have experienced an extraordinary degree of rearrangements and diversification compared to other classes of homeobox genes, including the emergence of a novel gene family in vertebrates. Investigating the functions and regulatory networks of CUT homeobox genes provides significant understanding into the molecular mechanisms underlying embryonic development and tissue homeostasis. Furthermore, aberrant expression or mutations in CUT homeobox genes have been associated with various human diseases, highlighting their relevance beyond developmental processes. This review will overview the well known roles of CUT homeobox genes in nervous system development, as well as their functions in other tissues across phylogeny.

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

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

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

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

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.

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.

Mutations and variants of ONECUT1 in diabetes

Genes involved in distinct diabetes types suggest shared disease mechanisms. Here we show that One Cut Homeobox 1 (ONECUT1) mutations cause monogenic recessive syndromic diabetes in two unrelated patients, characterized by intrauterine growth retardation, pancreas hypoplasia and gallbladder agenesis/hypoplasia, and early-onset diabetes in heterozygous relatives. Heterozygous carriers of rare coding variants of ONECUT1 define a distinctive subgroup of diabetic patients with early-onset, nonautoimmune diabetes, who respond well to diabetes treatment. In addition, common regulatory ONECUT1 variants are associated with multifactorial type 2 diabetes. Directed differentiation of human pluripotent stem cells revealed that loss of ONECUT1 impairs pancreatic progenitor formation and a subsequent endocrine program. Loss of ONECUT1 altered transcription factor binding and enhancer activity and NKX2.2/NKX6.1 expression in pancreatic progenitor cells. Collectively, we demonstrate that ONECUT1 controls a transcriptional and epigenetic machinery regulating endocrine development, involved in a spectrum of diabetes, encompassing monogenic (recessive and dominant) as well as multifactorial inheritance. Our findings highlight the broad contribution of ONECUT1 in diabetes pathogenesis, marking an important step toward precision diabetes medicine.

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.

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.

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.

Circadian Clock Regulation of Hepatic Energy Metabolism Regulatory Circuits

The liver is a critical organ of energy metabolism. At least 10% of the liver transcriptome demonstrates rhythmic expression, implying that the circadian clock regulates large programmes of hepatic genes. Here, we review the mechanisms by which this rhythmic regulation is conferred, with a particular focus on the transcription factors whose actions combine to impart liver- and time-specificity to metabolic gene expression.

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.

The Mouse Microbiome Is Required for Sex-Specific Diurnal Rhythms of Gene Expression and Metabolism

The circadian clock and associated feeding rhythms have a profound impact on metabolism and the gut microbiome. To what extent microbiota reciprocally affect daily rhythms of physiology in the host remains elusive. Here, we analyzed transcriptome and metabolome profiles of male and female germ-free mice. While mRNA expression of circadian clock genes revealed subtle changes in liver, intestine, and white adipose tissue, germ-free mice showed considerably altered expression of genes associated with rhythmic physiology. Strikingly, the absence of the microbiome attenuated liver sexual dimorphism and sex-specific rhythmicity. The resulting feminization of male and masculinization of female germ-free animals is likely caused by altered sexual development and growth hormone secretion, associated with differential activation of xenobiotic receptors. This defines a novel mechanism by which the microbiome regulates host metabolism.

CUX2 functions as an oncogene in papillary thyroid cancer

BACKGROUND

In recent years, the incidence of thyroid cancer (TC), the most common endocrine malignancy, has been increasing. Emerging evidence indicates that the CUT/CUX/CDP family of proteins can play an important role in tumor development and progression by regulating many cancer-related functions. However, the molecular functions of CUX2 in TC remain unknown.

METHODS

In this study, we used a series of loss-of-function experiments and Western blot analysis to investigate the function of CUX2 in TC and the mechanisms involved.

RESULTS

Our data revealed that CUX2 expression levels were upregulated in papillary thyroid cancer (PTC). Functionally, CUX2 silencing significantly inhibited PTC cell line (KTC-1 and BCPAP) proliferation, colony formation, migration, invasion, and apoptosis. Furthermore, CUX2 induced epithelial-mesenchymal transition (EMT) and influenced the phosphorylation of AKT and mTOR in the PI3K-AKT-mTOR pathways.

CONCLUSION

In summary, CUX2 may function as a tumor promoter in TC.

Stat5b Regulates Sexually Dimorphic Gene Expression in Zebrafish Liver

Sexual size dimorphism is an interesting phenomenon occurred in many fish species. Wildtype zebrafish exhibits a significant sexual dimorphism in body size at the adult stage. Previous studies indicated that sexual size dimorphism was eliminated in stat5b-mutated zebrafish. Herein, the comparative transcriptome analysis was conducted to observe the genes and pathways involved in sexual size dimorphism. The number of male-biased and female-biased genes was much less in the liver of stat5b mutant zebrafish than in wildtype. Gene ontology (GO) enrichment and Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis indicated that multiple pathways related to metabolism were affected upon loss of stat5b function. qRT-PCR results also validated that sexually dimorphic expression of a set of genes was lost when stat5b was mutated. Furthermore, the weighted correlation network analysis (WGCNA) detected many candidate genes related to the growth traits and stat5b function, such as greb1, lepr, and igf2b. Our data suggest that stat5b should regulate the sexually dimorphic gene expression in zebrafish liver and add in understanding of the molecular mechanisms underlying sexual size dimorphism in fish species.

The crux of Cux genes in neuronal function and plasticity

Cux1 and Cux2 are the vertebrate members of a family of homeodomain transcription factors (TF) containing Cut repeat DNA-binding sequences. Perturbation of their expression has been implicated in a wide variety of diseases and disorders, ranging from cancer to autism spectrum disorder (ASD). Within the nervous system, both genes are expressed during neurogenesis and in specific neuronal subpopulations. Their role during development and circuit specification is discussed here, with a particular focus on the cortex where their restricted expression in pyramidal neurons of the upper layers appears to be responsible for many of the specialized functions of these cells, and where their functions have been extensively investigated. Finally, we discuss how Cux TF represent a promising avenue for manipulating neuronal function and for reprogramming.

Metabolic and molecular framework for the enhancement of endurance by intermittent food deprivation

Evolutionary considerations suggest that the body has been optimized to perform at a high level in the food-deprived state when fatty acids and their ketone metabolites are a major fuel source for muscle cells. Because controlled food deprivation in laboratory animals and intermittent energy restriction in humans is a potent physiologic stimulus for ketosis, we designed a study to determine the impact of intermittent food deprivation during endurance training on performance and to elucidate the underlying cellular and molecular mechanisms. Male mice were randomly assigned to either ad libitum feeding or alternate-day food deprivation (ADF) groups, and half of the mice in each diet group were trained daily on a treadmill for 1 mo. A run to exhaustion endurance test performed at the end of the training period revealed superior performance in the mice maintained on ADF during training compared to mice fed ad libitum during training. Maximal O₂ consumption was increased similarly by treadmill training in mice on ADF or ad libitum diets, whereas respiratory exchange ratio was reduced in ADF mice on food-deprivation days and during running. Analyses of gene expression in liver and soleus tissues, and metabolomics analysis of blood suggest that the metabolic switch invoked by ADF and potentiated by exercise strongly modulates molecular pathways involved in mitochondrial biogenesis, metabolism, and cellular plasticity. Our findings demonstrate that ADF engages metabolic and cellular signaling pathways that result in increased metabolic efficiency and endurance capacity.-Marosi, K., Moehl, K., Navas-Enamorado, I., Mitchell, S. J., Zhang, Y., Lehrmann, E., Aon, M. A., Cortassa, S., Becker, K. G., Mattson, M. P. Metabolic and molecular framework for the enhancement of endurance by intermittent food deprivation.

Transcription factor activity rhythms and tissue-specific chromatin interactions explain circadian gene expression across organs

Temporal control of physiology requires the interplay between gene networks involved in daily timekeeping and tissue function across different organs. How the circadian clock interweaves with tissue-specific transcriptional programs is poorly understood. Here, we dissected temporal and tissue-specific regulation at multiple gene regulatory layers by examining mouse tissues with an intact or disrupted clock over time. Integrated analysis uncovered two distinct regulatory modes underlying tissue-specific rhythms: tissue-specific oscillations in transcription factor (TF) activity, which were linked to feeding-fasting cycles in liver and sodium homeostasis in kidney; and colocalized binding of clock and tissue-specific transcription factors at distal enhancers. Chromosome conformation capture (4C-seq) in liver and kidney identified liver-specific chromatin loops that recruited clock-bound enhancers to promoters to regulate liver-specific transcriptional rhythms. Furthermore, this looping was remarkably promoter-specific on the scale of less than 10 kilobases (kb). Enhancers can contact a rhythmic promoter while looping out nearby nonrhythmic alternative promoters, confining rhythmic enhancer activity to specific promoters. These findings suggest that chromatin folding enables the clock to regulate rhythmic transcription of specific promoters to output temporal transcriptional programs tailored to different tissues.

Feminization of Male Mouse Liver by Persistent Growth Hormone Stimulation: Activation of Sex-Biased Transcriptional Networks and Dynamic Changes in Chromatin States

Sex-dependent pituitary growth hormone (GH) secretory profiles-pulsatile in males and persistent in females-regulate the sex-biased, STAT5-dependent expression of hundreds of genes in mouse liver, imparting sex differences in hepatic drug/lipid metabolism and disease risk. Here, we examine transcriptional and epigenetic changes induced by continuous GH infusion (cGH) in male mice, which rapidly feminizes the temporal profile of liver STAT5 activity. cGH repressed 86% of male-biased genes and induced 68% of female-biased genes within 4 days; however, several highly female-specific genes showed weak or no feminization, even after 14 days of cGH treatment. Female-biased genes already in an active chromatin state in male liver generally showed early cGH responses; genes in an inactive chromatin state often responded late. Early cGH-responsive genes included those encoding two GH/STAT5-regulated transcriptional repressors: male-biased BCL6, which was repressed, and female-specific CUX2, which was induced. Male-biased genes activated by STAT5 and/or repressed by CUX2 were enriched for early cGH repression. Female-biased BCL6 targets were enriched for early cGH derepression. Changes in sex-specific chromatin accessibility and histone modifications accompanied these cGH-induced sex-biased gene expression changes. Thus, the temporal, sex-biased gene responses to persistent GH stimulation are dictated by GH/STAT5-regulated transcription factors arranged in a hierarchical network and by the dynamics of changes in sex-biased epigenetic states.

Activation of Male Liver Chromatin Accessibility and STAT5-Dependent Gene Transcription by Plasma Growth Hormone Pulses

Sex differences in pituitary growth hormone (GH) secretion (pulsatile in males vs near continuous/persistent in females) impart sex-dependent expression to hundreds of genes in adult mouse liver. Signal transducer and activator of transcription (STAT) 5, a GH-activated transcription factor that is essential for liver sexual dimorphism, is dynamically activated in direct response to each male plasma GH pulse. However, the impact of GH-induced STAT5 pulses on liver chromatin accessibility and downstream transcriptional events is unknown. In this study, we investigated the impact of a single pulse of GH given to hypophysectomized mice on local liver chromatin accessibility (DNase hypersensitive site analysis), transcription rates (heterogeneous nuclear RNA analysis), and gene expression (quantitative polymerase chain reaction and RNA sequencing) determined 30, 90, or 240 minutes later. The STAT5-dependent but sex-independent early GH response genes Igf1 and Cish showed rapid, GH pulse-induced increases in chromatin accessibility and gene transcription, reversing the effects of hypophysectomy. Rapid increases in liver chromatin accessibility and transcriptional activity were also induced in hypophysectomized male mice for some (Ces2b, Ugt2b38) but not for other liver STAT5-dependent male-biased genes (Cyp7b1). Moreover, in pituitary-intact male mice, Igf1, Cish, Ces2b, and Ugt2b38 all showed remarkable cycles of chromatin opening and closing, as well as associated cycles of induced gene transcription, which closely followed each endogenous pulse of liver STAT5 activity. Thus, the endogenous rhythms of male plasma GH pulsation dynamically open and then close liver chromatin at discrete, localized regulatory sites in temporal association with transcriptional activation of Igf1, Cish, and a subset of STAT5-dependent male-biased genes.

Regulation of drug metabolism and toxicity by multiple factors of genetics, epigenetics, lncRNAs, gut microbiota, and diseases: a meeting report of the 21st International Symposium on Microsomes and Drug Oxidations (MDO)

Variations in drug metabolism may alter drug efficacy and cause toxicity; better understanding of the mechanisms and risks shall help to practice precision medicine. At the 21^st International Symposium on Microsomes and Drug Oxidations held in Davis, California, USA, in October 2-6, 2016, a number of speakers reported some new findings and ongoing studies on the regulation mechanisms behind variable drug metabolism and toxicity, and discussed potential implications to personalized medications. A considerably insightful overview was provided on genetic and epigenetic regulation of gene expression involved in drug absorption, distribution, metabolism, and excretion (ADME) and drug response. Altered drug metabolism and disposition as well as molecular mechanisms among diseased and special populations were presented. In addition, the roles of gut microbiota in drug metabolism and toxicology as well as long non-coding RNAs in liver functions and diseases were discussed. These findings may offer new insights into improved understanding of ADME regulatory mechanisms and advance drug metabolism research.

Hepatic Long Intergenic Noncoding RNAs: High Promoter Conservation and Dynamic, Sex-Dependent Transcriptional Regulation by Growth Hormone

Long intergenic noncoding RNAs (lincRNAs) are increasingly recognized as key chromatin regulators, yet few studies have characterized lincRNAs in a single tissue under diverse conditions. Here, we analyzed 45 mouse liver RNA sequencing (RNA-Seq) data sets collected under diverse conditions to systematically characterize 4,961 liver lincRNAs, 59% of them novel, with regard to gene structures, species conservation, chromatin accessibility, transcription factor binding, and epigenetic states. To investigate the potential for functionality, we focused on the responses of the liver lincRNAs to growth hormone stimulation, which imparts clinically relevant sex differences to hepatic metabolism and liver disease susceptibility. Sex-biased expression characterized 247 liver lincRNAs, with many being nuclear RNA enriched and regulated by growth hormone. The sex-biased lincRNA genes are enriched for nearby and correspondingly sex-biased accessible chromatin regions, as well as sex-biased binding sites for growth hormone-regulated transcriptional activators (STAT5, hepatocyte nuclear factor 6 [HNF6], FOXA1, and FOXA2) and transcriptional repressors (CUX2 and BCL6). Repression of female-specific lincRNAs in male liver, but not that of male-specific lincRNAs in female liver, was associated with enrichment of H3K27me3-associated inactive states and poised (bivalent) enhancer states. Strikingly, we found that liver-specific lincRNA gene promoters are more highly species conserved and have a significantly higher frequency of proximal binding by liver transcription factors than liver-specific protein-coding gene promoters. Orthologs for many liver lincRNAs were identified in one or more supraprimates, including two rat lincRNAs showing the same growth hormone-regulated, sex-biased expression as their mouse counterparts. This integrative analysis of liver lincRNA chromatin states, transcription factor occupancy, and growth hormone regulation provides novel insights into the expression of sex-specific lincRNAs and their potential for regulation of sex differences in liver physiology and disease.