Next Article in Journal
The Organoid Era Permits the Development of New Applications to Study Glioblastoma
Previous Article in Journal
Immune Checkpoint Blockade in Gynecologic Cancers: State of Affairs
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Chronic Disruption of the Late Cholesterol Synthesis Leads to Female-Prevalent Liver Cancer

1
Centre for Functional Genomics and Bio-Chips, Institute of Biochemistry, Faculty of Medicine, University of Ljubljana, SI-1000 Ljubljana, Slovenia
2
Rudol-Schönheimer-Institute of Biochemistry, Divison of General Biochemistry, Faculty of Medicine, University of Leipzig, 04103 Leipzig, Germany
3
Medical Experimental Centre, Institute of Pathology, Faculty of Medicine, University of Ljubljana, SI-1000 Ljubljana, Slovenia
4
Institute of Pathology, Faculty of Medicine, University of Ljubljana, SI-1000 Ljubljana, Slovenia
*
Author to whom correspondence should be addressed.
Current Address: Section of Hepatology, Clinic and Polyclinic for Gastroenterology, Hepatology, Infectiology, Pneumology, University Clinic Leipzig, 04103 Leipzig, Germany.
Cancers 2020, 12(11), 3302; https://doi.org/10.3390/cancers12113302
Submission received: 14 September 2020 / Revised: 29 October 2020 / Accepted: 5 November 2020 / Published: 9 November 2020
(This article belongs to the Section Molecular Cancer Biology)

Abstract

:

Simple Summary

Hepatocellular carcinoma is a disease with a variety of molecular triggers and is usually reported to prevail in males. However, after the menopause, the disease is also increasing in the female population. Herein, we discovered that chronic depletion of cholesterol synthesis due to the knock-out of the gene Cyp51 from this pathway leads to female prevalent hepatocarcinogenesis in aging mice. There is a high similarity between our mouse model and the situation in humans. Multiple deregulated pathways of hepatocarcinogenesis are shared. A female-dependent metabolic reprogramming leading to this type of liver cancer is exposed for the first time and reflects on deregulated cholesterol synthesis as the metabolic trigger. These data are of crucial importance. Despite the higher overall prevalence of hepatocellular carcinoma in males, we need tools and biomarkers to further stratify patients and offer better diagnosis and treatment options to both sexes.

Abstract

While the role of cholesterol in liver carcinogenesis remains controversial, hepatocellular carcinoma generally prevails in males. Herein, we uncover pathways of female-prevalent progression to hepatocellular carcinoma due to chronic repression of cholesterogenic lanosterol 14α-demethylase (CYP51) in hepatocytes. Tumors develop in knock-out mice after year one, with 2:1 prevalence in females. Metabolic and transcription factor networks were deduced from the liver transcriptome data, combined by sterol metabolite and blood parameter analyses, and interpreted with relevance to humans. Female knock-outs show increased plasma cholesterol and HDL, dampened lipid-related transcription factors FXR, LXRα:RXRα, and importantly, crosstalk between reduced LXRα and activated TGF-β signalling, indicating a higher susceptibility to HCC in aging females. PI3K/Akt signalling and ECM-receptor interaction are common pathways that are disturbed by sex-specific altered genes. Additionally, transcription factors (SOX9)2 and PPARα were recognized as important for female hepatocarcinogenesis, while overexpressed Cd36, a target of nuclear receptor RORC, is a new male-related regulator of ECM-receptor signalling in hepatocarcinogenesis. In conclusion, we uncover the sex-dependent metabolic reprogramming of cholesterol-related pathways that predispose for hepatocarcinogenesis in aging females. This is important in light of increased incidence of liver cancers in post-menopausal women.

1. Introduction

Metabolic abnormalities result in metabolic associated fatty liver disease (MAFLD), previously termed (NAFLD), which affects over 30% of the western population and has no approved drug therapies [1]. MAFLD can further develop into hepatocellular carcinoma (HCC) and both diseases are predicted to increase in the following decades [2], not only in men but also in women [3].
Among the lipid factors associated with progressive liver damage is also cholesterol. The cellular cholesterol concentrations are fine-tuned by a balance of synthesis, uptake, and efflux [4]. Genes of cholesterol synthesis differ in their susceptibility to loss of function mutations and can affect cell metabolism in different manners, from promoting cell survival towards manifesting in malignancies [5,6]. Linking cholesterol synthesis and liver disease remains controversial. Cholesterol synthesis was shown to support the growth of HCC lesions after depletion of fatty acid synthesis [7,8], indicating crosstalk between fatty acid and cholesterol pathways in carcinogenesis. Overexpression of the squalene epoxidase enzyme of cholesterol synthesis was identified as a driver of NAFLD-induced HCC [9]. On the other hand, blocking cholesterol synthesis at the lanosterol 14α-demethylase (CYP51) step leads to the prepubertal onset of liver injury with ductular reaction and fibrosis and a pronounced male prevalent liver dysfunction before puberty [10]. In young adults, the liver damage was more prominent among female mice due to diminished cholesterol esters and elevated sterol substrates [11].
The sex disparity in liver pathologies is far from conclusive. While some studies do not describe sex as a cofounding factor [7,9,12,13], there is increasing evidence that after the menopause MAFLD occurs at a higher rate in women compared to men [14]. Female and male livers are metabolically distinct organs, which was demonstrated also in mathematical models of the liver with unique regulators of sex-specific metabolic outcomes [15]. A recent review concluded that clinical and epidemiological studies frequently fail to address sex differences appropriately [14].
This study aimed to address the long-term effect of disrupted liver cholesterol synthesis where the Cyp51 knock-out (KO) model of aging mice was applied in light of the relevance to humans. Surprisingly, the tumors developed spontaneously mainly in female aging mice. We identified novel metabolic pathways and transcription factors that can be linked to metabolism associated with female-prone hepatocarcinogenesis. We provide mechanistic insights into the metabolic reprogramming in this type of liver cancer, with special attention to that of females.

2. Results

2.1. Long-Term Hepatocyte Deletion of Cyp51 from Cholesterol Synthesis Results in HCC

We wanted to understand how the long-term depletion of cholesterol synthesis influences liver physiology. We expected that due to the damaged livers in young adulthood [16], the mice would have a shortened life span. In contrary to expectations, mice lived at least until 24 months (M). Their livers did not become cirrhotic but started to develop liver tumors between 12M and 24M of age. Hepatomegaly was more pronounced in females’ KOs (Table S1). The histopathological examination of livers identified moderate to severe ductular reaction bridging the adjacent portal tracts, accompanied by mild inflammation (predominantly mononuclear cells; arrows in Figure 1A), and fibrosis. Importantly, fibrosis was more pronounced in females (Sirius red staining (SR); Figure 1B). Yellow pigment likely representing lipofuscin (arrow in Figure S1) was observed in areas of ductular reaction, sometimes surrounded by other hepatic cells, evaluated as macrophages (stars in Figure 1A).
The incidence of tumors by sex end genotype is presented in Table S2. Liver tumors were observed initially in 12M female and 18M male KO mice (Figure 2Aa). They were defined as macroscopically visible nodules (>1 mm) of various colours. At 24M, the incidence of tumors was 77.8% for females and 50% for males and were observed only in KO mice. In wild-types (WT), a single tumor was found in a 2-year male, histologically described as adenoma (Table S2). At 24M, the KO mice showed a decreased body weight compared to controls, significantly higher liver weight, and elevated relative liver weights. This was especially true for female mice with liver tumors (Figure S2 and Table S1). Tumors were histologically classified as eosinophilic/clear cell nodules, hepatocellular adenomas, or HCCs according to published criteria [17]. None of the tumors showed histological features of biliary tumor, e.g., biliary adenoma or cholangiocarcinoma.
Histologic evaluation is often sufficient for distinguishing the two types of primary hepatocellular tumors (HCC and cholangiocarcinoma), except in the case of poorly differentiated tumors. For further molecular analyses, we selected only tumors histologically classified as HCC, with features such as the absence of portal tracts and hepatic lobules, thickened hepatic plates, solid growth, focal tubular or pseudoglandular structures, focal necrosis, numerous apoptotic cells, and presence of cellular and nuclear polymorphism with brisk mitotic activity (Figure 2Ab,c).
Histologic evaluation is the gold standard in distinguishing the two most common primary hepatocellular tumors, namely hepatocellular carcinoma (HCC) and cholangiocarcinoma, except in the case of poorly differentiated tumors. Immunohistochemistry is a useful tool in determining the histogenesis of tumours. However, no single marker has been found that is completely specific for HCC. Among many different HCC markers, we decided for glutamine synthetase (GS), as it was recognized as the single best immunostain for identifying hepatocellular tumors in mice in a recent study [18]. As expected for normal liver tissue, positive GS staining was observed in pericentral hepatocytes (PC), but not in midzonal (M) or periportal (PP) hepatocytes (Figure 2Ba). Different staining patterns were observed in tumor tissue of Cyp51 KO mice, from scattered focal (Figure 2Bb), to diffuse staining of neoplastic hepatocytes (Figure 2Bc).

2.2. Sex-Dependent Progression from Metabolism Associated Liver Disease towards HCC

To learn about molecular changes based on long-term metabolic insult—the ablation of Cyp51 in hepatocytes—we assessed the hepatic expression of particular genes linked to cholesterol homeostasis during aging and measured multiple metabolites and markers in the blood of females and males. The results are presented graphically in Figure 3, with the most statistically significant changes between the groups (male, female, WT, Cyp51 KO) over all time points (12M, 18M, 24M) in Table S3. At the expression level, deregulation of cholesterol synthesis genes (Table S3) was, during aging, more pronounced in female livers (Figure 3A: Cyp51 KO—violet line, WT—black line) compared to those of males (Figure 3A: Cyp51 KO—green line, WT—blue line). The key regulatory gene of cholesterol synthesis Hmgcr is, over time, overexpressed in the KOs, particularly in females at 12M. While at 12M Cyp51 is downregulated in the livers of Cyp51 KO, this is not the case anymore at older age. We were initially intrigued by this finding. However, the Cyp51 is knocked-out only in the hepatocytes and not in other cell types of the liver that seem to complement for the Cyp51 loss and increase its expression. If we look into the cholesterol synthesis gene expression in the tumors at 24M (Table S4), all genes except Cyp51 were upregulated, as expected for feedback regulation after the block in cholesterol synthesis, again significantly only in the female knock-outs. Investigation of bile acid (BA) synthesis pathway (Figure 3B, Table S3) showed the expression of Cyp8b1, Cyp7a1 and Cyp27a1 decreased in the KOs, with no significant sex differences. At 24M in tumors, Cyp7a1 is significantly downregulated in the female KOs (Table S4), indicating a drop in the classical pathway of BA synthesis.
A major sex disparity was observed when monitoring the blood plasma biochemical parameters during aging. Total plasma cholesterol, free fatty acids (FFA), HDL (high-density lipoprotein cholesterol), and triglycerides (TG) showed opposite trends between sexes. Cholesterol and HDL were higher in KO females and were increasing by age, while they were lower in KO males compared to male controls, and decreasing with age. A significant increase in ALT in female Cyp51 KO mice is an additional indication of chronic hepatic injury (Figure 3C and Table S3), which is more prominent in females.
To quantify the primary molecular insults in the livers, we measured sterol intermediates of cholesterol synthesis together with cholesterol (Figure 3D and Table S3). In the KO mice with inactivated CYP51A1 enzyme, its substrates lanosterol and dihydrolanosterol were highly elevated (~100-fold and ~500-fold, respectively), proving disrupted cholesterol synthesis. Dihydrolanosterol was higher than lanosterol, indicating that 24-dehydrocholesterol reductase (DHCR24) enzyme functions normally. Downstream the pathway, desmosterol was not significantly different between Cyp51 KO and WT mice; however, zymostenol was higher in Cyp51 KO livers. Small or no changes between the sexes were observed in intermediate sterols. In contrast to that, sex represents a statistically significant parameter for liver cholesterol in the interaction analysis between sex and age. Similarly as for plasma cholesterol, hepatic cholesterol is also increased in female mice during ageing, while this was not the case for males (Figure 3C,D and Table S3).

2.3. Transcriptome Analysis Identified Sex Differences in Gene Expression Linked to Hepatocarcinogenesis

We used Affymetrix microarrays for gene expression profiling of mouse liver samples from 24M Cyp51 KO mice (24M KO) and controls (24M WT). We performed an analysis of liver tumors and the surrounding tissue that is responsible for tumorigenesis. To monitor the progression of tumorigenesis, data was compared to precancerous stages at 19 weeks of age (19W KO versus 19W WT) [16] and discussed in view of the early processes towards carcinogenesis (Tables S5–S7).
Hierarchical clustering of differentially expressed (DE) genes showed a clear difference between 24M WT and 24M KO mice for both sexes (Figure 4A). Importantly, a clear sex imbalance in the number of DE genes was observed (Figure 4B), with more genes (63.4%) deregulated in the female livers, including Ctnnb1, Gsk3b, Nova1, Lmna, Gstp1, and Mup20 (Table S8). Among the DE genes deregulated only in males (11.6%) were Trp53, Ezh2, Foxa1, and Rb1. In both sexes, a majority of DE genes were positively enriched (56% in females and 68% in males) (Figure 4B). Using NetworkAnalyst [19], the functionally enriched networks of DE genes showed that positively enriched genes were mostly involved in pathways related to cancer and the extracellular matrix (ECM)-receptor interaction pathway, while genes from metabolic pathways were negatively enriched, especially among female KOs (Figure 4C).
Additionally, we selected the top 10 statistically sex-specific significant down- and up-regulated genes (Table S9), whose biological functions were examined. Lipid metabolism genes are among the most down-regulated, however, with some distinct differences in lipid related pathways in 24M KO females and KO 24M males. Most down-regulated genes in KO males are from metabolism of fatty acids, such as fatty acid elongation, and PPAR signalling. In the female KOs, we observed downregulation of steroid hormone biosynthesis and related metabolic pathways, bile secretion, and also fatty acid metabolism (arachidonic acids). Drug metabolism is also downregulated in KO females. Retinol and linoleic metabolism were among top 10 in both sexes.
Oppositely, top 10 up-regulated genes in male KO liver are mainly involved indifferent cancer pathways, PPAR signalling and ECM receptor interactions, while the female list includes other upregulated genes in PPAR signalling, Wnt signalling and diabetes pathways.

2.4. Enriched KEGG Pathways and Transcription Factors (TFs) as Triggers of Female Prevalent Hepatocarcinogenesis

From 323 KEGG pathways, we discovered 209 that were enriched in livers of precancerous 19W KO mice and 230 in 24M KO mice with tumors, compared to 24M WTs (Table 1 and Table S6). The enriched pathways include the activation of HCC, apoptosis, PI3K/Akt signalling, pathways in cancer, and proteoglycans in cancer. The ECM-receptor interaction pathway was observed to be the most positively enriched in 19W KO and 24M KO females and males (Table 1), which is in line with the histological findings. Chronic injury induces numerous molecular alterations in hepatocytes and causes a gradual accumulation of extracellular matrix components (ECM) that might lead to a complete architectural reconstruction of hepatocytes and tumor development.
We discovered multiple signalling pathways that describe the sex-dependent metabolic reprogramming in liver cancer. Among them, TGF-β signalling was significantly elevated only in Cyp51 KO female livers, both at the precancerous stage and also later (Table 1). Additionally, the Wnt signalling and circadian entrainment pathways were activated only in the cancers of KO females (Table 1). We further explored whether these pathways could serve as triggers for female hepatocarcinogenesis in the Cyp51 KO model. The selected protein markers of the TGF-β and Wnt signalling were evaluated by qPCR and IHC (Figures S3 and S4). The TGF-β showed a positive reaction, particularly in female Cyp51 KO mice, indicating this is a good biomarker to follow progression towards HCC, particularly in females.
Among lipid metabolic pathways, we found fatty acid degradation and primary bile acid synthesis negatively enriched in KO females at 19W and 24M, while this was the case only for 24M KO males (Table 1), suggesting that the metabolic reprogramming of fatty acid and BA metabolism is active earlier in female carcinogenesis.
Considering the general alteration of DE genes and KEGG pathways in hepatocarcinogenesis, we also investigated the transcriptional landscape by transcription factor (TF) enrichment analysis. This analysis indicated 321 enriched TFs in the precancerous 19W KO livers and 142 in tumor containing 24M KO livers, both in comparison to WT littermates (Table S7). The enriched TFs with specific functions are listed in Table 2. From the positively enriched transcriptional activity SP1, C-JUN, JUND, C-FOS, AP1 and C/EBPβ were enriched already in the precancerous stage in both sexes, while PPARα, (SOX9)2 and NFATC1 are positively enriched only in females, starting already at the precancerous stage. HIF1α was also enriched in females in the precancerous stage, and together with SOX9, also in tumors of both sexes. The list of negatively enriched TFs is shorter. Downregulated CLOCK is characteristic for HCC, irrespective of the sex. HNF4α and RORC are negatively enriched in females already at the precancerous stage and in both sexes in tumors. NR1B1 is the only TF downregulated in the precancerous male livers and later in tumors of both sexes. LXRα:RXRα is a female-specific HCC marker.
When we considered the top-10 sex-specific enriched TFs based on the fold-change (Table S10), additional information was gathered. PPARα and the oncogenic isoform (SOX9)2, which were exposed previously (Table 2), are among top 10 positively enriched TFs by fold-change in females. FXR is, by fold-change, the most downregulated female-specific TF in HCC that was not exposed in the overall TF deregulation of Table 2. LXRα:RXRα were, again, among the most negatively enriched TFs in female tumors (Table S10). SOX9 is a part of the Wnt signalling and is expressed in the progenitor/precursor cells during hepatocyte regeneration following liver injury [20]. Expression of heterodimer complex LXRα:RXRα correlates with the amount of the retinoic acid and cholesterol derivates [21], which is line with diminished activity in KO female mouse livers. It was shown that chronic administration of PPARα agonists leads to liver tumor development in mice [22]. However, it is still unclear how this receptor changes in the human liver cancer tissue. An investigation of PPARα mRNA expression in human HCC indicated an increased mRNA expression of PPARα and its target genes. These results point out the possible association of activated PPARα with human carcinogenesis [23]. NFATC1, which was not on the list of sex-specific top 10 TFs, was also indicated as female-specific TF (Table 2). The NFATC1 with calcineurin overexpression plays an important role in the proliferation of HCC [24]. NFATC1 transcription can be activated by AP1 (including FOS–JUN), the activation of which has been already described in T cells [25].
RORC (RORγ) was one of the TFs in our focus since sterols after lanosterol are ligands of this transcription factor [26]. Downregulation of RORC signalling is thus important due to lowered intermediate sterols upon the Cyp51 block.
To validate whether downregulation of RORC has effects in hepatocarcinogenesis, we evaluated the expression of downstream genes regulated by this transcription factor. We found a majority of RORC target genes decreased in Cyp51 KO mice (Table S11), since RORC serves as their activator. However, the RORC target Cd36 was upregulated since RORC is known to repress this gene [10]. The RORC target G6pc, which is involved in the regulation of the PI3K/Akt signalling pathway, was negatively enriched in female KO mice, representing a part of the female metabolic reprogramming (Table S11).
A summary of sex-dependent reprogramming of metabolic pathways and corresponding transcription factors is shown in Figure 5. Briefly, Figure 5 indicates connections between signalling pathways, TFs and cellular processes, which were deregulated in both sexes (top half of figure) and specific for female KOs (bottom half of figure), leading to the inhibition of basal metabolism and an increased possibility for tumor development.
The TF networks were visualized with STRING [27] to determine their interactions. Different colours of nodes represent biological processes in which enriched TFs of Cyp51 KO mice are involved. Surprisingly, SP1 and JUND deregulated in precancerous and tumor livers of both sexes represent nodules with most protein–protein associations linked to hepatocarcinogenesis (Figure S5).

2.5. Metabolic and Transcriptional Changes after Disrupted Cyp51 from Cholesterol Synthesis Align with Hepatocarcinogenesis in Humans

We performed a comparative analysis of sex-dependent metabolic adaptations upon inactivation of CYP51A1 in the hepatocytes in light of the human literature data. Figure 6 represents the proposed carcinogenesis mechanism in our mouse model with exposed signalling, their specific sex-dependent targets and possible transcription regulators that might offer clues for expanding the understanding of the human sex-dependent HCC carcinogenesis. In the mouse model, the chronic metabolic insult of disrupted cholesterol synthesis (accumulation of hepatic sterols lanosterol and dihydrolanosterol) resulted in a variety of sex-specific metabolic adaptations. The hepatic cholesterol was increased in KO females and decreased in KO males. This observation serves as one of the important features in cancer cells. As fast proliferating cells, cancer cells require high levels of cholesterol for membrane biogenesis and other functional needs. For example, the cholesterol-derived oncometabolite 6-oxo-cholestan-3β,5α-diol was enriched in patients with breast cancer and subsequently promotes tumor growth [6].
As shown in Figure 6, disturbance of cholesterol synthesis leads to the repression of selected nuclear receptors (FXR, LXRα:RXR, LRH1 and NR1B1:RXRα) in female KO hepatocytes, which are increasingly appreciated targets in the prevention of inflammation and fibrosis [6,28]. The positive enrichment of PPARα, observed only in females, was also confirmed in humans with NAFLD [21], but its role in HCC patients remains unclear [29]. FXR negative enrichment and downregulation is also in line with the human HCC progression. Negative enrichment of FXR in females and downregulation of its target gene Ndrg2 is in line with the human data in patients and confirmed tumor suppressor role of FXR and NDRG2 in human cell lines [30].
Aligned with the latest human literature, metabolic and transcriptional changes of disturbed cholesterol metabolism in females pinpoints on a new sex-dependent mechanism, which has not been revealed so far. We proposed that a higher incidence of HCC in females results from the crosstalk of downregulated transcriptional activity of LXRα and activation of TGFβ signalling. Positive expression of TGF-β signalling targets Gstp1, Tgfrb1, Tgfrb2 and Bmp7 was identified only in female KO mice. Among several human malignancies, TGF-β in humans regulates the progression of female HCC through the promotion of PI3K/Akt signalling [31] and ECM-receptor interaction [32]. PI3K/Akt signalling and ECM-receptor interaction are common pathways that are disturbed by altered gene expressions in human cancers [33,34]. In our mouse model, both pathways were activated, but through different positively enriched genes in females (Col5a2 and Mmp9) compared to males (Cd36). In the context of lipid homeostasis, Raichur S et al. showed the opposite expression patterns, where repressed RORC activity induces expression of Cd36. In both animal and human studies, increased Cd36 expression in metabolically active tissue has been implicated in the regulation of fatty acid metabolism and lipid accumulation [35,36]. In HCC patients, the number of fatty acid responsive TFs, including HIF1α and FOS, can transcriptionally upregulate Cd36, which is closely associated with induction of epithelial-mesenchymal transition (EMT) [37].
Moreover, Gck, which is known as a one of RORC targets, was downregulated in Rorc−/− mice [38], while in female Cyp51 KO livers, it was upregulated. The identity of all TFs regulating Gck promoter activity has not been fully elucidated. One of the possible co-regulations of the liver Gck expression may be through the insulin, which can increase Gck mRNA via PI3K/Akt signalling, regulating the HIF1α [39] that were both found to be activated in KO mice.
In humans, the expression of SOX9 correlates with the severity of liver fibrosis patients [40] and correlates with the progression to cirrhosis and HCC [41]. In our case, the activated TF (SOX9)2, and its targets Ctgf and Nfe2l1 (Table S11), were deregulated only in female KOs.

3. Discussion

The hallmark of cancer cells is metabolic reprogramming, causing new cellular demands in selective survival and growth [42]. Cholesterol is a precursor of steroid hormones and bile acids and an essential component of cellular membranes. Its circulating levels are under healthy homeostatic conditions regulated by a balance between cellular cholesterol synthesis, dietary intake, and removal of excess cholesterol [43].
CYP51A1 is a gene from the late (post-squalene) part of the pathway [44]. The disruption of this gene in hepatocytes starts with morphological alterations, such as ductular reaction accompanied by mild inflammation. We show herein that this can end with malignant liver tumors in ageing mice that are first observed at 12M in females and 18M in males, with the highest incidence in the 24M female Cyp51 KO mice. In addition, humans show sex-dependent differences in liver pathologies. Women more commonly present with acute liver failure, benign liver lesions, autoimmune hepatitis, and toxin-mediated hepatotoxicity. Oppositely, malignant liver tumors and viral hepatitis are less commonly observed in women [45]. However, non-alcoholic steatohepatitis (NASH), as a more serious form of MALD, more commonly progresses to HCC in women [46]. Generally, women are less prone to developing HCC than men (due to hormone-dependent resistance), while after the menopause a sharp increase in HCC incidence is witnessed [47].
Based on our data, we propose that the repression of CYP51A1 enzyme can contribute to HCC development. In humans, whole-body deletion of CYP51A1 is not reported, since this gene is essential and is embryonically lethal in the mouse [48]. Only two studies are reporting the CYP51A1 deleterious mutations, which resulted in infant liver failure and death [49,50]. Numerous polymorphisms in CYP51A1 have been associated with different phenotypes in humans [51], but the consequences of CYP51A1 heterozygosity in human hepatocarcinogenesis need yet to be established. CYP51A1 has a low mutation rate, lower than other cholesterol synthesis genes [52]. The CYP51A1 rs37417517 variant was found to be associated with the female-specific breast cancer and variant rs229188 with female genetics [53,54]. Additionally, the COSMIC database revealed the presence of somatic CYP51A1 point mutations in tumor samples, some also with a predicted deleterious effect on activity (Tables S12 and S13).
The synthesis of cholesterol is a housekeeping pathway, which includes more than 100 genes associated with cholesterol synthesis and its regulation [55]. The block in cholesterogenic Cyp51 in KO mice is in our case specific in hepatocytes. Thus, about 40% of Cyp51 gDNA remained unchanged in the livers of KO mice, which lead to an approximally 60% decrease in mRNA and an 80% decrease in the CYP51 protein. Lorbek et al. indicated the small foci of stained periportal hepatocytes that potentially originated from the oval cell compartment by immunohistochemical staining of CYP51. Thus, the high intensity of ductular proliferation in ageing mice may explain an increase in Cyp51 mRNA in female and male KO livers from 18M onward [16]. Oppositely, the drop in expression of other genes of cholesterol synthesis could be a result of hepatic inadequacy due to chronic liver disease, which was also shown in fibrotic patients [56].
The latest studies indicate that post-squalene intermediates of cholesterol synthesis are also signalling molecules. Inhibiting different enzymes from the pathway results in the accumulation of sterol intermediates with different cellular effects. Sterol products can activate RORC signalling [26], promote oligodendrocyte formation [57], or activate LXRα and inhibit EGFR signalling [58]. Inhibitors of lanosterol synthase were identified as potential anticancer drug targets [59] while the accumulation of lanosterol can cause degradation of 3-hydroxy-3-methyl-glutaryl-coenzyme A reductase (HMGCR) [60]. Our sterol measurements show significant accumulation of CYP51A1 substrates lanosterol and dihydrolanosterol, with no sex differences. We observed lower amounts of desmosterol and cholesterol in KO livers. The presence of desmosterol in the liver may be explained by compensatory effects on lower levels of DHCR24, as it was shown in the brain [61]. The liver cholesterol may be obtained from the diet, to compensate the defected de novo synthesis of cholesterol [62]. The increase in cholesterol in female KO livers may arise due to the cholesterol uptake from plasma, since cholesterol is required for tumor cell proliferation and growth [56,63].
The main difference between the genotypes is likely in the accumulation of sterol substrates rather than a lack of downstream sterols, while the difference between the sexes is shown by the sex-dependent transcriptome response to sterol imbalance. Due to the lack of sterols between lanosterol and zymosterol, which are proposed RORC natural ligands, we observe at 24M the repression of the RORC signalling pathway in both sexes. The RORC-dependent mechanisms in cancerogenesis were recently described. Oh et al. reported that RORC expression is decreased in basal-like subtype cancers and inversely correlated with histological grade and cancer drivers in breast cancer cohorts [64]. In breast cancer cell lines, RORC was shown to repress TGF-β/EMT-signalling and cancerogenic pathways, which indicates that, also in Cyp51 KO mice, decreased activity of RORC could contribute to increased TGF-β/EMT signalling. There is no data yet that would link SOX9 expression to sterols. (SOX9)2 is upregulated, particularly in KO females, and further research is needed to evaluate whether this is a direct effect of accumulating lanosterol and dihydrolanosterol.
Injured hepatocytes in Cyp51 KO mice produce different ECMs, exemplified by the enrichment of laminin and collagen genes (Col4a1, Lamc1, etc.) and a strongly positively enriched ECM-receptor interaction pathway in hepatocarcinogenesis of both sexes. In some way, these alterations could contribute to microenvironmental reorganization and induction of EMT that positively influences tumor progression [32]. A higher rate of cancer progression in female Cyp51 KO mice could be explained by the female-specific interplay of positively enriched TGF-β signalling. Activated TGF-β represents a central regulator of chronic liver disease and could contribute to the progression of all disease stages, both in animal models and in humans [65]. Interestingly, important crosstalk between TGF-β and Cyp51 gene expression has been shown in the Tgf-βfl/fl;Wnt-Cre mouse model [66].
Furthermore, the results of our analysis identified the key transcriptional regulators above the differentially expressed molecular pathways: positively enriched SOX9, SP1, and negatively enriched RORC. SOX9 is an important candidate for the cancer marker required for activation of TGF-β/Smad signalling [67]. SOX9 expression, especially its isoform 2 ((SOX9)2), can predict liver disease progression [40], but was so far not described as being associated more with the female type of carcinogenesis. Similarly to SOX9, the tumor progression driver SP1 is also a key factor for the induction of TGF-β during inflammation [68] and its overexpression in Cyp51 knock-outs relates to EMT. The impact of higher disease progression in females is reflected also in the interplay of transcription regulator SMAD3 on TGF-β signalling [69].
Increased plasma total cholesterol and HDL, and decreased plasma TG and FFA, are observed only in female KO mice. This is a possible effect of induced sex hormone-binding globulin (SHBG) by specific targets of an estrogen signalling pathway [70]. SHBG positively correlates with HDL and negatively with TGs and has been already associated with different cancers in ageing humans [71,72]. Supporting other hepatocarcinogenesis studies [73,74], the outcome of altered cholesterol and BA homeostasis was a chronic liver injury in our Cyp51 KO mouse model with widespread fibrosis, inflammation, and malignant transformation due to toxic sterol accumulation. Additionally, sterol accumulation, in our case the consequence of inhibition of CYP51A1 [16], could be a cause of the observed widespread increase in hepatic progenitors (also known as a ductular reaction) that may also differentiate into mature hepatocyte-like cells with oncogenic potential [75].

4. Materials and Methods

Detailed materials and methods are described in Appendix A. The experimental design of the study is shown in Figure A1.

4.1. Animal Study and Samples Collection

The generation of homozygous hepatocyte-specific Cyp51 knock-out (Cyp51 KO) mice has been reported previously [10,16]. Animals were bred under a 12 h light and dark cycle (7:00 a.m. until 7:00 p.m. light), at 22 ± 1 °C, humidity 55 ± 10%, on standard rodent chow (diet 1324, Altromin, Lage, Germany) and acidified tap water (pH = 3) ad libitum. They were housed in groups of 3–5 mice per cage. The results of routine health microbiological monitoring are shown in Table A1. Mice were sacrificed at different ages (12M, 18M and 24M) with cervical dislocation between 11 and 13 h (to minimize the circadian influence), after a 4–6 h of fasting (food withdrawal at 7 h a.m.). Blood was taken immediately after cervical dislocation from the right ventricle, collected into heparin-coated Vacuette MiniCollect® 1 mL Plasma Tubes (Greiner Bioone, Frickenhausen, Germany) and centrifuged for 4 °C, 3000× g for 15 min. Mice were weighed, organs (liver, kidney, spleen, heart, gonads) removed and patho-morphologically examined. For histology, the left lateral lobus was fixed in 4% formalin and embedded in paraffin for further histological analysis. Other parts of livers were cut into thin slices and snap-frozen in liquid nitrogen. The material was freshly frozen and stored at −80 °C for subsequent analysis.

4.2. Ethics Statement

Mouse experiments were approved by the Veterinary Administration of the Republic of Slovenia (license numbers 34401-31/2011/4 and 34401-52/2012/3) and were conducted in agreement with the National Institute of Health guidelines for work with laboratory animals, national legislation and Directive 2010/63/EU on the protection of animals used for scientific purposes.

4.3. Histological Analysis

For standard histological analyses, formalin-fixed paraffin-embedded tissue was used, sectioned to 4–5 µm and deparaffinized at 70 °C for 10 min and washed in two changes of Xylene, 100% ethanol, 95% ethanol, 70% ethanol and twice in acidified water. Hematoxylin and Eosin (H&E) stained liver sections of each mouse were scored for size, shape, and polymorphism of the hepatocytes, granulation of cytoplasm, presence of portal and parenchymal inflammation, ductular reaction, cholestasis, and presence of any type of nodule. The intensity of inflammation and ductular reaction was graded as: 0—normal, +—mild, ++—moderate and +++—severe. The extent of ductular reaction was assigned: p—proliferation around portal tracts, p-mz—proliferation extending from one portal tract to the middle zone, and p-p—bile ducts were extending from one portal field to another (bridging). For evaluation of liver fibrosis, Sirius red (SR) staining was performed by incubating deparaffinized sections in 0.1% Sirius Red solution for 1h (SR; 0.1% direct red 80, 1.2% picric acid in water), briefly destaining in diluted acetic acid, dehydrating in 70% ethanol, 95% ethanol, 100% ethanol and twice in Xylene and fixed with Roti Histokitt II (Carl Roth GmbH + Co. KG, Karlsruhe, Germany). Fibrosis in SR stained liver samples were graded: 0 as absent, + mild periportal, ++ moderate, and +++ as very strong fibrosis with bridging. Tumors, identified as macroscopically visible nodules (>1 mm) of various colors observed at autopsy, were histologically classified based on Mouse tumor classification criteria [17] as eosinophilic or clear cells nodules, adenomas, cholangiomas, cholangiocellular carcinoma or hepatocellular carcinoma (HCC).

4.4. Immunohistochemistry

To evaluate the expression of selected markers in liver samples submitted to the microarray experiment, immunolabeling for glutamine synthetase, β-catenin, and TGF-β1 was performed. Liver tissue sections were deparaffinized and subjected to degrading alcohol gradient and treated for antigen retrieval before staining. Afterwards, 3% hydrogen peroxide (H2O2) was used for endogenous peroxidase quenching, non-specific staining was blocked using 5% normal goat serum (Sigma-Aldrich, St. Louis, MO, USA) for 1 h at room temperature. Incubation with primary, anti-glutamine synthetase (BD Biosciences, San Jose, CA, USA; 1:1000 dilution), anti-β-catenin (BD Biosciences, San Jose, CA, USA; 1:1000 dilution) and anti-TGF-β1 (Promega, Madison, WI, USA; 1:500 dilution) antibodies in 1% goat serum in 0.1% Tris buffered saline with Tween 20 (TBST) was performed overnight at 4 °C. We used the DAKO EnVision Detection System (Agilent Technologies DAKO, Glostrup, Denmark) for antibody detection. The sections were subsequently counterstained with hematoxylin.

4.5. Analysis of Plasma Parameters

Total and HDL cholesterol, free fatty acids (FFA), triglycerides (TG), alanine aminotransferase (ALT), and aspartate aminotransferase (AST) were analyzed with the Architect ci8200 analyzer (Abbott Diagnostics, Abbott Park, IL, USA). The concentration of lipid parameters were measured as mmol/L and a linear regression model Age+Genotype+Sex+Age:Sex was used to fit the concentrations of individual lipids. The activity of ALT and AST were given in μkat/L and a linear regression model Age+Genotype+Sex+Genotype:Sex was used to fit their activity individually. ANOVA type III test was used to estimate the statistical significance of the effects and the Holm method was used to control the family-wise error rate (FWER) at α = 0.05.

4.6. Liver Sterol Analysis

Sterols isolation: Sterol intermediates were isolated with the protocol described previously [76]. Folch solution with samples corresponding to 80 mg of liver tissue was transfer to fresh tubes and 200 ng of internal standard Lathosterol-D7 (Avanti Polar Lipids, Alabaster, AL, USA) was added. Final samples were dissolved in 250 µL of methanol and transferred to HPLC vials for analysis. LC-MS analysis was performed as described [77]. Briefly, sterols were separated on combined columns of Luna® 3 µm PFP (2) 100 Å (100 mm and 150 mm length) on a Shimadzu Nexera XR HPLC. For mobile phase, 80% methanol, 10% water, 10% 1-propanol, and 0.05% Formic acid were used in the isocratic condition. HPLC was coupled with the Sciex Triple Quad 3500 mass spectrometer for detection. For ionization, APCI (Atmospheric-pressure chemical ionization) was used in positive mode and detection was made in MRM (Multiple reaction monitoring) mode. The sample concentration was normalized on internal standard (Lathosterol-d7) and sterol concentrations calculated to corresponding standards from Avanti Lipids. We measured the following sterol intermediates, with MRM in brackets: Zymosterol* (367/215), 24-dehydrolathosterol* (367/215), 7-dehydrodesmosterol* (365/199), Desmosterol (367/215), Zymostenol (369/215), Lathosterol-d7 (376/215), FF-MAS* (393/214), T-MAS* (395/243), Cholesterol (369/215), Lanosterol (409/191) and 24,25-dihydrolanosterol (411/191). *Sterol concentration was too low to quantify.

4.7. Gene Expression Profiling and Network Analysis

Transcriptome analysis of liver samples of Cyp51 WT and KO mice of both sexes was performed using Affymetrix GeneChip™ Mouse Gene 2.0 ST arrays and validated by qPCR (primers listed in Table S14; further details in Appendix A. The same mice and groups were used for plasma analysis (Table A2). For qPCR, a linear regression model Age+Genotype+Sex+Age:Genotype was used to fit normalized expression data of cholesterol-related genes Hmgcr, Sqle, Lss, Cyp51, Nsdhl and Dhcr7, and a reduced additive model Age+Genotype+Sex was used to fit data of fatty acid-related genes Cyp7a1, Cyp8b1 and Cyp27a1. ANOVA type III test was used to estimate the statistical significance of the effects and Holm method was used to control FWER at α = 0.05.
A total of 20 Affymetrix GeneChip™ Mouse Gene 2.0 ST Arrays (Affymetrix, Santa Clara, CA, USA) with liver samples were processed. From each of the Cyp51 KO livers, two samples were taken: tumor and surrounding tissue. Altogether, 6 tumor and 6 corresponding surrounding-tissue samples were obtained from Cyp51 KO mice (Table A3). The data were analysed using R/limma, controlling the false discovery rate at α = 0.05, and deposited to Gene Expression Omnibus (GEO) under accession number GSE127772. KEGG PATHWAY and TRANSFAC databases were used for functional enrichment studies. Gene expression data from 19 weeks (W) [16] WT and KO mice of both sexes were used to determine processes that contribute to tumor progression from 19W towards 24M. Network diagrams were created using NetworkAnalyst and STRING. The proposed mechanistic scheme was designed from data obtained from mice experiments and aligned to human literature data to increase their relevance.

5. Conclusions

The hallmark of cancer cells is metabolic reprogramming, causing new cellular demands in selective survival and growth. Metabolic reprogramming of hepatic cholesterol synthesis leads to female prevalent hepatocellular carcinoma during aging. Inactivation of Cyp51 affected multiple signalling pathways and transcription factors in a sex-dependent manner, resulting in metabolic-related HCC. The female-specific cholesterol-related mechanism of hepatocarcinogenesis is very similar to the human data and is in line with the fact that the liver is a sex-dimorphic metabolic organ and a major organ of cholesterol homeostasis. Since hepatocellular carcinoma is increasing in post-menopausal women, the female-specific aspects of liver cancer progression, diagnosis and treatment should become a central focus of research.

Supplementary Materials

The following are available online at https://www.mdpi.com/2072-6694/12/11/3302/s1, in files (1) Supplementary Figures and (2) Supplementary Tables. Figure S1: Patho-histological features of Cyp51 KO mice. Figure S2: The KO mice liver and body weights at different ages, both sexes. Figure S3: Immunohistochemical expression of TGF-β1 and β-catenin. Figure S4: qPCR expression profiles of TGF-β and Wnt signalling markers. Figure S5: TF networks of 24M KO mice. Table S1: List of gene sequences for qPCR validation. Table S2: Female Cyp51 KO mice had a higher incidence of liver tumors as males. Table S3: Differences between mouse liver/body/relative weights. Table S4: A heat-map of cholesterol and bile acid synthesis DNA-microarray gene expression in male and female KO mice. Table S5: The long-term ablation of Cyp51 in hepatocytes influence at gene and metabolic levels. Table S6: DNA-microarray gene expression tables. Table S7: KEGG pathways expression tables deduced from DNA-microarray gene expression pattern. Table S8: Transcription factors expression tables deduced from DNA-microarray gene expression pattern. Table S9: A clear sex imbalance in gene expression in 24M KOs. Table S10: List of sex-specific top 10 statistically significant down- and up-regulated genes (deduced from Table S6) in Cyp51 KO mice, aged 19 weeks and 24 months. Table S11: List of sex-specific top 10 statistically significant down- and up-regulared TFs (deduced from Table S8) in Cyp51 KO mice, aged 19 weeks and 24 months. Table S12: Target genes of TFs RORC and SOX9. Table S13: COSMIC database information of CYP51A1 point mutations in humans. Table S14: The positive data (mutated samples) for selected gene CYP51A1 from COSMIC database. The datasets generated during the current study are available in the Gene Expression Omnibus database (https://www.ncbi.nlm.nih.gov/geo/) under accession number GSE127772 and GSE58271.

Author Contributions

Conceptualization, D.R.; methodology, K.B.C., G.L., M.P., J.J., Ž.U. and C.S.; software, P.J.; validation, K.B.C. and C.S.; formal analysis, P.J.; investigation, K.B.C.; resources, D.R., and M.M.-S.; data curation, P.J.; writing—original draft preparation, K.B.C., M.P., J.J., T.R., and D.R.; writing—review and editing, all authors; visualization, K.B.C.; supervision, D.R.; project administration, D.R.; funding acquisition, D.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Slovenian Research Agency (ARRS) program grants P1-0104, P1-0390, and projects J1-9176 and BI-DE/17-19-8. Kaja Blagotinšek Cokan, Gregor Lorbek, and Žiga Urlep were supported by the graduate fellowship of ARRS.

Acknowledgments

We thank Ksenija Kodra, Helena Klavžar, Marjan Kužnik, and Doris Mahn for their help with experiments. We acknowledge Simon Horvat and Rok Keber for their previous contribution to the development of the Cyp51 knock-out mouse model. Thanks also to Nejc Nadižar for help in statistical data evaluation and figure preparation and Petra Ivanuša for help with sterol isolation. Thanks to collaborators Irena Vovk and Mitja Križman for LC-MS method development.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

CYP51Lanosterol 14α-demethylase
MAFLDMetabolic Associated Fatty Liver Disease
NAFLDNon-Alcoholic Fatty Liver Disease
HCCHepatocellular Carcinoma
WTWild-Type
KOKnock-Out
MMonth
WWeek
H&EHematoxylin and Eosin
SRSirius Red
FFAFatty Free Acids
TGTriglycerides
BABile Acids
CHOLCholesterol
HDLHigh-Density Lipoprotein
ALTAlanine Transaminase
ASTAspartate Transaminase
LC-MSLiquid Chromatography-Mass Spectrometry
HMGCR3-hydroxy-3-methyl-glutaryl-coenzyme A reductase
DHCR2424-dehydrocholesterol reductase
GEOGene Expression Omnibus
DEDifferentially Expression
ECMExtracellular Matrix
KEGGKyoto Encyclopedia of Genes and Genomes
EMTEpithelial to Mesenchymal Transition
TFTranscription Factor

Appendix A

Figure A1. Experimental design. Experiments were performed in female and male Cyp51 liver conditional knock-out mice at age of 12, 18, and 24 months, and the wild-type controls. Collected data from livers and blood samples were aligned to data of the 19 weeks old mice.
Figure A1. Experimental design. Experiments were performed in female and male Cyp51 liver conditional knock-out mice at age of 12, 18, and 24 months, and the wild-type controls. Collected data from livers and blood samples were aligned to data of the 19 weeks old mice.
Cancers 12 03302 g0a1
Table A1. List of negatively tested microorganisms in mice in routine health microbiological monitoring.
Table A1. List of negatively tested microorganisms in mice in routine health microbiological monitoring.
List of Negatively Tested Microorganisms
Citrobacter rodentiumMinute virus of mice
Corynebacterium kutscheriMouse hepatitis virus
Klebsiella oxytocaMouse parvovirus (rVP2)
Klebsiella pneumoniaeMouse rotavirus/EDIM
Pseudomonas aeruginosaMycoplasma pulmonis
Salmonela spp. Pneumonia virus of mice
Staphylococcus aureusReovirus type 3
Streptobacillus moniformisSendai virus
Streptococci beta-haemolytic Group ATheiler’s encephalomyelitis virus (GD VII)
Streptococci beta-haemolytic Group BMyobia musculi/Radfordia sp.
Streptococci beta-haemolytic Group COther ectoparasites
Streptococci beta-haemolytic Group GAspiculuris tetraptera
Streptococcus pneumoniaeCryptosporidium spp.
Adenovirus FLEntamoeba spp.
Adenovirus K87Giardia spp.
Clostridium piliformeHelicobacter bilis
Ectromelia virusSpironucleus spp.
General parvovirus (rNS-1)Syphacia obvelata
Lymphocytic choriomeningitisvirus

Plasma Analysis

Mice were placed into groups according to their age (12, 18 and 24 months, 24 mice per group), genotype (37 Cyp51 WT, 35 Cyp51 KO), and sex (40 female, 32 male). 72 mouse samples (4-7 samples per group) were analysed (Table A2).
Table A2. Number of mouse samples for plasma and RT-qPCR gene expression analysis.
Table A2. Number of mouse samples for plasma and RT-qPCR gene expression analysis.
Mouse SamplesAge
GenotypeSex12M18M24M
Cyp51 WTfemale777
male556
Cyp51 KOfemale577
male754

Gene Expression Analysis

The same mice and groups were used as for plasma analysis (Table A2). Total RNA was isolated from 30 µg of the liver using TRI Reagent (Sigma-Aldrich, St. Louis, MO, USA) procedure. RNA concentration and purity of each sample were assayed using NanoDrop 1000 Spectrophotometer (Thermo Fischer Scientific, Waltham, MA, USA). The RNA quality was checked with Agilent 2100 BioAnalyzer (Agilent Technologies, Santa Clara, CA, USA). Before cDNA synthesis, all liver RNA samples (2 µg) were treated with amplification grade DNase I (Roche, Basel, Switzerland), and reverse-transcribed with the Transcriptor Universal cDNA Master (Roche, Basel, Switzerland) according to manufacturer’s instructions.
Real-time quantitative reverse transcription polymerase chain reaction (RT-qPCR) was performed using a Roche LightCycler 480 (Roche, Basel, Switzerland). In all experiments, the PCR reaction consisted of 2.5 µL SYBR Green I Master (Roche, Basel, Switzerland), 0.75 µL of cDNA template, 0.6 µL of 2.5 µM random primer mix and 1.15 µL of RNAse-free water in a final volume of 5 µL. All experiments were carried out in three technical replicates for each of the sample using 384-well plates. PCR was performed with the following parameters: 10 min incubation at 95 °C followed by 45 cycles of 10 s at 95 °C, 20 s at 58 °C and 20 s at 72 °C. Utp6 and Hmbs were chosen as internal reference genes in mice liver for normalization of gene expression data using NormFinder [78]. The relative expression ratio was calculated using the ddCp as previously described [79]. The list of primer sequences used in RT-qPCR is provided under Table S1: oligonucleotide sequences.

Microarray-Based Gene Expression Profiling

Table A3. Number of 24M mouse samples for DNA-microarray analysis. ↔ indicates related samples, i.e., different tissue samples obtained from individual mice.
Table A3. Number of 24M mouse samples for DNA-microarray analysis. ↔ indicates related samples, i.e., different tissue samples obtained from individual mice.
24M Mouse SamplesGenotype
Cyp51 WTCyp51 KO
Tissue
NormalSurroundingHCC
SexFemale 4 4 4
Male4 22
The experiment was performed with 500 ng starting concentration of total RNA with a purity (absorbance ratio A260/280) >1.8 and quality (RIN) >4.5. The procedure was done as described previously [10,16] and considering the manufacturer’s instructions. Briefly, after overnight hybridization at 45 °C and 60 rpm/min microarrays were washed and stained using GeneChip Fluidics Station 450 and scanned on Affymetrix GeneChip Scanner 3000 7G. Image analysis and first quality check were done using the Affymetrix Expression Console™ version 1.3.
Further quality control and analysis were performed using R and Bioconductor software packages. RMA algorithm from package oligo [80] was used for the normalization of raw expression data. Package arrayQualityMetrics [81] was used to check additional quality control and assessment of potential outliers before and after normalization. Raw (CEL) as well as normalized data were deposited to GEO under accession number GSE127772 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE127772). Package limma [82] was used to fit individual normalized gene expression and gene set enrichment data using a linear regression model Genotype+Sex+Genotype:Sex. Empirical Bayes statistics was used to estimate the significance of the effects for each gene (gene set) and Benjamini-Hochberg procedure was used to control false discovery rate (FDR) at α=0.05 accounting for the number of tested genes (gene sets) and to infer differential expression of genes (enrichment of gene sets).
KEGG PATHWAY [82] and TRANSFAC database were used for functional enrichment studies [83]. Gene sets containing over 5 elements were constructed and tested for enrichment using the PGSEA package [84]. In the case of transcription factor enrichment, factors were merged based on their ID irrespective of their binding sites.
Gene expression data from 19 weeks (19W) old hepatocyte-specific Cyp51 KO mice (Cyp51flox/flox; Alb-Cre or LKO; i.e., KO) and their WT littermates (Cyp51flox/flox or LWT; i.e., WT) following standard rodent diet (low-fat no-cholesterol, LFnC) was analysed to determine processes that contribute to tumor progression from 19W towards 24M. Data is available from GEO under accession number GSE58271 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi) [16]. Data was preprocessed, annotation updated and analyzed using the same procedure as described above. A linear regression model Genotype+Sex+Genotype:Sex was used to fit individual normalized gene expression and gene set enrichment data from mouse samples shown in Table A4.
Table A4. Number of 19W mouse samples from GSE58271 that were included to determine processes that contribute to tumor progression from 19W towards 24M.
Table A4. Number of 19W mouse samples from GSE58271 that were included to determine processes that contribute to tumor progression from 19W towards 24M.
19W Mouse SamplesGenotype
WTKO
SexFemale33
Male33
Network diagrams of differences in gene expression in females and males (Figure 4C) were created using the NetworkAnalyst program [19]. For the generation and visualization of transcription factor networks (Figure S5), the list of female and male DEG was mapped to the ˝database˝ libraries of the STRING online software [27] using medium confidence score (0.4).

References

  1. Eslam, M.; Sanyal, A.J.; George, J.; International Consensus Panel. MAFLD: A consensus-driven proposed nomenclature for metabolic associated fatty liver disease. Gastroenterology 2020. [Google Scholar] [CrossRef]
  2. Nakagawa, H.; Hayata, Y.; Kawamura, S.; Yamada, T.; Fujiwara, N.; Koike, K. Lipid Metabolic Reprogramming in Hepatocellular Carcinoma. Cancers 2018, 10, 447. [Google Scholar] [CrossRef] [Green Version]
  3. Burra, P.; De Martin, E.; Gitto, S.; Villa, E. Influence of age and gender before and after liver transplantation. Liver Transplant. Off. Publ. Am. Assoc. Study Liver Dis. Int. Liver Transplant. Soc. 2013, 19, 122–134. [Google Scholar] [CrossRef]
  4. Chua, N.K.; Coates, H.W.; Brown, A.J. Cholesterol, cancer, and rebooting a treatment for athlete’s foot. Sci. Transl. Med. 2018, 10. [Google Scholar] [CrossRef]
  5. Horvat, S.; McWhir, J.; Rozman, D. Defects in cholesterol synthesis genes in mouse and in humans: Lessons for drug development and safer treatments. Drug Metab. Rev. 2011, 43, 69–90. [Google Scholar] [CrossRef]
  6. Huang, B.; Song, B.L.; Xu, C. Cholesterol metabolism in cancer: Mechanisms and therapeutic opportunities. Nat. Metab. 2020, 2, 132–141. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Che, L.; Chi, W.; Qiao, Y.; Zhang, J.; Song, X.; Liu, Y.; Li, L.; Jia, J.; Pilo, M.G.; Wang, J.; et al. Cholesterol biosynthesis supports the growth of hepatocarcinoma lesions depleted of fatty acid synthase in mice and humans. Gut 2020, 69, 177–186. [Google Scholar] [CrossRef] [Green Version]
  8. Zhu, W.W.; Lu, M.; Wang, X.Y.; Zhou, X.; Gao, C.; Qin, L.X. The fuel and engine: The roles of reprogrammed metabolism in metastasis of primary liver cancer. Genes Dis. 2020, 7, 299–307. [Google Scholar] [CrossRef]
  9. Liu, D.; Wong, C.C.; Fu, L.; Chen, H.; Zhao, L.; Li, C.; Zhou, Y.; Zhang, Y.; Xu, W.; Yang, Y.; et al. Squalene epoxidase drives NAFLD-induced hepatocellular carcinoma and is a pharmaceutical target. Sci. Transl. Med. 2018, 10. [Google Scholar] [CrossRef] [Green Version]
  10. Urlep, Z.; Lorbek, G.; Perse, M.; Jeruc, J.; Juvan, P.; Matz-Soja, M.; Gebhardt, R.; Bjorkhem, I.; Hall, J.A.; Bonneau, R.; et al. Disrupting Hepatocyte Cyp51 from Cholesterol Synthesis Leads to Progressive Liver Injury in the Developing Mouse and Decreases RORC Signalling. Sci. Rep. 2017, 7, 40775. [Google Scholar] [CrossRef] [PubMed]
  11. Lorbek, G.; Urlep, Z.; Rozman, D. Pharmacogenomic and personalized approaches to tackle nonalcoholic fatty liver disease. Pharmacogenomics 2016, 17, 1273–1288. [Google Scholar] [CrossRef] [Green Version]
  12. Wilson, M.A.; Buetow, K.H. Novel Mechanisms of Cancer Emerge When Accounting for Sex as a Biological Variable. Cancer Res. 2020, 80, 27–29. [Google Scholar] [CrossRef]
  13. Li, C.H.; Haider, S.; Shiah, Y.J.; Thai, K.; Boutros, P.C. Sex Differences in Cancer Driver Genes and Biomarkers. Cancer Res. 2018, 78, 5527–5537. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Lonardo, A.; Nascimbeni, F.; Ballestri, S.; Fairweather, D.; Win, S.; Than, T.A.; Abdelmalek, M.F.; Suzuki, A. Sex Differences in Nonalcoholic Fatty Liver Disease: State of the Art and Identification of Research Gaps. Hepatology 2019, 70, 1457–1469. [Google Scholar] [CrossRef]
  15. Cvitanovic Tomas, T.; Urlep, Z.; Moskon, M.; Mraz, M.; Rozman, D. LiverSex Computational Model: Sexual Aspects in Hepatic Metabolism and Abnormalities. Front. Physiol. 2018, 9, 360. [Google Scholar] [CrossRef]
  16. Lorbek, G.; Perse, M.; Jeruc, J.; Juvan, P.; Gutierrez-Mariscal, F.M.; Lewinska, M.; Gebhardt, R.; Keber, R.; Horvat, S.; Bjorkhem, I.; et al. Lessons from hepatocyte-specific Cyp51 knockout mice: Impaired cholesterol synthesis leads to oval cell-driven liver injury. Sci. Rep. 2015, 5, 8777. [Google Scholar] [CrossRef] [Green Version]
  17. Deschl, U.; Cattley, R.C.; Harada, T.; Kuttler, K.; Hailey, J.R.; Harting, F.; Leblanc, B.; Marsman, D.S.; Shirai, T. Liver, Gallbladder, and Exocrine Pancreas. In International Classification of Rodent Tumors: The Mouse; Mohr, U., Ed.; Springer: Berlin/Heidelberg Germany; New York, NY, USA, 2001; pp. 59–86. [Google Scholar]
  18. Salleng, K.J.; Revetta, F.L.; Deane, N.G.; Washington, M.K. The Applicability of a Human Immunohistochemical Panel to Mouse Models of Hepatocellular Neoplasia. Comp. Med. 2015, 65, 398–408. [Google Scholar]
  19. Xia, J.; Benner, M.J.; Hancock, R.E. NetworkAnalyst--integrative approaches for protein-protein interaction network analysis and visual exploration. Nucleic Acids Res. 2014, 42, W167–W174. [Google Scholar] [CrossRef] [Green Version]
  20. Leung, C.O.; Mak, W.N.; Kai, A.K.; Chan, K.S.; Lee, T.K.; Ng, I.O.; Lo, R.C. Sox9 confers stemness properties in hepatocellular carcinoma through Frizzled-7 mediated Wnt/beta-catenin signaling. Oncotarget 2016, 7, 29371–29386. [Google Scholar] [CrossRef] [Green Version]
  21. Steensels, S.; Qiao, J.; Ersoy, B.A. Transcriptional Regulation in Non-Alcoholic Fatty Liver Disease. Metabolites 2020, 10, 283. [Google Scholar] [CrossRef]
  22. Mello, T.; Materozzi, M.; Galli, A. PPARs and Mitochondrial Metabolism: From NAFLD to HCC. PPAR Res. 2016, 2016, 7403230. [Google Scholar] [CrossRef]
  23. Kurokawa, T.; Shimomura, Y.; Bajotto, G.; Kotake, K.; Arikawa, T.; Ito, N.; Yasuda, A.; Nagata, H.; Nonami, T.; Masuko, K. Peroxisome proliferator-activated receptor alpha (PPARalpha) mRNA expression in human hepatocellular carcinoma tissue and non-cancerous liver tissue. World J. Surg. Oncol. 2011, 9, 167. [Google Scholar] [CrossRef] [Green Version]
  24. Wang, S.; Kang, X.; Cao, S.; Cheng, H.; Wang, D.; Geng, J. Calcineurin/NFATc1 pathway contributes to cell proliferation in hepatocellular carcinoma. Dig. Dis. Sci. 2012, 57, 3184–3188. [Google Scholar] [CrossRef]
  25. Li, L.; Yu, J.; Duan, Z.; Dang, H.X. The effect of NFATc1 on vascular generation and the possible underlying mechanism in epithelial ovarian carcinoma. Int. J. Oncol. 2016, 48, 1457–1466. [Google Scholar] [CrossRef]
  26. Santori, F.R.; Huang, P.; van de Pavert, S.A.; Douglass, E.F., Jr.; Leaver, D.J.; Haubrich, B.A.; Keber, R.; Lorbek, G.; Konijn, T.; Rosales, B.N.; et al. Identification of natural RORgamma ligands that regulate the development of lymphoid cells. Cell Metab. 2015, 21, 286–298. [Google Scholar] [CrossRef] [Green Version]
  27. Szklarczyk, D.; Morris, J.H.; Cook, H.; Kuhn, M.; Wyder, S.; Simonovic, M.; Santos, A.; Doncheva, N.T.; Roth, A.; Bork, P.; et al. The STRING database in 2017: Quality-controlled protein-protein association networks, made broadly accessible. Nucleic Acids Res. 2017, 45, D362–D368. [Google Scholar] [CrossRef]
  28. Smeuninx, B.; Boslem, E.; Febbraio, M.A. Current and Future Treatments in the Fight Against Non-Alcoholic Fatty Liver Disease. Cancers 2020, 12, 1714. [Google Scholar] [CrossRef]
  29. Xiao, Y.B.; Cai, S.H.; Liu, L.L.; Yang, X.; Yun, J.P. Decreased expression of peroxisome proliferator-activated receptor alpha indicates unfavorable outcomes in hepatocellular carcinoma. Cancer Manag. Res. 2018, 10, 1781–1789. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Kovac, U.; Skubic, C.; Bohinc, L.; Rozman, D.; Rezen, T. Oxysterols and Gastrointestinal Cancers Around the Clock. Front. Endocrinol. 2019, 10, 483. [Google Scholar] [CrossRef] [Green Version]
  31. Moren, A.; Bellomo, C.; Tsubakihara, Y.; Kardassis, D.; Mikulits, W.; Heldin, C.H.; Moustakas, A. LXRalpha limits TGFbeta-dependent hepatocellular carcinoma associated fibroblast differentiation. Oncogenesis 2019, 8, 36. [Google Scholar] [CrossRef] [PubMed]
  32. Yang, J.D.; Nakamura, I.; Roberts, L.R. The tumor microenvironment in hepatocellular carcinoma: Current status and therapeutic targets. Semin. Cancer Biol. 2011, 21, 35–43. [Google Scholar] [CrossRef] [Green Version]
  33. Walker, C.; Mojares, E.; Del Rio Hernandez, A. Role of Extracellular Matrix in Development and Cancer Progression. Int. J. Mol. Sci. 2018, 19, 3028. [Google Scholar] [CrossRef] [Green Version]
  34. Jiang, N.; Dai, Q.; Su, X.; Fu, J.; Feng, X.; Peng, J. Role of PI3K/AKT pathway in cancer: The framework of malignant behavior. Mol. Biol. Rep. 2020, 47, 4587–4629. [Google Scholar] [CrossRef] [PubMed]
  35. Niavarani, S.R.; Lawson, C.; Bakos, O.; Boudaud, M.; Batenchuk, C.; Rouleau, S.; Tai, L.H. Lipid accumulation impairs natural killer cell cytotoxicity and tumor control in the postoperative period. BMC Cancer 2019, 19, 823. [Google Scholar] [CrossRef]
  36. Goldberg, I.J.; Eckel, R.H.; Abumrad, N.A. Regulation of fatty acid uptake into tissues: Lipoprotein lipase- and CD36-mediated pathways. J. Lipid Res. 2009, 50, S86–S90. [Google Scholar] [CrossRef] [Green Version]
  37. Nath, A.; Li, I.; Roberts, L.R.; Chan, C. Elevated free fatty acid uptake via CD36 promotes epithelial-mesenchymal transition in hepatocellular carcinoma. Sci. Rep. 2015, 5, 14752. [Google Scholar] [CrossRef] [Green Version]
  38. Takeda, Y.; Kang, H.S.; Freudenberg, J.; DeGraff, L.M.; Jothi, R.; Jetten, A.M. Retinoic acid-related orphan receptor gamma (RORgamma): A novel participant in the diurnal regulation of hepatic gluconeogenesis and insulin sensitivity. PLoS Genet. 2014, 10, e1004331. [Google Scholar] [CrossRef]
  39. Massa, M.L.; Gagliardino, J.J.; Francini, F. Liver glucokinase: An overview on the regulatory mechanisms of its activity. IUBMB Life 2011, 63, 1–6. [Google Scholar] [CrossRef] [PubMed]
  40. Athwal, V.S.; Pritchett, J.; Llewellyn, J.; Martin, K.; Camacho, E.; Raza, S.M.; Phythian-Adams, A.; Birchall, L.J.; Mullan, A.F.; Su, K.; et al. SOX9 predicts progression toward cirrhosis in patients while its loss protects against liver fibrosis. EMBO Mol. Med. 2017, 9, 1696–1710. [Google Scholar] [CrossRef]
  41. Athwal, V.S.; Pritchett, J.; Martin, K.; Llewellyn, J.; Scott, J.; Harvey, E.; Zaitoun, A.M.; Mullan, A.F.; Zeef, L.A.H.; Friedman, S.L.; et al. SOX9 regulated matrix proteins are increased in patients serum and correlate with severity of liver fibrosis. Sci. Rep. 2018, 8, 17905. [Google Scholar] [CrossRef]
  42. Infantino, V.; Dituri, F.; Convertini, P.; Santarsiero, A.; Palmieri, F.; Todisco, S.; Mancarella, S.; Giannelli, G.; Iacobazzi, V. Epigenetic upregulation and functional role of the mitochondrial aspartate/glutamate carrier isoform 1 in hepatocellular carcinoma. Biochim. Biophys. Acta Mol. Basis Dis. 2019, 1865, 38–47. [Google Scholar] [CrossRef]
  43. Cruz, P.M.; Mo, H.; McConathy, W.J.; Sabnis, N.; Lacko, A.G. The role of cholesterol metabolism and cholesterol transport in carcinogenesis: A review of scientific findings, relevant to future cancer therapeutics. Front. Pharmacol. 2013, 4, 119. [Google Scholar] [CrossRef] [Green Version]
  44. Acimovic, J.; Rozman, D. Steroidal triterpenes of cholesterol synthesis. Molecules 2013, 18, 4002–4017. [Google Scholar] [CrossRef] [Green Version]
  45. Guy, J.; Peters, M.G. Liver disease in women: The influence of gender on epidemiology, natural history, and patient outcomes. Gastroenterol. Hepatol. (N. Y.) 2013, 9, 633–639. [Google Scholar]
  46. Kohi, M.P. Gender-Related Differences in Hepatocellular Carcinoma: Does Sex Matter? J. Vasc. Interv. Radiol. JVIR 2016, 27, 1338–1341. [Google Scholar] [CrossRef]
  47. Hartwell, H.J.; Petrosky, K.Y.; Fox, J.G.; Horseman, N.D.; Rogers, A.B. Prolactin prevents hepatocellular carcinoma by restricting innate immune activation of c-Myc in mice. Proc. Natl. Acad. Sci. USA 2014, 111, 11455–11460. [Google Scholar] [CrossRef] [Green Version]
  48. Keber, R.; Motaln, H.; Wagner, K.D.; Debeljak, N.; Rassoulzadegan, M.; Acimovic, J.; Rozman, D.; Horvat, S. Mouse knockout of the cholesterogenic cytochrome P450 lanosterol 14alpha-demethylase (Cyp51) resembles Antley-Bixler syndrome. J. Biol. Chem. 2011, 286, 29086–29097. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Gillespie, R.L.; Urquhart, J.; Anderson, B.; Williams, S.; Waller, S.; Ashworth, J.; Biswas, S.; Jones, S.; Stewart, F.; Lloyd, I.C.; et al. Next-generation Sequencing in the Diagnosis of Metabolic Disease Marked by Pediatric Cataract. Ophthalmology 2016, 123, 217–220. [Google Scholar] [CrossRef]
  50. Patel, N.; Anand, D.; Monies, D.; Maddirevula, S.; Khan, A.O.; Algoufi, T.; Alowain, M.; Faqeih, E.; Alshammari, M.; Qudair, A.; et al. Novel phenotypes and loci identified through clinical genomics approaches to pediatric cataract. Hum. Genet. 2017, 136, 205–225. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Rezen, T.; Ogris, I.; Sever, M.; Merzel, F.; Golic Grdadolnik, S.; Rozman, D. Evaluation of Selected CYP51A1 Polymorphisms in View of Interactions with Substrate and Redox Partner. Front. Pharmacol. 2017, 8, 417. [Google Scholar] [CrossRef] [Green Version]
  52. Lewinska, M.; Zelenko, U.; Merzel, F.; Golic Grdadolnik, S.; Murray, J.C.; Rozman, D. Polymorphisms of CYP51A1 from cholesterol synthesis: Associations with birth weight and maternal lipid levels and impact on CYP51 protein structure. PLoS ONE 2013, 8, e82554. [Google Scholar] [CrossRef] [Green Version]
  53. Michailidou, K.; Lindstrom, S.; Dennis, J.; Beesley, J.; Hui, S.; Kar, S.; Lemacon, A.; Soucy, P.; Glubb, D.; Rostamianfar, A.; et al. Association analysis identifies 65 new breast cancer risk loci. Nature 2017, 551, 92–94. [Google Scholar] [CrossRef] [Green Version]
  54. Yashin, A.I.; Wu, D.; Arbeeva, L.S.; Arbeev, K.G.; Kulminski, A.M.; Akushevich, I.; Kovtun, M.; Culminskaya, I.; Stallard, E.; Li, M.; et al. Genetics of aging, health, and survival: Dynamic regulation of human longevity related traits. Front. Genet. 2015, 6, 122. [Google Scholar] [CrossRef] [Green Version]
  55. Haines, T.H. Do sterols reduce proton and sodium leaks through lipid bilayers? Prog. Lipid Res. 2001, 40, 299–324. [Google Scholar] [CrossRef]
  56. Jiang, J.T.; Xu, N.; Zhang, X.Y.; Wu, C.P. Lipids changes in liver cancer. J. Zhejiang Univ. Sci. B 2007, 8, 398–409. [Google Scholar] [CrossRef]
  57. Hubler, Z.; Allimuthu, D.; Bederman, I.; Elitt, M.S.; Madhavan, M.; Allan, K.C.; Shick, H.E.; Garrison, E.; Molly, T.K.; Factor, D.C.; et al. Accumulation of 8,9-unsaturated sterols drives oligodendrocyte formation and remyelination. Nature 2018, 560, 372–376. [Google Scholar] [CrossRef]
  58. Gabitova, L.; Restifo, D.; Gorin, A.; Manocha, K.; Handorf, E.; Yang, D.H.; Cai, K.Q.; Klein-Szanto, A.J.; Cunningham, D.; Kratz, L.E.; et al. Endogenous Sterol Metabolites Regulate Growth of EGFR/KRAS-Dependent Tumors via LXR. Cell Rep. 2015, 12, 1927–1938. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Phillips, R.E.; Yang, Y.; Smith, R.C.; Thompson, B.M.; Yamasaki, T.; Soto-Feliciano, Y.M.; Funato, K.; Liang, Y.; Garcia-Bermudez, J.; Wang, X.; et al. Target identification reveals lanosterol synthase as a vulnerability in glioma. Proc. Natl. Acad. Sci. USA 2019. [Google Scholar] [CrossRef] [Green Version]
  60. Chen, L.; Ma, M.Y.; Sun, M.; Jiang, L.Y.; Zhao, X.T.; Fang, X.X.; Man Lam, S.; Shui, G.H.; Luo, J.; Shi, X.J.; et al. Endogenous sterol intermediates of the mevalonate pathway regulate HMGCR degradation and SREBP-2 processing. J. Lipid Res. 2019, 60, 1765–1775. [Google Scholar] [CrossRef]
  61. Tint, G.S.; Pan, L.; Shang, Q.; Sharpe, L.J.; Brown, A.J.; Li, M.; Yu, H. Desmosterol in brain is elevated because DHCR24 needs REST for Robust Expression but REST is poorly expressed. Dev. Neurosci. 2014, 36, 132–142. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Russell, D.W. Cholesterol biosynthesis and metabolism. Cardiovasc. Drugs Ther. 1992, 6, 103–110. [Google Scholar] [CrossRef]
  63. Lu, M.; Hu, X.H.; Li, Q.; Xiong, Y.; Hu, G.J.; Xu, J.J.; Zhao, X.N.; Wei, X.X.; Chang, C.C.; Liu, Y.K.; et al. A specific cholesterol metabolic pathway is established in a subset of HCCs for tumor growth. J. Mol. Cell Biol. 2013, 5, 404–415. [Google Scholar] [CrossRef] [Green Version]
  64. Oh, T.G.; Wang, S.M.; Acharya, B.R.; Goode, J.M.; Graham, J.D.; Clarke, C.L.; Yap, A.S.; Muscat, G.E.O. The Nuclear Receptor, RORgamma, Regulates Pathways Necessary for Breast Cancer Metastasis. EBioMedicine 2016, 6, 59–72. [Google Scholar] [CrossRef] [Green Version]
  65. Dooley, S.; ten Dijke, P. TGF-beta in progression of liver disease. Cell Tissue Res. 2012, 347, 245–256. [Google Scholar] [CrossRef] [Green Version]
  66. Pelikan, R.C.; Iwata, J.; Suzuki, A.; Chai, Y.; Hacia, J.G. Identification of candidate downstream targets of TGFbeta signaling during palate development by genome-wide transcript profiling. J. Cell. Biochem. 2013, 114, 796–807. [Google Scholar] [CrossRef] [Green Version]
  67. Matheu, A.; Collado, M.; Wise, C.; Manterola, L.; Cekaite, L.; Tye, A.J.; Canamero, M.; Bujanda, L.; Schedl, A.; Cheah, K.S.; et al. Oncogenicity of the developmental transcription factor Sox9. Cancer Res. 2012, 72, 1301–1315. [Google Scholar] [CrossRef] [Green Version]
  68. Vizcaino, C.; Mansilla, S.; Portugal, J. Sp1 transcription factor: A long-standing target in cancer chemotherapy. Pharmacol. Ther. 2015, 152, 111–124. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Vervoort, S.J.; Lourenco, A.R.; Tufegdzic Vidakovic, A.; Mocholi, E.; Sandoval, J.L.; Rueda, O.M.; Frederiks, C.; Pals, C.; Peeters, J.G.C.; Caldas, C.; et al. SOX4 can redirect TGF-beta-mediated SMAD3-transcriptional output in a context-dependent manner to promote tumorigenesis. Nucleic Acids Res. 2018, 46, 9578–9590. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  70. Caldwell, J.D.; Jirikowski, G.F. Sex hormone binding globulin and aging. Horm. Metab. Res. Horm. Stoffwechs. Horm. Metab. 2009, 41, 173–182. [Google Scholar] [CrossRef]
  71. Hautanen, A. Synthesis and regulation of sex hormone-binding globulin in obesity. Int. J. Obes. Relat. Metab. Disord. J. Int. Assoc. Study Obes. 2000, 24 (Suppl. 2), S64–S70. [Google Scholar] [CrossRef] [Green Version]
  72. Vaidya, D. Sex hormones and circulating lipoprotein levels. Future Lipidol. 2008, 3, 603–606. [Google Scholar] [CrossRef]
  73. Cichoz-Lach, H.; Michalak, A. Oxidative stress as a crucial factor in liver diseases. World J. Gastroenterol. 2014, 20, 8082–8091. [Google Scholar] [CrossRef]
  74. Wu, X.Z.; Chen, D. Origin of hepatocellular carcinoma: Role of stem cells. J. Gastroenterol. Hepatol. 2006, 21, 1093–1098. [Google Scholar] [CrossRef]
  75. Fang, C.H.; Gong, J.Q.; Zhang, W. Function of oval cells in hepatocellular carcinoma in rats. World J. Gastroenterol. 2004, 10, 2482–2487. [Google Scholar] [CrossRef]
  76. Acimovic, J.; Lovgren-Sandblom, A.; Monostory, K.; Rozman, D.; Golicnik, M.; Lutjohann, D.; Bjorkhem, I. Combined gas chromatographic/mass spectrometric analysis of cholesterol precursors and plant sterols in cultured cells. J. Chromatogr. B Anal. Technol. Biomed. Life Sci. 2009, 877, 2081–2086. [Google Scholar] [CrossRef]
  77. Skubic, C.; Vovk, I.; Rozman, D.; Krizman, M. Simplified LC-MS Method for Analysis of Sterols in Biological Samples. Molecules 2020, 25, 4116. [Google Scholar] [CrossRef] [PubMed]
  78. Andersen, C.L.; Jensen, J.L.; Orntoft, T.F. Normalization of real-time quantitative reverse transcription-PCR data: A model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets. Cancer Res. 2004, 64, 5245–5250. [Google Scholar] [CrossRef] [Green Version]
  79. Kosir, R.; Acimovic, J.; Golicnik, M.; Perse, M.; Majdic, G.; Fink, M.; Rozman, D. Determination of reference genes for circadian studies in different tissues and mouse strains. BMC Mol. Biol. 2010, 11, 60. [Google Scholar] [CrossRef] [Green Version]
  80. Carvalho, B.S.; Irizarry, R.A. A framework for oligonucleotide microarray preprocessing. Bioinformatics 2010, 26, 2363–2367. [Google Scholar] [CrossRef]
  81. Kauffmann, A.; Gentleman, R.; Huber, W. arrayQualityMetrics—A bioconductor package for quality assessment of microarray data. Bioinformatics 2009, 25, 415–416. [Google Scholar] [CrossRef] [Green Version]
  82. Smyth, G.K. Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat. Appl. Genet. Mol. Biol. 2004, 3. [Google Scholar] [CrossRef] [PubMed]
  83. Matys, V.; Kel-Margoulis, O.V.; Fricke, E.; Liebich, I.; Land, S.; Barre-Dirrie, A.; Reuter, I.; Chekmenev, D.; Krull, M.; Hornischer, K.; et al. TRANSFAC and its module TRANSCompel: Transcriptional gene regulation in eukaryotes. Nucleic Acids Res. 2006, 34, D108–D110. [Google Scholar] [CrossRef] [Green Version]
  84. Furge, K.; Dykema, K. PGSEA: Parametric Gene Set Enrichment Analysis, R package version 1.44.0 ed; The R Foundation: Vienna, Austria, 2006. [Google Scholar]
Figure 1. Blocking cholesterol synthesis due to depletion of Cyp51 resulted in histological alterations of the liver, most prominent in females. (A) Histological examination revealed ductular reaction (*) and inflammation (→) (H&E, original magnification ×200), N = 5–10. (B) Sirius red stain (SR, original magnification ×200) showed more pronounced fibrosis in female transgenic mice, N = 3–5. M, month; KO, Cyp51 KO; WT, wild-type.
Figure 1. Blocking cholesterol synthesis due to depletion of Cyp51 resulted in histological alterations of the liver, most prominent in females. (A) Histological examination revealed ductular reaction (*) and inflammation (→) (H&E, original magnification ×200), N = 5–10. (B) Sirius red stain (SR, original magnification ×200) showed more pronounced fibrosis in female transgenic mice, N = 3–5. M, month; KO, Cyp51 KO; WT, wild-type.
Cancers 12 03302 g001
Figure 2. (A) Macroscopic and histological presentation of liver tumors. (a) Representative gross image of liver tumors in 24M KO mice. (b) Representative histologic image of HCC showing architectural abnormalities. (c) Higher magnification of (c), showing cellular and nuclear polymorphism with brisk mitotic activity. Hematoxylin and Eosin (H&E), original magnification A, B ×100, C ×400. (B) Immunohistochemical expression of GS marker of HCC in livers of WTs and KOs. Original magnification, (a) ×200, (b) ×100, (c) ×400. GS, glutamine synthetase; T, tumor; PP, periportal; PC, pericentral, M, midzonal; WT, wild-type; KO, Cyp51 KO; T, tumor.
Figure 2. (A) Macroscopic and histological presentation of liver tumors. (a) Representative gross image of liver tumors in 24M KO mice. (b) Representative histologic image of HCC showing architectural abnormalities. (c) Higher magnification of (c), showing cellular and nuclear polymorphism with brisk mitotic activity. Hematoxylin and Eosin (H&E), original magnification A, B ×100, C ×400. (B) Immunohistochemical expression of GS marker of HCC in livers of WTs and KOs. Original magnification, (a) ×200, (b) ×100, (c) ×400. GS, glutamine synthetase; T, tumor; PP, periportal; PC, pericentral, M, midzonal; WT, wild-type; KO, Cyp51 KO; T, tumor.
Cancers 12 03302 g002
Figure 3. Adaptation to hepatocyte block in cholesterol synthesis during aging. Measurements were performed in KO and WT mice at 12M, 18M and 24M of age. Expression of hepatic genes from (A) from cholesterol synthesis (B) and bile acid synthesis. N = 5 mice per sex/genotype/age group. (C) Plasma concentration of biochemical parameters, N = 5–10 mice per sex/genotype/age group. CHOL—cholesterol, HDL—high-density lipoproteins, FFA—free fatty acids, ALT—alanine aminotransferase, AST—aspartate aminotransferase. (D) Liver concentration of sterol intermediates and cholesterol, N = 4–6 mice per sex/genotype/age group. Light color bands—95% confidence interval. Dots—individual measurements. KO, Cyp51 KO; WT, Cyp51 WT, F = female; M = male.
Figure 3. Adaptation to hepatocyte block in cholesterol synthesis during aging. Measurements were performed in KO and WT mice at 12M, 18M and 24M of age. Expression of hepatic genes from (A) from cholesterol synthesis (B) and bile acid synthesis. N = 5 mice per sex/genotype/age group. (C) Plasma concentration of biochemical parameters, N = 5–10 mice per sex/genotype/age group. CHOL—cholesterol, HDL—high-density lipoproteins, FFA—free fatty acids, ALT—alanine aminotransferase, AST—aspartate aminotransferase. (D) Liver concentration of sterol intermediates and cholesterol, N = 4–6 mice per sex/genotype/age group. Light color bands—95% confidence interval. Dots—individual measurements. KO, Cyp51 KO; WT, Cyp51 WT, F = female; M = male.
Cancers 12 03302 g003
Figure 4. Target genes associated with blocked cholesterol synthesis. (A) Heat-map of top 50 DE transcripts that showed significant differences in at least one comparison: hierarchical clustering distinguished between the genotypes and sex of 24M KO (24M.KO~WT) mouse livers. (B) Stacked column plots with numbers of differentially up- and down-regulated sex-specific genes of 24M Cyp51 KO mice. (C) The functional enriched networks of up- and down-regulated DE genes in female and male 24M KO mice using NetworkAnalyst. The sizes of the nodes correspond to the number of DE genes within a particular molecular process. All pathway nodules are differentially expressed (p < 0.05), the gradation rising from red to yellow (most enriched) in the network of positively (left) and negatively (right) enriched molecular processes. Circles representing cancer pathways and ECM-receptor interaction are surrounded by green lines.
Figure 4. Target genes associated with blocked cholesterol synthesis. (A) Heat-map of top 50 DE transcripts that showed significant differences in at least one comparison: hierarchical clustering distinguished between the genotypes and sex of 24M KO (24M.KO~WT) mouse livers. (B) Stacked column plots with numbers of differentially up- and down-regulated sex-specific genes of 24M Cyp51 KO mice. (C) The functional enriched networks of up- and down-regulated DE genes in female and male 24M KO mice using NetworkAnalyst. The sizes of the nodes correspond to the number of DE genes within a particular molecular process. All pathway nodules are differentially expressed (p < 0.05), the gradation rising from red to yellow (most enriched) in the network of positively (left) and negatively (right) enriched molecular processes. Circles representing cancer pathways and ECM-receptor interaction are surrounded by green lines.
Cancers 12 03302 g004
Figure 5. Hepatocyte-specific Cyp51 KO resulted in the deregulation of multiple signalling pathways and TFs, leading to the development of female-prevalent liver cancer. Arrows indicate connections between KEGG signalling pathways, TFs, and cellular processes. KEGG pathways and TFs are violet if positively enriched and green if negatively enriched at 24M. Top half—selected KEGG signalling pathways and TFs that were deregulated in females and males, leading to the inhibition of basic metabolism and an increase in processes of cell death and growth, and cellular communication. MAPK and PI3K/Akt signalling pathways have a central role and are regulated by other signalling pathways and TFs, all directly or indirectly contributing to carcinogenesis. Bottom half—female-specific KEGG signalling pathways and TFs. TF NFATC1/2 is central to the female-specific phenotype and is connected by the TGF-β, mTOR and Wnt signalling pathways, positively enriched only in female Cyp51 KO livers.
Figure 5. Hepatocyte-specific Cyp51 KO resulted in the deregulation of multiple signalling pathways and TFs, leading to the development of female-prevalent liver cancer. Arrows indicate connections between KEGG signalling pathways, TFs, and cellular processes. KEGG pathways and TFs are violet if positively enriched and green if negatively enriched at 24M. Top half—selected KEGG signalling pathways and TFs that were deregulated in females and males, leading to the inhibition of basic metabolism and an increase in processes of cell death and growth, and cellular communication. MAPK and PI3K/Akt signalling pathways have a central role and are regulated by other signalling pathways and TFs, all directly or indirectly contributing to carcinogenesis. Bottom half—female-specific KEGG signalling pathways and TFs. TF NFATC1/2 is central to the female-specific phenotype and is connected by the TGF-β, mTOR and Wnt signalling pathways, positively enriched only in female Cyp51 KO livers.
Cancers 12 03302 g005
Figure 6. Sex-dependent metabolic and transcriptional changes after disrupted Cyp51 from cholesterol synthesis align with hepatocarcinogenesis in humans. Arrows indicate active connections between enriched genes, KEGG pathways, and TFs. Blocked arrows represent repression. Hatched arrows represent possible transcription regulation of the selected pathway or connection between selected pathways in KO mice. Metabolites, DE genes, KEGG pathways, and TFs are violet if positively enriched and green if negatively enriched. Icons were used from the Reactome icon library. All selected mouse signalling and metabolic pathways and TFs were aligned with human literature data.
Figure 6. Sex-dependent metabolic and transcriptional changes after disrupted Cyp51 from cholesterol synthesis align with hepatocarcinogenesis in humans. Arrows indicate active connections between enriched genes, KEGG pathways, and TFs. Blocked arrows represent repression. Hatched arrows represent possible transcription regulation of the selected pathway or connection between selected pathways in KO mice. Metabolites, DE genes, KEGG pathways, and TFs are violet if positively enriched and green if negatively enriched. Icons were used from the Reactome icon library. All selected mouse signalling and metabolic pathways and TFs were aligned with human literature data.
Cancers 12 03302 g006
Table 1. The most significantly enriched KEGG pathways in 19W KO (19W.KO~WT) and 24M KO (24M.KO~WT) mice. Parameters shown in bold are statistically significant. KEGG, Kyoto Encyclopaedia of Genes and Genomes; * p < 0.05, ** p < 0.01, *** p < 0.001.
Table 1. The most significantly enriched KEGG pathways in 19W KO (19W.KO~WT) and 24M KO (24M.KO~WT) mice. Parameters shown in bold are statistically significant. KEGG, Kyoto Encyclopaedia of Genes and Genomes; * p < 0.05, ** p < 0.01, *** p < 0.001.
KEGG Pathway19W.KO~WT24M.KO~WT
FemalesMalesFemalesMales
LogFCadj.
p-Value
LogFCadj.
p-Value
LogFCadj.
p-Value
LogFCadj.
p-Value
Cell cycle1.130.01 *1.060.0241.60<0.001 ***2.13<0.001 ***
Pathways in cancer1.090.001 **0.790.019 *1.89<0.001 ***1.84<0.001 ***
ECM-receptor interaction1.300.005 **1.500.005 **1.82<0.001 ***1.85<0.001 ***
Proteoglycans in cancer1.100.001 **0.930.005 **1.63<0.001 ***1.61<0.001 ***
Apoptosis1.160.001 **0.930.003 **1.65<0.001 ***1.49<0.001 ***
p53 signalling pathway0.670.011 *0.600.032 *1.26<0.001 ***1.33<0.001 ***
MicroRNAs in cancer1.050.001 **0.680.026 *0.87<0.001 ***0.95<0.001 ***
Hepatocellular carcinoma0.670.001 **0.320.1310.95<0.001 ***0.720.006 **
Circadian entrainment0.210.1230.130.3870.420.020 *0.270.197
TGF-β pathway0.610.002 **0.380.0510.450.007 **0.140.427
Wnt signalling pathway0.290.0510.100.5880.320.018 *0.240.118
PI3K/Akt signalling pathway1.000.006 **0.940.016 *1.61<0.001 ***1.51<0.001 ***
Fatty acid degradation−1.24<0.001 ***−0.490.053−1.180.001 **−1.120.009 **
Metabolic pathways−1.630.001 **−0.140.818−1.800.009 **−0.690.367
Peroxisome−0.690.007 **−0.570.033 *−1.82<0.001 ***−1.200.018 *
Primary bile acid biosynthesis−1.76<0.001 ***−0.600.103−1.290.002 **−1.400.004 **
Biosynthesis of unsaturated fatty acid−0.470.063−0.050.874−0.870.018 *−0.460.261
Table 2. The most significantly enriched TFs in 19W KO (19W.KO~WT) and 24M KO (24M.KO~WT). Parameters shown in bold are statistically significant. * p < 0.05, ** p < 0.01, *** p < 0.001; n.s.—not significant. SP-1: Specificity protein 1; AP1: Activator protein 1; C/EBPβ: DNA damage-inducible transcript 3; PPARα: Peroxisome proliferator-activated receptor alpha; HIF1α: Hypoxia-inducible factor 1-alpha; NFATC1: Nuclear factor of activated T-cells, cytoplasmic 1; CLOCK: Circadian locomotor output cycles kaput; LXRα: Liver X receptor alpha; RXRα: Retinoid X receptor alpha; HNF4α: Hepatocyte nuclear factor 4 alpha; RORα: RAR-related orphan receptor alpha, RORC: RAR-related orphan receptor gamma; NR1B1: Retinoic acid receptor alpha.
Table 2. The most significantly enriched TFs in 19W KO (19W.KO~WT) and 24M KO (24M.KO~WT). Parameters shown in bold are statistically significant. * p < 0.05, ** p < 0.01, *** p < 0.001; n.s.—not significant. SP-1: Specificity protein 1; AP1: Activator protein 1; C/EBPβ: DNA damage-inducible transcript 3; PPARα: Peroxisome proliferator-activated receptor alpha; HIF1α: Hypoxia-inducible factor 1-alpha; NFATC1: Nuclear factor of activated T-cells, cytoplasmic 1; CLOCK: Circadian locomotor output cycles kaput; LXRα: Liver X receptor alpha; RXRα: Retinoid X receptor alpha; HNF4α: Hepatocyte nuclear factor 4 alpha; RORα: RAR-related orphan receptor alpha, RORC: RAR-related orphan receptor gamma; NR1B1: Retinoic acid receptor alpha.
Transcription Factor19W.KO~WT24M.KO~WT
FemalesMalesFemalesMales
LogFCadj.
p-Value
LogFCadj.
p-Value
LogFCadj.
p-Value
LogFCadj.
p-Value
SP10.810.002 **0.850.001 **1.47<0.001 ***1.49<0.001 ***
C-JUN0.930.001 **0.750.005 **1.42<0.001 ***1.50<0.001 ***
JUND0.920.002 **0.720.013 *1.05<0.001 ***0.91<0.001 ***
C-FOS1.20<0.001 ***0.98<0.001 ***1.47<0.001 ***1.33<0.001 ***
AP11.00<0.001 ***0.660.003 **1.17<0.001 ***0.88<0.001 ***
C/EBPβ0.520.010 *0.76<0.001 ***0.91<0.001 ***0.870.004 **
PPARα0.510.04 *0.060.830.600.001 **0.220.27
SOX90.340.0820.300.1250.90<0.001 ***0.630.001 **
(SOX9)20.460.012 *0.330.0570.670.005 **0.500.062
HIF1α0.330.013 *0.230.0700.430.008 **0.580.004 **
NFATC10.340.042 *0.180.2840.360.035 *0.290.156
CLOCK−0.240.035−0.170.122−0.260.031 *−0.410.009 **
LXRα:RXRαn.s.n.s.n.s.n.s.−0.560.043 *−0.190.602
HNF4α−0.730.007 **−0.270.290−1.60<0.001 ***−1.170.019 *
RORα−0.72<0.001 ***−0.200.256n.s.n.s.n.s.n.s.
RORC−0.350.028 *−0.110.474−1.68<0.001 ***−2.10<0.001 ***
NR1B1−0.210.087−0.350.007 **−0.80<0.001 ***−0.59<0.001 ***
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cokan, K.B.; Urlep, Ž.; Lorbek, G.; Matz-Soja, M.; Skubic, C.; Perše, M.; Jeruc, J.; Juvan, P.; Režen, T.; Rozman, D. Chronic Disruption of the Late Cholesterol Synthesis Leads to Female-Prevalent Liver Cancer. Cancers 2020, 12, 3302. https://doi.org/10.3390/cancers12113302

AMA Style

Cokan KB, Urlep Ž, Lorbek G, Matz-Soja M, Skubic C, Perše M, Jeruc J, Juvan P, Režen T, Rozman D. Chronic Disruption of the Late Cholesterol Synthesis Leads to Female-Prevalent Liver Cancer. Cancers. 2020; 12(11):3302. https://doi.org/10.3390/cancers12113302

Chicago/Turabian Style

Cokan, Kaja Blagotinšek, Žiga Urlep, Gregor Lorbek, Madlen Matz-Soja, Cene Skubic, Martina Perše, Jera Jeruc, Peter Juvan, Tadeja Režen, and Damjana Rozman. 2020. "Chronic Disruption of the Late Cholesterol Synthesis Leads to Female-Prevalent Liver Cancer" Cancers 12, no. 11: 3302. https://doi.org/10.3390/cancers12113302

APA Style

Cokan, K. B., Urlep, Ž., Lorbek, G., Matz-Soja, M., Skubic, C., Perše, M., Jeruc, J., Juvan, P., Režen, T., & Rozman, D. (2020). Chronic Disruption of the Late Cholesterol Synthesis Leads to Female-Prevalent Liver Cancer. Cancers, 12(11), 3302. https://doi.org/10.3390/cancers12113302

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop