Overexpression of BUB1B contributes to progression of prostate cancer and predicts poor outcome in patients with prostate cancer

Transcriptome Network Analysis Identifies CXCL13-CXCR5 Signaling Modules in the Prostate Tumor Immune Microenvironment

The tumor immune microenvironment (TIME) consists of multiple cell types that contribute to the heterogeneity and complexity of prostate cancer (PCa). In this study, we sought to understand the gene-expression signature of patients with primary prostate tumors by investigating the co-expression profiles of patient samples and their corresponding clinical outcomes, in particular “disease-free months” and “disease reoccurrence”. We tested the hypothesis that the CXCL13-CXCR5 axis is co-expressed with factors supporting TIME and PCa progression. Gene expression counts, with clinical attributes from PCa patients, were acquired from TCGA. Profiles of PCa patients were used to identify key drivers that influence or regulate CXCL13-CXCR5 signaling. Weighted gene co-expression network analysis (WGCNA) was applied to identify co-expression patterns among CXCL13-CXCR5, associated genes, and key genetic drivers within the CXCL13-CXCR5 signaling pathway. The processing of downloaded data files began with quality checks using NOISeq, followed by WGCNA. Our results confirmed the quality of the TCGA transcriptome data, identified 12 co-expression networks, and demonstrated that CXCL13, CXCR5 and associated genes are members of signaling networks (modules) associated with G protein coupled receptor (GPCR) responsiveness, invasion/migration, immune checkpoint, and innate immunity. We also identified top canonical pathways and upstream regulators associated with CXCL13-CXCR5 expression and function.

Evaluation of core serous epithelial ovarian cancer genes as potential prognostic markers and indicators of the underlying molecular mechanisms using an integrated bioinformatics analysis

Ovarian cancer is a major cause of mortality in women. However, the molecular events underlying the pathogenesis of the disease are yet to be fully elucidated. In the present study, an integrated bioinformatics analysis was performed to identify core genes involved in serous epithelial ovarian cancer. A total of three expression datasets were downloaded from the Gene Expression Omnibus database, and included 46 serous epithelial ovarian cancer and 30 ovarian surface epithelium samples. The three datasets were merged, and batch normalization was performed. The normalized merged data were subsequently analyzed for differentially expressed genes (DEGs). In total, 2,212 DEGs were identified, including 1,300 upregulated and 912 downregulated genes. Gene Ontology analysis revealed that these DEGs were primarily involved in ‘regulation of cell cycle’, ‘mitosis’, ‘DNA packaging’ and ‘nucleosome assembly’. The main cellular components included ‘extracellular region part’, ‘chromosome’, ‘extracellular matrix’ and ‘condensed chromosome kinetochore’, whereas the molecular functions included ‘Calcium ion binding’, ‘polysaccharide binding’, ‘enzyme inhibitor activity’, ‘growth factor activity’, ‘cyclin-dependent protein kinase regulator activity’, ‘microtubule motor activity’ and ‘Wnt receptor activity’. Kyoto Encyclopedia of Genes and Genomes pathway analysis revealed that these DEGs were predominantly involved in ‘Wnt signaling pathway’, ‘pathways in cancer’, ‘PI3K-Akt signaling pathway’, ‘cell cycle’, ‘ECM-receptor interaction’, ‘p53 signaling pathway’ and ‘focal adhesion’. The 20 most significant DEGs were identified from the protein-protein interaction network, and Oncomine analysis of these core genes revealed that 13 were upregulated and two were downregulated in serous epithelial ovarian cancer. Survival analysis revealed that cyclin B1, polo like kinase 1, G protein subunit γ transducin 1 and G protein subunit γ 12 are key molecules that may be involved in the prognosis of serous epithelial ovarian cancer. These core genes may provide novel treatment targets, although their roles in the carcinogenesis and prognosis of serous epithelial ovarian cancer require further study.

CDK1, CCNB1, CDC20, BUB1, MAD2L1, MCM3, BUB1B, MCM2, and RFC4 May Be Potential Therapeutic Targets for Hepatocellular Carcinoma Using Integrated Bioinformatic Analysis

Hepatocellular carcinoma (HCC) is a malignant tumor with high mortality. The abnormal expression of genes is significantly related to the occurrence of HCC. The aim of this study was to explore the differentially expressed genes (DEGs) of HCC and to provide bioinformatics basis for the occurrence, prevention and treatment of HCC. The DEGs of HCC and normal tissues in GSE102079, GSE121248, GSE84402 and GSE60502 were obtained using R language. The GO function analysis and KEGG pathway enrichment analysis of DEGs were carried out using the DAVID database. Then, the protein–protein interaction (PPI) network was constructed using the STRING database. Hub genes were screened using Cytoscape software and verified using the GEPIA, UALCAN, and Oncomine database. We used HPA database to exhibit the differences in protein level of hub genes and used LinkedOmics to reveal the relationship between candidate genes and tumor clinical features. Finally, we obtained transcription factor (TF) of hub genes using NetworkAnalyst online tool. A total of 591 overlapping up-regulated genes were identified. These genes were related to cell cycle, DNA replication, pyrimidine metabolism, and p53 signaling pathway. Additionally, the GEPIA database showed that the CDK1, CCNB1, CDC20, BUB1, MAD2L1, MCM3, BUB1B, MCM2, and RFC4 were associated with the poor survival of HCC patients. UALCAN, Oncomine, and HPA databases and qRT-PCR confirmed that these genes were highly expressed in HCC tissues. LinkedOmics database indicated these genes were correlated with overall survival, pathologic stage, pathology T stage, race, and the age of onset. TF analysis showed that MYBL2, KDM5B, MYC, SOX2, and E2F4 were regulators to these nine hub genes. Overexpression of CDK1, CCNB1, CDC20, BUB1, MAD2L1, MCM3, BUB1B, MCM2, and RFC4 in tumor tissues predicted poor survival in HCC. They may be potential therapeutic targets for HCC.

Identification of key genes and microRNAs involved in kidney Wilms tumor by integrated bioinformatics analysis

Wilms tumor (WT) is one of the most common types of pediatric solid tumors; however, its molecular mechanisms remain unclear. The present study aimed to identify key genes and microRNAs (miRNAs), and to predict the underlying molecular mechanisms of WT using integrated bioinformatics analysis. Original gene expression profiles were downloaded from the Gene Expression Omnibus (GEO; accession, GSE66405) and The Cancer Genome Atlas (TCGA) databases. Similarly, miRNA expression patterns were downloaded from GEO (accession, GSE57370) and TCGA. R version 3.5.0 software was used to identify differentially expressed genes (DEGs) and differentially expressed miRNAs (DEMs) using the limma and edgeR packages. Kyoto Encyclopedia of Genes and Genomes pathway and Gene Ontology analyses were performed to examine the biological functions of the DEGs. Additionally, a protein-protein interaction (PPI) network was constructed to screen hub gene modules using Cytoscape software. By predicting target genes of the DEMs and integrating them with DEGs, the present study constructed a miRNA-mRNA regulatory network to predict the possible molecular mechanism of WT. Expression of hub genes was validated using the Oncomine database. A total of 613 genes and 29 miRNAs were identified to be differentially expressed in WT. By constructing a PPI network and screening hub gene modules, 5 upregulated genes, including BUB1 mitotic checkpoint serine/threonine kinase, BUB1B mitotic checkpoint serine/threonine kinase B, cell division cycle protein 45, cyclin B2 and pituitary tumor-transforming 1. These genes were identified to be associated with the cell cycle pathway, which suggested that these genes may serve important roles in WT. In addition, a miRNA-mRNA regulatory network was constructed and comprised 16 DEMs and 19 DEGs. In conclusion, key genes, miRNAs and the mRNA-miRNA regulatory network identified in the present study may improve understanding of the underlying molecular mechanisms in the occurrence and development of WT, and may aid the identification of potential biomarkers and therapeutic targets.

Integrative genomics analysis of hub genes and their relationship with prognosis and signaling pathways in esophageal squamous cell carcinoma

The main purpose of the present study was to recognize the integrative genomics analysis of hub genes and their relationship with prognosis and signaling pathways in esophageal squamous cell carcinoma (ESCC). The mRNA gene expression profile data of GSE38129 were downloaded from the Gene Expression Omnibus database, which included 30 ESCC and 30 normal tissue samples. The differentially expressed genes (DEGs) between ESCC and normal samples were identified using the GEO2R tool. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed to identify the functions and related pathways of the genes. The protein‑protein interaction (PPI) network of these DEGs was constructed with the Search Tool for the Retrieval of Interacting Genes and visualized with a molecular complex detection plug‑in via Cytoscape. The top five important modules were selected from the PPI network. A total of 928 DEGs, including ephrin‑A1 (EFNA1), collagen type IV α1 (COL4A1),  C‑X‑C chemokine receptor 2 (CXCR2), adrenoreceptor β2 (ADRB2), P2RY14, BUB1B, cyclin A2 (CCNA2), checkpoint kinase 1 (CHEK1), TTK, pituitary tumor transforming gene 1 (PTTG1) and COL5A1, including 498 upregulated genes, were mainly enriched in the 'cell cycle', 'DNA replication' and 'mitotic nuclear division', whereas 430 downregulated genes were enriched in 'oxidation‑reduction process', 'xenobiotic metabolic process' and 'cell‑cell adhesion'. The KEGG analysis revealed that 'ECM‑receptor interaction', 'cell cycle' and 'p53 signaling pathway' were the most relevant pathways. According to the degree of connectivity and adjusted P‑value, eight core genes were selected, among which those with the highest correlation were CHEK1, BUB1B, PTTG1, COL4A1 and CXCR2. Gene Expression Profiling Interactive Analysis in The Cancer Genome Atlas database for overall survival (OS) was applied among these genes and revealed that EFNA1 and COL4A1 were significantly associated with a short OS in 182 patients. Immunohistochemical results revealed that the expression of PTTG1 in esophageal carcinoma tissues was higher than that in normal tissues. Therefore, these genes may serve as crucial predictors for the prognosis of ESCC.

Insights into a Crucial Role of TRIP13 in Human Cancer

Thyroid Hormone Receptor Interacting Protein 13 (TRIP13) plays a key role in regulating mitotic processes, including spindle assembly checkpoint and DNA repair pathways, which may account for Chromosome instability (CIN). As CIN is a predominant hallmark of cancer, TRIP13 may act as a tumor susceptibility locus. Amplification of TRIP13 has been observed in various human cancers and implicated in several aspects of malignant transformation, including cancer cell proliferation, drug resistance and tumor progression. Here, we discussed the functional significance of TRIP13 in cell progression, highlighted the recent findings on the aberrant expression in human cancers and emphasized its significance for the therapeutic potential.

Identification of five hub genes as monitoring biomarkers for breast cancer metastasis in silico

BACKGROUND

Breast cancer is one of the most common endocrine cancers among females worldwide. Distant metastasis of breast cancer is causing an increasing number of breast cancer-related deaths. However, the potential mechanisms of metastasis and candidate biomarkers remain to be further explored.

RESULTS

The gene expression profiles of GSE102484 were downloaded from the Gene Expression Omnibus (GEO) database. Weighted gene co-expression network analysis (WGCNA) was used to screen for the most potent gene modules associated with the metastatic risk of breast cancer, and a total of 12 modules were identified based on the analysis. In the most significant module (R2 = 0.68), 21 network hub genes (MM > 0.90) were retained for further analyses. Next, protein-protein interaction (PPI) networks were used to further explore the biomarkers with the most interactions in gene modules. According to the PPI networks, five hub genes (TPX2, KIF2C, CDCA8, BUB1B, and CCNA2) were identified as key genes associated with breast cancer progression. Furthermore, the prognostic value and differential expression of these genes were validated based on data from The Cancer Genome Atlas (TCGA) and Kaplan-Meier (KM) Plotter. Receiver operating characteristic (ROC) curve analysis revealed that the mRNA expression levels of these five hub genes showed excellent diagnostic value for breast cancer and adjacent tissues. Moreover, these five hub genes were significantly associated with worse distant metastasis-free survival (DMFS) in the patient cohort based on KM Plotter.

CONCLUSION

Five hub genes (TPX2, KIF2C, CDCA8, BUB1B, and CCNA2) associated with the risk of distant metastasis were extracted for further research, which might be used as biomarkers to predict distant metastasis of breast cancer.

Integrated bioinformatics analysis to identify 15 hub genes in breast cancer

The aim of the present study was to identify the hub genes and provide insight into the tumorigenesis and development of breast cancer. To examine the hub genes in breast cancer, integrated bioinformatics analysis was performed. Gene expression profiles were obtained from the Gene Expression Omnibus (GEO) database and the differentially expressed genes (DEGs) were identified using the ‘limma’ package in R. Gene Ontology enrichment analysis and Kyoto Encyclopedia of Genes and Genomes pathway analysis was used to determine the functional annotations and potential pathways of the DEGs. Subsequently, a protein-protein interaction network analysis and weighted correlation network analysis (WGCNA) were conducted to identify hub genes. To confirm the reliability of the identified hub genes, RNA gene expression profiles were obtained from The Cancer Genome Atlas (TCGA)-breast cancer database, and WGCNA was used to screen for genes that were markedly correlated with breast cancer. By combining the results from the GEO and TCGA datasets, 15 hub genes were identified to be associated with breast cancer pathophysiology. Overall survival analysis was performed to examine the association between the expression of hub genes and the overall survival time of patients with breast cancer. Higher expression of all hub genes was associated with significantly shorter overall survival in patients with breast cancer compared with patients with lower levels of expression of the respective gene.

Bioinformatics analysis of a long non‑coding RNA and mRNA regulation network in rats with middle cerebral artery occlusion based on RNA sequencing

Long non‑coding RNAs (lncRNAs) have been proven to be critical gene regulators of development and disease. The main aim of the present study was to elucidate the lncRNA‑mRNA regulation network in ischemic stroke induced by middle cerebral artery occlusion (MCAO) using RNA sequencing (RNA‑seq) in rats. lncRNA expression profiles were screened in brain tissues to identify a number of differentially expressed lncRNAs (DELs) and genes (DEGs) by RNA‑seq. Reverse transcription‑quantitative polymerase chain reaction was performed to further confirm the lncRNA expression data. Furthermore, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were used to mine mRNA functions, and a lncRNA‑mRNA network was constructed. Additionally, cis‑ and trans‑regulatory gene analyses of DELs were predicted. A total of 134 DELs (fold change >2, false discovery rate <0.05) and 1,006 DEGs (fold change >2 and P<0.05) were identified. Eighteen lncRNAs were predicted to regulate heme oxygenase 1, mitotic checkpoint serine/threonine kinase B, chemokine ligand 2 and DNA Topoisomerase IIα, amongst other genes. These genes are all associated with a cellular response to inorganic substances, alkaloids, estradiol, reactive oxygen species, metal ions, oxidative stress, and are associated with metabolic pathways, chemokine signaling pathways, malaria, Parkinson's disease, the cell cycle and other GO and KEGG pathway enrichments. The present study identifies novel DELs and an lncRNA‑mRNA regulatory network that may allow for an improved understanding of the molecular mechanism of ischemic stroke induced by MCAO.

Identification of significant genes with poor prognosis in ovarian cancer via bioinformatical analysis

Ovarian cancer (OC) is the highest frequent malignant gynecologic tumor with very complicated pathogenesis. The purpose of the present academic work was to identify significant genes with poor outcome and their underlying mechanisms. Gene expression profiles of GSE36668, GSE14407 and GSE18520 were available from GEO database. There are 69 OC tissues and 26 normal tissues in the three profile datasets. Differentially expressed genes (DEGs) between OC tissues and normal ovarian (OV) tissues were picked out by GEO2R tool and Venn diagram software. Next, we made use of the Database for Annotation, Visualization and Integrated Discovery (DAVID) to analyze Kyoto Encyclopedia of Gene and Genome (KEGG) pathway and gene ontology (GO). Then protein-protein interaction (PPI) of these DEGs was visualized by Cytoscape with Search Tool for the Retrieval of Interacting Genes (STRING). There were total of 216 consistently expressed genes in the three datasets, including 110 up-regulated genes enriched in cell division, sister chromatid cohesion, mitotic nuclear division, regulation of cell cycle, protein localization to kinetochore, cell proliferation and Cell cycle, progesterone-mediated oocyte maturation and p53 signaling pathway, while 106 down-regulated genes enriched in palate development, blood coagulation, positive regulation of transcription from RNA polymerase II promoter, axonogenesis, receptor internalization, negative regulation of transcription from RNA polymerase II promoter and no significant signaling pathways. Of PPI network analyzed by Molecular Complex Detection (MCODE) plug-in, all 33 up-regulated genes were selected. Furthermore, for the analysis of overall survival among those genes, Kaplan–Meier analysis was implemented and 20 of 33 genes had a significantly worse prognosis. For validation in Gene Expression Profiling Interactive Analysis (GEPIA), 15 of 20 genes were discovered highly expressed in OC tissues compared to normal OV tissues. Furthermore, four genes (BUB1B, BUB1, TTK and CCNB1) were found to significantly enrich in the cell cycle pathway via re-analysis of DAVID. In conclusion, we have identified four significant up-regulated DEGs with poor prognosis in OC on the basis of integrated bioinformatical methods, which could be potential therapeutic targets for OC patients.

Overexpression of BUB1B, CCNA2, CDC20, and CDK1 in tumor tissues predicts poor survival in pancreatic ductal adenocarcinoma

Overexpressed genes in tumors usually contributed to aggressiveness in pancreatic ductal adenocarcinoma (PDAC). Using Gene Expression Omnibus (GEO) profiles including GSE46234, GSE71989, and GSE107610, we detected overexpressed genes in tumors with R program, which were enriched by Kyoto Encyclopedia of Genes and Genomes (KEGG), Gene ontology (GO), and Reactome pathway databases. Then, we performed a survival analysis of enriched genes based on TCGA profile. Our results revealed that high BUB1B, CCNA2, CDC20, and CDK1 expression in tumors was significantly associated with worse overall survival (OS) (Log rank P=0.00338, P=0.0447, P=0.00965, and P=0.00479, respectively), which was validated using a Kaplan–Meier plotter with a median cutoff (Log rank P=0.028, P=0.0035, P=0.039, and P=0.0033, respectively). Moreover, overexpression of BUB1B, CCNA2, CDC20, and CDK1 in tumor tissues was significantly associated with disease-free survival (DFS) in PDAC patients (Log rank P=0.00565, P=0.0357, P=0.00104, and P=0.00121, respectively). BUB1B, CCNA2, CDC20, and CDK1 were significantly overexpressed in deceased PDAC patients (all P<0.01) and in patients with recurrence/disease progression (all P<0.05). In addition, PDAC patients with neoplasms of histologic grade G3-4 had significantly higher BUB1B, CCNA2 and CDC20 levels (all P<0.05). In conclusion, the up-regulation of BUB1B, CCNA2, CDC20, CDK1, and WEE1 in tumor tissues are associated with worse OS and DFS in PDAC and is correlated with advanced tumor stage and tumor development.

Co-expression Network Analysis Identifies Four Hub Genes Associated With Prognosis in Soft Tissue Sarcoma

Background: Soft tissue sarcomas (STS) are heterogeneous tumors derived from mesenchymal cells that differentiate into soft tissues. The prognosis of patients who present with an STS is influenced by the regulation of a complex gene network. Methods: Weighted gene co-expression network analysis (WGCNA) was performed to identify gene modules associated with STS (Samples = 156). Results: Among the 11 modules identified, the black and blue modules were highly correlated with STS. However, using preservation analysis, the black module demonstrated low preservation, therefore the blue module was chosen as the module of interest. Furthermore, a total of 20 network hub genes were identified in the blue module, 12 of which were also hub nodes in the protein-protein interaction network of the module genes. Following additional verification, 4 of 12 genes (RRM2, BUB1B, CENPF, and KIF20A) demonstrated poorer overall survival and disease-free survival rate in the test datasets. In addition, gene set enrichment analysis (GSEA) demonstrated that samples with a high level of blue module eigengene (ME) were enriched in cell cycle and metabolism associated signaling pathways. Conclusion: In summary, co-expression network analysis identified four hub genes associated with prognosis for STS, which may diminish the prognosis by influencing cell cycle and metabolism associated signaling pathways.

Integrated analysis reveals key genes with prognostic value in lung adenocarcinoma

Background

Lung cancer is one of the most common malignant tumors. Despite advances in lung cancer therapies, prognosis of non-small-cell lung cancer is still unfavorable. The aim of this study was to identify the prognostic value of key genes in lung tumorigenesis.

Methods

Differentially expressed genes (DEGs) were screened out by GEO2R from three Gene Expression Omnibus cohorts. Common DEGs were selected for Kyoto Encyclopedia of Genes and Genomes pathway analysis and Gene Ontology enrichment analysis. Protein– protein interaction networks were constructed by the STRING database and visualized by Cytoscape software. Hub genes, filtered from the CytoHubba, were validated using the Gene Expression Profiling Interactive Analysis database, and their genomic alterations were identified by performing the cBioportal. Finally, overall survival analysis of hub genes was performed using Kaplan–Meier Plotter.

Results

From three datasets, 169 DEGs (70 upregulated and 99 downregulated) were identified. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses showed that upregulated DEGs were significantly enriched in cell cycle, p53 pathway, and extracellular matrix–receptor interactions; the downregulated DEGs were significantly enriched in PPAR pathway and tyrosine metabolism. The protein–protein interaction network consisted of 71 nodes and 305 edges, including 49 upregulated and 22 downregulated genes. The hub genes, including AURKB, BUB1B, KIF2C, HMMR, CENPF, and CENPU, were overexpressed compared with the normal group by Gene Expression Profiling Interactive Analysis analysis, and associated with reduced overall survival in lung cancer patients. In the genomic alterations analysis, two hotspot mutations (S2021C/F and E314K/V) were identified in Pfam protein domains.

Conclusion

DEGs, including AURKB, BUB1B, KIF2C, HMMR, CENPF, and CENPU, might be potential biomarkers for the prognosis and treatment of lung adenocarcinoma.

Aneuploidy: Cancer strength or vulnerability?

Aneuploidy is a very rare and tissue‐specific event in normal conditions, occurring in a low number of brain and liver cells. Its frequency increases in age‐related disorders and is one of the hallmarks of cancer. Aneuploidy has been associated with defects in the spindle assembly checkpoint (SAC). However, the relationship between chromosome number alterations, SAC genes and tumor susceptibility remains unclear. Here, we provide a comprehensive review of SAC gene alterations at genomic and transcriptional level across human cancers and discuss the oncogenic and tumor suppressor functions of aneuploidy. SAC genes are rarely mutated but frequently overexpressed, with a negative prognostic impact on different tumor types. Both increased and decreased SAC gene expression show oncogenic potential in mice. SAC gene upregulation may drive aneuploidization and tumorigenesis through mitotic delay, coupled with additional oncogenic functions outside mitosis. The genomic background and environmental conditions influence the fate of aneuploid cells. Aneuploidy reduces cellular fitness. It induces growth and contact inhibition, mitotic and proteotoxic stress, cell senescence and production of reactive oxygen species. However, aneuploidy confers an evolutionary flexibility by favoring genome and chromosome instability (CIN), cellular adaptation, stem cell‐like properties and immune escape. These properties represent the driving force of aneuploid cancers, especially under conditions of stress and pharmacological pressure, and are currently under investigation as potential therapeutic targets. Indeed, promising results have been obtained from synthetic lethal combinations exploiting CIN, mitotic defects, and aneuploidy‐tolerating mechanisms as cancer vulnerability.

Comprehensive analysis of multi Ewing sarcoma microarray datasets identifies several prognosis biomarkers

Ewing sarcoma (ES) is a common primary malignancy in children and adolescents. Progression of treatment methods hasn't contributed a lot to the imrovement of prognosis. To identify potential prognostic biomarkers, a meta-analysis pipeline of multi-gene expression datasets for ES from the Gene Expression Omnibus (GEO) was performed. Three datasets were screened and differential expression genes (DEGs) in ES samples compared with normal tissues were identified through limma package and subjected to network analysis. As a result, 1,470 DEGs were obtained which were mainly involved in biological processes associated with immune response and transcription regulation. Network analysis obtained 22 core genes with high network degree and fold change. Kaplan-Meier analysis based on ES datasets from The Cancer Genome Atlas identified five genes, including glycogen phosphorylase, muscle-associated, myocyte-specific enhancer factor 2C, tripartite motif containing 63, budding uninhibited by benzimidazoses1 and Ras GTPase-activating protein 1, whose altered expression profiles are significantly associated with survival. Changes of their expression values were further confirmed through RT-qPCR in ES cell and normal cell lines. Those genes may be considered as potential prognostic biomarkers of ES and should be helpful for its early diagnosis and treatment.

Dysfunction of Sister Chromatids Separation Promotes Progression of Hepatocellular Carcinoma According to Analysis of Gene Expression Profiling

Despite studying the various molecular mechanisms of hepatocellular carcinoma (HCC), effective drugs and biomarkers in HCC therapy are still scarce. The present study was designed to investigate dysregulated pathways, novel biomarkers and therapeutic targets for HCC. The gene expression dataset of GSE14520, which included 362 tumor and their paired non-tumor tissues of HCC, was extracted for processing by the Robust multi-array average (RMA) algorithm in the R environment. SAM methods were leveraged to identify differentially expressed genes (DEGs). Functional analysis of DEGs was performed using DAVID. The GeneMania and Cytohubba were used to construct the PPI network. To avoid individual bias, GSEA and survival analysis were employed to verify the results. The results of these analyses indicated that separation of sister chromatids was the most aberrant phase in the progression of HCC, and the most frequently involved genes, EZH2, GINS1, TPX2, CENPF, and BUB1B, require further study to be used as drug targets or biomarkers in diagnosis and treatment of HCC.

Overexpression of TPX2 is associated with progression and prognosis of prostate cancer

Targeting protein for Xenopus kinesin-like protein 2 (TPX2) activates Aurora kinase A during mitosis and targets its activity to the mitotic spindle, serving an important role in mitosis. It has been associated with different types of cancer and is considered to promote tumor growth. The aim of the present study was to explore the role of TPX2 in diagnosing prostate cancer (PCa). It was identified that TPX2 expression in PCa tissues was increased compared with benign prostate tissues. Microarray analysis demonstrated that TPX2 was positively associated with the Gleason score, tumor-node-metastasis (TNM) stage, clinicopathological stage, metastasis, overall survival and biochemical relapse-free survival. In vitro studies revealed that the high expression of TPX2 in PCa cells improved proliferative, invasive and migratory abilities, and repressed apoptosis of the PCa cells, without affecting tolerance to docetaxel. The results suggested that TPX2 serves as a tumorigenesis-promoting gene in PCa, and a potential therapeutic target for patients with PCa.

Upregulation of BUB1B, CCNB1, CDC7, CDC20, and MCM3 in Tumor Tissues Predicted Worse Overall Survival and Disease-Free Survival in Hepatocellular Carcinoma Patients

OBJECTIVE

To evaluate the association between upregulated differentially expressed genes (DEGs) and the outcomes of patients with hepatocellular carcinoma (HCC).

METHODS

Using Gene Expression Omnibus (GEO) datasets including GSE45436, GSE55092, GSE60502, GSE84402, and GSE17548, we detected upregulated DEGs in tumors. KEGG, GO, and Reactome enrichment analysis of the DEGs was conducted to clarify their function. The impact of the upregulated DEGs on patients' survival was analyzed based on TCGA profile.

RESULTS

161 shared upregulated DEGs were identified among GSE45436, GSE55092, GSE60502, and GSE84402 profiles. Cell cycle was the shared pathway/biological process in the gene sets investigation among databases of KEGG, GO, and Reactome. After being validated in GSE17548, 13 genes including BUB1B, CCNA2, CCNB1, CCNE2, CDC20, CDC6, CDC7, CDK1, CDK4, CDKN2A, CHEK1, MAD2L1, and MCM3 in cell cycle pathway were shared in the three databases for enrichment. The expression of BUB1B, CCNB1, CDC7, CDC20, and MCM3 was upregulated in HCC tissues when compared with adjacent normal tissues in 6.67%, 7.5%, 8.06%, 5.56%, and 9.72% of HCC patients, respectively. Overexpression of BUB1B, CCNB1, CDC7, CDC20, and MCM3 in HCC tissues accounted for poorer overall survival (OS) and disease-free survival (DFS) in HCC patients (all log rank P < 0.05). BUB1B, CCNB1, CDC7, CDC20, and MCM3 were all overexpressed in HCC patients with neoplasm histologic grade G3-4 compared to those with G1-2 (all P < 0.05). BUB1B, CCNB1, and CDC20 were significantly upregulated in HCC patients with vascular invasion (all P < 0.05). Additionally, levels of BUB1B, CCNB1, CDC7, and CDC20 were significantly higher in HCC patients deceased, recurred, or progressed (all P < 0.05).

CONCLUSION

Correlated with advanced histologic grade and/or vascular invasion, upregulation of BUB1B, CCNB1, CDC7, CDC20, and MCM3 in HCC tissues predicted worse OS and DFS in HCC patients. These genes could be novel therapeutic targets for HCC treatment.

Comprehensive analysis of differentially expressed genes associated with PLK1 in bladder cancer

BACKGROUND

The significance of PLK1 (polo-like kinase 1) has become increasingly essential as both a biomarker and a target for cancer treatment. Here, we aimed to determine the downstream genes of PLK1 and their effects on the carcinogenesis and progression of bladder cancer.

METHODS

Specific siRNA was utilized to silence the target gene expression. The cell proliferation, invasion and migration of bladder cancer cells by MTT assay, BrdU assay and transwell assay. The differential expression genes were identified using Affymetrix HTA2.0 Array. The KEGG, GO and STRING analysis were used to analyze the signaling pathway and protein-protein interaction. Spearman analysis was used to analyze the correlation between protein and protein, between protein and clincopathologic characteristics.

RESULTS

PLK1 siRNA hindered the proliferation, invasion and migration of bladder cancer cells, as determined by the MTT, BrdU and transwell assays. A total of 561 differentially expressed genes were identified using an Affymetrix HTA2.0 Array in PLK1 knockdown T24 cells. According to KEGG, GO and STRING analysis, five key genes (BUB1B, CCNB1, CDC25A, FBXO5, NDC80) were determined to be involved in cell proliferation, invasion and migration. PLK1 knockdown decreased BUB1B, CCNB1, CDC25A and NDC80 expressions but increased FBXO5 expression. BUB1B, CCNB1, CDC25A and NDC80 were positively correlated with cell proliferation, invasion, migration and PLK1 expression in tissues, but FBXO5 was negatively correlated with each of those factors. The results showed that the five genes expressions were significantly correlation with the PLK1 expression in normal bladder tissues and bladder cancer tissues. Four of them (BUB1B, CCNB1, CDC25A, NDC80) were obviously positive correlations with pT stage and metastasis. But FBXO5 was negative correlated with pT stage and metastasis. Furthermore, significant correlations were found between CCNB1 or CDC25A or NDC80 and histological grade; between BUB1B or NDC80 and recurrence.

CONCLUSION

Five downstream genes of PLK1 were associated with the regulation of cell proliferation, invasion and migration in bladder cancer. Furthermore, these genes may play important roles in bladder cancer and become important biomarkers and targets for cancer treatment.

The FOXM1/BUB1B signaling pathway is essential for the tumorigenicity and radioresistance of glioblastoma

Accumulating evidence indicates that mitotic checkpoint serine/threonine kinase B (BUB1B) plays a critical role in multiple types of cancer. However, the biological function and molecular regulatory mechanism of BUB1B in glioblastoma (GBM) remain unclear. In the present study, we identified that BUB1B expression was enriched in GBM tumors and was functionally required for tumor proliferation both in vitro and in vivo. Clinically, BUB1B expression was associated with poor prognosis in GBM patients and BUB1B-dependent radioresistance in GBM was decreased by targeting BUB1B via shRNAs. Mechanistically, forkhead box protein M1 (FOXM1) transcriptionally regulated BUB1B expression by binding to and then activating the BUB1B promoter. Therapeutically, we found that FOXM1 inhibitor attenuated tumorigenesis and radioresistance of GBM both in vitro and in vivo. Altogether, BUB1B promotes tumor proliferation and induces radioresistance in GBM, indicating that BUB1B could be a potential therapeutic target for GBM.