Skip to main content

A novel murine model of combined hepatocellular carcinoma and intrahepatic cholangiocarcinoma


Primary liver cancer (PLC) is a common gastrointestinal malignancy worldwide. While hepatocellular carcinoma (HCC) and intrahepatic cholangiocarcinoma (ICC) are two major pathologic types of PLC, combined HCC and ICC (cHCC-ICC) is a relatively rare subtype that shares both hepatocyte and cholangiocyte differentiation. However, the molecular feature of this unique tumor remains elusive because of its low incidence and lack of a suitable animal model. Herein, we generated a novel spontaneous cHCC-ICC model using a Sleeping Beauty-dependent transposon plasmid co-expressing oncogenic Myc and AKT1 and a CRISPR-Cas9 plasmid expressing single-guide RNA targeting p53 into mouse hepatocytes via in situ electroporation. The histological and transcriptional analysis confirmed that this model exhibits cHCC-ICC features and activates pathways committing cHCC-ICC formation, such as TGF-β, WNT, and NF-κB. Using this model, we further screened and identified LAMB1, a protein involved in cell adhesion and migration, as a potential therapeutic target for cHCC-ICC. In conclusion, our work presents a novel genetic cHCC-ICC model and provides new insights into cHCC-ICC.


Primary liver cancer (PLC), including hepatocellular carcinoma (HCC) and intrahepatic cholangiocarcinoma (ICC), represents the third cause of cancer-related death [1, 2]. Research has revealed that ICC and HCC share a monoclonal origin with bidirectional phenotype differentiation and may appear simultaneously [3]. Combined HCC and ICC (cHCC-ICC) is a rare tumor accounting for 0.4–14.2% of all PLC [4,5,6,7]. Since the diagnosis of cHCC-ICC relies on evidence of histological findings and patients who are not suitable for resection may be misdiagnosed with HCC or ICC alone, the actual incidence underestimates the current cHCC-ICC burden [5]. Although curative surgical resection or liver transplantation is considered the mainstay of clinical practice, up to 80% of patients relapse within 5 years due to lymph node metastasis and vascular invasion [8, 9]. Unfortunately, the 5-year survival in patients with unresectable cHCC-ICC does not exceed 30% due to inadequate response to current treatments [10, 11]. Taking these findings into account when making treatment decisions, a more focused understanding of the molecular pathology and identification of potential therapeutic targets of cHCC-ICC thus are urgently needed.

Compared to xenograft tumor models, genetically engineered tumor models develop de novo tumors that closely imitate the histopathological features of their human counterparts [12]. Hydrodynamic tail vein injection (HDVI) and in situ electroporation (Epo) are two ways to transfer foreign plasmid DNA directly into hepatocytes. HDVI can create a pressurized blood force that redirects the blood flow directly into the liver, leading to plasmid DNA entering the intracellular compartment of hepatocytes [13], while electroporation is an efficient way to introduce foreign genes into cultured cells and able to in situ transfer plasmids into hepatocytes [14].

Recently, high-throughput genomic studies have revealed a series of driver genes contributing to HCC or ICC tumorigenesis [15, 16]. c-Myc and TP53 are two top frequently mutated genes in HCC patients. Recently, co-delivery of c-Myc-encoding plasmid and CRISPR/Cas9-mediated p53 knockout via HDVI successfully developed spontaneous HCC in mice [17, 18]. In another study, Seehawer M et al. constructed a vector that co-expressed Myc and AKT1 to establish HCC in p19Arf−/− mice via HDVI [19]. Interestingly, the same vector led to ICC tumorigenesis by the approach of Epo, which could cause in situ necroptosis microenvironment, highlighting the hepatic microenvironment may contribute to lineage commitment during tumorigenesis [19].

In the current study, we applied a Sleeping Beauty-dependent transposon plasmid co-expressing oncogenic Myc and AKT1 in combination with a CRISPR-Cas9 plasmid expressing single-guide RNA targeting p53 to compare their tumorigenic capacity via either HDVI- or Epo-dependent hepatocyte delivery. Notably, we found transfection of these plasmids by Epo led to the cHCC-ICC formation. Taking advantage of this novel spontaneous model, we proposed that LAMB1 may serve as a therapeutic target for cHCC-ICC.

Materials and methods


A plasmid that co-expressing oncogenic Myc and AKT1 was a kind gift and has been constructed and described previously [19]. The SB13 transposase-encoding vector was kindly provided by Dr Yue Zhao. pX330 backbone expressing sgRNA targeting p53 was obtained from Tyler Jacks (Addgene plasmid #59910).

Animal studies

Male 4 to 6-week-old C57BL/6 J mice were purchased from Charles River (Shanghai, China), and all the animals used in the study were fed in a specific pathogen-free facility. All animal care and experimental protocols were approved by the Institutional Animal Care and Use Committee (IACUC) of Zhongshan Hospital, Fudan University.

Hydrodynamic tail vein injection and in situ electroporation

For hydrodynamic tail vein injection, 30 µg Myc + AKT1 co-expressing plasmid, 30 µg pX330 sg-p53 plasmid, and 10 µg SB13 transposase-encoding plasmid were prepared in 2 ml of sterile PBS and injected into a tail vein within 3–5 s per mouse. For in situ electroporation, 6-week-old wild-type C57BL/6 J mice were anesthetized, and the right lateral liver lobe was exposed after midline laparotomy. Plasmids described above were resolved in 50 µg sterile PBS and injected into the right lateral liver lobe using an insulin needle. In situ electroporation was performed with Squaure Wave Electroporator (Nepa Gene). The voltage and duration of electric pulse were 70 V and 75 ms, respectively. Two pulses were applied, and the interval was 500 ms.


Immunohistochemistry (IHC) was performed as previously described [20]. In brief, the sections of tumors were incubated with the following antibodies: HNF4α (ab201460, Abcam), CK-19 (ab52625, Abcam), PCNA (2586, Cell Signaling Technology), p-AKT (Ser473, 4060, Cell Signaling Technology), p-ERK (Thr202/Tyr204, 4370, Cell Signaling Technology), p-NF-κB (Ser536, 3033, Cell Signaling Technology), β-catenin (8480, Cell Signaling Technology), TGF-β (21898-1-AP, Proteintech) and FGFR2 (13042-1-AP, Proteintech).

RNA-sequencing and analysis

RNA-sequencing for whole transcriptome analysis was performed using Illumina NovaSeq 6000 platform according to the manufacturer’s protocol. Three biological replicates were applied for RNA-seq, and all RNA samples passed quality control with 5–8 Gb and Q20 ≥ 90. Hierarchical clustering of RNA-seq was performed using R language with the ‘pheatmap’ package. A Euclidean method was used to calculate distance measurements, while the ‘complete’ method was utilized to calculate the dissimilarity values for hierarchical clustering. For the heatmap visualization, gene expression values were normalized.

Bioinformatic analysis

Gene set enrichment analysis (GSEA) was performed using GSEA software version 4.1.0 with 1,000 permutations of the gene sets. The FPKM values from the RNA-seq were compared against the specific gene sets. Gene sets used in this study were downloaded from the MSigdb database (


The statistical analysis was performed using Prism Graphpad 7.0 software. Quantitative variables were analyzed by paired t-test. Kaplan-Meier analysis was used to compare OS between HDVI and Epo groups. Data were presented as mean ± standard deviation. p < 0.05 was considered statistically significant.


Generation of spontaneous liver cancer model through HDVI and Epo

Myc, AKT1, and TP53 are among the top frequent mutated genes in PLC. To assess their roles in liver cancer tumorigenesis, a vector expressing constitutively active oncogenic Myc and AKT1, and single-guild RNA targeting Tp53 (sgTp53) was transfected into mouse liver cells based on a Sleeping Beauty (SB) transposon system via either HDVI or Epo (Fig. 1A, B). Both methods induced liver cancer 3–4 weeks after vector delivery, showing comparable overall survival rates (Fig. 1C). Tumors generated by HDVI (termed as MAPHDVI) exhibited multifocal neoplasia, while tumors generated by Epo (termed as MAPEpo) exhibited a unilocular pattern (Fig. 1D). H&E staining revealed disorderly and mitotically activated epithelial cells with abnormal pleomorphic nuclei and loss of polarity in all MAPHDVI tumors (Fig. 1E). Also, bile duct-like tissues were observed in MAPEpo tumors rather than in MAPHDVI tumors (Fig. 1E). IHC staining confirmed that MAPHDVI were HCC, as evidenced by solid nuclear staining for hepatocyte nuclear factor 4-alpha (HNF4α), a liver-specific marker, but negative staining for cytokeratin 19 (CK-19), a marker of biliary differentiation. Interestingly, staining of both HNF4α and CK-19 in MAPEpo tumors with a clear boundary confirmed a combined HCC and ICC within the same tumor (Fig. 1E). These data indicate that induction of Myc and AKT1 and loss of TP53 by HDVI lead to spontaneous HCC formation while Epo leads to cHCC-ICC.

Fig. 1
figure 1

Generation of spontaneous primary liver cancer model using hydrodynamic tail vein injection (HDVI) and in situ electroporation (Epo). A Schematic of plasmids designed. B Schematic of intrahepatic delivery of plasmids through HDVI or Epo. C The overall survival rate of MAPHDVI and MAPEpo mice, respectively. D Representative images of livers from mice received HDVI- or Epo-dependent plasmids via intrahepatic transfection. E Representative hematoxylin and eosin (H&E), HNF4α, and CK-19 staining of liver tissues shown in Fig. 1D. Scale bar, 200 μm

Transcriptional analysis of the MAPHDVI and MAPEpo tumors

To further validate the morphological findings, RNA-seq analysis for MAPHDVI and MAPEpo was performed to investigate the transcriptome differences. Tumor tissues from MAPHDVI and MAPEpo, along with normal liver tissue from wild-type C57BL/6 J mice as control, were firstly assessed for hepatocyte and biliary differentiation markers. The results suggested that MAPHDVI tumors exhibited higher expression of hepatic markers, including Afp, Alb, Adh1, Hnf4a, Onecut1, F2, Aldob, and Fabp1. However, genes involved in biliary differentiation, including Krt19, Krt7, Ehf, Gprc5a, Nes, Tgfb1, Tgfb2, and Jag1, were highly expressed in MAPEpo tumors but not in MAPHDVI tumors [21] (Fig. 2A). Further GESA analysis confirmed that liver-specific gene signature was significantly enriched in normal liver tissues and MAPHDVI tumors rather than in MAPEpo tumors (Fig. 1B). A previous study transcriptionally profiled 20 cHCC-ICC and successfully identified a series of gene sets that could represent cHCC-ICC and HCC genome features [22]. cHCC-ICC-related gene set comprised 588 up-regulated genes, while the HCC-related gene set contained 656 up-regulated genes. Using these gene sets, Hierarchical clustering analysis revealed that MAPEpo tumors had higher expression of genes associated with cHCC-ICC. In contrast, MAPHDVI tumors had higher expression of genes involved in HCC. Collectively, these transcriptome data further support the successful generation of spontaneous cHCC-ICC through Epo.

Fig. 2
figure 2

Transcriptome analysis of MAPHDVI and MAPEpo tumors. A Expression of genes that involves in hepatocyte and biliary differentiation. B GESA analysis using transcriptome profiles of MAPHDVI and MAPEpo tumors. Enrichment plots were presented for liver-specific gene signature. C Hierarchal clustering analysis of HCC- and cHCC-ICC–related genes was performed in MAPHDVI and MAPEpo tumors

Signaling pathways that commit cHCC-ICC activates in MAPEpo tumors

Several key signaling pathways, including the AKT, RAS, NFκB, WNT-β-catenin, and TGF-β signaling pathways, have been identified as responsible for initiating and developing cHCC-ICC [23]. We, therefore, examined whether these pathways were also activated in the MAPEpo model. GESA analysis proved that genes in the RAS, NF-κB, WNT, and TGF-β signaling pathways were significantly enriched in MAPEpo tumors rather than in normal liver tissues (Fig. 3A). Although the enrichment score produced no statistical discrepancy, these pathways tended to be relatively enriched in MAPEpo tumor rather than in MAPHDVI tumor, probably because of the mixed HCC and ICC within the tumor (Fig. 3A). IHC analysis for HCC and ICC components was subsequently performed in MAPEpo tumors to dissect this phenotype better. A comparable PCNA staining pattern among MAPHDVI, HCC-MAPEpo, and ICC-MAPEpo tumors was observed, suggesting a similar proliferation in these tumor cells (Fig. 3B). Notably, IHC staining was strongly positive for p-AKT, p-ERK, p–NF-κB, β-catenin, and TGF-β in the ICC component of the MAPEpo tumors, where FGFR2 served as the positive control (Fig. 3B). Taken together, these data suggest that the signaling pathways that contributed to cHCC-ICC formation and development are activated in MAPEpo tumors.

Fig. 3
figure 3

Signaling pathways activated in MAPHDVI and MAPEpo tumors. A GESA analysis of gene signatures of indicated pathways in normal liver tissues, MAPHDVI and MAPEpo tumors. B Representative images of PCNA, p-AKT, p-ERK, p-NFκB, β-catenin, TGF-β, and FGFR2 staining in MAPHDVI tumors and HCC or ICC component of MAPEpo tumors, respectively. Scale bars, 50 μm

Identification of LAMB1 may serve as a potential therapy target for cHCC-ICC

As our preliminary data implied that MAPEpo tumor resembled human cHCC-ICC, we, therefore, utilized this model to explore its clinical relevance. Differentially expressed genes (Log2FC ≥ 1 or ≤ − 1, p < 0.05) between MAPEpo tumor and MAPHDVI tumor were profiled to find genes involved in cHCC-ICC but not in HCC. The results indicated that 552 genes were overexpressed while 568 genes were downregulated in MAPEpo tumor compared to that in MAPHDVI (Fig. 4A). To further verify the consistency of our findings, previously reported upregulated genes in cHCC-ICC and upregulated genes in MAPEpo tumors were compared. A total of 73 overlapping genes were identified (Fig. 4B).

Fig. 4
figure 4

LAMB1 was upregulated in cHCC-ICC. A Volcano plot analysis of 1120 differential genes expressed between MAPHDVI and MAPEpo tumors. The criteria of analysis were Log2FC ≥ 1 or ≤ -1, p < 0.05. Genes that upregulated in MAPEpo tumors were shown in blue plots, while genes that upregulated in MAPHDVI tumors were in red plots. B Venn gram presented 73 overlapping genes between 552 upregulated genes in MAPEpo tumors and 588 upregulated genes in the GSE35306 data set. C GO and D KEGG analysis of 73 genes obtained in Fig. 4B. E The expression of LAMB1 gene in human ICC, HCC, and cHCC-ICC from GSE35306, GSE84073, and GSE15765 datasets

Using these overlapping genes as input, Gene Oncology analysis revealed that biological processes (BP) were mainly enriched in tumor metastasis-related changes such as extracellular matrix organization (ECM), cell adhesion, and cell migration (Fig. 4C). Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis demonstrated that focal adhesion and ECM-receptor interaction were among the top two pathways enriched (Fig. 4D). These results were consistent with the invasive feature of cHCC-ICC.

In the 73 overlapping genes, Lamb1, which contributes to ECM activity and cell adhesion and is highly expressed in cHCC-ICC, was identified. We also found that LAMB1 expression in human cHCC-ICC (GSE35306) was significantly higher than in both ICC and HCC tissues. For further validation, we analyzed another two published datasets (GSE84073 and GSE15765) that contain transcriptome data of human ICC, HCC, and cHCC-ICC. It also exhibited that LAMB1 was upregulated in human cHCC-ICC compared with HCC (Fig. 4E). Considering minimal data on cHCC-ICC are available, no significant changes of Lamb1 were observed between cHCC-ICC and ICC. To sum up, these findings imply that Lamb1 may serve as a therapeutic target for cHCC-ICC, and further study is needed to dissect its role in cHCC-ICC.


cHCC-ICC is a rare type of PLC but has attracted increasing attention in recent years. Because the diagnosis of cHCC-ICC largely depends on histochemistry, the true incidence of cHCC-ICC is likely to be underestimated, making the knowledge and management of cHCC-ICC inaccessible. Thus, a tumor animal model that resembles clinical cHCC-ICC is essential for further investigation. In this study, we generated a spontaneous cHCC-ICC by induction of oncogenic Myc and AKT1 and loss of p53 through in situ electroporation. Evidence showed that this model shared similar transcriptome and oncogenic signaling pathways with human cHCC-ICC. Most importantly, we, for the first time, identified LAMB1 as a potential therapy target for cHCC-ICC.

Most solid tumors, especially HCC, develop in the context of chronic diseases and are composed of heterogeneous malignant cells. Importantly, the heterogeneous feature of tumor cells is a key reason for clinical drug resistance. Although traditional syngeneic or xenograft models are easy to perform, these models are unable to fully mimic specific human disease conditions and neglect the heterogeneous feature of tumor tissues. Thus, these common models have a limit to evaluating the drugs pre-clinically. It is well known that tumor cells are transformed from normal cells carrying oncogenic mutations, making it possible to induce PLC by genetic engineering. Accumulative studies have revealed specific mutations that lead to either HCC or ICC tumorigenesis. For example, c-Myc, CTNNB1 and TP53 are among top mutated genes for HCC patients [24]; while KRAS, a powerful oncogene involved in glandular malignant, is able to induce ICC when specifically expressed in mouse hepatocytes [25]. Currently, genetic engineered mouse model (GEMM) is an ideal tool for study HCC or ICC tumorigenesis and has clear genetic background resembling human disease, making it suitable for pre-clinical estimation. Unfortunately, few studies focus on cHCC-ICC as it’s not entirely clear how cHCC-ICC occurs. A recent study proved that an inflammatory tumor microenvironment directs lineage commitment of the PLC [19]. Here, we used electroporation to induce a necroptotic liver microenvironment and established cHCC-ICC formation by combing Myc and AKT1 knockin with p53 knockout. The morphological and genetic evidences proved successful induction of cHCC-ICC. Our work provides a simple GEMM for cHCC-ICC with direct clinical translational value. However, our study also has limitations. It remains elusive whether other gene combinations can lead to the cHCC-ICC formation in the same experimental setting. Moreover, our model does not exhibit the metastatic feature of cHCC-ICC, as evidenced by no lung metastatic lesion observed (data not shown). New gene combinations that can lead to metastasis should be tested in the future.

Whether cHCC-ICC is a unique or a subtype of HCC or ICC has long been controversial. cHCC-ICC can be further divided into three subtypes according to Allen and Lisa’s criteria that are separate type (HCC and ICC components physically separated), combined type (HCC and ICC component in the same tumor with clear boundaries), and mixed type (HCC and ICC component in the same tumor with no boundaries) [26]. Histologically, we found that MAPEpo tumors exhibited features with a combined type of cHCC-ICC. Recently, Xue et al. [3] comprehensively analyzed a total of 133 cHCC-ICC cases and revealed that combined type of cHCC-ICC acquired intense ICC-like landscapes, including high expression of KRT19 but a lower expression of HCC markers (including AFP and GPC3), which is also supported by our transcriptome data (Fig. 2A). Compared with MAPHDVI tumors, MAPEpo tumors tended to express markers of biliary differentiation, especially for Krt19, Krt7, Nes, Tgfb2, and Jag1. Consistently, Afp expression is much higher in MAPHDVI tumors than in MAPEpo tumors. These data suggested that MAPEpo tumors might exhibit more ICC-like characteristics.

Despite apparent molecular discrepancy among subtypes of cHCC-ICC, they all have poorer prognosis and more invasive features than HCC and are similar to ICC [27]. Interactions between cell adhesion or migration and ECM are vital factors that mediate tumor metastasis. In our study, we identified 73 overlapping genes that are highly upregulated in cHCC-ICC. Bioinformatic analysis indicated that these genes were mainly involved in ECM bioprocesses and cell adhesion. These findings further reinforced the invasive feature of cHCC-ICC.

Laminins, a family of extracellular matrix glycoproteins, are among the predominant component of ECM [28]. Evidence reported that laminins participated in tumor metastasis by promoting cell adhesion and migration, and their receptors expressed on the tumor cell surface [29,30,31]. We screened these 73 genes and identified a main differential gene Lamb1, a member of the laminin family. Taking advantage of several data sets reported previously, we established that Lamb1 was upregulated in cHCC-ICC compared to HCC, but its expression seemed comparable with ICC. It has been reported that LAMB1 overexpressed in several types of tumors and correlated with tumor metastasis and poor prognosis [32,33,34]. The prognosis value of LAMB1 was further screened tin TCGA database. The results suggested that LAMB1 is elevated in a series of tumors, but LAMB1 is negatively correlated with both OS and DFS for LIHC, CHOL, and COAD (Additional file 1: Fig. S1). Most importantly, KEGG and GSEA analysis showed that LAMB1 played a pivotal role in pro-metastatic processes, including focal adhesion, ECM-receptor interaction, and cellular junction in LIHC, CHOL, and COAD, highlighting that LAMB1 might also be a potential treating target for cHCC-ICC (Additional file 1: Fig. S2). Since clinical transcriptome data of patients with cHCC-ICC remains very limited, we were able only to verify the expression of Lamb1 in current published data; further larger-scale studies should be performed to assess its clinical value. Summarily, the present study established a human-resembling cHCC-ICC model via in situ electroporation. This novel preclinical model can be used to investigate the molecular feature of cHCC-ICC.

Availability of data and materials

The datasets used and analyzed during the current study are available from the corresponding author on reasonable request.


  1. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, Bray F. Global Cancer Statistics 2020: GLOBOCAN estimates of incidence and Mortality Worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71:209–49.

    Article  Google Scholar 

  2. Zhou J, Sun H, Wang Z, Cong W, Wang J, Zeng M, Zhou W, Bie P, Liu L, Wen T, et al. Guidelines for the diagnosis and treatment of hepatocellular carcinoma (2019 Edition). Liver Cancer. 2020;9:682–720.

    Article  Google Scholar 

  3. Xue R, Chen L, Zhang C, Fujita M, Li R, Yan SM, Ong CK, Liao X, Gao Q, Sasagawa S, et al. Genomic and transcriptomic profiling of combined hepatocellular and intrahepatic cholangiocarcinoma reveals distinct molecular subtypes. Cancer Cell. 2019;35:932-47.e938.

    Article  CAS  Google Scholar 

  4. Ramai D, Ofosu A, Lai JK, Reddy M, Adler DG. Combined hepatocellular cholangiocarcinoma: a population-based retrospective study. Am J Gastroenterol. 2019;114:1496–501.

    Article  Google Scholar 

  5. Garancini M, Goffredo P, Pagni F, Romano F, Roman S, Sosa JA, Giardini V. Combined hepatocellular-cholangiocarcinoma: a population-level analysis of an uncommon primary liver tumor. Liver Transpl. 2014;20:952–9.

    Article  Google Scholar 

  6. Jarnagin WR, Weber S, Tickoo SK, Koea JB, Obiekwe S, Fong Y, DeMatteo RP, Blumgart LH, Klimstra D. Combined hepatocellular and cholangiocarcinoma: demographic, clinical, and prognostic factors. Cancer. 2002;94:2040–6.

    Article  Google Scholar 

  7. Lee CC, Wu CY, Chen JT, Chen GH. Comparing combined hepatocellular-cholangiocarcinoma and cholangiocarcinoma: a clinicopathological study. Hepatogastroenterology. 2002;49:1487–90.

    Google Scholar 

  8. Yoon YI, Hwang S, Lee YJ, Kim KH, Ahn CS, Moon DB, Ha TY, Song GW, Jung DH, Lee JW, et al. Postresection outcomes of combined hepatocellular carcinoma-cholangiocarcinoma, hepatocellular carcinoma and intrahepatic cholangiocarcinoma. J Gastrointest Surg. 2016;20:411–20.

    Article  Google Scholar 

  9. Yin X, Zhang BH, Qiu SJ, Ren ZG, Zhou J, Chen XH, Zhou Y, Fan J. Combined hepatocellular carcinoma and cholangiocarcinoma: clinical features, treatment modalities, and prognosis. Ann Surg Oncol. 2012;19:2869–76.

    Article  Google Scholar 

  10. Yamashita YI, Aishima S, Nakao Y, Yoshizumi T, Nagano H, Kuroki T, Takami Y, Ide T, Ohta M, Takatsuki M, et al. Clinicopathological characteristics of combined hepatocellular cholangiocarcinoma from the viewpoint of patient prognosis after hepatic resection: high rate of early recurrence and its predictors. Hepatol Res. 2020;50:863–70.

    Article  CAS  Google Scholar 

  11. Magistri P, Tarantino G, Serra V, Guidetti C, Ballarin R, Di Benedetto F. Liver transplantation and combined hepatocellular-cholangiocarcinoma: feasibility and outcomes. Dig Liver Dis. 2017;49:467–70.

    Article  Google Scholar 

  12. Brown ZJ, Heinrich B, Greten TF. Mouse models of hepatocellular carcinoma: an overview and highlights for immunotherapy research. Nat Rev Gastroenterol Hepatol. 2018;15:536–54.

    Article  CAS  Google Scholar 

  13. Chen X, Calvisi DF. Hydrodynamic transfection for generation of novel mouse models for liver cancer research. Am J Pathol. 2014;184:912–23.

    Article  CAS  Google Scholar 

  14. Suzuki T, Shin BC, Fujikura K, Matsuzaki T, Takata K. Direct gene transfer into rat liver cells by in vivo electroporation. FEBS Lett. 1998;425:436–40.

    Article  CAS  Google Scholar 

  15. Dong LQ, Shi Y, Ma LJ, Yang LX, Wang XY, Zhang S, Wang ZC, Duan M, Zhang Z, Liu LZ, et al. Spatial and temporal clonal evolution of intrahepatic cholangiocarcinoma. J Hepatol. 2018;69:89–98.

    Article  CAS  Google Scholar 

  16. Gao Q, Zhu H, Dong L, Shi W, Chen R, Song Z, Huang C, Li J, Dong X, Zhou Y, et al. Integrated proteogenomic characterization of HBV-related hepatocellular carcinoma. Cell. 2019;179:561-77.e522.

    Article  CAS  Google Scholar 

  17. Ruiz de Galarreta M, Bresnahan E, Molina-Sánchez P, Lindblad KE, Maier B, Sia D, Puigvehi M, Miguela V, Casanova-Acebes M, Dhainaut M, et al. β-Catenin activation promotes Immune escape and resistance to Anti-PD-1 therapy in hepatocellular carcinoma. Cancer Discov. 2019;9:1124–41.

    Article  CAS  Google Scholar 

  18. Hu B, Yu M, Ma X, Sun J, Liu C, Wang C, Wu S, Fu P, Yang Z, He Y, et al. IFNα potentiates Anti-PD-1 efficacy by remodeling glucose metabolism in the hepatocellular carcinoma microenvironment. Cancer Discov. 2022;12:1718–41.

    Article  CAS  Google Scholar 

  19. Seehawer M, Heinzmann F, D’Artista L, Harbig J, Roux PF, Hoenicke L, Dang H, Klotz S, Robinson L, Dore G, et al. Necroptosis microenvironment directs lineage commitment in liver cancer. Nature. 2018;562:69–75.

    Article  CAS  Google Scholar 

  20. Sun J, Zhou C, Zhao Y, Zhang X, Chen W, Zhou Q, Hu B, Gao D, Raatz L, Wang Z, et al. Quiescin sulfhydryl oxidase 1 promotes sorafenib-induced ferroptosis in hepatocellular carcinoma by driving EGFR endosomal trafficking and inhibiting NRF2 activation. Redox Biol. 2021;41:101942.

    Article  CAS  Google Scholar 

  21. Lin YK, Fang Z, Jiang TY, Wan ZH, Pan YF, Ma YH, Shi YY, Tan YX, Dong LW, Zhang YJ, Wang HY. Combination of Kras activation and PTEN deletion contributes to murine hepatopancreatic ductal malignancy. Cancer Lett. 2018;421:161–9.

    Article  CAS  Google Scholar 

  22. Coulouarn C, Cavard C, Rubbia-Brandt L, Audebourg A, Dumont F, Jacques S, Just PA, Clément B, Gilgenkrantz H, Perret C, Terris B. Combined hepatocellular-cholangiocarcinomas exhibit progenitor features and activation of wnt and TGFβ signaling pathways. Carcinogenesis. 2012;33:1791–6.

    Article  CAS  Google Scholar 

  23. Wang G, Wang Q, Liang N, Xue H, Yang T, Chen X, Qiu Z, Zeng C, Sun T, Yuan W, et al. Oncogenic driver genes and tumor microenvironment determine the type of liver cancer. Cell Death Dis. 2020;11:313.

    Article  CAS  Google Scholar 

  24. Friemel J, Frick L, Unger K, Egger M, Parrotta R, Boge YT, Adili A, Karin M, Luedde T, Heikenwalder M, Weber A. Characterization of HCC Mouse Models: towards an etiology-oriented Subtyping Approach. Mol Cancer Res. 2019;17:1493–502.

    Article  CAS  Google Scholar 

  25. Tanaka M, Kunita A, Yamagishi M, Katoh H, Ishikawa S, Yamamoto H, Abe J, Arita J, Hasegawa K, Shibata T, Ushiku T. KRAS mutation in intrahepatic cholangiocarcinoma: linkage with metastasis-free survival and reduced E-cadherin expression. Liver Int. 2022;42:2329–40.

    Article  CAS  Google Scholar 

  26. Allen RA, Lisa JR. Combined liver cell and bile duct carcinoma. Am J Pathol. 1949;25:647–55.

    CAS  Google Scholar 

  27. Wakizaka K, Yokoo H, Kamiyama T, Ohira M, Kato K, Fujii Y, Sugiyama K, Okada N, Ohata T, Nagatsu A, et al. Clinical and pathological features of combined hepatocellular-cholangiocarcinoma compared with other liver cancers. J Gastroenterol Hepatol. 2019;34:1074–80.

    Article  CAS  Google Scholar 

  28. Patarroyo M, Tryggvason K, Virtanen I. Laminin isoforms in tumor invasion, angiogenesis and metastasis. Semin Cancer Biol. 2002;12:197–207.

    Article  CAS  Google Scholar 

  29. Liang Y, Chen X, Wu Y, Li J, Zhang S, Wang K, Guan X, Yang K, Bai Y. LncRNA CASC9 promotes esophageal squamous cell carcinoma metastasis through upregulating LAMC2 expression by interacting with the CREB-binding protein. Cell Death Differ. 2018;25:1980–95.

    Article  CAS  Google Scholar 

  30. Dai J, Cimino PJ, Gouin KH, Grzelak CA, Barrett A, Lim AR, Long A, Weaver S, Saldin LT, Uzamere A, et al. Astrocytic laminin-211 drives disseminated breast tumor cell dormancy in brain. Nat Cancer. 2022;3:25–42.

    Article  CAS  Google Scholar 

  31. Reuten R, Zendehroud S, Nicolau M, Fleischhauer L, Laitala A, Kiderlen S, Nikodemus D, Wullkopf L, Nielsen SR, McNeilly S, et al. Basement membrane stiffness determines metastases formation. Nat Mater. 2021;20:892–903.

    Article  CAS  Google Scholar 

  32. Shan H, Liu T, Gan H, He S, Deng J, Hu X, Li L, Cai L, He J, Long H, et al. RNA helicase DDX24 stabilizes LAMB1 to promote hepatocellular carcinoma progression. Cancer Res. 2022.

    Article  Google Scholar 

  33. Chen Z, Yang X, Bi G, Liang J, Hu Z, Zhao M, Li M, Lu T, Zheng Y, Sui Q, et al. Ligand-receptor interaction atlas within and between tumor cells and T cells in lung adenocarcinoma. Int J Biol Sci. 2020;16:2205–19.

    Article  CAS  Google Scholar 

  34. Meves A, Nikolova E, Heim JB, Squirewell EJ, Cappel MA, Pittelkow MR, Otley CC, Behrendt N, Saunte DM, Lock-Andersen J, et al. Tumor cell adhesion as a risk factor for sentinel lymph node metastasis in primary cutaneous melanoma. J Clin Oncol. 2015;33:2509–15.

    Article  CAS  Google Scholar 

Download references


Not applicable.


This study was partly supported by the National Natural Science Foundation of China (No. 81472673, and No. 81672720).

Author information

Authors and Affiliations



Concept, drafting and interpretation: RCX, FW, and JMZ. Execution and analysis: RCX, FW, JLS, WA, GCZ, ZYL, TTL, LD, XZS, and JMZ. All authors contributed to the article and approved the submitted version. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Ji-Min Zhu.

Ethics declarations

Ethics approval and consent to participate

All animal experiments were approved by the Animal Ethics and Welfare Committee of Zhongshan Hospital, Fudan University.

Consent for publication

All authors have approved the manuscript for submission.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1. Figure S1.

Multi-cancer analysis of the expression and prognostic role of LAMB1. A LAMB1 is significantly upregulated in multiple cancer types from TCGA data. B Kaplan-Meier analysis of the association between LAMB1 expression and OS or DFS in LIHC, CHOL, and COAD datasets. Figure S2. Multi-cancer analysis of biological functions and significant pathway of LAMB1. A Top 20 pathways enriched in the KEGG analysis in LIHC, CHOL, and COAD datasets. The red box represents pro-metastatic processes, including focal adhesion, ECM-receptor interaction, and cellular junction. B GESA analysis of gene signatures of indicated pathways in LIHC, CHOL, and COAD datasets.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Xu, RC., Wang, F., Sun, JL. et al. A novel murine model of combined hepatocellular carcinoma and intrahepatic cholangiocarcinoma. J Transl Med 20, 579 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Myc
  • AKT1
  • p53
  • LAMB1
  • Sleeping Beauty-dependent transposon system
  • In situ electroporation