BUB1 promotes proliferation of liver cancer cells by activating SMAD2 phosphorylation

BUB1/KIF14 complex promotes anaplastic thyroid carcinoma progression by inducing chromosome instability

Chromosome instability (CIN) is a common contributor driving the formation and progression of anaplastic thyroid cancer (ATC), but its mechanism remains unclear. The BUB1 mitotic checkpoint serine/threonine kinase (BUB1) is responsible for the alignment of mitotic chromosomes, which has not been thoroughly studied in ATC. Our research demonstrated that BUB1 was remarkably upregulated and closely related to worse progression-free survival. Knockdown of BUB1 attenuated cell viability, invasion, migration and induced cell cycle arrests, whereas overexpression of BUB1 promoted the cell cycle progression of papillary thyroid cancer cells. BUB1 knockdown remarkably repressed tumour growth and tumour formation of nude mice with ATC xenografts and suppressed tumour metastasis in a zebrafish xenograft model. Inhibition of BUB1 by its inhibitor BAY-1816032 also exhibited considerable anti-tumour activity. Further studies showed that enforced expression of BUB1 evoked CIN in ATC cells. BUB1 induced CIN through phosphorylation of KIF14 at serine1292 (Ser1292 ). Overexpression of the KIF14ΔSer1292 mutant was unable to facilitate the aggressiveness of ATC cells when compared with that of the wild type. Collectively, these findings demonstrate that the BUB1/KIF14 complex drives the aggressiveness of ATC by inducing CIN.

Expression of the checkpoint kinase BUB1 is a predictor of response to cancer therapies

The identification of clinically-relevant biomarkers is of upmost importance for the management of cancer, from diagnosis to treatment choices. We performed a pan-cancer analysis of the mitotic checkpoint budding uninhibited by benzimidazole 1 gene BUB1, in the attempt to ascertain its diagnostic and prognostic values, specifically in the context of drug response. BUB1 was found to be overexpressed in the majority of cancers, and particularly elevated in clinically aggressive molecular subtypes. Its expression was correlated with clinico-phenotypic features, notably tumour staging, size, invasion, hypoxia, and stemness. In terms of prognostic value, the expression of BUB1 bore differential clinical outcomes depending on the treatment administered in TCGA cancer cohorts, suggesting sensitivity or resistance, depending on the expression levels. We also integrated in vitro drug sensitivity data from public projects based on correlation between drug efficacy and BUB1 expression to produce a list of candidate compounds with differential responses according to BUB1 levels. Gene Ontology enrichment analyses revealed that BUB1 overexpression in cancer is associated with biological processes related to mitosis and chromosome segregation machinery, reflecting the mechanisms of action of drugs with a differential effect based on BUB1 expression.

Anti-tumor target screening of sea cucumber saponin Frondoside A: a bioinformatics and molecular docking analysis

Cancer remains the leading cause of death worldwide. In spite of significant advances in targeted and immunotherapeutic approaches, clinical outcomes for cancer remain poor. The aim of the present study was to investigate the potential mechanisms and therapeutic targets of Frondoside A for the treatment of liver, pancreatic, and bladder cancers. The data presented in our study demonstrated that Frondoside A reduced the viability and migration of HepG2, Panc02, and UM-UC-3 cancer cell in vitro. Moreover, we utilized the GEO database to screen and identify for differentially expressed genes (DEGs) in liver, pancreatic, and bladder cancers, which resulted in the identification of 714, 357, and 101 DEGs, respectively. Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway annotation were performed using the Metascape database for DEGs that were significantly associated with cancer development. The protein-protein interaction (PPI) networks of the identified DEGs in liver, pancreatic, and bladder cancers were analyzed using Cytoscape 3.9.0 software, and subsequently identified potential key genes that were associated with these networks. Subsequently, their prognostic values were assessed by gene expression level analysis and Kaplan-Meier survival analysis (GEPIA). Furthermore, we utilized TIMER 2.0 to investigate the correlation between the expression of the identified key gene and cancer immune infiltration. Finally, molecular docking simulations were performed to assess the affinity of Frondoside A and key genes. Our results showed a significant correlation between these DEGs and cancer progression. Combined, these analyses revealed that Frondoside A involves in the regulation of multiple pathways, such as drug metabolism, cell cycle in liver cancer by inhibiting the expression of CDK1, TOP2A, CDC20, and KIF20A, and regulates protein digestion and absorption, receptor interaction in pancreatic cancer by down-regulation of ASPM, TOP2A, DLGAP5, TPX2, KIF23, MELK, LAMA3, and ANLN. While in bladder cancer, Frondoside A regulates muscle contraction, complement and coagulation cascade by increase FLNC expression. In conclusion, the present study offers valuable insights into the molecular mechanism underlying the anticancer effects of Frondoside A, and suggests that Frondoside A can be used as a functional food supplement or further developed as a natural anti-cancer drug.

Inhibition of the serine/threonine kinase BUB1 reverses taxane resistance in prostate cancer

Men with incurable castration resistant prostate cancer (CRPC) are typically treated with taxanes; however, drug resistance rapidly develops. We previously identified a clinically relevant seven gene network in aggressive CRPC, which includes the spindle assembly checkpoint (SAC) kinase BUB1. Since SAC is deregulated in taxane resistant PC, we evaluated BUB1 and found that it was over-expressed in advanced PC patient datasets and taxane resistant PC cells. Treatment with a specific BUB1 kinase inhibitor re-sensitized resistant CRPC cells, including cells expressing constitutively active androgen receptor (AR) variants, to clinically used taxanes. Consistent with a role of AR variants in taxane resistance, ectopically expressed AR-V7 increased BUB1 levels and reduced sensitivity to taxanes. This work shows that disruption of BUB1 kinase activity reverted resistance to taxanes, which is essential to advancing BUB1 as a potential therapeutic target for intractable chemotherapy resistant CRPC including AR variant driven CRPC, which lacks durable treatment options.

Budding uninhibited by benzimidazoles 1 promotes cell proliferation, invasion, and epithelial-mesenchymal transition via the Wnt/β-catenin signaling in glioblastoma

The pathogenesis and progression of GBM (glioblastoma), as one of the most frequently occurring malignancies of the central nervous system, are regulated by several genes. BUB1 (budding uninhibited by benzimidazoles 1) is a mitotic checkpoint that plays an important role in chromosome segregation as well as in various tumors. However, its role in glioma is unknown. The current study discovered prominently elevated BUB1 in glioma and a significant relationship between BUB1 expression, a high World Health Organization grade, and a poor prognosis in glioma patients. Moreover, BUB1 triggered EMT (epithelial-mesenchymal transition) apart from promoting glioma cell proliferation, migration, and infiltration. Besides, BUB1 promoted EMT by activating the Wnt/β-catenin axis. As implied by our study, BUB1 probably has the potential as a target for GBM management.

Immunological role and prognostic value of the SKA family in pan-cancer analysis

Background

The spindle and kinetochore associated (SKA) complex, which plays important roles in proper chromosome segregation during mitosis by maintaining the stabilization of kinetochore-spindle microtubule attachment during mitosis, has recently been reported to exert regulatory effects on the initiation and progression of various human cancer types. Nevertheless, the prognostic significance and immune infiltration of the SKA family across cancers have not been well elucidated.

Methods

Using data from three large public datasets, including The Cancer Genome Atlas, Genotype-Tissue Expression, and Gene Expression Omnibus databases, a novel scoring system (termed the SKA score) was developed to quantify the SKA family level across cancers. We then evaluated the prognostic impact of the SKA score on survival and assessed the effect of the SKA score on immunotherapy at the pan-cancer level using multiomics bioinformatic analyses. The correlation of the SKA score and the tumor microenvironment (TME) was also explored in depth. Potential small molecular compounds and chemotherapeutic agents were assessed by CTRP and GDSC analyses. Immunohistochemistry was performed to verify the expression of the SKA family genes.

Results

Our results demonstrated a close correlation between the SKA score and tumor development and prognosis in multiple cancers. The SKA score was positively related to cell cycle pathways and DNA replication across cancers, such as E2F targets, the G2M checkpoint, MYC targets V1/V2, mitotic spindles and DNA repair. Additionally, the SKA score was negatively related to the infiltration of various immune cells with antitumor effects in the TME. In addition, the potential value of the SKA score was identified to predict immunotherapy response for melanoma and bladder cancer. We also demonstrated a correlation between SKA1/2/3 and the response to drug treatment across cancers and the promising potential of the SKA complex and its genes as therapeutic targets in cancer. Immunohistochemistry demonstrated that the expression differences of SKA1/2/3 were significant between the breast cancer group and the paracancerous group.

Conclusion

The SKA score plays a critical role in 33 cancer types and is highly related to tumor prognosis. Patients with elevated SKA scores have a clear immunosuppressive TME. The SKA score may serve as a predictor for patients receiving anti-PD-1/L1 therapy.

Effects of Sodium Chromate Exposure on Gene Expression Profiles of Primary Rat Hepatocytes (In Vitro)

Chromium exposure has adverse impacts on human health and the environment, whereas chromate-induced hepatotoxicity's detailed mechanism is still unclear. Therefore, the purpose of the current study was to reveal the crucial signaling pathways and genes linked to sodium chromate-induced hepatotoxicity. GSE19662, a gene expression microarray, was obtained from Gene Expression Omnibus (GEO). Six primary rat hepatocyte (PRH) samples from GSE19662 include sodium chromate-treated (n = 3) and the control PRH samples (n = 3). A total of 2,525 differentially expressed genes (DEGs) were obtained, especially 962, and 1,563 genes were up- and downregulated in sodium chromate-treated PRHs compared to the control. Gene ontology (GO) enrichment analysis suggested that those DEGs were involved in multiple biological processes, including the response to toxic substances, the positive regulation of apoptotic process, lipid and cholesterol metabolic process, and others. Signaling pathway enrichment analysis indicated that the DEGs were mainly enriched in MAPK, PI3K-Akt, PPAR, AMPK, cellular senescence, hepatitis B, fatty acid biosynthesis, etc. Moreover, many genes, including CYP2E1, CYP1A2, CYP2C13, CDK1, NDC80, and CCNB1, might contribute to sodium chromate-induced hepatotoxicity. Taken together, this study enhances our knowledge of the potential molecular mechanisms of sodium chromate-induced hepatotoxicity.

Comprehensive analysis of BUBs gene family in lung adenocarcinoma with immunological analysis

Lung adenocarcinoma (LUAD) is one of the most commonly malignant tumors, and major challenges remain in the treatment of LUAD. Budding uninhibited by benzimidazole 1/3 (BUB1/3) play significant roles in the process of spindle-assembly checkpoint (SAC) during mitosis. However, their roles in LUAD have not been established. Here, we performed an immunological analysis of BUB1/3 in LUAD using a comprehensive bioinformatics approach, quantitative real-time-PCR and Western blotting technique. Our results indicated that the expression levels of BUB1 and BUB3 in LUAD samples were higher than the expression levels in the control groups and were associated with some clinicopathologic parameters in patients with LUAD. BUB1/3 and their related genes were enriched in cell immune, and the immune infiltration analysis revealed that the BUB1/3 expression profile was significantly correlated with characteristics of immune cell infiltration. Survival analysis showed that the disease-free survival and overall survival of patients with LUAD decreased with an increase in the BUB1/3 expression levels. The mRNA and protein expression levels of BUB1 and BUB3 in each of the LUAD cell lines were upregulated to varying degrees. BUB1 and BUB3 are the potential immunological therapeutic intervention targets for patients with LUAD.

Screening of Hub Genes in Hepatocellular Carcinoma Based on Network Analysis and Machine Learning

Hepatocellular carcinoma (LIHC) is the fifth common cancer worldwide, and it requires effective diagnosis and treatment to prevent aggressive metastasis. The purpose of this study was to construct a machine learning-based diagnostic model for the diagnosis of liver cancer. Using weighted correlation network analysis (WGCNA), univariate analysis, and Lasso-Cox regression analysis, protein-protein interactions network analysis is used to construct gene networks from transcriptome data of hepatocellular carcinoma patients and find hub genes for machine learning. The five models, including gradient boosting, random forest, support vector machine, logistic regression, and integrated learning, were to identify a multigene prediction model of patients. Immunological assessment, TP53 gene mutation and promoter methylation level analysis, and KEGG pathway analysis were performed on these groups. Potential drug molecular targets for the corresponding hepatocellular carcinomas were obtained by molecular docking for analysis, resulting in the screening of 2 modules that may be relevant to the survival of hepatocellular carcinoma patients, and the construction of 5 diagnostic models and multiple interaction networks. The modes of action of drug-molecule interactions that may be effective against hepatocellular carcinoma core genes CCNA2, CCNB1, and CDK1 were investigated. This study is expected to provide research ideas for early diagnosis of hepatocellular carcinoma.

Upregulation of circ_0008812 and circ_0001583 predicts poor prognosis and promotes breast cancer proliferation

Background: Accumulating evidence suggests that circular RNAs (circRNAs) are highly correlated with tumor progression and pathogenesis in breast cancer. Whereas, their regulatory roles and corresponding mechanisms in breast cancer are still not exhaustive. Thus, we intended to establish circRNA-mediated competive endogenous RNA (ceRNA) network to uncover the possible roles and clinical implications of circRNAs in breast cancer.

Methods: Microarray and RNA-sequencing (RNA-seq) data were download from GEO and TCGA database to screen for differentially expressed RNAs (DEcircRNAs, DEmiRNAs, DEmRNAs) in breast cancer. By implementing online databases, we established ceRNA networks, performed gene set enrichment analysis, constructed protein-protein interaction (PPI) networks, and assessed the expression levels and prognostic significance of hub genes. Subsequently, we explored the functions of prognosis-related genes and constructed gene-drug interaction networks. Finally, the functional roles of DEcircRNAs in breast cancer were revealed via MTT and colony formation assay.

Results: Based on the identified 8 DEcircRNAs, 25 miRNAs and 216 mRNAs, a ceRNA regulatory network was established. Further analysis revealed that prominent enrichments were transcription factor binding, transforming growth factor-beta (TGF-β) and Apelin signaling pathway etc. PPI network and survival curves analysis showed that elevated levels of hub genes (RACGAP1 and KPNA2) were associated with poorer prognosis. They were found to be positively relevant to cell cycle and proliferation. Then a prognostic sub-network of ceRNA was constructed, consisting of 2 circRNAs, 4 miRNAs and 2 mRNAs. The gene-drug interaction network showed that numerous drugs could regulate the expression of these two prognosis-related genes. Functional experiments showed that depletion of circ_0008812 and circ_0001583 could significantly inhibit the proliferation of MCF-7 cells.

Conclusion: Our study constructed 4 prognostic regulatory axes that are significantly correlated with tumor prognosis in breast cancer patients, and uncover the roles of circ_0008812 and circ_0001583 in breast cancer, providing a new perspective into the molecular mechanisms of breast cancer pathogenesis.

DNA Damage Repair Classifier Defines Distinct Groups in Hepatocellular Carcinoma

Simple Summary

DNA repair pathways have been implicated in hepatocellular carcinoma outcomes. We found that hepatocellular carcinomas (HCC) could be separated into two groups (high and low) based on the overall expression of genes involved in DNA repair. Among the low repair group, there were three subgroups, one of which shared features of the high repair group. Given the important role of liver in metabolism and detoxification and its regenerative capacity, proliferation and DNA damage responses are critical in subdividing major biological categories of liver tumors. High repair samples showed more proliferative and regenerative signatures and had poorer outcomes versus the low repair that were more associated with the genes involved in normal liver biology. These biological groups suggest that dysregulation in endogenous liver processes promotes a pro-tumorigenic microenvironment that may facilitate tumor progression or identify tumors that require more substantial clinical intervention.

Abstract

DNA repair pathways have been associated with variability in hepatocellular carcinoma (HCC) clinical outcomes, but the mechanism through which DNA repair varies as a function of liver regeneration and other HCC characteristics is poorly understood. We curated a panel of 199 genes representing 15 DNA repair pathways to identify DNA repair expression classes and evaluate their associations with liver features and clinicopathologic variables in The Cancer Genome Atlas (TCGA) HCC study. We identified two groups in HCC, defined by low or high expression across all DNA repair pathways. The low-repair group had lower grade and retained the expression of classical liver markers, whereas the high-repair group had more clinically aggressive features, increased p53 mutant-like gene expression, and high liver regenerative gene expression. These pronounced features overshadowed the variation in the low-repair subset, but when considered separately, the low-repair samples included three subgroups: L1, L2, and L3. L3 had high DNA repair expression with worse progression-free (HR 1.24, 95% CI 0.81–1.91) and overall (HR 1.63, 95% CI 0.98–2.71) survival. High-repair outcomes were also significantly worse compared with the L1 and L2 groups. HCCs vary in DNA repair expression, and a subset of tumors with high regeneration profoundly disrupts liver biology and poor prognosis.

The cell cycle, cancer development and therapy

The process of cell division plays a vital role in cancer progression. Cell proliferation and error-free chromosomes segregation during mitosis are central events in life cycle. Mistakes during cell division generate changes in chromosome content and alter the balances of chromosomes number. Any defects in expression of TIF1 family proteins, SAC proteins network, mitotic checkpoint proteins involved in chromosome mis-segregation and cancer development. Here we discuss the function of organelles deal with the chromosome segregation machinery, proteins and correction mechanisms involved in the accurate chromosome segregation during mitosis.

Molecular Regulation of Yak Preadipocyte Differentiation and Proliferation by LncFAM200B and ceRNA Regulatory Network Analysis

The positive regulatory role of lncFAM200B in differentiation and lipid deposition in yak intramuscular preadipocytes has been demonstrated in our previous study. However, the regulatory mechanisms remain unclear. In this study, we aimed to produce complete mRNA and microRNA (miRNA) profiles after adenovirus-mediated lncFAM200B overexpression in yak preadipocytes using high-throughput sequencing. We constructed a competing endogenous RNA (ceRNA) network with lncFAM200B as the core and identified the functions of the selected target miRNA during cell proliferation and differentiation. We obtained 118 differentially expressed genes (DEGs) after lncFAM200B overexpression, 76 of which were up-regulated, including Notch signaling members NOTCH3, DTX3L, and HES4, and 42 DEGs were down-regulated, including genes related to the cell cycle (CCNA2, BUB1, CDC20, TOP2A, and KIF20A). Additionally, many ubiquitin-mediated proteolysis pathway members were also significantly up-regulated (BUA7, PML, TRIM21, and TRIM25). MiRNA sequencing showed that 13 miRNAs were significantly up-regulated, and 12 miRNAs were down-regulated. Among them, 29 targets of 10 differentially expressed miRNAs (DEMs) were differentially expressed, including miR-152-FBXO33, miR-6529a-TRIM21, miR-148c-NOTCH3, and the miR-6529b-HES4 axis. We further verified that overexpression and inhibition of miR-6529a can inhibit and promote, respectively, the proliferation and differentiation of preadipocytes. Taken together, our study not only revealed the regulatory network of lncFAM200B during yak preadipocytes differentiation but also laid a foundation for elucidating the cause for lower intramuscular fat content in yaks at the molecular level.

Transcriptomic and Functional Evidence for Differential Effects of MCF-7 Breast Cancer Cell-Secretome on Vascular and Lymphatic Endothelial Cell Growth

Vascular and lymphatic vessels drive breast cancer (BC) growth and metastasis. We assessed the cell growth (proliferation, migration, and capillary formation), gene-, and protein-expression profiles of Vascular Endothelial Cells (VECs) and Lymphatic Endothelial Cells (LECs) exposed to a conditioned medium (CM) from estrogen receptor-positive BC cells (MCF-7) in the presence or absence of Estradiol. We demonstrated that MCF-7-CM stimulated growth and capillary formation in VECs but inhibited LEC growth. Consistently, MCF-7-CM induced ERK1/2 and Akt phosphorylation in VECs and inhibited them in LECs. Gene expression analysis revealed that the LECs were overall (≈10-fold) more sensitive to MCF-7-CM exposure than VECs. Growth/angiogenesis and cell cycle pathways were upregulated in VECs but downregulated in LECs. An angiogenesis proteome array confirmed the upregulation of 23 pro-angiogenesis proteins in VECs. In LECs, the expression of genes related to ATP synthesis and the ATP content were reduced by MCF-7-CM, whereas MTHFD2 gene, involved in folate metabolism and immune evasion, was upregulated. The contrasting effect of MCF-7-CM on the growth of VECs and LECs was reversed by inhibiting the TGF-β signaling pathway. The effect of MCF-7-CM on VEC growth was also reversed by inhibiting the VEGF signaling pathway. In conclusion, BC secretome may facilitate cancer cell survival and tumor growth by simultaneously promoting vascular angiogenesis and inhibiting lymphatic growth. The differential effects of BC secretome on LECs and VECs may be of pathophysiological relevance in BC.

Identification of Key Prognostic Genes of Triple Negative Breast Cancer by LASSO-Based Machine Learning and Bioinformatics Analysis

Improved insight into the molecular mechanisms of triple negative breast cancer (TNBC) is required to predict prognosis and develop a new therapeutic strategy for targeted genes. The aim of this study is to identify key genes which may affect the prognosis of TNBC patients by bioinformatic analysis. In our study, the RNA sequencing (RNA-seq) expression data of 116 breast cancer lacking ER, PR, and HER2 expression and 113 normal tissues were downloaded from The Cancer Genome Atlas (TCGA). We screened out 147 differentially co-expressed genes in TNBC compared to non-cancerous tissue samples by using weighted gene co-expression network analysis (WGCNA) and differential gene expression analysis. Then, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were constructed, revealing that 147 genes were mainly enriched in nuclear division, chromosomal region, ATPase activity, and cell cycle signaling. After using Cytoscape software for protein-protein interaction (PPI) network analysis and LASSO feature selection, a total of fifteen key genes were identified. Among them, BUB1 and CENPF were significantly correlated with the overall survival rate (OS) difference of TNBC patients (p value < 0.05). In addition, BUB1, CCNA2, and PACC1 showed significant poor disease-free survival (DFS) in TNBC patients (p value < 0.05), and may serve as candidate biomarkers in TNBC diagnosis. Thus, our results collectively suggest that BUB1, CCNA2, and PACC1 genes could play important roles in the progression of TNBC and provide attractive therapeutic targets.

BUBs Are New Biomarkers of Promoting Tumorigenesis and Affecting Prognosis in Breast Cancer

Background

A tumor occurs because of abnormal cell multiplication caused by many variables like a significant disturbance in the regulation of cell growth and the instability of chromosome mitosis. Budding uninhibited by benzimidazoles 1 (BUB1), BUB1 mitotic checkpoint serine/threonine kinase B (BUB1B), and budding uninhibited by benzimidazoles 3 (BUB3) are key regulators of mitosis, and their abnormal expression is highly correlated with breast cancer (BrCa), sarcoma, hepatic carcinoma, and other malignant tumors. However, the occurrence of BUBs (BUB1, BUB1B, and BUB3) and the development of BrCa have not been systematically explained.

Methods

Find out the target gene by looking up literature on PubMed and CNKI. Using the R software, TCGA, GEO, Kaplan-Meier Plotter, TIMER, and other databases, we studied the level of transcription, genetic changes, and physiological functions of BUBs in BrCa patients and their relationship with the origin, development, prognosis, immunity, and drug resistance of BrCa patients. Findings. We found that the high expression level of BUBs in BrCa tissues proposed a poor prognosis. The multivariate Cox regression analysis suggested that BUB1B and BUB3 might be independent prognostic factors of BrCa. In addition, the Metascape functional enrichment analysis showed that BUBs may be involved in the composition of the spindle, chromosome, and other structures and play a role in mitosis, sister chromatid separation, and other processes. Pathway enrichment suggests that BUBs may affect the cell cycle and lead to abnormal proliferation. Meanwhile, we also found that BUB3 can negatively regulate B lymphocytes, and BUB1 and BUB1B inhibit immune responses by promoting the secretion level of checkpoint molecules of the immune system, leading to immune escape of tumor cells.

Conclusion

We speculate that BUB1, BUB1B, and BUB3 may be therapeutic targets for BrCa patients and also provide new therapeutic strategies for BrCa treatment.

Bioinformatics Analysis Highlight Differentially Expressed CCNB1 and PLK1 Genes as Potential Anti-Breast Cancer Drug Targets and Prognostic Markers

Breast cancer is one of the most common malignant tumors in women worldwide. Early diagnosis, treatment, and prognosis of breast cancer are global challenges. Identification of valid predictive diagnosis and prognosis biomarkers and drug targets are crucial for breast cancer prevention. This study characterizes differentially expressed genes (DEGs) based on the TCGA database by using DESeq2, edgeR, and limma. A total of 2032 DEGs, including 1026 up-regulated genes and 1006 down-regulated genes were screened. Followed with WGCNA, PPI analysis, GEPIA 2, and HPA database verification, thirteen hub genes including CDK1, BUB1, BUB1B, CDC20, CCNB2, CCNB1, KIF2C, NDC80, CDCA8, CENPF, BIRC5, AURKB, PLK1, MAD2L1, and CENPE were obtained, and they may serve as potential therapeutic targets of breast cancer. Especially, overexpression of CCNB1 and PLK1 are strongly associated with the low survival rate of breast cancer patients, demonstrating their potentiality as prognostic markers. Moreover, CCNB1 and PLK1 are highly expressed in all breast cancer stages, suggesting that they could be further studied as potential drug targets. Taken together, our study highlights CCNB1 and PLK1 as potential anti-breast cancer drug targets and prognostic markers.

The Interplay of Cohesin and RNA Processing Factors: The Impact of Their Alterations on Genome Stability

Cohesin, a multi-subunit protein complex, plays important roles in sister chromatid cohesion, DNA replication, chromatin organization, gene expression, transcription regulation, and the recombination or repair of DNA damage. Recently, several studies suggested that the functions of cohesin rely not only on cohesin-related protein-protein interactions, their post-translational modifications or specific DNA modifications, but that some RNA processing factors also play an important role in the regulation of cohesin functions. Therefore, the mutations and changes in the expression of cohesin subunits or alterations in the interactions between cohesin and RNA processing factors have been shown to have an impact on cohesion, the fidelity of chromosome segregation and, ultimately, on genome stability. In this review, we provide an overview of the cohesin complex and its role in chromosome segregation, highlight the causes and consequences of mutations and changes in the expression of cohesin subunits, and discuss the RNA processing factors that participate in the regulation of the processes involved in chromosome segregation. Overall, an understanding of the molecular determinants of the interplay between cohesin and RNA processing factors might help us to better understand the molecular mechanisms ensuring the integrity of the genome.

Prognostic Value of BUB1 for Predicting Non-Muscle-Invasive Bladder Cancer Progression

Non-muscle-invasive bladder cancer (NMIBC) is a common disease with a high recurrence rate requiring lifetime surveillance. Although NMIBC is not life-threatening, it can progress to muscle-invasive bladder cancer (MIBC), a lethal form of the disease. The management of the two diseases differs, and patients with MIBC require aggressive treatments such as chemotherapy and radical cystectomy. NMIBC patients at a high risk of progression benefit from early immediate cystectomy. Thus, identifying concordant markers for accurate risk stratification is critical to predict the prognosis of NMIBC. Candidate genetic biomarkers associated with NMIBC prognosis were screened by RNA-sequencing of 24 tissue samples, including 16 NMIBC and eight normal controls, and by microarray analysis (GSE13507). Lastly, we selected and investigated a mitotic checkpoint serine/threonine kinase, BUB1, that regulates chromosome segregation during the cell cycle. BUB1 gene expression was tested in 86 NMIBC samples and 15 controls by real-time qPCR. The performance of BUB1 as a prognostic biomarker for NMIBC was validated in the internal Chungbuk cohort (GSE13507) and the external UROMOL cohort (E-MTAB-4321). BUB1 expression was higher in NMIBC patients than in normal controls (p < 0.05), and the overexpression of BUB1 was correlated with NMIBC progression (log-rank test, p = 0.007). In in vitro analyses, BUB1 promoted the proliferation of bladder cancer cells by accelerating the G2/M transition of the cell cycle. Conclusively, BUB1 modulates the G2/M transition to promote the proliferation of bladder cancer cells, suggesting that it could serve as a prognostic marker in NMIBC.

Key signaling networks are dysregulated in patients with the adipose tissue disorder, lipedema

Objectives

Lipedema, a poorly understood chronic disease of adipose hyper-deposition, is often mistaken for obesity and causes significant impairment to mobility and quality-of-life. To identify molecular mechanisms underpinning lipedema, we employed comprehensive omics-based comparative analyses of whole tissue, adipocyte precursors (adipose-derived stem cells (ADSCs)), and adipocytes from patients with or without lipedema.

Methods

We compared whole-tissues, ADSCs, and adipocytes from body mass index–matched lipedema (n = 14) and unaffected (n = 10) patients using comprehensive global lipidomic and metabolomic analyses, transcriptional profiling, and functional assays.

Results

Transcriptional profiling revealed >4400 significant differences in lipedema tissue, with altered levels of mRNAs involved in critical signaling and cell function-regulating pathways (e.g., lipid metabolism and cell-cycle/proliferation). Functional assays showed accelerated ADSC proliferation and differentiation in lipedema. Profiling lipedema adipocytes revealed >900 changes in lipid composition and >600 differentially altered metabolites. Transcriptional profiling of lipedema ADSCs and non-lipedema ADSCs revealed significant differential expression of >3400 genes including some involved in extracellular matrix and cell-cycle/proliferation signaling pathways. One upregulated gene in lipedema ADSCs, Bub1, encodes a cell-cycle regulator, central to the kinetochore complex, which regulates several histone proteins involved in cell proliferation. Downstream signaling analysis of lipedema ADSCs demonstrated enhanced activation of histone H2A, a key cell proliferation driver and Bub1 target. Critically, hyperproliferation exhibited by lipedema ADSCs was inhibited by the small molecule Bub1 inhibitor 2OH-BNPP1 and by CRISPR/Cas9-mediated Bub1 gene depletion.

Conclusion

We found significant differences in gene expression, and lipid and metabolite profiles, in tissue, ADSCs, and adipocytes from lipedema patients compared to non-affected controls. Functional assays demonstrated that dysregulated Bub1 signaling drives increased proliferation of lipedema ADSCs, suggesting a potential mechanism for enhanced adipogenesis in lipedema. Importantly, our characterization of signaling networks driving lipedema identifies potential molecular targets, including Bub1, for novel lipedema therapeutics.

Hub Genes and Key Pathways of Intervertebral Disc Degeneration: Bioinformatics Analysis and Validation

Objective

To identify significant pathways and genes in intervertebral disc degeneration (IDD) based on bioinformatics analysis.

Design

The GEO database was used to download the GSE124272 dataset. Differentially expressed genes (DEGs) were analyzed using Limma package in R language. Then, gene ontologies (GO), Kyoto encyclopedia of genes and genomes (KEGG), and protein-protein interaction (PPI) networks were used to further identify hub genes. The mRNA expression levels of top six hub genes were verified.

Results

We found 563 DEGs, of which 214 were upregulated and 349 were downregulated. The top 5 GO terms and pathways were shown including immune response, cell cycle, and p53 pathway. Based on the PPI analysis, we verified the mRNA expression levels of 6 hub genes. The mRNA levels of CHEK1, CDCA2, SKA3, and KIF20A were upregulated in degenerative NP tissue than in healthy NP tissue. However, the mRNA level of BUB1 and SPC25 was downregulated.

Conclusions

This study may provide new biomarkers for the IDD and treatments to repair IDD related to CHEK1, CDCA2, SKA3, BUB1, KIF20A, and SPC25.

Dysregulation of KIF14 regulates the cell cycle and predicts poor prognosis in cervical cancer: a study based on integrated approaches

Cervical cancer (CC) is the most common malignant tumor in females. Although persistent high-risk human papillomavirus (HPV) infection is a leading factor that causes CC, few women with HPV infection develop CC. Therefore, many mechanisms remain to be explored, such as aberrant expression of oncogenes and tumor suppressor genes. To identify promising prognostic factors and interpret the relevant mechanisms of CC, the RNA sequencing profile of CC was downloaded from the Cancer Genome Atlas and the Gene Expression Omnibus databases. The GSE63514 dataset was analyzed, and differentially expressed genes (DEGs) were obtained by weighted coexpression network analysis and the edgeR package in R. Fifty-three shared genes were mainly enriched in nuclear chromosome segregation and DNA replication signaling pathways. Through a protein-protein interaction network and prognosis analysis, the kinesin family member 14 (KIF14) hub gene was extracted from the set of 53 shared genes, which was overexpressed and associated with poor overall survival (OS) and disease-free survival (DFS) of CC patients. Mechanistically, gene set enrichment analysis showed that KIF14 was mainly enriched in the glycolysis/gluconeogenesis signaling pathway and DNA replication signaling pathway, especially in the cell cycle signaling pathway. RT-PCR and the Human Protein Atlas database confirmed that these genes were significantly increased in CC samples. Therefore, our findings indicated the biological function of KIF14 in cervical cancer and provided new ideas for CC diagnosis and therapies.

In Silico Methods for the Identification of Diagnostic and Favorable Prognostic Markers in Acute Myeloid Leukemia

Acute myeloid leukemia (AML), the most common type of acute leukemia in adults, is mainly asymptomatic at early stages and progresses/recurs rapidly and frequently. These attributes necessitate the identification of biomarkers for timely diagnosis and accurate prognosis. In this study, differential gene expression analysis was performed on large-scale transcriptomics data of AML patients versus corresponding normal tissue. Weighted gene co-expression network analysis was conducted to construct networks of co-expressed genes, and detect gene modules. Finally, hub genes were identified from selected modules by applying network-based methods. This robust and integrative bioinformatics approach revealed a set of twenty-four genes, mainly related to cell cycle and immune response, the diagnostic significance of which was subsequently compared against two independent gene expression datasets. Furthermore, based on a recent notion suggesting that molecular characteristics of a few, unusual patients with exceptionally favorable survival can provide insights for improving the outcome of individuals with more typical disease trajectories, we defined groups of long-term survivors in AML patient cohorts and compared their transcriptomes versus the general population to infer favorable prognostic signatures. These findings could have potential applications in the clinical setting, in particular, in diagnosis and prognosis of AML.

Gene Expression Profiling of Olfactory Neuroblastoma Helps Identify Prognostic Pathways and Define Potentially Therapeutic Targets

Simple Summary

The gene expression profile of ONB defines a group of patients with a dismal prognosis and identifies potentially targetable pathways. Better prognostic stratification may offer new tailored approaches for the treatment and follow-up of ONB. The integration of new therapeutic agents with standard surgical and RT strategies may improve the outcomes in cases with worse prognoses. Furthermore, the ontogenesis of ONB in basal and neural subtypes is mirrored by different transcriptional pathways, paving the way towards different therapeutic approaches.

Abstract

Olfactory neuroblastoma (ONB) is a rare sinonasal neoplasm with a peculiar behavior, for which limited prognostic factors are available. Herein, we investigate the transcriptional pathways altered in ONB and correlate them with pathological features and clinical outcomes. We analyze 32 ONB patients treated with curative intent at two independent institutions from 2001 to 2019 for whom there is available pathologic and clinical data. We perform gene expression profiling on primary ONB samples and carry out functional enrichment analysis to investigate the key pathways associated with disease-free survival (DFS). The median age is 53.5 years; all patients undergo surgery and a pure endoscopic approach is adopted in the majority of cases (81.2%). Most patients have advanced disease (stages III–IV, 81.2%) and 84.4% undergo adjuvant (chemo)radiotherapy. The median follow-up is 35 months; 11 (26.8%) patients relapse. Clinical characteristics (gender, stage and Hyams’ grade) are not associated with the outcomes. In contrast, TGF-beta binding, EMT, IFN-alpha response, angiogenesis, IL2-STAT5 and IL6-JAK-STAT3 signaling pathways are enriched in patients experiencing recurrence, and significantly associated with shorter DFS. Clustering of transcriptional profiles according to pathological features indicates two distinct molecular groups, defined by either cytokeratin-positive or -negative immunostaining. Definition of the characterizing ONB transcriptomic pathways may pave the way towards tailored treatment approaches.

Expression and prognosis analyses of BUB1, BUB1B and BUB3 in human sarcoma

Budding Uninhibited By Benzimidazoles are a group of genes encoding proteins that play central roles in spindle checkpoint during mitosis. Improper mitosis may lead to aneuploidy which is found in many types of tumors. As a key mediator in mitosis, the dysregulated expression of BUBs has been proven to be highly associated with various malignancies, such as leukemia, gastric cancer, breast cancer, and liver cancer. However, bioinformatic analysis has not been applied to explore the role of the BUBs in sarcomas. Herein, we investigate the transcriptional and survival data of BUBs in patients with sarcomas using Oncomine, Gene Expression Profiling Interactive Analysis, Cancer Cell Line Encyclopedia, Kaplan-Meier Plotter, LinkedOmics, and the Database for Annotation, Visualization and Integrated Discovery. We found that the expression levels of BUB1, BUB1B and BUB3 were higher in sarcoma samples and cell lines than in normal controls. Survival analysis revealed that the higher expression levels of BUB1, BUB1B and BUB3 were associated with lower overall and disease-free survival in patients with sarcomas. This study implies that BUB1, BUB1B and BUB3 are potential treatment targets for patients with sarcomas and are new biomarkers for the prognosis of sarcomas.

Bioinformatics analysis identified hub genes in prostate cancer tumorigenesis and metastasis

OBJECTIVE

Prostate cancer (PCa) is the most frequent cancer found in males worldwide, and its mortality rate is increasing every year. To discover key molecular changes in PCa development and metastasis, we analyzed microarray data of localized PCa, metastatic PCa and normal prostate tissue samples from clinical specimens.

METHODS

Gene expression profiling datasets of PCa were analyzed online. The DAVID was used to perform GO functional and KEGG pathway enrichment analyses. CytoHubba in Cytoscape software was applied to identify hub genes. Survival data were downloaded from GEPIA. Gene expression data were obtained from ONCOMINE and UALCAN.

RESULTS

We obtained 4 sets of differentially expressed genes (DEGs), DEGs 1: a comparison of the gene expression between 4 normal prostate and 5 localized PCa samples in GSE27616; DEGs 2: a comparison of the gene expression between 6 normal prostate and 7 localized PCa samples in GSE3325; DEGs 3: a comparison of the gene expression between 5 localized PCa and 4 metastatic PCa samples in GSE27616; DEGs 4: a comparison of the gene expression between 7 localized PCa and 6 metastatic PCa samples in GSE3325. A comparison of these 4 sets of genes revealed 51 overlapped genes. GO function analysis revealed enrichment of the 51 DEGs in functions related to the proteinaceous extracellular matrix and centrosome, protein homodimerization activity and chromatin binding were the main functions of these genes, which participated in regulating cell division, mitotic nuclear division, proteinaceous extracellular matrix, cell adhesion and apoptotic process. KEGG pathway analysis indicated that these identified DEGs were mainly enriched in progesterone-mediated oocyte maturation, oocyte meiosis and cell cycle. We defined the 16 genes with the highest degree of connectivity as the hub genes in the 51 overlapped DEGs. Cox regression revealed TOP2A, CCNB2, BUB1, CDK1 and EZH2 were related to Disease-free survival (DFS). The expression levels of the 5 genes were 2.232-, 1.786-, 2.303-, 1.699-, and 1.986-fold higher in PCa than the levels in normal tissues, respectively (P < 0.05). We obtained 20 hub genes from DEGs by the comparison of normal prostate tissue vs. localized cancer tissue. Among them, KIF20A, CDKN3, PBK and CDCA2, were expressed higher in PCa than in normal tissues, and were associated with the DFS of PCa patients. Meanwhile, we obtained 20 hub genes from DEGs by the comparison of localized cancer tissue vs. metastatic cancer tissue. Cox regression revealed PLK1, CCNA2 and CDC20, were associated with both the DFS and overall survival of PCa patients.

CONCLUSIONS

The results suggested that the functions of KIF20A, CDKN3, PBK and CDCA2 may contribute to PCa development and the functions of PLK1, CCNA2 and CDC20 may contribute to PCa metastasis. Meanwhile, the functions of TOP2A, CCNB2, BUB1, CDK1 and EZH2 may contribute to both PCa development and metastasis.

Investigating the molecular control of deer antler extract on articular cartilage

BACKGROUND

Deer antler is considered as a precious traditional Chinese medicinal material and has been widely used to reinforce kidney's yang, nourish essence, and strengthen bone function. The most prominent bioactive components in deer antler are water-soluble proteins that play potential roles in bone formation and repair. The aim of this study was to explore the molecular control and therapeutic targets of deer antler extract (DAE) on articular cartilage.

METHODS

DAE was prepared as previously described. All rats were randomly divided into Blank group and DAE group (10 rats per group) after 7-day adaptive feeding. The rats in DAE group were orally administrated with DAE at a dose of 0.2 g/kg per day for 3 weeks, and the rats in Blank group were fed with drinking water. Total RNA was isolated from the articular cartilage of knee joints. RNA sequencing (RNA-seq) experiment combined with quantitative real-time polymerase chain reaction (qRT-PCR) verification assay was carried out to explore the molecular control and therapeutic targets of DAE on articular cartilage.

RESULTS

We demonstrated that DAE significantly increased the expression levels of functional genes involved in cartilage formation, growth, and repair and decreased the expression levels of susceptibility genes involved in the pathophysiology of osteoarthritis.

CONCLUSIONS

DAE might serve as a candidate supplement for maintaining cartilage homeostasis and preventing cartilage degeneration and inflammation. These effects were possibly achieved by accelerating the expression of functional genes involved in chondrocyte commitment, survival, proliferation, and differentiation and suppressing the expression of susceptibility genes involved in the pathophysiology of osteoarthritis. Thus, our findings will contribute towards deepening the knowledge about the molecular control and therapeutic targets of DAE on the treatment of cartilage-related diseases.

Identification and Analysis of Potential Key Genes Associated With Hepatocellular Carcinoma Based on Integrated Bioinformatics Methods

BACKGROUND

Hepatocellular carcinoma (HCC) is a type of primary liver tumor with poor prognosis and high mortality, and its molecular mechanism remains incompletely understood. This study aimed to use bioinformatics technology to identify differentially expressed genes (DEGs) in HCC pathogenesis, hoping to identify novel biomarkers or potential therapeutic targets for HCC research.

METHODS

The bioinformatics analysis of our research mostly involved the following two datasets: Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA). First, we screened DEGs based on the R packages (limma and edgeR). Using the DAVID database, the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of DEGs were carried out. Next, the protein-protein interaction (PPI) network of the DEGs was built in the STRING database. Then, hub genes were screened through the cytoHubba plug-in, followed by verification using the GEPIA and Oncomine databases. We demonstrated differences in levels of the protein in hub genes using the Human Protein Atlas (HPA) database. Finally, the hub genes prognostic values were analyzed by the GEPIA database. Additionally, using the Comparative Toxicogenomics Database (CTD), we constructed the drug-gene interaction network.

RESULTS

We ended up with 763 DEGs, including 247 upregulated and 516 downregulated DEGs, that were mainly enriched in the epoxygenase P450 pathway, oxidation-reduction process, and metabolism-related pathways. Through the constructed PPI network, it can be concluded that the P53 signaling pathway and the cell cycle are the most obvious in module analysis. From the PPI, we filtered out eight hub genes, and these genes were significantly upregulated in HCC samples, findings consistent with the expression validation results. Additionally, survival analysis showed that high level gene expression of CDC20, CDK1, MAD2L1, BUB1, BUB1B, CCNB1, and CCNA2 were connected with the poor overall survival of HCC patients. Toxicogenomics analysis showed that only topotecan, oxaliplatin, and azathioprine could reduce the gene expression levels of all seven hub genes.

CONCLUSION

The present study screened out the key genes and pathways that were related to HCC pathogenesis, which could provide new insight for the future molecularly targeted therapy and prognosis evaluation of HCC.