HJURP regulates cell proliferation and chemo-resistance via YAP1/NDRG1 transcriptional axis in triple-negative breast cancer

Comprehensive analysis of hub genes associated with cisplatin-resistance in ovarian cancer and screening of therapeutic drugs through bioinformatics and experimental validation

BACKGROUND

To identify key genes associated with cisplatin resistance in ovarian cancer, a comprehensive analysis was conducted on three datasets from the GEO database and through experimental validation.

METHODS

Gene expression profiles were retrieved from the GEO database. DEGs were identified by comparing gene expression profiles between cisplatin-sensitive and resistant ovarian cancer cell lines. The identified genes were further subjected to GO, KEGG, and PPI network analysis. Potential inhibitors of key genes were identified through methods such as LibDock nuclear molecular docking. In vitro assays and RT-qPCR were performed to assess the expression levels of key genes in ovarian cancer cell lines. The sensitivity of cells to chemotherapy and proliferation of key gene knockout cells were evaluated through CCK8 and Clonogenic assays.

RESULTS

Results showed that 12 genes influenced the chemosensitivity of the ovarian cancer cell line SKOV3, and 9 genes were associated with the prognosis and survival outcomes of ovarian cancer patients. RT-qPCR results revealed NDRG1, CYBRD1, MT2A, CNIH3, DPYSL3, and CARMIL1 were upregulated, whereas ERBB4, ANK3, B2M, LRRTM4, EYA4, and SLIT2 were downregulated in cisplatin-resistant cell lines. NDRG1, CYBRD1, and DPYSL3 knock-down significantly inhibited the proliferation of cisplatin-resistant cell line SKOV3. Finally, photofrin, a small-molecule compound targeting CYBRD1, was identified.

CONCLUSION

This study reveals changes in the expression level of some genes associated with cisplatin-resistant ovarian cancer. In addition, a new small molecule compound was identified for the treatment of cisplatin-resistant ovarian cancer.

Variations in fluid chemical potential induce fibroblast mechano-response in 3D hydrogels

Mechanical deformation of skin creates variations in fluid chemical potential, leading to local changes in hydrostatic and osmotic pressure, whose effects on mechanobiology remain poorly understood. To study these effects, we investigate the specific influences of hydrostatic and osmotic pressure on primary human dermal fibroblasts in three-dimensional hydrogel culture models. Cyclic hydrostatic pressure and hyperosmotic stress enhanced the percentage of cells expressing the proliferation marker Ki67 in both collagen and PEG-based hydrogels. Osmotic pressure also activated the p38 MAPK stress response pathway and increased the expression of the osmoresponsive genes PRSS35 and NFAT5. When cells were cultured in two-dimension (2D), no change in proliferation was observed with either hydrostatic or osmotic pressure. Furthermore, basal, and osmotic pressure-induced expression of osmoresponsive genes differed in 2D culture versus 3D hydrogels, highlighting the role of dimensionality in skin cell mechanotransduction and stressing the importance of 3D tissue-like models that better replicate in vivo conditions. Overall, these results indicate that fluid chemical potential changes affect dermal fibroblast mechanobiology, which has implications for skin function and for tissue regeneration strategies.

YAP1-induced RBM24 promotes the tumorigenesis of triple-negative breast cancer through the β-catenin pathway

Triple-negative breast cancer (TNBC) is the most aggressive breast cancer subtype and refractory to current treatments. RBM24 is an RNA-binding protein and shows the ability to regulate tumor progression in multiple cancer types. However, its role in TNBC is still unclear. In this study, we analyzed publicly available profiling data from TNBC tissues and cells. Loss- and gain-of-function experiments were performed to determine the function of RBM24 in TNBC cells. The mechanism for RBM24 action in TNBC was investigated. RBM24 was deregulated in TNBC tissues and TNBC cells with depletion of SIPA1, YAP1, or ARID1A, three key regulators of TNBC. Compared to MCF10A breast epithelial cells, TNBC cells had higher levels of RBM24. Knockdown of RBM24 inhibited TNBC cell proliferation, colony formation, and tumorigenesis, while overexpression of RBM24 promoted aggressive phenotype in TNBC cells. YAP1 overexpression induced the expression of RBM24 and the RBM24 promoter-driven luciferase reporter. YAP1 was enriched at the promoter region of RBM24. Overexpression of RBM24 increased β-catenin-dependent transcriptional activity. Most importantly, knockdown of CTNNB1 rescued RBM24 aggressive phenotype in TNBC cells. Collectively, the YAP1/RBM24/β-catenin axis plays a critical role in driving TNBC progression. RBM24 may represent a novel therapeutic target for TNBC treatment.

N-myc downstream-regulated gene 1 can promote vasculogenic mimicry and angiogenesis in urothelial carcinoma

Urothelial carcinoma (UC) of the bladder is a common cause of cancer-related death worldwide. Vasculogenic mimicry (VM) is a process by which the malignant cells can generate vascular-like structures formed of periodic acid-Schiff (PAS) positive/CD31 negative extracellular matrix independent of angiogenesis and thus promotes tumor progression. N-myc downstream-regulated gene 1 (NDRG1) is a protein that can modulate tumor angiogenesis; however, its role in regulating tumor angiogenesis and VM formation has not been previously investigated in UC. This study aims to evaluate the role of intra-tumor microvessel density (MVD) (as a surrogate measure of angiogenesis), VM, and NDRG1 in UC and their correlation with different clinicopathologic features, then assess the correlation between them in UC. Sixty specimens of UC of the bladder were included. PAS-CD31 immunohistochemical double staining method was used to evaluate the intra-tumor MVD and VM. Immunohistochemical expression of NDRG1 was also examined. VM and NDRG1 expression were detected in 41.7% and 83.3% of UC specimens respectively. The mean of intra-tumor MVD, VM area, and NDRG1 was significantly higher in tumors with higher grade, lymphovascular invasion, and higher T stage. NDRG1 expression was positively correlated with MVD and VM. We can suggest that MVD, VM, and NDRG1 may serve as poor prognostic markers for UC. The positive correlation between NDRG1 and both MVD and VM may provide the first evidence that NDRG1 can induce tumor angiogenesis and VM in UC which may offer a novel pathway for further therapeutic strategies.

Identification of a Group of Therapeutic Targets and Prognostic Biomarker for Triple Negative Breast Cancer

INTRODUCTION

Triple-negative breast cancer (TNBC) is a highly heterogeneous disease. Mining differentially expressed genes of TNBC is helpful to explore new therapeutic targets. This study aimed to investigate diagnostic biomarker genes in TNBC compared to normal tissue. Additionally, we explored the functions and prognostic value of these key genes as well as potential targeted drugs that could affect these genes.

METHODS

Differential gene expression analysis was conducted using the R software with data from the Gene Expression Omnibus (GEO) database. Then, the identified differentially expressed genes (DEGs) were used to construct a protein-protein interaction (PPI) network using the Search Tool for the Retrieval of Interacting Genes (STRING) database and Cytoscape software. The mRNA expression levels of key genes were analyzed using the UALCAN database with data from The Cancer Genome Atlas (TCGA). Enrichment and survival analyses were performed using R software. In addition, potential compounds showing sensitivity to key genes were identified by gene set cancer analysis (GSCA).

RESULTS

Compared with normal tissues, a total of 203 DEGs were upregulated in TNBC. These DEGs participated in various biological processes including nuclear division, microtubule binding, cell cycle, and the p53 signaling pathway. Through the PPI network analysis, ten key genes were identified, among which four genes showed significant correlation with poor progression-free interval (PFI) in patients with TNBC. Moreover, the four survival-related genes were found to act as sensitive therapeutic targets.

CONCLUSION

The identified four key genes were considered new biomarkers for diagnosis and prognosis and also potential therapeutic targets for TNBC.

HJURP is recruited to double-strand break sites and facilitates DNA repair by promoting chromatin reorganization

HJURP is overexpressed in several cancer types and strongly correlates with patient survival. However, the mechanistic basis underlying the association of HJURP with cancer aggressiveness is not well understood. HJURP promotes the loading of the histone H3 variant, CENP-A, at the centromeric chromatin, epigenetically defining the centromeres and supporting proper chromosome segregation. In addition, HJURP is associated with DNA repair but its function in this process is still scarcely explored. Here, we demonstrate that HJURP is recruited to DSBs through a mechanism requiring chromatin PARylation and promotes epigenetic alterations that favor the execution of DNA repair. Incorporation of HJURP at DSBs promotes turnover of H3K9me3 and HP1, facilitating DNA damage signaling and DSB repair. Moreover, HJURP overexpression in glioma cell lines also affected global structure of heterochromatin independently of DNA damage induction, promoting genome-wide reorganization and assisting DNA damage response. HJURP overexpression therefore extensively alters DNA damage signaling and DSB repair, and also increases radioresistance of glioma cells. Importantly, HJURP expression levels in tumors are also associated with poor response of patients to radiation. Thus, our results enlarge the understanding of HJURP involvement in DNA repair and highlight it as a promising target for the development of adjuvant therapies that sensitize tumor cells to irradiation.

Long noncoding RNA MEG3 regulates cell proliferation and apoptosis by disrupting microRNA-9-5p-mediated inhibition of NDRG1 in prostate cancer

BACKGROUND

Long noncoding RNA MEG3 has been described to be involved in the regulation of gene expression and cancer progression. However, the role of lncMEG3 in prostate cancer (PCa) remains largely uncharted.

METHODS

Differential expression of lncMEG3 was identified in PCa tissues using RNA-sequencing analysis. qRT-PCR was performed to examine the level of lncMEG3. Additionally, cellular fractionation and fluorescent in situ hybridization techniques were employed to determine the localization. Subsequently, functional assays were conducted to evaluate the impact of lncMEG3 and miR-9-5p on PCa proliferation and apoptosis in vitro and in vivo. The interaction between lncMEG3 and miR-9-5p was confirmed using RNA immunoprecipitation. Moreover, luciferase reporter assays were also utilized to investigate the relationship between miR-9-5p and NDRG1.

RESULTS

We observed downregulation of lncMEG3 in PCa cells and tissues. Patients with lower levels of lncMEG3 had a higher likelihood of experiencing biochemical recurrence. Overexpression of lncMEG3 resulted in the inhibition of PCa cell proliferation and the promotion of apoptosis. Moreover, lncMEG3 is competitively bound to miR-9-5p, preventing its inhibitory effect on the target gene NDRG1. This ultimately led to the inhibition of PCa cell proliferation and the promotion of apoptosis. Furthermore, increasing lncMEG3 levels also demonstrated inhibitory effects on PCa proliferation and promotion of apoptosis in vivo.

CONCLUSIONS

Our findings uncover a crucial role for lncMEG3 in inhibiting PCa proliferation and promoting apoptosis through disruption of miR-9-5p-mediated inhibition of NDRG1.

The role and application of vesicles in triple-negative breast cancer: Opportunities and challenges

Extracellular vesicles (EVs) carry DNA, RNA, protein, and other substances involved in intercellular crosstalk and can be used for the targeted delivery of drugs. Triple-negative breast cancer (TNBC) is rich in recurrent and metastatic disease and lacks therapeutic targets. Studies have proved the role of EVs in the different stages of the genesis and development of TNBC. Cancer cells actively secrete various biomolecules, which play a significant part establishing the tumor microenvironment via EVs. In this article, we describe the roles of EVs in the tumor immune microenvironment, metabolic microenvironment, and vascular remodeling, and summarize the application of EVs for objective delivery of chemotherapeutic drugs, immune antigens, and cancer vaccine adjuvants. EVs-based therapy may represent the next-generation tool for targeted drug delivery for the cure of a variety of diseases lacking effective drug treatment.

Unleashing the Power of Yes-Associated Protein in Ferroptosis and Drug Resistance in Breast Cancer, with a Special Focus on Therapeutic Strategies

The YAP protein is a critical oncogenic mediator within the Hippo signaling pathway and has been implicated in various cancer types. In breast cancer, it frequently becomes activated, thereby contributing to developing drug-resistance mechanisms. Recent studies have underscored the intricate interplay between YAP and ferroptosis within the breast tumor microenvironment. YAP exerts a negative regulatory effect on ferroptosis, promoting cancer cell survival and drug resistance. This review offers a concise summary of the current understanding surrounding the interplay between the YAP pathway, ferroptosis, and drug-resistance mechanisms in both bulk tumor cells and cancer stem cells. We also explore the potential of natural compounds alone or in combination with anticancer therapies for targeting the YAP pathway in treating drug-resistant breast cancer. This approach holds the promise of enhancing the effectiveness of current treatments and paving the way for developing novel therapeutics.

NFE2L1 restrains ferroptosis by transcriptionally regulating HJURP and participates in the progress of oral squamous cell carcinoma

Oral squamous cell carcinoma (OSCC) is a common head and neck malignancy with increasing mortality and high recurrence. In this work, we aim to explore the functional role of NFE2 like bZIP transcription factor 1 (NFE2L1) in OSCC progression. Based on databases analysis, we found that NFE2L1 was overexpressed in OSCC tumor tissues, and elevated NFE2L1 level induced poor prognosis of OSCC patients. Our results showed that NFE2L1 is upregulated in OSCC cells and overexpression of NFE2L1 promotes cell proliferation, and reduces the sensitivity of OSCC cells to erastin-induced ferroptosis. NFE2L1 upregulation decreased the levels of Fe2+, lipid reactive oxygen species and content of malondialdehyde, and increased the level of the key negative regulator of ferroptosis, GPX4 and SLC7A11. In NFE2L1 suppressed cells, these trends were reversed. Further results of dual luciferase reporter and chromatin immunoprecipitation assays confirmed that NFE2L1 could bind to the promoter of Holliday junction recognition protein (HJURP) to increase the transcriptional activity of HJURP, thus upregulating its expression. Inhibition of HJURP attenuated the proliferation and ferroptosis inhibition in NFE2L1 upregulated cells. In vivo tumorigenicity assay further proved that NFE2L1 promotes OSCC tumor growth. In summary, NFE2L1 restrains ferroptosis by transcriptionally regulating HJURP and participates in the progress of OSCC. Thus, NFE2L1 plays a key role in OSCC development and may be a promising therapeutic target for OSCC.

Identification of metabolism-related gene signature in lung adenocarcinoma

AIM

Lung cancer is one of the most common cancers in China and has a high mortality rate. Most patients who are diagnosed have lost the opportunity to undergo surgery. Aberrant metabolism is closely associated with tumorigenesis. We aimed to identify an effective metabolism-related prediction model for assessing prognosis based on the cancer genome atlas (TCGA) and GSE116959 databases.

METHODS

TCGA and GSE116959 datasets from Gene Expression Omnibus were used to obtain lung adenocarcinoma (LUAD) data. Additionally, we captured metabolism-related genes (MRGs) from the GeneCards database. First, we extracted differentially expressed genes using R to analyze the LUAD data. We then selected the same differentially expressed genes, including 168 downregulated and 77 upregulated genes. Finally, 218 differentially expressed MRGs (DEMRGs) were included to perform functional enrichment analysis and construct a protein-protein interaction network with the help of Cytoscape and Search Tool for the Retrieval of Interacting Genes database. Cytoscape was used to visualize the intensive intervals in the network. Then univariate and Least Absolute Shrinkage and Selection Operator Cox regression analyses, which assisted in identifying the overall survival (OS)-related DEMRGs and building a 10-DEMRG prognosis model, were performed. The prognostic values, tumor immunity relevance, and molecular mechanism were further investigated. A nomogram incorporating signature, age, gender, and TNM stage was established.

RESULTS

A 10-DEMRG model was established to forecast the OS of LUAD through Least Absolute Shrinkage and Selection Operator regression analysis. This prognostic signature stratified LUAD patients into low-risk and high-risk groups. The receiver operating characteristic curve and K-M analysis indicated good performance of the DEMRGs signature at predicting OS in the TCGA dataset. Univariate and multivariate Cox regression also revealed that the DEMRGs signature was an independent prognosis factor in LUAD. We noticed that the risk score was substantially related to the clinical parameters of LUAD patients, covering age and stage. Immune analysis results showed that risk score was associated with some immune cells and immune checkpoints. Nomogram also verified the clinical value of the DEMRGs signature.

CONCLUSION

In this study, we constructed a DEMRGs signature and established a prognostic nomogram that is robust and reliable to predict OS in LUAD. Overall, the findings could help with therapeutic customization and personalized therapies.

A review on the role of NDRG1 in different cancers

NDRG1 is a member of the α/β hydrolase superfamily that resides in the cytoplasm and participates in the stress responses, hormone response, cell growth, and differentiation. Several studies have pointed to the importance of NDRG1 in the carcinogenesis. This gene has been found to be up-regulated in an array of cancer types such as bladder, esophageal squamous cell carcinoma, endometrial, lung and liver cancers, but being down-regulated in other types of cancers such as colorectal, gastric and ovarian cancers. The current study summarizes the evidence on the role of NDRG1 in the carcinogenic processes in different types of tissues.

Advances in holliday junction recognition protein (HJURP): Structure, molecular functions, and roles in cancer

Oncogenes are increasingly recognized as important factors in the development and progression of cancer. Holliday Junction Recognition Protein (HJURP) is a highly specialized mitogenic protein that is a chaperone protein of histone H3. The HJURP gene is located on chromosome 2q37.1 and is involved in nucleosome composition in the mitotic region, forming a three-dimensional crystal structure with Centromere Protein A (CENP-A) and the histone 4 complex. HJURP is involved in the recruitment and assembly of centromere and kinetochore and plays a key role in stabilizing the chromosome structure of tumor cells, and its dysfunction may contribute to tumorigenesis. In the available studies HJURP is upregulated in a variety of cancer tissues and cancer cell lines and is involved in tumor proliferation, invasion, metastasis and immune response. In an in vivo model, overexpression of HJURP in most cancer cell lines promotes cell proliferation and invasiveness, reduces susceptibility to apoptosis, and promotes tumor growth. In addition, upregulation of HJURP was associated with poorer prognosis in a variety of cancers. These properties suggest that HJURP may be a possible target for the treatment of certain cancers. Various studies targeting HJURP as a prognostic and therapeutic target for cancer are gradually attracting interest and attention. This paper reviews the functional and molecular mechanisms of HJURP in a variety of tumor types with the aim of providing new targets for future cancer therapy.

Histone Chaperones and Digestive Cancer: A Review of the Literature

BACKGROUND

The global burden of digestive cancer is expected to increase. Therefore, crucial for the prognosis of patients with these tumors is to identify early diagnostic markers or novel therapeutic targets. There is accumulating evidence connecting histone chaperones to the pathogenesis of digestive cancer. Histone chaperones are now broadly defined as a class of proteins that bind histones and regulate nucleosome assembly. Recent studies have demonstrated that multiple histone chaperones are aberrantly expressed and have distinct roles in digestive cancers.

OBJECTIVE

The purpose of this review is to present the current evidence regarding the role of histone chaperones in digestive cancer, particularly their mechanism in the development and progression of esophageal, gastric, liver, pancreatic, and colorectal cancers. In addition, the prognostic significance of particular histone chaperones in patients with digestive cancer is discussed.

METHODS

According to PRISMA guidelines, we searched the PubMed, Embase, and MEDLINE databases to identify studies on histone chaperones and digestive cancer from inception until June 2022.

RESULTS

A total of 104 studies involving 21 histone chaperones were retrieved.

CONCLUSIONS

This review confirms the roles and mechanisms of selected histone chaperones in digestive cancer and suggests their significance as potential prognostic biomarkers and therapeutic targets. However, due to their non-specificity, more research on histone chaperones should be conducted in the future to elucidate novel strategies of histone chaperones for prognosis and treatment of digestive cancer.

Prognostic significance of YAP1 expression and its association with neuroendocrine markers in resected pulmonary large cell neuroendocrine carcinoma (LCNEC)

•

It is the first study to determine the prognostic relevance of YAP1 in pulmonary LCNEC. And we found YAP1 is a prognostic factor for worse survival, especially for DFS.

•

We investigated the relationship between YAP1 and NE markers (INSM1, DLL3, NeuroD1) and found that YAP1 expression was negatively correlated with INSM1 and DLL3, but not significantly correlated with NeuroD1.

•

Our sample size is large and the clinical data is complete.

•

The exploration of the prognostic mechanism of LCNEC is of great significance to its subtype classification and stratification of treatment and prognosis.

Background

YAP1 (Yes-associated protein 1), an important effector of the Hippo pathway, acts as an oncogene and is overexpressed in various malignant tumors. However, the function and expression pattern of YAP1 in pulmonary large cell neuroendocrine carcinoma (LCNEC) have not been systematically established. This study aimed to explore the relationship between YAP1 expression and neuroendocrine differentiation markers and their prognostic significance in LCNEC.

Materials and methods

YAP1 protein and neuroendocrine markers (INSM1, NeuroD1 and DLL3) expression were examined by immunohistochemical (IHC) staining in 80 resected pulmonary LCNEC cases. The possible association between these markers and clinicopathological features was evaluated and survival analyses were performed.

Results

YAP1 was highly expressed in 25% LCNECs (20/80) , especially at a relatively higher T stage (p = 0.015). YAP1 expression was negatively correlated with INSM1 (χ²=11.53, p = 0.001) and DLL3(χ²=8.55, p = 0.004), but not with NeuroD1 (p = 0.482). For survival analyses, YAP1 expression was associated with worse disease-free survival (DFS) and overall survival (OS) (median DFS: 13 months vs. not reached (NR), p = 0.0096; median OS: not reached, NR vs. NR, p = 0.038), and was an unfavorable prognostic factor for DFS (HR:3.285; 95%CI: 1.526-7.071, p = 0.002) and OS (HR: 2.864, 95% CI: 0.932-8.796, p = 0.066).

Conclusions

YAP1 was found to be conversely correlated with neuroendocrine markers and a prognostic factor for worse survival in resected LCNEC patients, and mechanisms need to be further investigated.

The expression, clinical relevance, and prognostic significance of HJURP in cholangiocarcinoma

Background

Cholangiocarcinoma (CCA) is the malignancy originating from the biliary epithelium, including intrahepatic (iCCA), perihilar (pCCA), and distal (dCCA) CCA. The prognosis of CCA is very poor, and the biomarkers of different CCA subsets should be investigated separately. Holliday junction recognition protein (HJURP) is a key component of the pre-nucleosomal complex, which is responsible for normal mitosis. The ectopic expression of HJURP has been reported in several cancers, but not CCA.

Materials and methods

In our study, we investigated the expression of HJURP in 127 CCA patients which were composed of 32 iCCAs, 71 pCCAs, and 24 dCCAs with immunohistochemistry and divided these patients into subgroups with a low or high expression of HJURP. With chi-square test and univariate and multivariate analyses, we evaluated the clinical relevance and prognostic significance of HJURP in iCCAs, pCCAs, and dCCAs.

Results

HJURP was ectopically upregulated in CCAs compared with the para-tumor tissues based on TCGA and other mRNA-seq databases. A high expression of HJURP was correlated with low overall survival rates of iCCA and pCCA, but not in dCCA. Moreover, HJURP was an independent prognostic biomarker in both iCCA and pCCA. Patients with high HJURP were more likely to suffer CCA-related death after operation.

Conclusions

HJURP was an independent prognostic biomarker in both iCCA and pCCA, but not in dCCA. Our results provide more evidence of the molecular features of different CCA subsets and suggest that patients with high HJURP are more high-risk, which can guide more precision follow-up and treatment of CCA.

Holliday Cross-Recognition Protein HJURP: Association With the Tumor Microenvironment in Hepatocellular Carcinoma and With Patient Prognosis

Background: Hepatocellular carcinoma is the most common type of primary liver cancer, and it is associated with poor prognosis. It often fails to respond to immunotherapy, highlighting the need to identify genes that are associated with the tumor microenvironment and may be good therapeutic targets. We and others have shown that the Holliday cross-recognition protein HJURP can promote the proliferation, migration, and invasion by hepatocellular carcinoma cells, and that HJURP overexpression is associated with poor survival. Here we explored the potential relationship between HJURP and the tumor microenvironment in hepatocellular carcinoma.

Methods: We used the Immuno-Oncology-Biological-Research (IOBR) software package to analyze the potential roles of HJURP in the tumor microenvironment. Using single-cell RNA sequencing data, we identified the cell clusters expressing abundant HJURP, then linked some of these clusters to certain bioprocesses using Gene Set Enrichment Analysis (GSEA). We validated the differential expression of HJURP in tumor-infiltrating CD8⁺ T cells, sorted by flow cytometry into populations based on the expression level of PD-1. We used weighted gene co-expression network analysis (WGCNA) to identify immunity-related genes whose expression strongly correlated with that of HJURP. The function of these genes was validated based on enrichment in Gene Ontology (GO) terms, and they were used to establish a prognosis prediction model.

Results: IOBR analysis suggested that HJURP is significantly related to the immunosuppressive tumor microenvironment and was significantly related to T cells, dendritic cells, and B cells. Based on single-cell RNA sequencing, HJURP was strongly expressed in T cells, erythrocytes, and B cells from normal liver tissues, as well as in CD8⁺ T cells, dendritic cells, and one cluster of hepatocytes in hepatocellular carcinoma tissues. Malignant hepatocytes strongly expressing HJURP were associated with the downregulation of immune bioprocesses. HJURP expression was significantly higher in CD8⁺ T cells strongly expressing PD-1 than in those expressing no or intermediate levels of PD1. WGCNA identified two module eigengenes (comprising 397 and 84 genes) related to the tumor microenvironment. We identified 24 hub genes and confirmed that they were related to immune regulation. A prognostic risk score model based on expression of HJURP, PPT1, PML, and CLEC7A showed moderate ability to predict survival.

Conclusion:HJURP is associated with tumor-infiltrating immune cells, immune checkpoints, and immune suppression in hepatocellular carcinoma. HJURP-related genes involved in immune responses may be useful for predicting patient prognosis.