High NUSAP1 expression predicts poor prognosis in colon cancer

Emerging role of RNA modification and long noncoding RNA interaction in cancer

RNA modification, especially N6-methyladenosine, 5-methylcytosine, and N7-methylguanosine methylation, participates in the occurrence and progression of cancer through multiple pathways. The function and expression of these epigenetic regulators have gradually become a hot topic in cancer research. Mutation and regulation of noncoding RNA, especially lncRNA, play a major role in cancer. Generally, lncRNAs exert tumor-suppressive or oncogenic functions and its dysregulation can promote tumor occurrence and metastasis. In this review, we summarize N6-methyladenosine, 5-methylcytosine, and N7-methylguanosine modifications in lncRNAs. Furthermore, we discuss the relationship between epigenetic RNA modification and lncRNA interaction and cancer progression in various cancers. Therefore, this review gives a comprehensive understanding of the mechanisms by which RNA modification affects the progression of various cancers by regulating lncRNAs, which may shed new light on cancer research and provide new insights into cancer therapy.

Identification of potential biomarkers for pancreatic ductal adenocarcinoma: a bioinformatics analysis

PDA is an aggressive cancer with a 5-year survival rate, which is very low. There is no effective prognosis or therapy for PDA because of the lack of target biomarkers. The objective of this article is to identify the target biomarkers for PDA using a bioinformatics approach. In this work, we have analysed the three microarray datasets from the NCBI GEO database. We used the Geo2R tool to analyse the microarray data with the Benjamini and Hochberg false discovery rate method, and the significance level cut-off was set to 0.05. We have identified 659 DEGs from the datasets. There are a total of 15 hub genes that were selected from the PPI network constructed using the STRING application. Furthermore, these 15 genes were evaluated on PDA patients using TCGA and GTEx databases in (GEPIA). The online tool DAVID was used to analyse the functional annotation information for the DEGs. The functional pathway enrichment was performed on the GO and KEGG. The hub genes were mainly enriched for cell division, chromosome segregation, protein binding and microtubule binding. Further, the gene alteration study was performed using the cBioportal tool and screened out six hub genes (ASPM, CENPF, BIRC5, TTK, DLGAP5, and TOP2A) with a high alteration rate in PDA samples. Furthermore, Kaplan-Meier survival analysis was performed on the six hub genes and identified poor-survival outcomes that may be involved in tumorigenesis and PDA development. So, this study concludes that, these six hub genes may be potential prognostic biomarkers for PDA.

Expression Profile and Gene Regulation Network of NUSAP1 in Pan Cancers Based on Integrated Bioinformatics Analysis

Background

Nucleolar and spindle-associated protein 1 (NUSAP1) plays key roles in microtubules and chromosomes in normal cells both structurally and functionally. In malignancies, NUSAP1 is frequently dysregulated and mutated. However, the expression profiles and biological functions of NUSAP1 in tumors remain unclear.

Methods

NUSAP1 expression in BALB/c mice and human normal or tumor tissues was examined using immunohistochemistry. Kaplan-Meier survival analysis was utilized to assess the prognostic significance of NUSAP1 in tumors, and principal component analysis and co-expression analysis were performed to explore the unique roles of NUSAP1. Gene ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses were performed with DAVID. The relevance between NUSAP1 and tumor-infiltrating immune cells was investigated using TIMER. A transcriptional regulation network was constructed using data from The Cancer Genome Atlas.

Results

NUSAP1 expression levels in various mice tissues were different. Compared with normal tissues, NUSAP1 was strongly expressed in several human tumor tissues. We believe that NUSAP1 distinctly impacts the prognosis of several cancers and plays various roles in thymoma and testicular germ cell tumors. Further, NUSAP1 expression levels were significantly positively associated with diverse infiltrating levels of immune cells, including B cells, CD4+ and CD8+ T cells, dendritic cells, and macrophages, in thymoma. The expression level of NUSAP1 demonstrated strong relevance with various immune markers in thymoma. Finally, the miR-1236-5p-NUSAP1 and TCF3-NUSAP1 network revealed the tumor-promoting role of NUSAP1 and pertinent underlying mechanisms in human liver hepatocellular carcinoma.

Conclusion

NUSAP1 may be regarded as a therapeutic target or potential prognostic biomarker for various cancer types.

Differentially expressed discriminative genes and significant meta-hub genes based key genes identification for hepatocellular carcinoma using statistical machine learning

Hepatocellular carcinoma (HCC) is the most common lethal malignancy of the liver worldwide. Thus, it is important to dig the key genes for uncovering the molecular mechanisms and to improve diagnostic and therapeutic options for HCC. This study aimed to encompass a set of statistical and machine learning computational approaches for identifying the key candidate genes for HCC. Three microarray datasets were used in this work, which were downloaded from the Gene Expression Omnibus Database. At first, normalization and differentially expressed genes (DEGs) identification were performed using limma for each dataset. Then, support vector machine (SVM) was implemented to determine the differentially expressed discriminative genes (DEDGs) from DEGs of each dataset and select overlapping DEDGs genes among identified three sets of DEDGs. Enrichment analysis was performed on common DEDGs using DAVID. A protein-protein interaction (PPI) network was constructed using STRING and the central hub genes were identified depending on the degree, maximum neighborhood component (MNC), maximal clique centrality (MCC), centralities of closeness, and betweenness criteria using CytoHubba. Simultaneously, significant modules were selected using MCODE scores and identified their associated genes from the PPI networks. Moreover, metadata were created by listing all hub genes from previous studies and identified significant meta-hub genes whose occurrence frequency was greater than 3 among previous studies. Finally, six key candidate genes (TOP2A, CDC20, ASPM, PRC1, NUSAP1, and UBE2C) were determined by intersecting shared genes among central hub genes, hub module genes, and significant meta-hub genes. Two independent test datasets (GSE76427 and TCGA-LIHC) were utilized to validate these key candidate genes using the area under the curve. Moreover, the prognostic potential of these six key candidate genes was also evaluated on the TCGA-LIHC cohort using survival analysis.

Development of necroptosis-related gene signature to predict the prognosis of colon adenocarcinoma

Colon adenocarcinoma (COAD) is a common malignancy and has a high mortality rate. However, the current tumor node metastasis (TNM) staging system is inadequate for prognostic assessment of COAD patients. Therefore, there is an urgent need to identify reliable biomarkers for the prognosis COAD patients. The aberrant expression of necroptosis-related genes (NRGs) is reported to be associated with tumorigenesis and metastasis. In the present work, we compared the expression profiles of NRGs between COAD patients and normal individuals. Based on seven differentially expressed NRGs, a risk score was defined to predict the prognosis of COAD patients. The validation results from both training and independent external cohorts demonstrated that the risk score is able to distinguish the high and low risk COAD patients with higher accuracies, and is independent of the other clinical factors. To facilitate its clinical use, by integrating the proposed risk score, a nomogram was built to predict the risk of individual COAD patients. The C-index of the nomogram is 0.75, indicating the reliability of the nomogram in predicting survival rates. Furthermore, two candidate drugs, namely dapsone and xanthohumol, were screed out and validated by molecular docking, which hold the potential for the treatment of COAD. These results will provide novel clues for the diagnosis and treatment of COAD.

Identification and clinical validation of NUSAP1 as a novel prognostic biomarker in ovarian cancer

Background

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.

Methods

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.

Results

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.

Conclusions

NUSAP1 may be a valuable prognostic marker, as well as a novel biomarker for evaluating the response to immunotherapy of patients with ovarian cancer.

Transcription factor HOXB2 upregulates NUSAP1 to promote the proliferation, invasion and migration of nephroblastoma cells via the PI3K/Akt signaling pathway

The transcription factor homeobox protein Hox-B2 (HOXB2) and its downstream factor nucleolar and spindle-associated protein 1 (NUSAP1) play important regulatory roles in cell proliferation, invasion and migration. However, their effects and specific mechanisms in nephroblastoma have not been previously investigated, to the best of our knowledge. Therefore, in the present study, the mRNA and protein expression levels of HOXB2 and NUSAP1 were determined in nephroblastoma cells using reverse transcription-quantitative PCR and western blot analyses, respectively. Furthermore, cell transfection experiments were carried out to knock down NUSAP1 and overexpress HOXB2 in nephroblastoma cell lines. The proliferative, invasive and migratory abilities of nephroblastoma cells were assessed by MTT, EdU, colony formation, wound healing and Transwell assays. In addition, the JASPAR website was used to predict the association between HOXB2 and NUSAP1, which was further verified by dual-luciferase reporter and chromatin immunoprecipitation assays. Finally, the expression levels of the PI3K/Akt signaling pathway-related proteins were measured by western blot analysis. The results showed that the expression of NUSAP1 was abnormally upregulated in nephroblastoma cell lines. However, NUSAP1 silencing attenuated the proliferation, invasion and migration abilities of nephroblastoma cells. The results also suggested that HOXB2 could transcriptionally activate NUSAP1. Therefore, HOXB2 overexpression abrogated the inhibitory effect of NUSAP1 silencing on the proliferation and metastasis of nephroblastoma cells, possibly via the PI3K/Akt signaling pathway. The aforementioned findings indicated that HOXB2 may upregulate NUSAP1 to promote the proliferation, invasion and migration of nephroblastoma cells via the PI3K/Akt signaling pathway.

Increased NUSAP1 expression is associated with lymph node metastasis and survival prognosis in bladder urothelial carcinoma

The main route of metastasis of bladder urothelial carcinoma is through lymph nodes; however, its exact mechanism remains unclear. In this study, we found an association of nucleolar and spindle associated protein 1 (NUSAP1) expression with BUC tissues along with lymph node metastasis and the survival prognosis. A total of 178 pathological specimens following radical bladder cancer resection were obtained. NUSAP1 expression was analyzed by immunohistochemistry. We evaluated the correlation between clinicopathological characteristics and NUSAP1 expression. Logistic regression was used to determine the independent variables that influenced lymph node metastasis. Uni- and multi-factorial Cox regression methods were used to determine the prognostic value of NUSAP1 expression in urothelial carcinoma of the bladder. High expression of NUSAP1 in BUC was not significantly related to the patient’s gender, age, or tumor number (p > 0.05), however was significantly associated with pathological grade, tumor diameter, pathological stage, and lymph node metastasis (p < 0.05). Lymph node metastasis was significantly correlated with pathological stage, pathological grade, tumor number, tumor diameter, and NUSAP1 expression (p < 0.05); only NUSAP1 expression was an independent predictor of lymph node metastasis in BUC (OR:1.786, 95% CI 1.229–2.596, p = 0.002). In addition, high NUSAP1 expression was an independent prognostic predictor for BUC. In BUC, NUSAP1 showed high expression and was significantly associated with lymph node metastasis, pathological stage, pathological grade, and tumor diameter. NUSAP1 was an independent predictor of lymph node metastasis and prognosis in BUC; higher expression indicated poorer prognosis of BUC patients.

microRNA-569 inhibits tumor metastasis in pancreatic cancer by directly targeting NUSAP1

MicroRNAs (miRNAs) are known to be involved in the development and progression of pancreatic cancer (PC). In this study, the prognostic significance and mechanistic role of microRNA-569 in PC were explored. Quantitative real-time PCR was used to detect the expression of microRNA-569 in PC tissues and cell lines. Scratch test and Transwell assay were conducted to detect migration and invasion ability. The xenograft nude mice model was used to determine tumor metastasis in vivo. The direct targets of microRNA-569 were determined by using bioinformatics analysis and a dual-luciferase reporter assay. The expression level of microRNA-569 was down-regulated in PC patients with a poor prognosis. In vitro and in vivo experiments indicated that over-expression of microRNA-569 inhibited the migration and invasion of PC cells. MicroRNA-569 negatively regulated NUSAP1 by directly binding its 3'-untranslated region. Further mechanism research implied that the ZEB1 pathway was involved in microRNA-569/NUSAP1 mediation of the biological behaviors in PC. These data demonstrated that microRNA-569 may exert a tumor-suppressing effect in PC and maybe a potential therapeutic target for PC patients.

Analyzing Roles of NUSAP1 From Clinical, Molecular Mechanism and Immune Perspectives in Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is one of the most common carcinomas worldwide. Our study aims to analyze how NUSAP1 affects progression of HCC from clinical, molecular mechanism and immune perspectives. Firstly, we downloaded GSE62232, GSE102079, GSE112790, and GSE121248 gene expression profile datasets from GEO database. R studio was used to screen DEGs of each dataset, and 86 overlapping DEGs of the four datasets were screened at last. Then, CytoHubba plug-in in Cytoscape software was used to screen out NUSAP1 from the 86 DEGs. Subsequently, survival analysis, clinical correlation analysis, independent prognostic analysis, and GSEA enrichment analysis of NUSAP1 were analyzed using HCC patients from GSE76427 dataset, ICGC database, and TCGA database. The results revealed that HCC patients with higher expression level of NUSAP1 had a worse prognosis. NUSAP1 was an independent prognostic factor of HCC, and it may promote HCC progress by regulating cell cycle. To further elucidate its underlying molecular mechanism, we used cBioProtal online data analysis tool to screen all co-expression genes of NUSAP1 and used top 300 co-expression genes to accomplish KEGG and GO enrichment analysis; the results confirmed that NUSAP1 accelerated progression of HCC by regulating cell cycle. We continued to draw KEGG pathway map of cell cycle using co-expression genes enriched in cell cycle pathway by KEGG online tool. The map depicted that most of co-expression genes of NUSAP1 were located in S phase and G2/M phase of the cell cycle, and they could regulate the genes in G1 phase. To further understand the mechanism of cell cycle, we also did qRT-PCR, Western blot, and flow cytometry; the results showed that NUSAP1 was closely associated with CDK4, CDK6, and cyclinD1, which could regulate G1 to S phase transition. Besides, we also analyzed correlation between NUSAP1 and immune cells using HCC patients from GSE76427 dataset, ICGC database, and TCGA database. NUSAP1 was associated with some immune cells, and we speculated that NUSAP1 could also promote HCC progression by influencing T cell CD4 memory resting and macrophage M0 through some underlying mechanism.

Four Novel Prognostic Genes Related to Prostate Cancer Identified Using Co-expression Structure Network Analysis

Prostate cancer (PCa) is one of the most common malignancies for males, but very little is known about its pathogenesis. This study aimed to identify novel biomarkers associated with PCa prognosis and elucidate the underlying molecular mechanism. First, The Cancer Genome Atlas (TCGA) RNA-sequencing data were utilized to identify differentially expressed genes (DEGs) between tumor and normal samples. The DEGs were then applied to construct a co-expression and mined using structure network analysis. The magenta module that was highly related to the Gleason score (r = 0.46, p = 3e-26) and tumor stage (r = 0.38, p = 2e-17) was screened. Subsequently, all genes of the magenta module underwent function annotation. From the key module, CCNA2, CKAP2L, NCAPG, and NUSAP1 were chosen as the four candidate genes. Finally, internal (TCGA) and external data sets (GSE32571, GSE70770, and GSE141551) were combined to validate and predict the value of real hub genes. The results show that the above genes are up-regulated in PCa samples, and higher expression levels show significant association with higher Gleason scores and tumor T stage. Moreover, receiver operating characteristic curve and survival analysis validate the excellent value of hub genes in PCa progression and prognosis. In addition, the protein levels of these four genes also remain higher in tumor tissues when compared with normal tissues. Gene set enrichment analysis and gene set variation analysis for a single gene reveal the close relation with cell proliferation. Meanwhile, 11 small molecular drugs that have the potential to treat PCa were also screened. In conclusion, our research identified four potential prognostic genes and several candidate molecular drugs for treating PCa.

The clinical relevance of gene expression based prognostic signatures in colorectal cancer

Colorectal cancer (CRC) is one of the most prevalent cancers, with more than one million new cases every year. In the last few decades, several advancements in therapeutic and preventative levels have reduced the mortality rate, but new biomarkers are required for improved prognosis. The alterations at the genetic and epigenetic level have been recognized as major players in tumorigenesis. The products of gene expression in the form of mRNA, microRNA, and long-noncoding RNA, have started to emerge as important regulatory molecules, playing an important role in cancer. Gene-expression based prognostic risk scores, which quantify and compare their expression, have emerged as promising biomarkers with enormous clinical value. These composite multi-gene models in which more than one gene is used to predict prognosis have been shown to be significantly effective in identifying patients with multiple clinico-pathological risks like overall mortality, response to chemotherapy, risk of metastasis, etc. The advent of microarray and advanced sequencing technologies have led to the generation of large datasets like TCGA (The Cancer Genome Atlas) and GEO (Gene Expression Omnibus), which have fueled the search for new biomarkers. Continuous evaluation of these candidate biomarkers in clinical settings is promising to improve the management of CRC. These composite gene signatures provide potential in identifying high-risk patients, which might help clinicians to better manage these patients and design appropriate personalized therapeutic interventions. In this review, we emphasize on composite prognostic scores from diverse resources with clinical utility in CRC.

EHD2 Overexpression Suppresses the Proliferation, Migration, and Invasion in Human Colon Cancer

Background: To investigate how EHD2 influences the development of colon cancer.Methods: Immunohistochemistry of 90 colon cancer tissue specimens were determined the expression of EHD2. The lentivirus-EHD2-transfected colon cancer cells were conducted to evaluate the biological behaviors.Results: EHD2 was closely associated with clinic pathological parameters (p < 0.001). EHD2 upregulation was relative with a longer overall survival. The results of the univariate and multivariate analyses indicated that EHD2 could be an independent prognosis marker. EHD2 overexpression suppressed cell invasion and proliferation, but enhanced cell apoptosis and cell cycle arrest.Conclusions: EHD2 might represent a therapeutic target of colon cancer.IMPACT STATEMENTWhat is already known on this subject? Membrane trafficking is crucial for cell proliferation, differentiation and apoptosis, especially tumorigenesis and development. EHD2 proteins play an important role in the regulation of membrane trafficking in endocytosis. EHD2 has been suggested to participate in the occurrence of some malignancies.What are the new findings? EHD2 could be an independent prognosis marker in colon cancer. EHD2 overexpression suppressed cell invasion and proliferation, but enhanced cell apoptosis and cell cycle arrest in vitro. EHD2 overexpression markedly increased the expression of EMT marker E-cadherin in colon cancer.How might it impact on clinical practice in the foreseeable future? EHD2 overexpression may inhibit tumorigenesis in colon cancer through the modulation of E-cadherin, the critical marker of EMT which is closely related to invasion and distant metastasis of tumor cells.

NUSAP1 Promotes Gastric Cancer Tumorigenesis and Progression by Stabilizing the YAP1 Protein

The Yes-associated protein (YAP1) is a main effector of the canonical Hippo pathway, which contributes greatly to tumor initiation, progression, and metastasis in multiple cancers, including gastric cancer (GC). Due to limited knowledge of YAP1 upregulation in cancer, it is a great challenge of therapeutic targets toward the Hippo-YAP1 pathway. Here, we identify nucleolar spindle-associated protein 1 (NUSAP1) as a novel binding partner of YAP1. The upregulation of NUSAP1 is associated with unfavorable clinical outcomes in GC patients, and NUSAP1 depletion impairs its oncogenic properties in vitro and in a xenograft model. Mechanistically, we discovered that NUSAP1 functions as a positive regulator of YAP1 protein stability, thereby inducing the transcription of Hippo pathway downstream target genes, such as CTGF and CYR61. More interestingly, we find that the cancer-promoting effects of NUSAP1 on GC cell growth, migration, and invasion are mainly mediated by YAP1. Furthermore, aberrant expression of NUSAP1 and YAP1 is highly correlated in GC cell lines and tissues. We herein clarify the role of the oncogenic NUSAP1-YAP1 axis in GC tumorigenesis and progression and, therefore, provide novel therapeutic targets for GC treatment.

Nucleolar and spindle-associated protein 1 accelerates cellular proliferation and invasion in nasopharyngeal carcinoma by potentiating Wnt/β-catenin signaling via modulation of GSK-3β

Nucleolar and spindle-associated protein 1 (NUSAP1) is a pivotal tumor-related protein that has been implicated in the progression of broad spectrum of tumors. However, no detailed study of the role of NUSAP1 in nasopharyngeal carcinoma (NPC) has been reported. The aim of this work is to enhance our understanding of NUSAP1 in the progression of NPC. By analyzing data available within the Oncomine database, we found that NUSAP1 expression was elevated in NPC relative to normal tissues. Further, we showed that NUSAP1 expression in clinical specimens of NPC and several NPC cell lines was elevated. Down-regulation of NUSAP1 by gene silencing markedly depleted the capacity of NPC cells to proliferate and invade. Contrastingly, overexpression of NUSAP1 potentiated the proliferative and invasive abilities of NPC cells. Further mechanistic research revealed that NUSAP1 knockdown decreased levels of Wnt/β-catenin signaling in NPC cells via a mechanism associated with downregulation of glycogen synthase kinase-3β (GSK-3β) phosphorylation. However, suppression of GSK-3β markedly abolished the inhibitory effect of NUSAP1 knockdown on Wnt/β-catenin signaling. Further, inhibition of Wnt/β-catenin signaling partially reversed NUSAP1-mediated tumor growth in NPC cells. In addition, NUSAP1 knockdown restrained tumorigenesis of NPC in vivo, and was associated with down-regulation of Wnt/β-catenin signaling. In conclusion, these findings demonstrate that NUSAP1 is capable of accelerating proliferation and invasion in NPC cells by potentiating Wnt/β-catenin signaling. Our study unveils a potential role of NUSAP1 in promoting NPC tumors and suggests that the protein is an attractive antitumor target for NPC treatment.

Identification of Potential Hub Genes and Therapeutic Drugs in Malignant Pleural Mesothelioma by Integrated Bioinformatics Analysis

INTRODUCTION

Malignant pleural mesothelioma (MPM) is closely linked to asbestos exposure and is an extremely aggressive tumor with poor prognosis.

OBJECTIVE

Our study aimed to elucidate hub genes and potential drugs in MPM by integrated bioinformatics analysis.

METHODS

GSE42977 was download from the Gene Expression Omnibus (GEO) database; the differentially expressed genes (DEGs) with adj.p value <0.05 and |logFC| ≥2 were identified. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed by DAVID database. The STRING database was used to construct a protein-protein interaction network, and modules analysis and hub genes acquisition were performed by Cytoscape. The Gene Expression Profiling Interactive Analysis (GEPIA) database was used to assess the impact of hub genes on the prognosis of MPM patients. The Drug-Gene Interaction database (DGIdb) was used to select the related drugs.

RESULTS

A total of 169 upregulated and 70 downregulated DEGs were identified. These DEGs are enriched in the pathway of extracellular matrix-receptor interaction, focal adhesion, PI3K-Akt signaling pathway, and PPAR signaling pathway. Finally, 10 hub genes (CDC20, CDK1, UBE2C, TOP2A, CCNB2, NUSAP1, KIF20A, AURKA, CEP55, and ASPM) were identified, which are considered to be closely related to the poor prognosis of MPM. In addition, 119 related drugs that may have a therapeutic effect on MPM were filtered out.

CONCLUSION

These discovered genes and small-molecule drugs provide some new ideas for further research on MPM.

Identification of Hub Genes in Anaplastic Thyroid Carcinoma: Evidence From Bioinformatics Analysis

Anaplastic thyroid carcinoma (ATC) is a rare type of thyroid cancer that results in fatal clinical outcomes; the pathogenesis of this life-threatening disease has yet to be fully elucidated. This study aims to identify the hub genes of ATC that may play key roles in ATC development and could serve as prognostic biomarkers or therapeutic targets. Two microarray datasets (GSE33630 and GSE53072) were obtained from the Gene Expression Omnibus database; these sets included 16 ATC and 49 normal thyroid samples. Differential expression analyses were performed for each dataset, and 420 genes were screened as common differentially expressed genes using the robust rank aggregation method. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses were conducted to explore the potential bio-functions of these differentially expressed genes (DEGs). The terms and enriched pathways were primarily associated with cell cycle, cell adhesion, and cancer-related signaling pathways. Furthermore, a protein-protein interaction network of DEG expression products was constructed using Cytoscape. Based on the whole network, we identified 7 hub genes that included CDK1, TOP2A, CDC20, KIF11, CCNA2, NUSAP1, and KIF2C. The expression levels of these hub genes were validated using quantitative polymerase chain reaction analyses of clinical specimens. In conclusion, the present study identified several key genes that are involved in ATC development and provides novel insights into the understanding of the molecular mechanisms of ATC development.

Identification of hub genes involved in the occurrence and development of hepatocellular carcinoma via bioinformatics analysis

Hepatocellular carcinoma (HCC) is a heterogeneous malignancy, which is a major cause of cancer morbidity and mortality worldwide. Thus, the aim of the present study was to identify the hub genes and underlying pathways of HCC via bioinformatics analyses. The present study screened three datasets, including GSE112790, GSE84402 and GSE74656 from the Gene Expression Omnibus (GEO) database, and downloaded the RNA-sequencing of HCC from The Cancer Genome Atlas (TCGA) database. The differentially expressed genes (DEGs) in both the GEO and TCGA datasets were filtered, and the screened DEGs were subsequently analyzed for functional enrichment pathways. A protein-protein interaction (PPI) network was constructed, and hub genes were further screened to create the Kaplan-Meier curve using cBioPortal. The expression levels of hub genes were then validated in different datasets using the Oncomine database. In addition, associations between expression and tumor grade, hepatitis virus infection status, satellites and vascular invasion were assessed. A total of 126 DEGs were identified, containing 70 upregulated genes and 56 downregulated genes from the GEO and TCGA databases. By constructing the PPI network, the present study identified hub genes, including cyclin B1 (CCNB1), cell-division cycle protein 20 (CDC20), cyclin-dependent kinase 1, BUB1 mitotic checkpoint serine/threonine kinase β (BUB1B), cyclin A2, nucleolar and spindle associated protein 1, ubiquitin-conjugating enzyme E2 C (UBE2C) and ZW10 interactor. Furthermore, upregulated CCNB1, CDC20, BUB1B and UBE2C expression levels indicated worse disease-free and overall survival. Moreover, a meta-analysis of tumor and healthy tissues in the Oncomine database demonstrated that BUB1B and UBE2C were highly expressed in HCC. The present study also analyzed the data of HCC in TCGA database using univariate and multivariate Cox analyses, and demonstrated that BUB1B and UBE2C may be used as independent prognostic factors. In conclusion, the present study identified several genes and the signaling pathways that were associated with tumorigenesis using bioinformatics analyses, which could be potential targets for the diagnosis and treatment of HCC.

MCM2 and NUSAP1 Are Potential Biomarkers for the Diagnosis and Prognosis of Pancreatic Cancer

Pancreatic cancer (PC) is one of the most malignant tumors. Despite considerable progress in the treatment of PC, the prognosis of patients with PC is poor. The aim of this study was to identify potential biomarkers for the diagnosis and prognosis of PC. First, the original data of three independent mRNA expression datasets were downloaded from the Gene Expression Omnibus and The Cancer Genome Atlas databases and screened for differentially expressed genes (DEGs) using the R software. Subsequently, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses of the DEGs were performed, and a protein-protein interaction (PPI) network was constructed to screen for hub genes. The hub genes were analyzed for genetic variations, as well as for survival, prognostic, and diagnostic value, using the cBioPortal and Gene Expression Profiling Interactive Analysis (GEPIA) databases and the pROC package. After screening for potential biomarkers, the mRNA and protein levels of the biomarkers were verified at the tissue and cellular levels using the Cancer Cell Line Encyclopedia, GEPIA, and the Human Protein Atlas. As a result, a total of 248 DEGs were identified. The GO terms enriched in DEGs were related to the separation of mitotic sister chromatids and the binding of the spindle to the extracellular matrix. The enriched pathways were associated with focal adhesion, ECM-receptor interaction, and phosphatidylinositol 3-kinase (PI3K)/AKT signaling. The top 20 genes were selected from the PPI network as hub genes, and based on the analysis of multiple databases, MCM2 and NUSAP1 were identified as potential biomarkers for the diagnosis and prognosis of PC. In conclusion, our results show that MCM2 and NUSAP1 can be used as potential biomarkers for the diagnosis and prognosis of PC. The study also provides new insights into the underlying molecular mechanisms of PC.

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

OBJECTIVE

We analysed the effect of expression of nucleolar spindle-associated protein 1 (NuSAP1) on the prognosis of breast cancer (BC) and investigated its potential mechanism of tumourigenicity in triple-negative breast cancer (TNBC) cell lines.

MATERIALS AND METHODS

We downloaded the RNA-seq breast cancer (BC) data from The Cancer Genome Atlas (TCGA) and screened for the NuSAP1 gene using R software. The clinical data for patients with BC were screened and analysed using R software. A survival curve was drawn using the Kaplan-Meier Plotter. Cell proliferation and invasion were verified by the Cell Counting Kit-8 and Transwell assays. Expression of NuSAP1, the Wnt/β-catenin pathway, and epithelial-mesenchymal-transition-related proteins in TNBC was detected using real-time quantitative polymerase chain reaction (qRT-PCR) and western blotting (WB).

RESULTS

Expression of NuSAP1 was upregulated in BC. The change in NuSAP1 expression levels was associated with multiple clinicopathological factors, and the higher the expression of NuSAP1 was, the shorter the survival time. In MDA-MB-231 and BT549 cells, knockdown of NuSAP1 expression resulted in a significant decrease in cell proliferation and invasion; a decrease in expression of cyclin D1, vimentin, Slug, Twist, wnt3a, and pβ-catenin; and an increase in expression of e-cadherin. The results of the sh-NuSAP1 + ov-NuSPA1 group were the opposite of the results of the sh-NuSAP1 group.

CONCLUSION

NuSAP1 is a carcinogen that facilitates progression of TNBC through the Wnt/β-catenin and epithelial-mesenchymal transition pathways.

RNA Sequencing and Bioinformatic Analysis on Retinoblastoma Revealing that Cell Cycle Deregulation Is a Key Process in Retinoblastoma Tumorigenesis

PURPOSE

Retinoblastoma (RB) is a primary pediatric ocular malignancy that can be fatal with inadequate treatment. While multimodal treatments are applied for eye salvage, vision loss and metastasis can occur in some patients. The present study aimed to explore key pathways and factors in RB pathogenesis, which could be potential targets for novel RB treatments.

METHODS

RNA sequencing was performed on three RB tissues and referenced with three normal retinas. Differentially expressed genes (DEGs) were identified from sequencing data and further analyzed with clustering analysis, function and pathway enrichment, protein-protein interaction (PPI), and data-mining analysis in order to screen for tumorigenic relevancy.

RESULTS

A total of 331 DEGs were identified by clustering analysis of RB tissues, and the expression patterns were significantly distinguishable from normal retinas. Function and pathway enrichment and PPI analysis together showed that cell cycle was the most prominently upregulated pathway found in RB tissues. Following comprehensive bioinformatic analyses, six key genes relevant to cell cycle regulation were identified, namely BUB1, RRM2, TPX2, UBE2C, NUSAP1, and DTL.

CONCLUSIONS

Cell cycle pathway and six relevant genes may be potential key factors in RB tumorigenesis and laying the foundation for prospective investigation on development of novel targeted therapies.

Decreased Expression of NUSAP1 Predicts Poor Overall Survival in Cervical Cancer

Background: Nucleolar and spindle-associated protein 1 (NUSAP1) was previously reported to be associated with poor prognosis in multiple cancers. In the present study, we comprehensively investigated the clinicopathological features and potential prognostic value of NUSAP1 in cervical squamous cell carcinoma and endocervical adenocarcinoma (CESC).

Methods: The expression profiles of the genes were extracted from Gene Expression Omnibus (GEO), The Cancer Genome Atlas (TCGA), International Cancer Genome Consortium (ICGC), Cancer Cell Line Encyclopedia (CCLE), Gene Expression Profiling Interactive Analysis (GEPIA), and The Human Protein Atlas databases. The association between clinicopathological characteristics and NUSAP1 was analyzed using logistic regression in TCGA patients and receiver operating characteristic (ROC) curve analysis for GSE7803, GSE9750, and GSE63514 datasets. The prognostic value of NUSAP1 in TCGA patients was evaluated using the Kaplan-Meier method and Cox regression. Gene set enrichment analysis (GSEA) was conducted using TCGA dataset.

Results: A total of 68 differentially expressed genes (DEGs) were identified in CESC. ROC analysis of NUSAP1 suggested that the area under the ROC curve was 0.968. Kaplan-Meier survival analysis indicated that CESC with low expression of NUSAP1 has a worse prognosis than CESC with high NUSAP1 expression (P = 0.005). The logistic regression revealed that low NUSAP1 expression in CESC was related to advanced tumor stage in TCGA database. Moreover, Cox regression analysis showed that NUSAP1 expression correlated significantly with prognosis in the case of patients in TCGA database. GSEA demonstrated that CESC patients with high expression of NUSAP1 were enriched in the G2M checkpoint, MYC targets, and breast cancer ZNF217.

Conclusion: The results suggest that identification of DEGs might enhance our understanding of the causes and molecular mechanisms underlying the development of CESC. Moreover, NUSAP1 may play an important role in CESC progression and prognosis and may serve as a valuable indicator of poor survival in CESC.

The role of nucleolar spindle-associated protein 1 in human ovarian cancer

Ovarian cancer (OC) is one of the deadliest malignancies of the female reproductive system. The present study focused on the role of Nucleolar spindle-associated protein 1 (NuSAP1) in OC. Relative expression of NuSAP1 was detected in OC tissues as well as cells. After knocking down NuSAP1 with lentivirus-mediated shRNA and verifying the knockdown efficiency via quantitative real-time polymerase chain reaction and Western blot assays, the cell proliferation, apoptosis, and cell cycle were determined using 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide, colony formation, and flow cytometry, respectively. Transwell assay was conducted to detect the migration and invasion of OC cells. It was showed that NuSAP1 was abundantly expressed in OC tissues and cell lines. After knocking down NuSAP1 in OC cells, in addition to significantly inhibiting proliferation and colony forming ability, it also promotes apoptosis and affects cell cycle distribution. Moreover, cells in the shNuSAP1 group showed significantly suppressed migration and invasion ability compared with that in the shCtrl group. In conclusion, NuSAP1 may act as an oncogenic factor in OC and therefore might serve as an indicator for prognosis and therapeutic target for OC treatment.

Nucleolar and Spindle Associated Protein 1 (NUSAP1) Promotes Bladder Cancer Progression Through the TGF-β Signaling Pathway

Purpose

NUSAP1 has been reported to be involved in the progression of several types of cancer. However, its expression and exact role in bladder cancer (BLCA) remains elusive. The aim of this study was to determine the expression and role of NUSAP1 in BLCA.

Methods

Tissue microarray, real-time PCR, Western blot and immunohistochemistry assays were carried out to determine NUSAP1 expression in BLCA tissues and cells. The biological roles of NUSAP1 were investigated using CCK-8, EdU labeling, flow cytometry, Transwell, and wound healing assays. Additionally, the effect of NUSAP1 on epithelial-mesenchymal transition (EMT) was investigated by Western blotting and real-time PCR.

Results

We found that NUSAP1 was upregulated in BLCA, and its expression was closely related to the poor prognosis of patients. Subsequently, we transfected 5637 and T24 cell lines with NUSAP1 siRNA and an NUSAP1 overexpression plasmid, respectively. NUSAP1 downregulation in 5637 cells inhibited cell proliferation, migration, and invasiveness and enhanced chemosensitivity to gemcitabine, while NUSAP1 overexpression in T24 cells resulted in the inverse effects. Moreover, NUSAP1 regulated EMT via the TGF-β signaling pathway, and when TGF-beta receptor 1 (TGFBR1) was inhibited with the inhibitor SB525334, the invasion and metastasis ability of BLCA cells was significantly suppressed, as well as p-Smad2/3 and vimentin expression.

Conclusion

Our above data demonstrate that NUSAP1 contributes to BLCA progression via the TGF-β signaling pathway.

Downregulation of NUSAP1 suppresses cell proliferation, migration, and invasion via inhibiting mTORC1 signalling pathway in gastric cancer

Gastric cancer (GC) is one of the most common causes of cancer-related death worldwide, and outstanding biomarkers for therapeutic targets or predicting GC survival are still lacking. Increasing evidence indicated that nucleolar and spindle associated protein 1 (NUSAP1) involved in regulating the progression of various cancers; however, its specific role in GC remained unclear. In this study, we found that NUSAP1 was upregulated in the GC tissues and cell lines via analysing data from The Cancer Genome Atlas (TCGA), gene expression omnibus (GEO), qRT-PCR, and western blot assays. Patients with high NUSAP1 expression levels showed shorter free-progression survival (FPS), larger tumour size, and higher lymphatic metastasis rate compared with those with low NUSAP1 expression. Further functional experiments revealed knockdown of NUSAP1 could inhibit the growth, migration, and invasion of GC cells in vitro and vivo. Additionally, silencing NUSAP1 induced G0/G1 phase arrest, apoptosis, and suppressed the epithelial-mesenchymal transition (EMT) process. Finally, we performed gene set enrichment analysis (GSEA) and observed NUSAP1 was positive with mTORC1 signalling pathway, which was verified by the subsequent immunoblotting. In conclusion, our findings suggested that NUSAP1 contributed to GC progression and may act as a potential therapeutic target for GC. SIGNIFICANCE OF THE STUDY: Our results firstly illuminated that NUSAP1 expression was significantly upregulated in GC tissues and predicted poor FPS. Silencing it could attenuate GC progression via inhibiting mTORC1 signalling pathway. Hence, NUSAP1 may act as a promising therapy target for GC.

iTRAQ‑based proteomics analysis of the therapeutic effects of combined anticancer bioactive peptides and oxaliplatin on gastric cancer cells

The combination of chemotherapeutic modalities may be more effective in treating gastric cancer compared with any modality alone. Previous studies have demonstrated that the combination of anticancer bioactive peptides (ACBP) and oxaliplatin (OXA) significantly inhibited the growth of the gastric cancer cell line MKN-45, promoted the apoptosis of MKN-45 cells, and caused an irreversible arrest of the MKN-45 cell cycle in the G2/M phase. In the present study, an isobaric tag for relative and absolute quantitation (iTRAQ)-based quantitative proteomics technique was used to determine the effect of ACBP-OXA treatment on the proteomics profile of MKN-45 cells. Notably, a total of 6,210 proteins were detected. Proteins with a >1.2-fold change in expression (either up- or downregulation) and P<0.05 were considered to be differentially expressed. A total of 256 differentially expressed proteins were identified through alignments with different groups. Compared with the control group, MKN-45 cells treated with ACBP, OXA and ACBP-OXA exhibited 17 (10 up- and 7 downregulated), 111 (27 up- and 84 downregulated) and 128 (53 up- and 75 downregulated) differentially expressed proteins, respectively. Of the 256 differentially expressed proteins, 6 (TPX2, NUSAP1, TOP2A, YAP, MKi-67 and GPC4) were verified by the parallel reaction monitoring method, which revealed that TPX2, NUSAP1, TOP2A, YAP, MKi-67 and GPC4 expression decreased with ACBP-OXA treatment. The cellular localization, functional annotation and biological pathways of differentially expressed proteins were examined by Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analysis. The results indicated that ACBP-OXA may act through the ribosome or the AMP-activated protein kinase (AMPK) signaling pathway, and the AMPK signaling pathway may be an important mediator of the inhibitory effects of ACBP-OXA on MKN-45 gastric cancer cells. In summary, iTRAQ-based proteomics analysis of the effect of ACBP-OXA on MKN-45 cells may guide future therapeutic strategies for gastric cancer. In addition, the present study may help provide new insights into the therapeutic role of combined ACBP and OXA in gastric cancer.

Screening and Identification of Key Biomarkers in Pancreatic Cancer: Evidence from Bioinformatic Analysis

Pancreatic cancer (PC) whose mortality is comparable to morbidity is a highly fatal disease. Early approaches of diagnosis and treatment for PC are quite limited, so it is of great urgency to figure out the exact tumorigenesis and development mechanism of PC. To identify the related molecular markers of pancreatic oncogenesis, we downloaded three microarray datasets (GSE63111, GSE101448, and GSE107610) from Gene Expression Omnibus (GEO) database. The common differentially expressed genes (DEGs) among them were identified, and the corresponding function enrichment analyses were accomplished. The protein-protein interaction network was conducted by Search Tool for the Retrieval of Interacting Genes (STRING), and the corresponding module analysis was accomplished by Cytoscape. There were 55 DEGs found in total. The molecular function and biological processes (BP) of these DEGs mainly include cytokinesis, mitotic nuclear division, cell division, cell proliferation, microtubule-based movement, and mineral absorption. Among the 55 DEGs, 14 hub genes were further confirmed and it was concluded that they mainly function in mitotic cytokinesis, microtubule-based movement, mitotic chromosome condensation, and mitotic spindle assembly from the BP analysis. The survival analysis showed that all the 14 hub genes, especially nucleolar and spindle associated protein 1 and abnormal spindle microtubule assembly, may involve in the tumorigenesis and development of PC. And they might be used as new biomarkers for auxiliary diagnosis and potential targets for immunotherapy of PC.

Identification of modules and functional analysis in CRC subtypes by integrated bioinformatics analysis

Colorectal cancer is one of the top three causes of cancer-related mortality globally, but no predictive molecular biomarkers are currently available for identifying the disease stage of colorectal cancer patients. Common molecular patterns in the disease, beyond superficial manifestations, can be significant in determining treatment choices. In this study, we used microarray data from colorectal cancer and adjacent normal tissue from the GEO database. These data were categorized into four consensus molecular subtypes based on distinct gene expression signatures. Weighted gene-based protein-protein interaction network analysis was performed for each subtype. NUSAP1, CD44, and COL4A1 modules were found to be statistically significant and present among all the subtypes and displayed though similar but not identical functional enrichment results. Reference of the characteristics of the subtypes to functional modules is necessary since the latter can stay resistant to platform changes and technique noise when compared with other analyses. The CMS4-mesenchymal group, which currently has a poor prognosis, was examined in the study. It is composed mainly of genes involved in immune and stromal expression, with modules focused on ECM dysregulation and chemokine biological processes. Hub genes detection and its' mapping into the protein-protein interaction network can be indicative of possible targets against specific modules. This approach identified subtypes using enrichment-oriented analysis in functional modules. Proper annotation of functional analysis of modules from different subtypes of CRC might be directive for finding extra options for treatment targets and guiding clinical routines.

The evolution of spindles and their mechanical implications for cancer metastasis

The mitotic spindle has long been known to play a crucial role in mitosis, orchestrating the segregation of chromosomes into two daughter cells during mitosis with high fidelity. Intracellular forces generated by the mitotic spindle are increasingly well understood, and recent work has revealed that the efficiency and the accuracy of mitosis is ensured by the scaling of mitotic spindle size with cell size. However, the role of the spindle in cancer progression has largely been ignored. Two recent studies point toward the role of mitotic spindle evolution in cancer progression through extracellular force generation. Cancer cells with lengthened spindles exhibit highly increased metastatic potential. Further, interpolar spindle elongation drives protrusive extracellular force generation along the mitotic axis to allow mitotic elongation, a morphological change that is required for cell division. Together, these findings open a new research area studying the role of the mitotic spindle evolution in cancer metastasis.

SNHG6 functions as a competing endogenous RNA to regulate E2F7 expression by sponging miR-26a-5p in lung adenocarcinoma

Increasing evidence has highlighted the pivotal roles of deregulated long non-coding RNAs (lncRNAs) in tumourigenesis. However, the biological functions and mechanisms of lncRNAs in human lung adenocarcinoma (LUAD) remain elusive. Small nucleolar RNA host gene 6 (SNHG6), a novel lncRNA, is aberrantly expressed in various cancers. In this study, SNHG6 was upregulated in LUAD tissues, and its upregulation was positively associated with advanced TNM stage, large tumour size and poor overall survival (OS) in LUAD patients. Gain- and loss-of-function experiments confirmed that SNHG6 promoted cell cycle progression, cell proliferation, migration and invasion, and epithelial-mesenchymal transition (EMT) in vitro. Animal experiments demonstrated that SNHG6 knockdown remarkably inhibited xenograft formation in vivo. Moreover, mechanistic experiments identified that SNHG6 functions as a competing endogenous RNA (ceRNA) through competitively sponging miR-26a-5p to regulate E2F7 expression, cell motility and EMT in LUAD cells. In summary, our findings reveal that SNHG6 may act as an oncogenic lncRNA in LUAD carcinogenesis by regulating the miR-26a-5p/E2F7 axis.