- Open Access
Custom gene expression panel for evaluation of potential molecular markers in hepatocellular carcinoma
BMC Medical Genomics volume 15, Article number: 235 (2022)
Hepatocellular carcinoma (HCC) is the second leading cause of cancer-related mortality worldwide. It is a highly heterogeneous disease with poor prognosis and limited treatment options, which highlights the need for reliable biomarkers. This study aims to explore molecular markers that allow stratification of HCC and may lead to better prognosis and treatment prediction.
Materials and methods
We studied 20 candidate genes (HCC hub genes, potential drug target genes, predominant somatic mutant genes) retrieved from literature and public databases with potential to be used as the molecular markers. We analysed expression of the genes by RT-qPCR in 30 HCC tumour and adjacent non-tumour paired samples from Vietnamese patients. Fold changes in expression were then determined using the 2−∆∆CT method, and unsupervised hierarchical clustering was generated using Cluster v3.0 software.
Clustering of expression data revealed two subtypes of tumours (proliferative and normal-like) and four clusters for genes. The expression profiles of the genes TOP2A, CDK1, BIRC5, GPC3, IGF2, and AFP were strongly correlated. Proliferative tumours were characterized by high expression of the c-MET, ARID1A, CTNNB1, RAF1, LGR5, and GLUL1 genes. TOP2A, CDK1, and BIRC5 HCC hub genes were highly expressed (> twofold) in 90% (27/30), 83% (25/30), and 83% (24/30) in the tissue samples, respectively. Among the drug target genes, high expression was observed in the GPC3, IGF2 and c-MET genes in 77% (23/30), 63% (19/30), and 37% (11/30), respectively. The somatic mutant Wnt/ß-catenin genes (CTNNB1, GLUL and LGR5) and TERT were highly expressed in 40% and 33% of HCCs, respectively. Among the HCC marker genes, a higher percentage of tumours showed GPC3 expression compared to AFP expression [73% (23/30) vs. 43% (13/30)].
The custom panel and molecular markers from this study may be useful for diagnosis, prognosis, biomarker-guided clinical trial design, and prediction of treatment outcomes.
Hepatocellular carcinoma (HCC) is the second most common cause of cancer-related mortality worldwide, with two-thirds of cases occurring in Asia . The hepatitis B (HBV) and hepatitis C (HCV) viruses are the major etiological agents for HCC globally , and in particular HBV-induced HCC is common with increasing numbers of cases in East Asia and Africa . Vietnam has a high prevalence of HBV infections, ranging from 10 to 20% in the general population and 20–40% among injecting drug users and HIV-positive patients . Prognosis of HCC remains poor because of the underlying chronic liver disease; late diagnosis, often at advanced stages of disease .
HCC is heterogeneous and classified into four Barcelona Clinic Liver Cancer (BCLC) stages (BCLC-A: early HCC, BCLC-B: intermediate HCC, BCLC-C: advanced HCC, BCLC-D: terminal HCC) for optimal clinical management using the staging depending on tumour burden, liver function (Child–Pugh score), and overall health [6, 7]. HCC is routinely diagnosed by blood tests (alanine transaminase [ALT], aspartate transaminase [AST], total bilirubin and direct bilirubin), liver biopsy, imaging, and serum alpha-fetoprotein (AFP) and albumin levels. HCC is heterogeneous, which is associated with poor prognostic outcomes because it is more difficult to diagnose accurately. Heterogeneity of HCC tumours are known to occur at different levels; firstly within the population and secondly within tumours of the same patient . Recently, rapid progress has been made in understanding the heterogeneity of HCC by molecular subclassification using molecular and genetic markers . HCC tumours are divided into two major phenotypic classes: the proliferation class and the non-proliferation class [10, 11], which have been identified on the basis of transcriptomic dysregulations and genetic alterations closely related to risk factors, pathologic features, and prognosis. This proliferation class is associated with HBV infection and has a poor clinical outcome. These are also subdivided into several molecular subclasses based on transcriptional profile, as discussed in extensively in several articles [11,12,13,14,15,16]. This classification would allow identification of predictive biomarkers and facilitate targeted treatments.
The pathogenesis of HCC is a multistep process involving progressive accumulation of genetic alterations in distinct oncogenes. Progression of HCC begins with chronic inflammation caused by viral (HBV, HCV) or non-viral factors that are known to alter the liver microenvironment. As a result, the liver cells increase the production of cytokines, reactive oxygen and nitrogen species that mediate liver injury and trigger the liver's regenerative response. This regeneration predisposes the hepatocyte cells to a variety of genetic alterations at the genomic and transcriptional levels leading to HCC [17,18,19,20]. The most commonly altered genetic variants in HCC include mutations of the TERT promoter, CTNNB1, TP53, RB1, CCNA2, CCNE1, PTEN, ARID1A, ARID2, RPS6KA3 or NFE2L2, CCND1, FGF19, VEGFA, MYC and MET [21, 22], leading to changes in TERT, Wnt/β-catenin, p53/p21 and RB1, AKT-mTOR, RAS-MAPK, VEGF/VEGFR, MET, IGF, ARID1A/ARID1B/ARID2 pathways. Understanding the pathways and genes in HCC plays an important role in cancer prognosis and prediction of treatment success through stratification of cases for personalized treatment.
Despite some recent success in the development of new drugs approved for treatment of advanced HCC , several clinical trials have failed or shown only modest improvement in overall survival . Further novel drug combinations and treatment regimens with first- and second-line drugs are currently under investigation . As HCC is clinically and molecularly very heterogeneous, the use of biomarkers may guide personalized treatment strategies in clinical trials, leading to more favourable therapeutic outcomes. Recent developments in biomarker-driven therapies show some promise. To identify biomarkers associated with heterogeneous HCC tumours, several studies and multiple insilico computational bioinformatics analysis [25,26,27,28,29] using publicly available high-throughput HCC datasets have identified critical genes and pathways [17,18,19,20]. Some of the recent markers identified in the studies include; Cell division cycle-associated protein-3 (CDCA3) , ribophorins (RPNs), ARID1A . These molecular markers of HCC enable to predict prognosis and to develop a new rationale for targeted therapeutic strategies. To date, several specific markers and key pathways involved in the HCC development have been identified potentially to assist early diagnosis, to predict prognosis and molecular targeted therapies of HCC tumours have been extensively reviewed. Recent developments in biomarker-driven therapies show some promises have been reviewed [32, 33]. For example, in a phase Ib/II HCC trial, patients with high MET mRNA expression showed a threefold increase in progression-free survival compared to patients with low c-MET expression when given a combination of the anti-VEGFR-2 mAb ramucirumab plus the anti-MET mAb emibetuzumab . The gene expression profiles are used to identify differentially expressed genes (DEGs) during cancer progression [35, 36], which enables the identification of biomarkers for the prognosis, diagnosis, and targeted therapy of tumours .
To understand the heterogeneity of HCC and identify appropriate molecular markers, the present study examined a panel of genes including HCC hub, drug target, and somatic mutation genes with the potential to serve as meaningful biomarkers in prognosis and prediction of therapeutic success.
Materials and methods
The clinical and diagnostic characteristics of the study group are summarized in Table 1. To determine gene expression, pairs of HCC tumour (T) and adjacent non-tumour (NT) tissue specimens were obtained from 30 HCC patients who underwent surgery at the Vietnam National Cancer Hospital, Tan Trieu, Ha Dong in Hanoi in 2018. All patients were negative for anti-HCV and anti-HIV antibodies, nor had they a history of drug abuse. Liver function tests including ALT, AST, total bilirubin and direct bilirubin, albumin and prothrombin were quantified. HCC patients were categorized according to the BCLC staging system [38, 39]. All blood and T and NT tissue samples were stored at − 80 °C until further use.
Serum levels assessment of AFP, CEA and CA 19-9
Serum levels of the markers AFP, CEA, and CA19-9 were measured using the ARCHITECT AFP reagent kit, Cat. No. 03P3625 (Abbott Ireland Diagnostics, Ireland), ARCHITECT CEA reagent kit, Cat. No. 07K6832 (Abbott Ireland Diagnostics, Ireland), and ARCHITECT CA 19-9XR Reagent Kit, Cat. No. 2K91 (Abbott GMBH & Co. Germany) respectively on the ARCHITECT i2000sr—an automated immunoassay analyser (Abbott Diagnostics, Abbott Park, IL, USA) according to the manufacturer's protocol.
Genes selection criteria for the qPCR panel
The custom qPCR panel included 20 genes (Table 2) retrieved from scientific literature with biomarker potential for diagnosis, prognosis, and prediction of treatment outcomes. Specifically, genes belonging to the categories of HCC hub genes (prognosis and diagnostic potential), drug target genes and somatic mutated genes were considered. Another inclusion criterion was that we considered genes whose expression was significantly upregulated in at least 30–40% of HCCs (www.proteinatlas.org).
RNA isolation and real-time qPCR
Total RNA was extracted from approximately 5 mm3 of HCC tumor or adjacent normal liver tissue using 1 mL TRIzol™ reagent (Thermo Fisher Scientific, Waltham, MA, USA). Next, the tissues were homogenized with a Dounce homogenizer, and RNA was isolated according to the manufacturer’s instructions. The RNA pellet was resuspended in 30–50 µL of RNA-free water, depending on the pellet's size. Next, quality and quantity were assessed using a Nanodrop spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA). For cDNA synthesis, one microgram of total RNA from each sample was reverse transcribed in 20 µL reaction using the RevertAid First Strand cDNA Synthesis Kit (Thermo Fisher Scientific, Waltham, MA, USA) according to the manufacturer’s instructions. The cDNA was diluted to 1:10 with nuclease-free water and subjected to quantitative real-time PCR.
The quantitative real-time PCR was performed with the LightCycler 480 Instrument II (Roche, Basel, Switzerland) in 96 well plates using a 10 µL final volume containing 0.3 µmol/L of forward and reverse primers, 2 µL of diluted cDNA and 1X SensiFAST™ SYBR® No-ROX Kit (Bioline, Memphis, TN, USA). PCR cycle parameters were initial denaturation at 94 °C for 2 min, 40 cycles at 95 °C for 10 s, and 60 °C for 25 s, with the fluorescence signal acquired at the end of each cycle. A melting curve analysis was performed at the end of the amplification. All 23 genes (20 genes under investigation plus 3 housekeeping genes) were quantified in duplicates in 96 well plates having a couple of non-template reactions with ACTB primers. The primer sequences used are given in the Additional file 1: Table S1.
The fold change gene expression between the tumour (T) and adjacent non-tumour (NT) pair was determined using the relative quantification 2−∆∆CT method . The fold change with ± two-fold to the transcript levels between T and adjacent NT was considered differentially regulated. For having a better control of the results, mean cycle threshold (Ct) of the three reference genes was taken for the normalization: Ct [ref] = mean (Ct [GAPDH], Ct [ACTB], Ct [TBP]). Correlations between genes were measured using Pearson correlation (r) and differences between groups were determined by either two-tailed t-test or ANOVA. The gene expression heatmap was generated using the Cluster v3.0 software (http://bonsai.hgc.jp/~mdehoon/software/cluster/software.html) . The normalized gene-expression profiles of 30 samples were subjected to unsupervised hierarchical clustering for similarities in expression data. Hierarchical clustering was applied to both rows and columns. The output file was visualized in JavaTreeview (https://sourceforge.net/projects/jtreeview/) .
Patients baseline characteristics and blood indices
Paired samples of T and NT tissues from 29 patients with BCLC-B classification and one patient with BCLC-A classification were analysed for the expression of the genes under investigation. The baseline characteristics of the patients are summarized in Table 1. Ninety percent were > 40 years of age at the time of HCC diagnosis, and the mean age of study participants was 57 ± 12 years; 83% were men. Chronic HBV infection was the predominant incidental factor of HCC development (87%). The remaining cases (13%) were due to unknown non-HBV and non-HCV factors. Among other risk factors, liver cirrhosis was found in 23% of patients, and smoking/alcohol consumption was observed in 57% of patients. Tumours in the right liver lobe (67%) were more frequent than in the left liver lobe (33%). Large (≥ 5 cm) tumours were detected in 50% of the cohort, an intermediate size (3–5 cm) in 47% and one patient (3%) had a tumour of < 3 cm diameter. Blood tests showed elevated ALT (> 40 U/L) and AST (> 56 U/L) in 30% and 13% of patients respectively, elevated total (> 17 µmol/L) and direct bilirubin (> 5 µmol/L) levels in 17% and lowered platelet counts (< 150 × 103/mL) in 20% of the cohort.
Identification of molecular HCC subtypes by gene expression profiling
Expression patterns of 20 genes were analysed in 30 paired tumour and non-tumour (T and NT) samples by RT-qPCR. Unsupervised hierarchical clustering using Manhattan distance and complete linkage revealed two subtypes of tumours (Fig. 1). Of the two subtypes, subtype-1 (proliferative) showed higher expression of genes compared to subtype-2 (normal-like). The proliferative subtype included 22 tumour samples compared eight for the normal-like subtypes.
Clustering of the genes revealed four clusters. The hub genes TOP2A, CDK1, and BIRC5 formed the top cluster with expression in 80–90% of the tumours. In addition to the hub genes, GPC3, IGF2, and AFP also showed uniform expression in both the proliferative, but not in the normal-like subtype. Expression of AFP showed a good correlation with the corresponding AFP serum levels (r:0.63, p < 0.001). The second cluster comprised of the PD-L1, VEGFA, IDH1, MCL1 and TP53 genes, which were diffusely upregulated in the proliferative subtype. The third cluster contained MDM4, ARID1A, c-MET, RAF, and CTNNB1, which were strongly upregulated in the proliferative subtype. The last cluster consisted of LGR5, GLUL1, TERT and CDC20, which were diffusely expressed in all subtypes. The bar graph represents the percentage of tumours from the highest to the lowest number of hub genes, drug target genes, and common somatic mutant genes that had > twofold gene expression (Fig. 2).
Correlation between genes
The expression profiles of genes in the first cluster across all tumours (TOP2A, CDK1, BIRC5, GPC3, IGF2 and AFP) were strongly correlated. The correlations between TOP2A and CDK1 (r:0.77, p < 0.0001), CDK1 and BIRC5 (r;0.77, p < 0.0001) and TOP2A and BIRC5 (r:0.65, p < 0.0001) were significant. The correlations between GPC3 and IGF2 (r:0.61, p < 0.0001) and between GPC3 and AFP (r:0.78, p < 0.00001) and between GPC3 and IGF2 (r:0.55, p < 0.0001) were also significant.
The correlation between TOP2A and CDK1 remained significant in both the tumour subtypes (r:0.58, p:0.004; r:0.79, p:0.02 in proliferative, and normal-like subtypes, respectively). The correlation between TOP2A and BIRC5 (r:0.57, p:0.006), TOP2A and GPC3 (r:0.48, p:0.02) was significant in the proliferative, but not in the normal-like subtype. Expression data for MDM4 were available for 24 of the 30 tumours and its correlation with TP53 was significant in proliferative (r:0.52, p:0.04), but not in normal-like (r:0.36, p:0.5). In the second gene cluster (IDH1, VEGFA, MCL1, PD-L1, TP53), statistically significant correlations were observed in proliferative subtype for VEGFA and IDH1 (r:0.42, p:0.05), VEGFA and MCL1 (r:0.43. p:0.05), VEGFA and TP53 (r:0.52, p:0.01), MCL1 and IDH1 (r:0.57, p < 0.01), PD-l1 and TP53 (r:0.57, p < 0.01). In the normal-like subtype, correlations were weak between genes and not statistically significant. In the third cluster, the correlations between c-MET and ARID1A (r:0.74, p < 0.01), c-MET and CTNNB1 (r:0.71, p < 0.01), CTNNB1 and ARID1A (r:0.51, p:0.03) and LGR5 and GLUL1 (r: 0.69, p:0.02) were only significant for the proliferative, but not in normal-like subtype. The correlation of TERT with CDC20 was significant in both the subtypes (r:0.60, p < 0.01; r:0.80, p:0.02 in proliferative and normal-like subtypes, respectively).
Comparison of tumour subtypes
The tumour sizes and serum levels of AFP, CEA and CA19-9 were compared between the proliferative and normal-like subtypes. Comparison of tumour sizes showed that the normal-like subtype showed larger tumour sizes in the normal-like (median [IQR]: 6.5 cm [3.8–8.5]) compared to the proliferative (4.8 cm [3–7]) and the difference was not statistically significant (t-test, two-tailed p:0.19).
Comparison of serum AFP levels revealed higher levels in the normal-like subtype (median 114 ng/ml [16–2462]; (log2-transformed median 6.6 [7–11]) compared with the proliferative (median 16 ng/ml [6–486]; (log2-transformed median 4.1 [3–8])). No statistical difference in serum AFP levels was observed between the groups (p:0.30) (Fig. 3A). The normal-type subtype had higher levels of the serum marker CEA (median 3.6 ng/ml [3–4]; (log2-transformed median 1.8 [1.6–2])) compared to the proliferative subtype (median 2.3 ng/ml [1.5–4]; (log2-transformed median 1.3 [0.6–2])), although the difference between groups was not statistically significant (p:0.35) (Fig. 3B). The proliferative subtype showed higher serum CA19-9 levels (median 13.4 U/L [9–31]; (log2-transformed median 3.7 [3–5])) compared with normal-like (median 7.6 U/L [1.5–14]; (log2-transformed median 2.9 [0.4–4])), but no statistical difference was found between groups (p:0.1) (Fig. 3C).
The intratumoural heterogeneity among the BCLC-B patients is a major impediment in effective HCC therapy and renders long-term prognoses difficult and largely unreliable. The expression profile of key mutated genes and molecular classification of HCC and a subclassification for BCLC-B into four stages (B1-B4) has been proposed for better patient management. However, contradictory results from validation studies have shown a greater need of improved scoring systems and clinically relevant biomarkers [61,62,63,64,65]. To address the issue of biomarkers, our study investigated the expression pattern of important genes in HCC tumours and assessed their utility as biomarkers.
By using hierarchical clustering of gene expression data our study revealed two tumour subtypes: a proliferative and normal-like subtype. The proliferative subtype showed mutually exclusive expression patterns, except for the hub genes GPC3, IGF2 and AFP. This observation is in agreement with findings published previously [44, 66, 67]. The proliferative subtype accounts for 50% of HCC cases and is associated with a highly aggressive phenotype and poor outcome . In our data, the proliferative subtype showed higher expression of the c-MET, ARID1A, CTNNB1, RAF1, MADM4, LGR5, GLUL1, TERT and CDC20 genes. Based on the expression profiles it is plausible that the RAS/MAPK, MET, AKT/mTOR and liver-Wnt pathways are responsible for cellular plasticity in the proliferative subtype [8, 44]. The normal-like subtype consisted of eight tumours and showed a scattered expression of hub genes (except TOP2A), drug-target genes and somatically mutated genes without any pattern, perhaps an indication that tissue analysed was mixed with normal tissue.
The higher expression of the hub genes TOP2A, CDK1 and BIRC5 in all samples indicates active cell division and proliferation. The selection of hub genes for the present study was based on their strong association with poor prognosis and overall survival in HCC patients [26, 27, 37, 40]. The strong correlation of TOP2A with CDK1 and BIRC5 as well as an earlier report showed that TOP2A could substitute other hub genes to be a good biomarker . Moreover, TOP2A is associated with tumour grades, HBV infection, and vascular invasion . In addition to expression of the hub genes, we noticed that GPC3 and IGF2 were expressed in 65–75% and AFP in more than 40% of the tumour samples. Our findings of GPC3 expression in 75% of HCC cases confirms earlier reports of its utility as a biomarker [42, 68]. GPC3 enhances cell proliferation through the Wnt/β-catenin pathway and is associated with a poor prognosis [42, 68]. IGF-2 expression is almost absent in adult hepatic cells, but its upregulation in HCC has been attributed to epigenetic mechanisms . GPC3 and IGF-2 are potential drug targets that have significantly reduced tumour growth and prolonged survival in Phase 1 clinical trials and in animal models respectively [70, 71]. Overall, we show that, in addition to HCC hub genes, GPC3 and IGF-2 may also serve as drug targets and early diagnostic markers for HCC.
The proliferative subtype in our study was characterised by high expression of the c-MET, ARID1A, CTNNB1, RAF1, MDM4, LGR5 and GLUL1 genes. Of these, c-MET along with hepatocyte growth factor (HGF) activates RAS-ERK and PI3K-AKT pathways, strengthening tumour aggressiveness with poor prognosis [72,73,74]. ARID1A has a context-dependent oncogenic and tumour suppressive function in HCC. A higher expression as seen in our study indicates oncogenic activity, supposedly mediated through cytochrome P450 and oxidative stress . RAF1 is a proto-oncogene that encodes MAP3 kinase and triggers cell proliferation through the ERK signalling pathway. Higher expression of RAF1 is associated with resistance to sorafenib, a tyrosine-kinase inhibitor and a standard drug for the treatment of advanced HCC . CTNNB1 encodes beta-catenin which is essential for the canonical Wnt pathway; higher expression of this gene in HCC is mediated through exonic mutations and epigenetic factors [77, 78]. LGR5 encodes a member of the G-protein coupled receptor superfamily and is believed to promote HCC metastasis formation by inducing epithelial-mesenchymal transition. GLUL1 catalyses the synthesis of glutamine from glutamate and ammonium in the liver tissue and is tightly controlled by Wnt/β-catenin signalling. Mutations in CTNNB1 activate the pathway leading to higher levels of GLUL1 and LGR5 [79, 80]. Overexpression of GLUL1 sensitizes HCCs to sorafenib, indicating the relevance of GLUL1 as a potential biomarker for the stratification of patients regarding treatment with sorafenib . The other two genes that showed expression in proliferative subtype tumours are TERT and CDC20. TERT is responsible for maintaining telomere length in tumours and re-activated in HCCs primarily through mutations in promoter region . Studies have shown a good concordance between mutations in CTNNB1 and the TERT promoter in HCCs, indicating that inhibitors targeting Wnt/β-catenin and TERT could be beneficial in HCC therapy . CDC20 regulates cell division through its interaction with anaphase-promoting complex/cyclosome (APC/C) and its overexpression is associated with poor prognosis .
Serum AFP, CA19-9, and CEA are used as preoperative tumour markers . Especially, AFP is commonly used in clinical practice to diagnose HCC and various tumours [86,87,88]. In combination with AFP, serum markers CA19-9 and CEA are being used to improve the diagnostic and prognostic performance of HCC patients. In this study, comparing serum AFP levels between the tumour subtypes showed that the normal-like subtype showed trend with higher AFP serum levels compared to the proliferative. However, it was statistically insignificant. This contrasts with published results, namely that the proliferative subtype has higher AFP levels . Interestingly, we observed a statistically significant positive correlation between AFP expression and AFP serum levels across all subtypes. Higher AFP levels promote HCC cell growth by activating the NF-κB pathway, in addition to suppressing the Fas/FADD-mediated apoptotic pathway [89, 90]. Higher AFP levels showed an association with higher metastasising activities and post-operative recurrence rates. In addition, we did not detect significant differences in serum levels of CA19-9 and CEA between proliferative and normal HCC subtypes.
Taken together our study has shown two main subtypes in BCLC-B classified tumours. We have demonstrated that a molecular classification of HCC can be achieved through a gene panel using RT-PCR. This approach enables patient stratification based on gene expression profiles for targeted personalized treatment. The genes c-MET, ARID1A, CTNNB1 and RAF1 showing an association with the proliferative subtype in our study may be used as molecular markers for subtype determination, and hub genes can be applied for HCC diagnosis. However, large prospective, well-designed follow-up studies are required to evaluate these marker genes for clinical applications. Besides being a retrospective study, with a rather small sample size our major limitation has been lacking longitudinal follow-ups of patients for a comprehensive survival analysis.
Availability of data and materials
The datasets generated and/or analysed during the current study are available in the Gene Expression Omnibus (GEO, http://www.ncbi.nlm.nih.gov/geo/), [GEO Submission (GSE191298)].
Hepatitis B virus
Hepatitis C virus
Real-time quantitative PCR
Barcelona Clinic Liver Cancer
Carbohydrate antigen 19-9
- CA 19-9:
Hepatocyte growth factor
Analysis of variance
Ferlay J, Colombet M, Soerjomataram I, Parkin DM, Pineros M, Znaor A, Bray F. Cancer statistics for the year 2020: an overview. Int J Cancer. 2021;149:778–89.
Yang JD, Hainaut P, Gores GJ, Amadou A, Plymoth A, Roberts LR. A global view of hepatocellular carcinoma: trends, risk, prevention and management. Nat Rev Gastroenterol Hepatol. 2019;16(10):589–604.
Park JW, Chen M, Colombo M, Roberts LR, Schwartz M, Chen PJ, Kudo M, Johnson P, Wagner S, Orsini LS, et al. Global patterns of hepatocellular carcinoma management from diagnosis to death: the BRIDGE Study. Liver Int. 2015;35(9):2155–66.
Nguyen VT. Hepatitis B infection in Vietnam: current issues and future challenges. Asia Pac J Public Health. 2012;24(2):361–73.
Park JW, Kim YJ, Kim DY, Bae SH, Paik SW, Lee YJ, Kim HY, Lee HC, Han SY, Cheong JY, et al. Sorafenib with or without concurrent transarterial chemoembolization in patients with advanced hepatocellular carcinoma: the phase III STAH trial. J Hepatol. 2019;70(4):684–91.
Llovet JM, Zucman-Rossi J, Pikarsky E, Sangro B, Schwartz M, Sherman M, Gores G. Hepatocellular carcinoma. Nat Rev Dis Primers. 2016;2:16018.
Llovet JM, Bru C, Bruix J. Prognosis of hepatocellular carcinoma: the BCLC staging classification. Semin Liver Dis. 1999;19(3):329–38.
Barcena-Varela M, Lujambio A. The endless sources of hepatocellular carcinoma heterogeneity. Cancers (Basel). 2021;13(11):2621.
Raja A, Haq F. Molecular classification of hepatocellular carcinoma: prognostic importance and clinical applications. J Cancer Res Clin Oncol. 2022;148(1):15–29.
Calderaro J, Ziol M, Paradis V, Zucman-Rossi J. Molecular and histological correlations in liver cancer. J Hepatol. 2019;71(3):616–30.
Rebouissou S, Nault JC. Advances in molecular classification and precision oncology in hepatocellular carcinoma. J Hepatol. 2020;72(2):215–29.
Lee JS, Heo J, Libbrecht L, Chu IS, Kaposi-Novak P, Calvisi DF, Mikaelyan A, Roberts LR, Demetris AJ, Sun Z, et al. A novel prognostic subtype of human hepatocellular carcinoma derived from hepatic progenitor cells. Nat Med. 2006;12(4):410–6.
Boyault S, Rickman DS, de Reynies A, Balabaud C, Rebouissou S, Jeannot E, Herault A, Saric J, Belghiti J, Franco D, et al. Transcriptome classification of HCC is related to gene alterations and to new therapeutic targets. Hepatology. 2007;45(1):42–52.
Chiang DY, Villanueva A, Hoshida Y, Peix J, Newell P, Minguez B, LeBlanc AC, Donovan DJ, Thung SN, Sole M, et al. Focal gains of VEGFA and molecular classification of hepatocellular carcinoma. Cancer Res. 2008;68(16):6779–88.
Hoshida Y, Nijman SM, Kobayashi M, Chan JA, Brunet JP, Chiang DY, Villanueva A, Newell P, Ikeda K, Hashimoto M, et al. Integrative transcriptome analysis reveals common molecular subclasses of human hepatocellular carcinoma. Cancer Res. 2009;69(18):7385–92.
Cancer Genome Atlas Research Network. Electronic address wbe, Cancer Genome Atlas Research N: Comprehensive and integrative genomic characterization of hepatocellular carcinoma. Cell. 2017;169(7):1327–41.
Kumar M, Zhao X, Wang XW. Molecular carcinogenesis of hepatocellular carcinoma and intrahepatic cholangiocarcinoma: one step closer to personalized medicine? Cell Biosci. 2011;1(1):5.
Thorgeirsson SS, Grisham JW. Molecular pathogenesis of human hepatocellular carcinoma. Nat Genet. 2002;31(4):339–46.
Hussain SP, Schwank J, Staib F, Wang XW, Harris CC. TP53 mutations and hepatocellular carcinoma: insights into the etiology and pathogenesis of liver cancer. Oncogene. 2007;26(15):2166–76.
Yeh SH, Chen PJ, Shau WY, Chen YW, Lee PH, Chen JT, Chen DS. Chromosomal allelic imbalance evolving from liver cirrhosis to hepatocellular carcinoma. Gastroenterology. 2001;121(3):699–709.
Dhanasekaran R, Nault JC, Roberts LR, Zucman-Rossi J. Genomic medicine and implications for hepatocellular carcinoma prevention and therapy. Gastroenterology. 2019;156(2):492–509.
Tornesello ML, Buonaguro L, Izzo F, Buonaguro FM. Molecular alterations in hepatocellular carcinoma associated with hepatitis B and hepatitis C infections. Oncotarget. 2016;7(18):25087–102.
Llovet JM, Montal R, Villanueva A. Randomized trials and endpoints in advanced HCC: role of PFS as a surrogate of survival. J Hepatol. 2019;70(6):1262–77.
Huang A, Yang XR, Chung WY, Dennison AR, Zhou J. Targeted therapy for hepatocellular carcinoma. Signal Transduct Target Ther. 2020;5(1):146.
Yoon JK, Choi JY, Rhee H, Park YN. MRI features of histologic subtypes of hepatocellular carcinoma: correlation with histologic, genetic, and molecular biologic classification. Eur Radiol. 2022;32(8):5119–33.
Chen Z, Chen J, Huang X, Wu Y, Huang K, Xu W, Xie L, Zhang X, Liu H. Identification of potential key genes for hepatitis B virus-associated hepatocellular carcinoma by bioinformatics analysis. J Comput Biol. 2019;26(5):485–94.
Li N, Li L, Chen Y. The identification of core gene expression signature in hepatocellular carcinoma. Oxid Med Cell Longev. 2018;2018:3478305.
Gao X, Wang X, Zhang S. Bioinformatics identification of crucial genes and pathways associated with hepatocellular carcinoma. Biosci Rep. 2018;38(6).
Li L, Lei Q, Zhang S, Kong L, Qin B. Screening and identification of key biomarkers in hepatocellular carcinoma: evidence from bioinformatic analysis. Oncol Rep. 2017;38(5):2607–18.
Wang Z, Chen S, Wang G, Li S, Qin X. CDCA3 is a novel prognostic biomarker associated with immune infiltration in hepatocellular carcinoma. Biomed Res Int. 2021;2021:6622437.
Feng Y, Tang X, Li C, Su Y, Wang X, Li N, Zhang A, Jiang F, Wu C. ARID1A is a prognostic biomarker and associated with immune infiltrates in hepatocellular carcinoma. Can J Gastroenterol Hepatol. 2022;2022:3163955.
Altundag O. Recent advances in systemic therapy for hepatocellular carcinoma. Exp Clin Transpl. 2022.
Cucarull B, Tutusaus A, Rider P, Hernaez-Alsina T, Cuno C, Garcia de Frutos P, Colell A, Mari M, Morales A. Hepatocellular carcinoma: molecular pathogenesis and therapeutic advances. Cancers (Basel). 2022;14(3):621.
Harding JJ, Zhu AX, Bauer TM, Choueiri TK, Drilon A, Voss MH, Fuchs CS, Abou-Alfa GK, Wijayawardana SR, Wang XA, et al. A phase Ib/II study of ramucirumab in combination with emibetuzumab in patients with advanced cancer. Clin Cancer Res. 2019;25(17):5202–11.
Shi YX, Zhu T, Zou T, Zhuo W, Chen YX, Huang MS, Zheng W, Wang CJ, Li X, Mao XY, et al. Prognostic and predictive values of CDK1 and MAD2L1 in lung adenocarcinoma. Oncotarget. 2016;7(51):85235–43.
Liang B, Li C, Zhao J. Identification of key pathways and genes in colorectal cancer using bioinformatics analysis. Med Oncol. 2016;33(10):111.
Zhu Q, Sun Y, Zhou Q, He Q, Qian H. Identification of key genes and pathways by bioinformatics analysis with TCGA RNA sequencing data in hepatocellular carcinoma. Mol Clin Oncol. 2018;9(6):597–606.
Diaz-Gonzalez A, Reig M, Bruix J. Treatment of hepatocellular carcinoma. Dig Dis. 2016;34(5):597–602.
Cillo U, Vitale A, Grigoletto F, Farinati F, Brolese A, Zanus G, Neri D, Boccagni P, Srsen N, D’Amico F, et al. Prospective validation of the Barcelona Clinic Liver Cancer staging system. J Hepatol. 2006;44(4):723–31.
Wei Z, Liu Y, Qiao S, Li X, Li Q, Zhao J, Hu J, Wei Z, Shan A, Sun X, et al. Identification of the potential therapeutic target gene UBE2C in human hepatocellular carcinoma: an investigation based on GEO and TCGA databases. Oncol Lett. 2019;17(6):5409–18.
Yang H, Zhang X, Cai XY, Wen DY, Ye ZH, Liang L, Zhang L, Wang HL, Chen G, Feng ZB. From big data to diagnosis and prognosis: gene expression signatures in liver hepatocellular carcinoma. PeerJ. 2017;5:e3089.
Guo M, Zhang H, Zheng J, Liu Y. Glypican-3: A New Target for Diagnosis and Treatment of Hepatocellular Carcinoma. J Cancer. 2020;11(8):2008–21.
Neureiter D, Stintzing S, Kiesslich T, Ocker M. Hepatocellular carcinoma: therapeutic advances in signaling, epigenetic and immune targets. World J Gastroenterol. 2019;25(25):3136–50.
Zucman-Rossi J, Villanueva A, Nault JC, Llovet JM. Genetic Landscape and biomarkers of hepatocellular carcinoma. Gastroenterology. 2015;149(5):1226–39.
de Bono J, Lin CC, Chen LT, Corral J, Michalarea V, Rihawi K, Ong M, Lee JH, Hsu CH, Yang JC, et al. Two first-in-human studies of xentuzumab, a humanised insulin-like growth factor (IGF)-neutralising antibody, in patients with advanced solid tumours. Br J Cancer. 2020;122(9):1324–32.
Wang H, Rao B, Lou J, Li J, Liu Z, Li A, Cui G, Ren Z, Yu Z. The function of the HGF/c-Met axis in hepatocellular carcinoma. Front Cell Dev Biol. 2020;8:55.
Okoye-Okafor UC, Bartholdy B, Cartier J, Gao EN, Pietrak B, Rendina AR, Rominger C, Quinn C, Smallwood A, Wiggall KJ, et al. New IDH1 mutant inhibitors for treatment of acute myeloid leukemia. Nat Chem Biol. 2015;11(11):878–86.
Hird AW, Tron AE. Recent advances in the development of Mcl-1 inhibitors for cancer therapy. Pharmacol Ther. 2019;198:59–67.
Jochemsen AG. Reactivation of p53 as therapeutic intervention for malignant melanoma. Curr Opin Oncol. 2014;26(1):114–9.
Adnane L, Trail PA, Taylor I, Wilhelm SM. Sorafenib (BAY 43-9006, Nexavar), a dual-action inhibitor that targets RAF/MEK/ERK pathway in tumor cells and tyrosine kinases VEGFR/PDGFR in tumor vasculature. Methods Enzymol. 2006;407:597–612.
Liu L, Qin S, Zheng Y, Han L, Zhang M, Luo N, Liu Z, Gu N, Gu X, Yin X. Molecular targeting of VEGF/VEGFR signaling by the anti-VEGF monoclonal antibody BD0801 inhibits the growth and induces apoptosis of human hepatocellular carcinoma cells in vitro and in vivo. Cancer Biol Ther. 2017;18(3):166–76.
An P, Xu J, Yu Y, Winkler CA. Host and viral genetic variation in HBV-related hepatocellular carcinoma. Front Genet. 2018;9:261.
Huang W, Skanderup AJ, Lee CG. Advances in genomic hepatocellular carcinoma research. Gigascience. 2018;7(11):giy135.
Ruden M, Puri N. Novel anticancer therapeutics targeting telomerase. Cancer Treat Rev. 2013;39(5):444–56.
Pez F, Lopez A, Kim M, Wands JR, Caron de Fromentel C, Merle P. Wnt signaling and hepatocarcinogenesis: molecular targets for the development of innovative anticancer drugs. J Hepatol. 2013;59(5):1107–17.
Niu ZS, Niu XJ, Wang WH. Genetic alterations in hepatocellular carcinoma: an update. World J Gastroenterol. 2016;22(41):9069–95.
Schulze K, Nault JC, Villanueva A. Genetic profiling of hepatocellular carcinoma using next-generation sequencing. J Hepatol. 2016;65(5):1031–42.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method. Methods. 2001;25(4):402–8.
de Hoon MJ, Imoto S, Nolan J, Miyano S. Open source clustering software. Bioinformatics. 2004;20(9):1453–4.
Saldanha AJ. Java Treeview—extensible visualization of microarray data. Bioinformatics. 2004;20(17):3246–8.
Muller L, Hahn F, Mahringer-Kunz A, Stoehr F, Gairing SJ, Foerster F, Weinmann A, Galle PR, Mittler J, Pinto Dos Santos D, et al. Immunonutritive scoring in patients with hepatocellular carcinoma undergoing transarterial chemoembolization: Prognostic nutritional index or controlling nutritional status score? Front Oncol. 2021;11:696183.
Golfieri R, Bargellini I, Spreafico C, Trevisani F. Patients with Barcelona Clinic Liver Cancer stages B and C hepatocellular carcinoma: time for a subclassification. Liver Cancer. 2019;8(2):78–91.
Bolondi L, Burroughs A, Dufour JF, Galle PR, Mazzaferro V, Piscaglia F, Raoul JL, Sangro B. Heterogeneity of patients with intermediate (BCLC B) Hepatocellular Carcinoma: proposal for a subclassification to facilitate treatment decisions. Semin Liver Dis. 2012;32(4):348–59.
Ha Y, Shim JH, Kim SO, Kim KM, Lim YS, Lee HC. Clinical appraisal of the recently proposed Barcelona Clinic Liver Cancer stage B subclassification by survival analysis. J Gastroenterol Hepatol. 2014;29(4):787–93.
Weinmann A, Koch S, Sprinzl M, Kloeckner R, Schulze-Bergkamen H, Duber C, Lang H, Otto G, Worns MA, Galle PR. Survival analysis of proposed BCLC-B subgroups in hepatocellular carcinoma patients. Liver Int. 2015;35(2):591–600.
Nault JC, Villanueva A. Biomarkers for hepatobiliary cancers. Electronic 1527–3350.
Desert R, Nieto N, Musso O. Dimensions of hepatocellular carcinoma phenotypic diversity. World J Gastroenterol. 2018;24(40):4536–47.
Nishida T, Kataoka H. Glypican 3-targeted therapy in hepatocellular carcinoma. Cancers (Basel). 2019;11(9):1339.
Thomas H. Liver cancer: IGF2—an epigenetic oncodriver in HCC. Nat Rev Gastroenterol Hepatol. 2016;13(11):625.
Shi D, Shi Y, Kaseb AO, Qi X, Zhang Y, Chi J, Lu Q, Gao H, Jiang H, Wang H, et al. Chimeric antigen receptor-glypican-3 T-cell therapy for advanced hepatocellular carcinoma: results of phase I trials. Clin Cancer Res. 2020;26(15):3979–89.
Martinez-Quetglas I, Pinyol R, Dauch D, Torrecilla S, Tovar V, Moeini A, Alsinet C, Portela A, Rodriguez-Carunchio L, Sole M, et al. IGF2 is up-regulated by epigenetic mechanisms in hepatocellular carcinomas and is an actionable oncogene product in experimental models. Gastroenterology. 2016;151(6):1192–205.
Garcia-Vilas JA, Medina MA. Updates on the hepatocyte growth factor/c-Met axis in hepatocellular carcinoma and its therapeutic implications. World J Gastroenterol. 2018;24(33):3695–708.
Organ SL, Tsao MS. An overview of the c-MET signaling pathway. Ther Adv Med Oncol. 2011;3(1 Suppl):S7–19.
Fujita R, Blot V, Wong E, Stewart C, Lieuw V, Richardson R, Banah A, Villicana J, Timmer A, Coronella J, et al. A novel non-agonist c-Met antibody drug conjugate with superior potency over a c-Met tyrosine kinase inhibitor in c-Met amplified and non-amplified cancers. Cancer Biol Ther. 2020;21(6):549–59.
Sun X, Wang SC, Wei Y, Luo X, Jia Y, Li L, Gopal P, Zhu M, Nassour I, Chuang JC, et al. Arid1a has context-dependent oncogenic and tumor suppressor functions in liver cancer. Cancer Cell. 2017;32(5):574–89.
Tian N, Wu D, Tang M, Sun H, Ji Y, Huang C, Chen L, Chen G, Zeng M. RAF1 expression is correlated with HAF, a parameter of liver computed tomographic perfusion, and may predict the early therapeutic response to sorafenib in advanced hepatocellular carcinoma patients. Open Med (Wars). 2020;15:167–74.
Khalaf AM, Fuentes D, Morshid AI, Burke MR, Kaseb AO, Hassan M, Hazle JD, Elsayes KM. Role of Wnt/beta-catenin signaling in hepatocellular carcinoma, pathogenesis, and clinical significance. J Hepatocell Carcinoma. 2018;5:61–73.
Javanmard D, Najafi M, Babaei MR, Karbalaie Niya MH, Esghaei M, Panahi M, Safarnezhad Tameshkel F, Tavakoli A, Jazayeri SM, Ghaffari H, et al. Investigation of CTNNB1 gene mutations and expression in hepatocellular carcinoma and cirrhosis in association with hepatitis B virus infection. Infect Agent Cancer. 2020;15:37.
Lachenmayer A, Alsinet C, Savic R, Cabellos L, Toffanin S, Hoshida Y, Villanueva A, Minguez B, Newell P, Tsai HW, et al. Wnt-pathway activation in two molecular classes of hepatocellular carcinoma and experimental modulation by sorafenib. Clin Cancer Res. 2012;18(18):4997–5007.
Desert R, Rohart F, Canal F, Sicard M, Desille M, Renaud S, Turlin B, Bellaud P, Perret C, Clement B, et al. Human hepatocellular carcinomas with a periportal phenotype have the lowest potential for early recurrence after curative resection. Hepatology. 2017;66(5):1502–18.
Sohn BH, Park IY, Shin JH, Yim SY, Lee JS. Glutamine synthetase mediates sorafenib sensitivity in beta-catenin-active hepatocellular carcinoma cells. Exp Mol Med. 2018;50(1):e421.
Jager K, Walter M. Therapeutic Targeting of Telomerase. Genes (Basel). 2016;7(7):39.
Kotiyal S, Evason KJ. Exploring the interplay of telomerase reverse transcriptase and beta-catenin in hepatocellular carcinoma. Cancers (Basel). 2021;13(16):4202.
Li J, Gao JZ, Du JL, Huang ZX, Wei LX. Increased CDC20 expression is associated with development and progression of hepatocellular carcinoma. Int J Oncol. 2014;45(4):1547–55.
Edoo MIA, Chutturghoon VK, Wusu-Ansah GK, Zhu H, Zhen TY, Xie HY, Zheng SS. Serum biomarkers AFP, CEA and CA19-9 combined detection for early diagnosis of hepatocellular carcinoma. Iran J Public Health. 2019;48(2):314–22.
Montal R, Andreu-Oller C, Bassaganyas L, Esteban-Fabro R, Moran S, Montironi C, Moeini A, Pinyol R, Peix J, Cabellos L, et al. Molecular portrait of high alpha-fetoprotein in hepatocellular carcinoma: implications for biomarker-driven clinical trials. Br J Cancer. 2019;121(4):340–3.
Ahn JC, Teng PC, Chen PJ, Posadas E, Tseng HR, Lu SC, Yang JD. Detection of circulating tumor cells and their implications as a biomarker for diagnosis, prognostication, and therapeutic monitoring in hepatocellular carcinoma. Hepatology. 2021;73(1):422–36.
Su F, Weiss NS, Beste LA, Moon AM, Jin GY, Green P, Berry K, Ioannou GN. Screening is associated with a lower risk of hepatocellular carcinoma-related mortality in patients with chronic hepatitis B. J Hepatol. 2021;74(4):850–9.
Chen T, Dai X, Dai J, Ding C, Zhang Z, Lin Z, Hu J, Lu M, Wang Z, Qi Y, et al. AFP promotes HCC progression by suppressing the HuR-mediated Fas/FADD apoptotic pathway. Cell Death Dis. 2020;11(10):822.
Li QT, Qiu MJ, Yang SL, Fang X, He XX, Wang MM, Li YN, Xiong ZF, Huang S. Alpha-fetoprotein regulates the expression of immune-related proteins through the NF-kappaB (P65) pathway in hepatocellular carcinoma cells. J Oncol. 2020;2020:9327512.
We thank all study subjects for their participation. The authors acknowledge Mr. Bui Dinh Tung, and Ms. Bui Thuy Linh for their experimental support.
Open Access funding enabled and organized by Projekt DEAL. NXH acknowledges a Grant from the Vietnam National Foundation for Science and Technology Development (NAFOSTED, Grant Number: 108.02-2018.315). TPV acknowledges Grants from the German Federal Ministry of Education and Research (BMBF01DP19006A, BMBF01DP17047).
Ethics approval and consent to participate
The study was approved by the institutional review board of the 108 Military Central Hospital and the 103 Military Hospital of the Vietnam Military Medical University, Hanoi, Vietnam. All experiments were performed in accordance with relevant guidelines and regulations. Informed written consent was obtained from the participants in the study.
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.
About this article
Cite this article
Pallerla, S.R., Hoan, N.X., Rachakonda, S. et al. Custom gene expression panel for evaluation of potential molecular markers in hepatocellular carcinoma. BMC Med Genomics 15, 235 (2022). https://doi.org/10.1186/s12920-022-01386-7
- Hepatocellular carcinoma
- Hepatitis B virus