Sex-specific mouse liver gene expression: genome-wide analysis of developmental changes from pre-pubertal period to young adulthood

Genetic hypogonadal (Gnrh1hpg) mouse model uncovers influence of reproductive axis on maturation of the gut microbiome during puberty

The gut microbiome plays a key role in human health and gut dysbiosis is linked to many sex-specific diseases including autoimmune, metabolic, and neurological disorders. Activation of the hypothalamic-pituitary-gonadal (HPG) axis during puberty leads to sexual maturation and development of sex differences through the action of gonadal sex steroids. While the gut microbiome also undergoes sex differentiation, the mechanisms involved remain poorly understood. Using a genetic hypogonadal (hpg) mouse model, we sampled the fecal microbiome of male and female wild-type and hpg mutant mice before and after puberty to determine how microbial taxonomy and function are influenced by age, sex, and the HPG axis. We showed that HPG axis activation during puberty is required for sexual maturation of the gut microbiota composition, community structure, and metabolic functions. We also demonstrated that some sex differences in taxonomic composition and amine metabolism developed independently of the HPG axis, indicating that sex chromosomes are sufficient for certain sex differences in the gut microbiome. In addition, we showed that age, independent of HPG axis activation, led to some aspects of pubertal maturation of the gut microbiota community composition and putative functions. These results have implications for microbiome-based treatments, indicating that sex, hormonal status, and age should be considered when designing microbiome-based therapeutics.

Endo-reduplication in mouse liver after conditional mutation of ORC2 and combined mutation of ORC1 and ORC2

The six subunit Origin Recognition Complex (ORC) is essential for loading MCM2-7 at origins of DNA replication to promote initiation of DNA replication in organisms ranging from S. cerevisiae to humans. In rare instances, as in cancer cell-lines in culture with mutations in ORC1 , ORC2 or ORC5 , or in endo-reduplicating mouse hepatocytes in vivo without ORC1 , DNA replication has been observed in the virtual absence of individual ORC subunits. Although ORC1 is dispensable in the mouse liver for endo-reduplication, because of the homology of ORC1 with CDC6, it could be argued that CDC6 was substituting for ORC1 to restore functional ORC. Here, we have created mice with a conditional deletion of ORC2 , to demonstrate that mouse embryo fibroblasts require ORC2 for proliferation, but that the mouse hepatocytes can carry out DNA synthesis in vitro and endo-reduplicate in vivo , despite the deletion of ORC2 . Combining the conditional mutation of ORC1 and ORC2 revealed that the mouse liver can still carry out endo-reduplication despite the deletion of the two genes, both during normal development and after partial hepatectomy. Since endo-reduplication, like normal S phase replication, requires the presence of MCM2-7 on the chromatin, these results suggest that in primary hepatocytes there is a mechanism to load sufficient MCM2-7 to carry out effective DNA replication despite the virtual absence of two subunits of ORC.

Pharmacological activation of constitutive androstane receptor induces female-specific modulation of hepatic metabolism

Background & Aims

The constitutive androstane receptor (CAR) is a nuclear receptor that binds diverse xenobiotics and whose activation leads to the modulation of the expression of target genes involved in xenobiotic detoxification and energy metabolism. Although CAR hepatic activity is considered to be higher in women than in men, its sex-dependent response to an acute pharmacological activation has seldom been investigated.

Methods

The hepatic transcriptome, plasma markers, and hepatic metabolome, were analysed in Car+/+ and Car-/- male and female mice treated either with the CAR-specific agonist 1,4-bis[2-(3,5-dichloropyridyloxy)]benzene (TCPOBOP) or with vehicle.

Results

Although 90% of TCPOBOP-sensitive genes were modulated in a sex-independent manner, the remaining 10% showed almost exclusive female liver specificity. These female-specific CAR-sensitive genes were mainly involved in xenobiotic metabolism, inflammation, and extracellular matrix organisation. CAR activation also induced higher hepatic oxidative stress and hepatocyte cytolysis in females than in males. Hepatic expression of flavin monooxygenase 3 (Fmo3) was almost abolished and was associated with a decrease in hepatic trimethylamine-N-oxide (TMAO) concentration in TCPOBOP-treated females. In line with a potential role in the control of TMAO homeostasis, CAR activation decreased platelet hyper-responsiveness in female mice supplemented with dietary choline.

Conclusions

More than 10% of CAR-sensitive genes are sex-specific and influence hepatic and systemic responses such as platelet aggregation. CAR activation may be an important mechanism of sexually-dimorphic drug-induced liver injury.

Impact and implications

CAR is activated by many drugs and pollutants. Its pharmacological activation had a stronger impact on hepatic gene expression and metabolism in females than in males, and had a specific impact on liver toxicity and trimethylamine metabolism. Sexual dimorphism should be considered when testing and/or prescribing xenobiotics known to activate CAR.

HNF4α isoforms regulate the circadian balance between carbohydrate and lipid metabolism in the liver

Hepatocyte Nuclear Factor 4α (HNF4α), a master regulator of hepatocyte differentiation, is regulated by two promoters (P1 and P2) which drive the expression of different isoforms. P1-HNF4α is the major isoform in the adult liver while P2-HNF4α is thought to be expressed only in fetal liver and liver cancer. Here, we show that P2-HNF4α is indeed expressed in the normal adult liver at Zeitgeber time (ZT)9 and ZT21. Using exon swap mice that express only P2-HNF4α we show that this isoform orchestrates a distinct transcriptome and metabolome via unique chromatin and protein-protein interactions, including with different clock proteins at different times of the day leading to subtle differences in circadian gene regulation. Furthermore, deletion of the Clock gene alters the circadian oscillation of P2- (but not P1-)HNF4α RNA, revealing a complex feedback loop between the HNF4α isoforms and the hepatic clock. Finally, we demonstrate that while P1-HNF4α drives gluconeogenesis, P2-HNF4α drives ketogenesis and is required for elevated levels of ketone bodies in female mice. Taken together, we propose that the highly conserved two-promoter structure of the Hnf4a gene is an evolutionarily conserved mechanism to maintain the balance between gluconeogenesis and ketogenesis in the liver in a circadian fashion.

Corticosteroid-Binding Globulin (SERPINA6) Consolidates Sexual Dimorphism of Adult Rat Liver

Produced by the liver, corticosteroid-binding globulin (CBG) regulates the plasma distribution and actions of glucocorticoids. A sex difference in pituitary growth hormone secretion patterns established during puberty in rats results in increased hepatic CBG production and 2-fold higher plasma corticosterone levels in females. Glucocorticoids control hepatic development and metabolic activities, and we have therefore examined how disrupting the SerpinA6 gene encoding CBG influences plasma corticosterone dynamics, as well as liver gene expression in male and female rats before and after puberty. Comparisons of corticosterone plasma clearance and hepatic uptake in adult rats, with or without CBG, indicated that CBG limits corticosterone clearance by reducing its hepatic uptake. Hepatic transcriptomic profiling revealed minor sex differences (207 differentially expressed genes) and minimal effect of CBG deficiency in 30-day-old rats before puberty. While liver transcriptomes in 60-day-old males lacking CBG remained essentially unchanged, 2710 genes were differentially expressed in wild-type female vs male livers at this age. Importantly, ∼10% of these genes lost their sexually dimorphic expression in adult females lacking CBG, including those related to cholesterol biosynthesis, inflammation, and lipid and amino acid catabolism. Another 203 genes were altered by the loss of CBG specifically in adult females, including those related to xenobiotic metabolism, circadian rhythm, and gluconeogenesis. Our findings reveal that CBG consolidates the sexual dimorphism of the rat liver initiated by sex differences in growth hormone secretion patterns and provide insight into how CBG deficiencies are linked to glucocorticoid-dependent diseases.

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/ .

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

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

NUDT7 regulates total hepatic CoA levels and the composition of the intestinal bile acid pool in male mice fed a Western diet

Nudix hydrolase 7 (NUDT7) is an enzyme that hydrolyzes CoA species, is highly expressed in the liver, and resides in the peroxisomes. Peroxisomes are organelles where the preferential oxidation of dicarboxylic fatty acids occurs and where the hepatic synthesis of the primary bile acids cholic acid and chenodeoxycholic acid is completed. We previously showed that liver-specific overexpression of NUDT7 affects peroxisomal lipid metabolism but does not prevent the increase in total liver CoA levels that occurs during fasting. We generated Nudt7^-/- mice to further characterize the role that peroxisomal (acyl-)CoA degradation plays in the modulation of the size and composition of the acyl-CoA pool and in the regulation of hepatic lipid metabolism. Here, we show that deletion of Nudt7 alters the composition of the hepatic acyl-CoA pool in mice fed a low-fat diet, but only in males fed a Western diet does the lack of NUDT7 activity increase total liver CoA levels. This effect is driven by the male-specific accumulation of medium-chain dicarboxylic acyl-CoAs, which are produced from the β-oxidation of dicarboxylic fatty acids. We also show that, under conditions of elevated synthesis of chenodeoxycholic acid derivatives, Nudt7 deletion promotes the production of tauromuricholic acid, decreasing the hydrophobicity index of the intestinal bile acid pool and increasing fecal cholesterol excretion in male mice. These findings reveal that NUDT7-mediated hydrolysis of acyl-CoA pathway intermediates in liver peroxisomes contributes to the regulation of dicarboxylic fatty acid metabolism and the composition of the bile acid pool.

Ontogenetic Change in Male Expression of Testosterone-Responsive Genes Contributes to the Emergence of Sex-Biased Gene Expression in Anolis sagrei

Sex differences in gene expression tend to increase with age across a variety of species, often coincident with the development of sexual dimorphism and maturational changes in hormone levels. However, because most transcriptome-wide characterizations of sexual divergence are framed as comparisons of sex-biased gene expression across ages, it can be difficult to determine the extent to which age-biased gene expression within each sex contributes to the emergence of sex-biased gene expression. Using RNAseq in the liver of the sexually dimorphic brown anole lizard (Anolis sagrei), we found that a pronounced increase in sex-biased gene expression with age was associated with a much greater degree of age-biased gene expression in males than in females. This pattern suggests that developmental changes in males, such as maturational increases in circulating testosterone, contribute disproportionately to the ontogenetic emergence of sex-biased gene expression. To test this hypothesis, we used four different experimental contrasts to independently characterize sets of genes whose expression differed as a function of castration and/or treatment with exogenous testosterone. We found that genes that were significantly male-biased in expression or upregulated as males matured tended to be upregulated by testosterone, whereas genes that were female-biased or downregulated as males matured tended to be downregulated by testosterone. Moreover, the first two principal components describing multivariate gene expression indicated that exogenous testosterone reversed many of the feminizing effects of castration on the liver transcriptome of maturing males. Collectively, our results suggest that developmental changes that occur in males contribute disproportionately to the emergence of sex-biased gene expression in the Anolis liver, and that many of these changes are orchestrated by androgens such as testosterone.

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.

Maternal obesity (MO) programs morphological changes in aged rat offspring small intestine in a sex dependent manner: Effects of maternal resveratrol supplementation

Maternal obesity (MO) leads to offspring metabolic problems. The mechanisms involved are multifactorial. The small intestine plays an important role in the absorption of nutrients and is modified as we age. Few studies have explored MO programming effects on offspring (F1) small intestine morphology. The aim of this study was to investigate MO effects on old adult F1 intestinal morphology, and whether any F1 intestinal changes due to MO were modified by maternal resveratrol supplementation. From weaning throughout pregnancy and lactation, female Wistar rats (F0) ate standard chow (controls, C: 5%-fat) or high-fat diet (MO: 25%-fat). One month before mating at postnatal day (PND) 120 through lactation half of each group received 20 mg/kg/day of resveratrol orally (Cres or MOres). After weaning F1 were fed with chow diet until the end of the study at PND 650. Body weight, percent of fat, glucose, cholesterol and triglyceride serum concentrations were determined. F1 small intestinal samples were collected for histological analysis. Male F1 body weight was higher in MO and MOres compared with C and Cres. Female F1 body weight and percent of fat was higher in MO than C and MOres. Triglyceride concentrations were higher in MO and MOres male F1 compared with C and Cres. There were no differences among groups in female triglyceride concentrations. Male F1 duodenal villus height was smaller in MO compared with MOres. Female F1 duodenal and jejunal crypt depth was smaller in MO compared with C and was greater compared with MOres. Female F1 villus height in jejunum was greater in MO compared with MOres. In conclusion, exposure to the developmental challenge of MO changed the aged F1 intestinal morphological and metabolic profiles. Maternal resveratrol supplementation ameliorated these effects in an F1 sex dependent manner.

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

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

Global analysis of expression, maturation and subcellular localization of mouse liver transcriptome identifies novel sex-biased and TCPOBOP-responsive long non-coding RNAs

BACKGROUND

While nuclear transcription and RNA processing and localization are well established for protein coding genes (PCGs), these processes are poorly understood for long non-coding (lnc)RNAs. Here, we characterize global patterns of transcript expression, maturation and localization for mouse liver RNA, including more than 15,000 lncRNAs. PolyA-selected liver RNA was isolated and sequenced from four subcellular fractions (chromatin, nucleoplasm, total nucleus, and cytoplasm), and from the chromatin-bound fraction without polyA selection.

RESULTS

Transcript processing, determined from normalized intronic to exonic sequence read density ratios, progressively increased for PCG transcripts in going from the chromatin-bound fraction to the nucleoplasm and then on to the cytoplasm. Transcript maturation was similar for lncRNAs in the chromatin fraction, but was significantly lower in the nucleoplasm and cytoplasm. LncRNA transcripts were 11-fold more likely to be significantly enriched in the nucleus than cytoplasm, and 100-fold more likely to be significantly chromatin-bound than nucleoplasmic. Sequencing chromatin-bound RNA greatly increased the sensitivity for detecting lowly expressed lncRNAs and enabled us to discover and localize hundreds of novel regulated liver lncRNAs, including lncRNAs showing sex-biased expression or responsiveness to TCPOBOP a xenobiotic agonist ligand of constitutive androstane receptor (Nr1i3).

CONCLUSIONS

Integration of our findings with prior studies and lncRNA annotations identified candidate regulatory lncRNAs for a variety of hepatic functions based on gene co-localization within topologically associating domains or transcription divergent or antisense to PCGs associated with pathways linked to hepatic physiology and disease.

Embryonic Onset of Sexually Dimorphic Heart Rates in the Viviparous Fish, Gambusia holbrooki

In fish, little is known about sex-specific differences in physiology and performance of the heart and whether these differences manifest during development. Here for the first time, the sex-specific heart rates during embryogenesis of Gambusia holbrooki, from the onset of the heart rates (HRs) to just prior to parturition, was investigated using light cardiogram. The genetic sex of the embryos was post-verified using a sex-specific genetic marker. Results reveal that heart rates and resting time significantly increase (p < 0.05) with progressive embryonic development. Furthermore, both ventricular and atrial frequencies of female embryos were significantly higher (p < 0.05) than those of their male sibs at the corresponding developmental stages and remained so at all later developmental stages (p < 0.05). In concurrence, the heart rate and ventricular size of the adult females were also significantly (p < 0.05) higher and larger respectively than those of males. Collectively, the results suggest that the cardiac sex-dimorphism manifests as early as late-organogenesis and persists through adulthood in this species. These findings suggest that the cardiac measurements can be employed to non-invasively sex the developing embryos, well in advance of when their phenotypic sex is discernible. In addition, G. holbrooki could serve as a better model to study comparative vertebrate cardiovascular development as well as to investigate anthropogenic and climatic impacts on heart physiology of this species, that may be sex influenced.

Sex-biased DNA methylation in papillary thyroid cancer

Background: Women have a higher risk of developing papillary thyroid cancer (PTC) than men. DNA methylation is known to differ between men and women. Materials & methods: Using the human methylation 450 BeadChip and RNA-sequence, we profiled the genome-wide DNA methylation patterns of papillary thyroid carcinoma patients and para-carcinoma tissue. Results: We first identified 398 different expression genes (DEGs) between males and females PTC. Then we analyzed the relationship between differentially methylated positions (DMPs) and DEGs at gene promoter regions and identified 39 genes and explored DMP-DEGs' correlation with immune cell infiltration and tumor purity. We also analyzed the relationship between genomic regions and enhancers. Conclusion: Our study identified 39 DMP-DEGs providing some new insights into the mechanisms of methylation-mediated gender differences in PTC.

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.

Tissue and sex-specific programming of DNA methylation by perinatal lead exposure: implications for environmental epigenetics studies

Early developmental environment can influence long-term health through reprogramming of the epigenome. Human environmental epigenetics studies rely on surrogate tissues, such as blood, to assess the effects of environment on disease-relevant but inaccessible target tissues. However, the extent to which environment-induced epigenetic changes are conserved between these tissues is unclear. A better understanding of this conservation is imperative for effective design and interpretation of human environmental epigenetics studies. The Toxicant Exposures and Responses by Genomic and Epigenomic Regulators of Transcription (TaRGET II) consortium was established by the National Institute of Environmental Health Sciences to address the utility of surrogate tissues as proxies for toxicant-induced epigenetic changes in target tissues. We and others have recently reported that perinatal exposure to lead (Pb) is associated with adverse metabolic outcomes. Here, we investigated the sex-specific effects of perinatal exposure to a human environmentally relevant level of Pb on DNA methylation in paired liver and blood samples from adult mice using enhanced reduced-representation bisulphite sequencing. Although Pb exposure ceased at 3 weeks of age, we observed thousands of sex-specific differentially methylated cytosines in the blood and liver of Pb-exposed animals at 5 months of age, including 44 genomically imprinted loci. We observed significant tissue overlap in the genes mapping to differentially methylated cytosines. A small but significant subset of Pb-altered genes exhibit basal sex differences in gene expression in the mouse liver. Collectively, these data identify potential molecular targets for Pb-induced metabolic diseases, and inform the design of more robust human environmental epigenomics studies.

Effect of endurance exercise training on liver gene expression in male and female mice

Chronic endurance exercise is a therapeutic strategy in the treatment of many chronic diseases in humans, including the prevention and treatment of metabolic diseases such as diabetes mellitus. Metabolic, cardiorespiratory, and endocrine pathways targeted by chronic endurance exercise have been identified. In the liver, however, the cellular and molecular pathways that are modified by exercise and have preventive or therapeutic relevance to metabolic disease need to be elucidated. The mouse model used in the current study allows for the quantification of a human-relevant exercise "dosage". In this study we show hepatic gene expression differences between sedentary female and sedentary male mice and that chronic exercise modifies the transcription of hepatic genes related to metabolic disease and steatosis in both male and female mice. Chronic exercise induces molecular pathways involved in glucose tolerance, glycolysis, and gluconeogenesis while producing a decrease in pathways related to insulin resistance, steatosis, fibrosis, and inflammation. Given these findings, this mouse exercise model has potential to dissect the cellular and molecular hepatic changes following chronic exercise with application to understanding the role that chronic exercise plays in preventing human diseases. Novelty: Exercise modifies the hepatic gene expression and hepatic pathways related to metabolic disease in male and female mice. Sex differences were seen in hepatic gene expression between sedentary and exercised mice. The mouse exercise model used in this study allows for application and evaluation of exercise effects in human 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.

A reporter mouse for non-invasive detection of toll-like receptor ligands induced acute phase responses

The acute phase response (APR) is a systemic first-line defense against challenges including infection, trauma, stress, and neoplasia. Alteration of acute phase protein (APP) levels in plasma is the most important change during acute phase response. C-reactive protein (CRP), which increases dramatically during inflammation onset, is an indicator of inflammation. To monitor the process of APR, we generated human CRP promoter-driven luciferase transgenic (hCRP-Luc) mice to quantify the hCRP promoter activation in vivo. The naïve female hCRP-Luc mice express low basal levels of liver bioluminescence, but the naïve male hCRP-Luc mice do not. Thus, female hCRP-Luc mice are suitable for monitoring the process of APR. The liver bioluminescence of female hCRP-Luc mice can be induced by several toll-like receptor (TLR) ligands. The expression of liver bioluminescence was highly sensitive to endotoxin stimulation in a dose-dependent manner. On-off-on bioluminescence response was noted in female hCRP-Luc mice upon two endotoxin stimulations one month apart. The LPS-induced bioluminescence of the female hCRP-Luc mice was IL-6-mediated and associated with APP alpha-1-acid glycoprotein expression. In conclusion, the female hCRP-Luc mouse is a non-invasive, sensitive and reusable reporter tool for APR.

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.

Inducible Loss of the Aryl Hydrocarbon Receptor Activates Perigonadal White Fat Respiration and Brown Fat Thermogenesis via Fibroblast Growth Factor 21

The aryl hydrocarbon receptor (AHR) is a ligand-activated transcription factor highly expressed in hepatocytes. Researchers have employed global and liver-specific conditional Ahr knockout mouse models to characterize the physiological roles of the AHR; however, the gestational timing of AHR loss in these models can complicate efforts to distinguish the direct and indirect effects of post-gestational AHR deficiency. Utilizing a novel tamoxifen-inducible AHR knockout mouse model, we analyzed the effects of hepatocyte-targeted AHR loss in adult mice. The data demonstrate that AHR deficiency significantly reduces weight gain and adiposity, and increases multilocular lipid droplet formation within perigonadal white adipose tissue (gWAT). Protein and mRNA expression of fibroblast growth factor 21 (FGF21), an important hepatokine that activates thermogenesis in brown adipose tissue (BAT) and gWAT, significantly increases upon AHR loss and correlates with a significant increase of BAT and gWAT respiratory capacity. Confirming the role of FGF21 in mediating these effects, this phenotype is reversed in mice concomitantly lacking AHR and FGF21 expression. Chromatin immunoprecipitation analyses suggest that the AHR may constitutively suppress Fgf21 transcription through binding to a newly identified xenobiotic response element within the Fgf21 promoter. The data demonstrate an important AHR-FGF21 regulatory axis that influences adipose biology and may represent a “druggable” therapeutic target for obesity and its related metabolic disorders.

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.

Establishment and analysis of a mouse model that regulates sex-related differences in liver drug metabolism

The adult liver performs many metabolic functions for maintaining homeostasis. There are several sex differences in liver function and disease pathogenesis. One important function of the liver is drug metabolism, where cytochrome p450s (CYPs) in hepatocytes are the main enzymes involved. The toxicity of various drugs and chemicals differs with sex due to differences in hepatocytic CYP expression. However, the molecular mechanism regulating sex-related differences in drug metabolism remains unknown. In this study, we identified transcriptional regulator B-cell lymphoma 6 (Bcl6) as an important factor in sex-biased differential CYP expression. Microarray analysis of livers derived from liver-specific Bcl6-knockout mice showed that Bcl6 is required for sex-biased CYP expression patterns in the liver. Additionally, quantitative PCR analysis revealed that hepatocytic expression of male-biased genes, such as Cyp2d9, Cyp2u1, Cyp4a12a/12b, and Cyp7b1, in liver-specific Bcl6-knockout male mice significantly decreased to levels similar to those observed in wild-type female mice. Conversely, hepatocytic expression of female-biased genes, such as Cyp2a4/2a5, Cyp2b9, Cyp3a41, and Cyp17a1, significantly increased in liver-specific Bcl6-knockout male mice. Deletion of Bcl6 caused female-like expression of CYPs in male livers. These results suggest that Bcl6 is a key regulator of sex-related differential regulation of drug metabolism. Moreover, serum sex hormone levels and fertility did not change in liver-specific, Bcl6-knockout mice. Hepatocytic Bcl6 regulates sex-related differential CYP expression in the liver without changing the sex of the whole body. Thus, this mouse model is useful for analyzing liver-specific sex-dependent regulation of drug metabolism and pathogenesis.

Genomics of sex hormone receptor signaling in hepatic sexual dimorphism

The liver plays a crucial role in a variety of physiological processes. Sexual dimorphism is markedly defined in liver disorders, such as fatty liver diseases and liver cancer, but barely addressed in the normal liver. Distinct sex hormone signaling between male and female livers is the major driving factor for hepatic sexual dimorphism. Over 6000 genes are differently expressed between male and female livers in mice. Here we address how sex hormone receptors, estrogen receptor alpha (ERα) and androgen receptor (AR), mediate sexually dimorphic gene expression in mouse livers. We identified 5192 ERα target genes and 4154 AR target genes using ChIP-Seq. Using liver-specific ERα or AR knockout mice, we further identified direct and functional target genes of ERα (123 genes) and AR (151 genes) that contribute to hepatic sexual dimorphism. We also found that the most significant sexually dimorphic gene expression was initiated at birth by comparing hepatic gene expression data from the embryonic stage E10.5 to the postnatal stage P60 during liver development. Overall, our study indicates that sex hormone receptor signaling drives sexual dimorphism of hepatic gene expression throughout liver development.

The molecular pathways underlying early gonadal development

The body of knowledge surrounding reproductive development spans the fields of genetics, anatomy, physiology and biomedicine, to build a comprehensive understanding of the later stages of reproductive development in humans and animal models. Despite this, there remains much to learn about the bi-potential progenitor structure that the ovary and testis arise from, known as the genital ridge (GR). This tissue forms relatively late in embryonic development and has the potential to form either the ovary or testis, which in turn produce hormones required for development of the rest of the reproductive tract. It is imperative that we understand the genetic networks underpinning GR development if we are to begin to understand abnormalities in the adult. This is particularly relevant in the contexts of disorders of sex development (DSDs) and infertility, two conditions that many individuals struggle with worldwide, with often no answers as to their aetiology. Here, we review what is known about the genetics of GR development. Investigating the genetic networks required for GR formation will not only contribute to our understanding of the genetic regulation of reproductive development, it may in turn open new avenues of investigation into reproductive abnormalities and later fertility issues in the adult.

Evaluating mice lacking serum carboxylesterase as a behavioral model for nerve agent intoxication

Mice and other rodents are typically utilized for chemical warfare nerve agent research. Rodents have large amounts of carboxylesterase in their blood, while humans do not. Carboxylesterase nonspecifically binds to and detoxifies nerve agent. The presence of this natural bioscavenger makes mice and other rodents poor models for studies identifying therapeutics to treat humans exposed to nerve agents. To obviate this problem, a serum carboxylesterase knockout (Es1 KO) mouse was created. In this study, Es1 KO and wild type (WT) mice were assessed for differences in gene expression, nerve agent (soman; GD) median lethal dose (MLD) values, and behavior prior to and following nerve agent exposure. No expression differences were detected between Es1 KO and WT mice in more than 34 000 mouse genes tested. There was a significant difference between Es1 KO and WT mice in MLD values, as the MLD for GD-exposed WT mice was significantly higher than the MLD for GD-exposed Es1 KO mice. Behavioral assessments of Es1 KO and WT mice included an open field test, a zero maze, a Barnes maze, and a sucrose preference test (SPT). While sex differences were observed in various measures of these tests, overall, Es1 KO mice behaved similarly to WT mice. The two genotypes also showed virtually identical neuropathological changes following GD exposure. Es1 KO mice appear to have an enhanced susceptibility to GD toxicity while retaining all other behavioral and physiological responses to this nerve agent, making the Es1 KO mouse a more human-like model for nerve agent research.

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.

Gene expression profiling of puberty-associated genes reveals abundant tissue and sex-specific changes across postnatal development

The timing of human puberty is highly variable, sexually dimorphic, and associated with adverse health outcomes. Over 20 genes carrying rare mutations have been identified in known pubertal disorders, many of which encode critical components of the hypothalamic-pituitary-gonadal (HPG) axis. Recent genome-wide association studies (GWAS) have identified more than 100 candidate genes at loci associated with age at menarche or voice breaking in males. We know little about the spatial, temporal or postnatal expression patterns of the majority of these puberty-associated genes. Using a high-throughput and sensitive microfluidic quantitative PCR strategy, we profiled the gene expression patterns of the mouse orthologs of 178 puberty-associated genes in male and female mouse HPG axis tissues, the pineal gland, and the liver at five postnatal ages spanning the pubertal transition. The most dynamic gene expression changes were observed prior to puberty in all tissues. We detected known and novel tissue-enhanced gene expression patterns, with the hypothalamus expressing the largest number of the puberty-associated genes. Notably, over 40 puberty-associated genes in the pituitary gland showed sex-biased gene expression, most of which occurred peri-puberty. These sex-biased genes included the orthologs of candidate genes at GWAS loci that show sex-discordant effects on pubertal timing. Our findings provide new insight into the expression of puberty-associated genes and support the possibility that the pituitary plays a role in determining sex differences in the timing of puberty.

Catsper1 promoter is bidirectional and regulates the expression of a novel lncRNA

The Catsper1 gene, whose expression is restricted to male germ cells, has great importance in reproductive biology because of its function in sperm motility and fertilization. We previously reported that the promoter of this gene has transcriptional activity in either direction in a heterologous system. In the present study, we found that the Catsper1 promoter has in vitro transcriptional activity in either orientation in GC-1 spg mouse spermatogonial cells. The results also showed that this promoter regulates the expression of a new divergent Catsper1 gene named Catsper1au (Catsper1 antisense upstream transcript). Catsper1au is expressed in adult male mouse testis and liver tissues but not in female mouse liver or ovary tissues. In the testis, Catsper1au is expressed in embryos at 11.5 days post-coitum and from newborns to adults. This gene is also expressed in 1- to 3-week postnatal hearts and in 1-week to adult stage livers. The analysis of the 1402 bp whole genome sequence revealed that Catsper1au is an intronless and polyadenylated lncRNA, located in the nuclei of Sertoli and spermatogenic cells from adult testis. These data indicate that Catsper1au is divergently expressed from the Catsper1 promoter and could regulate gene expression during spermatogenesis.

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.

Macro- and microelements in the rat liver, kidneys, and brain tissues; sex differences and effect of blood removal by perfusion in vivo

Concentrations of macro- and microelements in animal organs indicate the animal health status and represent reference data for animal experiments. Their levels in blood and tissues could be different between sexes, and could be different with and without blood in tissues. To test these hypotheses, in adult female and male rats the concentrations of various elements were measured in whole blood, blood plasma, and tissues from blood-containing (nonperfused) and blood-free liver, kidneys, and brain (perfused in vivo with an elements-free buffer). In these samples, 6 macroelements (Na, Mg, P, S, K, Ca) and 14 microelements (Fe, Mn, Co, Cu, Zn, Se, I, As, Cd, Hg, Pb, Li, B, Sr) were determined by inductively coupled plasma mass spectrometry following nitric acid digestion. In blood and plasma, female- or male-dominant sex differences were observed for 6 and 5 elements, respectively. In nonperfused organs, sex differences were observed for 3 (liver, brain) or 9 (kidneys) elements, whereas in perfused organs, similar differences were detected for 9 elements in the liver, 5 in the kidneys, and none in the brain. In females, perfused organs had significantly lower concentrations of 4, 5, and 2, and higher concentrations of 10, 4, and 7 elements, respectively, in the liver, kidneys, and brain. In males, perfusion caused lower concentrations of 4, 7, and 2, and higher concentrations of 1, 1, and 7 elements, respectively, in the liver, kidneys, and brain. Therefore, the residual blood in organs can significantly influence tissue concentrations of various elements and their sex-dependency.

Disrupting Hepatocyte Cyp51 from Cholesterol Synthesis Leads to Progressive Liver Injury in the Developing Mouse and Decreases RORC Signalling

Development of mice with hepatocyte knockout of lanosterol 14α-demethylase (H^Cyp51−/−) from cholesterol synthesis is characterized by the progressive onset of liver injury with ductular reaction and fibrosis. These changes begin during puberty and are generally more aggravated in the knockout females. However, a subgroup of (pre)pubertal knockout mice (runts) exhibits a pronounced male prevalent liver dysfunction characterized by downregulated amino acid metabolism and elevated Casp12. RORC transcriptional activity is diminished in livers of all runt mice, in correlation with the depletion of potential RORC ligands subsequent to CYP51 disruption. Further evidence for this comes from the global analysis that identified a crucial overlap between hepatic Cyp51^−/− and Rorc^−/− expression profiles. Additionally, the reduction in RORA and RORC transcriptional activity was greater in adult H^Cyp51−/− females than males, which correlates well with their downregulated amino and fatty acid metabolism. Overall, we identify a global and sex-dependent transcriptional de-regulation due to the block in cholesterol synthesis during development of the Cyp51 knockout mice and provide in vivo evidence that sterol intermediates downstream of lanosterol may regulate the hepatic RORC activity.

Metabolic Profiling Reveals Effects of Age, Sexual Development and Neutering in Plasma of Young Male Cats

Neutering is a significant risk factor for obesity in cats. The mechanisms that promote neuter-associated weight gain are not well understood but following neutering, acute changes in energy expenditure and energy consumption have been observed. Metabolic profiling (GC-MS and UHPLC-MS-MS) was used in a longitudinal study to identify changes associated with age, sexual development and neutering in male cats fed a nutritionally-complete dry diet to maintain an ideal body condition score. At eight time points, between 19 and 52 weeks of age, fasted blood samples were taken from kittens neutered at either 19 weeks of age (Early Neuter (EN), n = 8) or at 31 weeks of age (Conventional Neuter (CN), n = 7). Univariate and multivariate analyses were used to compare plasma metabolites (n = 370) from EN and CN cats. Age was the primary driver of variance in the plasma metabolome, including a developmental change independent of neuter group between 19 and 21 weeks in lysolipids and fatty acid amides. Changes associated with sexual development and its subsequent loss were also observed, with differences at some time points observed between EN and CN cats for 45 metabolites (FDR p<0.05). Pathway Enrichment Analysis also identified significant effects in 20 pathways, dominated by amino acid, sterol and fatty acid metabolism. Most changes were interpretable within the context of male sexual development, and changed following neutering in the CN group. Felinine metabolism in CN cats was the most significantly altered pathway, increasing during sexual development and decreasing acutely following neutering. Felinine is a testosterone-regulated, felid-specific glutathione derivative secreted in urine. Alterations in tryptophan, histidine and tocopherol metabolism observed in peripubertal cats may be to support physiological functions of glutathione following diversion of S-amino acids for urinary felinine secretion.

Deep sexual dimorphism in adult medaka fish liver highlighted by multi-omic approach

Sexual dimorphism describes the features that discriminate between the two sexes at various biological levels. Especially, during the reproductive phase, the liver is one of the most sexually dimorphic organs, because of different metabolic demands between the two sexes. The liver is a key organ that plays fundamental roles in various physiological processes, including digestion, energetic metabolism, xenobiotic detoxification, biosynthesis of serum proteins, and also in endocrine or immune response. The sex-dimorphism of the liver is particularly obvious in oviparous animals, as the female liver is the main organ for the synthesis of oocyte constituents. In this work, we are interested in identifying molecular sexual dimorphism in the liver of adult medaka fish and their sex-variation in response to hepatotoxic exposures. By developing an integrative approach combining histology and different high-throughput omic investigations (metabolomics, proteomics and transcriptomics), we were able to globally depict the strong sexual dimorphism that concerns various cellular and molecular processes of hepatocytes comprising protein synthesis, amino acid, lipid and polysaccharide metabolism, along with steroidogenesis and detoxification. The results of this work imply noticeable repercussions on the biology of oviparous organisms environmentally exposed to chemical or toxin issues.

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.

Cyp1b1 affects external control of mouse hepatocytes, fatty acid homeostasis and signaling involving HNF4α and PPARα

Cytochrome P450 1b1 (Cyp1b1) is expressed in endothelia, stellate cells and pre-adipocytes, but not hepatocytes. Deletion alters liver fatty acid metabolism and prevents obesity and hepatic steatosis. This suggests a novel extra-hepatocyte regulation directed from cells that express Cyp1b1. To characterize these mechanisms, microarray gene expression was analyzed in livers of normal and congenic Cyp1b1-ko C57BL/6 J mice fed either low or high fat diets. Cyp1b1-ko gene responses indicate suppression of endogenous PPARα activity, a switch from triglyceride storage to mitochondrial fatty acid oxidation and decreased oxidative stress. Many gene responses in Cyp1b1-ko are sexually dimorphic and correspond to increased activity of growth hormone mediated by HNF4α. Male responses stimulated by GH pulses are enhanced, whereas responses that decline exhibit further suppression, including Cyp regulation by PPARα, CAR and PXR. These effects of Cyp1b1 deletion overlap with effects caused by deletion of the small heterodimeric partner, a suppressor of these nuclear factors. Redirection of gene expression associated with liver fat homeostasis in Cyp1b1-ko mice that directs hypothalamic control of GH and leptin. Cyp1b1-ko suppresses neonatal Scd1 and delays adult maturation of dimorphic GH/HNF4α signaling. Alternatively, deletion may diminish hypothalamic metabolism of estradiol, which establishes adult GH regulation.

Gender Differences in Response to Prolonged Every-Other-Day Feeding on the Proliferation and Apoptosis of Hepatocytes in Mice

Intermittent fasting decreases glucose and insulin levels and increases insulin sensitivity and lifespan. Decreased food intake influences the liver. Previous studies have shown gender differences in response to various types of caloric restriction, including every-other-day (EOD) feeding, in humans and rodents. Our goal was to show the influence of prolonged EOD feeding on the morphology, proliferation and apoptosis of livers from male and female mice. After nine months of an EOD diet, the livers from male and female mice were collected. We examined their morphology on histological slides using the Hematoxilin and Eosine (H_E) method and Hoechst staining of cell nuclei to evaluate the nuclear area of hepatocytes. We also evaluated the expression of mRNA for proto-oncogens, pro-survival proteins and apoptotic markers using Real Time Polimerase Chain Reaction (PCR). We noted increased lipid content in the livers of EOD fed female mice. EOD feeding lead to a decrease of proliferation and apoptosis in the livers of female and male mice, which suggest that tissue maintenance occurred during EOD feeding. Our experiment revealed sex-specific expression of mRNA for proto-oncogenes and pro-survival and pro-apoptotic genes in mice as well as sex-specific responses to the EOD treatment.

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(-4), 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.

Sexually dimorphic gene expression emerges with embryonic genome activation and is dynamic throughout development

Background

As sex determines mammalian development, understanding the nature and developmental dynamics of the sexually dimorphic transcriptome is important. To explore this, we generated 76 genome-wide RNA-seq profiles from mouse eight-cell embryos, late gestation and adult livers, together with 4 ground-state pluripotent embryonic (ES) cell lines from which we generated both RNA-seq and multiple ChIP-seq profiles. We complemented this with previously published data to yield 5 snap-shots of pre-implantation development, late-gestation placenta and somatic tissue and multiple adult tissues for integrative analysis.

Results

We define a high-confidence sex-dimorphic signature of 69 genes in eight-cell embryos. Sex-chromosome-linked components of this signature are largely conserved throughout pre-implantation development and in ES cells, whilst the autosomal component is more dynamic. Sex-biased gene expression is reflected by enrichment for activating and repressive histone modifications. The eight-cell signature is largely non-overlapping with that defined from fetal liver, neither was it correlated with adult liver or other tissues analysed. The number of sex-dimorphic genes increases throughout development. We identified many more dimorphic genes in adult compared to fetal liver. However, approximately two thirds of the dimorphic genes identified in fetal liver were also dimorphic in adult liver. Sex-biased expression differences unique to adult liver were enriched for growth hormone-responsiveness. Sexually dimorphic gene expression in pre-implantation development is driven by sex-chromosome based transcription, whilst later development is characterised by sex dimorphic autosomal transcription.

Conclusion

This systematic study identifies three distinct phases of sex dimorphism throughout mouse development, and has significant implications for understanding the developmental origins of sex-specific phenotypes and disease in mammals.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1506-4) contains supplementary material, which is available to authorized users.

Lessons from hepatocyte-specific Cyp51 knockout mice: impaired cholesterol synthesis leads to oval cell-driven liver injury

We demonstrate unequivocally that defective cholesterol synthesis is an independent determinant of liver inflammation and fibrosis. We prepared a mouse hepatocyte-specific knockout (LKO) of lanosterol 14α-demethylase (CYP51) from the part of cholesterol synthesis that is already committed to cholesterol. LKO mice developed hepatomegaly with oval cell proliferation, fibrosis and inflammation, but without steatosis. The key trigger was reduced cholesterol esters that provoked cell cycle arrest, senescence-associated secretory phenotype and ultimately the oval cell response, while elevated CYP51 substrates promoted the integrated stress response. In spite of the oval cell-driven fibrosis being histologically similar in both sexes, data indicates a female-biased down-regulation of primary metabolism pathways and a stronger immune response in males. Liver injury was ameliorated by dietary fats predominantly in females, whereas dietary cholesterol rectified fibrosis in both sexes. Our data place defective cholesterol synthesis as a focus of sex-dependent liver pathologies.

Sex-related differences in murine hepatic transcriptional and proteomic responses to TCDD

2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) is an environmental contaminant that produces myriad toxicities in most mammals. In rodents alone, there is a huge divergence in the toxicological response across species, as well as among different strains within a species. But there are also significant differences between males and females animals of a single strain. These differences are inconsistent across model systems: the severity of toxicity is greater in female rats than males, while male mice and guinea pigs are more sensitive than females. Because the specific events that underlie this difference remain unclear, we characterized the hepatic transcriptional response of adult male and female C57BL/6 mice to 500μg/kg TCDD at multiple time-points. The transcriptional profile diverged significantly between the sexes. Female mice demonstrated a large number of altered transcripts as early as 6h following treatment, suggesting a large primary response. Conversely, male animals showed the greatest TCDD-mediated response 144h following exposure, potentially implicating significant secondary responses. Nr1i3 was statistically significantly induced at all time-points in the sensitive male animals. This mRNA encodes the constitutive androstane receptor (CAR), a transcription factor involved in the regulation of xenobiotic metabolism, lipid metabolism, cell cycle and apoptosis. Surprisingly though, changes at the protein level (aside from the positive control, CYP1A1) were modest, with only FMO3 showing clear induction, and no genes with sex-differences. Thus, while male and female mice show transcriptional differences in their response to TCDD, their association with TCDD-induced toxicities remains unclear.

Cytochrome P450 1B1: An unexpected modulator of liver fatty acid homeostasis

Cytochrome P450 1b1 (Cyp1b1) expression is absent in mouse hepatocytes, but present in liver endothelia and activated stellate cells. Increased expression during adipogenesis suggests a role of Cyp1b1 metabolism in fatty acid homeostasis. Wild-type C57BL/6j (WT) and Cyp1b1-null (Cyp1b1-ko) mice were provided low or high fat diets (LFD and HFD, respectively). Cyp1b1-deletion suppressed HFD-induced obesity, improved glucose tolerance and prevented liver steatosis. Suppression of lipid droplets in sinusoidal hepatocytes, concomitant with enhanced glycogen granules, was a consistent feature of Cyp1b1-ko mice. Cyp1b1 deletion altered the in vivo expression of 560 liver genes, including suppression of PPARγ, stearoyl CoA desaturase 1 (Scd1) and many genes stimulated by PPARα, each consistent with this switch in energy storage mechanism. Ligand activation of PPARα in Cyp1b1-ko mice by WY-14643 was, nevertheless, effective. Seventeen gene changes in Cyp1b1-ko mice correspond to mouse transgenic expression that attenuated diet-induced diabetes. The absence of Cyp1b1 in mouse hepatocytes indicates participation in energy homeostasis through extra-hepatocyte signaling. Extensive sexual dimorphism in hepatic gene expression suggests a developmental impact of estrogen metabolism by Cyp1b1. Suppression of Scd1 and increased leptin turnover support enhanced leptin participation from the hypothalamus. Cyp1b1-mediated effects on vascular cells may underlie these changes.