Unbiased, genome-wide in vivo mapping of transcriptional regulatory elements reveals sex differences in chromatin structure associated with sex-specific liver gene expression

HDI-STARR-seq: Condition-specific enhancer discovery in mouse liver in vivo

Background STARR-seq and other massively-parallel reporter assays are widely used to discover functional enhancers in transfected cell models, which can be confounded by plasmid vector-induced type-I interferon immune responses and lack the multicellular environment and endogenous chromatin state of complex mammalian tissues. Results Here, we describe HDI-STARR-seq, which combines STARR-seq plasmid library delivery to the liver, by hydrodynamic tail vein injection (HDI), with reporter RNA transcriptional initiation driven by a minimal Albumin promoter, which we show is essential for mouse liver STARR-seq enhancer activity assayed 7 days after HDI. Importantly, little or no vector-induced innate type-I interferon responses were observed. Comparisons of HDI-STARR-seq activity between male and female mouse livers and in livers from males treated with an activating ligand of the transcription factor CAR ( Nr1i3 ) identified many condition-dependent enhancers linked to condition-specific gene expression. Further, thousands of active liver enhancers were identified using a high complexity STARR-seq library comprised of ~ 50,000 genomic regions released by DNase-I digestion of mouse liver nuclei. When compared to stringently inactive library sequences, the active enhancer sequences identified were highly enriched for liver open chromatin regions with activating histone marks (H3K27ac, H3K4me1, H3K4me3), were significantly closer to gene transcriptional start sites, and were significantly depleted of repressive (H3K27me3, H3K9me3) and transcribed region histone marks (H3K36me3). Conclusions HDI-STARR-seq offers substantial improvements over current methodologies for large scale, functional profiling of enhancers, including condition-dependent enhancers, in liver tissue in vivo, and can be adapted to characterize enhancer activities in a variety of species and tissues by selecting suitable tissue- and species-specific promoter sequences.

HDI-STARR-seq: Condition-specific enhancer discovery in mouse liver in vivo

STARR-seq and other massively-parallel reporter assays are widely used to discover functional enhancers in transfected cell models, which can be confounded by plasmid vector-induced type-I interferon immune responses and lack the multicellular environment and endogenous chromatin state of complex mammalian tissues. Here, we describe HDI-STARR-seq, which combines STARR-seq plasmid library delivery to the liver, by hydrodynamic tail vein injection (HDI), with reporter RNA transcriptional initiation driven by a minimal Albumin promoter, which we show is essential for mouse liver STARR-seq enhancer activity assayed 7 days after HDI. Importantly, little or no vector-induced innate type-I interferon responses were observed. Comparisons of HDI-STARR-seq activity between male and female mouse livers and in livers from males treated with an activating ligand of the transcription factor CAR (Nr1i3) identified many condition-dependent enhancers linked to condition-specific gene expression. Further, thousands of active liver enhancers were identified using a high complexity STARR-seq library comprised of ~50,000 genomic regions released by DNase-I digestion of mouse liver nuclei. When compared to stringently inactive library sequences, the active enhancer sequences identified were highly enriched for liver open chromatin regions with activating histone marks (H3K27ac, H3K4me1, H3K4me3), were significantly closer to gene transcriptional start sites, and were significantly depleted of repressive (H3K27me3, H3K9me3) and transcribed region histone marks (H3K36me3). HDI-STARR-seq offers substantial improvements over current methodologies for large scale, functional profiling of enhancers, including condition-dependent enhancers, in liver tissue in vivo, and can be adapted to characterize enhancer activities in a variety of species and tissues by selecting suitable tissue- and species-specific promoter sequences.

Sex differences in brain cell-type specific chromatin accessibility in schizophrenia

Our understanding of the sex-specific role of the non-coding genome in serious mental illness remains largely incomplete. To address this gap, we explored sex differences in 1,393 chromatin accessibility profiles, derived from neuronal and non-neuronal nuclei of two distinct cortical regions from 234 cases with serious mental illness and 235 controls. We identified sex-specific enhancer-promoter interactions and showed that they regulate genes involved in X-chromosome inactivation (XCI). Examining chromosomal conformation allowed us to identify sex-specific cis- and trans-regulatory domains (CRDs and TRDs). Co-localization of sex-specific TRDs with schizophrenia common risk variants pinpointed male-specific regulatory regions controlling a number of metabolic pathways. Additionally, enhancers from female-specific TRDs were found to regulate two genes known to escape XCI, (XIST and JPX), underlying the importance of TRDs in deciphering sex differences in schizophrenia. Overall, these findings provide extensive characterization of sex differences in the brain epigenome and disease-associated regulomes.

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.

Multiomics Analyses Reveal Sex Differences in Mouse Renal Proximal Subsegments

SIGNIFICANCE STATEMENT

Sex-dependent differences in kidney function are recognized but the underlying molecular mechanisms are largely unexplored. Advances in genomics and proteomic technologies now allow extensive characterization of differences between the same cell types of males and females. Multiomics integrating RNA-seq, ATAC-seq, and proteomics data to investigate differences in gene expression, chromatin accessibility, and protein expression in proximal tubules of male and female mice identified many sex-biased genes and proteins associated with kidney functions, including metabolic and transport processes. Sex differences may also arise from variations of the interaction between transcription factors and accessible chromatin regions. A comprehensive web resource is provided to advance understanding of sex differences in cells of the proximal tubule.

BACKGROUND

Sex differences have been increasingly recognized as important in kidney physiology and pathophysiology, but limited resources are available for comprehensive interrogation of sex differences.

METHODS

RNA-seq and ATAC-seq of microdissected mouse proximal tubules and protein mass spectrometry of homogenized perfused mouse kidneys reveal differences in proximal tubule cells of males and females.

RESULTS

The transcriptomic data indicated that the major differences in the proximal tubules between the sexes were in the S2/S3 segments, and most of the sex-biased transcripts mapped to autosomes rather than to the sex chromosomes. Many of the transcripts exhibiting sex-biased expression are involved in monocarboxylic acid metabolic processes, organic anion transport, and organic acid transport. The ATAC-seq method on microdissected tubules captured chromatin accessibility. Many of the more than 7000 differentially accessible DNA regions identified were in distal regions. Motif analyses revealed a lack of direct involvement of estrogen receptors or the androgen receptor (absence of canonical hormone response elements), suggesting an indirect regulatory role of sex hormones. Instead, analyses identified several transcription factors (TFs) ( Tead1 , Nfia/b , and Pou3f3 ) whose interplay with proximal tubule-specific TFs ( e.g. , Hnf1b , Hnf4a ) may contribute to sex differences. Finally, the whole-kidney proteome was correlated with the transcriptome, and many sex-biased proteins ( e.g. , Cyp2e1, Acsm2/3) were identified.

CONCLUSIONS

Sex-dependent cis-regulatory elements interact with TFs in ways that lead to sex-biased gene expression in proximal tubule cells. These data are provided as a user-friendly web page at https://esbl.nhlbi.nih.gov/MRECA/PT/ .

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.

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.

High Risk-Human Papillomavirus in HNSCC: Present and Future Challenges for Epigenetic Therapies

Head and Neck Squamous Cell Carcinoma (HNSCC) is a highly heterogeneous group of tumors characterized by an incidence of 650,000 new cases and 350,000 deaths per year worldwide and a male to female ratio of 3:1. The main risk factors are alcohol and tobacco consumption and Human Papillomavirus (HPV) infections. HNSCC cases are divided into two subgroups, the HPV-negative (HPV−) and the HPV-positive (HPV+) which have different clinicopathological and molecular profiles. However, patients are still treated with the same therapeutic regimens. It is thus of utmost importance to characterize the molecular mechanisms underlying these differences to find new biomarkers and novel therapeutic targets towards personalized therapies. Epigenetic alterations are a hallmark of cancer and can be exploited as both promising biomarkers and potential new targets. E6 and E7 HPV oncoviral proteins besides targeting p53 and pRb, impair the expression and the activity of several epigenetic regulators. While alterations in DNA methylation patterns have been well described in HPV+ and HPV− HNSCC, accurate histone post-translational modifications (hPTMs) characterization is still missing. Herein, we aim to provide an updated overview on the impact of HPV on the hPTMs landscape in HNSCC. Moreover, we will also discuss the sex and gender bias in HNSCC and how the epigenetic machinery could be involved in this process, and the importance of taking into account sex and/or gender also in this field.

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.

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.

Hepatic steatosis and steatohepatitis: a functional meta-analysis of sex-based differences in transcriptomic studies

Background

Previous studies have described sex-based differences in the epidemiological and clinical patterns of non-alcoholic fatty liver disease (NAFLD); however, we understand relatively little regarding the underlying molecular mechanisms. Herein, we present the first systematic review and meta-analysis of NAFLD transcriptomic studies to identify sex-based differences in the molecular mechanisms involved during the steatosis (NAFL) and steatohepatitis (NASH) stages of the disease.

Methods

Transcriptomic studies in the Gene Expression Omnibus database were systematically reviewed following the PRISMA statement guidelines. For each study, NAFL and NASH in premenopausal women and men were compared using a dual strategy: gene-set analysis and pathway activity analysis. Finally, the functional results of all studies were integrated into a meta-analysis.

Results

We reviewed a total of 114 abstracts and analyzed seven studies that included 323 eligible patients. The meta-analyses identified significantly altered molecular mechanisms between premenopausal women and men, including the overrepresentation of genes associated with DNA regulation, vinculin binding, interleukin-2 responses, negative regulation of neuronal death, and the transport of ions and cations in premenopausal women. In men, we discovered the overrepresentation of genes associated with the negative regulation of interleukin-6 and the establishment of planar polarity involved in neural tube closure.

Conclusions

Our meta-analysis of transcriptomic data provides a powerful approach to identify sex-based differences in NAFLD. We detected differences in relevant biological functions and molecular terms between premenopausal women and men. Differences in immune responsiveness between men and premenopausal women with NAFLD suggest that women possess a more immune tolerant milieu, while men display an impaired liver regenerative response.

Graphical abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13293-021-00368-1.

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.

Genome-Wide Sex and Gender Differences in Cancer

Despite their known importance in clinical medicine, differences based on sex and gender are among the least studied factors affecting cancer susceptibility, progression, survival, and therapeutic response. In particular, the molecular mechanisms driving sex differences are poorly understood and so most approaches to precision medicine use mutational or other genomic data to assign therapy without considering how the sex of the individual might influence therapeutic efficacy. The mandate by the National Institutes of Health that research studies include sex as a biological variable has begun to expand our understanding on its importance. Sex differences in cancer may arise due to a combination of environmental, genetic, and epigenetic factors, as well as differences in gene regulation, and expression. Extensive sex differences occur genome-wide, and ultimately influence cancer biology and outcomes. In this review, we summarize the current state of knowledge about sex-specific genetic and genome-wide influences in cancer, describe how differences in response to environmental exposures and genetic and epigenetic alterations alter the trajectory of the disease, and provide insights into the importance of integrative analyses in understanding the interplay of sex and genomics in cancer. In particular, we will explore some of the emerging analytical approaches, such as the use of network methods, that are providing a deeper understanding of the drivers of differences based on sex and gender. Better understanding these complex factors and their interactions will improve cancer prevention, treatment, and outcomes for all individuals.

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.

Age-related hearing loss: Why we need to think about sex as a biological variable

It has long been known that age-related hearing loss (ARHL) is more common, more severe, and with an earlier onset in men compared to women. Even in the absence of confounding factors such as noise exposure, these sexdifferences in susceptibility to ARHL remain. In the last decade, insight into the pleiotrophic nature by which estrogen signaling can impact multiple signaling mechanisms to mediate downstream changes in gene expression and/or elicit rapid changes in cellular function has rapidly gathered pace, and a role for estrogen signaling in the biological pathways that confer neuroprotection is becoming undeniable. Here I review the evidence why we need to consider sex as a biological variable (SABV) when investigating the etiology of ARHL. Loss of auditory function with aging is frequency-specific and modulated by SABV. Evidence also suggests that differences in cochlear physiology between women and men are already present from birth. Understanding the molecular basis of these sex differences in ARHL will accelerate the development of precision medicine therapies for ARHL.

Unravelling the complex genetics of common kidney diseases: from variants to mechanisms

Genome-wide association studies (GWAS) have identified hundreds of loci associated with kidney-related traits such as glomerular filtration rate, albuminuria, hypertension, electrolyte and metabolite levels. However, these impressive, large-scale mapping approaches have not always translated into an improved understanding of disease or development of novel therapeutics. GWAS have several important limitations. Nearly all disease-associated risk loci are located in the non-coding region of the genome and therefore, their target genes, affected cell types and regulatory mechanisms remain unknown. Genome-scale approaches can be used to identify associations between DNA sequence variants and changes in gene expression (quantified through bulk and single-cell methods), gene regulation and other molecular quantitative trait studies, such as chromatin accessibility, DNA methylation, protein expression and metabolite levels. Data obtained through these approaches, used in combination with robust computational methods, can deliver robust mechanistic inferences for translational exploitation. Understanding the genetic basis of common kidney diseases means having a comprehensive picture of the genes that have a causal role in disease development and progression, of the cells, tissues and organs in which these genes act to affect the disease, of the cellular pathways and mechanisms that drive disease, and of potential targets for disease prevention, detection and therapy.

The impact of sex differences on genomic research

Sex and gender differences affect all dimensions of human health ranging from the biological basis of disease to therapeutic access, choice and response. Genomics research has long ignored the role of sex differences as potential modulators and the concept is gaining more attention only recently. In the present review we summarize the current knowledge of the impact of sex differences on genomic and epigenomic research, the potential interaction of genomics and gender and the role of these differences in disease etiopathogenesis. Sex differences can emerge from differences in the sex chromosomes themselves, from their interaction with the genome and from the influence of hormones on genomic processes. The impact of these processes on the incidence of autoimmune and oncologic disease is well documented. The growing field of systems biology, which aims at integrating information from different networks of the human body, could also greatly benefit from this approach. In the present review we summarize the current knowledge and provide recommendations for the future performance of sex-sensitive genomics research.

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.

The role of sex in the genomics of human complex traits

Nearly all human complex traits and disease phenotypes exhibit some degree of sex differences, including differences in prevalence, age of onset, severity or disease progression. Until recently, the underlying genetic mechanisms of such sex differences have been largely unexplored. Advances in genomic technologies and analytical approaches are now enabling a deeper investigation into the effect of sex on human health traits. In this Review, we discuss recent insights into the genetic models and mechanisms that lead to sex differences in complex traits. This knowledge is critical for developing deeper insight into the fundamental biology of sex differences and disease processes, thus facilitating precision medicine.

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

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

Impact of CAR Agonist Ligand TCPOBOP on Mouse Liver Chromatin Accessibility

Activation of the nuclear receptor and transcription factor CAR (Nr1i3) by its specific agonist ligand TCPOBOP (1, 4-bis[2-(3, 5-dichloropyridyloxy)]benzene) dysregulates hundreds of genes in mouse liver and is linked to male-biased hepatocarcinogenesis. To elucidate the genomic organization of CAR-induced gene responses, we investigated the distribution of TCPOBOP-responsive RefSeq coding and long noncoding RNA (lncRNA) genes across the megabase-scale topologically associating domains (TADs) that segment the genome, and which provide a structural framework that functionally constrains enhancer-promoter interactions. We show that a subset of TCPOBOP-responsive genes cluster within TADs, and that TCPOBOP-induced genes and TCPOBOP-repressed genes are often found in different TADs. Further, using DNase-seq and DNase hypersensitivity site (DHS) analysis, we identified several thousand genomic regions (ΔDHS) where short-term exposure to TCPOBOP induces localized changes (increases or decreases) in mouse liver chromatin accessibility, many of which cluster in TADs together with TCPOBOP-responsive genes. Sites of chromatin opening were highly enriched nearby genes induced by TCPOBOP and chromatin closing was highly enriched nearby genes repressed by TCPOBOP, consistent with TCPOBOP-responsive ΔDHS serving as enhancers and promoters that positively regulate CAR-responsive genes. Gene expression changes lagged behind chromatin opening or closing for a subset of TCPOBOP-responsive ΔDHS. ΔDHS that were specifically responsive to TCPOBOP in male liver were significantly enriched for genomic regions with a basal male bias in chromatin accessibility; however, the male-biased response of hepatocellular carcinoma-related genes to TCPOBOP was not associated with a correspondingly male-biased ΔDHS response. These studies elucidate the genome-wide organization of CAR-responsive genes and of the thousands of associated genomic sites where TCPOBOP exposure induces both rapid and persistent changes in chromatin accessibility.

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

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

Loss of growth hormone-mediated signal transducer and activator of transcription 5 (STAT5) signaling in mice results in insulin sensitivity with obesity

Growth hormone (GH) has an important function as an insulin antagonist with elevated insulin sensitivity evident in humans and mice lacking a functional GH receptor (GHR). We sought the molecular basis for this sensitivity by utilizing a panel of mice possessing specific deletions of GHR signaling pathways. Metabolic clamps and glucose homeostasis tests were undertaken in these obese adult C57BL/6 male mice, which indicated impaired hepatic gluconeogenesis. Insulin sensitivity and glucose disappearance rate were enhanced in muscle and adipose of mice lacking the ability to activate the signal transducer and activator of transcription (STAT)5 via the GHR (Ghr-391^-/-) as for GHR-null (GHR^-/-) mice. These changes were associated with a striking inhibition of hepatic glucose output associated with altered glycogen metabolism and elevated hepatic glycogen content during unfed state. The enhanced hepatic insulin sensitivity was associated with increased insulin receptor β and insulin receptor substrate 1 activation along with activated downstream protein kinase B signaling cascades. Although phosphoenolpyruvate carboxykinase (Pck)-1 expression was unchanged, its inhibitory acetylation was elevated because of decreased sirtuin-2 expression, thereby promoting loss of PCK1. Loss of STAT5 signaling to defined chromatin immunoprecipitation targets would further increase lipogenesis, supporting hepatosteatosis while lowering glucose output. Finally, up-regulation of IL-15 expression in muscle, with increased secretion of adiponectin and fibroblast growth factor 1 from adipose tissue, is expected to promote insulin sensitivity.-Chhabra, Y., Nelson, C. N., Plescher, M., Barclay, J. L., Smith, A. G., Andrikopoulos, S., Mangiafico, S., Waxman, D. J., Brooks, A. J., Waters, M. J. Loss of growth hormone-mediated signal transducer and activator of transcription 5 (STAT5) signaling in mice results in insulin sensitivity with obesity.

Insulin signaling and reduced glucocorticoid receptor activity attenuate postprandial gene expression in liver

Hepatic circadian gene transcription is tightly coupled to feeding behavior, which has a profound impact on metabolic disorders associated with diet-induced obesity. Here, we describe a genomics approach to uncover mechanisms controlling hepatic postprandial gene expression. Combined transcriptomic and cistromic analysis identified hundreds of circadian-regulated genes and enhancers controlled by feeding. Postprandial suppression of enhancer activity was associated with reduced glucocorticoid receptor (GR) and Forkhead box O1 (FOXO1) occupancy of chromatin correlating with reduced serum corticosterone levels and increased serum insulin levels. Despite substantial co-occupancy of feeding-regulated enhancers by GR and FOXO1, selective disruption of corticosteroid and/or insulin signaling resulted in dysregulation of specific postprandial regulated gene programs. In combination, these signaling pathways operate a major part of the genes suppressed by feeding. Importantly, the feeding response was disrupted in diet-induced obese animals, which was associated with dysregulation of several corticosteroid- and insulin-regulated genes, providing mechanistic insights to dysregulated circadian gene transcription associated with obesity.

The liver is an essential organ regulating metabolic homeostasis in response to fluctuations of metabolites induced by daily rhythms of food intake. Homeostasis is maintained by precise dynamic regulation of signaling pathways controlling a wealth of enzymatic reactions involving lipid, bile acid, amino acid and glucose synthesis, storage, and redistribution in hepatocytes. Precise temporal expression of hepatic enzymes is crucial for metabolic homeostasis; a major part of circadian hepatic protein expression is regulated by precisely timed gene transcription. Here, we use a genomics approach to identify genes and regulatory regions of the genome involved in feeding-regulated gene expression. We find that transcription factors acting downstream of glucocorticoid and insulin signaling are enriched at regulatory regions repressed by feeding. Importantly, insulin and glucocorticoid signaling operate to cooperatively control the majority of feeding-mediated gene repression, and these signaling pathways are dysregulated in diet-induced obesity impacting dynamic hepatic gene expression.

Sex differences in lipid metabolism are affected by presence of the gut microbiota

Physiological processes are differentially regulated between men and women. Sex and gut microbiota have each been demonstrated to regulate host metabolism, but it is unclear whether both factors are interdependent. Here, we determined to what extent sex-specific differences in lipid metabolism are modulated via the gut microbiota. While male and female Conv mice showed predominantly differential expression in gene sets related to lipid metabolism, GF mice showed differences in gene sets linked to gut health and inflammatory responses. This suggests that presence of the gut microbiota is important in sex-specific regulation of lipid metabolism. Further, we explored the role of bile acids as mediators in the cross-talk between the microbiome and host lipid metabolism. Females showed higher total and primary serum bile acids levels, independent of presence of microbiota. However, in presence of microbiota we observed higher secondary serum bile acid levels in females compared to males. Analysis of microbiota composition displayed sex-specific differences in Conv mice. Therefore, our data suggests that bile acids possibly play a role in the crosstalk between the microbiome and sex-specific regulation of lipid metabolism. In conclusion, our data shows that presence of the gut microbiota contributes to sex differences in lipid metabolism.

Personalized therapy when tackling nonalcoholic fatty liver disease: a focus on sex, genes, and drugs

INTRODUCTION

Nonalcoholic fatty liver disease (NAFLD) is the most frequent liver disease in the world. It describes a term for a group of hepatic diseases including steatosis, fibrosis, and cirrhosis that can finally lead to hepatocellular carcinoma. There are many factors influencing NAFLD initiation and progression, such as obesity, dyslipidemia, insulin resistance, genetic factors, and hormonal changes. However, there is also lean-NAFLD which is not associated with obesity. NAFLD is considered to be a sexually dimorphic disease. In most cases, men have a higher prevalence for the disease compared to premenopausal women. Areas covered: In this review, we first summarize the NAFLD disease epidemiology, pathology, and diagnosis. We describe NAFLD progression with the focus on sexual and genetic differences for disease development and pharmacological treatment. Personalized treatment for multifactorial NAFLD is discussed in consideration of different factors, including genetics, gender and sex. Expert opinion: The livers of female and male NAFLD patients have different metabolic capacities which influence the metabolism of all drugs applied to such patients. This aspect is not yet sufficiently taken into account. The liver computational models might quicken the pace toward assessing personalized disease progression and treatment options.

Computational prediction of CTCF/cohesin-based intra-TAD loops that insulate chromatin contacts and gene expression in mouse liver

CTCF and cohesin are key drivers of 3D-nuclear organization, anchoring the megabase-scale Topologically Associating Domains (TADs) that segment the genome. Here, we present and validate a computational method to predict cohesin-and-CTCF binding sites that form intra-TAD DNA loops. The intra-TAD loop anchors identified are structurally indistinguishable from TAD anchors regarding binding partners, sequence conservation, and resistance to cohesin knockdown; further, the intra-TAD loops retain key functional features of TADs, including chromatin contact insulation, blockage of repressive histone mark spread, and ubiquity across tissues. We propose that intra-TAD loops form by the same loop extrusion mechanism as the larger TAD loops, and that their shorter length enables finer regulatory control in restricting enhancer-promoter interactions, which enables selective, high-level expression of gene targets of super-enhancers and genes located within repressive nuclear compartments. These findings elucidate the role of intra-TAD cohesin-and-CTCF binding in nuclear organization associated with widespread insulation of distal enhancer activity.

Implications of STAT3 and STAT5 signaling on gene regulation and chromatin remodeling in hematopoietic cancer

STAT3 and STAT5 proteins are oncogenic downstream mediators of the JAK–STAT pathway. Deregulated STAT3 and STAT5 signaling promotes cancer cell proliferation and survival in conjunction with other core cancer pathways. Nuclear phosphorylated STAT3 and STAT5 regulate cell-type-specific transcription profiles via binding to promoter elements and exert more complex functions involving interaction with various transcriptional coactivators or corepressors and chromatin remodeling proteins. The JAK–STAT pathway can rapidly reshape the chromatin landscape upon cytokine, hormone, or growth factor stimulation and unphosphorylated STAT proteins also appear to be functional with respect to regulating chromatin accessibility. Notably, cancer genome landscape studies have implicated mutations in various epigenetic modifiers as well as the JAK–STAT pathway as underlying causes of many cancers, particularly acute leukemia and lymphomas. However, it is incompletely understood how mutations within these pathways can interact and synergize to promote cancer. We summarize the current knowledge of oncogenic STAT3 and STAT5 functions downstream of cytokine signaling and provide details on prerequisites for DNA binding and gene transcription. We also discuss key interactions of STAT3 and STAT5 with chromatin remodeling factors such as DNA methyltransferases, histone modifiers, cofactors, corepressors, and other transcription factors.

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.

Functional Roles of Sex-Biased, Growth Hormone-Regulated MicroRNAs miR-1948 and miR-802 in Young Adult Mouse Liver

Sex-specific temporal patterns of pituitary growth hormone (GH) secretion determine the sex-biased transcription of hundreds of genes in the liver and impart important sex differences in liver physiology, metabolism, and disease. Sex differences in hepatic gene expression vary widely, ranging from less than twofold to >1000-fold in the mouse. Here, we use small RNA sequencing to discover 24 sex-biased mouse liver microRNAs (miRNAs), and then investigate the roles of two of these miRNAs in GH-regulated liver sex differences. Studies in prepubertal and young adult mice, and in mice in which pituitary hormones are ablated or where sex-specific hepatic GH signaling is dysregulated, demonstrated that the male-biased miR-1948 and the female-biased miR-802 are both regulated by sex-specific pituitary GH secretory patterns, acquire sex specificity at puberty, and are dependent on the GH-activated transcription factor STAT5 for their sex-specific expression. Both miRNAs are within genomic regions characterized by sex-biased chromatin accessibility. miR-1948, an uncharacterized miRNA, has essential features for correct Drosha/Dicer processing, generates a bona fide mature miRNA with strong strand bias for the 5p arm, and is bound by Argonaute in liver tissue, as is miR-802. In vivo studies using inhibitory locked nucleic acid sequences revealed that miR-1948-5p preferentially represses female-biased messenger RNAs (mRNAs) and induces male-biased mRNAs in male liver; conversely, miR-802-5p preferentially represses male-biased mRNAs and increases levels of female-biased mRNAs in female liver. Cytochrome P450 mRNAs were strongly enriched as targets of both miRNAs. Thus, miR-1948-5p and miR-802-5p are functional components of the GH regulatory network that shapes sex-differential gene expression in mouse liver.

Regulation of circadian clock transcriptional output by CLOCK:BMAL1

The mammalian circadian clock relies on the transcription factor CLOCK:BMAL1 to coordinate the rhythmic expression of 15% of the transcriptome and control the daily regulation of biological functions. The recent characterization of CLOCK:BMAL1 cistrome revealed that although CLOCK:BMAL1 binds synchronously to all of its target genes, its transcriptional output is highly heterogeneous. By performing a meta-analysis of several independent genome-wide datasets, we found that the binding of other transcription factors at CLOCK:BMAL1 enhancers likely contribute to the heterogeneity of CLOCK:BMAL1 transcriptional output. While CLOCK:BMAL1 rhythmic DNA binding promotes rhythmic nucleosome removal, it is not sufficient to generate transcriptionally active enhancers as assessed by H3K27ac signal, RNA Polymerase II recruitment, and eRNA expression. Instead, the transcriptional activity of CLOCK:BMAL1 enhancers appears to rely on the activity of ubiquitously expressed transcription factors, and not tissue-specific transcription factors, recruited at nearby binding sites. The contribution of other transcription factors is exemplified by how fasting, which effects several transcription factors but not CLOCK:BMAL1, either decreases or increases the amplitude of many rhythmically expressed CLOCK:BMAL1 target genes. Together, our analysis suggests that CLOCK:BMAL1 promotes a transcriptionally permissive chromatin landscape that primes its target genes for transcription activation rather than directly activating transcription, and provides a new framework to explain how environmental or pathological conditions can reprogram the rhythmic expression of clock-controlled genes.

Sex-biased microRNA expression in mammals and birds reveals underlying regulatory mechanisms and a role in dosage compensation

Sexual dimorphism depends on sex-biased gene expression, but the contributions of microRNAs (miRNAs) have not been globally assessed. We therefore produced an extensive small RNA sequencing data set to analyze male and female miRNA expression profiles in mouse, opossum, and chicken. Our analyses uncovered numerous cases of somatic sex-biased miRNA expression, with the largest proportion found in the mouse heart and liver. Sex-biased expression is explained by miRNA-specific regulation, including sex-biased chromatin accessibility at promoters, rather than piggybacking of intronic miRNAs on sex-biased protein-coding genes. In mouse, but not opossum and chicken, sex bias is coordinated across tissues such that autosomal testis-biased miRNAs tend to be somatically male-biased, whereas autosomal ovary-biased miRNAs are female-biased, possibly due to broad hormonal control. In chicken, which has a Z/W sex chromosome system, expression output of genes on the Z Chromosome is expected to be male-biased, since there is no global dosage compensation mechanism that restores expression in ZW females after almost all genes on the W Chromosome decayed. Nevertheless, we found that the dominant liver miRNA, miR-122-5p, is Z-linked but expressed in an unbiased manner, due to the unusual retention of a W-linked copy. Another Z-linked miRNA, the male-biased miR-2954-3p, shows conserved preference for dosage-sensitive genes on the Z Chromosome, based on computational and experimental data from chicken and zebra finch, and acts to equalize male-to-female expression ratios of its targets. Unexpectedly, our findings thus establish miRNA regulation as a novel gene-specific dosage compensation mechanism.

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.

TFOS DEWS II Sex, Gender, and Hormones Report

One of the most compelling features of dry eye disease (DED) is that it occurs more frequently in women than men. In fact, the female sex is a significant risk factor for the development of DED. This sex-related difference in DED prevalence is attributed in large part to the effects of sex steroids (e.g. androgens, estrogens), hypothalamic-pituitary hormones, glucocorticoids, insulin, insulin-like growth factor 1 and thyroid hormones, as well as to the sex chromosome complement, sex-specific autosomal factors and epigenetics (e.g. microRNAs). In addition to sex, gender also appears to be a risk factor for DED. "Gender" and "sex" are words that are often used interchangeably, but they have distinct meanings. "Gender" refers to a person's self-representation as a man or woman, whereas "sex" distinguishes males and females based on their biological characteristics. Both gender and sex affect DED risk, presentation of the disease, immune responses, pain, care-seeking behaviors, service utilization, and myriad other facets of eye health. Overall, sex, gender and hormones play a major role in the regulation of ocular surface and adnexal tissues, and in the difference in DED prevalence between women and men. The purpose of this Subcommittee report is to review and critique the nature of this role, as well as to recommend areas for future research to advance our understanding of the interrelationships between sex, gender, hormones and DED.

Transcriptional role of androgen receptor in the expression of long non-coding RNA Sox2OT in neurogenesis

The complex architecture of adult brain derives from tightly regulated migration and differentiation of precursor cells generated during embryonic neurogenesis. Changes at transcriptional level of genes that regulate migration and differentiation may lead to neurodevelopmental disorders. Androgen receptor (AR) is a transcription factor that is already expressed during early embryonic days. However, AR role in the regulation of gene expression at early embryonic stage is yet to be determinate. Long non-coding RNA (lncRNA) Sox2 overlapping transcript (Sox2OT) plays a crucial role in gene expression control during development but its transcriptional regulation is still to be clearly defined. Here, using Bicalutamide in order to pharmacologically inactivated AR, we investigated whether AR participates in the regulation of the transcription of the lncRNASox2OTat early embryonic stage. We identified a new DNA binding region upstream of Sox2 locus containing three androgen response elements (ARE), and found that AR binds such a sequence in embryonic neural stem cells and in mouse embryonic brain. Our data suggest that through this binding, AR can promote the RNA polymerase II dependent transcription of Sox2OT. Our findings also suggest that AR participates in embryonic neurogenesis through transcriptional control of the long non-coding RNA Sox2OT.

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.

Hyper- and hypo- nutrition studies of the hepatic transcriptome and epigenome suggest that PPARα regulates anaerobic glycolysis

Diet plays a crucial role in shaping human health and disease. Diets promoting obesity and insulin resistance can lead to severe metabolic diseases, while calorie-restricted (CR) diets can improve health and extend lifespan. In this work, we fed mice either a chow diet (CD), a 16 week high-fat diet (HFD), or a CR diet to compare and contrast the effects of these diets on mouse liver biology. We collected transcriptomic and epigenomic datasets from these mice using RNA-Seq and DNase-Seq. We found that both CR and HFD induce extensive transcriptional changes, in some cases altering the same genes in the same direction. We used our epigenomic data to infer transcriptional regulatory proteins bound near these genes that likely influence their expression levels. In particular, we found evidence for critical roles played by PPARα and RXRα. We used ChIP-Seq to profile the binding locations for these factors in HFD and CR livers. We found extensive binding of PPARα near genes involved in glycolysis/gluconeogenesis and uncovered a role for this factor in regulating anaerobic glycolysis. Overall, we generated extensive transcriptional and epigenomic datasets from livers of mice fed these diets and uncovered new functions and gene targets for PPARα.

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.

Mouse genotypes drive the liver and adrenal gland clocks

Circadian rhythms regulate a plethora of physiological processes. Perturbations of the rhythm can result in pathologies which are frequently studied in inbred mouse strains. We show that the genotype of mouse lines defines the circadian gene expression patterns. Expression of majority of core clock and output metabolic genes are phase delayed in the C56BL/6J line compared to 129S2 in the adrenal glands and the liver. Circadian amplitudes are generally higher in the 129S2 line. Experiments in dark - dark (DD) and light - dark conditions (LD), exome sequencing and data mining proposed that mouse lines differ in single nucleotide variants in the binding regions of clock related transcription factors in open chromatin regions. A possible mechanisms of differential circadian expression could be the entrainment and transmission of the light signal to peripheral organs. This is supported by the genotype effect in adrenal glands that is largest under LD, and by the high number of single nucleotide variants in the Receptor, Kinase and G-protein coupled receptor Panther molecular function categories. Different phenotypes of the two mouse lines and changed amino acid sequence of the Period 2 protein possibly contribute further to the observed differences in circadian gene expression.

Sexually Dimorphic Expression of eGFP Transgene in the Akr1A1 Locus of Mouse Liver Regulated by Sex Hormone-Related Epigenetic Remodeling

Sexually dimorphic gene expression is commonly found in the liver, and many of these genes are linked to different incidences of liver diseases between sexes. However, the mechanism of sexually dimorphic expression is still not fully understood. In this study, a pCAG-eGFP transgenic mouse strain with a specific transgene integration site in the Akr1A1 locus presented male-biased EGFP expression in the liver, and the expression was activated by testosterone during puberty. The integration of the pCAG-eGFP transgene altered the epigenetic regulation of the adjacent chromatin, including increased binding of STAT5b, a sexually dimorphic expression regulator, and the transformation of DNA methylation from hypermethylation into male-biased hypomethylation. Through this de novo sexually dimorphic expression of the transgene, the Akr1A1^eGFP mouse provides a useful model to study the mechanisms and the dynamic changes of sexually dimorphic gene expression during either development or pathogenesis of the liver.

Disruption of STAT5b-Regulated Sexual Dimorphism of the Liver Transcriptome by Diverse Factors Is a Common Event

Signal transducer and activator of transcription 5b (STAT5b) is a growth hormone (GH)-activated transcription factor and a master regulator of sexually dimorphic gene expression in the liver. Disruption of the GH hypothalamo-pituitary-liver axis controlling STAT5b activation can lead to metabolic dysregulation, steatosis, and liver cancer. Computational approaches were developed to identify factors that disrupt STAT5b function in a mouse liver gene expression compendium. A biomarker comprised of 144 STAT5b-dependent genes was derived using comparisons between wild-type male and wild-type female mice and between STAT5b-null and wild-type mice. Correlations between the STAT5b biomarker gene set and a test set comprised of expression datasets (biosets) with known effects on STAT5b function were evaluated using a rank-based test (the Running Fisher algorithm). Using a similarity p-value ≤ 10⁻⁴, the test achieved a balanced accuracy of 99% and 97% for detection of STAT5b activation or STAT5b suppression, respectively. The STAT5b biomarker gene set was then used to identify factors that activate (masculinize) or suppress (feminize) STAT5b function in an annotated mouse liver and primary hepatocyte gene expression compendium of ~1,850 datasets. Disruption of GH-regulated STAT5b is a common phenomenon in liver in vivo, with 5% and 29% of the male datasets, and 11% and 13% of the female datasets, associated with masculinization or feminization, respectively. As expected, liver STAT5b activation/masculinization occurred at puberty and suppression/feminization occurred during aging and in mutant mice with defects in GH signaling. A total of 70 genes were identified that have effects on STAT5b activation in genetic models in which the gene was inactivated or overexpressed. Other factors that affected liver STAT5b function were shown to include fasting, caloric restriction and infections. Together, these findings identify diverse factors that perturb the hypothalamo-pituitary-liver GH axis and disrupt GH-dependent STAT5b activation in mouse liver.

A Genome-Wide Perspective on Metabolism

Mammals have at least 210 histologically diverse cell types (Alberts, Molecular biology of the cell. Garland Science, New York, 2008) and the number would be even higher if functional differences are taken into account. The genome in each of these cell types is differentially programmed to express the specific set of genes needed to fulfill the phenotypical requirements of the cell. Furthermore, in each of these cell types, the gene program can be differentially modulated by exposure to external signals such as hormones or nutrients. The basis for the distinct gene programs relies on cell type-selective activation of transcriptional enhancers, which in turn are particularly sensitive to modulation. Until recently we had only fragmented insight into the regulation of a few of these enhancers; however, the recent advances in high-throughput sequencing technologies have enabled the development of a large number of technologies that can be used to obtain genome-wide insight into how genomes are reprogrammed during development and in response to specific external signals. By applying such technologies, we have begun to reveal the cross-talk between metabolism and the genome, i.e., how genomes are reprogrammed in response to metabolites, and how the regulation of metabolic networks is coordinated at the genomic level.