Identification of novel lipid biomarkers in xmrk- and Myc-induced models of hepatocellular carcinoma in zebrafish
Cancer & Metabolism volume 10, Article number: 7 (2022)
Hepatocellular carcinoma (HCC) is the predominant form of liver cancer and is accompanied by complex dysregulation of lipids. Increasing evidence suggests that particular lipid species are associated with HCC progression. Here, we aimed to identify lipid biomarkers of HCC associated with the induction of two oncogenes, xmrk, a zebrafish homolog of the human epidermal growth factor receptor (EGFR), and Myc, a regulator of EGFR expression during HCC.
We induced HCC in transgenic xmrk, Myc, and xmrk/Myc zebrafish models. Liver specimens were histologically analyzed to characterize the HCC stage, Oil-Red-O stained to detect lipids, and liquid chromatography/mass spectrometry analyzed to assign and quantify lipid species. Quantitative real-time polymerase chain reaction was used to measure lipid metabolic gene expression in liver samples. Lipid species data was analyzed using univariate and multivariate logistic modeling to correlate lipid class levels with HCC progression.
We found that induction of xmrk, Myc and xmrk/Myc caused different stages of HCC. Lipid deposition and class levels generally increased during tumor progression, but triglyceride levels decreased. Myc appears to control early HCC stage lipid species levels in double transgenics, whereas xmrk may take over this role in later stages. Lipid metabolic gene expression can be regulated by either xmrk, Myc, or both oncogenes. Our computational models showed that variations in total levels of several lipid classes are associated with HCC progression.
These data indicate that xmrk and Myc can temporally regulate lipid species that may serve as effective biomarkers of HCC progression.
The adult liver can regulate lipid synthesis and degradation allowing it to play a major role in lipid metabolism . Dysregulation of lipid production occurs in human hepatocellular carcinoma (HCC) [2, 3], a form of liver cancer that is the third leading cause of cancer-related death in the world according to the World Health Organization (https://www.who.int/news-room/fact-sheets/detail/cancer). Because of the societal impact of HCC, there is considerable interest in the development of new models that can act as analytical platforms for the study of this disease. The zebrafish model is readily adaptable for studying human cancer as the zebrafish genome exhibits a high degree of sequence conservation with human oncogenes, exhibits similar tumor development and physiology, and is readily adaptable to genetic and chemical screening [4, 5]. Further, cancer cells expressed in zebrafish have very similar genetic and genomic characteristics relative to their human counterparts . Recently, transgenic adult zebrafish models have been successfully used to study oncogenes, immune physiology, gender-based effects, hormonal signaling, metastasis, drug efficacy, and tumor progression and regression in HCC [7,8,9,10,11,12,13].
We recently generated a doxycycline (DOX)-inducible liver tumor zebrafish model that transgenically expresses the myelocytomatosis (Myc) and the Xiphophorus melanoma receptor tyrosine kinase (xmrk activated epidermal growth factor receptor (EGFR) homolog) oncogenes [10, 11]. These models allowed us to control the temporality and location of oncogene induction via DOX exposure through regulation of the liver fatty acid–binding protein (fabp10a) promoter which is only expressed in hepatocytes . Here, we used the xmrk and Myc transgenic fish models to identify potential lipid biomarkers of HCC progression. Histological analysis demonstrated that Myc fish developed hyperplasia, xmrk fish induced HCC stage I, and some xmrk/Myc fish acquired the characteristics of HCC stage II. Single and double transgenic fish had less lipid accumulation in the early stages of tumor progression but began to exhibit markedly increased lipid deposition in the later stages. This increased lipid accumulation is similar to conditions observed during nonalcoholic fatty liver disease (NAFLD), a spectrum of disorders that comprise nonalcoholic steatohepatitis, cirrhosis, and HCC , although specific lipid levels can be either decreased or increased in NAFLD .
Comparative lipidomic analysis of the different transgenic models revealed that the levels of most lipid species in cancer cells increased except for triglycerides. Examination of the lipid species regulation profiles suggests that Myc may drive certain changes in lipid levels at earlier time points and that xmrk may take over this role at some later stage before the onset of mortality. Analyses of the relative expression of lipogenic transcription factors, lipogenic enzymes, and lipid β-oxidation genes indicate that xmrk and Myc may regulate distinct metabolic genes and may even counteract each other’s effects. Univariate and multifactorial modeling showed that several lipid classes increased during HCC progression, but that increased levels of other classes were correlated with restoration of the normal phenotype. Thus, xmrk and Myc signaling may act through specific lipids and lipid metabolic genes which may have utility as either diagnostic or prognostic biomarkers.
Zebrafish maintenance and transgene induction
Zebrafish embryos were raised as previously described . Adult zebrafish were euthanized by adding 1× tricaine (Sigma, St. Louis, USA) to the maintenance solution. Generation of TO(xmrk) (x+m−); TO(Myc) (x−m+) and TO(xmrk/Myc) (x+m+) transgenic models under the control of the lfabp promoter has been described previously [8, 11]. All zebrafish studies were approved by the Deakin University Animal Welfare Committee (AWC G17-2015) and the Institutional Animal Care and Use Committee of the National University of Singapore (Protocol 079/07). Transgenes were induced by exposing adult zebrafish (90 days old) to 40 μg/ml doxycycline (DOX) [10, 11].
Histology sample preparation and analysis
Liver samples taken from 1.5, 3, 4.5, 6, and 7.5 days post-DOX treatment (dpt) fish were slowly dehydrated using a series of increasing ethanol solution concentrations (70%, 90%, 95%, and 100%). Specimens were then embedded in paraffin using a Leica EG1120. Sectioning at 5 μm was performed using a Reichert-Jung 2030 microtome. Hematoxylin (H) (Vector Laboratories, Burlingame, CA, #H3404) and Eosin (E; Sigma, #HT110232) were used for H and E staining. Oil-Red-O (Sigma) staining of zebrafish liver sections was performed as previously described . Identification and classification of tumor types were based on previously established criteria .
To normalize samples, three to seven zebrafish livers were dissected and homogenized using a Pro200 homogenizer (Pro Scientific, Oxford, CT). Protein concentration was then quantified using a Pierce™ BCA protein assay kit (Life Technologies, Carlsbad, CA). An aliquot of each sample normalized to 30 μg was used for lipid extraction. Samples and standards used in this analysis were as previously described . Briefly, samples from 1.5, 3, 4.5, 6, and 7.5 dpt were lyophilized to remove all liquid. Samples were then reconstituted prior to extraction in 10 μl of water and 10 μl of the internal standard mix was added to each sample. Lipids were extracted by adding 200 μl chloroform/methanol (2:1) followed by sonication for 30 min before the supernatant was transferred to a 96-well plate and dried under vacuum in a SpeedVac Concentrator (Thermo Scientific, Waltham, MA). Samples were then reconstituted with 50 μl water-saturated butanol and 50 μl methanol with 10 mM ammonium formate, and analyzed by LC ESI-MS/MS, a technique that allows consistent measurement of various lipid species and classes in zebrafish over different experimental time points , using an Agilent 1200 LC system (Santa Clara, CA) and an AB Sciex 4000 Qtrap mass spectrometer (Darmstadt, Germany). Data were analyzed using MultiQuant 2.1 software (AB Sciex).
Quantitative real-time PCR analysis
Total RNA was extracted from 3, 4.5, and 6 dpt liver samples using the TRIzol® (Invitrogen #15596-018) and chloroform method per the manufacturer’s protocol. RNA was reverse-transcribed into cDNA using Transcriptor First Strand cDNA Synthesis Kit (Roche Applied Science, Penzberg, Germany). qRT-PCR was carried out using SYBR Green I Master Mix (Roche Applied Science) and a LightCycler® 480 machine (Roche Life Science). Cycle numbers (Ct) of triplicate samples were averaged and then normalized to β-actin expression to obtain a ΔCt value. Final gene expression levels were calculated as 2−ΔΔCt (see Supplemental Table 1 for a complete list of primers used in this study).
Lipid disease stage progression modeling analysis
As a final step, we constructed a model that would allow us to graphically depict the correlation between lipid levels and disease stage [21, 22]. 24 univariate ordinal logistic regression models relating each lipid to disease stage, i.e., normal, hyperplasia/adenoma, HCC1, and HCC2 were constructed using data from all time points (1.5, 3, 4.5, 6, 7.5 dpt). To facilitate ease of interpretation, each lipid was scaled by dividing the value by roughly ¼ of its standard deviation. Then, we put all 24 logistic regression models into a multivariate fractional polynomial (MFP) ordinal logistic model with backward selection variable procedures enabled. This approach allowed us to construct a concise model where only random error and no systematic bias is present in lipid measurement and where additive effects that might propagate are accounted for by removing lipids that contribute to multicollinearity and only those total lipid classes that display independence in each HCC stage are considered [20,21,22]. All variables remained linear after the MFP procedure, with the exception of GM3, which was given a square root transformation.
Lipidomic samples were statistically evaluated (GraphPAD PRISM, La Jolla, CA) using an ANOVA with Tukey’s post hoc test. qRT-PCR samples were statistically evaluated using an ANOVA. All statistical analyses for the logistic modeling analysis were performed with Stata v16.1 (College Station, TX). A probability level of p < 0.05 was deemed significant throughout.
Rapid oncogenic transformation of hepatocytes into HCC cells in xmrk/Myc transgenic zebrafish
Three DOX-inducible zebrafish transgenic oncogene lines, xmrk (x+m−), Myc (x−m+), and xmrk/Myc (x+m+) were evaluated for their ability to transform hepatocytes into HCC cells. After 1.5 days post-treatment (dpt), all the transgenic groups (N = 10 per group) showed similar histology characteristics compared to the non-transgenic DOX-exposed control group (x−m−) (Fig. 1), but at 3 dpt, single (x+m− and x−m+) and double transgenic zebrafish (x+m+) exhibited hyperplasia with increased cell number (Fig. 1B–D). Cellular transformation to HCC Grade I was first detected at 4.5 dpt in x+m− and x+m+ samples with evident increased apoptosis (Fig. 1B, D, arrows), and multiple nuclei in one cell (Fig. 1D, black box). By 7.5 dpt, two x+m+ individuals showed grade II HCC phenotype, and x+m+ fish started to die after 7.5 dpt; therefore, 7.5 dpt was chosen as the end point for analysis. Unlike the other samples, x−m+ livers remained hyperplastic throughout all time points (Fig. 1C).
Lipid class and species levels in zebrafish livers during tumor progression
Oil-red-O (ORO) staining was performed on liver samples to detect lipids. We found that x−m− samples had red staining lipid droplets present across all time points (1.5 to 7.5 dpt), while in transgenic samples, droplets were also detected and became particularly noticeable at 7.5 dpt (Fig. 2). Interestingly, lipid droplets are detected in greater number in control x−m− shortly after DOX induction (1.5 dpt to 3 dpt); however, after 6 and 7.5 dpt, lipid droplets are present in excess in all three transgene models compared to controls (Fig. 2). We analyzed liver sample lipid levels from 1.5 to 7.5 dpt using LC ESI-MS/MS. We identified many lipid classes whose levels exhibited temporal changes during these intervals, e.g., triglycerides (TG), dihexosylceramide (DHC), alkylphosphatidylcholine (PC(O)), monosialodihexosyl-ganglioside (GM3), dihydroceramides (dhCer), ceramides (Cer), monohexosylceramide (MHC), trihexosylceramide (THC), sphingomyelin (SM), lyso-platelet-activating factor (LPAF), lysophosphatidylcholine (LPC), phosphatidylcholine (PC), alkenylphosphatidyl-choline (PC(P)), phosphatidylethanolamine (PE), alkylphosphatidylethanolamine (PE(O)), alkenylphosphatidylethanolamine (PE(P)), lysophosphatidylethanolamine (LPE), phosphatidylinositol (PI), lysophosphatidylinositol (LPI), phosphatidylserine (PS), phosphatidylglycerol (PG), cholesteryl ester (CE), free cholesterol (COH), and diacylglycerol (DG) (Fig. 3, supplementary Fig. 1). DHC and PC(O) levels in x+m− fish were higher at 7.5 dpt (Fig. 3A). GM3 lipid levels in x+m+ samples were increased at 6 dpt (Fig. 3A). PI lipid levels in x+m− and x−m+ samples were significantly higher at 7.5 dpt (Fig. 3A). At 6 dpt, x+m+ livers expressed higher levels of PE (Fig. 3B), and at 6 and 7.5 dpt, x+m− samples had higher levels of PS (Fig. 3B). TG levels exhibited a downward trend across all time points (Fig. 3B). Heat-map analysis revealed that relative lipid class abundance in the transgenic tissues compared to control generally increased between 1.5 to 7.5 dpt (Fig. 4).
The LC ESI-MS/MS analysis also showed that inducing xmrk and Myc caused altered levels of numerous lipid species (Table 1, supplementary Fig. 2). In x−m+ samples, we found the following levels at various time points: 1.5 dpt, several decreased PEs, PE 20.4-0:0 increased; 3 dpt, PC(O)s, PEs, and PI 20:4-0-0 increased; 4.5 dpt PC(O) 35:4 was increased; 6 dpt, PC(O) 20:0-0:0, PEs, and two PIs increased; 7.5 dpt, several TGs were decreased, while several PC(O)s, PEs, PIs, and PSs were all increased (see Table 1 for details). In x+m− samples, we found that xmrk induced regulation of different lipid species when compared to Myc transgenic samples. At 1.5 dpt, x+m− fish did not have decreased lipid species, but x+m− fish had increased levels of DHC 20:0, GM3 16:0, and multiple PC(O)s, PEs, PIs and TGs (Table 1). At 3 dpt, x+m− fish livers showed decreased levels of lipid species, e.g., the TGs, while only one PC(O), 35:4, was increased in x+m− samples (Table 1). Xmrk only transgenic samples at 4.5 dpt also showed increased PC(O) 35:4 and many other lipid species, e.g., several other PC(O)s, DHC 22:0, several PEs, PIs, and two PS species (Table 1). At 6 dpt, x+m− fish also showed increased PC(O) 20:0-0:0, but x+m− fish showed increased levels of DHCs, GM3 24:1, additional PC(O)s, PEs, PIs, and several PSs (Table 1). At 7.5 dpt, x+m− fish samples did not show decreased TGs, but did show increased levels of PC(O)s, PEs, PIs, and PSs, and xmrk induction alone increased DHCs (Table 1).
The double transgenics also had unique lipid species characteristics. At 1.5 dpt, we found that x+m+ fish also had lipid species decreased in x−m+ samples and PC(O)s, additional PEs, PI 38:4 and PS 38.4, while unlike x+m−, x+m+ had increased PI 34.1 and decreased PC(O) 36:4 and PE 36:3 (Table 1). Double transgenics at 3 dpt exhibited similar lipid species regulation to x−m+ samples; however, x+m+ fish also had several TGs increased (Table 1). At 4.5 dpt, x+m+ lipid species levels in liver tissues resembled those in x+m− samples, but not to the same extent, and expressed increased GM3 (22:0 and 24:1) (Table 1). At 6 dpt, x+m+ liver profiles were similar to those observed in x+m−, but x+m+ samples had increased levels of several TGs (Table 1). Data from 7.5 dpt time points showed that double transgenics exhibited similar PC(O), PE, and PI lipid species upregulation as observed in the single transgenics (Table 1). However, x+m+ fish did not have significant levels of PS species, but only they expressed GM3 species (Table 1).
Expression of lipid metabolism genes in liver tumors
We next measured the expression of key metabolic genes in DOX-treated liver tissue specimens at selected stages (3, 4.5, and 6 dpt) corresponding to the rapid period of transformation into hyperplasia and HCC (Fig. 1B–D). These genes included: lipogenic factors (peroxisome proliferator-activated receptor gamma [pparg], sterol regulatory element-binding transcription factor 1 [srebf1], CCAAT/enhancer-binding protein alpha [cebpa]), lipogenic enzymes (fatty acid synthase [fasn], phosphatidic acid phosphatase [pap], diacylglycerol O-acyltransferase 2 [dgat2]), and lipid β-oxidation genes (peroxisome proliferator-activated receptor alpha b [pparab], carnitine palmitoyltransferase 1 [cpt1], cytochrome P450 monooxygenase [cyp4a10], acyl-coenzyme A oxidase 3 [acox3] )[23,24,25,26,27,28,29,30,31,32]. The majority of these genes exhibited downregulation at 3 dpt compared to the controls, with the exception of cpt1, which was upregulated in x−m+ samples (Fig. 5A, B). At 4.5 dpt, in x−m+ samples, cebpa, dgat2, fasn, and srebf1 were upregulated and pparb was downregulated, in x+m− samples, dgat2 and srebf1 were downregulated, while acox3, cebpa, cyp4a10, fasn, pap, and pparab showed increased expression, and in x+m+, cyp4a10, dgat2, pparab, and srebf1 expression was decreased, while cebpa and fasn were upregulated (Fig. 5A, B). At 6 dpt, x−m+ transgenics showed upregulated acox3, cpt1, cyp4a10, fasn, pparg, and srebf1 expression, and decreased pap expression; x+m− showed downregulated dgat2, pap, and pparg expression, and increased acox3, cebpa, fasn, and pparab expression, while x+m+ samples exhibited upregulation of acox3, cebpa, cpt1, fasn, and srebf1, and downregulation of cyp4a10 and dgat2 expression with pap and pparg expression being unchanged (Fig. 5A, B).
Association of lipid class regulation with HCC disease progression
To identify lipid classes that may contribute to HCC disease progression, we related each lipid class to HCC disease stage using univariate ordinal logistic regression models. Univariate analysis showed that total DHC, GM3, PC(O), and PI levels were correlated with increased promotion of HCC II (Fig. 6A), while increasing total PS, PE, and TG levels were associated with the normal, non-HCC phenotype (Fig. 6B, see supplementary Figs. 3, 4, 5 and 6 for additional univariate data sets). However, since univariate analyses of individual lipid classes may not accurately reflect their potential systemic role in HCC progression, we also used multivariate fractional polynomial (MFP) analysis to interrelate the potential effect of each lipid class on HCC development (Table 2). The MFP analysis confirmed that DHC, PC(O), GM(3), PS, PE, PI, and TG levels were associated with the progression of HCC during oncogenic induction (Table 2, see supplementary Fig. 1 for excluded lipid class data sets). Specifically, four lipid sub-classes displayed increased odds of disease severity: DHC [OR = 2.17; (1.30–3.61); p = 0.003], GM3 [OR = 1.64; (1.28–2.11); p < 0.001], PC(O) [OR = 4.32; (1.87–10.0); p = 0.001], and PI [OR = 3.22; (1.85–5.62); p < 0.001] had increased odds of disease severity per 40, 10 square root, 630, and 1530 units, respectively (Table 2). The opposite relationship was observed for 3 lipid sub-classes: PE [OR = 0.52; (0.29–0.95); p = 0.033], PS [OR = 0.29; (0.14–0.61); p = 0.001], and TG [OR = 0.73; (0.60–0.89); p = 0.001] that had decreased odds of disease severity per 2140, 350, and 1810 units, respectively (Table 2). The MFP analysis produced scaling factors of 40 (DHC), 10 (GM3), 630 (PC(O)), 2140 (PE), 1530 (PI), 350 (PS), and 1810 (TG) units, respectively. Also, the final model had the most superior Akaike information criterion  values (88.70) and R2 (0.55) among examined models.
The results of the univariate ordinal logistic regression models were placed into a multivariate fractional polynomial analysis with backward selection variable procedures enabled to determine which lipid classes significantly modulated the HCC disease stage. All results were linear except that total GM3 used a square root nonlinear transformation (AIC =88.7 and R2 = 0.55). OR, odds ratio; p < 0.05.
Lipids and their metabolites can play roles as energy, signaling, and/or prognostic biomarkers in HCC [34, 35], but oncogene-mediated regulation of specific lipids and related metabolomic genes during liver cancer progression have not been well studied. Therefore, we utilized transgenic zebrafish models of HCC expressing the liver cancer–inducing oncogenes, xmrk and Myc [8,9,10,11] to analyze lipidomic profiles and lipid metabolic genes at different HCC stages. Xmrk is a fish homolog of the human epidermal growth factor receptor (EGFR), a transmembrane receptor tyrosine kinase that regulates multiple signaling pathways controlling cell proliferation, migration, and inhibition . EGFR plays a complex role in liver cancer where it is typically upregulated . Some forms of Myc act as transcription factors that can modulate EGFR gene expression , and Myc amplification is correlated with HCC development and lipogenesis [39, 40]. Further, EGFR stimulation can upregulate Myc expression in HCC . We found that x+m−, x−m+, and x+m+ liver samples exhibited characteristics of HCC oncogenesis  within 7.5 dpt (Fig. 1) indicating that these oncogenes modulate phenotypic attributes of HCC. These data suggest that combined transgene induction may enhance disease progression perhaps via xmrk-mediated self-potentiation by a signaling loop integrating Myc.
The presence of increased lipid levels at later time points in the single and double transgenics (Fig. 2) suggests that xmrk and Myc might promote HCC oncogenesis via specific types of lipids, but surprisingly, few of the lipid classes that we detected were altered (Fig. 3, supplementary Fig. 1). During HCC, lipid metabolites and lipid-related gene regulation can exhibit extensive up- and downregulation , and we observed altered regulation of some lipid classes, e.g., DHC, PC(O), GM3, PI, PE, and PS (Fig. 3A, B), where a general upward trend in levels was found beginning at 4.5 or 6 dpt, and TGs, where there is a general downward trend across all time points (Fig. 3B). The lipid classes PE and PI also can be dysregulated and have been correlated with HCC disease and its progression [43, 44]. Further, the overall increasing relative levels of lipid classes across the time points (Fig. 4) may suggest ongoing lipid metabolic reprogramming which is observed in human HCC as well, where it has been associated with enhanced lipid synthesis . In summary, altered lipid levels could represent a response to the different metabolic and/or signaling demands encountered during HCC progression.
Analysis of specific lipid species showed that induction of xmrk and Myc alone in HCC may regulate largely distinct sets of lipid species suggesting that each transgene may target lipid-based metabolism and signaling mechanisms (Table 1). In human HCC cell lines, the expression of EGFR and other genes has been correlated with lipid species regulation . Further, the Myc oncogene can modulate some lipid species during HCC, particularly members of the PG class . Recently, EGFR and Myc signals have been shown to be integrated into a signaling axis that modulates oncogenesis . Therefore, expressing both oncogenes simultaneously may result in the utilization of different lipid species when compared to the single transgene models. Interestingly, when we compared the lipid species associated with the single transgenics with lipid species associated with the double transgenics at the same time points, we found that at 1.5 and 3 dpt, Myc appeared to drive lipid species levels in the double transgenics, but that from 4.5 to 6 dpt, xmrk may supplant this role in x+m+ HCC liver samples (Table 1). These results suggest that shifts in the metabolic or signaling demands associated with particular HCC disease stages may be responsible for the transposition of xmrk and Myc regulation of lipid species in x+m+. However, neither oncogene has a primary role in modulating lipid species regulation in x+m+ samples at 7.5 dpt, when mortality begins to appear (Table 1). Lipid-based mechanisms can regulate cell death mechanisms in cancer , and some lipid classes, e.g., GM3, which exhibits increased levels in late HCC stage x+m+ samples (Fig. 3A), can inhibit EGFR function  suggesting that the levels of some lipids in the transgenic models could be associated with late HCC stage cell death.
Distinct xmrk and Myc temporal regulation of lipids could mean that these oncogenes modulate different lipid metabolic genes. Measuring the expression profiles of these genes over various time points allows assessment of their status as biomarkers without necessitating the analysis of the corresponding enzyme profiles or their mechanistic role in regulating lipid levels. The expression profiles for the lipogenic factors, pparg and srebf1, and for the enzyme, dgat2, in x+m+ samples are generally similar to those observed in xmrk only model (Fig. 5A), suggesting that xmrk may control these genes; although, Myc may be counteracting xmrk regulation of srebf1 by 6 dpt (Fig. 5A). Activation of pparg causes growth inhibition and cell death in HCC cells  suggesting that xmrk-mediated reduction of pparg expression (Fig. 5A) may promote oncogenesis in HCC. Increased expression of srebf1 is associated with cell proliferation in liver cancer  indicating that Myc may contribute to HCC oncogenesis at 6 dpt (Fig. 5A). As suppression of the lipogenic enzyme gene, dgat2, is associated with increased HCC cell proliferation , xmrk induction also could be driving disease progression by targeting dgat2 function (Fig. 5B). Expression of the lipid β-oxidation genes exhibited more complex regulation with pparab correlated with Myc, cpt1 with both xmrk and Myc, and cyp4a10 with xmrk (Fig. 5C). As reduced expression of pparab is associated with increased HCC development , Myc could be promoting liver cancer onset by suppressing lipid breakdown via pparab function. Further, as increased cpt1 is associated with HCC progression , these data support the suggestion that xmrk and Myc induction together promote the development of liver cancer via enhanced cpt1 expression. Interestingly, cyp4a10 expression is increased in liver tumor development in mice , but we found that cyp4a10 expression in the x+m+ samples was reduced suggesting that xmrk may counteract HCC progression through this gene (Fig. 5C). Therefore, xmrk and Myc may modulate distinct subsets of lipid metabolic genes that may have some utility as biomarkers of HCC progression.
In support of this concept, altered lipid levels may act as an indicator of HCC disease progression as lipid levels change during HCC [42, 48]. DHC is associated with the regulation of differentiation, proliferation, and programmed cell death , and its elevated levels (Fig. 6A, Table 2) may be an indicator of oncogenic progression. GM3 can suppress cancer progression ; therefore, its elevation during HCC progression (Fig. 6A, Table 2) could be a biomarker of a cellular response to counteract oncogenesis. PC(O) is a phospholipid, and some phospholipid metabolites can mediate proliferative growth and programmed cell death in cancer ; therefore, PC(O) elevation (Fig. 6A, Table 2) may be a marker of advanced HCC disease stage. PI species can participate in cell signaling in cancer  suggesting that increased PI levels during HCC (Fig. 6A, Table 2) could indicate either pro- or antioncogenic responses. PE phospholipids exhibit increased levels during cancer where they play a role in cell division and death . Increased PE (Fig. 6B, Table 2) during HCC progression could represent activation of a cell death mechanism that counteracts liver cancer. PS phospholipid levels may regulate phagocytosis, immunosuppression, tumor growth, and metastasis . Therefore, increased PS levels (Fig. 6B, Table 2) could indicate phagocytic anticancer activity. TG levels in the liver are regulated by uptake, secretory, and metabolic mechanisms , and higher TG levels (Fig. 6B, Table 2) could represent reduced mobilization of TG for driving mechanisms that promote HCC oncogenesis. Thus, altered levels of some lipid classes and their constituent species may act as diagnostic biomarkers of HCC disease progression.
When performing liver lipid analysis in HCC, it is important to bear in mind that changes in the diet can affect liver lipid content . For this reason, our control and DOX-exposed fish were fed the same diet and amount to rule out inducing any effect from qualitative or quantitative variations in food intake. Moreover, in rats, dietary change can affect liver lipid content in non-alcoholic fatty liver disease (NAFLD) models, but the main lipid classes that displayed a strong variability in the liver were TGs and the lysophospholipids, LPC and LPE, which were higher in high-fat diet (HFD) rat livers . This observation is consistent with our data as neither LPC nor LPE were identified as biomarkers of cancer progression (Figure 6). Further, unlike in many cases of HFD, where TGs levels are dramatically increased in the liver [57,58,59,60], in the data we derived, decreased TG levels indicate a non-cancerous liver (Figure 6B). Therefore, the lipid species identified in this study are not responding to diet variations and, instead, are hepatocyte cancer state specific. As zebrafish are lecithotrophic organisms, i.e., the embryo receives no nutrition other than what the yolk sac contains originally, a difference in the lipid content deposited in the yolk sac by the mother during oocyte formation might affect the future lipid content of the liver in adults. However, this effect is unlikely to have occurred in this study, because reproductive females were from similar batches of fish and were all fed the same quality and quantity of food, so a variation in yolk sac lipid content should not be present or responsible for the lipid variation observed in the HCC fish models.
These studies indicate that induction of xmrk, Myc, and xmrk/Myc causes temporal alteration of lipid class levels and metabolic gene expression at different stages of HCC. Increased lipid class levels are associated with HCC progression or restoration of the normal phenotype and may act as predictive biomarkers of HCC.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Santos CR, Schulze A. Lipid metabolism in cancer. FEBS J. 2012;279:2610–23.
Bian X, Liu R, Meng Y, Xing D, Xu D, Lu Z. Lipid metabolism and cancer. J Exp Med. 2020;218:e20201606.
Hu B, Lin JZ, Yang XB, Sang XT. Aberrant lipid metabolism in hepatocellular carcinoma cells as well as immune microenvironment: a review. Cell Prolif. 2020;53:e12772.
Hason M, Bartůněk P. Zebrafish models of cancer-new insights on modeling human cancer in a non-mammalian vertebrate. Genes (Basel). 2019;10:935.
Gamble JT, Elson DJ, Greenwood JA, Tanguay RL, Kolluri SK. The zebrafish xenograft models for investigating cancer and cancer therapeutics. Biology (Basel). 2021;10:252.
Lam SH, Wu YL, Vega VB, Miller LD, Spitsbergen J, Tong Y, et al. Conservation of gene expression signatures between zebrafish and human liver tumors and tumor progression. Nat Biotechnol. 2006;24:73–5.
Lu JW, Hsia Y, Tu HC, Hsiao YC, Yang WY, Wang HD, et al. Liver development and cancer formation in zebrafish. Birth Defects Res C Embryo Today. 2011;93:157–72.
Li Z, Huang X, Zhan H, Zeng Z, Li C, Spitsbergen JM, et al. Inducible and repressable oncogene-addicted hepatocellular carcinoma in Tet-on xmrk transgenic zebrafish. J Hepatol. 2012;56:419–25.
Li Z, Luo H, Li C, Huo X, Yan C, Huang X, et al. Transcriptomic analysis of a transgenic zebrafish hepatocellular carcinoma model reveals a prominent role of immune responses in tumour progression and regression. Int J Cancer. 2014;135:1564–73.
Zheng W, Li Z, Nguyen AT, Li C, Emelyanov A, Gong Z. Xmrk, kras and myc transgenic zebrafish liver cancer models share molecular signatures with subsets of human hepatocellular carcinoma. PLoS One. 2014;9:e91179.
Li Z, Zheng W, Li H, Li C, Gong Z. Synergistic induction of potential Warburg effect in zebrafish hepatocellular carcinoma by co-transgenic expression of Myc and xmrk oncogenes. PLoS One. 2015;10:e0132319.
Li H, Lu JW, Huo X, Li Y, Li Z, Gong Z. Effects of sex hormones on liver tumor progression and regression in Myc/xmrk double oncogene transgenic zebrafish. Gen Comp Endocrinol. 2019;277:112–21.
Nakayama J, Makinoshima H. Zebrafish-based screening models for the identification of anti-metastatic drugs. Molecules. 2020;25:2407.
Zhang X, Li C, Gong Z. Development of a convenient in vivo hepatotoxin assay using a transgenic zebrafish line with liver-specific DsRed expression. PLoS One. 2014;9:e91874.
Pei K, Gui T, Kan D, Feng H, Jin Y, Yang Y, et al. An overview of lipid metabolism and nonalcoholic fatty liver disease. Biomed Res Int. 2020;2020:4020249.
Puri P, Baillie RA, Wiest MM, Mirshahi F, Choudhury J, Cheung O, et al. A lipidomic analysis of nonalcoholic fatty liver disease. Hepatology. 2007;46:1081–90.
Westerfield M. The Zebrafish Book. A guide for the laboratory use of zebrafish (Danio rerio). 5th Edition. Eugene: University of Oregon Press; 2007.
Li C, Li P, Tan YM, Lam SH, Chan EC, Gong Z. Metabolomic characterizations of liver injury caused by acute arsenic toxicity in zebrafish. PLoS One. 2016;11:e0151225.
Weir JM, Wong G, Barlow CK, Greeve MA, Kowalczyk A, Almasy L, et al. Plasma lipid profiling in a large population-based cohort. J Lipid Res. 2013;54:2898–908.
Fraher D, Sanigorski A, Mellett NA, Meikle PJ, Sinclair AJ, Gibert Y. Zebrafish embryonic lipidomic analysis reveals that the yolk cell is metabolically active in processing lipid. Cell Rep. 2016;14:1317–29.
Sauerbrei W, Royston P. Building multivariable prognostic and diagnostic models: transformation of the predictors by using fractional polynomials. J R Statist Soc A. 1999;162:71–94.
Zhang Z. Multivariable fractional polynomial method for regression model. Ann Transl Med. 2016;4:174.
Hsu HT, Chi CW. Emerging role of the peroxisome proliferator-activated receptor-gamma in hepatocellular carcinoma. J Hepatocell Carcinoma. 2014;1:127–35.
Lu GD, Ang YH, Zhou J, Tamilarasi J, Yan B, Lim YC, et al. CCAAT/enhancer binding protein α predicts poorer prognosis and prevents energy starvation-induced cell death in hepatocellular carcinoma. Hepatology. 2015;61:965–78.
Qu Q, Zeng F, Liu X, Wang QJ, Deng F. Fatty acid oxidation and carnitine palmitoyltransferase I: emerging therapeutic targets in cancer. Cell Death Dis. 2016;7:e2226.
Kuwata K, Inoue K, Ichimura R, Takahashi M, Kodama Y, Yoshida M. Constitutive active/androstane receptor, peroxisome proliferator-activated receptor α, and cytotoxicity are involved in oxadiazon-induced liver tumor development in mice. Food Chem Toxicol. 2016;88:75–86.
Xiao YB, Cai SH, Liu LL, Yang X, Yun JP. Decreased expression of peroxisome proliferator-activated receptor alpha indicates unfavorable outcomes in hepatocellular carcinoma. Cancer Manag Res. 2018;10:1781–9.
Zhao X, Zhao L, Yang H, Li J, Min X, Yang F, et al. Pyruvate kinase M2 interacts with nuclear sterol regulatory element-binding protein 1a and thereby activates lipogenesis and cell proliferation in hepatocellular carcinoma. J Biol Chem. 2018;293:6623–34.
Che L, Paliogiannis P, Cigliano A, Pilo MG, Chen X, Calvisi DF. Pathogenetic, prognostic, and therapeutic role of fatty acid synthase in human hepatocellular carcinoma. Front Oncol. 2019;9:1412.
Li Y, Li T, Jin Y, Shen J. Dgat2 reduces hepatocellular carcinoma malignancy via downregulation of cell cycle-related gene expression. Biomed Pharmacother. 2019;115:108950.
Xu A, Wang B, Fu J, Qin W, Yu T, Yang Z, et al. Diet-induced hepatic steatosis activates Ras to promote hepatocarcinogenesis via CPT1α. Cancer Lett. 2019;442:40–52.
Brohée L, Crémer J, Colige A, Deroanne C. Lipin-1, a versatile regulator of lipid homeostasis, is a potential target for fighting cancer. Int J Mol Sci. 2021;22:4419.
Portet S. A primer on model selection using the Akaike Information Criterion. Infect Dis Model. 2020;5:111–28.
Hayes CN, Zhang P, Chayama K. The role of lipids in hepatocellular carcinoma. In: Tirnitz-Parker JEE, editor. Hepatocellular carcinoma. Brisbane: Codon Publications; 2019. p. 95–110.
Sasaki N, Toyoda M, Ishiwata T. Gangliosides as signaling regulators in cancer. Int J Mol Sci. 2021;22:5076.
Monroe JD, Basheer F, Gibert Y. Xmrks the spot: fish models for investigating epidermal growth factor receptor signaling in cancer research. Cells. 2021;10:1132.
Komposch K, Sibilia M. EGFR signaling in liver diseases. Int J Mol Sci. 2015;17:30.
Perini G, Diolaiti D, Porro A, Della VG. In vivo transcriptional regulation of N-Myc target genes is controlled by E-box methylation. Proc Natl Acad Sci U S A. 2005;102:12117–22.
Takahashi Y, Kawate S, Watanabe M, Fukushima J, Mori S, Fukusato T. Amplification of c-myc and cyclin D1 genes in primary and metastatic carcinomas of the liver. Pathol Int. 2007;57:437–42.
Gouw AM, Margulis K, Liu NS, Raman SJ, Mancuso A, Toal GG, et al. The MYC oncogene cooperates with sterol-regulated element-binding protein to regulate lipogenesis essential for neoplastic growth. Cell Metab. 2019;30:556–572.e5.
Pan Z, Liu C, Zhi Y, Xie Z, Wu L, Jiang M, et al. LIMK1 nuclear translocation promotes hepatocellular carcinoma progression by increasing p-ERK nuclear shuttling and by activating c-Myc signalling upon EGF stimulation. Oncogene. 2021;40:2581–95.
Budhu A, Roessler S, Zhao X, Yu Z, Forgues M, Ji J, et al. Integrated metabolite and gene expression profiles identify lipid biomarkers associated with progression of hepatocellular carcinoma and patient outcomes. Gastroenterology. 2013;144:1066–75.
Li Z, Guan M, Lin Y, Cui X, Zhang Y, Zhao Z, et al. Aberrant lipid metabolism in hepatocellular carcinoma revealed by liver lipidomics. Int J Mol Sci. 2017;18:2550.
Ismail IT, Elfert A, Helal M, Salama I, El-Said H, Fiehn O. Remodeling lipids in the transition from chronic liver disease to hepatocellular carcinoma. Cancers (Basel). 2020;13:88.
Nakagawa H, Hayata Y, Kawamura S, Yamada T, Fujiwara N, Koike K. Lipid metabolic reprogramming in hepatocellular carcinoma. Cancers (Basel). 2018;10:447.
Rivas Serna IM, Romito I, Maugeri A, Lo Re O, Giallongo S, Mazzoccoli G, et al. A lipidomic signature complements stemness features acquisition in liver cancer cells. Int J Mol Sci. 2020;21:8452.
Huang C, Freter C. Lipid metabolism, apoptosis and cancer therapy. Int J Mol Sci. 2015;16:924–49.
Sangineto M, Villani R, Cavallone F, Romano A, Loizzi D, Serviddio G. Lipid metabolism in development and progression of hepatocellular carcinoma. Cancers (Basel). 2020;12:1419.
Ogretmen B. Sphingolipid metabolism in cancer signalling and therapy. Nat Rev Cancer. 2018;18:33–50.
Zheng C, Terreni M, Sollogoub M, Zhang Y. Ganglioside GM3 and its role in cancer. Curr Med Chem. 2019;26:2933–47.
Ridgway ND. The role of phosphatidylcholine and choline metabolites to cell proliferation and survival. Crit Rev Biochem Mol Biol. 2013;48:20–38.
Beloribi-Djefaflia S, Vasseur S, Guillaumond F. Lipid metabolic reprogramming in cancer cells. Oncoscience. 2014;2:15–30.
Tan LT, Chan KG, Pusparajah P, Lee WL, Chuah LH, Khan TM, et al. Targeting membrane lipid a potential cancer cure? Front Pharmacol. 2017;8:12.
Chang W, Fa H, Xiao D, Wang J. Targeting phosphatidylserine for cancer therapy: prospects and challenges. Theranostics. 2020;10:9214–29.
Alves-Bezerra M, Cohen DE. Triglyceride metabolism in the liver. Compr Physiol. 2017;8:1–8.
Hodson L, Rosqvist F, Parry SA. The influence of dietary fatty acids on liver fat content and metabolism. Proc Nutr Soc. 2020;79:30–41.
Yaligar J, Gopalan V, Kiat OW, Sugii S, Shui G, Lam BD, et al. Evaluation of dietary effects on hepatic lipids in high fat and placebo diet fed rats by in vivo MRS and LC-MS techniques. PLoS One. 2014;9:e91436.
Handayani D, Meyer BJ, Chen J, Brown SH, Mitchell TW, Huang XF. A high-dose Shiitake mushroom increases hepatic accumulation of triacylglycerol in rats fed a high-fat diet: underlying mechanism. Nutrients. 2014;6:650–62.
Siersbæk M, Varticovski L, Yang S, Baek S, Nielsen R, Mandrup S, et al. High fat diet-induced changes of mouse hepatic transcription and enhancer activity can be reversed by subsequent weight loss. Sci Rep. 2017;7:40220.
Li X, Wang Z, Klaunig JE. The effects of perfluorooctanoate on high fat diet induced non-alcoholic fatty liver disease in mice. Toxicology. 2019;416:1–14.
The authors would like to acknowledge support from the UMMC Cancer Center and Research Institute.
YG is supported by a National Institute of Health grant (P20 GM104357), by a COBRE/MS CEPR National Institute of Health grant (P20GM121334) and by a National Institute of Health grant (R01DE029803) and is supported by the Cancer Center and Research Institute. SL is partially supported by the Mississippi Center for Clinical and Translational Research and the Mississippi Center of Excellence in Perinatal Research COBRE funded by the National Institute of General Medical Sciences of the National Institutes of Health under Award Numbers 5U54GM115428 and P20GM121334. ZG is supported by grants from the Ministry of Education, Singapore (R154000B70114 and R154000B88112).
Ethics approval and consent to participate
All zebrafish studies were approved by the Deakin University Animal Welfare Committee (AWC G17-2015) and the Institutional Animal Care and Use Committee of the National University of Singapore (Protocol 079/07).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Expression profiles for lipid classes not associated with HCC disease progression. Additional lipid classes were detected in HCC zebrafish control and transgenic livers over the 1.5 to 7.5 dpt time course but were not associated with HCC disease progression by the MFP analysis. A. dhCer, B. Cer, C. MHC, D. THC, E. SM, F. LPAF, G. LPC, H. PC, I. PC(P), J. PE(O), K. PE(P), L. LPE, M. LPI, N. PG, O. CE, P. COH, Q. DG. Color key: x-m-, brown column; x-m+, blue column; x+m-, red column; x+m+, green column. “*” = p < 0.05; “**” = p < 0.01
Expression profiles for lipid species after transgene induction. A-C. Total number of lipid species showing a decrease (dark gray) or an increase (light gray) in all transgene types from 1.5 to 7.5 dpt. D. Heat map showing all lipid species detected at 7.5 dpt in all transgene types compared to DOX only (no transgene) control liver tissues. A red band shows an increase of a particular lipid species in cancer cells, a green band shows a decrease of a particular lipid species, and a dark band indicated no change. Data expressed as log2 transformation of the value of lipid species/value of control lipid species.
Univariate modelling analysis for dhCer, Cer, MHC, THC and SM lipid classes. Several lipid classes analyzed with univariate modelling analysis were not associated with HCC disease progression by MFP analysis. A. dhCer, B. Cer, C. MHC, D. THC, E. SM. X-axis units: pmol/mg. Key: “OR” (odds ratio); dark blue lines (normal); dark red lines (hyperplasia/ adenoma); green lines (HCC stage I); orange lines (HCC stage II).
Univariate modelling analysis for DG, PC(P), LPAF, LPC and PC lipid classes. Several lipid classes analyzed with univariate modelling analysis were not associated with HCC disease progression by MFP analysis. A. DG, B. PC(P), C. LPAF, D. LPC, E. PC. X-axis units: pmol/mg. Key: “OR” (odds ratio); dark blue lines (normal); dark red lines (hyperplasia/ adenoma); green lines (HCC stage I); orange lines (HCC stage II).
Univariate modelling analysis for LPE, PE(P), LPI, PE(O) and PG lipid classes. Several lipid classes analyzed with univariate modelling analysis were not associated with HCC disease progression by MFP analysis. A. LPE, B. PE(P), C. LPI, D. PE(O), E. PG. X-axis units: pmol/mg. Key: “OR” (odds ratio); dark blue lines (normal); dark red lines (hyperplasia/ adenoma); green lines (HCC stage I); orange lines (HCC stage II).
Univariate modelling analysis for CE and COH lipid classes. Some lipid classes analyzed with univariate modelling analysis were not associated with HCC disease progression by MFP analysis. A. CE, B. COH. X-axis units: pmol/mg. Key: “OR” (odds ratio); dark blue lines (normal); dark red lines (hyperplasia/ adenoma); green lines (HCC stage I); orange lines (HCC stage II).
List of qPCR primers.
About this article
Cite this article
Monroe, J.D., Fraher, D., Huang, X. et al. Identification of novel lipid biomarkers in xmrk- and Myc-induced models of hepatocellular carcinoma in zebrafish. Cancer Metab 10, 7 (2022). https://doi.org/10.1186/s40170-022-00283-y