Abstract
Hepatocellular carcinoma (HCC), the most common liver cancer, exhibits a higher incidence in males. Here, we report that mice lacking the bile acid regulators, Farnesoid X Receptor (FXR) and Small Heterodimer Partner (SHP), recapitulate the sex difference in liver cancer risk. Since few therapeutic options are available, we focused on understanding the intrinsic protection afforded to female livers. Transcriptomic analysis in control and FXR and SHP double knockout livers identified female-specific changes in metabolism, including amino acids, lipids and steroids. We examined if the obtained transcriptomic signatures correlate with the survival outcomes for HCC patients to assess the translational potential of this murine HCC model. Gene signature that is unique to the knockout females correspond with low-grade tumors and better survival. Ovariectomy blunts the metabolic changes in female livers and promotes tumorigenesis that, intriguingly, coincides with increases in serum bile acid (BA) levels. Despite similar genetics, we found higher serum BA concentrations in males, whereas female knockout mice excreted more BAs. Decreasing enterohepatic BA recirculation using cholestyramine, an FDA-approved resin, dramatically reduced the liver cancer burden in male mice. Overall, we reveal that sex-specific BA metabolism leading to lower circulating BA concentration protects female livers from developing cancer. Thus, targeting BA excretion may be a promising therapeutic strategy against HCC.
Significance
We show that female-specific gene profiles identified in Fxr-/-, Shp-/- double knockout (DKO) mice correlate with better outcomes for HCC patients and uncover sex difference in circulation and excretion of bile acid. Overall, we demonstrate that increasing (ovariectomy or chemical injury) or decreasing (pharmacologically with FDA approved resin) serum bile acids, not hepatic bile acids, promoted or alleviated liver cancer burden.
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).
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].
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.
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].
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).
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).
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.
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.)
References
- 1.Epidemiology of viral hepatitis and hepatocellular carcinomaGastroenterology 142:1264–1273
- 2.Incidence of Hepatocellular Carcinoma in All 50 United States, From 2000 Through 2012Gastroenterology 152:812–820
- 3.Global, regional, and national age-sex specific all-cause and cause-specific mortality for 240 causes of death, 1990-2013: a systematic analysis for the Global Burden of Disease Study 2013Lancet 385:117–71
- 4.Hepatocellular carcinomaLancet 379:1245–55
- 5.Hepatocellular carcinomaNat Rev Dis Primers 7
- 6.Age-specific sex difference in the incidence of hepatocellular carcinoma in the United StatesOncotarget 8:68131–68137
- 7.Role of reproductive factors in hepatocellular carcinoma: Impact on hepatitis B- and C-related riskHepatology 38:1393–400
- 8.Androgen receptor is a new potential therapeutic target for the treatment of hepatocellular carcinomaGastroenterology 135:947–55
- 9.Reproductive factors, exogenous hormone use and risk of hepatocellular carcinoma among US women: results from the Liver Cancer Pooling ProjectBr J Cancer 112:1266–72
- 10.Foxa1 and Foxa2 are essential for sexual dimorphism in liver cancerCell 148:72–83
- 11.Gender disparity in liver cancer due to sex differences in MyD88-dependent IL-6 productionScience 317:121–4
- 12.The Emerging Hallmarks of Cancer MetabolismCell Metab 23:27–47
- 13.Bile acids activate YAP to promote liver carcinogenesisCell Rep 5:1060–9
- 14.Spontaneous hepatocarcinogenesis in farnesoid X receptor-null miceCarcinogenesis 28:940–6
- 15.Orphan receptor small heterodimer partner suppresses tumorigenesis by modulating cyclin D1 expression and cellular proliferationHepatology 48:289–98
- 16.Spontaneous development of liver tumors in the absence of the bile acid receptor farnesoid X receptorCancer Res 67:863–7
- 17.Mutations in the nuclear bile acid receptor FXR cause progressive familial intrahepatic cholestasisNat Commun 7
- 18.Functional variants of the central bile acid sensor FXR identified in intrahepatic cholestasis of pregnancyGastroenterology 133:507–16
- 19.Mutations in the small heterodimer partner gene are associated with mild obesity in Japanese subjectsProc Natl Acad Sci U S A 98:575–80
- 20.Farnesoid X receptor deficiency induces nonalcoholic steatohepatitis in low-density lipoprotein receptor-knockout mice fed a high-fat dietJ Pharmacol Exp Ther 328:116–22
- 21.Increased activation of the Wnt/beta-catenin pathway in spontaneous hepatocellular carcinoma observed in farnesoid X receptor knockout miceJ Pharmacol Exp Ther 338:12–21
- 22.Systems-level analysis of gene expression data revealed NR0B2/SHP as potential tumor suppressor in human liver cancerMol Cells 30:485–91
- 23.Epigenetic inhibition of nuclear receptor small heterodimer partner is associated with and regulates hepatocellular carcinoma growthGastroenterology 134:793–802
- 24.Pathogenesis of primary sclerosing cholangitis and advances in diagnosis and managementGastroenterology 145:521–36
- 25.Severe bile salt export pump deficiency: 82 different ABCB11 mutations in 109 familiesGastroenterology 134:1203–14
- 26.Genotype correlates with the natural history of severe bile salt export pump deficiencyJ Hepatol 73:84–93
- 27.Bile acids: natural ligands for an orphan nuclear receptorScience 284:1365–8
- 28.The enzymes, regulation, and genetics of bile acid synthesisAnnu Rev Biochem 72:137–74
- 29.Targeting bile-acid signalling for metabolic diseasesNat Rev Drug Discov 7:678–93
- 30.Endogenous bile acids are ligands for the nuclear receptor FXR/BARMol Cell 3:543–53
- 31.The role of bile salt export pump mutations in progressive familial intrahepatic cholestasis type IIJ Clin Invest 110:965–72
- 32.Combined deletion of Fxr and Shp in mice induces Cyp17a1 and results in juvenile onset cholestasisJ Clin Invest 121:86–95
- 33.Bile acids as endogenous etiologic agents in gastrointestinal cancerWorld J Gastroenterol 15:3329–40
- 34.Bile acids promote diethylnitrosamine-induced hepatocellular carcinoma via increased inflammatory signalingAm J Physiol Gastrointest Liver Physiol 311:G91–G104
- 35.Dysregulated hepatic bile acids collaboratively promote liver carcinogenesisInt J Cancer 139:1764–75
- 36.ClueGO: a Cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networksBioinformatics 25:1091–3
- 37.Involvement of p21(WAF1/Cip1), p27(Kip1), and p18(INK4c) in troglitazone-induced cell-cycle arrest in human hepatoma cell linesHepatology 33:1087–97
- 38.Androgen receptor roles in hepatocellular carcinoma, fatty liver, cirrhosis and hepatitisEndocr Relat Cancer 21:R165–82
- 39.Peroxisome proliferator-activated receptor gamma inhibition prevents adhesion to the extracellular matrix and induces anoikis in hepatocellular carcinoma cellsCancer Res 65:2251–9
- 40.Nrf2 is a potential prognostic marker and promotes proliferation and invasion in human hepatocellular carcinomaBMC Cancer 15
- 41.FOXO1 inhibits the invasion and metastasis of hepatocellular carcinoma by reversing ZEB2-induced epithelial-mesenchymal transitionOncotarget 8:1703–1713
- 42.An 18-year-old woman with a 15-cm liver mass and an ammonia level of 342Gastroenterology 143
- 43.Biopsy-proven Hepatocellular Carcinoma in a 53-year-old Woman With Arginase DeficiencyPediatr Dev Pathol 20:517–521
- 44.Diversion of aspartate in ASS1-deficient tumours fosters de novo pyrimidine synthesisNature 527:379–383
- 45.Hepatocellular carcinoma in a research subject with ornithine transcarbamylase deficiencyMol Genet Metab 105:263–5
- 46.Recombinant human arginase inhibits proliferation of human hepatocellular carcinoma by inducing cell cycle arrestCancer Lett 277:91–100
- 47.Discovery of a Carbamoyl Phosphate Synthetase 1-Deficient HCC Subtype With Therapeutic Potential Through Integrative Genomic and Experimental AnalysisHepatology 74:3249–3268
- 48.An Essential Role for Liver ERalpha in Coupling Hepatic Metabolism to the Reproductive CycleCell Rep 15:360–71
- 49.Ethinylestradiol treatment induces multiple canalicular membrane transport alterations in rat liverJ Clin Invest 91:2714–20
- 50.Drug- and estrogen-induced cholestasis through inhibition of the hepatocellular bile salt export pump (Bsep) of rat liverGastroenterology 118:422–30
- 51.Estrogen receptor alpha mediates 17alpha-ethynylestradiol causing hepatotoxicityJ Biol Chem 281:16625–31
- 52.Estrogen Replacement Reduces Risk and Increases Survival Times of Women With Hepatocellular CarcinomaClin Gastroenterol Hepatol 15:1791–1799
- 53.Hormone replacement therapy is associated with reduced hepatocellular carcinoma risk and improved survival in postmenopausal women with hepatitis B: A nationwide long-term population-based cohort studyPLoS One 17
- 54.Multihormonal regulation of hepatic sinusoidal Ntcp gene expressionAm J Physiol Gastrointest Liver Physiol 287:G782–94
- 55.Effects of sulfation patterns on intestinal transport of bile salt sulfate estersAm J Physiol 238:G34–9
- 56.Mechanisms of gender-specific regulation of mouse sulfotransferases (Sults)Xenobiotica 41:187–97
- 57.Quantitative estimation of the hydrophilic-hydrophobic balance of mixed bile salt solutionsJ Lipid Res 30:719–30
- 58.Sulfation and renal excretion of bile salts in patients with cirrhosis of the liverGastroenterology 68:534–44
- 59.Promotive effects of steroids and bile acids on hepatocarcinogenesis initiated by diethylnitrosamineCancer Res 42:2426–8
- 60.Farnesoid X Receptor Activation Promotes Hepatic Amino Acid Catabolism and Ammonium Clearance in MiceGastroenterology 152:1462–1476
- 61.Alternate pathways of bile acid synthesis in the cholesterol 7alpha-hydroxylase knockout mouse are not upregulated by either cholesterol or cholestyramine feedingJ Lipid Res 42:1594–603
- 62.Global gene profiling of spontaneous hepatocellular carcinoma in B6C3F1 mice: similarities in the molecular landscape with human liver cancerToxicol Pathol 39:678–99
- 63.Toxicogenomics directory of chemically exposed human hepatocytesArch Toxicol 88:2261–87
- 64.Identification of prognostic biomarkers for patients with hepatocellular carcinoma after hepatectomyOncol Rep 41:1586–1602
- 65.Identification of a bile acid and bile salt metabolism-related lncRNA signature for predicting prognosis and treatment response in hepatocellular carcinomaSci Rep 13
- 66.Integrated analysis of microbiome and host transcriptome reveals correlations between gut microbiota and clinical outcomes in HBV-related hepatocellular carcinomaGenome Med 12
- 67.Association between Pre-Diagnostic Serum Bile Acids and Hepatocellular Carcinoma: The Singapore Chinese Health StudyCancers (Basel 13
- 68.Prevention of spontaneous hepatocarcinogenesis in farnesoid X receptor-null mice by intestinal-specific farnesoid X receptor reactivationHepatology 61:161–70
- 69.Classification and prediction of survival in hepatocellular carcinoma by gene expression profilingHepatology 40:667–76
- 70.Application of comparative functional genomics to identify best-fit mouse models to study human cancerNat Genet 36:1306–11
- 71.A novel prognostic subtype of human hepatocellular carcinoma derived from hepatic progenitor cellsNat Med 12:410–6
- 72.Sixty-five gene-based risk score classifier predicts overall survival in hepatocellular carcinomaHepatology 55:1443–52
- 73.Yes-associated protein 1 and transcriptional coactivator with PDZ-binding motif activate the mammalian target of rapamycin complex 1 pathway by regulating amino acid transporters in hepatocellular carcinomaHepatology 63:159–72
- 74.A unique metastasis gene signature enables prediction of tumor relapse in early-stage hepatocellular carcinoma patientsCancer Res 70:10202–12
- 75.Genome-wide survey of recurrent HBV integration in hepatocellular carcinomaNat Genet 44:765–9
- 76.Neoangiogenesis-related genes are hallmarks of fast-growing hepatocellular carcinomas and worst survival. Results from a prospective studyGut 65:861–9
- 77.Electronic address, w.b.e. and N. Cancer Genome Atlas Research, Comprehensive and Integrative Genomic Characterization of Hepatocellular CarcinomaCell 169:1327–1341
- 78.Distinct properties of cell-type-specific and shared transcription factor binding sitesMol Cell 52:25–36
- 79.Research resource: Aorta- and liver-specific ERalpha-binding patterns and gene regulation by estrogenMol Endocrinol 28:1337–51
- 80.Changes in Gene Expression and Estrogen Receptor Cistrome in Mouse Liver Upon Acute E2 TreatmentMol Endocrinol 30:709–32
Article and author information
Author information
Version history
- Sent for peer review:
- Preprint posted:
- Reviewed Preprint version 1:
Copyright
© 2024, Patton et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
- views
- 268
- downloads
- 12
- citations
- 0
Views, downloads and citations are aggregated across all versions of this paper published by eLife.