Comparative Transcriptional Signature Analysis of Peripheral Blood Mononuclear Cells in Early Stage of Hepatitis B-related Hepatocellular Carcinoma

authors:

avatar Omid Gholizadeh 1 , 2 , avatar Hossein Bannazadeh Baghi ORCID 1 , 2 , avatar Mahin Ahangar Oskouee 2 , avatar Nader Mohammadzadeh 3 , avatar Behrooz Naghili 2 , avatar Narges Eslami 2 , avatar Mohsen Nasiri-toosi 4 , avatar Alka Hasani 2 , avatar Mohammad Ahngarzadeh Rezaei 2 , avatar Tina Taikandi 5 , avatar Vahdat Poortahmasebi ORCID 1 , 2 , *

Infectious and Tropical Diseases Research Center, Tabriz University of Medical Sciences, Tabriz, Iran
Department of Bacteriology and Virology, School of Medicine, Tabriz University of Medical Sciences, Tabriz, Iran
Central Laboratory of East Azerbaijan Province, Tabriz University of Medical Sciences, Tabriz, Iran
Liver Transplantation Research center, Tehran University of Medical Sciences, Tehran, Iran
Department of Biology, Faculty of Life Sciences, Denesh Alborz University, Abyek, Iran

how to cite: Gholizadeh O, Bannazadeh Baghi H, Ahangar Oskouee M, Mohammadzadeh N, Naghili B, et al. Comparative Transcriptional Signature Analysis of Peripheral Blood Mononuclear Cells in Early Stage of Hepatitis B-related Hepatocellular Carcinoma. Hepat Mon. 2022;22(1):e130862. https://doi.org/10.5812/hepatmon-130862.

Abstract

Background:

Hepatocellular carcinoma (HCC) is a prevalent and life-threatening tumor with high morbidity and mortality. Proper prediction and prognosis are incredibly stressed to diagnose HCC and increase patient survival.

Objectives:

This research aims to evaluate gene expression levels of pre-differentiated transcripts for those suffering from chronic hepatitis B (CHB) and HCC.

Methods:

To examine the previously analyzed peripheral blood mononuclear cells (PBMCs) transcriptomic array data, we selected seven differentially expressed genes (DEGs) in normal versus CHB and CHB versus HCC (CD44, SP3, USP8, E2F2, UFM1, IFN regulative factor binding protein 2 (IRF2BP2), and T-cell intracellular antigen 1 (TIA1)). The study included individuals with treatment-naïve CHB (n = 30) and primary HCC (n = 25) and healthy controls (n = 15). Subsequently, the expression of genes was assayed using qRT-PCR. A phylogenetic evaluation was performed using direct sequencing of HBsAg.

Results:

In HCC patients, 60% (n = 15) were HBeAg-positive. HBeAg was negative in all CHB patients, but all were anti-HBe-positive. The hepatitis B virus (HBV) load of HCC patients was more than that of CHB subjects. All patients were of the Iranian race and HBV D genotype. The expression of five transcriptional markers (CD44, SP3, USP8, E2F2, and UFM1) was higher in HCC patients than in CHB and healthy subjects, which was similar to the initial microarray data analysis.

Conclusions:

Transcriptional signatures may be related to the pathogenesis of HCC and used as diagnostic biological markers for the initial monitoring and prediction of HCC.

1. Background

Hepatitis B virus (HBV) possesses hepatic tropism and annular double-stranded DNA. This virus causes acute and chronic hepatitis B (CHB) in humans. Chronic hepatitis B (CHB) is the ultimate central stimulus for worldwide cirrhosis, hepatocellular carcinoma (HCC), and hepatic failure (1). Although effective antiviral drugs and vaccines are available, HBV disease remains a well-being trouble. Furthermore, the virus still infects approximately 350 million people and contributes to the morbidity and mortality associated with viral hepatitis (2). Therefore, discovering new methods is necessary to develop effective early diagnostic and treatment strategies. A minimally invasive procedure, such as obtaining blood, can be an attractive alternative to liver biopsy (3).

Peripheral blood mononuclear cells (PBMCs) are an essential immune system element. Interestingly, examining critical genes in PBMCs has shown the relationship between altered gene expression and diseases. In different diseases, PBMCs are used for monitoring gene expression. It has been proposed that changed gene expression in PBMCs might be employed as a potential clinical tool for collecting data since it can give insight into the diagnosis or prognosis of the disease (4). The identification of new biomarkers may assist in disease management and interventions to decrease mortality rates. For example, examining the transcriptional profile modifications in HCC and CHB patients’ PBMCs might help predict and determine HCC (5).

Our previous study acquired comprehensive transcription analyses from the Gene Expression Omnibus (GEO) database with code GSE58208 (6). This databank includes PBMCs to evaluate the transcript in three subjects: CHB, HCC, and healthy controls. We showed that most differentially expressed genes (DEGs) in CHB carriers and HCC patients are related to the immune response. Also, 208 genes were differentially expressed among normal and CHB patients (normal-CHB). Moreover, 186 genes were differentially expressed among CHB and HCC patients (normal-CHB). Our analysis found that the normal-CHB and CHB-HCC networks only comprised 201 and 179 DEGs. Seven common genes were found in both groups, including CD44, SP3, USP8, E2F2, UFM1, IFN regulative factor binding protein 2 (IRF2BP2), and T-cell intracellular antigen 1 (TIA1) (Figure 1). Of interest, six of these seven DEGs were upregulated in CHB-HCC compared with normal-CHB. However, an independent experimental method should evaluate gene expression data to ensure the obtained results. Accordingly, our study examined the expression of genes mentioned above using qRT-PCR in PBMC specimens from CHB and HCC.

“Normal- chronic hepatitis B (CHB)” and “CHB-hepatocellular carcinoma (HCC) states” have seven overlapping differentially expressed genes (DEGs). CD44, E2F2, USP8, IRF2BP2, and SP3 are six DEGs that were upregulated in CHB-HCC and downregulated in normal-CHB.
“Normal- chronic hepatitis B (CHB)” and “CHB-hepatocellular carcinoma (HCC) states” have seven overlapping differentially expressed genes (DEGs). CD44, E2F2, USP8, IRF2BP2, and SP3 are six DEGs that were upregulated in CHB-HCC and downregulated in normal-CHB.

2. Objectives

This investigation analyzed the gene expression pattern of CD44, SP3, USP8, E2F2, UFM1, IRF2BP2, and TIA1 obtained from transcriptome analysis.

3. Methods

3.1. Patients

This study was conducted between August 2020 and February 2022. Eighty-five whole blood samples were collected from healthy subjects (n = 30), treatment-naïve CHB patients (n = 30), and primary HCC patients (n = 25). Demographic information, such as gender and age, and laboratory information, like serological and biochemical tests, were obtained from all patients. Patients with immunosuppression, hepatitis C virus (HCV), and hepatitis D virus were excluded. Healthy controls without HBV, HCV, or human immunodeficiency virus, were selected; they also had normal liver function tests and no history of liver disease. In addition, HBsAg positivity, levels of serum alanine transaminase (ALT) and aspartate transaminase (AST) for at least six months, HBsAg/Anti-HBe status, HBV viral load, and histopathological evidence were considered criteria for inactive CHB carriers. All participants gave informed consent, and the Tabriz University of Medical Sciences Ethics Committee authorized the study.

3.2. Laboratory Tests, Hepatitis B Virus DNA Extraction, and Quantitative Real-time PCR

The ALT and AST levels were determined in all 85 blood samples. Furthermore, HBV serological markers such as HBsAg, HBeAg, and anti-HBe were assessed using commercialized enzyme-linked immunosorbent assay (ELISA) kits (Acon, USA) as maintained in the factory’s protocol. Hepatitis B virus DNA was sequestered from a 200 μL aliquot of serum utilizing a High Pure Viral Nucleic Acid Kit (Roche, Germany) as reported in the manufacturer’s protocol. Next, to elute the DNA, 50 μL of elution buffer was utilized, which was subsequently stored at -70°C. Regardless of serological results, HBV viral load was determined in all specimens using a qRT-PCR kit (Fast-Track Diagnostics, Luxembourg) as reported in the manufacturer’s protocol. The kit protocol specifies a lower limit of detection (LOD) of 500 copies/mL.

3.3. Nested PCR, Sequencing, and Phylogenetic Analysis of HBsAg

The complete sequences of the HBsAg gene were amplified for all positive samples using specific primers. As depicted, two rounds of nested PCR (N-PCR) were used to expand the HBsAg area within the virus genome (7). The whole HBsAg segment of the HBV genome was sequenced by ABI 3130 Hereditary Analyzer (Hereditary Analyzer ABI-3130 DNA Sequencer, Cultivate City, CA, USA). Phylogenetic analysis was used for genotyping HBV genotypes A through H utilizing reference sequences. A maximum likelihood (ML) phylogenetic tree was constructed using the alignment of HBsAg region of patients by rooting with a woodchuck hepatitis virus (WHV) sequence. The MEGA X program was used to generate the phylogenetic tree, and the Kimura two-parameter matrix was used to measure genetic distance (8). Bootstrap resampling and rebuilding were performed 1000 times to validate the phylogenetic tree’s reliability (9).

3.4. Peripheral Blood Mononuclear Cell Isolation, RNA Extraction, and Complementary DNA Synthesis

Ficoll density gradient centrifugation (Ambion, USA) extracted PBMCs from 5 mL of fresh whole EDTA-treated blood for 20 minutes at 20°C at 1,000x g. PBMCs were pelleted in phosphate-buffered saline (PBS). Washed and pellet cells were put in 1 mL of RPMI before being taken to the RNA extraction step. The whole RNA was isolated from PBMCs utilizing TRIzol reagent (Ambion, USA), following the manufacturer’s instructions. RNase-free DNase I (Qiagen, Germany) was used to annihilate genomic DNA impurity from RNA samples. The concentration of total RNA extracted was measured using a nanodrop spectrophotometer (Boeco, Germany). Before being stored at -70°C, the isolated RNA was adjusted to a final concentration of 1 mg/ml. Finally, complementary DNA (cDNA) synthesis was carried out utilizing the cDNA synthesis kit (Qiagen, Hilden, Germany) based on the manufacturer’s protocols.

3.5. Relative Quantitative Real-time PCR

The whole RNA was extracted from PBMCs utilizing TRIzol reagent (Ambion, USA) in conformity with the producer’s recommendations. The relative expression of mRNAs (CD44, SP3, USP8, E2F2, UFM1, IRF2BP2, and TIA1) turned into measured using Real Q Plus 2x Master Mix Green (Ampliqon, Denmark). PCR protocols were as follows: initial denaturation at 95°C for 10 min, followed by 40 cycles of PCR at 95°C for 15 s and 57°C for 30 s. Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was used as a reference gene to evaluate the linearity and accuracy of RT-PCR. Table 1 lists the primers used in this assay. Based on the 2-ΔΔCT method, mRNA expression templates were investigated (10). Relative values of each gene were expressed as fold change to compare mRNA levels between patients.

Table 1.

Sequences of Primers Used in qRT-PCR

Genes and DirectionPrimer Sequence (5’ to 3’)
CD44
ForwardCTGCAGGTATGGGTTCATAG
ReverseATATGTGTCATACTGGGAGGTG
UFM1
ForwardCCGTTCACAGCAGTGCTAAA
ReverseCCGTGCTTCAGGAAAACA
E2F2
ForwardCCTTGGAGGCTACTGACAGC
ReverseCCACAGGTAGTCGTCCTGGT
TIA1
ForwardGTGTCCTGAGAGAGCAAGCA
ReverseCCACTCGGGCATCTCTCAAA
SP3
ForwardCCAGGATGTGGTAAAGTCTA
ReverseCTCCATTGTCTCATTTCCAG
IRF2BP2
ForwardAGTTCTGCTTCCCTTGCTCC
ReverseAGGATGGTTGCAATTTCCCC
USP8
ForwardGACGCCACCTGCATCTATAGAAG
ReverseGGAAAGTAAAACTGTCCTGCGCAA
GAPDH
ForwardCCCTGAACGCTATCACTGCA
ReverseCCATCACGCCACAGTTTC

3.6. Statistical Analysis

The data were expressed as percentage and mean ± standard deviation. The differences of categorical variables were investigated using chi-square test. The normal distribution was estimated using the Kolmogorov-Smirnov (K-S) test. The Mann-Whitney U test was used for non-normal continuous variables. IBM SPSS Statistics 24 (SPSS Inc., Chicago, Illinois, USA) was used for statistical analysis. P-values lower than 0.05 were supposed to be significant for all tests.

4. Results

4.1. General Characteristics of Patients

Demographic and clinical specifications of all 85 subjects, including healthy controls (n = 30), CHB patients (n = 30), and HCC patients (n = 25), are summarized in Table 2. Of the 85 participants, 44 (51.8%) and 41 (48.2%) were male and female. There was no significant difference between the studied groups in gender (P-value = 0.213). The average age of CHB and HCC patients was 45.86 ± 11.02 and 58.60 ± 6.92, respectively. Statistically, there was a significant difference between CHB and HCC patients in age (P-value < 0.001). In terms of liver injury enzymes (ALT and AST), we found statistical differences between HCC and CHB patients (P value < 0.001). Also, the HBV viral load was higher in HCC patients compared to the CHB patients (P-value < 0.0001).

Table 2.

Clinical Characteristics of Participants a

FactorsAll Subjects (n = 85)Healthy Controls (n = 30)Chronic Hepatitis B Virus (n = 30)Hepatocellular Carcinoma (n = 25)P-Value
Sex0.213
Male44 (51.8)16 (53.3)12 (40)16 (51.8)
Female41 (48.2)14 (46.7)18 (60)9 (48.2)
Age (y)50.22 ± 10.6047.07 ± 8.1945.86 ± 11.0258.60 ± 6.92< 0.001
ALT (IU/L)42.70 ± 22.3826.80 ± 4.0234.57 ± 13.5066.20 ± 20.61< 0.001
AST (IU/L)42.39 ± 23.3724.93 ± 5.4736.57 ± 9.4665.27 ± 25.37< 0.001
Viral load (Log2 copies/mL)15.12 ± 3.37-12.42 ± 1.1017.83 ± 2.57< 0.001
HBeAg< 0.001
Positive20 (23.5)0 (0)0 (0)15 (60)
Negative65 (76.5)30 (100)30 (100)10 (10)

In HCC patients, 60% (n = 15) were positive for HBeAg. All CHB patients were negative for HBeAg, but all were positive for anti-HBe. A comparison was made between HCC and CHB subjects. Statistical analysis revealed a significant difference between HCC and CHB patients in HBeAg status (P-value < 0.001). HBV DNA was quantified in serum samples. The HBV load was significantly higher in HCC patients than in CHB subjects (P-value < 0.001).

4.2. Phylogenetic and Mutational Analysis

The sufferers were of the Iranian race and genotype D for HBV (Figure 2). Several mutations were observed in the “a” determinant location of HBsAg (amino acid residues 124 to 147). Six (3.6%) affected individuals had an amino acid spurt in the “a” designation. According to the number of mutations/sites, the average nucleotide mutation frequency (dN/dS) across all sequences was less than one. A dN/dS ratio of less than one indicated negative selection pressure or sequence stability for the HBsAg gene.

Most extreme probability phylogenetic tree development based on the ClustalW alignment of the HBsAg arrangements (681-bp) obtained from Iranian patients
Most extreme probability phylogenetic tree development based on the ClustalW alignment of the HBsAg arrangements (681-bp) obtained from Iranian patients

4.3. Distinct mRNA Signature in Hepatocellular Carcinoma Patients

We evaluated the relative expression of a panel of genes (CD44, SP3, USP8, E2F2, UFM1, IRF2BP2, and TIA1) and detected a specific mRNA signature in HCC sufferers compared to CHB sufferers and healthy subjects (Figure 3). PCR results showed that five out of seven genes (CD44, SP3, USP8, E2F2, and UFM1) were significantly upregulated in the PBMCs of HCC patients, which was almost similar to the results of the initial microarray data analysis. Patients with HCC were comparable to normal subjects, and the levels were identical to the previous results (data not shown). These findings indicate that the abovementioned genes might be associated with HCC progression. Our last microarray surveys showed that the mRNA level of IRF2BP2 was considerably upregulated in HCC patients. Also, the mRNA level of TIA1 was downregulated in CHB sufferers. However, qRT-PCR data showed that the dysregulation of IRF2BP2 and TIA1 genes in the patient group was mild and statistically insignificant.

The relative expression alterations of the seven differentially expressed genes identified by PCR were compared between hepatocellular carcinoma patients, chronic hepatitis Bpatients, and healthy Controls. * indicates statistical significance.
The relative expression alterations of the seven differentially expressed genes identified by PCR were compared between hepatocellular carcinoma patients, chronic hepatitis Bpatients, and healthy Controls. * indicates statistical significance.

5. Discussion

Hepatocellular carcinoma is a leading cause of death in individuals with cirrhosis and is one of the significant cancers worldwide (11). In parts of Asia, the higher rate of HCC resulted from a high prevalence of hepatitis B and C virus infections (12). Hepatocellular carcinoma incidence has been on the rise for the past 40 years, and the diagnosis of HCC can often be challenging (13). Despite significant advances in the diagnosis and treatment of HCC over a long period, the prognosis of HCC remains poor (14). To improve the prognosis of HCC, diagnosing the tumor in the early stages is crucial when patients are asymptomatic and their liver function is preserved. This situation enables effective therapies with significant survival benefits to be applied (15). As a result, investigating the biomarkers for the early detection, treatment, and prognosis of HCC is important and fundamental.

For a long time, inventive bioinformatics and computational procedures have been successfully connected to a wide assortment of oncology inquiries (16). Microarray and high-throughput innovations have given an effective method to filter thousands of genes that might be implicated in the occurrence and advancement of HCC. Hence, these genes can serve as potential targets for diagnosing and treating HCC (17). Furthermore, the microarray technique provides an effective tool to explore gene regulation patterns and molecular mechanisms associated with the oncogenesis and progression of HCC. Numerous studies from other researchers and a couple of our previous investigations have merely been perpetrated based on a pure system biology approach and a bioinformatics attitude for gene expression analysis. Our previous study was a computational analysis and network-based investigation (6) but not an experimental examination. However, further validation tests may lead to an improvement in our credibility. Also, others have used PBMC transcriptome analysis to elucidate pathogenesis and detect specific biomarkers, for instance, in individuals suffering from hepatitis B-related acute-on-chronic hepatic failure (18).

In our prior investigation, we conducted a network-based study for the transcriptome of PBMCs in three subjects, including healthy controls and CHB and HCC patients. CD44, SP3, USP8, E2F2, UFM1, IRF2BP2, and TIA1 were normal in patients with CHB-CHB and CHB-HCC systems. Except for TIA1, the expression of other genes was upregulated in CHB-HCC and downregulated in normal-CHB subject (6). In this investigation, we evaluated the gene expression data of these seven DEGs by an experimental method to validate the results already obtained. To our knowledge, no study has detailed the expression levels or prognostic value of the previously mentioned genes in the PBMCs of HCC patients.

Our investigation indicated that five mRNAs, including CD44, SP3, USP8, E2F2, and UFM1, could remarkably discriminate between patients with HCC and CHB, which is consistent with our previous computational study. CD44 involves several physiological processes, including organ development, immune function, and hematopoiesis. This gene controls various pathways that alter cancer development, invasion, metastasis, and resistance to treatment (19). Several studies have reported tissue expression of CD44 as a favorable predictor in HCC patients (20). Furthermore, in our study, CD44 was upregulated entirely in the PBMCs of HCC patients. E2Fs are a family of essential translation factors that interact in numerous tissues and organs. Cell proliferation, sequestration, DNA repair, cell cycle direction, and apoptosis are all molecular capacities of E2Fs (21). Reports have detailed the expression or function of E2Fs in several human malignancies. According to studies, E2Fs may be a biomarker for tumor prediction (22, 23). In our research, E2F2 was generally upregulated in the PBMCs of HCC patients. Huang et al., using a network-based analysis, demonstrated that tissue expression of E2Fs had a predictive value in various stages and pathological grades of HCC (24). Evaluating the E2F expression in the PBMCs of HCC patients and its predictive value has not been investigated. However, other studies have demonstrated that the blood E2F expression in patients with lung and prostate cancer can be examined by qRT-PCR (25, 26). SP3, in the particular protein/Krüppel-like factor transcription factor family (27), is significantly upregulated in the PBMCs of HCC patients. Microarray data analysis from the GSE58208 dataset revealed that CD44, SP3, and E2F2 genes were upregulated in patients with HCC (6). In our investigation, we observed that these genes were considerably overexpressed in HCC patients compared to CHB patients. In any case, no striking contrast between the two mRNAs was found in CHB patients. Ubiquitin-fold modifier 1 (UFM1) is a small molecule ubiquitin protein with a similar function to ubiquitin. In several pathophysiological forms, UFM1 and its modification framework are embroiled (28).

USP8 has been recognized as an independent prognostic marker for cervical squamous cell carcinoma (CSCC) in various studies. Therefore, it may act as a diagnostic and therapeutic target in CSCC patients (29, 30). Notably, UFM1 and USP8 were generally upregulated in the PBMCs of HCC patients. T-cell intracellular antigen 1 is involved in several biological processes, including RNA metabolism, both within the cytoplasm and in the nucleus (31). IFN regulative factor binding protein 2 is an IRF-2-dependent transcriptional repressor involved in cell apoptosis, survival, and differentiation (32). However, qRT-PCR results revealed slight and statistically insignificant dysregulation of the TIA1 and IRF2BP2 genes in HCC and CHB patients. Our findings are the first to demonstrate that HCC PBMCs have considerably different transcriptional signatures than CHB PBMCs and that higher gene expression is associated with HCC progression. In this research, our gene expression results were obtained from Singaporean patients (33), but we investigated nominee genes in Iranian individuals with CHB and HCC. Ethnicity emerges not because of contrasts in genetic makeup but because of differences in the expression of shared genes between ethnic groups (34). As a result, various hosts and viral variables, including as high viral load, gender, and HBsAg levels, are significant in the evolution of HCC. No serological tests for HBeAg, anti-HBe, and HBV genotypes are performed in the GSE58208 dataset (33). In our research, the direct sequencing of 55 HBsAg sequences showed that all of them were for HBV genotype D. Subsequent studies in different Iranian districts revealed that genotype D is the same genotype of HBV that has spread throughout the country (35). The HBV viral load and genotyping are now essential indicators of the development of antiviral drug resistance and the severity of HBV disease (36, 37).

5.1. Conclusions

Our study showed that different HBV disease states are associated with differentially expressed of several mRNAs signatures. The present transcriptional signature may be employed for the early detection of HBV-related HCC and might play an essential role in the pathogenesis of HBV-induced HCC.

Acknowledgements

References

  • 1.

    Yuen MF, Chen DS, Dusheiko GM, Janssen HLA, Lau DTY, Locarnini SA, et al. Hepatitis B virus infection. Nat Rev Dis Primers. 2018;4:18035. [PubMed ID: 29877316]. https://doi.org/10.1038/nrdp.2018.35.

  • 2.

    Polaris Observatory Collaborators. Global prevalence, treatment, and prevention of hepatitis B virus infection in 2016: a modelling study. Lancet Gastroenterol Hepatol. 2018;3(6):383-403. [PubMed ID: 29599078]. https://doi.org/10.1016/S2468-1253(18)30056-6.

  • 3.

    Chin JL, Pavlides M, Moolla A, Ryan JD. Non-invasive Markers of Liver Fibrosis: Adjuncts or Alternatives to Liver Biopsy? Front Pharmacol. 2016;7:159. [PubMed ID: 27378924]. [PubMed Central ID: PMC4913110]. https://doi.org/10.3389/fphar.2016.00159.

  • 4.

    Li X, Wang Y, Han C, Li P, Zhang H. Colorectal cancer progression is associated with accumulation of Th17 lymphocytes in tumor tissues and increased serum levels of interleukin-6. Tohoku J Exp Med. 2014;233(3):175-82. [PubMed ID: 24990228]. https://doi.org/10.1620/tjem.233.175.

  • 5.

    Lei B, Zhou J, Xuan X, Tian Z, Zhang M, Gao W, et al. Circular RNA expression profiles of peripheral blood mononuclear cells in hepatocellular carcinoma patients by sequence analysis. Cancer Med. 2019;8(4):1423-33. [PubMed ID: 30714679]. [PubMed Central ID: PMC6488130]. https://doi.org/10.1002/cam4.2010.

  • 6.

    Poortahmasebi V, Alavian SM, Nasiri-Toosi M, Norouzi M, Hosseini M, Jazayeri SM. Transcriptome analysis of peripheral blood mononuclear cells from chronic hepatitis B and hepatocellular carcinoma patients: a network-based attitude. Future Virol. 2017;12(12):729-37. https://doi.org/10.2217/fvl-2017-0087.

  • 7.

    Jazayeri SM, Basuni AA, Sran N, Gish R, Cooksley G, Locarnini S, et al. HBV core sequence: definition of genotype-specific variability and correlation with geographical origin. J Viral Hepat. 2004;11(6):488-501. [PubMed ID: 15500549]. https://doi.org/10.1111/j.1365-2893.2004.00534.x.

  • 8.

    Kimura M. A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol. 1980;16(2):111-20. [PubMed ID: 7463489]. https://doi.org/10.1007/BF01731581.

  • 9.

    Efron B, Halloran E, Holmes S. Bootstrap confidence levels for phylogenetic trees. Proc Natl Acad Sci U S A. 1996;93(23):13429-34. [PubMed ID: 8917608]. [PubMed Central ID: PMC24110]. https://doi.org/10.1073/pnas.93.23.13429.

  • 10.

    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. [PubMed ID: 11846609]. https://doi.org/10.1006/meth.2001.1262.

  • 11.

    Ferlay J, Soerjomataram I, Dikshit R, Eser S, Mathers C, Rebelo M, et al. Cancer incidence and mortality worldwide: sources, methods and major patterns in GLOBOCAN 2012. Int J Cancer. 2015;136(5):E359-86. [PubMed ID: 25220842]. https://doi.org/10.1002/ijc.29210.

  • 12.

    Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394-424. [PubMed ID: 30207593]. https://doi.org/10.3322/caac.21492.

  • 13.

    European Association For The Study Of The Liver; European Organisation For Research; Treatment Of Cancer. EASL-EORTC clinical practice guidelines: management of hepatocellular carcinoma. J Hepatol. 2012;56(4):908-43. [PubMed ID: 22424438]. https://doi.org/10.1016/j.jhep.2011.12.001.

  • 14.

    Okajima W, Komatsu S, Ichikawa D, Miyamae M, Ohashi T, Imamura T, et al. Liquid biopsy in patients with hepatocellular carcinoma: Circulating tumor cells and cell-free nucleic acids. World J Gastroenterol. 2017;23(31):5650-68. [PubMed ID: 28883691]. [PubMed Central ID: PMC5569280]. https://doi.org/10.3748/wjg.v23.i31.5650.

  • 15.

    Ayuso C, Rimola J, Vilana R, Burrel M, Darnell A, Garcia-Criado A, et al. Diagnosis and staging of hepatocellular carcinoma (HCC): current guidelines. Eur J Radiol. 2018;101:72-81. [PubMed ID: 29571804]. https://doi.org/10.1016/j.ejrad.2018.01.025.

  • 16.

    Xu X, Wang X, Fu B, Meng L, Lang B. Differentially expressed genes and microRNAs in bladder carcinoma cell line 5637 and T24 detected by RNA sequencing. Int J Clin Exp Pathol. 2015;8(10):12678-87. [PubMed ID: 26722457]. [PubMed Central ID: PMC4680402].

  • 17.

    Xie S, Jiang X, Zhang J, Xie S, Hua Y, Wang R, et al. Identification of significant gene and pathways involved in HBV-related hepatocellular carcinoma by bioinformatics analysis. PeerJ. 2019;7. e7408. [PubMed ID: 31392101]. [PubMed Central ID: PMC6677124]. https://doi.org/10.7717/peerj.7408.

  • 18.

    Zhou Q, Ding W, Jiang L, Xin J, Wu T, Shi D, et al. Comparative transcriptome analysis of peripheral blood mononuclear cells in hepatitis B-related acute-on-chronic liver failure. Sci Rep. 2016;6:20759. [PubMed ID: 26861114]. [PubMed Central ID: PMC4748289]. https://doi.org/10.1038/srep20759.

  • 19.

    Xu H, Niu M, Yuan X, Wu K, Liu A. CD44 as a tumor biomarker and therapeutic target. Exp Hematol Oncol. 2020;9(1):36. [PubMed ID: 33303029]. [PubMed Central ID: PMC7727191]. https://doi.org/10.1186/s40164-020-00192-0.

  • 20.

    Luo Y, Tan Y. Prognostic value of CD44 expression in patients with hepatocellular carcinoma: meta-analysis. Cancer Cell Int. 2016;16:47. [PubMed ID: 27330410]. [PubMed Central ID: PMC4912706]. https://doi.org/10.1186/s12935-016-0325-2.

  • 21.

    Attwooll C, Lazzerini Denchi E, Helin K. The E2F family: specific functions and overlapping interests. EMBO J. 2004;23(24):4709-16. [PubMed ID: 15538380]. [PubMed Central ID: PMC535093]. https://doi.org/10.1038/sj.emboj.7600481.

  • 22.

    Chen HZ, Tsai SY, Leone G. Emerging roles of E2Fs in cancer: an exit from cell cycle control. Nat Rev Cancer. 2009;9(11):785-97. [PubMed ID: 19851314]. [PubMed Central ID: PMC3616489]. https://doi.org/10.1038/nrc2696.

  • 23.

    Emmrich S, Putzer BM. Checks and balances: E2F-microRNA crosstalk in cancer control. Cell Cycle. 2010;9(13):2555-67. [PubMed ID: 20581444]. https://doi.org/10.4161/cc.9.13.12061.

  • 24.

    Huang YL, Ning G, Chen LB, Lian YF, Gu YR, Wang JL, et al. Promising diagnostic and prognostic value of E2Fs in human hepatocellular carcinoma. Cancer Manag Res. 2019;11:1725-40. [PubMed ID: 30863181]. [PubMed Central ID: PMC6388971]. https://doi.org/10.2147/CMAR.S182001.

  • 25.

    Pipinikas CP, Nair SB, Kirby RS, Carter ND, Fenske CD. Measurement of blood E2F3 mRNA in prostate cancer by quantitative RT-PCR: a preliminary study. Biomarkers. 2007;12(5):541-57. [PubMed ID: 17701752]. https://doi.org/10.1080/13547500701391353.

  • 26.

    Al Ahmed HA, Nada O. E2F3 transcription factor: A promising biomarker in lung cancer. Cancer Biomark. 2017;19(1):21-6. [PubMed ID: 28269748]. https://doi.org/10.3233/CBM-160196.

  • 27.

    Suske G, Bruford E, Philipsen S. Mammalian SP/KLF transcription factors: bring in the family. Genomics. 2005;85(5):551-6. [PubMed ID: 15820306]. https://doi.org/10.1016/j.ygeno.2005.01.005.

  • 28.

    Ha BH, Jeon YJ, Shin SC, Tatsumi K, Komatsu M, Tanaka K, et al. Structure of ubiquitin-fold modifier 1-specific protease UfSP2. J Biol Chem. 2011;286(12):10248-57. [PubMed ID: 21228277]. [PubMed Central ID: PMC3060479]. https://doi.org/10.1074/jbc.M110.172171.

  • 29.

    Yan M, Zhao C, Wei N, Wu X, Cui J, Xing Y. High Expression of Ubiquitin-Specific Protease 8 (USP8) Is Associated with Poor Prognosis in Patients with Cervical Squamous Cell Carcinoma. Med Sci Monit. 2018;24:4934-43. [PubMed ID: 30010158]. [PubMed Central ID: PMC6067021]. https://doi.org/10.12659/MSM.909235.

  • 30.

    Poondla N, Chandrasekaran AP, Kim KS, Ramakrishna S. Deubiquitinating enzymes as cancer biomarkers: new therapeutic opportunities? BMB Rep. 2019;52(3):181-9. [PubMed ID: 30760385]. [PubMed Central ID: PMC6476481]. https://doi.org/10.5483/BMBRep.2019.52.3.048.

  • 31.

    Sanchez-Jimenez C, Izquierdo JM. T-cell intracellular antigens in health and disease. Cell Cycle. 2015;14(13):2033-43. [PubMed ID: 26036275]. [PubMed Central ID: PMC4614917]. https://doi.org/10.1080/15384101.2015.1053668.

  • 32.

    Ramalho-Oliveira R, Oliveira-Vieira B, Viola JPB. IRF2BP2: A new player in the regulation of cell homeostasis. J Leukoc Biol. 2019;106(3):717-23. [PubMed ID: 31022319]. https://doi.org/10.1002/JLB.MR1218-507R.

  • 33.

    Shi M, Chen MS, Sekar K, Tan CK, Ooi LL, Hui KM. A blood-based three-gene signature for the non-invasive detection of early human hepatocellular carcinoma. Eur J Cancer. 2014;50(5):928-36. [PubMed ID: 24332572]. https://doi.org/10.1016/j.ejca.2013.11.026.

  • 34.

    Spielman RS, Bastone LA, Burdick JT, Morley M, Ewens WJ, Cheung VG. Common genetic variants account for differences in gene expression among ethnic groups. Nat Genet. 2007;39(2):226-31. [PubMed ID: 17206142]. [PubMed Central ID: PMC3005333]. https://doi.org/10.1038/ng1955.

  • 35.

    Norouzi M, Ramezani F, Khedive A, Karimzadeh H, Alavian SM, Malekzadeh R, et al. Hepatitis B Virus Genotype D is the Only Genotype Circulating in Iranian Chronic Carriers, the Unique Pattern of Genotypic Homogeneity. J Gastroenterol Hepatol Res. 2014;3(9):1238-43. https://doi.org/10.17554/j.issn.2224-3992.2014.03.418.

  • 36.

    Garson JA, Grant PR, Ayliffe U, Ferns RB, Tedder RS. Real-time PCR quantitation of hepatitis B virus DNA using automated sample preparation and murine cytomegalovirus internal control. J Virol Methods. 2005;126(1-2):207-13. [PubMed ID: 15847939]. https://doi.org/10.1016/j.jviromet.2005.03.001.

  • 37.

    Soriano V, Puoti M, Peters M, Benhamou Y, Sulkowski M, Zoulim F, et al. Care of HIV patients with chronic hepatitis B: updated recommendations from the HIV-Hepatitis B Virus International Panel. AIDS. 2008;22(12):1399-410. [PubMed ID: 18614862]. https://doi.org/10.1097/QAD.0b013e3282f8b46f.