Identification of Rad51 as a prognostic biomarker correlated with immune infiltration in hepatocellular carcinoma

The DNA Damage Response and Inflammation in Cancer

Genomic stability in normal cells is crucial to avoid oncogenesis. Accordingly, multiple components of the DNA damage response (DDR) operate as bona fide tumor suppressor proteins by preserving genomic stability, eliciting the demise of cells with unrepairable DNA lesions, and engaging cell-extrinsic oncosuppression via immunosurveillance. That said, DDR sig-naling can also favor tumor progression and resistance to therapy. Indeed, DDR signaling in cancer cells has been consistently linked to the inhibition of tumor-targeting immune responses. Here, we discuss the complex interactions between the DDR and inflammation in the context of oncogenesis, tumor progression, and response to therapy.

SIGNIFICANCE

Accumulating preclinical and clinical evidence indicates that DDR is intimately connected to the emission of immunomodulatory signals by normal and malignant cells, as part of a cell-extrinsic program to preserve organismal homeostasis. DDR-driven inflammation, however, can have diametrically opposed effects on tumor-targeting immunity. Understanding the links between the DDR and inflammation in normal and malignant cells may unlock novel immunotherapeutic paradigms to treat cancer.

RAD51 Inhibition Shows Antitumor Activity in Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC), the major type of liver cancer, causes a high annual mortality worldwide. RAD51 is the critical recombinase responsible for homologous recombination (HR) repair in DNA damage. In this study, we identified that RAD51 was upregulated in HCC and that RAD51 silencing or inhibition reduced the proliferation, migration, and invasion of HCC cells and enhanced cell apoptosis and DNA damage. HCC cells with the combinatorial treatments of RAD51 siRNA or inhibitor and sorafenib demonstrated a synergistic effect in inhibiting HCC cell proliferation, migration, and invasion, as well as inducing cell apoptosis and DNA damage. Single RAD51 silencing or sorafenib reduced RAD51 protein expression and weakened HR efficiency, and their combination almost eliminated RAD51 protein expression and inhibited HR efficiency further. An in vivo tumor model confirmed the RAD51 inhibitor's antitumor activity and synergistic antitumor activity with sorafenib in HCC. RNA-Seq and gene set enrichment analysis (GSEA) in RAD51-inactivated Huh7 cells indicated that RAD51 knockdown upregulated cell apoptosis and G1/S DNA damage checkpoint pathways while downregulating mitotic spindle and homologous recombination pathways. Our findings suggest that RAD51 inhibition exhibits antitumor activities in HCC and synergizes with sorafenib. Targeting RAD51 may provide a novel therapeutic approach in HCC.

RAD54L promotes progression of hepatocellular carcinoma via the homologous recombination repair pathway

Hepatocellular carcinoma (HCC) is a malignant tumor with high incidence worldwide. The underlying mechanisms remain poorly understood. The DNA metabolic process of homologous recombination repair (HRR) has been linked to a high probability of tumorigenesis and drug resistance. This study aimed to determine the role of HRR in HCC and identify critical HRR-related genes that affect tumorigenesis and prognosis. A total of 613 tumor and 252 para-carcinoma tissue samples were collected from The Cancer Genome Atlas (TCGA) and the International Cancer Genome Consortium (ICGC) to obtain differentially expressed genes (DEGs). HRR-related genes were assessed using gene enrichment and pathway analyses. Survival analysis was performed using the Kaplan-Meier method in the Gene Expression Profiling Interactive Analysis portal. The levels of RAD54L in the HRR pathway were detected by RT-qPCR and western blotting in para-carcinoma and HCC tissues and in L02 normal human liver cells and Huh7 HCC cells. Immunohistochemistry (IHC) was performed on the clinical specimens to determine the connection between gene expression and clinical features. Bioinformatics analysis revealed that the HRR pathway was enriched in HCC tissues. Upregulation of HRR pathway DEGs in HCC tissues was positively correlated with tumor pathological staging and negatively associated with patient overall survival. RAD54B, RAD54L, and EME1 genes in the HRR pathway were screened as markers for predicting HCC prognosis. RT-qPCR identified RAD54L as the most significantly expressed of the three genes. Western blotting and IHC quantitative analyses further demonstrated that RAD54L protein levels were higher in HCC tissues. IHC analysis of 39 pairs of HCC and para-carcinoma tissue samples also revealed an association between RAD54L and Edmondson-Steiner grade and the proliferation-related gene Ki67. The collective findings positively correlate RAD54L in the HRR signaling pathway with HCC staging and implicate RAD54L as a marker to predict HCC progression.

Pan-cancer analysis of the prognostic and immunological role of Fanconi anemia complementation group E

Fanconi anemia (FA) genes contribute to tumorigenesis by regulating DNA repair. Despite its importance for assembly and functionality of the FA core complex, no pan-cancer analysis of FANCE was performed. We aimed to provide a comprehensive understanding of the role of FANCE in cancers. Based on The Cancer Genome Atlas (TCGA), Cancer Cell Line Encyclopedia (CCLE), Genotype Tissue-Expression (GTEx), Human Protein Atlas (HPA), Gene Expression Omnibus (GEO), and Cancer Single-cell Atlas (CancerSEA) databases, we investigated the carcinogenicity of FANCE using various bioinformatics methods, including FANCE expression and prognosis, immune invasion, tumor mutation burden, microsatellite instability, and neoantigens. We monitored Fance mutations in mice that caused tumorigenesis. FANCE expression and activity scores were upregulated in 15 and 21 cancers. High expression of FANCE affected shorter overall survival (OS) in seven cancers and longer overall survival in three cancers. It was correlated with shorter overall survival and progression-free interval (PFI) in endometrial cancer and longer overall survival and PFI in cervical cancer. FANCE expression negatively correlated with stromal/immune scores in 21 cancers including cervical cancer, endometrial cancer, and ovarian cancer. FANCE expression negatively correlated with CD8 T cells in endometrial cancer and positively correlated with M1 macrophages in cervical cancer, possibly related to cancer prognosis. FANCE positively correlated with immune checkpoint inhibitors PD-1, PD-L1, and CTLA4 in endometrial cancer and ovarian cancer. FANCE expression positively correlated with microsatellite instability, tumor mutational burden, and neoantigens in 7, 22, and five cancers, especially in endometrial cancer, potentially increasing the effectiveness of immunotherapy. Single-cell sequencing data showed FANCE was primarily expressed in cancer cells in cervical and ovarian cancer, and in fibroblasts in endometrial cancer. Fance heterozygous mutant mice had increased tumor incidences and shorter overall survival and tumor-free survival (TFS) than Fance homozygous mutant mice and wild-type mice. Conclusively, FANCE potential to serve as a biomarker for cancer prognosis and may predict cancer immunotherapy responses. Fance heterozygous mutant resulted in increased tumorigenesis and poor prognosis in mice.

Prognostic biomarkers correlated with immune infiltration in non-small cell lung cancer

Lung cancer is the leading cause of cancer-related mortality in men and women globally. Non-small cell lung cancer (NSCLC) is the most prevalent subtype, accounting for 85-90% of all cancers. Although there have been dramatic advances in therapeutic approaches in recent decades, the recurrence and metastasis rates of NSCLC are as high as 30-40% with the 5-year overall survival rate being less than 15%. Therefore, it is necessary to explore the pathogenesis of NSCLC at the genetic level and identify prognostic biomarkers and novel therapeutic targets. Here, we aimed to identify mutated genes with high frequencies in Chinese NSCLC patients using next-generation sequencing and to investigate their relationships with the tumor mutation burden (TMB) and tumor immune microenvironment. A total of 110 NSCLC patients were enrolled to profile the genetic variations. Mutations in EGFR (62.37%), TP53 (61.29%), LRP1B (13.98%), FAT1 (12.90%), KMT2D (11.83%), CREBBP (10.75%), and RB1 (9.68%) were most prevalent. TP53, LRP1B, KMT2D, and CREBBP mutations were all significantly associated with high TMB (P < 0.05 or P < 0.01). The infiltrating levels of immune cells and immune molecules were enriched significantly in the LRP1B mutation group. LRP1B mutations significantly correlated with stimulating and inhibitory immunoregulators. Gene set enrichment analysis revealed that cell cycle, the Notch signaling pathway, the insulin signaling pathway, and the mTOR signaling pathway are related to LRP1B mutations in the immune system. LRP1B mutations may be of clinical importance in enhancing the anti-tumor immune response and may be a promising biomarker for predicting immunotherapy responsiveness.

Screening for Biomarkers for Progression from Oral Leukoplakia to Oral Squamous Cell Carcinoma and Evaluation of Diagnostic Efficacy by Multiple Machine Learning Algorithms

The aim of the study is to identify key genes during the progression from oral leukoplakia (OL) to oral squamous cell carcinoma (OSCC) and predict effective diagnoses. Weighted gene co-expression network analysis (WGCNA) and differential expression analysis were performed to identify seven genes associated with the progression from OL to OSCC. Twelve machine learning algorithms including k-nearest neighbor (KNN), neural network (NNet), and extreme gradient boosting (XGBoost) were used to construct multi-gene models, which revealed that each model had good diagnostic efficacy. The functional mechanism or the pathways associated with these genes were evaluated using enrichment analysis, subtype clustering, and immune infiltration analysis. The enrichment analysis revealed that the genes enriched were associated with the cell cycle, cell division, and intracellular energy metabolism. The immunoassay results revealed that the genes primarily affected the infiltration of proliferating T cells and macrophage polarization. Finally, a nomogram and Kaplan-Meier survival analysis were used to predict the prognostic efficacy of key genes in OSCC patients. The results showed that genes could predict the prognosis of the patients, and patients in the high-risk group had a poor prognosis. Our study identified that the seven key genes, including DHX9, BCL2L12, RAD51, MELK, CDC6, ANLN, and KIF4A, were associated with the progression from OL to OSCC. These genes had good diagnostic efficacy and could be used as potential biomarkers for the prognosis of OSCC patients.

Targeting DNA repair to enhance the efficacy of sorafenib in hepatocellular carcinoma

The multityrosine kinase inhibitor sorafenib remains an important systemic treatment option for hepatocellular carcinoma (HCC). Signaling pathways, which are targeted by sorafenib, are involved in checkpoint and DNA repair response, RAD51 being a candidate protein. Here, we aim to evaluate the effect of the human RAD51 inhibitor B02 in combination with sorafenib in human HCC cells. Impact of RAD51 expression on HCC patient survival was evaluated by an in silico approach using Human Protein Atlas dataset. Cell viability of HUH7, AKH12, AKH13, and 3P was assessed by neutral red assay. To measure the cytotoxicity, we quantified loss of membrane integrity by lactate dehydrogenase release. We also employed colony formation assay and hanging drop method to assess clonogenic and invasive ability of HCC cell lines upon sorafenib and B02 treatment. Cell cycle distribution and characterization of apoptosis was evaluated by flow cytometry. In silico approach revealed that HCC patients with higher expression of RAD51 messenger RNA had a significantly shorter overall survival. The RAD51 inhibitor B02 alone and in combination with sorafenib significantly reduced viability, colony formation ability, and invasion capacity of HCC cells. Cell cycle analysis revealed that the combination of both agents reduces the proportion of cells in the G2/M phase while leading to an accumulating in the subG1 phase. The RAD51 inhibitor B02 seems to be a promising agent for HCC treatment and enhances the antitumor effects of sorafenib in vitro.

A Pan-Cancer Bioinformatic Analysis of RAD51 Regarding the Values for Diagnosis, Prognosis, and Therapeutic Prediction

Background

RAD51, a critical protein for DNA repairment, has been found to associate with multiple cancer types, but, so far, a systematic pan-cancer analysis of RAD51 has not been done yet.

Methods

Data were obtained from multiple open databases and genetic alteration, gene expression, survival association, functional enrichment, stemness, mutation association, immunity association, and drug therapy association of RAD51were analyzed. A prognostic model of RAD51 for overall glioma was constructed as an example application of RAD51 as a biomarker.

Results

RAD51 was overexpressed in 28 types of cancers and was associated with worse overall survival in 11 cancer types. RAD51 correlated genes were enriched in cell cycle terms. RAD51 was associated with cancer stemness, tumor mutational burden, and multiple immunomodulators in different cancer types. RAD51 expression was different across immune subtypes in 11 cancer types. RAD51 was closely associated with cancer immune microenvironments in some cancer types. Proliferating T cells was the cell type that expressed highest RAD51 across most of the cancer samples analyzed. RAD51 expression had an AUC of over 0.5 in 12 of the 23 ICB subcohorts. The Tumor Immune Dysfunction and Exclusion of 9 cancer types were different between RAD51 high and low groups. RAD51 expression showed negative correlations with the sensitivity of most drugs. A prognostic nomogram was constructed with a high confidence.

Conclusion

RAD51 is a clinical valuable biomarker for multiple cancer types, regarding its potential power for diagnosis, prognosis, and therapeutic prediction.

Ataxia-televangelist mutated (ATM)/ ATR serine/threonine kinase (ATR)-mediated RAD51 recombinase (RAD51) promotes osteogenic differentiation and inhibits osteoclastogenesis in osteoporosis

Osteoporosis is a metabolic bone disease that significantly affects the quality of life and can even lead to death. In this study, we aimed to investigate the role of RAD51 recombinase (RAD51) in osteoblast and osteoclast differentiation. We analyzed differentially expressed genes using microarray analysis. The osteogenic differentiation capability was analyzed by alkaline phosphatase (ALP) staining and alizarin red staining assays. Osteogenesis and osteoclast related genes expression was detected using quantitative real-time PCR (qPCR) and Western blotting. The phosphorylation of Ataxia-telangiectasia mutated (ATM) and ATR serine/threonine kinase (ATR) was tested using Western blotting. The effect of RAD51 on osteoporosis was also explored in vivo. The results showed that RAD51 was downregulated in osteoporosis, but upregulated in differentiated osteoblasts. Overexpression of RAD51 enhanced the differentiation of osteoblasts and suppressed the formation of osteoclasts. Furthermore, p-ATM and p-ATR levels were upregulated in osteoblasts and downregulated in osteoclasts. RAD51 expression was reduced by the ATM/ATR pathway inhibitor AZ20. AZ20 treatment inhibited osteoblastogenesis and promoted osteoclastogenesis, whereas RAD51 reversed the effects induced by AZ20. Moreover, RAD51 improved bone microarchitecture in vivo. Taken together, ATM/ATR signaling-mediated RAD51 promoted osteogenic differentiation and suppressed osteoclastogenesis. These findings reveal a critical role for RAD51 in osteoporosis.

Solute carrier family 12 member 8 (SLC12A8) is a potential biomarker and related to tumor immune cell infiltration in bladder cancer

The solute carrier family has been reported to play critical roles in the progression of several cancers; however, the relationship between solute carrier family 12 member 8 (SLC12A8) and bladder cancer (BC) has not been clearly confirmed. This study explores the prognostic value of SLC12A8 for BC and its correlation with immune cell infiltration. We found that the expression of SLC12A8 mRNA was significantly overexpressed in BC tissues compared with noncancerous tissues in multiple public databases, and the result was validated using real-time PCR and immunohistochemistry (IHC). The Kaplan-Meier method and Cox proportional hazards models were used to evaluate the prognostic value of SLC12A8 for BC. The high expression of SLC12A8 led to a shorter overall survival time and was an unfavorable prognostic biomarker for BC. The mechanisms of SLC12A8 promoting tumorigenesis were investigated by Gene Set Enrichment Analysis (GSEA). Moreover, the correlations of SLC12A8 expression with the tumor-infiltrating immune cells (TICs) in BC were explored using TIMER 2.0 and CIBERSORT. SLC12A8 was associated with CD4+ T cells, dendritic cells, neutrophils, and macrophages infiltration. The expression of SLC12A8 was positively correlated with crucial immune checkpoint molecules. In conclusion, SLC12A8 might be an unfavorable prognostic biomarker in BC related to tumor immune cell infiltration.

The long noncoding RNA noncoding RNA activated by DNA damage (NORAD)-microRNA-496-Interleukin-33 axis affects carcinoma-associated fibroblasts-mediated gastric cancer development

Carcinoma-associated fibroblasts (CAFs) are one of the crucial parts of in the tumor microenvironment and contribute to tumor progression. Interleukin-33 (IL-33), a tissue-derived nuclear cytokine from the IL-1 family, has been found abnormally expressed in tumor cells and Fibroblast. However, the role and mechanism of IL-33 in the interaction between gastric cancer (GC) cells and CAFs need investigation. Presently, we inquire into the function of lncRNA NORAD-miR-496 axis-mediated IL-33 in modulating the GC-CAFs interaction. Real-time reverse transcription-polymerase chain reaction (RT-PCR) was adopted to gauge the expression of NORAD, miR-496, and IL-33 in GC tissues and cells, and gain- or loss-of-function assays were conducted to investigate the role of them in GC. A GC cell-CAFs co-culture model was established to explore the interaction between CAFs and GCs. As exhibited, NORAD was up-regulated in GC tissues and cells, while miR-496 was remarkably down-regulated. Overexpressing NORAD substantially promoted the proliferation, migration, invasion, and EMT of GC cells and repressed cell death, while overexpressing miR-496 had the opposite effects. Additionally, NORAD enhanced the IL-33 expression and the release of IL-33 from GC cells. The dual-luciferase reporter assay confirmed that miR-496 was a target of NORAD and targeted IL-33. CAFs aggravated the malignant behaviors of GC cells as indicated by both experiments. However, NORAD knockdown in CAFs reversed CAFs-mediated promotive effects on GC cells. In conclusion, NORAD enhanced the promotive effect of CAFs in GC cells by up-regulating IL-33 and targeting miR-496, which provided new insights into the microenvironment of GC cells and CAFs.Abbreviation ANOVA: Analysis of Variance; BCA:Bicinchoninic acid; CAFs: carcinoma-associated fibroblasts; CCK-8: cell counting kit-8; ceRNA: competing endogenous RNA; DAPI: 4',6-diamidino-2-phenylindole; DMEM: Dulbecco's minimal essential medium/Ham's; ECL: enhanced chemiluminiscent; ELISA: Enzyme-Linked Immunosorbent Assay; EMT: epithelial-mesenchymal transition; FBS: fetal bovine serum; FISH:Fluorescence in situ hybridization; FITC:fluorescein isothiocyanate; FSP:fibroblast-specific protein; GAPDH: glyceraldehyde-3-phosphate dehydrogenase; GC: gastric cancer; IHC: immunohistochemistry; IL: Interleukin; lncRNA: long Noncoding RNA; miR-496: microRNA-496; MMP-14:matrix metalloproteinase-14; MUT:mutant; MYH9: myosin heavy chain 9; NFs: normal fibroblasts; NORAD: Noncoding RNA activated by DNA damage; ORF: open reading frame; PBS: phosphate-buffered saline; PMSF: Phenylmethylsulfonyl fluoride; PVDF: polyvinylidene difluoride; RIPA: Radio-Immunoprecipitation Assay; RT-PCR: Real-time reverse transcription polymerase chain reaction; S100A4:S100 calcium binding protein A4; SDS-PAGE: sodium dodecyl sulfate-polyacrylamide gel electrophoresis; sh-NC: short-hairpin RNA negative control; sh-NORAD: short-hairpin RNA of NORAD; α-SMA: α-smooth muscle actin; TBST: Tris-buffered saline with Tween-20; TGF-β1: Transforming growth factor β1; TUNEL: TdT-mediated dUTP Nick-End Labeling; TWIST1: the twist-related protein 1; VEGF-C: vascular endothelial growth factor C; WT: Wildtype.

Identification of a novel glycolysis-related signature to predict the prognosis of patients with breast cancer

Background

Breast cancer (BC) has a high incidence and mortality rate in females. Its conventional clinical characteristics are far from accurate for the prediction of individual outcomes. Therefore, we aimed to develop a novel signature to predict the survival of patients with BC.

Methods

We analyzed the data of a training cohort from the Cancer Genome Atlas (TCGA) database and a validation cohort from the Gene Expression Omnibus (GEO) database. After the applications of Gene Set Enrichment Analysis (GSEA) and Cox regression analyses, a glycolysis-related signature for predicting the survival of patients with BC was developed; the signature contained AK3, CACNA1H, IL13RA1, NUP43, PGK1, and SDC1. Furthermore, on the basis of expression levels of the six-gene signature, we constructed a risk score formula to classify the patients into high- and low-risk groups. The receiver operating characteristic (ROC) curve and the Kaplan-Meier curve were used to assess the predicted capacity of the model. Later, a nomogram was developed to predict the outcomes of patients with risk score and clinical features over a period of 1, 3, and 5 years. We further used Human Protein Atlas (HPA) database to validate the expressions of the six biomarkers in tumor and sample tissues, which were taken as control.

Results

We constructed a six-gene signature to predict the outcomes of patients with BC. The patients in the high-risk group showed poor prognosis than those in the low-risk group. The area under the curve (AUC) values were 0.719 and 0.702, showing that the prediction performance of the signature is acceptable. Additionally, Cox regression analysis revealed that these biomarkers could independently predict the prognosis of BC patients with BC without being affected by clinical factors. The expression levels of all six biomarkers in BC tissues were higher than that in normal tissues; however, AK3 was an exception.

Conclusion

We developed a six-gene signature to predict the prognosis of patients with BC. Our signature has been proved to have the ability to make an accurate prediction and might be useful in expanding the hypothesis in clinical research.

Hydroxygenkwanin Increases the Sensitivity of Liver Cancer Cells to Chemotherapy by Inhibiting DNA Damage Response in Mouse Xenograft Models

Molecules involved in DNA damage response (DDR) are often overexpressed in cancer cells, resulting in poor responses to chemotherapy and radiotherapy. Although treatment efficacy can be improved with the concomitant use of DNA repair inhibitors, the accompanying side effects can compromise the quality of life of patients. Therefore, in this study, we identified a natural compound that could inhibit DDR, using the single-strand annealing yeast-cell analysis system, and explored its mechanisms of action and potential as a chemotherapy adjuvant in hepatocellular carcinoma (HCC) cell lines using comet assay, flow cytometry, Western blotting, immunofluorescence staining, and functional analyses. We developed a mouse model to verify the in vitro findings. We found that hydroxygenkwanin (HGK) inhibited the expression of RAD51 and progression of homologous recombination, thereby suppressing the ability of the HCC cell lines to repair DNA damage and enhancing their sensitivity to doxorubicin. HGK inhibited the phosphorylation of DNA damage checkpoint proteins, leading to apoptosis in the HCC cell lines. In the mouse xenograft model, HGK enhanced the sensitivity of liver cancer cells to doxorubicin without any physiological toxicity. Thus, HGK can inhibit DDR in liver cancer cells and mouse models, making it suitable for use as a chemotherapy adjuvant.

DNA Damage Repair Profiles Alteration Characterize a Hepatocellular Carcinoma Subtype With Unique Molecular and Clinicopathologic Features

Hepatocellular carcinoma (HCC) is one of the most common malignancies and displays high heterogeneity of molecular phenotypes. We investigated DNA damage repair (DDR) alterations in HCC by integrating multi-omics data. HCC patients were classified into two heterogeneous subtypes with distinct clinical and molecular features: the DDR-activated subtype and the DDR-suppressed subtype. The DDR-activated subgroup is characterized by inferior prognosis and clinicopathological features that result in aggressive clinical behavior. Tumors of the DDR-suppressed class, which have distinct clinical and molecular characteristics, tend to have superior survival. A DDR subtype signature was ultimately generated to enable HCC DDR classification, and the results were confirmed by using multi-layer date cohorts. Furthermore, immune profiles and immunotherapy responses are also different between the two DDR subtypes. Altogether, this study illustrates the DDR heterogeneity of HCCs and is helpful to the understanding of personalized clinicopathological and molecular mechanisms responsible for unique tumor DDR profiles.

Immune implication of FAM83D gene in hepatocellular carcinoma

FAM83D has been demonstrated to contribute to tumorigenesis. However, its immune effects in hepatocellular carcinoma (HCC) have not been reported. This study aimed to identify the immune role of FAM83D in HCC. FAM83D was over-expressed in HCC and contributed to poor prognosis according to the results of data analysis based on The Cancer Genome Atlas (TCGA). Afterward, the levels of immune cells infiltration were found to be correlated with the expression level of FAM83D in HCC. Through TISIDB and cBioPortal network tools, a total of 82 FAM83D-associated genes were screened out, including 12 immunoinhibitors, 20 immunostimulators and 50 tightly co-expressed genes. TCGA cohort was divided into train set and test set on the basis of the proportion of 7:3. According to FAM83D-associated immunomodulators, a four gene predicted model was established using train set via the Cox regression analysis. Survival analysis demonstrated that the overall survival (OS) of high-risk HCC patients was poor compared with the patients in low-risk group. The reliability and predicted power of the risk-score model were identified by a receiver operating characteristic (ROC) curve. A risk-score based nomogram as well as a calibration curve, which were created could be used to anticipate patient’s 1-year, 3-year and 5-year survival probabilities. The test set was used to validate these results. Our findings showed that the FAM83D gene was related with HCC immunity. The immune marker chosen could be a promising biomarker for HCC prognosis.