Identification and clinical validation of NUSAP1 as a novel prognostic biomarker in ovarian cancer
BMC Cancer volume 22, Article number: 690 (2022)
Nucleolar and spindle-associated protein 1 (NUSAP1) was shown to be involved in cell cycle regulation in cancer. However, its prognostic value and underlying mechanism in ovarian cancer remain unclear.
Oncomine, TCGA, CCLE, and UALCAN databases were used to analyze the expression level of NUSAP1 in ovarian cancer. The Kaplan–Meier plotter database was used to evaluate its prognostic value. The results from these analyses were further validated using immunohistochemical assay. The potential molecular mechanism of NUSAP1 in ovarian cancer was assessed with respect to homologous recombination repair, mismatch repair, and immunology using different databases.
Database analyses and experimental results demonstrated that NUSAP1 was highly expressed in ovarian cancer, its levels being correlated with the FIGO stage. High NUSAP1 expression was an independent risk factor affecting the prognosis of patients with epithelial ovarian cancer. Moreover, NUSAP1 was associated with cell cycle, DNA replication, homologous recombination, and p53 signaling pathway. A positive correlation was identified between the expression of NUSAP1 and BRCA1/2 in ovarian cancer. In addition, NUSAP1 was associated with the expression of DNA mismatch repair genes and immune cell infiltration.
NUSAP1 may be a valuable prognostic marker, as well as a novel biomarker for evaluating the response to immunotherapy of patients with ovarian cancer.
Ovarian cancer is a common malignant tumor of the female reproductive system. The Global Cancer Observatory 2020 statistics reported that ovarian cancer causes approximately 313,959 new cases and 207,252 deaths worldwide each year . Surgery and platinum-based chemotherapy are routine therapy options for ovarian cancer [2, 3]. Although the prognosis of early-stage ovarian cancer is favorable with treatment, nearly 75% of women present with International Federation of Gynecology and Obstetrics (FIGO) stage III–IV, which is associated to reduced 5-year survival rates . Therefore, it is important to clarify the mechanism of ovarian tumorigenesis and explore new molecular markers.
The nucleolar and spindle-associated protein 1 (NUSAP1) is a kind of microtubule- and chromatin-binding protein involved in the regulation of spindle formation and stability, chromosome segregation, and cytokinesis [5, 6]. Its expression peaks when cells transition from the G2 phase of the cell cycle to the mitotic phase . NUSAP1 is highly expressed in multiple malignant tumors, such as breast cancer , hepatocellular carcinoma , esophageal squamous cell carcinoma , and glioblastoma . High NUSAP1 expression promotes malignant biological behaviors, such as tumor invasion and metastasis, and is associated with poor prognosis in patients with astrocytoma , breast cancer , colon cancer , and cervical cancer . In addition, NUSAP1 maintains resistance to the toxic effects of anti-tubulin chemotherapeutics . Therefore, NUSAP1 has extensive research prospects in cancer. However, its prognostic value and underlying specific molecular involvement in ovarian cancer have not yet been clarified.
Our study aimed to explore the clinical significance of NUSAP1 and provide a new research direction for finding biomarkers of ovarian cancer. This study is the first to report the prognostic value of NUSAP1, as well as its involvement in signaling pathways and correlation with immune cell infiltration in ovarian cancer.
Data extraction from Oncomine, TCGA and CCLE databases
The Oncomine database (http://www.oncomine.org) was used to validate the NUSAP1 expression in cancer . The criteria were set as previously reported . Furthermore, we explored the expression of NUSAP1 in multiple cancers from The Cancer Genome Atlas (TCGA) program (https://www.cancer.gov/tcga).
The Cancer Cell Line Encyclopedia (CCLE) project (https://www.broadinstitute.org/ccle) includes 1457 cell lines and 84,434 genes, and provides information on the genetic and pharmacologic characterization of cancer models . It was used to validate the expression profiles of NUSAP1 in multiple cancer cell lines, and the correlation between the expression of NUSAP1 and BRCA1/2.
Data extraction from UALCAN and HPA databases
UALCAN (http://ualcan.path.uab.edu) uses data from Clinical Proteomic Tumor Analysis Consortium (CPTAC) to analyze protein expression . CPTAC integrates the data of genomic and proteomic to improve the understanding of the molecular basis of cancer. The database was used to validate the protein levels of NUSAP1 in ovarian cancer.
The Human Protein Atlas (HPA) database (https://www.proteinatlas.org) provides tissue and cell distribution information of proteins . The protein expression of NUSAP1 in ovarian cancer was validated using immunohistochemistry assays.
Kaplan–Meier plotter analysis
The Kaplan–Meier plotter (https://kmplot.com/analysis) can be used to assess the prognosis information of 54,675 genes in different cancer types . In this study, the Kaplan–Meier plotter was used to analyze the prognostic value of NUSAP1 by calculating its respective hazard ratios (HRs) and 95% confidence intervals (CIs). Log-rank P < 0.05 indicated statistically significant differences.
Sample source and clinical data
In this study, 114 ovarian tissues obtained between 2008 and 2012 from the Shengjing Hospital of China Medical University were embedded in paraffin. Patients included in the study did not receive radiation therapy, chemotherapy, or hormone therapy before surgery. Comprehensive staging operation and cytoreductive surgeries for ovarian cancer were performed in early and advanced stage cases, respectively. Complete clinicopathological data was available. All pathological sections were re-examined and clearly diagnosed by pathologists using the diagnostic criteria stated in the World Health Organization 2014 classification before follow-up. Cystadenomas or cystadenofibromas with > 10% borderline histology were classified as borderline ovarian tumor. The 114 samples included 78 epithelial ovarian cancer cases, 14 epithelial ovarian borderline tumors, 10 epithelial ovarian benign tumors, and 12 normal ovarian tissues. No significant difference in age among the four groups was observed (P > 0.05): the malignant tumor group had a median age of 52 years (range: 27–83 years), for the borderline tumor group was of 43 years (19–84), for the benign tumor group was of 48 years (28–61), and for the normal ovarian group was of 43 years (32–65). Within the malignant tumor group, 23 cases were grade 1, 20 cases were grade 2, and 35 cases were grade 3. According to the 2009 FIGO standards, 34 patients had stage I-II and 44 had stage III-IV ovarian cancer. Evaluation of lymph node metastases revealed metastases in 18 cases, no metastasis in 33 cases, and no lymphadenectomy in 27 cases. Within the borderline tumor group, there were seven samples of serous borderline tumor and seven samples of mucinous borderline tumor. Within the benign tumor group, there was one sample of benign serous tumor and nine samples of benign mucinous tumor.
Immunohistochemistry (IHC) assay
IHC assay and scoring method were carried out according to previously published methods [18, 23]. The samples were dewaxed with xylene and rehydrated with sequential solutions of decreasing ethanol concentrations. Thereafter, the immunostaining permeation solution was used to permeate the membrane. After inhibiting the endogenous peroxidase activity, antigen retrieval was performed with heat treatment. After cooling, the samples were blocked with goat serum for 30 min at 37 °C. A working concentration of 1:100 polyclonal antibody against NUSAP1 (12024–1-AP; Proteintech, Wuhan, China) was used to evaluate NUSAP1 expression in ovarian cancer. After overnight incubation at 4 °C, biotinylated secondary antibodies were added for 30 min at 37 °C. 3,3′-Diaminobenzidine was used to visualize the antibody complexes. The sections were then stained with hematoxylin, dehydrated, and examined. Positive control and negative control were set for each batch. All sections were visualized under a Nikon Eclipse Ci microscope (Tokyo, Japan) and captured by NIS-Elements F software. The final scores resulted from multiplying the staining and percentage of stained cells scores: 0–2 scores (–), 3–4 scores ( +), 5–8 scores (+ +), and 9–12 scores (+ + +). Two senior pathologists scored each tissue slice independently to reduce errors.
Gene set enrichment analysis (GSEA)
According to the median expression of NUSAP1, the TCGA-OV data set was divided into high and low expression groups. GSEA version 3.0 was used to analyze the data . We downloaded the C2.cp.kegg.v6.1.symbols.gmt data cluster from the Molecular Signatures Database, and performed enrichment analysis. The number of random combinations was set to 1,000.
Metascape database analysis
Metascape (http://metascape.org) is an online platform that provides gene enrichment, biological process annotation, and protein interaction network information . Metascape was used to evaluate the functional and pathway enrichment of NUSAP1 and its co-expressed genes.
cBioPortal database analysis
The cBioPortal database (https://www.cbioportal.org) is an online platform that provides information on DNA copy numbers, DNA methylation, mRNA and microRNA expression, and nonsynonymous mutations . The database was used to analyze the correlation between NUSAP1 and BRCA1/2 in ovarian cancer.
Immune infiltration analysis
The TIMER database (http://timer.cistrome.org) contains the information of immune cell infiltration in different cancer types . We evaluated the correlation of NUSAP1 expression with immune infiltrates, and the effect of NUSAP1 copy number variations on the infiltration of immune cells. Furthermore, the TISIDB database (http://cis.hku.hk/TISIDB), which integrates multiple public databases to provide the information about tumor-immune system interactions, was employed to explore the distribution of NUSAP1 expression across molecular subtypes and immune subtypes in ovarian cancer .
SPSS 22.0 software (IBM Corporation, Armonk, NY, USA) was used to analyze the data. Chi-squared and Fisher’s exact probability tests were performed to analyze the count data, and Student’s t-test was used to analyze the measurement data. Survival curves were analyzed via Kaplan–Meier and log-rank tests. The Cox regression model was used to analyze the relationships between NUSAP1 expression and clinicopathological parameters. P < 0.05 was the cutoff value established for statistical significance.
Validation of NUSAP1 expression level in different databases
We initially evaluated NUSAP1 expression in multiple cancers through the Oncomine database (Fig. 1a). Of these studies, 92 showed that the expression of NUSAP1 was significantly upregulated, including three studies in ovarian cancer. The difference in expression between multiple cancer and normal tissues was evaluated through the TCGA database (Fig. 1b). Ovarian cancer datasets showed that NUSAP1 expression was significantly upregulated in ovarian cancer. The data of each tumor cell line were downloaded from the CCLE database and similar results were obtained (Fig. 1c). Furthermore, the UALCAN database was used to analyze the expression of NUSAP1 in ovarian cancer and its correlation with clinical parameters. NUSAP1 expression was significantly increased in ovarian cancer compared with that of normal tissues, especially in FIGO stage III (Fig. 1d, e). In addition, the level of NUSAP1 was higher in Caucasian people and patients aged between 61 and 80 years old compared with other subgroups (Fig. 1f, g). The HPA database evaluated the level of NUSAP1 in ovarian serous carcinoma, mucinous carcinoma, endometrioid carcinoma, and normal tissues using immunohistochemical data. There were different degrees of staining in ovarian cancer tissues, but no staining in normal tissues. Representative photomicrographs of various tissue samples were selected for protein expression validation (Fig. 1h).
Relationship between NUSAP1 expression level and ovarian cancer prognosis
The correlation between NUSAP1 expression and patient prognosis was analyzed using the Kaplan–Meier plotter database. Results showed that high NUSAP1 expression was correlated with poor overall survival (OS) (Fig. 2a). However, the analysis of progression-free survival (PFS) did not yield similar results (Fig. 2b). NUSAP1 upregulation indicated poor OS in patients with FIGO stage I+II (Fig. 2c, d). In addition, NUSAP1 upregulation was correlated with poor OS in patients with grades 3 and 4 (Fig. 2e, f, g, h). A significant correlation was observed between upregulation of NUSAP1 and poor OS in patients harboring TP53 mutations (Fig. 2i, j). Upregulation of NUSAP1 was correlated with poor OS in ovarian serous cancer, but not in ovarian endometrioid cancer (Fig. 2k, l). Further analysis showed that high NUSAP1 expression conferred poor OS in high-grade serous ovarian cancer (HGSC) (Fig. 2m), especially in the early-stage (Fig. 2n, o).
We further analyzed the correlation between the expression of NUSAP1 and the prognosis of chemotherapy after optimal surgery. For all patients with ovarian cancer, NUSAP1 could not predict the prognosis of receiving platinum-paclitaxel combination or platinum chemotherapy (Fig. 3a, b). Surprisingly, serous ovarian cancer patients with high NUSAP1 expression had a poor prognosis after receiving these chemotherapy regimens (Fig. 3c, d). Similar results were obtained when patients were treated with docetaxel (Fig. 3e). However, NUSAP1 had no effect on the prognosis of patients with serous ovarian cancer treated with gemcitabine, topotecan, or avastin (Fig. 3f, g, h). As for patients with HGSC, upregulation of NUSAP1 was correlated with poor prognosis after platinum or docetaxel therapy (Fig. 3i, j, k). Therefore, NUSAP1 upregulation may be related to platinum and docetaxel-based chemotherapy resistance in patients with serous ovarian cancer.
NUSAP1 is highly expressed in ovarian epithelial malignancies
NUSAP1 expression in our tissues cohort was evaluated by IHC assay. Results showed that NUSAP1 was mainly expressed in the cytoplasm and nucleus. The positive rate and high-positive rate of the malignant tumor group were significantly higher than those of the borderline, benign, and normal group. No significant differences in NUSAP1 expression were observed among borderline, benign and normal tissues (Fig. 4a; Table 1).
NUSAP1 expression in epithelial ovarian malignant tumors is related to the FIGO stage
Then, we evaluated the association between NUSAP1 expression and clinicopathological characteristics (Table 2). Seventy-eight patients with ovarian epithelial malignant tumors were divided into two groups according to the NUSAP1 expression level. High NUSAP1 expression had a notable correlation with FIGO stage. The high-positive rate of NUSAP1 expression was significantly higher in the patients with FIGO III-IV compared with that in the patients with FIGO I-II. NUSAP1 showed the highest positive rate in ovarian clear cell carcinoma. In addition, the high-positive rate of NUSAP1 in HGSC was higher than that in low-grade serous ovarian cancer. However, NUSAP1 level was not significantly correlated with tumor grade, lymph node metastasis, residual tumor size, or pathological type.
High NUSAP1 expression is an independent risk factor for the prognosis of patients with ovarian epithelial malignant tumors
We followed-up 78 patients with ovarian epithelial malignant tumors (last follow-up time was September 30, 2017), 65 of whom (83.33%) had complete clinical information, whereas 13 cases (16.67%) were lost to follow-up. Kaplan–Meier and log-rank tests showed that the 5-year survival rates of patients with high NUSAP1 expression, FIGO III-IV, and residual tumor size > 1 cm were significantly decreased (Fig. 4b). In the multivariate Cox′s proportional hazards model analysis, NUSAP1 expression level, FIGO stage, and residual tumor size were found to be independent prognostic factors (Table 3).
To investigate the significance of NUSAP1 in epithelial ovarian cancer of different stages, we analyzed the prognostic significance stratified by FIGO stage. Both Kaplan–Meier survival analysis and Cox regression analysis showed that high expression of NUSAP1 was associated with shorter OS in early-stage, but not in advanced-stage epithelial ovarian cancer (Fig. 4b, Table 4). Therefore, the prognostic significance of NUSAP1 mainly exists in early-stage.
Ovarian serous cancer is the most common pathological type of ovarian cancer. Therefore, we focused on the prognostic value of NUSAP1 in the specific pathological subtype. Kaplan–Meier survival analysis showed that for patients with serous and non-serous ovarian cancer, high expression of NUSAP1 was associated with poor OS (Fig. 4b). By univariate analysis, NUSAP1 expression and residual tumor size predicted poor OS in serous ovarian cancer. However, for patients with HGSC, the Cox regression model only suggested the correlation between residual tumor size and prognosis (Table 5).
Molecular mechanism of NUSAP1 involved in ovarian cancer development
To explore the relevant pathways in which NUSAP1 is involved, we conducted gene set enrichment analysis. The results indicated that the cell cycle, DNA replication, homologous recombination, and p53 signaling pathway were enriched in NUSAP1 high-expressing samples (Fig. 5a). Calcium signaling pathway, extracellular matrix-receptor interaction, focal adhesion, and the Wnt signaling pathway were enriched in NUSAP1 low-expressing samples (Fig. 5b). The Metascape database also confirmed that NUSAP1 and its co-expressed genes were involved in the regulation of DNA repair (GO: 0006281), replication fork (GO: 0005657), and base excision repair (ko03410) (Fig. 5c).
Homologous recombination repair (HRR) is a key repair mechanism for DNA double-strand breaks, in which BRCA1/2 participates as key tumor suppressor genes. The correlation between NUSAP1 and BRCA1/2 was analyzed through the CCLE and cBioPortal databases. Both databases showed a positive correlation between NUSAP1 and BRCA1/2 expression (Fig. 5d, e). In addition, studies confirmed that DNA methyltransferase (DNMT) inhibitors induced a BRCAness phenotype. We further analyzed the correlation between NUSAP1 and methyltransferase, and found that NUSAP1 was positively correlated with DNMT1, DNMT3A, and DNMT3B (Fig. 5f).
The mismatch repair (MMR) function is responsible for repairing DNA base mismatches that occur during replication. Mutation of MMR-related genes will cause DNA replication errors, thereby resulting in genetic instability and higher somatic mutations. We analyzed the correlation between NUSAP1 and key MMR genes, among them, NUSAP1 was positively correlated with MSH2 and MSH6 (Fig. 5g).
Correlation between NUSAP1 expression and immune cells in ovarian cancer
In order to explore the predictive value of NUSAP1 in ovarian cancer immunotherapy response, we analyzed the relationship between NUSAP1 and immune infiltration. The Timer database showed that high expression of NUSAP1 indicated high infiltrating abundances of neutrophils, macrophages, B cells, CD4 + T cells, and myeloid dendritic cells (Fig. 6a). We further compared the tumor infiltration level in ovarian cancer with different somatic copy number variations of NUSAP1 (Fig. 6b). Compared with arm-level gain, normal copy number of NUSAP1 was correlated with increased infiltration of CD8 + T cells, although the correlation between NUSAP1 expression and CD8 + T cells infiltration level was not significant.
Furthermore, multivariable hazards models were used to evaluate the impact of NUSAP1 expression in the presence of varying immune cells. NUSAP1 had 1.2-fold higher risk on PFS (Fig. 6c).
Moreover, we explored the distribution of NUSAP1 expression across molecular subtypes in ovarian cancer using the TISIDB database. TCGA ovarian samples were classified into differentiated, immunoreactive, mesenchymal, and proliferative subtypes. The results showed that the expression of NUSAP1 was relatively high in immunoreactive and proliferative subtypes (Fig. 6d). Further analysis of immune subtypes showed that the expression of NUSAP1 was the lowest in the inflammatory (C3) subtype (Fig. 6e). Hence, NUSAP1 could be a potential biomarker for immunotherapy response prediction.
Ovarian cancer currently lacks effective early screening and diagnostic methods and is prone to recurrence, metastasis, and drug resistance, resulting in no significant improvement in the 5-year OS of affected patients. Therefore, finding valuable tumor markers and exploring their contribution for ovarian cancer development are crucial for the diagnosis, treatment, and prognosis of patients with ovarian cancer. In this study, the expression of NUSAP1 was validated though multiple databases. IHC assay confirmed that NUSAP1 is highly expressed in ovarian cancer and is an independent risk factor for patient survival prognosis. Furthermore, database analyses indicated molecular mechanisms through which NUSAP1 may be involved in ovarian cancer development.
Numerous studies have confirmed that NUSAP1 is upregulated in breast cancer , hepatocellular carcinoma , and other malignant tumors, and is closely related to malignant biological behaviors, such as cell cycle regulation, invasion, and metastasis. NUSAP1 silencing inhibited DNMT1 expression, leading to the inhibition of liver and colorectal cancer progression . NUSAP1 promotes the nuclear translocation of GLI family zinc finger 1 and activates the Hedgehog signaling pathway, resulting in increased aggressiveness of astrocytoma . Li et al. found that NUSAP1 promotes metastasis and induces cancer stem cell-like properties of cervical cancer cells by activating Wnt/β-catenin signals . Moreover, Xu et al. demonstrated that inhibition of NUSAP1 suppresses cell growth and metastasis by regulating the BTG2/PI3K/Akt signals in non-small-cell lung cancer . These studies indicate that NUSAP1 is closely related to the occurrence of various tumors, and its high expression may promote the malignant biological behavior of tumors. Recently, Shen et al. proposed six hub genes, including NUSAP1, to predict the development and prognosis of ovarian cancer based on the Gene Expression Omnibus data. However, this study did not conduct experiments to verify the results of the bioinformatics analysis and evaluate the contribution of each hub gene separately . Zhang et al. showed that knocking down NUSAP1 in ovarian cancer cells promotes apoptosis and affects the cell cycle distribution . However, its specific mechanism of function in ovarian cancer has not been clarified. Herein, through the joint analysis of Oncomine, TCGA, and CCLE databases, we found that NUSAP1 is highly expressed in multiple malignant tumors, and its protein level is correlated with FIGO stage. Results of IHC assay showed that the positive rate of NUSAP1 in ovarian malignant tumors is significantly higher than that in other groups. In addition, a significant correlation between NUSAP1 level and FIGO stage was observed. Further analysis of the expression of NUSAP1 in different histological types of ovarian cancer showed that NUSAP1 had the highest positive rate in ovarian clear cell carcinoma. However, due to the limitation of the number of clear cell carcinoma cases, it is necessary to expand the sample size to further verify the conclusion.
High NUSAP1 expression predicts poor prognosis in astrocytoma , breast cancer , esophageal squamous cell carcinoma , and colon cancer . However, opposite results were obtained in the prognostic analysis of cervical cancer [15, 32]. In addition, studies have suggested a correlation between NUSAP1 and chemotherapeutic resistance. In oral squamous cell carcinoma, NUSAP1 downregulation enhances the anti-tumor effect of paclitaxel by activating apoptotic pathways . Moreover, Zhang et al. reported that knockdown of NUSAP1 increased susceptibility to epirubicin in invasive breast cancer . Recent investigations showed that NUSAP1 stabilizes the ATR by sumoylation via its SAP domain, thereby promoting chemotherapeutic resistance to temozolomide and doxorubicin . In this study, we followed up cancer patients for more than 5 years to explore the prognostic value of NUSAP1. The result of Cox regression analysis, which was consistent with the results of bioinformatics analysis, showed that high NUSAP1 expression conferred poor prognosis of patients with ovarian cancer. In addition, NUSAP1 showed an even more significant correlation with prognosis in patients with early-stage. Ovarian serous cancer is the most common pathological type of ovarian cancer, and HGSC is a leading cause of ovarian cancer-related death . Five years after diagnosis, advanced-stage HGSC cases had a survival rate of 32.1% . The results of bioinformatics analysis and immunohistochemistry statistics indicated that high expression of NUSAP1 was correlated with poor OS of patients with serous ovarian cancer. However, for the patients with HGSC, our clinical sample analysis results were not consistent with the database, which may be limited by the number of cases and the heterogeneity of the study cohort.
We also explored the molecular mechanism of NUSAP1 contributing for the tumorigenesis and progression of ovarian cancer. Database analysis showed that NUSAP1 is involved in the regulation of the DNA repair, homologous recombination, base excision repair, replication fork, cell cycle, DNA replication, and p53 signaling pathway. Studies have confirmed that these pathways are closely related to ovarian cancer. NUSAP1 has been confirmed to help generate microtubules near chromatin, leading to chromatin-induced spindle formation, which is essential for cell proliferation . In particular, Mills et al. identified NUSAP1 as a substrate for SCF-type E3 ubiquitin ligase during S/G2 phase of the cell cycle and observed a cell cycle regulatory interaction between NUSAP1 and a SUMO (small ubiquitin-like modifier) E3 ligase complex through mass spectrometry-based proteomics [16, 39]. In addition, inhibition of NUSAP1 expression reduced cyclin A2 and cyclin B1 levels, and induced cell cycle arrest [10, 15, 40].
Several types of DNA damage repair (DDR) pathways, including nucleotide excision repair, base excision repair, MMR, homologous recombination, and non-homologous end joining, can detect DNA damage and initiate the process of cell repair. Hence, alterations in DDR pathways lead to genomic instability and cancer development . As an error-free repair mechanism, HRR is a key repair method for DNA double-strand breaks. Epigenetic modifications and mutations in HRR-related genes can cause homologous recombination deficiency, thereby enhancing the sensitivity of patients to PARP inhibitors . Studies showed that NUSAP1 knockdown reduces the sumoylation level of ATR and cause DNA damage in glioblastoma multiforme cells . In addition, NUSAP1, similar to BRCA1, plays a role in DNA double-stand break repair via the homologous recombination and single-strand annealing pathways in breast cancer, and it protects BRCA1 from proteasome-mediated degradation . The combination of NUSAP1 and BRCA1 improves the prognostic power in triple-negative breast cancer compared with the calculation based on age, menstrual status, and lymph node status . We found that NUSAP1 is correlated with DNA repair, replication fork and homologous recombination in ovarian cancer. In addition, both the CCLE and cBioPortal database showed a positive correlation between the expression of NUSAP1 and BRCA1 in ovarian cancer. Therefore, NUSAP1 in ovarian cancer may participate in HRR by regulating the levels of ATR and BRCA1/2; thus, low expression of NUSAP1 in ovarian cancer implies sensitivity to PARP inhibitors. Studies have shown that the acquisition of a drug resistant phenotype is related to the accumulation of epigenetic variations in tumor suppressor and DNA repair genes . DNMTs are a class of epigenetic regulatory enzymes, of which DNMT1 is involved in maintaining the methylation status after DNA synthesis and DNMT3A/B are involved in the de novo methylation of DNA . It was reported that DNMT inhibitors induce a BRCAness phenotype through the downregulation of homologous recombination and NHEJ genes, thereby sensitizing non-small cell lung cancer cells to PARP inhibitors . In addition, combination of DNMT inhibitors and PARP inhibitors promotes the accumulation of reactive oxygen species (ROS), thereby sensitizing breast and ovarian cancer cells to PARP inhibitors in a ROS-cAMP/PKA-dependent manner . Therefore, we further evaluated the correlation between NUSAP1 and methyltransferases, and found that NUSAP1 was positively correlated with DNMT1, DNMT3A and DNMT3B. Therefore, patients with high expression of NUSAP1 can be treated with the combination of DNMT inhibitors and PARP inhibitors to enhance treatment efficacy and overall survival.
The MMR system is responsible for repairing DNA base mismatches ; thus, MMR gene mutations and hypermethylation in the promoter region will result in microsatellite instability (MSI), which are mostly caused by inactivation of MLH1, PMS2, MSH2, or MSH6 . Solid tumors with deficient MMR usually have immunogenicity and extensive T-cell infiltration, which make them more responsive to the treatment with immune checkpoint inhibitors (ICIs) [48, 49]. Although the response rate to ICIs in ovarian cancer is usually lower than in other cancer types, a proportion of ovarian cancers exhibiting MSI show an increased number of CD3 + and CD8 + tumor infiltrating lymphocytes [50, 51]. Herein, we analyzed the correlation between NUSAP1 and key MMR genes, which revealed that NUSAP1 is correlated with MSH2 and MSH6. Furthermore, we analyzed the correlation between NUSAP1 and immune infiltration level. NUSAP1 showed a strongest correlation with neutrophil tumor infiltration. Neutrophils are the most abundant white blood cells and are closely related to inflammation; however, the role and characteristics of neutrophils in cancer are controversial. Studies showed that neutrophils have more anti-tumor phenotype during early stages of tumor progression, whereas become more tumorigenic during aggressive stages [52, 53]. Neutrophils maintain tumor growth by suppressing T cell activation, accelerating gene mutation, tumor cell proliferation, angiogenesis and metastasis . It has been observed that high infiltration with tumor-associated neutrophils is associated with poor prognosis in multiple cancers . Comprehensive meta-analysis showed that high neutrophil-to-lymphocyte ratio is associated with shorter PFS and OS in ovarian cancer . In addition, the infiltration of CD8 + T cells in tumors was found to be associated with favorable prognosis in ovarian cancer . Therefore, high NUSAP1 expression may cause increased infiltration of neutrophils and decreased infiltration of CD8 + T cells, leading to the progression and poor prognosis of ovarian cancer. Altogether, NUSAP1 may represent a new biomarker for predicting the efficacy of immunotherapy in ovarian cancer.
In summary, database analyses and IHC assay confirmed that NUSAP1 is highly expressed in ovarian cancer. High NUSAP1 level is an independent risk factor affecting the survival and prognosis of patients with ovarian cancer, especially for patients with early-stage or serous ovarian cancer. Moreover, NUSAP1 expression is closely related to response to chemotherapy. Through multi-omics analysis, we confirmed that high NUSAP1 expression is closely related to signaling pathways that influence the development of cancer, such as cell cycle, DNA replication, homologous recombination, and p53 signals. In addition, a correlation was found between the expression of NUSAP1 and immune cell infiltration. Therefore, this study highlights the significant value of NUSAP1 for predicting prognosis and response to chemotherapy and immunotherapy in ovarian cancer. Further investigation using in vitro and in vivo models are needed to validate the biological significance of NUSAP1 in different pathological subtypes of ovarian cancer and to elucidate the molecular mechanism of NUSAP1 in regulating DNA repair. In addition, the ability of NUSAP1 as a biomarker to predict the efficacy of PARP inhibitors alone or in combination with DNMT inhibitions warrants further evaluation, which will assist in the development of future individualized therapeutic strategies for ovarian cancer.
Availability of data and materials
The datasets generated by the TCGA Research Network during the current study are available on the website, https://www.cancer.gov/tcga. Other data generated or analysed during this study are included in this published article.
Cancer Cell Line Encyclopedia
Clinical Proteomic Tumor Analysis Consortium
DNA damage repair
International Federation of Gynecology and Obstetrics
Gene set enrichment analysis
High-grade serous ovarian cancer
Human Protein Atlas
Homologous recombination repair
Immune checkpoint inhibitor
DNA mismatch repair
Non-homologous end joining
Nucleolar and spindle-associated protein 1
Reactive oxygen species
Small ubiquitin-like modifier
The Cancer Genome Atlas
Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J Clin. 2021;71(3):209–49.
Bristow RE, Tomacruz RS, Armstrong DK, Trimble EL, Montz FJ. Survival effect of maximal cytoreductive surgery for advanced ovarian carcinoma during the platinum era: a meta-analysis. J Clin Oncol. 2002;20:1248–59.
Kyrgiou M, Salanti G, Pavlidis N, Paraskevaidis E, Ioannidis JP. Survival benefits with diverse chemotherapy regimens for ovarian cancer: meta-analysis of multiple treatments. J Natl Cancer Inst. 2006;98:1655–63.
Lheureux S, Braunstein M, Oza AM. Epithelial ovarian cancer: Evolution of management in the era of precision medicine. CA Cancer J Clin. 2019;69:280–304.
Ribbeck K, Raemaekers T, Carmeliet G, Mattaj IW. A role for NuSAP in linking microtubules to mitotic chromosomes. Curr Biol. 2007;17:230–6.
Ribbeck K, Groen AC, Santarella R, Bohnsack MT, Raemaekers T, Kocher T, et al. NuSAP, a mitotic RanGTP target that stabilizes and cross-links microtubules. Mol Biol Cell. 2006;17:2646–60.
Raemaekers T, Ribbeck K, Beaudouin J, Annaert W, Van Camp M, Stockmans I, et al. NuSAP, a novel microtubule-associated protein involved in mitotic spindle organization. J Cell Biol. 2003;162:1017–29.
Kotian S, Banerjee T, Lockhart A, Huang K, Catalyurek UV, Parvin JD. NUSAP1 influences the DNA damage response by controlling BRCA1 protein levels. Cancer Biol Ther. 2014;15:533–43.
Wang Y, Ju L, Xiao F, Liu H, Luo X, Chen L, et al. Downregulation of nucleolar and spindle-associated protein 1 expression suppresses liver cancer cell function. Exp Ther Med. 2019;17:2969–78.
Guan C, Liu Z, Lu C, Xiao M, Shi H, Ni R, et al. Nucleolar spindle-associated protein 1 promotes tumorigenesis and predicts poor prognosis in human esophageal squamous cell carcinoma. J Cell Biochem. 2019. https://doi.org/10.1002/jcb.28452.
Qian Z, Li Y, Ma J, Xue Y, Xi Y, Hong L, et al. Prognostic value of NUSAP1 in progression and expansion of glioblastoma multiforme. J Neurooncol. 2018;140:199–208.
Wu X, Xu B, Yang C, Wang W, Zhong D, Zhao Z, et al. Nucleolar and spindle associated protein 1 promotes the aggressiveness of astrocytoma by activating the Hedgehog signaling pathway. J Exp Clin Cancer Res. 2017;36:127.
Chen L, Yang L, Qiao F, Hu X, Li S, Yao L, et al. High Levels of Nucleolar Spindle-Associated Protein and Reduced Levels of BRCA1 Expression Predict Poor Prognosis in Triple-Negative Breast Cancer. PLoS ONE. 2015;10:e0140572.
Liu Z, Guan C, Lu C, Liu Y, Ni R, Xiao M, et al. High NUSAP1 expression predicts poor prognosis in colon cancer. Pathol Res Pract. 2018;214:968–73.
Li H, Zhang W, Yan M, Qiu J, Chen J, Sun X, et al. Nucleolar and spindle associated protein 1 promotes metastasis of cervical carcinoma cells by activating Wnt/β-catenin signaling. J Exp Clin Cancer Res. 2019;38:33.
Emanuele MJ, Elia AE, Xu Q, Thoma CR, Izhar L, Leng Y, et al. Global Identification of Modular Cullin-Ring Ligase Substrates. Cell. 2011;147:459–74.
Rhodes DR, Kalyana-Sundaram S, Mahavisno V, Varambally R, Yu J, Briggs BB, et al. Oncomine 3.0: Genes, pathways, and networks in a collection of 18,000 cancer gene expression profles. Neoplasia. 2007;9:166–80.
Gou R, Zhu L, Zheng M, Guo Q, Hu Y, Li X, et al. Annexin A8 can serve as potential prognostic biomarker and therapeutic target for ovarian cancer: based on the comprehensive analysis of Annexins. J Transl Med. 2019;17:275.
Barretina J, Caponigro G, Stransky N, Venkatesan K, Margolin AA, Kim S, et al. The Cancer Cell Line Encyclopedia enables predictive modelling of anticancer drug sensitivity. Nature. 2012;483:603–7.
Chen F, Chandrashekar DS, Varambally S, Creighton CJ. Pan-cancer molecular subtypes revealed by mass-spectrometry-based proteomic characterization of more than 500 human cancers. Nat Commun. 2019;10:5679.
Uhlen M, Zhang C, Lee S, Sjöstedt E, Fagerberg L, Bidkhori G, et al. A pathology atlas of the human cancer transcriptome. Science. 2017;357:eaan2507.
Gyorffy B, Lánczky A, Szállási Z. Implementing an online tool for genome-wide validation of survival-associated biomarkers in ovarian cancer using microarray data from 1287 patients. Endocr Relat Cancer. 2012;19:197–208.
Gou R, Hu Y, Liu O, Dong H, Gao L, Wang S, et al. PGK1 Is a Key Target for Anti-Glycolytic Therapy of Ovarian Cancer: Based on the Comprehensive Analysis of Glycolysis-Related Genes. Front Oncol. 2021;11:682461.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102:15545–50.
Zhou Y, Zhou B, Pache L, Chang M. Metascape provides a biologist-oriented resource for the analysis 600 of systems-level datasets. Nat Commun. 2019;10:1523.
Gao J, Aksoy BA, Dogrusoz U, Dresdner G, Gross B, Sumer SO, et al. Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Sci Signal. 2013;6:pl1.
Li T, Fan J, Wang B, Traugh N, Chen Q, Liu JS, et al. TIMER: A Web Server for Comprehensive Analysis of Tumor-Infiltrating Immune Cells. Cancer Res. 2017;77:e108–10.
Ru B, Wong CN, Tong Y, Yi Zhong J, Wa Zhong SS, Wu WC, et al. TISIDB: an integrated repository portal for tumor-immune system interactions. Bioinformatics. 2019;35:4200–2.
Xu Z, Wang Y, Xiong J, Cui F, Wang L, Peng H. NUSAP1 knockdown inhibits cell growth and metastasis of non-small-cell lung cancer through regulating BTG2/PI3K/Akt signaling. J Cell Physiol. 2020;235:3886–93.
Shen J, Yu S, Sun X, Yin M, Fei J, Zhou J. Identification of key biomarkers associated with development and prognosis in patients with ovarian carcinoma: evidence from bioinformatic analysis. J Ovarian Res. 2019;12:110.
Zhang Y, Huang K, Cai H, Chen S, Sun D, Jiang P. The role of nucleolar spindle-associated protein 1 in human ovarian cancer. J Cell Biochem. 2020;121:4397–405.
Xie Q, Ou-Yang W, Zhang M, Wang H, Yue Q. Decreased Expression of NUSAP1 Predicts Poor Overall Survival in Cervical Cancer. J Cancer. 2020;11:2852–63.
Okamoto A, Higo M, Shiiba M, Nakashima D, Koyama T, Miyamoto I, et al. Down-Regulation of Nucleolar and Spindle-Associated Protein 1 (NUSAP1) Expression Suppresses Tumor and Cell Proliferation and Enhances Anti-Tumor Effect of Paclitaxel in Oral Squamous Cell Carcinoma. PLoS ONE. 2015;10:e0142252.
Zhang X, Pan Y, Fu H, Zhang J. Nucleolar and Spindle Associated Protein 1 (NUSAP1) Inhibits Cell Proliferation and Enhances Susceptibility to Epirubicin In Invasive Breast Cancer Cells by Regulating Cyclin D Kinase (CDK1) and DLGAP5 Expression. Med Sci Monit. 2018;24:8553–64.
Zhao Y, He J, Li Y, Lv S, Cui H. NUSAP1 potentiates chemoresistance in glioblastoma through its SAP domain to stabilize ATR. Signal Transduct Target Ther. 2020;5:44.
Pan J, Hu Y, Sun S, Chen L, Schnaubelt M, Clark D, et al. Glycoproteomics-based signatures for tumor subtyping and clinical outcome prediction of high-grade serous ovarian cancer. Nat Commun. 2020;11:6139.
Peres LC, Cushing-Haugen KL, Köbel M, Harris HR, Berchuck A, Rossing MA, et al. Invasive Epithelial Ovarian Cancer Survival by Histotype and Disease Stage. J Natl Cancer Inst. 2019;111:60–8.
Vanden Bosch A, Raemaekers T, Denayer S, Torrekens S, Smets N, Moermans K, et al. NuSAP is essential for chromatin-induced spindle formation during early embryogenesis. J Cell Sci. 2010;123:3244–55.
Mills CA, Suzuki A, Arceci A, Mo JY, Duncan A, Salmon ED, et al. Nucleolar and spindle-associated protein 1 (NUSAP1) interacts with a SUMO E3 ligase complex during chromosome segregation. J Biol Chem. 2017;292:17178–89.
Roy S, Hooiveld GJ, Seehawer M, Caruso S, Heinzmann F, Schneider AT, et al. microRNA 193a–5p Regulates Levels of Nucleolar- and Spindle-Associated Protein 1 to Suppress Hepatocarcinogenesis. Gastroenterology. 2018;155:1951–1966.e26.
McMullen M, Karakasis K, Madariaga A, Oza AM. Overcoming Platinum and PARP-Inhibitor Resistance in Ovarian Cancer. Cancers (Basel). 2020;12:1607.
Gou R, Dong H, Lin B. Application and reflection of genomic scar assays in evaluating the efficacy of platinum salts and PARP inhibitors in cancer therapy. Life Sci. 2020;261:118434.
Zhang J, Yang C, Wu C, Cui W, Wang L. DNA Methyltransferases in Cancer: Biology, Paradox, Aberrations, and Targeted Therapy. Cancers (Basel). 2020;12:2123.
Abbotts R, Topper MJ, Biondi C, Fontaine D, Goswami R, Stojanovic L, et al. DNA methyltransferase inhibitors induce a BRCAness phenotype that sensitizes NSCLC to PARP inhibitor and ionizing radiation. Proc Natl Acad Sci U S A. 2019;116:22609–18.
Pulliam N, Fang F, Ozes AR, Tang J, Adewuyi A, Keer H, et al. An Effective Epigenetic-PARP Inhibitor Combination Therapy for Breast and Ovarian Cancers Independent of BRCA Mutations. Clin Cancer Res. 2018;24:3163–75.
Modrich P. Mismatch repair, genetic stability, and cancer. Science. 1994;266:1959–60.
Pećina-Šlaus N, Kafka A, Salamon I, Bukovac A. Mismatch Repair Pathway, Genome Stability and Cancer. Front Mol Biosci. 2020;7:122.
Kubota Y, Kawazoe A, Sasaki A, Mishima S, Sawada K, Nakamura Y, et al. The Impact of Molecular Subtype on Efficacy of Chemotherapy and Checkpoint Inhibition in Advanced Gastric Cancer. Clin Cancer Res. 2020;26:3784–90.
Pasanen A, Loukovaara M, Bützow R. Clinicopathological significance of deficient DNA mismatch repair and MLH1 promoter methylation in endometrioid endometrial carcinoma. Mod Pathol. 2020;33:1443–52.
Fraune C, Rosebrock J, Simon R, Hube-Magg C, Makrypidi-Fraune G, Kluth M, et al. High homogeneity of MMR deficiency in ovarian cancer. Gynecol Oncol. 2020;156:669–75.
Xiao X, Dong D, He W, Song L, Wang Q, Yue J, et al. Mismatch repair deficiency is associated with MSI phenotype, increased tumor-infiltrating lymphocytes and PD-L1 expression in immune cells in ovarian cancer. Gynecol Oncol. 2018;149:146–54.
Mishalian I, Bayuh R, Levy L, Zolotarov L, Michaeli J, Fridlender ZG. Tumor-associated neutrophils (TAN) develop pro-tumorigenic properties during tumor progression. Cancer Immunol Immunother. 2013;62:1745–56.
Ghoneum A, Almousa S, Warren B, Abdulfattah AY, Shu J, Abouelfadl H, et al. Exploring the clinical value of tumor microenvironment in platinum-resistant ovarian cancer. Semin Cancer Biol. 2021;S1044–579X(21):00003–1.
Jaillon S, Ponzetta A, Di Mitri D, Santoni A, Bonecchi R, Mantovani A. Neutrophil diversity and plasticity in tumour progression and therapy. Nat Rev Cancer. 2020;20:485–503.
Shaul ME, Fridlender ZG. Tumour-associated neutrophils in patients with cancer. Nat Rev Clin Oncol. 2019;16:601–20.
Chen S, Zhang L, Yan G, Cheng S, Fathy AH, Yan N, et al. Neutrophil-to-Lymphocyte Ratio Is a Potential Prognostic Biomarker in Patients with Ovarian Cancer: A Meta-Analysis. Biomed Res Int. 2017;2017:7943467.
Sato E, Olson SH, Ahn J, Bundy B, Nishikawa H, Qian F, et al. Intraepithelial CD8+ tumor-infiltrating lymphocytes and a high CD8+/regulatory T cell ratio are associated with favorable prognosis in ovarian cancer. Proc Natl Acad Sci U S A. 2005;102:18538–43.
This work was supported by the Key R&D Guidance Plan Project in Liaoning Province (2019JH8/10300022), and the Beijing Kanghua Foundation for the Development of Traditional Chinese and Western Medicine Gynecological Oncology Special Research Fund (KH-2021-LLZX-010). The funding bodies did not play a role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Ethics approval and consent to participate
The study was approved by the Research Ethics Board at Shengjing Hospital of China Medical University, and an informed consent exemption was approved. All methods were carried out in accordance with relevant guidelines and regulations.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Gou, R., Zheng, M., Hu, Y. et al. Identification and clinical validation of NUSAP1 as a novel prognostic biomarker in ovarian cancer. BMC Cancer 22, 690 (2022). https://doi.org/10.1186/s12885-022-09753-4
- Ovarian cancer
- Prognostic biomarker
- Immune infiltration