The differential expression of NEDD1 in HCC
We assessed the differential expression analysis of NEDD1 in HCC patients through the Wilcoxon rank-sum. The results indicated that NEDD1 expression was obviously higher in LIHC, LUSC, et al. than in normal tissue by comparing NEDD1 expression in normal samples of Genotype-Tissue Expression (GTEx) and TCGA, and tumor samples of TCGA (Figs. 1A,1C). We further demonstrated our findings in paired samples and found that NEDD1 expression also higher in CHOL, KICH, LIHC et al. (Fig. 1B). The results showed that NEDD1 was closely associated with cancer.
The relationship between NEDD1 expression and HCC prognosis
Kaplan-Meier survival curve was conducted to explore the relevance of NEDD1 expression and OS of patients in pan cancer, which indicated that NEDD1 was relevant to poor OS in HCC(P=0.002), Kidney renal clear cell carcinoma(P=0.046) and lung adenocarcinoma (P=0.020) (Fig. 2). Subsequently, forest plot regarding OS of the univariate and multivariate Cox regression analysis in HCC demonstrated that high NEDD1 expression was an independent risk indicator for HCC patients (HR 1.643, 95%CI 1.125–2.398; P = 0.01) (Fig. 3). The finding indicated that NEDD1 was significantly related to HCC prognosis and probably be a novel prognostic marker for HCC.
Immune infiltration analysis of NEDD1 in HCC patients
Using TIMER, we assessed the association of NEDD1 expression with the abundance of immune infiltration in HCC. The results reveled that tumor purity is not directly associated with the NEDD1 expression, but a positive relevance existed between NEDD1 expression and the abundance of infiltrating B cells(r = 0.36, P = 6.08e-12), CD8+T cells(r = 0.334, P = 2.41e-10), CD4+ T cells (r =0.471, P = 2.27e-20), Macrophages (r = 0.548, P = 3.90e-28), Neutrophils (r = 0.579, P = 2.62e-32), and DCs (r = 0.492, P = 3.93e-22) in HCC(Fig. 4A). Meanwhile, the results reveled that the immune infiltrating levels were not directly relevant to OS of HCC patients (Fig. 4B). Subsequently, Spearman correlation was performed to explore the relevance of NEDD1 with the immune cells’ infiltration levels quantified by ssGSEA method. Our results showed that there were positive association of NEDD1 expression with Th2 cells and Th cells, and the negative correlation of NEDD1 with cytotoxic cells, pDC and DC was found in Fig. 4C. As show in Fig. 4D and Fig. 4E, the level of NEDD1 expression were significantly positively related to Th cells (r=0.440, P<0.001) and Th2 cells (r =0.510, P<0.001).
Previous studies suggested that immune checkpoint inhibitors are a promising treatment modality for the effective treatment of cancer (23). We used GEPIA to explore the relevance of NEED1 with the immune checkpoint genes, which indicated that NEDD1 was closely associated with immune checkpoint genes in pan cancer (Fig. 5A). Meanwhile, our findings reveled that NEDD1 expression was related to the expression of CTLA4 (r=0.260, P<0.001), CD274(PD-L1) (r=0.400, P<0.001) and PTCD1(PD-1) (r=0.270, P<0.001). (Fig. 5B,5C)
GO and KEGG analyses of NEDD1
To analyze the biological function of NEDD1, we employed the LinkedOmics database to explore the co-expressed genes of NEDD1 in LIHC from TCGA, firstly (Fig. 6A). It indicated that the expression of PPP1R12A, GTE3, et al were positively associated with NEDD1 in Fig. 6B. On the contrary, ZNHIT1, DCI, et al were negatively associated with NEDD1 in Fig. 6C. Subsequently, based on GSEA in LinkedOmics, we conducted GO analysis. The results demonstrated that NEDD1 was significantly associated with the Cell cycle, Homologous recombination and MicroRNAs in cancer signaling pathways (Fig. 6D). Meanwhile, NEDD1 plays an important part in the process of histone binding, helicase activity and is located in chromosomal region and microtubule organizing center part (Figs. 6E,6F).
In addition, on the basis of GSEA, we performed KEGG analysis demonstrating that NEDD1 expression is positively associated with the following signaling pathways, including Cell cycle, MicroRNAs in cancer and Ribosome. It also showed that NEDD1 was negatively correlated with Ribosome pathway (Figs. 7A-7D). These results demonstrated that NEDD1 played a significant impact in HCC development.
Construction and validation of prognostic nomogram in HCC patients
Further analysis revealed the diagnostic value of NEED1 in OS and disease-specific survival (DSS) of HCC. The results indicated that the AUC values of 1-year OS, 3-year OS and 5-year OS were 0.688, 0.629 and 0.628, respectively. It also suggested that the AUC values of 1-year DSS, 3-year DSS and 5-year DSS were 0.722, 0.658 and 0.634, respectively (Figs. 8C,8D). Thus, the nomogram was built to estimate the 1-, 3-, and 5-year survival probability by integrating NEDD1 expression with relevant variables, including age, gender, tumor status and other clinical factors (Fig. 8A). There are 373 original data, 208 cases with missing variable information, and the final number of enrolled samples is 165. The consistency test (Concordance, C-index) of the nomogram was 0.712 (0.663-0.761). As well, the calibration curve illustrated that there was excellent agreement between the predicted and actual survival in the prognostic nomogram, suggesting that it had good predictive value. (Fig. 8B).
Subsequently, primary liver cancer data from ICGC was employed as an external validation of the nomogram. As shown Figs.9A-9C, it demonstrated that in the ICGC database, this modal has good predictive power and can clearly distinguish between low, medium and high risks HCC patients.