Introduction

Liver cancer, a leading cause for cancer-related death, has diverse etiologies, and displays sex-difference with reduced risk in females compared to males [1–5]. Since current therapies for liver cancer fall short, we posit that understanding molecular mechanisms functioning in the female livers will reveal new therapeutic targets. Earlier studies have reported the role of sex hormones [6–9], transcription factors FoxA1/A2 [10], and cytokine Il6 signaling [11] in regulating the sex difference in hepatocellular carcinoma (HCC), but the role of metabolic pathways remains poorly understood.

Rewiring of cellular metabolism enables the tumor cells to maintain viability and grow disproportionately [12]. We previously showed that the combined deletion of nuclear receptors, Farnesoid X Receptor (FXR) and Small Heterodimer Partner (SHP) resulted in spontaneous liver cancer in the year-old male mice [13]. In this study we report that unlike the males, female Fxr^-/-, Shp^-/-double knockout (DKO) mice exhibit protection against tumorigenesis, and thus mimic the sexual dimorphism in liver cancer incidence observed in clinics. Although 15-month-old individual Fxr knockout and individual Shp knockout mice were previously shown to develop liver cancer, but unlike the DKO mice, their incidence does not show 100 percent penetrance nor sex differences [14–16].

Mutations and reduction in Fxr, and Shp transcript levels have been noted in cholestasis (reduced bile flow and subsequent increase in hepatic and serum bile acids (BA)), fatty liver disease and liver cancer [17–23]. Moreover, individuals with chronic cholestasis exhibit increased risk for HCC [24–26]. Typically, BA levels are tightly controlled via receptor signaling, including FXR and SHP [27–31].

Consistently, combined loss of Fxr and Shp in mice results in juvenile onset of cholestasis that progresses to HCC [32]. We and others have shown that excessive accumulation and dysregulation of BA homeostasis is directly linked with liver cancer risk [13, 26, 33–35]. However, whether BAs are contributing factors to the sex differences seen in HCC prevalence has not been evaluated.

Therefore, we performed transcriptomic analysis to identify distinct gene profiles from both sexes of control and DKO mice. Then, using five separate human clinical HCC cohorts, we tested the clinical utility of the identified gene signatures from our mouse model. Next, we investigated the role for the endogenous estrogen signaling in the DKO mice by performing ovariectomy. We measured hepatic, serum, urine, and fecal BAs from male and female mice to understand their homeostasis.

Finally, we manipulated the circulating BA levels in the DKO mice either with chemical challenge or BA binding resins and examined its consequence on hepatic tumorigenesis. Overall, our data uncover that the differential BA homeostasis between the two sexes can orchestrate the observed gender differences in HCC burden in clinics.

Results

Fxr^-/-, Shp^-/- double knockout (DKO) mice phenocopy clinical features of HCC

Here, we report that DKO mice exhibit the sexually dimorphic incidence of HCC observed in the clinic. Despite the loss of BA homeostatic machinery, one-year-old DKO female mice did not develop liver tumors but showed modest fat accumulation and mild fibrosis. On the contrary, DKO male livers revealed HCC and well-defined adenomas, along with robust steatosis, and fibrosis (Fig. 1A-F). Even at six months of age, female livers were smaller and displayed reduced hepatic fibrosis compared to males (Supplementary Fig. 1). Sex-difference in tumor burden was reflected in the gross liver to body weight ratio which is significantly higher in DKO males than the DKO females (Fig. 1G). Serum ALT and AST were elevated in the DKO animals compared to WT (Fig. 1H), which correlates with the cholestatic phenotype of the mice. More importantly, liver cancer patients show reduction in Fxr and Shp transcript expression (Fig. 1I).

Figure 1:

Sex-specific metabolic programs regulate liver tumorigenesis

To identify transcriptional mechanisms that can contribute towards the sex differences in the incidence of hepatic tumorigenesis, we analyzed one-year-old male and female livers. DKO males and females displayed striking differences in hepatic gene expression profile (Supplementary Table 1, Supplementary Fig. 2 and GEO GSE151524) with DKO males showing enrichment of endoplasmic reticulum stress, unfolded protein response, and immune function (Fig. 2A-B). Additionally, network analysis with ClueGO [36] revealed interactions between drug metabolism, inflammation, ERK signaling, and steroid metabolism in DKO males (Fig. 2C). On the contrary, DKO females displayed pathway enrichment of steroid metabolism and clustering of lipid, glucose, amino acid, and steroid metabolism, along with increased sulfotransferase activity (Fig. 2D). Next, we parsed the sex-specific upregulated gene sets to identify unique transcription factor motifs. Overrepresented motifs in DKO males, including AR, FOXA1, FOXA2, NRF2 and PPARγ [10, 37–40] correlate with tumor-promoting functions (Supplementary Table 2). In contrast, in DKO females, FOXO1, E2F, and ERα (Supplementary Table 3) were dominant motifs and are associated with regulating metabolic function during liver carcinogenesis [41].

Figure 2:

Transcriptomic signature of the DKO mice correlates with poor overall survival in the clinical datasets

To investigate the clinical relevance of the DKO mouse model, we analyzed the WT and DKO murine transcriptomic signatures in a sex specific manner and compared these to five separate clinical HCC datasets. Using class prediction, the patient data were sorted based on their similarity to one-year-old DKO gene signatures (Supplementary Fig. 2-3). Computational prediction scores (BCCP: 1 represents complete match and 0 represents no match) using the patient samples revealed that the DKO male signature matches the later stages (2-5) of liver cancer, whereas the DKO female-specific signature matches with stage 0 (Supplementary Fig. 3).

Of the 1100 patient data, we found approximately (∼45%) showed similar transcriptomic signature to that of either DKO male or DKO female which corresponded to lower overall survival (OS), but not recurrence-free survival (RFS). WT gene signatures were used as controls (Fig. 3A-B). Although DKO female mice do not develop liver cancer, it is pertinent to note that these mice do lack Fxr and Shp expression, display chronic cholestasis similar to their male counterparts, and hence the global gene changes associate with poor OS. On the contrary, when we focused on the gene signature that was distinctly changed only in the DKO female livers not the DKO males, we found that patients (∼54.71%) who displayed this subset of gene signature had better OS as well as RFS (Fig. 3C). These findings underscore a high potential for clinical translation of data generated from the DKO mouse model. Moreover, by focusing on specific transcript changes in the DKO female livers, we uncovered a subset of metabolic genes that correspond to better survival and might be responsible for their protection against cancer.

Figure 3:

We first focused on the pathways pertaining to amino acid metabolism and ureagenesis since individuals with mutations in the urea cycle disorder exhibit high risk for liver cancer incidence [42–45]. Moreover, amino acids are known to play a central role in tumorigenesis. We found overexpression of genes regulating amino acid metabolism and ureagenesis, including carbamoyl phosphate synthetase (Cps1), ornithine transcarbamylase (Otc), argininosuccinate synthetase (Ass1), argininosuccinate lyase (Asl), and arginase (Arg1) in DKO female livers (Supplementary Fig. 4A). The upregulation of genes which encode these enzymes correlated well with the protection afforded to the DKO female livers since loss of function mutations of these genes are linked to HCC [45–47]. Additionally, our analysis showed that patients with increased expression of urea cycle genes (DKO-UreaCycle) exhibited a better clinical outcome (Supplementary Fig. 4B).

Estrogen signaling controls amino acid and bile acid metabolism in the liver

Since estrogen signaling was previously shown to regulate amino acid metabolism [48], we examined its role in controlling the expression of urea cycle genes in the DKO female livers. To do this, we ovariectomized (OVX) DKO mice and found that indeed the hepatic expression of all these genes, Cps1, Asl1, Ass, Otc and Arg1 were significantly blunted in the absence of endogenous estrogen signal (Supplementary Fig. 4C). But when we measured the urea cycle metabolites, we did not find any significant change in the intermediate nor urea production except for a decrease in ornithine levels (Supplementary Fig. 4D), in DKO females compared to the DKO males. We reason that static measurements may not reflect the flux into the urea cycle.

Besides amino acid metabolism, estrogen signaling has been shown to affect BA homeostasis and cause cholestasis [49–51]. So, we anticipated that OVX will conversely lead to lower BA levels in DKO female mice. Instead, we found that OVX led to liver cancer development in otherwise resistant year old DKO female mice (Fig. 4A-B) and the serum BA levels doubled in these mice (Fig. 4C). This finding is consistent with the tumorigenic role for BAs, and coincides with clinical data wherein post-menopausal women are susceptible to developing HCC whereas their resistance to liver cancer can be restored upon hormone replacement therapy [52, 53].

Figure 4:

To overcome the confounding effects of ageing and cancer, we examined young WT and DKO female mice with and without OVX. In addition, we also challenged these mice with BA excess (Fig. 4D). As expected, OVX resulted in reduction of basal hepatic Erα gene expression in both WT and DKO mice (Fig. 4E). In the DKO mice, which display high basal levels of BA synthesis and sulphation genes, we found dramatic induction of Erα gene upon BA treatment (Fig. 4E).

Importantly, the increase in Erα gene coincided with reduced expression of Cyp7a1, Cyp8b1 and sulphotransferase known to sulphate estrogen and BAs, Sult2a1 transcripts, indicating a role for estrogen signaling in regulating BA homeostasis in the DKO livers. OVX in WT mice led to lower basal levels of Cyp7a1 and Sult2a1 but not Cyp8b1 whereas all these three genes were significantly reduced in the DKO livers (Fig. 4E-4H). Unlike the OVX WT, which maintained CA-mediated suppression of BA synthetic genes, consistent with intact FXR signaling, DKO+OVX mice did not alter their expression (Fig. 4F-G).

We next examined if the recruitment of ERα to BA synthesis genes exhibit any sex difference in WT and DKO livers by ChIP-PCR. We find that ERα was preferentially recruited to Cyp7a1 in a sex specific manner (Figure 4I). Cyp8b1 showed a similar trend but not Sult2a1. Also, we did not find any sex-specific patterns in ERα occupancy in Ldlr, and Pgr genes, which were used as positive controls for ERα ChIP assays (Figure 4I). These data along with increased BAs upon OVX suggest that ERα signaling is pertinent to control BA synthesis especially in the absence of FXR, as seen in the CA-fed sham DKO mice.

DKO mice display sexual dimorphism in BA homeostasis

WT mice do not show overt changes in serum BAs between the two sexes, however, genetically identical DKO mice displayed dramatically lower serum BAs in females compared to males (Supplementary Table 4). Nonetheless, the serum BA concentration in DKO females was higher than WT females. This was intriguing.

So, we analyzed the expression of genes involved in BA synthesis, transport, and metabolism in both sexes of DKO mice. Consistent with Fxr and Shp deletion that results in the loss of negative feedback on BA biosynthesis, both sexes of DKO mice have significantly higher expression of Cyp7a1 and Cyp8b1 genes that are involved in classical BA synthesis (Fig. 5A). The male dominant expression of Cyp7b1 in the WT is lost in the DKO mice. On the other hand, Cyp27a1, which initiates alternative BA synthesis, was increased in a female-specific manner (Fig. 5A).

Figure 5:

Next, we examined BA transport. We found that hepatic transcript levels of the key BA efflux pump, bile salt export pump (Bsep, Abcb11) was reduced in both sexes of DKO mice, consistent with loss of Fxr (24) (Fig. 5B), while the expression of canalicular efflux transporters, Abcb1 (Mdr1) and Abcc2 (Mrp2) were unchanged (Fig. 5B). Also, the BA uptake transporter, sodium taurocholate co-transporting polypeptide (Ntcp, Slc10a1) showed lower transcript levels in females (Fig. 5B), which is in line with previous findings that estradiol can downregulate Slc10a1 expression [54].

We then investigated the transcript expression of Sult2a1, which contributes to BA sulfation—a modification that can reduce enterohepatic recirculation [55]. As expected, hepatic Sult2a1 expression was predominant in females irrespective of the genotype (Fig. 5C) [56]. Sulphated BAs are excreted in urine to eliminate excess BA during cholestasis [57, 58]. Total urine BA levels were higher in DKO males, reflecting a larger circulating BA pool than in DKO females (Supplementary Table 5). However, DKO female mice exhibited higher percent sulphated BAs (Fig. 5D, Supplementary Table 5) that corroborates with the high Sult2a1 expression in females.

BA compositional analysis was performed in the serum, hepatic, urine, and feces of DKO males and females (Fig. 5E-H). Both sexes of DKO mice showed abundant muricholates in the serum but there were modest differences in the composition, indicating a slightly hydrophilic BAs in the DKO females (Fig. 5E). Moreover, hepatic BA composition was indifferent between the two sexes (Fig. 5F). These results indicated that rather than synthesis or transport, excretion may be different between DKO male and females. Notably, we found that both urine and fecal levels and composition between male and female DKO were distinct (Fig. 5G-H).

Since urinary BA excretion is a small proportion and cannot explain the 50% decrease in circulating BAs in DKO females, we performed untargeted metabolomics using the fecal samples. BAs accounted for 20% of the fecal samples in the males whereas in the females it was double the amount indicative of twice the amount being excreted in DKO females. The female DKOs may be protected against detrimental tumor promoting BA signaling due to their higher BA excretion.

Increasing fecal BA excretion is sufficient to reduce liver cancer risk

Finally, to test this, we promoted fecal BA excretion in DKO males by using cholestyramine (CHR), a resin that binds BAs. We fed nine-month-old DKO male mice with a 2% CHR-containing diet since, by this age, they have already developed tumor nodules. The CHR diet was continued until one year of age, mimicking a therapeutic intervention strategy (Fig. 6A). CHR diet as expected led to a 50% reduction in circulating BA levels and altered BA composition in DKO males (Fig. 6B-C). DKO males fed chow exhibit severe hepatic tumorigenesis, whereas CHR-fed DKO males have a drastically lower tumor burden with only small liver nodules and are protected from developing aggressive carcinomas. Histological analysis revealed that CHR treatment lowered the number of nodules and dysplastic changes, but increased steatosis in DKO males (Fig. 6D).

Figure 6:

Unbiased cluster analysis of serum BA composition between the two cohorts of DKO mice revealed that the BA profiles of CHR-fed DKO male mice grouped with DKO females (Supplementary Fig. 5). Conversely, increasing circulating BAs by bile duct injury with 3,5-diethoxycarbonyl −1,4-dihydrocollidine (DDC) in DKO females resulted in the development of large liver tumors in DKO females (Supplemental Fig. 6).

Overall, these findings demonstrate that modulating circulating BAs is sufficient to change the liver cancer outcome with lowering their levels leading to subsequent protection and vice versa.

Discussion

Here, we demonstrate that the sex-differences in BA homeostasis confers the sexual disparity noted in HCC risk. Importantly, elevated BA concentrations are reported in patients with HCC [13, 26, 33–35, 59]. Using a genetic mouse model of excess BAs that develop spontaneous HCC, we uncovered distinct transcriptional control of metabolism between the two sexes. Both Fxr and Shp transcript levels were downregulated in HCC patients. Moreover, differential gene expression specifically of the DKO female correlated well with better survival, highlighting the translational relevance for our model. Thus, these gene signatures could be utilized as a potential prognostic marker for HCC progression and survival.

Both BA homeostasis, and amino acid metabolism were altered between the two sexes. Of note, genes controlling ureagenesis were higher in the DKO females, and consistent with previous findings, we were able to recapitulate estrogen-mediated regulation of some these genes signaling [48]. In line with these findings, patients with urea cycle enzyme deficiencies have a 200x higher incidence of HCC, highlighting the importance of amino acid metabolism in hepatic tumorigenesis [42–45]. Also, BAs have been shown to promote amino acid catabolic machinery [60], which indicates that BAs may be a central node in liver cancer. Intriguingly, hepatic urea analysis did not reveal any difference between the DKO male and female mice. A caveat being we are measuring steady state urea levels rather than the flux of this pathway.

We examined and found estrogen signaling can regulate the expression of BA synthesis and sulphation genes. DKO female mice challenged with CA diet showed robust increase in hepatic Erα transcript coincided with the suppression in BA synthesis in the absence of Fxr and Shp. Consistently higher recruitment of ER to the classical BA synthetic genes was noted in female livers. DKO OVX mice that have blunted Erα gene expression, also showed lower transcript levels of Cyp7a1, and Cyp8b1 that were unaltered post-CA challenge. Currently we do not know how to reconcile this data other than indicating a potential ER-independent mechanism.

Nonetheless, we confirmed the known sex differences in BA synthesis, such as a female-dominant Cyp27a1 expression and male-dominant Cyp7b1 pattern in WT mice. Loss of Fxr and Shp altered many of the gene expressions irrespective of the sex. For instance, Cyp2c70 expression was reduced per se irrespective of the sex and the male dominance of Cyp7b1 was lost in the DKO mice.

Of note, OVX of DKO females increased the serum BA levels and lost their protection against the development of liver tumorigenesis. This finding fully recapitulates the clinical data, wherein post-menopausal women are equally prone to HCC incidence as males.

BA analysis show that DKO female mice have a hydrophilic composition and excrete BA proportions. So, we tested and demonstrate a potent therapeutic utility of reducing BA levels in serum using a generic FDA-approved BA binding resin, Cholestyramine (CHR) in dramatically reducing the tumor burden. This study highlights that lowering enterohepatic recirculation is a beneficial strategy in modulating liver cancer. Although Cyp7a1 expression is reported to be induced in CHR fed mice [61], long-term CHR feeding in DKO mice lowered Cyp7a1 expression but induced Cyp8b1 transcripts. Conversely, DDC fed DKO females that develop hepatic tumors show a corresponding decrease in Cyp8b1 transcript. Also, patients with HCC exhibit a reduction in Cyp8b1 expression [62–64], which promotes a more hydrophilic ratio of BA composition.

Although species differences in BAs between mice and human are a limitation, several fundamental understandings have been gained from mouse experiments. Another caveat is that the DKO mouse model mimics the progression of cholestasis to HCC and not all the etiologies, so the observed sex differences in circulating BAs may be limited to these subsets of HCC. Nevertheless, elevated BA concentrations are seen in various liver disease conditions and particularly inherited disorders of cholestasis predisposes to HCC onset. More recently clinical studies support the utilization of BAs as prognostic markers [65–67].

Overall, our findings demonstrate that female cholestatic mice exhibit increased excretion and lower serum BAs compared to males. But hepatic BAs were not different between the sexes. These data highlight a fundamental difference in circulating BA that contributes to the sex-differences seen in HCC incidence. Accordingly, we show that lowering enterohepatic BA recirculation is beneficial in our model. Our results are in line with previous findings that had implicated intestinal FXR signaling being crucial rather than hepatic FXR to prevent liver tumorigenesis [68]. Taken together, these results reveal that drugs inhibiting intestinal reabsorption of BAs (Asbt inhibitor, IBAT inhibitor) that are on clinical trials for NASH and cholestasis can be evaluated as potential therapeutics to combat HCC.

Methods Experimental

Design

This study was designed to elucidate the role for bile acids (BAs) in the sexually dimorphic incidence of HCC and assess the therapeutic benefits of reducing circulating BA levels on HCC development. Fxr^-/-, Shp^-/- (DKO) mice were maintained on a C57BL/6 background at the University of Illinois, Urbana-Champaign. Male and female mice were sacrificed at 8-12 weeks, or 6 and 12 months after birth. Mice were housed on a standard 12-hour light/dark cycle and fed normal chow and water ad libitum. To study estrogen signaling, bilateral ovariectomies were performed on WT and DKO females at 8-10 weeks old, followed by one week of recovery. 2% Cholestyramine-supplemented chow was fed to 9-month-old DKO male mice for a period of 3 months while 0.1% DDC (3,5-diethoxycarbonyl-1, 4-dihydrocollidine) was fed to 10-month-old DKO female mice for 3 months. Urea cycle studies were performed on mice after overnight fasting. All studies were carried out as outlined in the Guide for the Care and Use of Laboratory Animals prepared by the National Academy of Sciences and published by the National Institutes of Health (National Institutes of Health publication 86-23, revised 1985). For biological harvesting, mice were anaesthetized and euthanized as described by IUCAC. Tissue was flash frozen in liquid nitrogen and blood was collected and spun down for serum.

Serum Chemistry

Blood was collected by retro-orbital bleeding and centrifuged at 8000g x 10 minutes to separate serum. Serum ALT and AST were measured using Infinity ALT and Infinity AST kits (Thermo Scientific). Calorimetric measurement of serum and hepatic bile acids was performed with the Total Bile Acid (NBT method) kit (Genway Biotech).

Bile Acid analysis

Serum and urine from DKO male and female mice fed chow was analyzed for composition of bile acids and their sulfated metabolites at University of Nebraska Medical Center. Liquid chromatographic-mass spectrometry analysis was performed with a Waters ACQUITY column (Milford, MA). Bile acids and internal standards were prepared in methanol and analyzed. Serum from DKO male and female mice fed chow, DKO males fed CHR, and DKO females fed DDC was analyzed for bile acid composition at Baylor College of Medicine Metabolomics Core. Liquid chromatographic-mass spectrometry analysis was performed with a Waters ACQUITY UPLC BEH C18 column (Milford, MA). Bile acids were detected in negative mode, with L-Zeatine added to each sample as an internal standard.

Metabolite Profiling

Liver tissue was weighed and sonicated in 70% methanol, followed by centrifugation. Supernatant was flash frozen and used for subsequent LC-MS analysis for urea cycle metabolites. Tissue lysate was used for BCA assay to determine protein concentration of each sample. All metabolite concentrations were normalized to protein concentration of the lysate.

Untargeted Metabolomics

Fecal samples were weighed into microcentrifuge tubes and homogenized in 50% MeOH/H2O solution with a 1:10 w/v ratio, for 5 minutes at 5 Hz. The samples were centrifuged at 14000 rpm for 15 minutes, then a 200uL aliquot of each supernatant was transferred to a 96-well plate and dried under centrifugal vacuum. The dried extracts were covered and stored at −80 °C until analysis, at which time the samples were resuspended in 200uL of 50% MeOH/H2O solution with 1uM sulfadimethoxine as internal standard and diluted three-fold for analysis. Untargeted LC-MS/MS was performed on a Thermo Vanquish UPLC system coupled to a Q-Exactive Orbitrap mass spectrometer (ThermoFisher Scientific, Bremen, Germany). A polar C18 column (Kinetex polar C18, 100 x 2.1 mm, 2.6 μm particle size, 100 A pore size; Phenomenex, Torrance, CA USA) was used as the stationary phase, and a high-pressure binary gradient pump was used to deliver the mobile phase, which consisted of solvent A (100% H2O + 0.1 % formic acid (FA)) and solvent B (100% acetonitrile (ACN) + 0.1 % FA). The flow rate was set to 0.5mL/min and the injection volume for each sample was 5uL. Following injection, samples were eluted with the following gradient: 0-1.0 min, 5% B; 1.0-1.1 min, 25%; 6.0 min, 70%; 7.0 min, 100%; 7.5-8.0 min, 5%. MS data was collected in positive mode and electrospray ionization (ESI) parameters were set to 53 L/min for sheath gas, 14 L/min for auxiliary gas, 0 L/min for spare gas, and 400°C for auxiliary gas temperature. The spray voltage was set to 3500 V, the capillary temperature to 320 °C, and the S-Lens radio frequency level to 50 V. MS1 data were collected from 150-1500 m/z with a resolution of 35,000 at m/z 200 with one micro scan. The maximum ion injection time was set to 100 ms with an automatic gain control (AGC) target of 1.0E6.

MS/MS spectra were collected using data dependent acquisition (DDA), where the top 5 most abundant ions in the MS1 scan were selected for fragmentation. Normalized collision energies were increased stepwise from 20, 30, to 40. MS2 data were collected with a resolution of 17,500 at m/z 200 with one micro scan and an AGC of 5.0E5. All untargeted LC-MS/MS data are publicly available from the MassIVE data repository under accession number MSV000089715.

MS1 feature detection and MS/MS pairing was performed using MZmine 2.37corr17.7_kai_merge (55, 56). An intensity threshold of 5E4 and 1E3 were set for MS1 and MS2 detection, respectively, with centroid data. MS1 chromatogram construction was performed using the ADAP chromatogram builder, where the minimum group size was set to 5, group intensity threshold was 5E4, minimum highest intensity was 1.5E5, and mass tolerance was 0.005 m/z or 10 ppm. Chromatogram deconvolution was then performed using a local minimum search algorithm with a chromatographic threshold of 80%, a search minimum in retention time (RT) range of 0.2 min, minimum relative height of 1%, minimum absolute threshold height of 1.5E5, minimum ratio for top/edge of 1, and a peak duration of 0.05-2.0 min. Pairing between MS1 and MS2 was performed with a mass tolerance of 0.005 m/z or 10 ppm and RT range of 0.2 min. Isotope peaks were grouped, then features from different samples were aligned using the same mass and RT tolerances; alignment was performed by placing a weight of 75 on m/z and 25 on RT. A peak area feature table was exported as a .csv file and consensus MS/MS spectral data were exported in mgf format. Feature-based molecular networking and MolNetEnhancer workflows were then performed with this data using GNPS (gnps.ucsd.edu). The corresponding jobs can be found at: https://gnps.ucsd.edu/ProteoSAFe/status.jsp?task=d697d44ec18440d29d0771f84ba7cccd andhttps://gnps.ucsd.edu/ProteoSAFe/status.jsp?task=e3edba56efba4a27b073a9031c60b5e5, respectively.

Histology

Liver samples were fixed in 10% neutral-buffered formalin, sections were cut at a thickness of ∼5 microns and used for hematoxylin and eosin, and sirius red staining.

RNA extraction and quantitative PCR analysis

Total RNA from liver was prepared according to the TRIzol (Invitrogen) protocol. cDNA was synthesized using Maxima Reverse Transcriptase (Thermo Scientific), according to manufacturer’s protocol. q-RTPCR was performed on an Illumina Eco Platform. For qRT-PCR analysis, 50 ng of cDNA was added to each SYBR green-based reaction. qRT-PCR primers are provided in Supplementary Table 6.

Microarray

Microarray was performed by Dr. Ju-Seog Lee’s laboratory at the MD Anderson Cancer Center. Liver samples from 12-month-old male and female WT and DKO mice were collected and snap-frozen. Total RNA was isolated, labeled, and hybridized to BeadChip Array MouseWG-6 (Illumina). Bead chips were scanned with an Illumina BeadArray Reader. Microarray analysis was performed on the Illumina mouseRefseq-8 Expression platform. Upregulated gene sets were generated from genes with fold change > 1.3 (p<0.0001) compared to control group (i.e. DKO males vs. DKO females). These gene sets were then used for downstream analyses with DAVID Bioinformatics Resources Analysis Software and ClueGO [36].

Transcription Factor Motif Analysis

GeneXplain software was used to identify enriched transcription factor binding sites (TFBS) using the upregulated gene sets generated from the microarray. Analysis included regions from −1000 to 100 bp relative to transcription start site. Only those TFBS enriched with p≤ 0.01 were included in tables.

Extraction of transcriptomic signature

Multiple transcriptomic signatures were extracted from the microarray data of the DKO mouse model (Table 1). DKO_All signature was generated from the comparison between wild type (WT) male and female mouse, and DKO_Male and DKO_Female signature from WT male and female mouse, respectively.

DKO_FvsM, DKO_Estrogen, DKO_BA, and DKO_Urea signatures were made from the comparison between DKO male and female mouse. Using the gene expression dataset after normalization, signature genes were selected by T-test and logFC (p<0.001 and log2FC>1 or <-1).

Gene expression data from HCC tumors

Gene expression data from the National Cancer Institute (NCI) cohort were generated in earlier studies [69–71] and the data are publicly available from the NCBI’s GEO database (GSE1898 and GSE4024). Gene expression data from Korea, Samsung, Modena, and Fudan cohorts have been described previously and are available from the NCBI’s GEO database (accession numbers, GSE14520, GSE16757, GSE43619, GSE36376, and GSE54236) [72–76]. TCGA RNA sequencing data for HCC was downloaded from the University of California, Santa Cruz, Genomics Institute (https://xenabrowser.net/) 12. FPKM-normalized data were log-transformed.

Tumor specimens and clinical data were obtained from HCC patients who had undergone hepatectomy as a primary treatment for HCC at multiple institutes as described in their original study. Except for TCGA cohort, patients and tissues were collected based on availability of high quality of frozen tissues for genomic studies. For TCGA cohort [77], surgical resection of biopsy biospecimens were collected from patients diagnosed with HCC and had not received prior treatment for their disease (ablation, chemotherapy, or radiotherapy). Institutional review boards at each tissue source site reviewed protocols and consent documentation and approved submission of cases to TCGA. Hematoxylin and eosin (H&E) stained sections from each sample were subjected to independent pathology review to confirm that the tumor specimen was histologically consistent with the allowable HCC. Each case was reviewed independently by at least 3 liver pathologists, with no clinical or molecular information.

Prediction model

To predict the class similar with DKO signature in human HCC cohort, we used a classification algorithm based on Bayesian compound covariate predictor (BCCP). After the integration of the signature matrix with human HCC dataset, the Bayesian probability for each human HCC sample were calculated by using class prediction procedure in BRB Arraytools. We set 0.5 as the cut-off of Bayesian probability for each signature.

Validation of clinical relevance

Prognostic significance was evaluated rigorously for overall and recurrence-free survival in human HCC cohort based on the predicted class calculated by the BCCP algorithm using multiple DKO-signatures. A total of 5 human HCC transcriptomic cohorts were used in this study (Fudan, Korea, Samsung, TCGA, Modena). All DKO-signatures were evaluated in each human HCC cohort and meta-cohort. To identify the gender difference in human HCC cohort, we did subgroup analysis for gender, age in meta-cohort. BCCP scores (BCCP probability) were compared in all population, and gender subgroup. The analysis for potential correlation between the class predicted by DKO-signature and staging system of HCC in terms of TNM, BCLC, CLIP classification were performed in meta-cohort.

ERα ChIP-Seq Analysis

ERα-ChIP assay was performed in both sexes of WT and DKO mice. ERα-F10 antibody (sc-8002, Santa Cruz) was used to perform the pulldown followed by qPCR. We also analyzed for BED files for ERα ChIP-Seq from three independent studies [78–80] were obtained from Cistrome DB and visualized using UCSC genome browser on mouse GRCm38/mm10 assembly. For tracks 3 and 4, track size was set to auto-size.

Statistical Analysis

All statistical tests were performed using GraphPad Prism software. Data are presented as means ± SEM. Multiple group comparisons were analyzed by one-way and two-way ANOVA with the post hoc Bonferroni test. Unpaired t-test was used for comparison between two groups. P values ≤0.05 were determined to be significant, unless otherwise noted in legends.

Data Availability

The gene expression data generated and used in this publication have been deposited in NCBI’s Gene Expression Omnibus and are accessible through GEO Series accession number GSE151524. (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE151524).

Acknowledgements

We thank the Systems biology laboratory at the University of Texas, MD Anderson Cancer Center for initial analysis and performing the microarray studies. We would also like to thank the metabolite analysis core at Baylor College of Medicine for performing bile acid composition analysis. We thank Dr. Bhoomika Mathur for initial help with the DDC experiment. The authors also thank Drs. Auinash Kalsotra and Stephanie Ceman for their comments and critiques during the preparation of this manuscript. We would like to thank Ms. Angela Major at Texas Children’s’ Hospital for histological preparation and analysis. The authors thank Dr. Lucas Li at the Roy Carver Metabolomics Core at University of Illinois at Urbana-Champaign for Urea Cycle metabolite analysis.

Author contributions

M.E.P., S.K., A.E.D and S.A. contributed to conception, experimental design, data acquisition, analysis and interpretation, and drafting the article. L.J.T. helped with urea cycle expression analysis, R.N.T. and Y.A. analyzed bile acid composition and interpretation from serum and urine, E.G., M.P, P.D performed untargeted metabolomics, A.E.D, Q.Z, Z.M.E analyzed ERa ChIP-seq data, H.S.L and J.S.L performed all the HCC patient cohort analysis, M.J.F. evaluated and interpreted the liver histology.

Competing interests

None

Disclosures

PCD is an advisor and holds equity in Cybele and Sirenas and a Scientific co-founder, advisor and holds equity to Ometa, Enveda, and Arome with prior approval by UC-San Diego. PCD also consulted for DSM animal health in 2023.

Funding

This work was supported in part by NIDDK grant, R01 DK113080 (S. A), Research Scholar Grant ACS 132640-RSG (S.A) and UIUC start-up funds (S.A.)