Homeobox Genes and Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is the sixth most common type of cancer, and is the third leading cause of cancer-related deaths each year. It involves a multi-step progression and is strongly associated with chronic inflammation induced by the intake of environmental toxins and/or viral infections (i.e., hepatitis B and C viruses). Although several genetic dysregulations are considered to be involved in disease progression, the detailed regulatory mechanisms are not well defined. Homeobox genes that encode transcription factors with homeodomains control cell growth, differentiation, and morphogenesis in embryonic development. Recently, more aberrant expressions of Homeobox genes were found in a wide variety of human cancer, including HCC. In this review, we summarize the currently available evidence related to the role of Homeobox genes in the development of HCC. The objective is to determine the roles of this conserved transcription factor family and its potential use as a therapeutic target in future investigations.


Epidemiology
Hepatocellular carcinoma (HCC) is the most common type of primary liver cancer and is ranked as the second leading cause of cancer-related mortality worldwide [1]. The global annual death as a result of HCC was reported to be 700,000 patients [2]. Intriguingly, its incidence worldwide differs from heterogeneous prevalence of risk factors. The highest incidence of HCC has been reported in East/Southeast Asia and Africa and the lowest in South/Central Asia and Europe [3]. Commonly mentioned risk factors for HCC include chronic viral hepatitis (i.e., hepatitis B virus (HBV) and hepatitis C virus (HCV)), chronic alcohol use, environment pollutants, obesity, and diabetes mellitus. Chronic hepatitis B and exposure to aflatoxin are major risk factors for the occurrence of HCC in sub-Saharan Africa and East Asia, whereas chronic hepatitis C is the major risk factor in the USA, Europe, and Japan [4]. In recent years, nonalcoholic fatty liver disease (NAFLD) has gradually emerged as a leading cause of HCC in Western and Asian populations [5].

Etiologies
Various possible mechanisms, which link these risk factors to hepatocarcinogenesis, have been proposed. HBV-an enveloped DNA virus-belongs to the Hepadnaviridae family and includes eight genotypes (i.e., A to H), which have their respective geographical distribution [6]. Studies have shown that the HBV X protein (HBx) is a 154-amino acid polypeptide that plays an essential role in the development of HCC. HBx may directly promote hepatocytes transforming into tumor-initiating cells through the activation of Wnt/β-catenin signaling [7]. HCV is a single-stranded positive RNA virus belonging to the Hepacivirus genus in the family Flavivaridae, including seven major genotypes [8]. HCV-induced progressive liver cirrhosis is a well-known risk factor for the development of HCC. Of note, HCC can occur more than 10 years after eradication of HCV, with an annual rate of 1% [9]. Liver cirrhosis is an established risk factor for HCC; it represents the final stage of liver fibrosis and usually develops in response to chronic liver injuries [10]. Chronic alcohol consumption and consequent liver cirrhosis play a causative role in the development of HCC. The consumption of contaminated animal and plant products may expose individuals to aflatoxins, another common risk factor for the development of HCC. Aflatoxin B1 (AFB1) is the most potent liver carcinogen among the four aflatoxins (i.e., B1, B2, G1, and G2). p53 gene mutations are associated with high exposure to AFB1. These mutations, such as codon 249 transversion, appeared in 50% of HCC cases [11]. NAFLD encompasses a spectrum of pathological changes characterized by different degrees of fat accumulation in the hepatocytes. This condition is attributable to overnutrition and is strongly associated with metabolic syndrome. Nonalcoholic steatohepatitis (NASH) is a severe subtype of NAFLD, with the histologic features of lobular inflammation and hepatocyte ballooning. Patients with NASH are predisposed to liver fibrosis, cirrhosis, and HCC [12]. Several mechanisms, including increased levels of tumor necrosis factor-α (TNF-α), interleukin 6 (IL-6), and leptin, have been correlated with carcinogenesis from NASH.

Treatments
In the past, HCC was usually diagnosed at an advanced stage, following the development of symptoms and impairment of liver function. At that point, treatment was often futile with poor median survival rates (i.e., <3 months) [13]. Currently, a substantial proportion of HCC patients continues to have a poor liver reserve and/or compromised portal vein flow. Thus, these untreated cases of HCC are associated with poor prognosis. With advancements in early HCC detection technology and surveillance programs, the curative treatment has improved the five-year survival rates, ranging from 50% to 75% [14]. Despite the availability of several therapeutic options for HCC (i.e., hepatic resection, liver transplantation, locoregional therapies, and systemic therapies), the treatment strategy must be individualized for each patient. The Barcelona Clinic Liver Cancer staging system is widely used worldwide to establish the prognosis and most appropriate treatment strategy for patients at different stages [15]. Although the so-called curative treatments (i.e., surgical resection, liver transplantation, and radiofrequency ablation) have greatly improved the outcomes of HCC, disease recurrence and intrahepatic metastasis continue to pose challenges in the treatment of these patients. The Barcelona Clinic Liver Cancer algorithm suggests systemic treatment for advanced HCC. Sorafenib-an oral multikinase inhibitor of cell proliferation through a strong inhibition of the serine/threonine kinase RAF-is the first approved systemic medication for the treatment of advanced HCC [16,17]. The efficacy of sorafenib has been demonstrated in several clinical studies [18]. However, although sorafenib is currently considered the best option for treating advanced HCC, it only increases life expectancy by a few months. In the recent years, certain limitations of radiotherapy for HCC treatment, such as inability to deliver a tumoricidal dose, have been overcome. In patients who are not candidates for orthotopic liver transplantation or resection, the tumor can be precisely targeted to deliver the appropriate dose through modern liver-directed radiotherapy, including three-dimensional conformal radiotherapy, charged particle radiotherapy, and stereotactic body radiotherapy [19]. HCC develops several immunosuppressive mechanisms to evade the immunological surveillance system and progress further. The immune checkpoint regulation and its associated molecules have led to advancements in cancer treatment. Among the different types of molecules involved in the immune checkpoints, programmed death-ligand 1 (PD-L1), found on the surface of cancer cells and stromal cells; programmed cell death 1 (PD-1); and cytotoxic T lymphocyte-associated protein 4, found on the surface of T cells have been shown to participate in the crucial steps for the suppression of T-cell function by cancer cells [20,21]. Active efforts toward immunotherapy for HCC include the development of monoclonal antibodies against molecules of the immune checkpoint [22].

Homeobox Genes
Homeobox genes, which are master regulatory genes controlling the development of each segment, were firstly discovered in the fruit fly Drosophila melanogaster [23]. In flies, three Homeobox genes (Ubx, Abd-A, and Abd-B) belong to the bithorax complex and the other five Homeobox genes (Lab, Pb, Dfd, Scr, and Antp) belong to the antennapedia complex. Homeobox genes specify the regional identity from the anterior to the posterior of the body segments of the fly [24]. The concept that a cluster of regulatory genes may control the development of segments is conserved in a wide range of organisms (from Caenorhabditis elegans to humans) [25]. In humans, the Homeobox genes are classified into four clusters (i.e., A to D) and located at 7p15, 17q21.2, 12q13, and 2q31. On the basis of sequence similarity and relative position of the cluster, each cluster is identified with 13 paralog groups with 9-11 protein-coding genes [26]. The Homeobox gene has a consensus element with 180 bps, and its sequence encodes approximately 60 amino acids conserved homeodomain, which is a DNA-binding globular domain. This domain contains a stable three-helix bundle preceded by a flexible N-terminal arm [27]. On the basis of the evolutionary tree, Homeobox genes can be classified into 11 classes (i.e., ANTP, PRD, LIM, POU, HF, SINE, TALE, CUT, PROS, ZF, and CERS), containing approximately 250 Homeobox genes in humans [26]. The increasing numbers of Homeobox genes are as a result of two extra rounds of genome duplication, and subsequent loss of paralogs has been discovered in humans [23]. The Homeobox genes are involved in various processes ranging from the earliest stages of development, embryonic stem cells [28], to patterning (particularly the Homeobox genes). Notably, mutations in the Homeobox genes cause developmental defects [29]. An increasing number of studies have demonstrated the aberrant expression of Homeobox genes in tumorigenesis, suggesting that, besides developmental regulation, these genes play a critical role in the development of cancer [30]. In this review, we will focus on the role of the Homeobox genes in the generation of HCC, including tumor-initiating stem-like cells (TICs), epithelial to mesenchymal transition (EMT), immunotolerance, and viral infection (Figure 1), and list the promoting or repressing function of Homeobox genes in HCC (Tables 1 and 2).

Homeobox Genes
Homeobox genes, which are master regulatory genes controlling the development of each segment, were firstly discovered in the fruit fly Drosophila melanogaster [23]. In flies, three Homeobox genes (Ubx, Abd-A, and Abd-B) belong to the bithorax complex and the other five Homeobox genes (Lab, Pb, Dfd, Scr, and Antp) belong to the antennapedia complex. Homeobox genes specify the regional identity from the anterior to the posterior of the body segments of the fly [24]. The concept that a cluster of regulatory genes may control the development of segments is conserved in a wide range of organisms (from Caenorhabditis elegans to humans) [25]. In humans, the Homeobox genes are classified into four clusters (i.e., A to D) and located at 7p15, 17q21.2, 12q13, and 2q31. On the basis of sequence similarity and relative position of the cluster, each cluster is identified with 13 paralog groups with 9-11 protein-coding genes [26]. The Homeobox gene has a consensus element with 180 bps, and its sequence encodes approximately 60 amino acids conserved homeodomain, which is a DNA-binding globular domain. This domain contains a stable three-helix bundle preceded by a flexible N-terminal arm [27]. On the basis of the evolutionary tree, Homeobox genes can be classified into 11 classes (i.e., ANTP, PRD, LIM, POU, HF, SINE, TALE, CUT, PROS, ZF, and CERS), containing approximately 250 Homeobox genes in humans [26]. The increasing numbers of Homeobox genes are as a result of two extra rounds of genome duplication, and subsequent loss of paralogs has been discovered in humans [23]. The Homeobox genes are involved in various processes ranging from the earliest stages of development, embryonic stem cells [28], to patterning (particularly the Homeobox genes). Notably, mutations in the Homeobox genes cause developmental defects [29]. An increasing number of studies have demonstrated the aberrant expression of Homeobox genes in tumorigenesis, suggesting that, besides developmental regulation, these genes play a critical role in the development of cancer [30]. In this review, we will focus on the role of the Homeobox genes in the generation of HCC, including tumor-initiating stem-like cells (TICs), epithelial to mesenchymal transition (EMT), immunotolerance, and viral infection (Figure 1), and list the promoting or repressing function of Homeobox genes in HCC ( Table 1; Table 2).   EMT-epithelial to mesenchymal transition.

Homeobox Genes (Oct4 and Nanog) in TICs of HCC
TICs have been known to be capable of recapitulating the primary tumor heterogeneity owing to their self-renewal, differentiation, and tumor-initiation capacities [98,99]. One of the main properties of TICs is their self-renewal ability, which has been identified through performing sphere formation assays in vitro and through transplantation experiments with limiting dilution in mice [100,101]. Several cell-surface antigens have been discovered on liver TICs, including PROM1 (CD133), THY1 (CD90), epithelial cell adhesion molecule (EpCAM), and CD24 [102]. Increasing evidence indicates that liver TICs play a critical role in hepatocarcinogenesis [66,103]. In Homeobox genes, octamer-binding transcription factor 4 (Oct4) and Nanog reportedly control the self-renewal and pluripotency of pluripotent stem cells [104,105]. Moreover, both of them were identified as critical non-cell-surface markers in TICs of various types of tumors [106,107]. Therefore, it is essential to investigate the mechanism through which the Homeobox genes are regulated in TICs. The elucidation of this mechanism may assist treating physicians in achieving improved clinical outcomes for patients with HCC. MiR-429 has been identified as an oncogene in HCC by promoting the self-renewal, tumorigenicity, and chemoresistance of EpCAM-positive TICs. MiR-429 promotes the transcription of Oct4 by reducing the expression of retinoblastoma-binding protein 4, which otherwise inhibits the E2F1 transactivation of Oct4 transcription [68]. The HCC cells with high Nanog expression, which have stem-cell-like characteristics, are characterized by a high capacity for tumor invasion and metastasis and are resistant to treatment with sorafenib and cisplatin [62]. These cells are required for Nanog to promote the expression of insulin-like growth factor (IGF) 2 and IGF 1 receptor (IGF1R) for self-renewal [62], and the NODAL/SMAD3 signaling pathway for tumor invasion [63]. It has been reported that EMT-associated transcription factors Snail and Slug or the androgen/androgen receptor axis directly regulate the expression of Nanog under transforming growth factor-β (TGF-β) signaling in TICs [58,59,64]. Interestingly, chromatin immunoprecipitation-sequencing analyses of Nanog showed that it regulates the expression of mitochondrial metabolic genes required in CD133-positive TICs. Moreover, it represses mitochondrial oxidative phosphorylation genes to prevent the reactive oxygen species (ROS) induction and induce oxidation of fatty acid for self-renewal and drug resistance of TICs [66]. Nevertheless, in murine nonalcoholic steatohepatitis model, Nanog contribute to the reprogramming of hepatic progenitor cells in driving these cells to TICs [60].
Increasing studies have also focused on signaling pathways that can manipulate the Nanog expression to alter the properties of TICs. Receptor for activated C kinase1 (RACK1) was first discovered through interactions between the activated form of protein kinase C, which positively regulates self-renewal and chemoresistance of CD13, and CD133 double-positive HCC cells by directly binding with Nanog. Its binding abrogated the recruitment of E3 ubiquitin ligase FBXW8 and degradation of proteasome to promote stemness and chemoresistance [108]. Myeloid cell leukemia-1 (Mcl-1) promotes HCC stemness and self-renewal by regulating the level of Nanog, Sox2, and KLF4 [109]. Transforming acidic coiled-coil protein 3 (TACC3), which is involved in cell mitosis and transcriptional activity, is required for Wnt/β-catenin and PI3K/AKT signaling pathways to maintain stem cell transcription factors, including Bmi1, c-Myc, and Nanog [110]. In addition, negative regulators of Nanog expression in TICs have been identified; for example, the GATA transcription factor 5 (GATA5) negatively regulates the expression of Oct4, Nanog, Klf4, c-myc, and EpCAM in TICs by binding with β-catenin and inhibiting β-catenin nuclear translocation [111]. Levels of Atonal homolog 8 (ATOH8), a basic-helix-loop-helix (bHLH) transcription factor, are reduced in HCC patients, and loss of ATOH8 increases the transcription of OCT4, NANOG, and CD133 and reduces chemo-sensitivity to 5-fluorouracil and cisplatin [112].

Homeobox in EMT of HCC
EMT is described as a conversion process of adherent epithelial cells into migratory mesenchymal cells with highly invasive properties. In HCC, EMT are endowed with HCC progression by inhibiting apoptosis and senescence, escaping immune reactions, and gaining chemoresistance [113][114][115], and are also involved in the early stages of tumor transformation into aggressive malignancies by increasing the potential for invasiveness and metastasis. Typically, cells undergoing EMT are characterized by decreased levels of E-cadherin, increased levels of N-cadherin and vimentin, and translocation of β-catenin from the membrane to the nucleus. An increasing number of studies have shown that Homeobox genes can positively or negatively regulate EMT in HCC as follows.
Zinc finger E-box binding Hox 1 (ZEB1) and ZEB2 are master transcription factors that positively regulate invasion and metastasis by promoting EMT in cancer cells [116]. Both ZEB1 and ZEB2 comprise two zinc finger domains at the N-and C-terminals and the Smad interaction domain, homeodomain, and CtBP interaction domain at the central region. In breast cancer, ZEB1 and ZEB2 have been shown to directly bind to the E-box located in the E-cadherin promoter [117], recruiting the CtBP transcriptional co-repressors [118] and/or the SWI/SNF chromatin remodeling protein BRG1 [119]. This results in repression of E-cadherin expression and promotion of EMT. In HCC patients, ZEB1 expression has been associated with low expression of E-cadherin, venous invasion, and tumor/node/metastasis (TNM)stage. Patients expressing high levels of ZEB1 and low levels of E-cadherin are associated with poorer prognosis [78].
In recent years, several regulations have been identified upstream of ZEB1 or ZEB2. Depletion of thrombomodulin, a natural anticoagulation factor, induces HCC cell migration by induction of ZEB1 and reduction of E-cadherin [76]. Moreover, 14-3-3e-a protein belonging to the highly conserved in eukaryotic cells 14-3-3 protein family-suppresses the expression of E-cadherin via regulation of ZEB1 [77]. It has been shown that this regulation increases the risk of metastasis and decreases the survival rates in HCC patients. Previous studies have shown that MYC-associated zinc finger protein (MAZ) promotes the expression of c-Myc, Ras, VEGF, and podoplanin and represses that of p53, Sp4, and endothelial nitric oxide synthase in tumor development. MAZ was also shown to promote EMT and metastasis in HCC patients [85]. Depletion of the minus-end-directed motor protein kinesin family member C1 [87]-also termed HSET-downregulates ZEB1 to reduce EMT in HCC. Numerous studies have shown that mutations of the liver kinase B1 cause cancer (i.e., lung, mammary gland, ovarian, melanoma, and HCC). In addition, it is involved in liver carcinogenesis by promoting ZEB1 expression-associated EMT [89]. Forkhead box Q1, a forkhead transcription factor, induces EMT in HCC by up-regulating the expression of ZEB2 [92].
Regulatory proteins can manipulate the levels of ZEB1 and ZEB2 in EMT. However, there are numerous microRNAs (miRNAs) and long non-coding RNAs (lncRNAs) that participate in ZEB1-mediated EMT [120,121]. For example, miR-139-5p, which suppresses EMT in HCC by binding to ZEB1 and ZEB2 mRNA, was downregulated in HCC tissue [94]. In addition, MiR-101 binds to the 3 -untranslated region of ZEB1 to silence and disrupt the EMT in HCC. However, this regulation can be reversed by small nucleolar host gene 6 transcript, acting as a competing endogenous RNA [81,84]. MiR-211-5p and miR-154 targeting ZEB2 suppress HCC metastasis and tumor growth [93,97]. Gα12 gep is an oncogene that generates G-protein-coupled receptors sensing the increasing levels of ligands in tumor microenvironments. The activated Gα12 promotes EMT through the induction of ZEB1 and reduction of miRNAs (i.e., miR-192, miR-215, and miR-200a), which target ZEB1 and ZEB2 [80,90]. Long non-coding RNA activated by TGF-β (lncRNA-ATB) promotes EMT under TGF-β stimulation by binding with the miR-200 family, which directly targets ZEB1 and ZEB2 mRNA [79]. Another lncRNA ZFAS1 also up-regulated the expression of ZEB1 by binding competitively to miR-150, consequently inducing EMT and invasion [82]. An increasing number of studies have shown that upstream antisense transcription controls the transcription of the corresponding genes [122]. Interestingly, a non-coding antisense transcript, ZEB1 antisense 1, is located in physical contiguity with ZEB1. ZEB1 antisense 1 promotes HCC invasion, metastasis, and EMT by targeting ZEB1 [83]. A similar regulation mechanism has also been found between ZEB2 and ZEB2 antisense 1 [95]. Recent studies showed that lncRNAs competitively bind with miRNAs that maintain the expression of ZEB1 and EMT, including lncRNA MALAT1 and miR-143-3p, lncRNA TUG1, and Mir-142-3p [86,88].

Homeobox Genes and Immunotolerance in HCC
Liver sinusoids are involved in the central immunological functions of the liver [123]. With this system, the liver has the capacity to remove different microbes, microbe-associated molecules, and DAMPs, which continuously circulate from the gut to the liver. The diverse innate and adaptive immune cells residing in the liver to detect and clear bloodborne infections include Kupffer, natural killer, natural killer T, and CD4+/8+ T cells. These immune cells respond to large and diverse cell-surface ligands expressed by infected, damaged, or transformed cells, leading to the changes in the innate and adaptive immune reactions through the production of different potent cytokines [124]. Sinusoidal endothelial cells, stellate cells, and hepatocytes participate in the maintenance of balance between immunotolerance and immune activation to prevent liver damage from non-pathological or continuous inflammatory stimuli or systemic immunotolerance [125]. Deregulation of the precisely regulated immunological network promotes the development of HCC, possibly owing to chronic infection (e.g., infection with HBV or HCV), fat accumulation (i.e., NASH), or DAMPs generated as a result of toxic liver damage (alcoholic liver disease) [126].
Rare reports focus on Homeobox gene and immune suppression in tumor cells. One example is that of the intestine-specific Homeobox transcription factor (ISX), which has a consensus homeodomain (Figure 1), and is a proto-oncogene involved in HCC development [55][56][57]. Recently, our lab demonstrated that ISX is involved in a positive feedback loop, including inflammation, tryptophan catabolism, and immune suppression. IL-6 induces the transcriptional activation of ISX to promote the production of tryptophan catabolic enzymes tryptophan 2,3-dioxygenase and indoleamine 2,3-dioxygenase 1 in HCC. Both enzymes increase the level of tryptophan catabolite, kynurenine, and aryl hydrocarbon receptor (AhR) and activate the kynurenine/AHR axis. Its activation promotes a positive feedback mechanism to increase ISX associated proliferation, tumorigenesis, and immunotolerance. Apart from the AHR-dependent immunotolerance, overexpression of ISX induced level of genes encoding the immune modulators CD86 (B7-2) and PD-L1 and presented a repressed CD8 + T-cell response.

Homeobox Genes in HBV-and HCV-Associated HCC
Clinical and epidemiological studies have linked chronic hepatic inflammation to the pathogenesis of HCC [127]. The chronic HBV-and HCV-induced inflammation contributes to the development of HCC [128]. Notably, this tumorigenesis was associated with the regulation of numerous Homeobox genes [30]. The Hox-A13 protein is induced in liver stem-like cells and HBV or HCV-infected HCC, but not in hepatocytes and bile duct epithelia [129]. The function of the prospero-related Hox protein 1 in the progression of HCC is debatable [73,[130][131][132][133]. However, it represses HBV antigen expression and genome replication through repression of the enhancer II/core promoter, preS1 promoter, and enhancer I/X promoter of HBV [134]. In HBV-infected HCC cells, IL-6 increases IGF1R and result in the stemness-related properties that evaluate the Oct4/Nanog, which confers poor prognosis [69]. Furthermore, the HBx encoded by the HBV X gene has been shown to be associated with HCC development. The HBX with C-terminally truncation causes more malignant HCC by promoting metastasis and tumorigenicity. The reason is that HBx-∆C1 is involved in the regulation of the properties of liver cancer stem cells by up-regulating the expression of Homeobox and NANOG through the stat3 pathway [135]. The regulation of HCV by Homeobox promotes TICs. Overexpression of HCV non-structural protein NS5A and induction Toll-like receptor 4 (TLR4) by alcohol-induced endotoxemia in hepatocyte synergistically generates liver damage and tumor development, resulting in the generation of Nanog-positive TICs from the mice model [65,136]. Of note, TLR4/Nanog-dependent TICs are also found in HCC patients. Research has identified a TLR4/Nanog-mediated activation of YAP1 and IGF2BP3, which are novel molecules responsible for the inhibition of the TGF-β pathway and the development of chemoresistance [102]. Using a high-cholesterol/high-fat diet, it was shown that a higher proportion of hepatocyte-specific NS5A transgenic mice developed liver tumors containing TICs. This was attributed to the activation of TLR4-Nanog and pSTAT3 signaling pathways through an exaggerated EMT via the induction of Twist1 [65].

Conclusions
Dysregulated expression of Homeobox genes is wildly identified in different aspects of HCC development, including HBV and HCV infection, TICs, EMT, and immunotolerance. Homeobox genes are discovered in both positive and negative regulation in HCC progression; however, transcriptional regulatory networks of Homeobox genes in HCC remain unclear. Integrating single cell RNA sequencing in HCC to build a systemic interaction network of Homeobox genes in HCC is noteworthy. It can provide more information for the cell type specific mechanism of the regulation of Homeobox genes in the tumor microenvironment and for the design of more efficacious therapy of HCC patients.
In recent years, Homeobox genes have been identified at different stages of the progression of HCC; thus, Homeobox genes are the potential therapeutic targets, but some questions remain to be answered. One of the strategies is regulating the expression or stability of Homeobox genes, but more effort is required to identify the cohesive signaling regulation for using a specific inhibitor and activator or to develop the drug for specific targets. Second, the structure of the homeodomain is not suitable for the development of a small molecular drug to block the DNA-binding ability of Homeobox. The conservation of homeodomain increases the difficulties regarding the specificity of drug and needs to be solved by structure identification. Third, it has been known that the expressions of Homeobox genes in HCC patients are diverse; thus, development of a high-throughput screening method for detecting aberrant expressions of Homeobox genes would be valuable to develop a precise medicine tailored to individual needs. Therefore, more and more basic and clinical studies are required to enhance the knowledge regarding Homeobox genes in HCC patients.

Conflicts of Interest:
The authors declare no conflict of interest.