Huaping Chen1Junrong Wu3Liuyi Lu1Zuojian Hu1Xi Li1Li Huang1Xiaolian Zhang1Mingxing Chen1
Xue Qin1*Li Xie2*- 1Department of Clinical Laboratory, First Affiliated Hospital of Guangxi Medical University, Nanning, China
- 2Department of Clinical Laboratory, Second Affiliated Hospital of Guangxi Medical University, Nanning, China
- 3Department of Clinical Laboratory, Affiliated Tumor Hospital of Guangxi Medical University, Nanning, China
Aims: In the cancer-related research field, there is currently a major need for a greater number of valuable biomarkers to predict the prognosis of hepatocellular carcinoma (HCC). In this study, we aimed to screen hub genes related to immune cell infiltration and explore their prognostic value for HCC.
Methods: We analyzed five datasets (GSE46408, GSE57957, GSE74656, GSE76427, and GSE87630) from the Gene Expression Omnibus database to screen the differentially expressed genes (DEGs). A protein–protein interaction network of the DEGs was constructed using the Search Tool for the Retrieval of Interacting Genes; then, the hub genes were identified. Functional enrichment of the genes was performed on the Metascape website. Next, the expression of these hub genes was validated in several databases, including Oncomine, Gene Expression Profiling Interactive Analysis 2 (GEPIA2), and Human Protein Atlas. We explored the correlations between the hub genes and infiltrated immune cells in the TIMER2.0 database. The survival curves were generated in GEPIA2, and the univariate and multivariate Cox regression analyses were performed using TIMER2.0.
Results: The top ten hub genes [DNA topoisomerase II alpha (TOP2A), cyclin B2 (CCNB2), protein regulator of cytokinesis 1 (PRC1), Rac GTPase-activating protein 1 (RACGAP1), aurora kinase A (AURKA), cyclin-dependent kinase inhibitor 3 (CDKN3), nucleolar and spindle-associated protein 1 (NUSAP1), cell division cycle-associated 5 (CDCA5), abnormal spindle microtubule assembly (ASPM), and non-SMC condensin I complex subunit G (NCAPG)] were identified in subsequent analysis. These genes are most markedly enriched in cell division, suggesting their close association with tumorigenesis. Multi-database analyses validated that the hub genes were upregulated in HCC tissues. All hub genes positively correlated with several types of immune infiltration, including B cells, CD4+ T cells, macrophages, and dendritic cells. Furthermore, these hub genes served as independent prognostic factors, and the expression of these hub genes combing with the macrophage levels could help predict an unfavorable prognosis of HCC.
Conclusion: In sum, these hub genes (TOP2A, CCNB2, PRC1, RACGAP1, AURKA, CDKN3, NUSAP1, CDCA5, ASPM, and NCAPG) may be pivotal markers for prognostic prediction as well as potentially work as targets for immune-based intervention strategies in HCC.
Introduction
Hepatocellular carcinoma (HCC), the second leading cause of cancer-related death in the world, is a commonly fatal cancer with an unfavorable prognosis due to its complex genetics and clinical features (Siegel et al., 2015; Villanueva, 2019). To a certain extent, high heterogeneity contributes to the low survival rate of HCC despite excision, transplantation, transcatheter arterial chemoembolization, and radiofrequency ablation, among others, having been widely used in HCC treatment (Imamura et al., 2003; Forner et al., 2018; Caruso and Nault, 2019; Yin et al., 2019). Timely and effective intervention for HCC patients can improve not only their quality of life but also their survival rate (Deng et al., 2020; Kim et al., 2020). Therefore, the identification of new prognostic biomarkers and therapeutic targets plays a crucial role in HCC therapy.
Several prognostic biomarkers have been widely applied in HCC, such as alpha-fetoprotein and des-gamma-carboxyprothrombin (Fox et al., 2014; Abe et al., 2017). Nevertheless, these markers depend on the significant burthen of a tumor, which has resulted in their often limited application and inconsistent performance assessments (Park and Park, 2013). Numerous studies have demonstrated that valuable prognostic molecules can be identified from the bioinformatics analysis of high-throughput data, such as functional genes (Ma et al., 2020; Wang et al., 2020; Xie et al., 2020). From this, it can be determined that immune-associated genes may play a crucial role in HCC outcomes and targeted therapies on immune cells; thus, related genes have been developed for the reactivation of adaptive and innate immune systems and the creation of a strong antitumoral immune response. For instance, some researchers have found that inhibitors of programed death-1, programed death-ligand 1, and cytotoxic T lymphocyte-associated antigen 4 produce antitumoral effects on HCC cells (Langhans et al., 2019; Wang J. et al., 2019; Zhang Y. et al., 2020). Unfortunately, only sectional HCC patients with determinate tumor types react to the current immunotherapies; therefore, there is an urgent need to identify more underlying immune targets.
In the present study, differentially expressed genes (DEGs) were screened from the Gene Expression Omnibus (GEO) database. Then, we performed enrichment and protein–protein interaction (PPI) analyses of these genes to comprehend the functions of DEGs and identify the top ten hub genes in HCC. Next, we explored the potential correlations between each of the hub genes and infiltrated immune cells in the TIMER2.0 database. Furthermore, we visualized the prognostic landscape of candidate hub genes using several databases, including Oncomine (Rhodes et al., 2007), Gene Expression Profiling Interactive Analysis 2 (GEPIA2) (Tang et al., 2019), Human Protein Atlas (HPA) (Ponten et al., 2008), and TIMER2.0 (Li et al., 2020).
Materials and Methods
Data Source
The gene expression datasets (GSE46408, GSE57957, GSE74656, GSE76427, and GSE87630) of HCC were obtained from GEO1. All of the data included in the present study was available online. Information on these five datasets is summarized in Table 1.
Table 1. Details of GEO HCC data.
DEGs Processing
GEO2R2, an interactive online tool that can compare two different groups in a GEO dataset, was applied to screen the DEGs (Davis and Meltzer, 2007). The DEGs were defined as different expression genes between tumor and tumor-adjacent tissues in HCC patients with an adjusted p value < 0.05 and an absolute log fold-change (FC) > 1. Accordingly, to decrease the false discovery rate, the p values were adjusted using the Benjamini and Hochberg method. The overlapping up- and downregulated DEGs from these five datasets were identified using TBtools software (Chen et al., 2020).
Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) Pathway Analyses
Functional enrichment analyses played a crucial role in the identification of biological characteristics in transcriptome data. In this study, the Metascape database was used to conduct the Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses in overlapping DEGs and identified hub genes (Zhou Y. et al., 2019). In order to choose the more remarkable terms within each cluster, a p value < 0.01 was considered statistically significant. We used bar charts and bubble diagrams, respectively, to visualize the results of the GO and KEGG analyses.
PPI Network Construction and Module Analysis
The Search Tool for the Retrieval of Interacting Genes (STRING) is an interacting gene database designed to analyze PPI information (Szklarczyk et al., 2019). The overlapping DEGs were mapped in STRING to generate a network with functional interactions; then, this PPI network was visualized using Cytoscape software. Next, we employed the Molecular Complex Detection (MCODE) plugin to determine the most significant gene modules. Moreover, the cytoHubba plugin was applied to identify the hub genes using the Maximal Clique Centrality (MCC) method, and the top ten hub genes with the highest MCC scores were subjected to the subsequent analyses.
Validation for mRNA and Protein Levels of Hub Genes in Multi-Databases
The Oncomine database3 is a publicly available cancer database that facilitates the analysis of genome-wide expression in multifarious cancers. In the present study, the Oncomine was utilized to overview the mRNA expression of candidate genes with a p value < 0.0001 and |FC| > 1.5.
The GEPIA2 database4, which includes TCGA and GTEx data, was applied to analyze the differential expression of the hub genes in the HCC and normal groups, and the cutoff values were set as |log2FC| = 1.0 and p value = 0.01.
Furthermore, immunohistochemistry analysis was conducted online to observe the distribution and protein level of the candidate hub genes in the HPA database5.
Survival Analysis of Hub Genes
We employed GEPIA2 to perform overall survival (OS) and relapse-free survival (RFS) analyses, with a median group cutoff in 362 HCC patients. The survival curves with the calculated hazard rate (HR) and the log-rank p value were presented on the charts. Additionally, TIMER2.0 was utilized to conduct a univariate Cox analysis to validate the results from the GEPIA2 analyses.
TIMER2.0 Database Analysis
TIMER2.06 is a comprehensive source that can explore the relationship between two genes, and the correlations of mRNA expression and immune infiltration. In this study, we analyzed the correlations of hub gene expression and several tumor-associated immune cells, including B cells, CD8+ T cells, CD4+ T cells, macrophages, neutrophils, and dendritic cells (DCs). A p value less than 0.01 was considered statistically significant to identify the more prominent correlation between the hub genes and immune cells.
Further OS analyses were performed with macrophage and single hub gene expression. Moreover, we constructed ten multivariate Cox proportional hazard models, each of which comprised seven variables, including age, tumor stage, gender, race, tumor purity, macrophage level, and expression of the single candidate hub gene. The survival curves, featuring patterns of single gene expression, and macrophage levels were shown on the diagram. The association between each macrophage and OS was displayed as the low or high expression of a single hub gene.
Statistical Analysis
The Student’s t test or non-parametric Mann–Whitney test was utilized to compare the two independent groups, as appropriate. The correlations between the candidate hub genes as well as the relationship of these genes and immune cell infiltration were assessed using Spearman’s correlation. The log-rank test was used to calculate the HR and log-rank p value to compare the survival curves. Univariate and multivariate Cox regression models were applied to calculate the HR and Cox p value. If not specifically stated, p values < 0.05 were considered statistically significant.
Results
Identification and Enrichment Analysis of DEGs
After screening the DEGs according to the criteria, 1,414, 417, 454, 493, and 1,163 genes were identified from GSE46408, GSE57957, GSE74656, GSE76427, and GSE87630, respectively (Table 1). These genes shared 107 DEGs among these five datasets, of which, there were 18 upregulated and 89 downregulated genes (Figures 1A,B and Supplementary Figure 1).
Figure 1. Identification of the DEGs between liver tumor and non-cancerous tissues in GSE46408, GSE57957, GSE74656, GSE76427, and GSE87630. Venn diagram of (A) upregulated and (B) downregulated DEGs based on the five GEO datasets. The overlapping areas represent the commonly altered DEGs. The t-test was used to analyze DEGs, with the cutoff criteria of |log FC| > 1.0 and adjusted p < 0.05. DEGs, differentially expressed genes; GEO, Gene Expression Omnibus; log FC, log fold change.
As shown in Figure 2A, the results of the GO analysis suggested that the overlapping DEGs were principally enriched in the biological process (BP), especially in the detoxification of copper ion, the organic acid catabolic process, and the monocarboxylic acid metabolic process. In terms of cellular component (CC), the DEGs were enriched in plasma lipoprotein particles and extracellular matrix. In regard to molecular function (MF), the identified DEGs were significantly enriched in the aspects of peptide hormone binding, insulin-like growth factor binding, and carboxylic acid binding. Also, KEGG analysis showed that the overlapping DEGs were dramatically concentrated in mineral absorption and tryptophan metabolism (Figure 2B).
Figure 2. GO annotation and KEGG pathway enrichment analysis of DEGs. The significantly enriched (A) GO terms and (B) KEGG pathways. GO, gene ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes; DEGs, differentially expressed genes; BP, biological process; CC, cellular component; MF, molecular function. GO:0010273, detoxification of copper ion; GO:0016054, organic acid catabolic process; GO:0032787, monocarboxylic acid metabolic; GO:0034754, cellular hormone metabolic process; GO:1901615, organic hydroxy compound metabolic process; GO:0071466, cellular response to xenobiotic stimulus; GO:0001867, complement activation, lectin pathway; GO:0008209, androgen metabolic process; GO:0006544, glycine metabolic process; GO:0007076, mitotic chromosome condensation; GO:0009617, response to bacterium; GO:0048589, developmental growth; GO:0044262, cellular carbohydrate metabolic process; GO:0031214, biomineral tissue development; GO:0034358, plasma lipoprotein particle; GO:0031012, extracellular matrix; GO:0017046, peptide hormone binding; GO:0005520, insulin-like growth factor binding; GO:0031406, carboxylic acid-binding.
PPI Network Construction and Hub Gene Identification
To seek the interactions of the overlapping DEGs, a PPI network, which included 84 nodes and 191 edges, was constructed and visualized in the Cytoscape (Figure 3A). As Figures 3B–D shows, the three most prominent subnetworks were identified using the MCODE plugin according to the connective degrees. Moreover, the CytoHubba plugin was used to determine the top ten hub genes based on their MCC scores (Figure 3E). Interestingly, the hub-gene network from the CytoHubba analysis completely coincided with the highest score module from the MCODE analysis, both of which included the same 10 nodes and 45 edges. Notably, DNA topoisomerase II alpha (TOP2A) was the most significant gene, with the highest MCC score of 362,912, followed by cyclin B2 (CCNB2) (MCC score = 362,906), protein regulator of cytokinesis 1 (PRC1) (MCC score = 362,904), Rac GTPase-activating protein 1 (RACGAP1) (MCC score = 362,904), aurora kinase A, (AURKA) (MCC score = 362,888), cyclin-dependent kinase inhibitor 3 (CDKN3) (MCC score = 362,886), nucleolar and spindle-associated protein 1 (NUSAP1) (MCC score = 362,880), cell division cycle-associated 5 (CDCA5) (MCC score = 362,880), abnormal spindle microtubule assembly (ASPM) (MCC score = 362,880), and non-SMC condensin I complex subunit G (NCAPG) (MCC score = 362,880). All the hub genes were upregulated in HCC tissues. In addition, the analyses of the correlations between candidate hub genes on the mRNA level were conducted on the TIMER2.0. Of these ten genes, every two genes showed a significant correlation (Spearman’s rho value > 0.7; p < 0.05) with or without tumor purity adjustment (Supplementary Table 1).
Figure 3. PPI network and clusters identification. (A) The PPI network of DEGs was constructed using Cytoscape. Upregulated genes are marked in red, and downregulated genes are marked in green. (B) Cluster 1 with 10 nodes and 45 edges. (C) Cluster 2 with 7 nodes and 18 edges. (D) Cluster 3 with 8 nodes and 14 edges. MCODE plugin was employed for the detection of clusters. Upregulated genes are marked in red, and downregulated genes are marked in green. (E) The top ten hub genes were identified using CytoHubba and ranked by the MCC score. The hub-gene network from CytoHubba analysis completely coincided with the Cluster 1 from MCODE analysis, both of which included the same 10 nodes and 45 edges. PPI, protein–protein interaction; DEGs, differentially expressed genes; MCODE, molecular complex detection; MCC, maximal clique centrality.
The high scores of these candidate hub genes indicated that there would be potential biological effects in the hub genes network; thus, we further determined the functional enrichment of these genes. The BP analysis proved that the hub genes were dramatically enriched in terms of cell division, the positive regulation of mitotic nuclear division, female gamete generation, and mitotic cell cycle phase transition. Besides, these genes were significantly enriched in terms of the mitotic spindle in the CC analysis and in terms of protein kinase binding in the MP analysis (Supplementary Figure 2). These results implied that these hub genes are closely associated with tumorigenesis.
The Expression of Hub Genes Was Upregulated in Multi-Databases
To verify the dependability of the results from the bioinformatics analysis, we next determined the mRNA expression levels of the hub genes in the Oncomine and GEPIA2 databases. The results showed that the transcriptional levels of TOP2A, CCNB2, PRC1, RACGAP1, AURKA, CDKN3, NUSAP1, CDCA5, ASPM, and NCAPG were significantly overexpressed in HCC tissue when compared with the normal controls (Figures 4A–K), indicating their potential oncogenic effects.
Figure 4. The expression of the hub genes was validated in HCC tissues. (A) The mRNA expression patterns of hub genes in the Oncomine database. This figure shows the numbers of datasets with statistically significant mRNA overexpression (red) or under-expression (blue) of the target gene (cancer vs. normal). The p value threshold is 0.0001. The number in each cell represents the number of analyses that met the threshold within the analyses and cancer types. The mRNA expression levels of (B) TOP2A, (C) CCNB2, (D) PRC1, (E) RACGAP1, (F) AURKA, (G) CDKN3, (H) NUSAP1, (I) CDCA5, (J) ASPM, and (K) NCAPG in LIHC tissues and normal liver tissues using GEPIA2. *p < 0.01 was considered statistically significant. HCC, hepatocellular carcinoma; LIHC, liver hepatocellular carcinoma. The “Yellow Box” specified the mRNA expression of these hub genes in liver cancer.
In the HPA database analysis, we found that the protein levels of TOP2A, PRC1, RACGAP1, AURKA, NUSAP1, and CDCA5 were significantly higher in the HCC tissues when compared to normal liver tissues (Figures 5A–F).
Figure 5. Representative immunohistochemistry images in HCC and normal liver tissues derived from the HPA database. (A) TOP2A incubating with HPA006458 antibody, (B) PRC1 incubating with HPA034521 antibody, (C) RACGAP1 incubating with HPA043912 antibody, (D) AURKA incubating with HPA002636 antibody, (E) NUSAP1 incubating with HPA042904 antibody, (F) CDCA5 incubating with HPA023691 antibody. HCC, hepatocellular carcinoma; HPA, Human Protein Atlas.
The mRNA Levels of Hub Genes Are Positively Associated With Immune Infiltration in HCC
Numerous studies have demonstrated that the infiltration of tumor-associated immune cells and immune-related genes is correlated with the development and prognosis of HCC (Duan et al., 2019; Huang et al., 2020; Tang et al., 2020). Remarkably, targeting on these immune cells and/or genes has been a prospective workable approach in HCC therapy (Jayant et al., 2020). In this study, we attempted to explore the relationship between the mRNA expression of hub genes and immune cell infiltration using TIMER2.0. As Table 2 presents, each of the hub genes correlated with tumor purity in HCC tissues. Notably, we observed that these ten genes presented significant associations with infiltrating levels of B cells, CD4+ T cells, macrophages, and DCs, out of which these genes most strongly correlated with B cells [correlation coefficient (COR), 0.363 to 0.451; p < 0.01], macrophages (COR 0.262 to 0.384; p < 0.01), and DCs (COR, 0.449 to 0.548; p < 0.01), indicating that these hub genes were positively related to tumor-associated B cells, macrophages, and DCs in the HCC microenvironment.
Table 2. Correlation analysis between candidate hub genes and immune cells in the TIMER2.0 database.
The Overexpression of Hub Genes Predicts Poor Prognosis in HCC
In the GEPIA2 analysis, the KM plotter analyses showed that the upregulated mRNA levels of TOP2A, PRC1, RACGAP1, AURKA, NUSAP1, CDCA5, ASPM, and NCAPG were correlated with worse OS in HCC patients; however, there was no significant correlation between CCNB2/CDKN3 and HCC prognosis (Figures 6A–J). For the RFS analyses, the overexpression of all ten hub genes could predict unfavorable prognosis in HCC (Figures 7A–J). Additionally, we performed a univariate Cox regression analysis of the candidate hub genes in TIMER2.0. Table 3 summarizes the validation of the prognostic values of these hub genes in the TIMER2.0 database, suggesting that each of the hub genes may be an independent risk factor in HCC.
Figure 6. Overall survival analysis of ten hub genes in HCC patients. (A) TOP2A, (B) CCNB2, (C) PRC1, (D) RACGAP1, (E) AURKA, (F) CDKN3, (G) NUSAP1, (H) CDCA5, (I) ASPM, and (J) NCAPG. Log-rank p value < 0.05 was considered statistically significant. HCC, hepatocellular carcinoma; OS, overall survival; HR, hazard ratio.
Figure 7. Disease-free survival analysis of ten hub genes in HCC patients. (A) TOP2A, (B) CCNB2, (C) PRC1, (D) RACGAP1, (E) AURKA, (F) CDKN3, (G) NUSAP1, (H) CDCA5, (I) ASPM, and (J) NCAPG. Log-rank p value < 0.05 was considered statistically significant. HCC, hepatocellular carcinoma; DFS, disease-free survival; HR, hazard ratio.
Table 3. Univariate Cox proportional hazards analyses of overall survival in the TIMER2.0.
The Overexpression of Hub Genes Accompanied by a High Level of Infiltrated Macrophages Predicts Poor Prognosis in HCC
Tumor-associated macrophages (TAMs) majorly facilitate tumor angiogenesis, invasion, and metastasis and lead to a poor prognosis in HCC (Zhang et al., 2019; Zhao et al., 2020). There have been studies substantiating that targeting TAMs becomes a potential individualized precision single or combined therapy (Li et al., 2019). Therefore, the identification of TAM-related genes will contribute to providing more potential targets of the individualized precision treatment and improve the prognosis of HCC. In the present study, we evaluated the prognostic efficiency of the combination of TAMs and expression patterns for the single hub gene. As shown in Figures 8A–J, there was no significant relationship between the TAMs and prognosis under the low expression level of TOP2A/CCNB2/PRC1/RACGAP1/AURKA/CDKN3/NUSAP1/CD CA5/ASPM/NCAPG. However, under high CCNB2 expression, higher macrophage levels had a worse outcome in HCC [HR = 1.6, p = 0.0366]. Similarly, the high macrophage level predicted unfavorable prognosis under the high expression of RACGAP1 (HR = 1.79, p = 0.0131), AURKA (HR = 1.64, p = 0.0283), CDKN3 (HR = 1.94, p = 0.0066), ASPM (HR = 1.72, p = 0.0177), and NCAPG (HR = 1.66, p = 0.0304).
Figure 8. Overall survival analysis for combining the expression of single hub gene and macrophage in HCC patients. (A) TOP2A, (B) CCNB2, (C) PRC1, (D) RACGAP1, (E) AURKA, (F) CDKN3, (G) NUSAP1, (H) CDCA5, (I) ASPM, and (J) NCAPG. P value < 0.05 was considered statistically significant. HCC, hepatocellular carcinoma; HR, hazard ratio.
Moreover, we further established ten multivariate Cox regression analyses, each of which included seven variables: macrophage level, age, stage, gender, race, tumor purity, and expression of a single candidate gene (Figures 9A–J). The results showed that, after adjustments of age, stage, gender, race, and tumor purity, there was still no statistical correlation between TAM and prognosis with the low expression of TOP2A/CCNB2/PRC1/RACGAP1/AURKA/CDKN3/NUSAP1/CD CA5/ASPM/NCAPG; nevertheless, the lower level of TAM could predict favorable prognosis under the high expression of TOP2A (HR = 2.1, p = 0.0078)/CCNB2 (HR = 2.2, p = 0.0046)/PRC1 (HR = 1.92, p = 0.0221)/RACGAP1 (HR = 2.09, p = 0.0091)/AURKA (HR = 2.08, p = 0.0061)/ CDKN3 (HR = 2.7, p = 0.0009)/NUSAP1 (HR = 2.25, p = 0.0043)/CDCA5 (HR = 2.05, p = 0.0086)/ASPM (HR = 2.18, p = 0.0053)/NCAPG (HR = 2.25, p = 0.0043) (Figures 10A–J). These results suggested that each of the hub genes was an independent unfavorable prognostic biomarker and that combining their respective expression levels with the TAMs would help them play a more effective role in the prognosis prediction of HCC.
Figure 9. The univariate and multivariate Cox regression analysis of these hub genes. (A) TOP2A, (B) CCNB2, (C) PRC1, (D) RACGAP1, (E) AURKA, (F) CDKN3, (G) NUSAP1, (H) CDCA5, (I) ASPM, and (J) NCAPG. P value < 0.05 was considered statistically significant.
Figure 10. Overall survival analysis for combining the expression of single hub gene and macrophage in HCC patients after adjusting five confounding factors, including age, stage, gender, race, and tumor purity. (A) TOP2A, (B) CCNB2, (C) PRC1, (D) RACGAP1, (E) AURKA, (F) CDKN3, (G) NUSAP1, (H) CDCA5, (I) ASPM, and (J) NCAPG. Log-rank p value < 0.05 was considered statistically significant. HCC, hepatocellular carcinoma; HR, hazard ratio.
Discussion
Hepatocellular carcinoma, one of the malignant cancers with high heterogeneity, is frequently diagnosed in its middle and advanced stages (Friemel et al., 2015; Buczak et al., 2018). Surgical resection remains the most crucial technique for HCC treatment; however, its therapeutic effects are always unsatisfactory (Abe et al., 2017; Deng et al., 2020). Thus, there is a need to screen novel carcinoma biomarkers and therapeutic targets. In the present study, comprehensive and bioinformatics analyses of multi-databases were applied to determine ten hub genes that appeared to be correlated with infiltrated immune cells in HCC. These genes were identified as independent prognostic factors in HCC patients.
In our study, five GEO datasets (GSE46408, GSE57957, GSE74656, GSE76427, and GSE87630) shared 107 common DEGs in HCC tissues; following this, the PPI network was constructed based on these genes. The results revealed a most significant module using the MCODE analysis that completely coincided with the subnetwork identification from the CytoHubba analysis. Particularly, ten hub genes in the module were upregulated in HCC tissues: TOP2A, CCNB2, PRC1, RACGAP1, AURKA, CDKN3, NUSAP1, CDCA5, ASPM, and NCAPG, respectively. The enrichment analyses presented that these hub genes were dramatically enriched in several terms of the BP analysis, including cell division, the positive regulation of mitotic nuclear division, female gamete generation, and the mitotic cell cycle phase transition. These genes were also significantly enriched in terms of the mitotic spindle in the CC analysis and in terms of the protein kinase binding in the MP analysis, suggesting that there is a close association between the hub genes and tumorigenesis. The validation in Oncomine and GEPIA2 confirms that the mRNA levels of TOP2A, CCNB2, PRC1, RACGAP1, AURKA, CDKN3, NUSAP1, CDCA5, ASPM, and NCAPG were significantly overexpressed in the HCC tissues, and at this validation, the p < 0.0001 and p < 0.01 were set in Oncomine and GEPIA2, respectively to more accurately identify the expression pattern of hub genes between HCC and normal tissues. Further HPA analysis also demonstrated that, compared to normal liver tissues, HCC tissues had significantly higher protein levels of TOP2A, PRC1, RACGAP1, AURKA, NUSAP1, and CDCA5, while we could not obtain the protein expression of CCNB2, CDKN3, ASPM, and NCAPG of HCC from the HPA website. These hub genes were validated to be closely correlated with infiltrated immune cells using the TIMER2.0 database. Both survival curves and univariate Cox regression analyses suggested that these candidate hub genes have a strong predictive ability for HCC. Previous studies demonstrated that TAMs extremely facilitate tumor angiogenesis and lead to a detrimental prognosis in HCC (Zhang et al., 2019; Zhao et al., 2020). The identification of TAM-related genes will facilitate providing more potential targets of the individualized precision treatment and improve the prognosis of HCC. Thus, we explored the prognostic value for the combination of candidate gene and TAM expression patterns in HCC and found that there was no significant relationship between TAMs and prognosis under a low expression level of the single hub gene. Meanwhile, under the high expression of CCNB2/RACGAP1/AURKA/CDKN3/ASPM/NCAPG, high TAM levels predicted unfavorable prognosis. Furthermore, the multivariate Cox regression models indicated that all candidate hub genes were independent predictors and that combining their respective expression levels with TAM will help them play a more effective role in the prognosis prediction of HCC.
In the enrichment analysis of candidate hub genes, we observed that nine out of ten genes were significantly associated with cell division: TOP2A, CCNB2, PRC1, RACGAP1, AURKA, NUSAP1, CDCA5, ASPM, and NCAPG. The remaining gene, CDKN3, has a dual function in the regulation of the cell cycle. On one hand, CDKN3 serves as a cyclin-dependent kinase inhibitor while interacting with and dephosphorylating CDK2 kinase, thereby restraining its activation (Hannon et al., 1994; Johnson et al., 2002); on the other hand, CDKN3 can act as a MDM2-binding protein that forms a complex with MDM2 and P53, thus suppressing the production of P21, leading to the acceleration of cell cycle progression (Okamoto et al., 2006). These candidate hub genes have been demonstrated to work as oncogenes and are associated with clinical prognosis in numerous solid neoplasms, particularly in HCC (Roy et al., 2018; Gong et al., 2019; Lin J. et al., 2019). The TOP2A gene encoded a DNA topoisomerase, which controls and changes the topological status of DNA in the process of transcription and functions as a target for some antitumor agents (Delgado et al., 2018; Kitdumrongthum et al., 2020). Other bioinformatic analyses showed that TOP2A was related to development in cancers of the liver, esophagus, stomach, cervix, and lung, among others (Zeng et al., 2019; Kou et al., 2020; Zhang T. et al., 2020; Zhou et al., 2020). CCNB2 is an important element for the process of cell cycle regulation. U2AF homology motif kinase 1 facilitates the nuclear enrichment of MYB proto-oncogene like 2 by affecting the expression of CCNB2 to regulate cell cycle and proliferation (Wei et al., 2019), and reduced transmembrane protein 9 can contribute to decreased CCNB2 levels and then promote cell cycle arrest in HCC cells (Zhang et al., 2016). PRC1, RACGAP1, and CDCA5 were identified as the crucial genes in the pathological progression from cirrhosis to HCC, and their hypomethylation may drive the high expression of these genes (Lin Y. et al., 2019). AURKA can induce the metastasis of irradiated residuary HCC while promoting an epithelial-mesenchymal transition and cancer stem cell properties (Chen et al., 2017). MYC proto-oncogene and AURKA regulate the expression of each other at a mRNA level identified as a MYC-AURKA feedback loop (Lu et al., 2015). CDKN3 overexpression can shorten the survival of HCC cells and shift sensitivity to chemotherapeutic drugs across the AKT/P53/P21 signaling pathway; besides, CDKN3 has been shown to be downregulated in advanced tumor stages (Dai et al., 2016). On the contrary, Chunyang et al. presented that the upregulation of CDKN3 might facilitate cell proliferation via the stimulation of the G1-S transition (Xing et al., 2012). NUSAP1 is a target for mir193A-5p, and evidence has shown that mir193A-5p might block the tumorigenesis of HCC through reducing NUSAP1 (Roy et al., 2018). ASPM is associated with the development of HCV-related cirrhosis via the regulation of tumor-associated phosphorylation (Wang et al., 2017). ASPM is also considered a prognostic biomarker that predicts the increased possibility of invasive or metastatic HCC (Lin et al., 2008). NCAPG down-regulation indicates the suppression of HCC progression, possibly via the PI3K-AKT signaling pathway (Gong et al., 2019; Wang Y. et al., 2019). However, among these candidate hub genes, only PRC1 has been reported to be an immune-related gene in a weighted gene co-expression network analysis (Wang et al., 2020).
Previous bioinformatics analyses have revealed that some TOP2A, CCNB2, PRC1, RACGAP1, AURKA, NUSAP1, CDCA5, ASPM, and NCAPG can be identified as key genes basing on different screening rules (Cai et al., 2019; Wang M. et al., 2019; Zhou Z. et al., 2019; Song et al., 2020). Comparing these previous studies, ours has the following advantages: First, this study included five GEO datasets, while others included two or three gene expression microarrays. In general, a greater number of included samples indicate more credible results in integrated research. Second, we constructed ten gene-macrophage Cox regression models. Inevitably, there were still several limitations of the present study. The included datasets came from different platforms, which might lead to an uncertain systematic bias. Furthermore, TIMER2.0 is a visual website based on tumor tissue information from the Cancer Genome Atlas database (Li et al., 2020). Thus, although tumor purity adjustment was performed in the correlation analyses between the immune cell and candidate genes, there was still systematic bias. To overcome this issue, the application of single-cell RNA sequencing at a higher resolution should be conducted (Papalexi and Satija, 2018). Finally, future experiments in vivo/in vitro should be performed to verify the results of this bioinformatics analysis.
Numerous studies have revealed that immune cell infiltration, TAMs, for instance, can serve as a biomarker for the diagnosis and prognosis of various cancers (Ali et al., 2016; Zhou R. et al., 2019). Thus, we assessed the prognostic value of the combination of TAMs and expression patterns for each of the hub genes. The results showed that the high TAM level predicted unfavorable prognosis under the condition of the high expression of the hub genes, while there was no significant correlation between the TAM and prognosis under the condition of the low expression level of these genes. These results suggested that the combination of hub-gene expression and the TAM levels would play a more effective role in the prognosis prediction of HCC.
In summary, we identified ten genes with a positive correlation with infiltrated-immune cells. These candidate genes present the marked prognostic value of HCC and act as independent prognosis factors for patients with HCC. Moreover, these genes may function in the progression of HCC. Furthermore, we determined that combining the expression of these genes and TAMs can provide a more efficient HCC prognosis prediction. Overall, these findings suggest that these hub genes may be potential targets of immune therapy.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author Contributions
XQ and LX: conceptualization. XL and LH: GEO searches and analyses. JW, ZH, and MC: transcriptome data visualization and validation. HC: writing. HC, JW, and LL: proofreading. All authors have read and approved the manuscript.
Funding
This research was supported by the National Natural Science Foundation of China (81960390) and the National Natural Science Foundation of Guangxi (2017GXNSFBA198068).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2020.575762/full#supplementary-material
Abbreviations
HCC, hepatocellular carcinoma; DEGs, differentially expressed genes; GEO, Gene Expression Omnibus; PPI, protein–protein interaction; GEPIA2, Gene Expression Profiling Interactive Analysis 2; HPA, human protein atlas; GO, Gene Ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes; STRING, Search Tool for the Retrieval of Interacting Genes; BP, biological process; CC, cellular component; MF, molecular function; TOP2A, DNA topoisomerase II alpha; CCNB2, cyclin B2; PRC1, protein regulator of cytokinesis 1; RACGAP1, Rac GTPase-activating protein 1; AURKA, aurora kinase A; CDKN3, cyclin-dependent kinase inhibitor 3; NUSAP1, nucleolar and spindle-associated protein 1; CDCA5, cell division cycle-associated 5; ASPM, abnormal spindle microtubule assembly; NCAPG, non-SMC condensin I complex subunit G; DC, dendritic cell; COR, correlation coefficient; OS, overall survival; RFS, relapse-free survival; TAM, tumor-associated macrophage; HR, hazard rate; FC, fold change; MCODE, Molecular Complex Detection; MCC, Maximal Clique Centrality.
Footnotes
- ^ http://www.ncbi.nlm.nih.gov/geo
- ^ https://www.ncbi.nlm.nih.gov/geo/geo2r/
- ^ https://www.oncomine.org/
- ^ http://gepia2.cancer-pku.cn
- ^ https://www.proteinatlas.org/
- ^ http://timer.cistrome.org/
References
Abe, T., Tashiro, H., Kobayashi, T., Hattori, M., Kuroda, S., and Ohdan, H. (2017). Glasgow prognostic score and prognosis after hepatectomy for hepatocellular carcinoma. World J. Surg. 41, 1860–1870. doi: 10.1007/s00268-017-3909-7
Ali, H. R., Chlon, L., Pharoah, P. D., Markowetz, F., and Caldas, C. (2016). Patterns of immune infiltration in breast cancer and their clinical implications: a gene-expression-based retrospective Study. PLoS Med. 13:e1002194. doi: 10.1371/journal.pmed.1002194
Buczak, K., Ori, A., Kirkpatrick, J. M., Holzer, K., Dauch, D., Roessler, S., et al. (2018). Spatial tissue proteomics quantifies inter- and intratumor heterogeneity in hepatocellular carcinoma (HCC). Mol. Cell. Proteom. MCP 17, 810–825. doi: 10.1074/mcp.RA117.000189
Cai, C., Wang, W., and Tu, Z. (2019). Aberrantly DNA methylated-differentially expressed genes and pathways in hepatocellular carcinoma. J. Cancer 10, 355–366. doi: 10.7150/jca.27832
Caruso, S., and Nault, J. C. (2019). A dive into the deep heterogeneity of hepatocellular carcinoma. Gastroenterology 157, 1477–1479. doi: 10.1053/j.gastro.2019.10.003
Chen, C., Chen, H., Zhang, Y., Thomas, H. R., Frank, M. H., He, Y., et al. (2020). TBtools - an integrative toolkit developed for interactive analyses of big biological data. Mol. Plant [Epub ahead of print]. doi: 10.1016/j.molp.2020.06.009
Chen, C., Song, G., Xiang, J., Zhang, H., Zhao, S., and Zhan, Y. (2017). AURKA promotes cancer metastasis by regulating epithelial-mesenchymal transition and cancer stem cell properties in hepatocellular carcinoma. Biochem. Biophys. Res. Commun. 486, 514–520. doi: 10.1016/j.bbrc.2017.03.075
Dai, W., Miao, H., Fang, S., Fang, T., Chen, N., and Li, M. (2016). CDKN3 expression is negatively associated with pathological tumor stage and CDKN3 inhibition promotes cell survival in hepatocellular carcinoma. Mol. Med. Rep. 14, 1509–1514. doi: 10.3892/mmr.2016.5410
Davis, S., and Meltzer, P. S. (2007). GEOquery: a bridge between the gene expression omnibus (GEO) and BioConductor. Bioinformatics 23, 1846–1847. doi: 10.1093/bioinformatics/btm254
Delgado, J. L., Hsieh, C. M., Chan, N. L., and Hiasa, H. (2018). Topoisomerases as anticancer targets. Biochem. J. 475, 373–398. doi: 10.1042/BCJ20160583
Deng, Y., Zhu, J., Liu, Z., Huang, M., Chang, D. W., and Gu, J. (2020). Elevated systemic inflammatory responses, factors associated with physical and mental quality of life, and prognosis of hepatocellular carcinoma. Aging 12, 4357–4370. doi: 10.18632/aging.102889
Duan, M., Goswami, S., Shi, J. Y., Wu, L. J., Wang, X. Y., Ma, J. Q., et al. (2019). Activated and exhausted MAIT cells foster disease progression and indicate poor outcome in hepatocellular carcinoma. Clin. Cancer Res. 25, 3304–3316. doi: 10.1158/1078-0432.CCR-18-3040
Forner, A., Reig, M., and Bruix, J. (2018). Hepatocellular carcinoma. Lancet 391, 1301–1314. doi: 10.1016/S0140-6736(18)30010-2
Fox, R., Berhane, S., Teng, M., Cox, T., Tada, T., Toyoda, H., et al. (2014). Biomarker-based prognosis in hepatocellular carcinoma: validation and extension of the BALAD model. Br. J. Cancer 110, 2090–2098. doi: 10.1038/bjc.2014.130
Friemel, J., Frick, L., and Weber, A. (2015). Intratumor heterogeneity in HCC. Aging 7, 350–351. doi: 10.18632/aging.100760
Gong, C., Ai, J., Fan, Y., Gao, J., Liu, W., Feng, Q., et al. (2019). NCAPG promotes the proliferation of hepatocellular carcinoma through PI3K/AKT signaling. OncoTargets Ther. 12, 8537–8552. doi: 10.2147/OTT.S217916
Hannon, G. J., Casso, D., and Beach, D. (1994). KAP: a dual specificity phosphatase that interacts with cyclin-dependent kinases. Proc. Natl. Acad. Sci. U.S.A. 91, 1731–1735. doi: 10.1073/pnas.91.5.1731
Huang, R., Chen, Z., Li, W., Fan, C., and Liu, J. (2020). Immune systemassociated genes increase malignant progression and can be used to predict clinical outcome in patients with hepatocellular carcinoma. Int. J. Oncol. 56, 1199–1211. doi: 10.3892/ijo.2020.4998
Imamura, H., Matsuyama, Y., Tanaka, E., Ohkubo, T., Hasegawa, K., Miyagawa, S., et al. (2003). Risk factors contributing to early and late phase intrahepatic recurrence of hepatocellular carcinoma after hepatectomy. J. Hepatol. 38, 200–207. doi: 10.1016/s0168-8278(02)00360-4
Jayant, K., Habib, N., Huang, K. W., Podda, M., Warwick, J., and Arasaradnam, R. (2020). Immunological basis of genesis of hepatocellular carcinoma: unique challenges and potential opportunities through immunomodulation. Vaccines 8:247. doi: 10.3390/vaccines8020247
Johnson, L. N., De Moliner, E., Brown, N. R., Song, H., Barford, D., Endicott, J. A., et al. (2002). Structural studies with inhibitors of the cell cycle regulatory kinase cyclin-dependent protein kinase 2. Pharmacol. Ther. 93, 113–124. doi: 10.1016/s0163-7258(02)00181-x
Kim, S., Rovgaliyev, B., Lee, J. M., Lee, K. W., Hong, S. K., Cho, J. H., et al. (2020). Clinical significance of de novo malignancy after liver transplant: a single-center study. Transpl. Proc. [Epub ahead of print]. doi: 10.1016/j.transproceed.2020.02.148
Kitdumrongthum, S., Reabroi, S., Suksen, K., Tuchinda, P., Munyoo, B., Mahalapbutr, P., et al. (2020). Inhibition of topoisomerase IIalpha and induction of DNA damage in cholangiocarcinoma cells by altholactone and its halogenated benzoate derivatives. Biomed. Pharmacother. Biomed. Pharmacother. 127:110149. doi: 10.1016/j.biopha.2020.110149
Kou, F., Sun, H., Wu, L., Li, B., Zhang, B., Wang, X., et al. (2020). TOP2A promotes lung adenocarcinoma cells’ malignant progression and predicts poor prognosis in lung adenocarcinoma. J. Cancer 11, 2496–2508. doi: 10.7150/jca.41415
Langhans, B., Nischalke, H. D., Kramer, B., Dold, L., Lutz, P., Mohr, R., et al. (2019). Role of regulatory T cells and checkpoint inhibition in hepatocellular carcinoma. Cancer Immunol. Immunother. CII 68, 2055–2066. doi: 10.1007/s00262-019-02427-4
Li, T., Fu, J., Zeng, Z., Cohen, D., Li, J., Chen, Q., et al. (2020). TIMER2.0 for analysis of tumor-infiltrating immune cells. Nucleic Acids Res. 48, W509–W514. doi: 10.1093/nar/gkaa407
Li, Z., Wu, T., Zheng, B., and Chen, L. (2019). Individualized precision treatment: targeting TAM in HCC. Cancer Lett. 458, 86–91. doi: 10.1016/j.canlet.2019.05.019
Lin, J., Hou, Y., Huang, S., Wang, Z., Sun, C., Wang, Z., et al. (2019). Exportin-T promotes tumor proliferation and invasion in hepatocellular carcinoma. Mol. Carcinogenesis 58, 293–304. doi: 10.1002/mc.22928
Lin, S. Y., Pan, H. W., Liu, S. H., Jeng, Y. M., Hu, F. C., Peng, S. Y., et al. (2008). ASPM is a novel marker for vascular invasion, early recurrence, and poor prognosis of hepatocellular carcinoma. Clin. Cancer Res. 14, 4814–4820. doi: 10.1158/1078-0432.CCR-07-5262
Lin, Y., Liang, R., Ye, J., Li, Q., Liu, Z., Gao, X., et al. (2019). A twenty gene-based gene set variation score reflects the pathological progression from cirrhosis to hepatocellular carcinoma. Aging 11, 11157–11169. doi: 10.18632/aging.102518
Lu, L., Han, H., Tian, Y., Li, W., Zhang, J., Feng, M., et al. (2015). Aurora kinase A mediates c-Myc’s oncogenic effects in hepatocellular carcinoma. Mol. Carcinogenesis 54, 1467–1479. doi: 10.1002/mc.22223
Ma, X., Zhou, L., and Zheng, S. (2020). Transcriptome analysis revealed key prognostic genes and microRNAs in hepatocellular carcinoma. PeerJ 8:e8930. doi: 10.7717/peerj.8930
Okamoto, K., Kitabayashi, I., and Taya, Y. (2006). KAP1 dictates p53 response induced by chemotherapeutic agents via Mdm2 interaction. Biochem. Biophys. Res. Commun. 351, 216–222. doi: 10.1016/j.bbrc.2006.10.022
Papalexi, E., and Satija, R. (2018). Single-cell RNA sequencing to explore immune cell heterogeneity. Nat. Rev. Immunol. 18, 35–45. doi: 10.1038/nri.2017.76
Park, H., and Park, J. Y. (2013). Clinical significance of AFP and PIVKA-II responses for monitoring treatment outcomes and predicting prognosis in patients with hepatocellular carcinoma. BioMed Res. Int. 2013:310427. doi: 10.1155/2013/310427
Ponten, F., Jirstrom, K., and Uhlen, M. (2008). The Human Protein Atlas–a tool for pathology. J. Pathol. 216, 387–393. doi: 10.1002/path.2440
Rhodes, D. R., Kalyana-Sundaram, S., Mahavisno, V., Varambally, R., Yu, J., Briggs, B. B., et al. (2007). Oncomine 3.0: genes, pathways, and networks in a collection of 18,000 cancer gene expression profiles. Neoplasia 9, 166–180. doi: 10.1593/neo.07112
Roy, S., Hooiveld, G. J., Seehawer, M., Caruso, S., Heinzmann, F., Schneider, A. T., et al. (2018). microRNA 193a-5p regulates levels of nucleolar- and spindle-associated protein 1 to suppress hepatocarcinogenesis. Gastroenterology 155, 1951.e26–1966.e26. doi: 10.1053/j.gastro.2018.08.032
Siegel, R. L., Miller, K. D., and Jemal, A. (2015). Cancer statistics, 2015. CA Cancer J. Clin. 65, 5–29. doi: 10.3322/caac.21254
Song, X., Du, R., Gui, H., Zhou, M., Zhong, W., Mao, C., et al. (2020). Identification of potential hub genes related to the progression and prognosis of hepatocellular carcinoma through integrated bioinformatics analysis. Oncol. Rep. 43, 133–146. doi: 10.3892/or.2019.7400
Szklarczyk, D., Gable, A. L., Lyon, D., Junge, A., Wyder, S., Huerta-Cepas, J., et al. (2019). STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 47, D607–D613. doi: 10.1093/nar/gky1131
Tang, Y., Zeng, Z., Wang, J., Li, G., Huang, C., Dong, X., et al. (2020). Combined signature of nine immune-related genes: a novel risk score for predicting prognosis in hepatocellular carcinoma. Am. J. Transl. Res. 12, 1184–1202.
Tang, Z., Kang, B., Li, C., Chen, T., and Zhang, Z. (2019). GEPIA2: an enhanced web server for large-scale expression profiling and interactive analysis. Nucleic Acids Res. 47, W556–W560. doi: 10.1093/nar/gkz430
Villanueva, A. (2019). Hepatocellular carcinoma. New Engl. J. Med. 380, 1450–1462. doi: 10.1056/NEJMra1713263
Wang, D., Liu, J., Liu, S., and Li, W. (2020). Identification of crucial genes associated with immune cell infiltration in hepatocellular carcinoma by weighted gene co-expression network analysis. Front. Genet. 11:342. doi: 10.3389/fgene.2020.00342
Wang, F., Chang, Y., Li, J., Wang, H., Zhou, R., Qi, J., et al. (2017). Strong correlation between ASPM gene expression and HCV cirrhosis progression identified by co-expression analysis. Digestive Liver Dis. 49, 70–76. doi: 10.1016/j.dld.2016.10.017
Wang, J., Meng, J., Ran, W., Lee, R. J., Teng, L., Zhang, P., et al. (2019). Hepatocellular carcinoma growth retardation and PD-1 blockade therapy potentiation with synthetic high-density lipoprotein. Nano Lett. 19, 5266–5276. doi: 10.1021/acs.nanolett.9b01717
Wang, M., Wang, L., Wu, S., Zhou, D., and Wang, X. (2019). Identification of key genes and prognostic value analysis in hepatocellular carcinoma by integrated bioinformatics analysis. Int. J. Genomics 2019:3518378. doi: 10.1155/2019/3518378
Wang, Y., Gao, B., Tan, P. Y., Handoko, Y. A., Sekar, K., Deivasigamani, A., et al. (2019). Genome-wide CRISPR knockout screens identify NCAPG as an essential oncogene for hepatocellular carcinoma tumor growth. FASEB J. 33, 8759–8770. doi: 10.1096/fj.201802213RR
Wei, T., Weiler, S. M. E., Toth, M., Sticht, C., Lutz, T., Thomann, S., et al. (2019). YAP-dependent induction of UHMK1 supports nuclear enrichment of the oncogene MYBL2 and proliferation in liver cancer cells. Oncogene 38, 5541–5550. doi: 10.1038/s41388-019-0801-y
Xie, W., Wang, B., Wang, X., Hou, D., Su, H., and Huang, H. (2020). Nine hub genes related to the prognosis of HBV-positive hepatocellular carcinoma identified by protein interaction analysis. Ann. Transl. Med. 8:478. doi: 10.21037/atm.2020.03.94
Xing, C., Xie, H., Zhou, L., Zhou, W., Zhang, W., Ding, S., et al. (2012). Cyclin-dependent kinase inhibitor 3 is overexpressed in hepatocellular carcinoma and promotes tumor cell proliferation. Biochem. Biophys. Res. Commun. 420, 29–35. doi: 10.1016/j.bbrc.2012.02.107
Yin, Z., Dong, C., Jiang, K., Xu, Z., Li, R., Guo, K., et al. (2019). Heterogeneity of cancer-associated fibroblasts and roles in the progression, prognosis, and therapy of hepatocellular carcinoma. J. Hematol. Oncol. 12:101. doi: 10.1186/s13045-019-0782-x
Zeng, L., Fan, X., Wang, X., Deng, H., Zhang, K., Zhang, X., et al. (2019). Bioinformatics analysis based on multiple databases identifies hub genes associated with hepatocellular carcinoma. Curr. Genomics 20, 349–361. doi: 10.2174/1389202920666191011092410
Zhang, J., Chang, L., Zhang, X., Zhou, Z., and Gao, Y. (2019). Meta-analysis of the prognostic and clinical value of tumor-associated macrophages in hepatocellular carcinoma. J. Investig. Surg. [Epub ahead of print]. doi: 10.1080/08941939.2019.1631411
Zhang, T., Wang, B. F., Wang, X. Y., Xiang, L., Zheng, P., Li, H. Y., et al. (2020). Key genes associated with prognosis and tumor infiltrating immune cells in gastric cancer patients identified by cross-database analysis. Cancer Biother. Radiopharm. 35, 696–710. doi: 10.1089/cbr.2019.3423
Zhang, Y., Ran, Y., Xiong, Y., Zhong, Z. B., Wang, Z. H., Fan, X. L., et al. (2016). Effects of TMEM9 gene on cell progression in hepatocellular carcinoma by RNA interference. Oncol. Rep. 36, 299–305. doi: 10.3892/or.2016.4821
Zhang, Y., Zhang, L., Xu, Y., Wu, X., Zhou, Y., and Mo, J. (2020). Immune-related long noncoding RNA signature for predicting survival and immune checkpoint blockade in hepatocellular carcinoma. J. Cell. Physiol. 235, 9304–9316. doi: 10.1002/jcp.29730
Zhao, X., Wang, X., You, Y., Wen, D., Feng, Z., Zhou, Y., et al. (2020). Nogo-B fosters HCC progression by enhancing Yap/Taz-mediated tumor-associated macrophages M2 polarization. Exp. Cell Res. 391:111979. doi: 10.1016/j.yexcr.2020.111979
Zhou, R., Zhang, J., Zeng, D., Sun, H., Rong, X., Shi, M., et al. (2019). Immune cell infiltration as a biomarker for the diagnosis and prognosis of stage I-III colon cancer. Cancer Immunol. Immunother. CII 68, 433–442. doi: 10.1007/s00262-018-2289-7
Zhou, W., Wu, J., Liu, X., Ni, M., Meng, Z., Liu, S., et al. (2020). Identification of crucial genes correlated with esophageal cancer by integrated high-throughput data analysis. Medicine 99:e20340. doi: 10.1097/MD.0000000000020340
Zhou, Y., Zhou, B., Pache, L., Chang, M., Khodabakhshi, A. H., Tanaseichuk, O., et al. (2019). Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat. Commun. 10:1523. doi: 10.1038/s41467-019-09234-6
Keywords: hepatocellular carcinoma, immune infiltration, prognosis, tumor-associated macrophage, biomarker
Citation: Chen H, Wu J, Lu L, Hu Z, Li X, Huang L, Zhang X, Chen M, Qin X and Xie L (2021) Identification of Hub Genes Associated With Immune Infiltration and Predict Prognosis in Hepatocellular Carcinoma via Bioinformatics Approaches. Front. Genet. 11:575762. doi: 10.3389/fgene.2020.575762
Received: 24 June 2020; Accepted: 07 December 2020;
Published: 11 January 2021.
Edited by:
Yanni Sun, City University of Hong Kong, Hong KongReviewed by:
Manal Said Fawzy, Suez Canal University, EgyptKa-Chun Wong, City University of Hong Kong, Hong Kong
Copyright © 2021 Chen, Wu, Lu, Hu, Li, Huang, Zhang, Chen, Qin and Xie. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xue Qin, qinxue919@126.com; Li Xie, drlixie@163.com