Skip to main content

CTNNB1 mutations, TERT polymorphism and CD8+ cell densities in resected hepatocellular carcinoma are associated with longer time to recurrence



Hepatocellular carcinoma (HCC) is a fatal disease characterized by early genetic alterations in telomerase reverse transcriptase promoter (TERTp) and β-catenin (CTNNB1) genes and immune cell activation in the tumor microenvironment. As a novel approach, we wanted to assess patient survival influenced by combined presence of mutations and densities of CD8+ cytotoxic T cells.


Tissue samples were obtained from 67 HCC patients who had undergone resection. We analysed CD8+ T cells density, TERTp mutations, rs2853669 polymorphism, and CTNNB1 mutations. These variables were evaluated for time to recurrence (TTR) and disease free survival (DFS).


TERTp mutations were found in 75.8% and CTNNB1 mutations in 35.6% of the patients. TERTp mutations were not associated with survival but polymorphism rs2853669 in TERTp was associated with improved TTR and DFS. CTNNB1 mutations were associated with improving TTR. High density of CD8+ T-lymphocytes in tumor center and invasive margin correlated with longer TTR and DFS. Combined genetic and immune factors further improved survival showing higher predictive values. E.g., combining CTNNB1 mutations and high density of CD8+ T-lymphocytes in tumor center yielded HRs of 0.12 (0.03–0.52), p = 0.005 for TTR and 0.25 (0.09–0.74), p = 0.01 for DFS.


The results outline a novel integrative approach for prognostication through combining independent predictive factors from genetic and immune cell profiles. However, larger studies are needed to explore multiple cell types in the tumor microenvironment.

Peer Review reports


Primary livers cancer are the sixth most frequent diagnosed tumors globally and the third cause of cancer related death worldwide. Hepatocellular carcinoma (HCC) encompasses up to 85% of all primary liver cancers [1]. It develops as a result of a chronic liver disease which stepwise progresses from chronic inflammation to cirrhosis, fibrosis and HCC [2]. Malignant transformation of hepatocytes is a complex molecular pro-cess comprising numerous genetic and epigenetic alterations. These changes increase genetic diversity and accelerate evolution of tumor cells. However, increasing diversity stimulates immune system to action against cells presenting neoepitopes [3]. CD8+ cytotoxic T cells play a pivotal role in antitumor immune response. They are also one of many tumor-infiltrating lymphocytes (TIL), which are vital in predicting clinical outcome (overall survival (OS) and disease free-survival (DFS)) in primary and metastatic liver cancer. Evidence shows that high densities of CD8+ cells in HCC are associated with better OS and DFS [4]. These at-tributes point to the role of intratumor CD8+ cells as a potential prognostic biomarker [4, 5].

In cancer cells, reactivation of telomerase prevents telomere shortening and may lead to unlimited cellular proliferation, enabling subsequent transformation. Almost 90% of human cancers show reactivation of telomerase [6]. Telomerase reverse transcriptase (TERT) is the catalytic subunit of telomerase, which additionally contains the telomerase RNA component (TERC) [7]. In HCC, TERT promoter (TERTp) mutations are the most frequent and earliest somatic alterations with a prevalence of about 60% [8, 9]. Until to-date, two recurrent somatic mutations in TERTp (− 124 and − 146 upstream of the ATG start site) have been associated with many types of cancer, including HCC. They create de novo binding sites for the ETS family of transcription factors and increase TERT expression [10]. TERT single nucleotide polymorphism (SNP) rs2853669 disrupts an existing ETS2 binding site, which results in decreased telomerase activity [11]. It was shown that SNP rs2853669 correlates with patient survival with variable effects in various cancers [12,13,14]. On the other hand, the same SNP was associated with poor prognosis in HCC [15]. Hot-spot mutations present in exon 3 of the β-catenin (CTNNB1) gene are critical for hepatocarcinogenesis and are among the most frequent alterations (20–40% of cases) [16]. They predominantly occur at the phosphorylation sites, induce accumulation of nuclear β-catenin and lead to activation of the Wnt signaling pathway [17]. Moreover, mutations in exon 3 of CTNNB1tend to be associated with favourable prognosis in HCC [18, 19].

Research in the past decade has improved our understanding of hepatocarcinogenesis highlighting the multitude of genetic changes in the tumor but also the repertoire of interacting cells in the tumor microenvironments [20,21,22]. The aim of present study is to examine the interplay of most frequent driver mutations and one well-known immune marker. As a novel hypothesis, we posit that these two environments may act independently in influencing survival but when combined they may act additively supporting each other’s action. We chose to examine the interplay between genetic factors (TERTp and CTNNB1) in the tumor and immune factors (cytotoxic CD8+ cells) in the tumor microenvironment and their prognostic value in terms of time to recurrence (TTR), DFS and OS among resected HCC patients. Although the present patient series of 67 patients is relatively small it was homogenous as to the treatment by resection in a single center and precise follow-up thus providing a valuable addition to HCC patients.


Patient selection and characteristic

The study cohort consisted of 67 HCC patients (16 women and 51 men) of median age 69 years (range 24–86) who underwent primary tumor resection in Pilsen University hospital between 1997 and 2019. Patients were not subjected to neo-adjuvant therapy before operation, nor had distant metastasis. Clinical characteristics of the enrolled patients are presented in Table 1. The study was conducted in accordance with the Declaration of Helsinki and was approved by the Ethical Board of the Faculty of Medicine and University hospital in Pilsen (118/2021, March 11 2021).

Table 1 Clinical characteristics of recruited patients

DNA extraction, sanger sequencing

DNA was extracted from FFPE (formaldehyde fixed paraffin embedded) tumor tissue with RecoverAll™ Total Nucleic Acid Isolation Kit for FFPE according to the manufacturer’s protocol. Isolated DNA was used to obtain TERTp and CTNNB1 exon 3 mutation profiles. Mutations were analyzed in tumor and non-tumor adjacent tissue. PCR was carried out in 25-μl volume containing 40 ng DNA,0.1 mM dNTP, 2.5 mM Mg2+ and 0.11 μM of each primer and 1 unit Taq polymerase (Taq DNA polymerase, Top-Bio) and 5% glycerol for TERTp. For CTNNB1, 20 ng of DNA, 0.2 mM dNTP, 2.5 mM Mg2+ and 0.2 μM of each primer and 1 unit Taq polymerase was used. TERTp and CTNNB1 amplification PCR was performed under following conditions: 95 °C for 5 min/3 min, next 33 cycles of 95 °C for 30 sec, 68 °C/60 °C for 1 min, 72 °C for 45 sec respectively. Primer pairs were described previously [23, 24]. Sanger sequencing was performed in both directions. Chromatograms were analyzed by visual comparison and mutation was called when the height of the over-lapping peak encompassed at least 33% of the primary peak.

Polymorphism rs2853669 genotyping

SNP rs2853669 at TERTp was genotyped with TaqMan® SNP Genotyping Assays (Applied Bioscience) and 20 ng of isolated DNA on CFX96 Touch Deep Well Real-Time PCR Detection System (Bio-Rad). Genotyping conditions were as followed: 95 °C for 10 min, next 50 cycles of 95 °C for 15 sec and 60 °C for 1 min.

Histology and immunohistochemistry

For each patient FFPE tissue blocks containing center of the tumor (TC) and invasive margin (IM) and, when available, separate block with non-tumor liver, 2- to 3-cm distant from the tumor site, were retrieved from the pathology archive; the detailed procedure are described elsewhere [25]. One or two tis-sue sections of 4-μm thickness were cut and mounted onto BOND Plus microscope slides (K8020, Leica Bio-systems Newcastle Ltd., United Kingdom). Immunoperoxidase detection of CD8+ T cells was performed using fully automated BOND-III IHC/ISH stainer with anti-CD8 monoclonal primary antibodies (clone 4B11) from Leica Biosystems (Newcastle Ltd., United Kingdom) and horseradish peroxidase (HRP)-linker antibody conjugate system (Bond™ Polymer Refine Detection). Sections were counterstained with Mayer’s hema-toxylin and embedded into Micromount mounting medium (Leica). Appropriate positive (tonsils) and negative tissue control samples were used. All sections were examined under Olympus CX41 microscope (Olympus, Tokyo, Japan) by two pathologists (AT and LC). Sections were microanatomically divided into TC and IM and non-tumor liver. The IM was defined as a 1000 μm-region centered on the border separating the malignant cell nests from the host tissue (Hendry et al. 2017), whereas TC represents the remaining tumor area [26]. Eight fields of view were selected by systematic uniform random sampling from the TC and IM using the objective 20×. Pictures were captured by the PromiCam 3-3CP digital camera (Promicra, Prague, Czech Republic), coupled with the QuickPhoto Industrial 3.2 software (Promicra, Prague, Czech Republic). The stereological analysis was performed by AT, who was blinded to the clinical outcome, using the computer assisted stereology software Ellipse (ViDiTo, Kosice, Slovak Republic). CD8+ nucleated cell profiles were counted using a probe consisting of a set of 2D unbiased counting frames (UCF). The number of CD8+ profiles per mm2 of area (QA = numerical density) was calculated as the number of positive lymphocytic profiles divided by the total, investigated counting frame area. To eliminate skewness in the distribution, the raw density of CD8+ nucleated cell profiles CD8+ (further denoted as density) was converted into corresponding percentile values and categorized into Low (below 25th percentile) vs Intermediate-High (25th–100th percentile) (further denoted as high).

Statistical analysis

Statistical analysis was performed in the R environment (v.4.1.1). Spearman correlation was used to evaluate relations between variables. Survival analysis was performed with the Finalfit package [27]. Influence of individual and combined variables on TTR, DFS, OS was assessed by univariate and multivariable Cox regression. TTR was defined as period from tumor resection to recurrence/metastasis diagnosis. DFS was described as time from tumor resection until recurrence/metastasis detection or death due to any cause. OS was determined as the time from tumor resection to death due to any cause. Kaplan-Meier analysis was performed with the survival package and plots were generated with the survminer package [28, 29]. Differences between groups were assessed with Log-rank test. CD8+ densities in TC and IM were compared with the Wilcoxon test. Results were considered statistically significant when p-value < 0.05.


TERT promoter and CTNNB1 mutations and polymorphism rs2853669

We sequenced 67 HCC cases for TERTp and CTNNB1 mutation analysis. Good quality sequences were obtained for 58 and 59 cases, respectively. We detected TERTp mutation in 44 (75.8%) patients. TERTp mutations were located in two mutually exclusive hotspots − 124 bp in 41 and − 146 bp in three cases, all of them were heterozygous changes. CTNNB1 mutations were identified in 21 patients (35.6%): 20 heterozygous and 1 homozygous. Of them 16 patients had one missense mutation, one had a missense mutation accompanied by a synonymous change, one had three missense mutations and two has only a synonymous change (Supplementary Table 1). Mutational changes in TERTp and CTNNB1 were not observed in adjacent non-tumor tissue. Additionally, all patients were successfully genotyped for polymorphism rs2853669. Overall frequencies for AA, AG and GG genotypes were 44.8, 32.8 and 22.4%, respectively. For most of our genetic and clinical data, we did not find any correlations. However, correlations were found between − 146 bp TERTp mutation and concentration of alfa-fetoprotein (p = 0.004), and between CTNNB1 mutations and diabetes mellitus (p = 0.004) and extent of encapsulation (p = 0.03) (Supplementary Table 2). We did not detect any correlations between mutations and the polymorphism. Clinical data were discussed in more detail in our previous work (25); only TNM stage (HR = 1.5. p = 0.032) and age (HR = 0.96, p = 0.007) had impact on patients TTR.

CD8+ distribution

IM displayed significantly higher density of CD8+ cells than TC (p < 0.0001). CD8+ cell densities in non-tumor liver were higher compared with TC (p < 0.0001), and lower when compared with IM. (Fig. 1., Table 2). Significant correlations between clinical data and CD8+ cells are presented in Supplement (Supplementary Table 3).

Fig. 1
figure 1

Immunoperoxidase staining for CD8+ lymphocytes in HCC. Regions of interest are denoted: TC (tumor center) and IM (invasive margin). The IM was defined as a 1-mm region centered on the border separating the malignant cell nests from the host tissue. The TC represented the remaining tumor area. Scale bar 1000 μm

Table 2 Density of CD8+ cell (number of nucleated cell profiles/mm2)

Prognostic values of CTNNB1, TERTp mutation, rs2853669 and CD8+ cell densities

We performed univariate Cox regression and Kaplan-Meier analysis to investigate the association between genetic, immune factors and TTR, DFS, OS. To confirm differences between groups we performed Log-rank test (Supplementary Table 4). None of the variables were associated with OS. We did not detect any associations between TERTp mutations and TTR or DFS. Presence of a CTNNB1 mutation was associated with longer TTR but did not affect DFS (Table 3, Fig. 2). Presence of homozygote GG in rs2853669 polymorphism was associated with longer TTR and longer DFS, whereas heterozygotes were not associated (Table 3, Supplementary Fig. 1). High CD8+ cell densities in TC and IM were associated with lower risk of recurrence and longer DFS. (Table 3, Supplementary Fig. 2).

Table 3 Univariable analysis of genetic and immune variables associated with TTR and DFS
Fig. 2
figure 2

Kaplan-Meier analysis for TTR (A) and DFS (B) according to presence or absence of CTNNB1 mutation

Combined genetic and immune variables and their prognostic value

Next, we analyzed associations of combined genetic and immune factors with TTR, DFS and OS. The variables were grouped as follows: presence or absence CTNNB1 mutation, presence of AA, AG, or GG genotype for rs2853669 with low or high density of CD8+ cells. Altought TERTp mutations were not statistically significant in our individual analysis, due to previous reports [15, 30] we decided to combine them with polymorphism rs2853669 and evaluate their impact on patient’s outcomes. Individual analyses of CTNNB1 mutation and CD8+ cells density showed that the presence of mutation and high cell densities (HD) could be considered a favorable prognostic factor. Thus, we decided to analyze their combined influence on TTR and DFS. Cases were divided in 3 groups: presence of two good prognostic factors (CTNNB1(+)/CD8 + (HD)), presence of only one (CTNNB1(+/−)/CD8(LD/HD)) and presence of adverse prognostic factor (CTNNB1(−)/CD8(LD)). Presence of mutation in CTNNB1 and high CD8+ cells densities in TC and IM were associated with lower risk of recurrence (HR = 0.12 p = 0.005 and HR = 0.17 p = 0.012 respectively) and longer disease-free survival (HR = 0.25, p = 0.005 and HR = 0.34 p = 0.025 respectively) (Table 4 and Fig. 3).

Table 4 Analysis for TTR and DFS for combined genetic and immune variables
Fig. 3
figure 3

Kaplan-Meier analysis for TTR according to presence/absence of the CTNNB1mutation and high/low CD8+ densities in TC (A) and IM (C) and DFS in TC (B) and IM (D)

Similar to mutation in CTNNB1, the GG genotype in SNP rs2853669 can also be considered a “good” prognostic factor (Table 3). Therefore, we combined it with CD8+ cell densities in TC and IM. We observed statistically significant association between combination of density CD8+ with the GG genotype with TTR in TC and IM (HR = 0.10, p = 0.013 and HR = 0.17, p = 0.015, respectively) and with DFS in TC and IM (HR = 0.19, p = 0.011 and HR = 0.22, p = 0.003, respectively). Additionally, combination of CD8+ cell densities and the AG heterozygote in IM was associated with longer DFS survival (HR =0.36, p = 0.037). These findings were con-firmed in Kaplan-Meier analysis. (Fig. 3 and Supplementary Fig. 3) Combination of TERTp mutation (− 124 bp position) and polymorphism rs2853669 (where GG genotype can be considered as “good” prognostic factor) revealed that there was statistically significant association between GG genotype and lack of TERTp mutation and longer TTR (HR = 0.08, p = 0.024) and longer DFS survival (HR = 0.23, p = 0.024). Similar pattern was observed for presence of TERTp mutation and GG genotype, showing longer TTR and longer DFS (0.15, p = 0.009 and 0.20, p = 0.007, respectively).

Multivariable analysis for individual and combined variables

Multivariable Cox regression models for each individual and combined parameters adjusted for TNM staging and age were built to assess their independent ability to predict TTR (Supplementary Table 5 and 6). Among the individual parameters, CTNNB1 mutation and rs2853669 were independently associated with TTR. Immune markers (CD8+ T cell densities in TC and IM) lost predictive capabilities after adjustment for TNM stage and age. Among combined parameters, combinations between rs2853669 and cell densities in TC and IM were associated with TTR. Moreover, combination of TERTp mutation and rs2853669 retain their predictive value after adjustment for TNM stage and age and was associated with longer TTR.


The novelty of the present study was that for the first time results were presented on the interplay be-tween genetic (CTNNB1 mutation and TERTp SNP rs2853669) and immune (CD8+ cells density) factors, their prognostic value and present usefulness of TTR in HCC prognostication. The GG genotype of SNPrs2853669, mutation in CTNNB1 and high CD8+ cell density were individually associated with improved TTR and DFS. Presence of CTNNB1 mutation and high CD8+ cell densities in both region of interest (TC or IM) strongly predicted longer TTR and DFS. Similarly, high CD8 + cell densities in TC or IM and presence of GG genotype were also associated with longer TTR and DFS. Combining positive genetic and immune factors additively improved each other’s influence.

TERTp mutations are a well-known cause of telomerase reactivation in many cancers [31]. We confirmed high frequency of TERTp mutation in the present HCC, not associated with viral hepatitis patients (75.8%). Our mutation frequency is somewhat higher than has been reported for another European dataset (56.6%) for no-viral HCC [32]. In concordance with available data, the nucleotide change in − 124 bp was the most common TERTp mutation (93.2%). TERTp mutation and polymorphism rs2853669 did not show any correlation with clinical and pathological data in our analysis. Neither did we observe any associations with TTR, DFS or OS for the TERTp mutation; the relation with prognosis is rather controversial, some literatureshowing association of TERTp mutation with poor DSF and OS [33] other does not [34, 35]. Our data showed a statistically significant association of GG genotype of SNP rs2853669 with longer TTR (HR = 0.29, p = 0.02) and longer DFS (HR = 0.42, p = 0.025). The associations were not significant for the AG heterozygotes. Presence of G allele modulated effect of TERTp mutations because the SNP destroys the pre-existing binding site for transcription factors EST [36]. Our data on the association of this SNP with favourable patient outcome are in agreement with prognostic data for various cancers [37]. Nevertheless, even deviant results have been reported for unknown reasons in an Asian population [15]. Previous reports which evaluate impact of TERTp mutation and rs2853669 and their combined effect on patient survival are not in agreement showing no correlation with survival or correlation with poor prognosis [15, 30]. Our data show opposite pattern where combination of GG allele and lack/presence of TERTp mutation is significantly correlated with longer TTR and longer DFS. We speculate firstly, that our patient etiology differed from the previous reports on viral etiology. Secondly, our data show a lesser role for a TERTp mutational status compared to higher predictive values for the rs2885669 GG allele.

The CTNNB1 gene is a key partner of the Wnt/β-catenin signaling pathway, which is activated in 30–50% of HCC cases [2]. This gene is one of the most frequently mutated cancer driver genes in HCC [37]. CTNNB1mutations are missense mutations, most of which are located at the exon 3 phosphorylation sites and they increase the expression of β-catenin protein in the cytoplasm and nucleus [17, 38]. According to the literature, nuclear expression is associated with improved prognosis, opposite to cytoplasmic and membranous expression, which are associates with poorer prognosis [39, 40]. CTNNB1mutations have been identified preferentially in early stages and HCV-related HCC [41, 42]. Our observed frequency of CTNNB1 mutations (35.6%) and the affected nucleotide sequences (Supplement Table 1) agreed with the literature (17). Available data show associations of CTNNB1 mutation with favorable prognosis (especially with longer OS) [19]. Our data showed no association with OS or DFS, but we observe longer TTR (HR =0.31, p = 0.02). Moreover, we found no association between CTNNB1 mutations and clinical or pathological data, which may be related to the non-viral etiology of the present patients compared to the majority of data originating from Asian populations with viral etiology.

TILs (T cells, B cells, natural killer (NK) cells) represent one of the elements of host immune response to tumor [5]. They may manifest pro- or anti-tumor characteristics dependent on lymphocyte subsets and their phenotypic orientation [43]. The role of TILs has been extensively explored in HCC [4, 5, 43, 44]. These results show that patients with high numbers of NK cells, B-cells, T-lymphocytes and dendritic cells had a favorable prognosis while poor prognosis was associated with high numbers of monocytes, neutrophils and Tregs [44]. Role of high CD8+ cell densities has been explored in many cancers, including HCC, usually showing associations with improved prognosis [4, 45,46,47]. In HCC, CD8+ cells were among the prognostic factors considered in the meta-analysis by Xu et al., showing associations of CD8+ densities in tumor and its margin with higher OS and DFS [4]. As high density of CD8+ cells predicts a favorable OS and DFS they may be used as biomarker of patient prognosis. On par with available data, our CD8+ analysis also showed high predictive value for DFS (HR = 0.41, p = 0.009), but not OS [48]. In our research, we also observed influence of high density of CD8+ cells on longer TTR (HR = 0.32, p = 0.005).

Many analyses have been published showing association with DFS or OS, but we did not find any data considering TTR. We believe that among the survival parameters analysed by us, TTR is highly relevant for prognostication. TTR is a tumor-specific outcome in contrast to less specific endpoints DFS and OS.

Previous prognostic studies have either used genetic or immune cell markers, but we are not aware of previous studies combining the two. There is a scientific rationale for the combination, because mutations are present in the tumor and immune response is mounted by the tumor microenvironment. These may be viewed as independent predictors of outcome although they interact, and mutation frequency may generally stimulate immune response [49, 50]. We decided to explore combined effect of CTNNB1 mutation and CD8+ cells in TC and IM as the first scenario. The second scenario was to combine three observed genotypes in rs2853669 with CD8+ cell densities. In both scenarios, certain combinations of CD8+ cell densities with genetic alteration refined prognostic value of our results. Cases presenting β-catenin mutation and high CD8+ cell densities in TC/IM presented much longer TTR (HR = 0.12, p = 0.005/HR = 0.17, p = 0.005 respectively) and longer DFS (HR = 0.25, p = 0.01/HR = 0.34, p = 0.02) than cases presenting only one of these “good” prognostic factors (Table 3 and Fig. 2). Similar to the first scenario, the GG genotype in SNP rs2853669 combined with high density of CD8+ cells presented longer TTR and longer DFS in both tumor regions of interest (Table 3). We could demonstrate that the combination of genetic and immune variable presents a much better prognostic value than individual variables. Multivariable Cox regression was performed for both individual and combined factors. Some individual and combined factors lost their significance after adjustment. We believe that was provoked by size of our group thus, our multivariate results are not reliable.

Our study presents few limitations, the foremost being the low case number. The case numbers limit particularly our analysis for combined factors. Thus, our present data cover only a small segment of genetic and immune landscape of tumor microenvironment. Nevertheless, we believe that a follow up study performed on greater number of cases and exploring different genetic alterations and different types of immune cells may further refine the prognostic tools. These data could be invaluable in treatment decision making.


Our results show that factors both in the tumor and its environment influence HCC survival. The best predictive values for TTR were obtained by combining data for CTNNB1 mutations and CD8+ cell densities. Although we were able to analyze the most common somatic mutations in HCC, we were able to use only one immune cell marker, CD8+ cells, in this first analysis validating the integrative approach. In the future, many more cell types of the tumor microenvironment need to be tested in combination with the genetic markers. Moreover, it seems that combination of TERTp mutation (position − 124 bp) and rs2853669 may be a separate predictive factor for patients with non – viral etiology of HCC. The results, although based on small numbers of resected patients, are encouraging for study further integrative approaches for HCC risk stratification.

Availability of data and materials

The datasets generated and analyzed during the current study are available in the Open Science Forum repository. Data can be accessed openly from the link:


  1. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global Cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71(3):209–49.

    Article  PubMed  Google Scholar 

  2. Llovet JM, Kelley RK, Villanueva A, Singal AG, Pikarsky E, Roayaie S, et al. Hepatocellular carcinoma. Nat Rev Dis Prim. 2021;7:1.

    Google Scholar 

  3. Chen DS, Mellman I. Elements of cancer immunity and the cancer-immune set point. Nature. 2017;541(7637):321–30.

    Article  CAS  PubMed  Google Scholar 

  4. Xu X, Tan Y, Qian Y, Xue W, Wang Y, Du J, et al. Clinicopathologic and prognostic significance of tumor-infiltrating CD8+ T cells in patients with hepatocellular carcinoma: a meta-analysis. Medicine (Baltimore). 2019;98(2):e13923.

    Article  CAS  Google Scholar 

  5. Ding W, Xu X, Qian Y, Xue W, Wang Y, Du J, et al. Prognostic value of tumor-infiltrating lymphocytes in hepatocellular carcinoma: a meta-analysis. Medicine (Baltimore). 2018;97(50):e13301.

    Article  Google Scholar 

  6. Bell RJA, Rube HT, Xavier-Magalhães A, Costa BM, Mancini A, Song JS, et al. Understanding TERT promoter mutations: a common path to immortality. Mol Cancer Res. 2016;14(4):315–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Daniel M, Peek GW, Tollefsbol TO. Regulation of the human catalytic subunit of telomerase (hTERT). Gene. 2012;498(2):135–46.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Dratwa M, Wysoczanska B, Lacina P, Kubik T, Bogunia-Kubik K. TERT-regulation and roles in Cancer formation. Front Immunol. 2020;11:589929.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Chakravarti D, LaBella KA, DePinho RA. Telomeres: history, health, and hallmarks of aging. Cell. 2021;184(2):306–22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Heidenreich B, Kumar R. TERT promoter mutations in telomere biology. Mutat Res Rev Mutat Res. 2017;771:15–31.

    Article  CAS  PubMed  Google Scholar 

  11. Hsu CP, Hsu NY, Lee LW, Ko JL. Ets2 binding site single nucleotide polymorphism at the hTERT gene promoter--effect on telomerase expression and telomere length maintenance in non-small cell lung cancer. Eur J Cancer. 2006;42(10):1466–74.

    Article  CAS  PubMed  Google Scholar 

  12. Nagore E, Rachakonda S, Kumar R. TERT promoter mutations in melanoma survival. Oncotarget. 2019;10(16):1546–8.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Rachakonda PS, Hosen I, de Verdier PJ, Fallah M, Heidenreich B, Ryk C, et al. TERT promoter mutations in bladder cancer affect patient survival and disease recurrence through modification by a common polymorphism. Proc Natl Acad Sci U S A. 2013;110(43):17426–31.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Simon M, Hosen I, Gousias K, Rachakonda S, Heidenreich B, Gessi M, et al. TERT promoter mutations: a novel independent prognostic factor in primary glioblastomas. Neuro-Oncology. 2015;17(1):45–52.

    Article  CAS  PubMed  Google Scholar 

  15. Ko E, Seo HW, Jung ES, Kim BH, Jung G. The TERT promoter SNP rs2853669 decreases E2F1 transcription factor binding and increases mortality and recurrence risks in liver cancer. Oncotarget. 2016;7(1):684–99.

    Article  PubMed  Google Scholar 

  16. Tornesello ML, Buonaguro L, Tatangelo F, Botti G, Izzo F, Buonaguro FM. Mutations in TP53, CTNNB1 and PIK3CA genes in hepatocellular carcinoma associated with hepatitis B and hepatitis C virus infections. Genomics. 2013;102(2):74–83.

    Article  CAS  PubMed  Google Scholar 

  17. Gao C, Wang Y, Broaddus R, Sun L, Xue F, Zhang W. Exon 3 mutations of CTNNB1 drive tumorigenesis: a review. Oncotarget. 2018;9(4):5492–508.

    Article  PubMed  Google Scholar 

  18. Mo Z, Wang Y, Cao Z, Li P, Zhang S. An integrative analysis reveals the underlying association between CTNNB1 mutation and immunotherapy in hepatocellular carcinoma. Front Oncol. 2020;10:853.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Wang Z, Sheng YY, Gao XM, Wang CQ, Wang XY, Lu XU, et al. beta-catenin mutation is correlated with a favorable prognosis in patients with hepatocellular carcinoma. Mol Clin Oncol. 2015;3(4):936–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Galluzzi L, Chan TA, Kroemer G, Wolchok JD, Lopez-Soto A. The hallmarks of successful anticancer immunotherapy. Sci Transl Med. 2018;10(459).

  21. Galon J, Bruni D. Tumor immunology and tumor evolution: intertwined histories. Immunity. 2020;52(1):55–81.

    Article  CAS  PubMed  Google Scholar 

  22. Shalapour S, Karin M. Pas de Deux: control of anti-tumor immunity by Cancer-associated inflammation. Immunity. 2019;51(1):15–26.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Friemel J, Rechsteiner M, Frick L, Böhm F, Struckmann K, Egger M, et al. Intratumor heterogeneity in hepatocellular carcinoma. Clin Cancer Res. 2015;21(8):1951–61.

    Article  CAS  PubMed  Google Scholar 

  24. Yuan X, Cheng G, Yu J, Zheng S, Sun C, Sun Q, et al. The TERT promoter mutation incidence is modified by germline TERT rs2736098 and rs2736100 polymorphisms in hepatocellular carcinoma. Oncotarget. 2017;8(14):23120–9.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Trailin A, Cervenkova L, Ambrozkiewicz F, Ali E, Kasi P, Palek R, Hošek P, Třeška V, Daum O, Tonar Z, Liška V, Hemminki K. T- and B-cells in the inner invasive margin of hepatocellular carcinoma after resection associate with favorable prognosis. Cancers. 2022;14(3):604.

  26. Hendry S, Salgado R, Gevaert T, Russell PA, John T, Thapa B, et al. Assessing tumor-infiltrating lymphocytes in solid tumors: a practical review for pathologists and proposal for a standardized method from the international Immunooncology biomarkers working group: part 1: assessing the host immune response, TILs in invasive breast carcinoma and ductal carcinoma in situ, metastatic tumor deposits and areas for further research. Adv Anat Pathol. 2017;24(5):235–51.

    Article  PubMed  PubMed Central  Google Scholar 

  27. finalfit: Quickly Create Elegant Regression Results Tables and Plots when Modelling.

  28. survminer: Drawing Survival Curves using 'ggplot2'.

  29. A Package for Survival Analysis in R. R package version 3.2–13.

  30. Pezzuto F, Izzo F, Buonaguro L, Annunziata C, Tatangelo F, Botti G, et al. Tumor specific mutations in TERT promoter and CTNNB1 gene in hepatitis B and hepatitis C related hepatocellular carcinoma. Oncotarget. 2016;7(34):54253–62.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Nault JC, Calderaro J, Di Tommaso L, Balabaud C, Zafrani ES, Bioulac-Sage P, et al. Telomerase reverse transcriptase promoter mutation is an early somatic genetic alteration in the transformation of premalignant nodules in hepatocellular carcinoma on cirrhosis. Hepatology. 2014;60(6):1983–92.

    Article  CAS  PubMed  Google Scholar 

  32. Pezzuto F, Buonaguro L, Buonaguro FM, Tornesello ML. Frequency and geographic distribution of TERT promoter mutations in primary hepatocellular carcinoma. Infect Agent Cancer. 2017;12:27.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  33. Pezzuto F, Izzo F, De Luca P, Biffali E, Buonaguro L, Tatangelo F, Buonaguro FM, Tornesello ML. Clinical significance of telomerase reverse-transcriptase promoter mutations in hepatocellular carcinoma. Cancers. 2021;13(15):3771.

  34. Chen YL, Jeng YM, Chang CN, Lee HJ, Hsu HC, Lai PL, et al. TERT promoter mutation in resectable hepatocellular carcinomas: a strong association with hepatitis C infection and absence of hepatitis B infection. Int J Surg. 2014;12(7):659–65.

    Article  CAS  PubMed  Google Scholar 

  35. Lee HW, Park TI, Jang SY, Park SY, Park WJ, Jung SJ, et al. Clinicopathological characteristics of TERT promoter mutation and telomere length in hepatocellular carcinoma. Medicine (Baltimore). 2017;96(5):e5766.

    Article  CAS  Google Scholar 

  36. Rachakonda S, Hoheisel JD, Kumar R. Occurrence, functionality and abundance of the TERT promoter mutations. Int J Cancer. 2021;149(11):1852–62.

    Article  CAS  PubMed  Google Scholar 

  37. Guichard C, Amaddeo G, Imbeaud S, Ladeiro Y, Pelletier L, Maad IB, et al. Integrated analysis of somatic mutations and focal copy-number changes identifies key genes and pathways in hepatocellular carcinoma. Nat Genet. 2012;44(6):694–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Waisberg J, Saba GT. Wnt−/−beta-catenin pathway signaling in human hepatocellular carcinoma. World J Hepatol. 2015;7(26):2631–5.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Wong CM, Fan ST, Ng IO. beta-catenin mutation and overexpression in hepatocellular carcinoma: clinicopathologic and prognostic significance. Cancer. 2001;92(1):136–45.

    Article  CAS  PubMed  Google Scholar 

  40. Mao TL, Chu JS, Jeng YM, Lai PL, Hsu HC. Expression of mutant nuclear beta-catenin correlates with non-invasive hepatocellular carcinoma, absence of portal vein spread, and good prognosis. J Pathol. 2001;193(1):95–101.

    Article  CAS  PubMed  Google Scholar 

  41. Hsu HC, Jeng YM, Mao TL, Chu JS, Lai PL, Peng SY. Beta-catenin mutations are associated with a subset of low-stage hepatocellular carcinoma negative for hepatitis B virus and with favorable prognosis. Am J Pathol. 2000;157(3):763–70.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Khalaf AM, Fuentes D, Morshid AI, Burke MR, Kaseb AO, Hassan M, et al. Role of Wnt/beta-catenin signaling in hepatocellular carcinoma, pathogenesis, and clinical significance. J Hepatocell Carcinoma. 2018;5:61–73.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Hiraoka N. Tumor-infiltrating lymphocytes and hepatocellular carcinoma: molecular biology. Int J Clin Oncol. 2010;15(6):544–51.

    Article  PubMed  Google Scholar 

  44. Hsiao Y-W, Chiu L-T, Chen C-H, Shih W-L, Lu T-P. Tumor-infiltrating leukocyte composition and prognostic power in hepatitis B- and hepatitis C-related hepatocellular carcinomas. Genes. 2019;10(8):630.

  45. Zheng X, Song X, Shao Y, Xu B, Chen L, Zhou Q, et al. Prognostic role of tumor-infiltrating lymphocytes in gastric cancer: a meta-analysis. Oncotarget. 2017;8(34):57386–98.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Huang Y, Ma C, Zhang Q, Ye J, Wang F, Zhang Y, et al. CD4+ and CD8+ T cells have opposing roles in breast cancer progression and outcome. Oncotarget. 2015;6(19):17462–78.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Mei Z, Liu Y, Liu C, Cui A, Liang Z, Wang G, et al. Tumour-infiltrating inflammation and prognosis in colorectal cancer: systematic review and meta-analysis. Br J Cancer. 2014;110(6):1595–605.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Sun C, Xu J, Song J, Liu C, Wang J, Weng C, et al. The predictive value of Centre tumour CD8(+) T cells in patients with hepatocellular carcinoma: comparison with Immunoscore. Oncotarget. 2015;6(34):35602–15.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Bi F, Chen Y, Yang Q. Significance of tumor mutation burden combined with immune infiltrates in the progression and prognosis of ovarian cancer. Cancer Cell Int. 2020;20:373.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Melo CM, Vidotto T, Chaves LP, Lautert-Dutra W, Reis RBD, Squire JA. The role of somatic mutations on the immune response of the tumor microenvironment in prostate Cancer. Int J Mol Sci. 2021;22(17):9550.

Download references


Not applicable.


This project has received funding from the European Union’s Horizon 2020 research and innovation pro-gramme under grant agreement N°856620, a grant of Ministry of Health of the Czech Republic AZV NU21–03-00506, Czech Science Foundation, project no. 21-14082S, project of the Ministry of Education, Youth and Sports of the Czech Republic (INTER-ACTION no. LTA-USA-19032) and the programme Cooperatio no. 43 – Surgical Disciplines.

Author information

Authors and Affiliations



Conceptualization: K.H. and V.L., methodology: F.A., A.T., L.C., R.V., O.D., Z.T., validation: F.A., A.T., R.V. and K.H., formal analysis: F. A, R.V. and A.T., investigation: F.A., A.T., V.H., M.A., resources: V.T., R.P., and V.L., data curation: F.A., A.T., R.P., V.T., L.C., writing original draft: F.A. and K.H., review and editing: A.T., K.H., R.V., V.L., O.D., supervision: K.H., project administration: V.L., funding acquisition: K.H. and V.L. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Filip Ambrozkiewicz.

Ethics declarations

Ethics approval and consent to participate

The Ethical approval for the study was given by the Ethical Board of Faculty of Medicine and University hospital in Pilsen (118/2021, 11 March 2021). The need for informed consent was waived by the Ethical Board of Faculty of Medicine and University hospital in Pilsen (118/2021, 11 March 2021) because of the retrospective nature of the study. This study was carried out in accordance with the ethical standards laid down in the Declaration of Helsinki (2013 version).

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

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

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Ambrozkiewicz, F., Trailin, A., Červenková, L. et al. CTNNB1 mutations, TERT polymorphism and CD8+ cell densities in resected hepatocellular carcinoma are associated with longer time to recurrence. BMC Cancer 22, 884 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: