- Open Access
Evaluate the diagnostic and prognostic value of NUSAP1 in papillary thyroid carcinoma and identify the relationship with genes, proteins, and immune factors
World Journal of Surgical Oncology volume 20, Article number: 207 (2022)
Nucleolar spindle-associated protein 1 (NUSAP1) is reported to be a useful diagnostic and prognostic marker for a variety of cancers, but relevant studies are lacking in papillary thyroid carcinoma (PTC).
The relationship between NUSAP1 expression and the overall survival (OS) of pan-cancer was examined by GEPIA and KMplot. We explored the relationship between NUSAP1 and clinical PTC data based on the THCA dataset of TCGA and the GEO dataset of NCBI; GO, KEGG analysis, and ceRNA networks were performed on co-expressed genes through LinkedOmics and Starbase. We assessed the relevance between NUSAP1 and the tumor microenvironment using ESTIMATE, correlations between NUSAP1 and immune cells with TIMER, the relationship between NUSAP1 and immunotherapy by TCIA, and small-molecule drugs targeting NUSAP1 that can be discovered using the CMap database.
Higher expression of NUSAP1 in pan-cancer tissues was correlated with shorter OS. NUSAP1 was also significantly expressed in PTC tissues and was an independent prognostic risk factor. Compared to the NUSAP1 low expression group, the NUSAP1 high expression group was more likely to also have lymph node metastasis, pathological PTC type, shorter progression-free survival (PFS), and higher scores for immune checkpoint inhibitor treatment. The genes associated with NUSAP1 were mostly involved in the cell cycle, immune-related pathways, and AITD. Ten lncRNAs (GAS5, SNHG7, UCA1, SNHG1, HCP5, DLEU2, HOTAIR, TP53TG1, SNHG12, C9orf106), eleven miRNAs (hsa-miR-10a-5p, hsa-miR-10b-5p, hsa-miR-18a-5p, hsa-miR-18b-5p, hsa-miR-128-3p, hsa-miR-214-3p, hsa-miR-219a-2-3p, hsa-miR-339-5p, hsa-miR-494-3p, hsa-miR-545-3p, hsa-miR-769-5p), and one mRNA (NUSAP1) were constructed. NUSAP1 participated in the formation of the tumor microenvironment. CMap predicted the 10 most important small molecules about NUSAP1.
In PTC, NUSAP1 shows good diagnostic value and prognostic value; NUSAP1 impacts the cell cycle, immune-related pathways, and AITD and has a complex effect on the tumor microenvironment in PTC.
Thyroid cancer is a common endocrine malignancy. According to the World Health Organization, papillary thyroid carcinoma (PTC) is the most common histological type of thyroid cancer . The incidence of thyroid cancer has increased in recent years, with the most prominent changes occurring in the incidence of PTC . Previous studies  have indicated that the recurrence rate of PTC reaches 28% and is related to age, tumor type, and disease stage. Relapse is more likely with increased age, follicular papillary thyroid carcinoma (FVPTC), T4 stage, and cervical lymph node involvement with distant metastasis. Therefore, it is necessary to study and identify high-risk patients in order to ensure appropriate treatment . PTC has been shown to be more prevalent among people with autoimmune diseases. In a prospective cohort study of 9851 subjects, Morais et al.  demonstrated that the risk of malignant nodules increased significantly when the patient had autoimmune thyroiditis (AIT) as compared to those who did not have the disease. Briseis et al.  conducted a prospective study of about 50,000 patients with follow-up for 10 years. Compared with people without diabetes, those with diabetes had a significantly higher risk of developing thyroid cancer. Furthermore, diabetes has been linked to differentiated thyroid cancer. Although previous studies have demonstrated that autoimmune diseases increase the risk of PTC, there are few molecular mechanisms to explain this. Therefore, it is necessary to study the mechanism of PTC in-depth to better understand this condition.
Nucleolar spindle-associated protein 1  (NUSAP1) is a protein-related to cell proliferation. This protein is selectively expressed during cell proliferation and reaches a peak at the transition from interphase to mitosis and is then degraded during cell division. It was first reported in 2007 that expression of NUSAP1 was increased in melanoma compared to normal tissue . A number of recent studies have shown the involvement of NUSAP1 in cancer and have explored its mechanism of action in various diseases, including glioblastoma , bladder cancer , nasopharyngeal carcinoma , non-small-cell lung cancer , liver cancer , osteosarcoma , ovarian cancer , breast cancer , and gastric cancer . The report by Guo et al.  demonstrated that NUSAP1 induces the activation of the Hippo signaling pathway by promoting the stability of Yes1-associated transcriptional regulator (YAP1) protein, thereby promoting the growth and invasion of gastric cancer cells and tumor growth. It has also been reported that silencing NUSAP1 reduces the Wnt/β-catenin signaling pathway in nasopharyngeal cancer by downregulating phosphorylation of glycogen synthase kinase 3 beta (GSK-3β), reducing cell proliferation, and reducing invasion . In 2011, Salmena et al . proposed that competing endogenous RNAs (ceRNAs) can regulate gene expression by binding to the same miRNA. A growing number of researchers demonstrated that ceRNAs were involved in the development of various tumors. However, few studies have been conducted on the role of the NUSAP1-related ceRNA regulatory network in PTC.
Immunotherapy has brought a huge revolution to the treatment of cancer. Recent developments in immune checkpoint blockade therapy (ICB) have resulted in increased success for patients with CD274 (also known as PDL1), PDCD1 (also known as PD1), and CTLA4 . Mehnert et al.  found through a clinical trial on PTC (NCT02054806) that pembrolizumab (PD-1 inhibitor) not only exerts anti-tumor activity, but also has controllable safety. French  concluded that a large amount of data is currently accumulating on the use of single-agent ICB therapy and ICB combination therapy to treat thyroid cancer. Therefore, it is necessary to find some new small-molecule drugs.
Extensive research has been conducted on NUSAP1, but the research value of NUSAP1 in PTC remains unclear. Therefore, in the current study, we used bioinformatics to conduct functions and therapy studies related to NUSAP1 in PTC to lay a theoretical foundation for subsequent research.
Materials and methods
We downloaded The Cancer Genome Atlas Thyroid Cancer (TCGA THCA) dataset, including miRNA expression, clinical data, and gene expression, from the University of California at Santa Cruz (UCSC) Xena browserFootnote 1. The dataset included 59 normal thyroid tissues and 497 PTC tissues.
For evaluating the diagnostic value of NUSAP1, we downloaded data from the National Center for Biotechnology Information (NCBI) Gene Expression Omnibus (GEO) databaseFootnote 2. The download included normal thyroid tissues and PTC tissues. The dataset included GSE27155 (N = 4, T = 51), GSE33630 (N = 45, T = 49), GSE58545 (N = 18, T = 27), and GSE60542 (N = 30, T = 33).
The expression of NUSAP1 in pan-carcinoma
The Tumor Immune Estimation Resource (TIMER) website was used to analyze the expression of NUSAP1 in other cancers and its correlation with the abundance of six types of tumor immune cells. The Gene Expression Profiling Interactive Analysis (GEPIA)Footnote 3  includes RNAseq data and clinical data from TCGA and the Genotype-Tissue Expression project (GTEx), which can be used to analyze the relationship between NUSAP1 expression and overall survival (OS) in other cancers. The KMplot database was used to validate the prognosis of NUSAP1 in other cancers .
The expression of NUSAP1 in different datasets of PTC
First, NUSAP1 mRNA levels were compared in normal thyroid tissue and PTC tissue through the THCA dataset. The diagnostic value of NUSAP1 was evaluated by receiver operating characteristic (ROC). The GSE27155, GSE33630, GSE58545, and GSE60542 datasets were used to further analyze the mRNA level and diagnostic value of NUSAP1. According to the upper quartile of NUSAP1 in the THCA dataset, samples were categorized into the L-NUSAP1 group or H-NUSAP1 group. Samples below the upper quartile were placed in the L-NUSAP1 group, and samples higher than the upper quartile were defined as the H-NUSAP1 group. Kaplan-Meier (K-M) survival curves were evaluated by the log-rank test. The endpoint of this study was progression-free survival (PFS). A heat map was used to visualize the differences in the clinical data between the L-NUSAP1 and H-NUSAP1 groups.
Analysis of the underlying mechanism of NUSAP1
The LinkedOmics websiteFootnote 4 is a multi-omics database that includes 32 types of cancer, with data from TCGA and the Clinical Proteomics Tumor Analysis Consortium (CPTAC) . We used this website to analyze the molecules that were co-expressed with NUSAP1 and visualized expression using volcano maps and heat maps. According to the Pearson correlation coefficient, the first 1000 genes positively correlated with NUSAP1 and the first 1000 genes negatively correlated with NUSAP1 were selected; a PPI network was constructed using StringFootnote 5 (version 11.0; interaction value = 0.99). A Cytoscape was used to visualize (version 3.8.0) and display the core genes at the same time. Finally, GO analysis and KEGG analysis were performed by gene set enrichment analysis (GSEA).
Analysis of NUSAP1 and tumor microenvironment
ESTIMATE  was based on the ssGSEA algorithm to score the matrix and immune gene sets. Immune cell infiltration (immuneScore), stromal content (stromalScore), and tumor purity were calculated for each THCA sample through the “estimate” R package. The TIMERFootnote 6 website  is a tumor-infiltrating immune cell component analysis website for analyzing the correlation between the expression of NUSAP1 and the abundance of six tumor-infiltrating immune cells. The p-values and correlation values (cor) were obtained using Spearman’s rank correlation test after purity adjustment. Scatter plots were used to display NUSAP1 and immune cell biomarkers [B cell (CD19 and CD79A), CD8+ T cell (CD8A and CD8B), CD4+ T cell (CD4), M1 macrophage (NOS2, IRF5, and PTGS2), and M2 macrophage (CD163, VSIG4, and MS4A4A), neutrophil (CEACAM8, ITGAM, and CCR7), and dendritic cell (HLA-DPB1, HLA-DQB1, HLA-DRA, HLA-DPA1, CD1C, NRP1, and ITGAX)].
Construction of the ceRNA network
StarBase v2. 0 searches for micorRNA targets through high-throughput CLIP-Seq experimental data and degradome experimental data was used to predict mRNAs. Intersecting genes predicted to bind to NUSAP1 and differentially expressed in PTC were considered as potential miRNAs and lncRNAs in PTC. Finally, a lncRNA-miRNA-mRNA network was constructed, and the data were visualized using the Cytoscape software 3.8.0.
The value of NUSAP1 in PTC immunotherapy
The correlation between NUSAP1 and immune checkpoints (CD274, CTLA4, and PDCD1) was analyzed through the GEPIA website.
TCIAFootnote 7  contains medical images of many tumors along with clinical information. The difference between PD1 and CTLA4 treatments in the L-NUSAP1 group and the H-NUSAP1 group was analyzed and visualized through the TCIA website. The “RColorBrewer” R package was used to analyze the differences between the L-NUSAP1 and H-NUSAP1 groups in different immunotypes (C1, C2, C3, C4, C5, and C6).
Small-molecule drugs identified by connectivity map
Genes related to NUSAP1 were identified by LinkedOmics. CMap  query was performed on inputting the top 150 positive and top 150 negatively related genes to find small molecules with the same effect as NUSAP1. Negatively related drugs that reverse the effect of NUSAP1 were found in the top 10.
Differences in NUSAP1 expression in normal thyroid tissue and PTC tissue were analyzed using the Mann-Whitney test. Survival analysis was performed by the log-rank test. The difference in clinical information between the L-NUSAP1 and H-NUSAP1 groups was analyzed using the chi-squared test, and the correlation was analyzed using the Spearman correlation. Data were processed through R (v4.2.2) and SPSS version 25.0. Data were visualized with R (v4.2.2) and GraphPad Prism V.9.0 software.
The expression of NUSAP1 and pan-cancer
In order to explore the expression of NUSAP1 in cancer, we analyzed its expression level, prognostic value, and relationship with the six main immune cells. First, we analyzed the expression of NUSAP1 in pan-cancers and visualized the expression of NUSAP1 in 32 cancers and 2 cancer subtypes through box plots. As shown in Fig. 1a, the expression level of NUSAP1 was increased significantly in 16 cancers, highlighted by a gray background color. Metastatic skin melanoma cancer was more common than skin melanoma cancer. But there was no significant change in 16 cancers, so only the expression of NUSAP1 in cancer tissues is shown. For high-expressing cancers, the expression in different subtypes is further analyzed. Among them, the HPV-positive group of head and neck squamous cell carcinoma is higher than the HPV-negative group, and the difference in the three subtypes of breast cancer is not very significant. The prognostic value of NUSAP1 in these 16 cancers was analyzed by log-rank test. The results showed that NUSAP1 was related to five cancers; prostate adenocarcinoma (PRAD), kidney chromophobe (KICH), kidney renal papillary cell carcinoma (KIRP), liver hepatocellular carcinoma (LIHC), and lung adenocarcinoma (LUAD) were associated with shorter OS (Fig. 1b–f). Out of 9499 cases of pan-cancer, cancer risks in the high-expressing NUSAP1 group were 1.7 times that of the low-expressing NUSAP1 group (p < 0.0001) (Fig. 1g). Immune cells in LIHC were most closely associated with NUSAP1, with B cells, DC cells, and macrophages being the most obvious (Fig. 1h). To further investigate the prognostic potential of NUSAP1 in cancer, we used the KMplot database to study NUSAP1 in cancer (Additional file 1: Fig. S1). We found that high NUSAP1 expression was also closely related to the poorer prognosis of OS in KIRC, KIRP, LIHC, and LUAD. Besides, pancreatic ductal adenocarcinoma (PDAC) had the same result.
The expression and diagnostic value of NUSAP1 in PTC
The THCA dataset was used to examine the NUSAP1 mRNA levels in normal thyroid tissue and PTC. As illustrated in Fig. 2a, the expression of NUSAP1 was higher in PTC compared with normal thyroid tissue (p < 0.0001). To determine the diagnostic value, we plotted the ROC curve and calculated its area as 0.7221 (p < 0.0001), indicating that it had a good diagnostic value (Fig. 2b). For the purpose of excluding the influence of individual differences, the THCA dataset was utilized to conduct a paired analysis. The results showed that NUSAP1 expression was increased in 59 paired PTC samples (Fig. 2c, p < 0.0001]. In order to verify the above results, we download the GSE27155, GSE33630, GSE58545, and GSE60542 datasets from the GEO database. As shown in Fig. 2d in the GSE27155 dataset, the expression of NUSAP1 in PTC was significantly higher than that in normal thyroid tissue (p < 0.0076), which was consistent with the results from other datasets [GSE33630 (p < 0.0001), GSE58545 (p <0.0012), GSE60542 (p < 0.0001)] (Fig. 2e–g). The areas under the ROC curve were 0.8824 (p < 0.0115), 0.8063 (p < 0.0001), 0.7809 (p < 0.0018), and 0.7758 (p < 0.0002); all findings support the above results (Fig. 2h–k).
The prognostic value and clinical value of NUSAP1
In order to detect the prognostic value of NUSAP1 in PTC, the THCA dataset was got from the TCGA database and divided into L-NUSAP1 (n = 373) and H-NUSAP1 (n = 124), according to the NUSAP1 expression level. KM survival analysis showed that the progression-free survival (PFS) of the H-NUSAP1 group was shorter (p < 0.0001) compared with the L-NUSAP1 group (Fig. 3a). By regression analysis, we determined that age (p < 0.008), tumor stage (p < 0.001), tumor mutation burden (TMB, p < 0.004), and NUSAP1 expression (p < 0.005) were different from PFS (Table 1). Multivariate Cox regression analysis revealed that NUSAP1 may be an independent prognostic risk factor (HR = 1.5328, 95%CI 1.1171–2.1031, p = 0.0081) after adjusting for gender, age, stage, and TMB (Fig. 3b). In Cox regression, for dichotomous variables, we used males and < 55 years as the reference group, respectively. To determine the clinical value of NUSAP1 in PTC, the chi-squared test was applied to explore the diversity in various clinical data between the L-NUSAP1 and H-NUSAP1 groups (Table 2). The results showed that the H-NUSAP1 group was more prone to lymph node metastasis (p < 0.013) and more typical PTC types (p < 0.05) (Fig. 3c).
Analysis of related genes and functions of NUSAP1
Next, we aimed to study the function of NUSAP1 and related genes. A total of 501 cases of PTC in the THCA dataset were analyzed through the LinkedOmics website. The results of the study showed that a total of 4281 genes were negatively correlated with NUSAP1 expression, and 6195 genes were positively correlated (p < 0.01); these data are visualized as a volcano graph (Fig. 4a). Heat maps were used to visualize the top 50 positively correlated and top 50 negatively correlated genes that were most closely connected to NUSAP1 (Fig. 4b, c). In order to analyze the interaction with NUSAP1-related proteins, 2000 proteins closely linked to NUSAP1 were visualized through String and Cytoscape, and the first 30 core proteins were displayed using R language, including MCM4, CDC20, PLK1, CDK1, MCM6, CDC45, CDK2 MCM5, BUB1, CDC6, and (Fig. 5a, b). We next analyzed the role of NUSAP1 in PTC; GO and KEGG analysis were performed by GSEA. The results showed a positive correlation with biological processes (BP) (Fig. 6a) related to mitotic cell cycle phase transition and adaptive immune response; a negative correlation was found for BP related to protein transmembrane transport and mitochondrial respiratory chain complex assembly. For the cell component (CC) (Fig. 6b), the chromosomal region was positively correlated; negatively related BP included motility and mitochondrial protein complex. For molecular function (MF) (Fig. 6c), BP primarily associated with nucleic acid receptor activity and cyclin-dependent protein kinase activity were positively correlated. For KEGG (Fig. 6d), positive correlations were identified for cytokine-cytokine receptor interaction, cell cycle, type I diabetes mellitus, NF-kappa B signaling pathway, p53 signaling pathway, DNA replication, autoimmune thyroid diseases, and cell adhesion molecules (CAMs); a negative correlation was identified for oxidative phosphorylation. To learn more about the NUSAP1 function in the H-NUSAP1 group, GO and KEGG analysis were performed by GSEA (Fig. 6e). The most likely related KEGG pathways (Fig. 6f) included “CELL_CYCLE” and “P53_SIGNALING_PATHWAY.”
The relationship between NUSAP1 and tumor microenvironment
To investigate the association between NUSAP1 and the tumor microenvironment, we first calculated stromalScore, immuneScore, and tumor purity in PTC tissues using the ESTIMATE method. The stromalScore and immuneScore were much higher in the H-NUSAP1 group than in the L-NUSAP1 group (both p-values were < 0.0001) (Fig. 7a, b). Moreover, the tumor purity level in the L-NUSAP1 group was substantially lower than that in the H-NUSAP1 group (Fig. 7c). We used a histogram to see the proportion of immune cells in PTC tissue (Additional file 1: Fig. S2). In order to analyze the relationship between NUSAP1 and immune cells, the correlation between NUSAP1 and B cells, dendritic cells, CD4+ T cells, CD8+ T cells, macrophages, and neutrophils in PTC was analyzed by TIMER. The correlation coefficients of NUSAP1 and the abovementioned immune cells, respectively, were 0.508, 0.608, 0.35, 0.24, 0.517, and 0.352 (p < 0.0001, Fig. 7d). Heat maps were used to show the interaction between different immune cells in PTC (Additional file 1: Fig. S3). To further analyze the relationship between NUSAP1 and immune cells, we first needed to better understand the changes in immune cell molecular markers. The results (Fig. 7e) revealed positive associations for all, including B cell molecular markers (CD19 and CD79A), CD8+ T cell molecular markers (CD8A, CD8B), CD4+ T cell molecular markers (CD4), M1 macrophage molecular markers (IRF5 and PTGS2), M2 macrophage molecular markers (CD163, VSIG4, and MS4A4A), neutrophil molecular markers (CEACAM8, ITGAM, and CCR7), and dendritic cell molecular markers (HLA-DPB1, HLA-DQB1, HLA-DRA, HLA-DPA1, CD1C, and ITGAX). NUSAP1 was not correlated with the M1 macrophage molecular marker NOS2 or the dendritic cell molecular marker NRP1.
Construction of ceRNA network
To better explore the different roles in NUSAP1-related miRNAs and lncRNAs of papillary thyroid carcinoma, we constructed mRNA–miRNA–lncRNA networks. By exploring the miRNAs that bind to NUSAP1 through the starbase website, we found that after removing duplicate values, 24 miRNAs had binding sites for NUSAP1 (Fig. 8a). Since miRNAs negatively regulate NUSAP1, we performed the miRNA expression quantification data in the TCGA to search for miRNAs that are underexpressed in PTC. Eleven miRNAs were identified. To further find related lncRNAs, we analyzed the lncRNAs bound to miRNAs again by starbase and ensured that lncRNAs were highly expressed in PTC then found GAS5, SNHG7, UCA1, SNHG1, HCP5, DLEU2, HOTAIR, TP53TG1, SNHG12, and C9orf106. Based on the above findings, we constructed a ceRNA network by cytoscape (Fig. 8b).
The value of NUSAP1 in therapy
The GEPIA website was utilized to evaluate the association between NUSAP1 and immune checkpoints (CD274, CTLA4, and PDCD1) in order to further study the effect of NUSAP1 on PTC immunotherapy. The results showed that the R value of NUSAP1 and CD274 was 0.36 (Fig. 9a, p < 0.0001), the R value of CTLA4 was 0.43 (Fig. 9b, p < 0.0001), and the R value of PDCD1 was 0.26 (Fig. 9c, p < 0.0001). On this premise, the treatment scores of the immune checkpoints were examined, and there were no statistical differences between the L-NUSAP1 group and the H-NUSAP1 group without CTLA4 and PD1 treatment (Fig. 9d). PD1 treatment alone (Fig. 9e) and CTLA4 treatment alone (Fig. 9f) showed higher immune scores in the H-NUSAP1 group than in the L-NUSAP1 group, and the difference was statistically significant (p < 0.0001). With combined PD1 and CTLA4 treatment (Fig. 9g), the immunotherapy score in the H-NUSAP1 group still had a higher immune score (p < 0.0001), but there were almost no statistically significant changes with respect to immune phenotyping (Fig. 9h).
To explore drugs related to NUSAP1, genes significantly associated with NUSAP1 were entered into the Cmap database to find the gene expression profiles of small molecules matching these genes. The top 10 related drugs/molecules are shown in Table 3 and sorted by score. The top negative enrichment scores belong to CDK inhibitors, including palbociclib and purvalanol-a.
The rapid development of sequencing technology has promoted our understanding of pan-cancer, especially for thyroid cancer [31, 32]. In this analysis, we first identified the differential expression of NUSAP1 in pan-cancer using TIMER and validated it in KMplot. High NUSAP1 expression is associated with poor disease prognosis, and NUSAP1 is interrelated to immune cell abundance. Using public databases, we found that the level of NUSAP1 expression in PTC tissues was significantly higher than in normal tissue, showed good diagnostic value, was closely related to the prognosis of PTC patients, and was an independent risk factor for PTC. Therefore, we completed an in-depth study of the potential functions of NUSAP1 and discovered that it was related to cell proliferation and immune response. The analysis of NUSAP1 and the tumor microenvironment lays the foundation, and the analysis results were consistent with expectations. In the H-NUSAP1group, stromal and immune cell infiltration is likely to be greater, and tumor purity is possible to be less. We also constructed a lncRNA–miRNA–mRNA network. Considering that NUSAP1 is closely relevant to immune infiltration, it is speculated that NUSAP1 is related to the efficacy of immunosuppressive agents. We used TCIA data to analyze that PD1 immunotherapy and CTLA4 immunotherapy may be better in the H-NUSAP1group and screened out several potential small molecule drugs that can reverse NUSAP1.
NUSAP1 is highly abundant in many cancers, is related to a short OS period, and has a significant positive correlation with six kinds of immune cells. Therefore, it is conjectured that NUSAP1 also performs certain biological functions in thyroid cancer. Li et al.  found that the expression of NUSAP1 was notably increased in anaplastic thyroid carcinoma (ATC). Therefore, we used the THCA dataset to analyze the expression level of NUSAP1 in PTC. The results showed that NUSAP1 was also increased significantly in PTC tissues and in four GEO databases. The same results are seen in all datasets. In recent years, studies have reported that NUSAP1 was involved in malignant pleural mesothelioma (MPM) , liver cancer , prostate cancer , and other cancers, with diagnostic value. We used the THCA dataset and analyzed the GSE27155, GSE33630, GSE58545, and GSE60542 datasets to demonstrate that NUSAP1 has a better diagnostic value in PTC. In PTC, patients with high NUSAP1 levels have shorter PFS and are prone to relapse. In NSCLC , LIHC , and other cancers, NUSAP1 is associated with a poor prognosis of the disease. We demonstrated that in PTC, NUSAP1 is an independent risk factor. It was reported that central lymph node metastasis (CLNM) was a risk factor for poor prognosis and recurrence of papillary thyroid microcarcinoma . In addition, patients with high NUSAP1 expression are more likely to develop lymph node metastasis, and these patients primarily have typical pathological types of PTC, which further illustrates the potential of NUSAP1 to serve as a prognostic marker.
We found that NUSAP1 interacts with cell cycle-related proteins, such as MCM4, CDC20, and PLK1. In addition, we used GSEA analysis to determine that NUSAP1 may affect the cell cycle, DNA replication, cell adhesion, the NF-KB signaling pathway, and autoimmune thyroid disease (AITD); high expression of NUSAP1 does indeed participate in the cell cycle pathway. In addition, it was also discovered that high expression of NUSAP1 is linked to ubiquitination-mediated protein degradation. A recent protein genomics study on PTC reported that the E3 ubiquitin ligase HUWE1 and DUB were upregulated in tumors and metastatic disease . According to reports, silencing NUSAP1 inhibited proliferation, migration, and invasion of bladder cancer cells . Zhao et al.  found that NUSAP1 was involved in GBM cell proliferation, apoptosis, and DNA damage in patients with malignant glioma (GBM) and was related to GBM chemotherapy resistance. Autoimmunity and inflammation are risk factors for PTC. Zhou et al. and Li et al. [42, 43] reported that HT was a protective factor for PTC, and high TPOAb concentrations were less prone to CLNM. However, it has also been suggested that high urinary iodine, as a risk factor for HT, promotes CLNM of PTC . At present, the exact mechanism of the association between AITD and PTC remains unclear . We observed that NUSAP1 participates in AITD, which provides a theoretical basis for exploring the relationship between AITD and PTC.
As we all know, immune cells that infiltrate a tumor affect the tumor microenvironment. The tumor microenvironment has a complex impact on the cancer cells with respect to invasion, proliferation, and migration. In liver cancer, Chen et al. and Huang et al. [38, 46] reported that NUSAP1 expression was closely related to immune cells. Our study also found that NUSAP1 was positively correlated with dendritic cells, B cells, macrophages, CD8+ T cells, CD4+ T cells, and neutrophils. However, high expression of NUSAP1 is related to the poor prognosis of PTC, and it is speculated that NUSAP1 may be related to immune escape. This indicates that NUSAP1 plays a complex change in tumor immunity.
mRNAs act as a sponge for miRNAs and are regulated by miRNAs . Starbase is the most common tool for predicting mRNA-binding miRNAs. We revealed a total of 24 miRNAs with binding sites to NUSAP1, and hsa-miR-10a-5p, hsa-miR-10b-5p, hsa-miR-18a-5p, hsa-miR-18b-5p, hsa-miR-128-3p, hsa-miR-214-3p, hsa-miR-219a-2-3p, hsa-miR-339-5p, hsa-miR-494-3p, hsa-miR-545-3p, and hsa-miR-769-5p were lowly expressed in PTC. Chen et al.  verified the binding site of hsa-miR-769-5p and NUSAP1 by luciferase gene reporter gene detection, and NUSAP1 can reverse the cell growth, migration, invasion, and apoptosis caused by overexpression of hsa-miR-769-5p in bladder cancer. NUSAP1 is also regulated by hsa-miR-18b-5p and affects the proliferation of hepatocellular carcinoma . We constructed a ceRNA network by finding lncRNAs that bind to miRNAs and identified a total of 10 lncRNAs. It is found that lncRNAs can bind to multiple miRNAs, among which GAS5 binds the most miRNAs. GAS5, SNHG7, UCA1, SNHG1, HCP5, DLEU2, HOTAIR, and SNHG12 were reported to be differentially expressed in PTC and involved in malignant behaviors [49,50,51,52,53,54,55,56].
Next, we analyzed the interaction between NUSAP1 and three common immune checkpoint molecules and revealed that NUSAP1 was mainly related to CD274 and CTLA4, with a low correlation to PDCD1. However, the expression level of NUSAP1 is related to PD1 immunotherapy scores and CTLA4 immunotherapy scores. The treatment scores were all positively correlated, and it was also positively correlated with PD1 and CTLA4 combined treatment. Anti-cancer immunotherapy , especially immune checkpoint inhibitors, can release the immune system and activate cytotoxic lymphocytes to kill cancer cells. Consistent with the research of Thorsson et al. , our research also observed that PTC was mainly C3 (inflammatory type), but we found that high expression of NUSAP1 was not different from the immunophenotyping of PTC. Based on the above results, we further predicted small-molecule drugs that might reverse the effect of NUSAP1, and finally discovered 10 small-molecule drugs. Similar to the above results, CDK inhibitors played an important role in reversing the effects of NUSAP1. In this class of drugs, two small-molecule drugs have been screened, including palbociclib and purvalanol-a. The antitumor effect of palbociclib had been widely reported. In addition, selumetinib increased radioiodine affinity in radioactive iodine (RAI)-refractory tumors, which exerted its effectiveness in advanced thyroid cancer .
Our analysis shows that NUSAP1 can be used as a diagnostic marker of PTC and can affect prognosis. Overexpression of NUSAP1 is associated with lymph node metastasis and recurrence of PTC. NUSAP1 may have an effect on tumor cell growth by affecting the cell cycle and DNA damage repair. Altering cell adhesion impacts the migration and invasion capabilities of tumor cells. This impacts PTC through altering AITD and has a wide range of prognostic values in pan-cancer. We created a NUSAP1-specific ceRNA network and predicted NUSAP1-related small-molecule drugs. However, the role of NUSAP1 in the tumor microenvironment is complex. It is positively correlated with a variety of immune cells and is also positively correlated with immune checkpoint molecules. It is speculated that NUSAP1 is related to immune escape. Therefore, further studies are required to determine the role of NUSAP1 in PTC to provide a basis for the user of NUSAP1 as a therapeutic target of PTC.
Through a multifaceted analysis of the value and function of NUSAP1 in PTC, our study predicts that NUSAP1 may become a potential diagnostic and prognostic marker of PTC, as well as a predictive molecule for the efficacy of therapy. NUSAP1 is related to the development and recurrence of PTC and also provides a link between AITD and PTC. At the same time, NUSAP1 may be involved in immune escape and affect the tumor microenvironment, so immune checkpoint inhibitor therapy may have a better effect. However, further experiments are required to verify the above results.
Availability of data and materials
The data used to support this study are available publicly in the University of California at Santa Cruz (UCSC) database and Gene Expression Omnibus (GEO).
Nucleolar spindle-associated protein 1
Papillary thyroid carcinoma
Follicular papillary thyroid carcinoma
Yes1-associated transcriptional regulator
Glycogen synthase kinase 3 beta
Immune checkpoint blockade
- TCGA THCA:
The Cancer Genome Atlas Thyroid Cancer
University of California at Santa Cruz
National Center for Biotechnology Information
Gene Expression Omnibus
Tumor Immune Estimation Resource
Gene Expression Profiling Interactive Analysis
Genotype-Tissue Expression Project
Receiver operating characteristic
Clinical Proteomics Tumor Analysis Consortium
Gene set enrichment analysis
Kidney renal papillary cell carcinoma
Liver hepatocellular carcinoma
Tumor mutation burden
Cell adhesion molecules
Anaplastic thyroid carcinoma
Malignant pleural mesothelioma
Autoimmune thyroid disease
Fagin JA, Wells SA Jr. Biologic and clinical perspectives on thyroid cancer. N Engl J Med. 2016;375(11):1054–67. https://doi.org/10.1056/NEJMra1501993.
Roman BR, Morris LG, Davies L. The thyroid cancer epidemic, 2017 perspective. Curr Opin Endocrinol Diabetes Obes. 2017;24(5):332–6. https://doi.org/10.1097/MED.0000000000000359.
Grogan RH, Kaplan SP, Cao H, Weiss RE, Degroot LJ, Simon CA, et al. A study of recurrence and death from papillary thyroid cancer with 27 years of median follow-up. Surgery. 2013;154(6):1436–46; discussion 46–7. https://doi.org/10.1016/j.surg.2013.07.008.
Raue F, Frank-Raue K. Thyroid cancer: risk-stratified management and individualized therapy. Clin Cancer Res. 2016;22(20):5012–21. https://doi.org/10.1158/1078-0432.CCR-16-0484.
Silva de Morais N, Stuart J, Guan H, Wang Z, Cibas ES, Frates MC, et al. The impact of Hashimoto thyroiditis on thyroid nodule cytology and risk of thyroid cancer. J Endocr Soc. 2019;3(4):791–800. https://doi.org/10.1210/js.2018-00427.
Aschebrook-Kilfoy B, Sabra MM, Brenner A, Moore SC, Ron E, Schatzkin A, et al. Diabetes and thyroid cancer risk in the National Institutes of Health-AARP Diet and Health Study. Thyroid. 2011;21(9):957–63. https://doi.org/10.1089/thy.2010.0396.
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(6):1017–29. https://doi.org/10.1083/jcb.200302129.
Ryu B, Kim DS, Deluca AM, Alani RM. Comprehensive expression profiling of tumor cell lines identifies molecular signatures of melanoma progression. PLoS One. 2007;2(7):e594. https://doi.org/10.1371/journal.pone.0000594.
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(1):44. https://doi.org/10.1038/s41392-020-0137-7.
Chen Y, Zhang W, Kadier A, Zhang H, Yao X. MicroRNA-769-5p suppresses cell growth and migration via targeting NUSAP1 in bladder cancer. J Clin Lab Anal. 2020;34(5):e23193. https://doi.org/10.1002/jcla.23193.
Zhang L, Dang Y, Wang Y, Fan X. Nucleolar and spindle-associated protein 1 accelerates cellular proliferation and invasion in nasopharyngeal carcinoma by potentiating Wnt/β-catenin signaling via modulation of GSK-3β. J Bioenerg Biomembr. 2020;52(6):441–51. https://doi.org/10.1007/s10863-020-09860-6.
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(4):3886–93. https://doi.org/10.1002/jcp.29282.
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(4):2969–78. https://doi.org/10.3892/etm.2019.7314.
Wang H, Liu Z, Wu P, Wang H, Ren W. NUSAP1 accelerates osteosarcoma cell proliferation and cell cycle progression via upregulating CDC20 and cyclin A2. Onco Targets Ther. 2021;14:3443–54. https://doi.org/10.2147/ott.S295818.
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(11):4397–405. https://doi.org/10.1002/jcb.29661.
Sun L, Shi C, Liu S, Zhang E, Yan L, Ji C, et al. Overexpression of NuSAP1 is predictive of an unfavourable prognosis and promotes proliferation and invasion of triple-negative breast cancer cells via the Wnt/β-catenin/EMT signalling axis. Gene. 2020;747:144657. https://doi.org/10.1016/j.gene.2020.144657.
Ge Y, Li Q, Lin L, Jiang M, Shi L, Wang B, et al. Downregulation of NUSAP1 suppresses cell proliferation, migration, and invasion via inhibiting mTORC1 signalling pathway in gastric cancer. Cell Biochem Funct. 2020;38(1):28–37. https://doi.org/10.1002/cbf.3444.
Guo H, Zou J, Zhou L, Zhong M, He Y, Huang S, et al. NUSAP1 promotes gastric cancer tumorigenesis and progression by stabilizing the YAP1 protein. Front. Oncol. 2020;10:591698. https://doi.org/10.3389/fonc.2020.591698.
Salmena L, Poliseno L, Tay Y, Kats L, Pandolfi PP. A ceRNA hypothesis: the Rosetta Stone of a hidden RNA language? Cell. 2011;146(3):353–8. https://doi.org/10.1016/j.cell.2011.07.014.
Yao L, Jia G, Lu L, Bao Y, Ma W. Factors affecting tumor responders and predictive biomarkers of toxicities in cancer patients treated with immune checkpoint inhibitors. Int Immunopharmacol. 2020;85:106628. https://doi.org/10.1016/j.intimp.2020.106628.
Mehnert JM, Varga A, Brose MS, Aggarwal RR, Lin CC, Prawira A, et al. Safety and antitumor activity of the anti-PD-1 antibody pembrolizumab in patients with advanced, PD-L1-positive papillary or follicular thyroid cancer. BMC Cancer. 2019;19(1):196. https://doi.org/10.1186/s12885-019-5380-3.
French JD. Immunotherapy for advanced thyroid cancers - rationale, current advances and future strategies. Nat Rev Endocrinol. 2020;16(11):629–41. https://doi.org/10.1038/s41574-020-0398-9.
Tang Z, Li C, Kang B, Gao G, Li C, Zhang Z. GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 2017;45(W1):W98–W102. https://doi.org/10.1093/nar/gkx247.
Lanczky A, Gyorffy B. Web-based survival analysis tool tailored for medical research (KMplot): development and implementation. J Med Internet Res. 2021;23(7):e27633. https://doi.org/10.2196/27633.
Vasaikar SV, Straub P, Wang J, Zhang B. LinkedOmics: analyzing multi-omics data within and across 32 cancer types. Nucleic Acids Res. 2018;46(D1):D956–D63. https://doi.org/10.1093/nar/gkx1090.
Yoshihara K, Shahmoradgoli M, Martinez E, Vegesna R, Kim H, Torres-Garcia W, et al. Inferring tumour purity and stromal and immune cell admixture from expression data. Nat Commun. 2013;4:2612. https://doi.org/10.1038/ncomms3612.
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(21):e108–e10. https://doi.org/10.1158/0008-5472.CAN-17-0307.
Li JH, Liu S, Zhou H, Qu LH, Yang JH. starBase v2.0: decoding miRNA-ceRNA, miRNA-ncRNA and protein-RNA interaction networks from large-scale CLIP-Seq data. Nucleic Acids Res. 2014;42(Database issue):D92–7. https://doi.org/10.1093/nar/gkt1248.
Charoentong P, Finotello F, Angelova M, Mayer C, Efremova M, Rieder D, et al. Pan-cancer immunogenomic analyses reveal genotype-immunophenotype relationships and predictors of response to checkpoint blockade. Cell Rep. 2017;18(1):248–62. https://doi.org/10.1016/j.celrep.2016.12.019.
Subramanian A, Narayan R, Corsello SM, Peck DD, Natoli TE, Lu X, et al. A next generation connectivity map: L1000 platform and the first 1,000,000 profiles. Cell. 2017;171(6):1437–52 e17. https://doi.org/10.1016/j.cell.2017.10.049.
Li C, Yuan Q, Xu G, Yang Q, Hou J, Zheng L, et al. A seven-autophagy-related gene signature for predicting the prognosis of differentiated thyroid carcinoma. World J Surg Oncol. 2022;20(1):129. https://doi.org/10.1186/s12957-022-02590-6.
He Y, Yu X, Zhang M, Guo W. Pan-cancer analysis of m(5)C regulator genes reveals consistent epigenetic landscape changes in multiple cancers. World J Surg Oncol. 2021;19(1):224. https://doi.org/10.1186/s12957-021-02342-y.
Li L, Zhu M, Huang H, Wu J, Meng D. Identification of hub genes in anaplastic thyroid carcinoma: evidence from bioinformatics analysis. Technol Cancer Res Treat. 2020;19:1533033820962135. https://doi.org/10.1177/1533033820962135.
Lin P, Zhao Y, Li X, Liang Z. KIAA0101 in malignant pleural mesothelioma: a potential diagnostic and prognostic marker. Comb Chem High Throughput Screen. 2021. https://doi.org/10.2174/1386207324666210707105634.
Gao S, Gang J, Yu M, Xin G, Tan H. Computational analysis for identification of early diagnostic biomarkers and prognostic biomarkers of liver cancer based on GEO and TCGA databases and studies on pathways and biological functions affecting the survival time of liver cancer. BMC Cancer. 2021;21(1):791. https://doi.org/10.1186/s12885-021-08520-1.
Feng T, Wei D, Li Q, Yang X, Han Y, Luo Y, et al. Four novel prognostic genes related to prostate cancer identified using co-expression structure network analysis. Front Genet. 2021;12:584164. https://doi.org/10.3389/fgene.2021.584164.
Zhao D, Mu HJ, Shi HB, Bi HX, Jiang YF, Liu GH, et al. Identification of therapeutic targets and mechanisms of tumorigenesis in non-small cell lung cancer using multiple-microarray analysis. Medicine (Baltimore). 2020;99(44):e22815. https://doi.org/10.1097/md.0000000000022815.
Huang R, Liu J, Li H, Zheng L, Jin H, Zhang Y, et al. Identification of hub genes and their correlation with immune infiltration cells in hepatocellular carcinoma based on GEO and TCGA databases. Front Genet. 2021;12:647353. https://doi.org/10.3389/fgene.2021.647353.
Wen X, Jin Q, Cen X, Qiu M, Wu Z. Clinicopathologic predictors of central lymph node metastases in clinical node-negative papillary thyroid microcarcinoma: a systematic review and meta-analysis. World J Surg Oncol. 2022;20(1):106. https://doi.org/10.1186/s12957-022-02573-7.
Krishnan A, Berthelet J, Renaud E, Rosigkeit S, Distler U, Stawiski E, et al. Proteogenomics analysis unveils a TFG-RET gene fusion and druggable targets in papillary thyroid carcinomas. Nat Commun. 2020;11(1):2056. https://doi.org/10.1038/s41467-020-15955-w.
Gao S, Yin H, Tong H, Zhan K, Yang G, Hossain MA, et al. Nucleolar and spindle associated protein 1 (NUSAP1) promotes bladder cancer progression through the TGF-β signaling pathway. Onco Targets Ther. 2020;13:813–25. https://doi.org/10.2147/ott.S237127.
Li X, Zhang H, Zhou Y, Cheng R. Risk factors for central lymph node metastasis in the cervical region in papillary thyroid carcinoma: a retrospective study. World J Surg Oncol. 2021;19(1):138. https://doi.org/10.1186/s12957-021-02247-w.
Zhou SL, Guo YP, Zhang L, Deng T, Xu ZG, Ding C, et al. Predicting factors of central lymph node metastasis and BRAF(V600E) mutation in Chinese population with papillary thyroid carcinoma. World J Surg Oncol. 2021;19(1):211. https://doi.org/10.1186/s12957-021-02326-y.
Zeng Z, Li K, Wang X, Ouyang S, Zhang Z, Liu Z, et al. Low urinary iodine is a protective factor of central lymph node metastasis in papillary thyroid cancer: a cross-sectional study. World J Surg Oncol. 2021;19(1):208. https://doi.org/10.1186/s12957-021-02302-6.
Ferrari SM, Fallahi P, Elia G, Ragusa F, Ruffilli I, Paparo SR, et al. Thyroid autoimmune disorders and cancer. Semin Cancer Biol. 2020;64:135–46. https://doi.org/10.1016/j.semcancer.2019.05.019.
Chen H, Wu J, Lu L, Hu Z, Li X, Huang L, et al. Identification of hub genes associated with immune infiltration and predict prognosis in hepatocellular carcinoma via bioinformatics approaches. Front Genet. 2020;11:575762. https://doi.org/10.3389/fgene.2020.575762.
Hansen TB, Jensen TI, Clausen BH, Bramsen JB, Finsen B, Damgaard CK, et al. Natural RNA circles function as efficient microRNA sponges. Nature. 2013;495(7441):384–8. https://doi.org/10.1038/nature11993.
Yang Z, Li J, Feng G, Wang Y, Yang G, Liu Y, et al. Hepatitis B virus X protein enhances hepatocarcinogenesis by depressing the targeting of NUSAP1 mRNA by miR-18b. Cancer Biol Med. 2019;16(2):276–87. https://doi.org/10.20892/j.issn.2095-3941.2018.0283.
Zuo Y, Qu C, Tian Y, Wen Y, Xia S, Ma M. The HIF-1/SNHG1/miR-199a-3p/TFAM axis explains tumor angiogenesis and metastasis under hypoxic conditions in breast cancer. Biofactors. 2021;47(3):444–60. https://doi.org/10.1002/biof.1702.
Guo L, Lu J, Gao J, Li M, Wang H, Zhan X. The function of SNHG7/miR-449a/ACSL1 axis in thyroid cancer. J Cell Biochem. 2020;121(10):4034–42. https://doi.org/10.1002/jcb.29569.
Liang L, Xu J, Wang M, Xu G, Zhang N, Wang G, et al. LncRNA HCP5 promotes follicular thyroid carcinoma progression via miRNAs sponge. Cell Death Dis. 2018;9(3):372. https://doi.org/10.1038/s41419-018-0382-7.
Liu Y, Li YF, Liu J, Deng ZG, Zeng L, Zhou WB. Long noncoding RNA GAS5 targeting miR-221-3p/cyclin-dependent kinase inhibitor 2B axis regulates follicular thyroid carcinoma cell cycle and proliferation. Pathobiology. 2021;88(4):289–300. https://doi.org/10.1159/000513338.
Yang J, Huang Y, Dong B, Dai Y. Long noncoding RNA DLEU2 drives the malignant behaviors of thyroid cancer through mediating the miR-205-5p/TNFAIP8 axis. Endocr. Connect. 2021;10(4):471–83. https://doi.org/10.1530/ec-21-0046.
Li L, Wang J, Li Z, Qiu S, Cao J, Zhao Y, et al. Diagnostic value of serum lncRNA HOTAIR combined with galectin-3 in benign and papillary thyroid carcinoma. Cancer Manag Res. 2021;13:6517–25. https://doi.org/10.2147/cmar.S312784.
Liu J, Tang X, Lv J, Peng X, Zhang K, Li C, et al. LncRNAs SNHG12 and LINC00152 were associated with progression of patients with papillary thyroid carcinoma. Future Oncol. 2019;15(36):4167–79. https://doi.org/10.2217/fon-2019-0016.
Gao H, Yang JY, Tong LX, Jin H, Liu CZ. Long noncoding RNA UCA1 promotes proliferation and metastasis of thyroid cancer cells by sponging miR-497-3p. Eur Rev Med Pharmacol Sci. 2020;24(14):7555. https://doi.org/10.26355/eurrev_202007_22203.
Ferrari SM, Fallahi P, Galdiero MR, Ruffilli I, Elia G, Ragusa F, et al. Immune and inflammatory cells in thyroid cancer microenvironment. Int J Mol Sci. 2019;20(18). https://doi.org/10.3390/ijms20184413.
Thorsson V, Gibbs DL, Brown SD, Wolf D, Bortone DS, Ou Yang TH, et al. The immune landscape of cancer. Immunity. 2018;48(4):812–30 e14. https://doi.org/10.1016/j.immuni.2018.03.023.
Ho AL, Dedecjus M, Wirth LJ, Tuttle RM, Inabnet WB 3rd, Tennvall J, et al. Selumetinib plus adjuvant radioactive iodine in patients with high-risk differentiated thyroid cancer: a phase III, randomized, placebo-controlled trial (ASTRA). J Clin Oncol. 2022:JCO2100714. https://doi.org/10.1200/JCO.21.00714.
We thank International Science Editing for editing this manuscript.
Ethics approval and consent to participate
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.
Kaplan-Meier survival curves comparing the high and low expression of NUSAP1 in 5 types of cancer. Survival curves of OS in pancreatic ductal adenocarcinoma (PDAC), Kidney Renal Clear Cell Carcinoma (KIRC), Kidney renal papillary cell carcinoma (KIRP), Liver hepatocellular carcinoma (LIHC), Lung adenocarcinoma (LUAD). Figure S2. The proportion of immune cells in the PTC tissue of the THCA data set. The x-axis shows the sample name, and the y-axis shows the proportion of immune cells. Figure S3. Correlation between immune cells in PTC tissue of THCA data set. The content in the lower left corner displays the correlation coefficient, and the upper right corner visualizes the correlation coefficient. Blue represents positive correlation, red represents negative correlation, and the darker the color, the higher the correlation.
About this article
Cite this article
Gao, T., Zhao, L., Zhang, F. et al. Evaluate the diagnostic and prognostic value of NUSAP1 in papillary thyroid carcinoma and identify the relationship with genes, proteins, and immune factors. World J Surg Onc 20, 207 (2022). https://doi.org/10.1186/s12957-022-02652-9
- Papillary thyroid carcinoma
- Prognostic markers
- Tumor immunity