Super-Enhancer-Associated Hub Genes In Chronic Myeloid Leukemia Identified Using Weighted Gene Co-Expression Network Analysis

KIF2C: An important factor involved in signaling pathways, immune infiltration, and DNA damage repair in tumorigenesis

BACKGROUNDS

Poorly regulated mitosis and chromosomal instability are common characteristics in malignant tumor cells. Kinesin family member 2 C (KIF2C), also known as mitotic centromere-associated kinesin (MCAK) is an essential component during mitotic regulation. In recent years, KIF2C was shown to be dysregulated in several tumors and was involved in many aspects of tumor self-regulation. Research on KIF2C may be a new direction and target for anti-tumor therapy.

OBJECT

The article aims at reviewing current literatures and summarizing the research status of KIF2C in malignant tumors as well as the oncogenic signaling pathways associated with KIF2C and its role in immune infiltration.

RESULT

In this review, we summarize the KIF2C mechanisms and signaling pathways in different malignant tumors, and briefly describe its involvement in pathways related to classical chemotherapeutic drug resistance, such as MEK/ERK, mTOR, Wnt/β-catenin, P53 and TGF-β1/Smad pathways. KIF2C upregulation was shown to promote tumor cell migration, invasion, chemotherapy resistance and inhibit DNA damage repair. It was also highly correlated with microRNAs, and CD4 +T cell and CD8 +T cell tumor immune infiltration.

CONCLUSION

This review shows that KIF2C may function as a new anticancer drug target with great potential for malignant tumor treatment and the mitigation of chemotherapy resistance.

Identification of module genes and functional pathway analysis in septic shock subtypes by integrated bioinformatics analysis

BACKGROUND

The present study aimed to identify the module genes and key gene functions and biological pathways of septic shock (SS) through integrated bioinformatics analysis.

METHODS

In the study, we performed batch correction and principal component analysis on 282 SS samples and 79 normal control samples in three datasets, GSE26440, GSE95233 and GSE57065, to obtain a combined corrected gene expression matrix containing 21,654 transcripts. Patients with SS were then divided into three molecular subtypes according to sample subtyping analysis.

RESULTS

By analyzing the demographic characteristics of the different subtypes, we found no statistically significant differences in gender ratio and age composition among the three groups. Then, three subtypes of differentially expressed genes (DEGs) and specific upregulated DEGs (SDEGs) were identified by differential gene expression analysis. We found 7361 DEGs in the type I group, 5594 DEGs in the type II group, and 7159 DEGs in the type III group. There were 1698 SDEGs in the type I group, 2443 in the type II group, and 1831 in the type III group. In addition, we analyzed the correlation between the expression data of 5972 SDEGs in the three subtypes and the gender and age of 227 patients, constructed a weighted gene co-expression network, and identified 11 gene modules, among which the module with the highest correlation with gender ratio was MEgrey. The modules with the highest correlation with age composition were MEgrey60 and MElightyellow. Then, by analyzing the differences in module genes among different subgroups of SS, we obtained the differential expression of 11 module genes in four groups: type I, type II, type III and the control group. Finally, we analyzed the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment of all module DEGs, and the GO function and KEGG pathway enrichment of different module genes were different.

CONCLUSIONS

Our findings aim to identify the specific genes and intrinsic molecular functional pathways of SS subtypes, as well as further explore the genetic and molecular pathophysiological mechanisms of SS.

Advances in lung adenocarcinoma: A novel perspective on prognoses and immune responses of CENPO as an oncogenic superenhancer

Lung adenocarcinoma (LUAD) is the most prevalent form of lung cancer globally, and its treatment remains a significant challenge. Therefore, it is crucial to comprehend the microenvironment to improve therapy and prognosis urgently. In this study, we utilized bioinformatic methods to analyze the transcription expression profile of patient samples with complete clinical information from the TCGA-LUAD datasets. To validate our findings, we also analyzed the Gene Expression Omnibus (GEO) datasets. The super-enhancer (SE) was visualized using the peaks of the H3K27ac and H3K4me1 ChIP-seq signal, which were identified by the Integrative Genomics Viewer (IGV). To further investigate the role of Centromere protein O (CENPO) in LUAD, we conducted various assays including Western blot, qRT-PCR, flow cytometry, wound healing and transwell assays to assess the cell functions of CENPO in vitro. The overexpression of CENPO is linked to a poor prognosis in patients with LUAD. Strong signal peaks of H3K27ac and H3K4me1 were also observed near the predicted SE regions of CENPO. CENPO was found to be positively associated with the expression levels of immune checkpoints and drug IC50 value (Roscovitine and TGX221), but negatively associated with the fraction levels of several immature cells and drug IC50 value (CCT018159, GSK1904529A, Lenaildomide, and PD-173074). Additionally, CENPO-associated prognostic signature (CPS) was identified as an independent risk factor. The high-risk group for LUAD is identified based on CPS enrichment, which involved not only endocytosis that transfers mitochondria to promote cell survival in response to chemotherapy but also cell cycle promotion that leads to drug resistance. The removal of CENPO significantly suppressed metastasis and induced arrest and apoptosis of LUAD cells. The involvement of CENPO in the immunosuppression of LUAD provides a prognostic signature for LUAD patients.

Super-Enhancer–Associated nine-gene prognostic score model for prediction of survival in chronic lymphocytic leukemia patients

Chronic lymphocytic leukemia (CLL) is a type of highly heterogeneous mature B-cell malignancy with various disease courses. Although a multitude of prognostic markers in CLL have been reported, insights into the role of super-enhancer (SE)–related risk indicators in the occurrence and development of CLL are still lacking. A super-enhancer (SE) is a cluster of enhancers involved in cell differentiation and tumorigenesis, and is one of the promising therapeutic targets for cancer therapy in recent years. In our study, the CLL-related super-enhancers in the training database were processed by LASSO-penalized Cox regression analysis to screen a nine-gene prognostic model including TCF7, VEGFA, MNT, GMIP, SLAMF1, TNFRSF25, GRWD1, SLC6AC, and LAG3. The SE-related risk score was further constructed and it was found that the predictive performance with overall survival and time-to-treatment (TTT) was satisfactory. Moreover, a high correlation was found between the risk score and already known prognostic markers of CLL. In the meantime, we noticed that the expressions of TCF7, GMIP, SLAMF1, TNFRSF25, and LAG3 in CLL were different from those of healthy donors (p < 0.01). Moreover, the risk score and LAG3 level of matched pairs before and after treatment samples varied significantly. Finally, an interactive nomogram consisting of the nine-gene risk group and four clinical traits was established. The inhibitors of mTOR and cyclin-dependent kinases (CDKs) were considered effective in patients in the high-risk group according to the pRRophetic algorithm. Collectively, the SE-associated nine-gene prognostic model developed here may be used to predict the prognosis and assist in the risk stratification and treatment of CLL patients in the future.

CENPO is Associated with Immune Cell Infiltration and is a Potential Diagnostic and Prognostic Marker for Hepatocellular Carcinoma

Purpose

To examine the expression, clinical significance, and potential regulatory mechanism of centromere protein O (CENPO) in hepatocellular carcinoma (HCC).

Methods

CENPO expression in pan-cancer was studied using the TCGA-GTEx database, in HCC and normal liver tissues using the GEO and TCGA databases, and in clinical HCC samples by RT-qPCR. The diagnostic value of CENPO was assessed using receiver operating characteristic curves. Univariate and multivariate regression analyses of factors associated with HCC prognosis were performed. CENPO function and its mechanism in HCC were explored using GO, KEGG, and GSEA analyses of differentially expressed genes (DEGs). Association of CENPO expression with immune cell infiltration and immune checkpoint-associated molecules was conducted using TCGA data and the TIMER2.0 database. Relationships between CENPO expression and DNA methylation were analyzed using the UALCAN and cBioPortal databases. CENPO expression in HCC cell lines was detected using RT-qPCR.

Results

CENPO is upregulated in most cancers, including HCC and cell lines, and is a potential biomarker for HCC diagnosis (AUC = 0.936, 95% CI: 0.911–0.960). Higher CENPO expression was associated with poorer outcomes in patients with HCC (OS, p = 0.004; DSS, p = 0.002; PFI, p < 0.001), and CENPO was an independent predictor of factors influencing overall survival in HCC. DEGs between samples with high and low CENPO levels were enriched in various biological processes, including activation of the G2M checkpoint and other signaling pathways, while CENPO expression correlated with HCC immune cell infiltration and immune checkpoint-associated molecules, as well as CENPO promoter methylation (p < 0.001).

Conclusion

In HCC and cell lines, CENPO is overexpressed, a potential diagnostic marker and an indicator of poor prognosis. CENPO may regulate HCC development by influencing nuclear division and tumor immune infiltration and is regulated by methylation, making it a potential target for HCC immunotherapy.

Super enhancers as master gene regulators in the pathogenesis of hematologic malignancies

Transcriptional deregulation of multiple oncogenes, tumor suppressors and survival pathways is a cancer cell hallmark. Super enhancers (SE) are long stretches of active enhancers in close linear proximity that ensure extraordinarily high expression levels of key genes associated with cell lineage, function and survival. SE landscape is intrinsically prone to changes and reorganization during the course of normal cell differentiation. This functional plasticity is typically utilized by cancer cells, which remodel their SE landscapes to ensure oncogenic transcriptional reprogramming. Multiple recent studies highlighted structural genetic mechanisms in non-coding regions that create new SE or hijack already existing ones. In addition, alterations in abundance/activity of certain SE-associated proteins or certain viral infections can elicit new super enhancers and trigger SE-driven transcriptional changes. For these reasons, SE profiling emerged as a powerful tool for discovering the core transcriptional regulatory circuits in tumor cells. This, in turn, provides new insights into cancer cell biology, and identifies main nodes of key cellular pathways to be potentially targeted. Since SEs are susceptible to inhibition, their disruption results in exponentially amassing 'butterfly' effect on gene expression and cell function. Moreover, many of SE elements are druggable, opening new therapeutic opportunities. Indeed, SE targeting drugs have been studied preclinically in various hematologic malignancies with promising effects. Herein, we review the unique features of SEs, present different cis- and trans-acting mechanisms through which hematologic tumor cells acquire SEs, and finally, discuss the potential of SE targeting in the therapy of hematologic malignancies.

Network analysis with biological data of cancer patients: A scoping review

BACKGROUND & OBJECTIVE

Network Analysis (NA) is a mathematical method that allows exploring relations between units and representing them as a graph. Although NA was initially related to social sciences, the past two decades was introduced in Bioinformatics. The recent growth of the networks' use in biological data analysis reveals the need to further investigate this area. In this work, we attempt to identify the use of NA with biological data, and specifically: (a) what types of data are used and whether they are integrated or not, (b) what is the purpose of this analysis, predictive or descriptive, and (c) the outcome of such analyses, specifically in cancer diseases.

METHODS & MATERIALS

The literature review was conducted on two databases, PubMed & IEEE, and was restricted to journal articles of the last decade (January 2010 - December 2019). At a first level, all articles were screened by title and abstract, and at a second level the screening was conducted by reading the full text article, following the predefined inclusion & exclusion criteria leading to 131 articles of interest. A table was created with the information of interest and was used for the classification of the articles. The articles were initially classified to analysis studies and studies that propose a new algorithm or methodology. Each one of these categories was further screened by the following clustering criteria: (a) data used, (b) study purpose, (c) study outcome. Specifically for the studies proposing a new algorithm, the novelty presented in each one was detected.

RESULTS

& Conclusions: In the past five years researchers are focusing on creating new algorithms and methodologies to enhance this field. The articles' classification revealed that only 25% of the analyses are integrating multi-omics data, although 50% of the new algorithms developed follow this integrative direction. Moreover, only 20% of the analyses and 10% of the newly developed methodologies have a predictive purpose. Regarding the result of the works reviewed, 75% of the studies focus on identifying, prognostic or not, gene signatures. Concluding, this review revealed the need for deploying predictive and multi-omics integrative algorithms and methodologies that can be used to enhance cancer diagnosis, prognosis and treatment.

An interactive nomogram based on clinical and molecular signatures to predict prognosis in multiple myeloma patients

Although novel drugs and treatments have been developed and improved, multiple myeloma (MM) is still recurrent and difficult to cure. In the present study, the magenta module containing 400 hub genes was determined from the training dataset of GSE24080 through weighted gene co-expression network analysis (WGCNA). Then, using the least absolute shrinkage and selection operator (Lasso) analysis, a fifteen-gene signature was firstly selected and the predictive performance for overall survival (OS) was favorable, which was identified by Receiver Operating Characteristic (ROC) curves. The risk score model was constructed based on survival-associated fifteen genes from the Lasso model, which classified MM patients into high-risk and low-risk groups. Areas under the curve (AUC) of ROC curve and log-rank test showed that the high-risk group was correlated to the dismal survival outcome of MM patients, which was also identified in testing dataset of GSE9782. The calibration plot, the AUC value of the ROC curve and Concordance-index showed that the interactive nomogram with risk score could favorably predict the probability of multi-year OS of MM patients. Therefore, it may help clinicians make a precise therapeutic decision based on the easy-to-use tool of the nomogram.

Construction of a Five-Super-Enhancer-Associated-Genes Prognostic Model for Osteosarcoma Patients

Osteosarcoma is a malignant tumor most commonly arising in children and adolescents and associated with poor prognosis. In recent years, some prognostic models have been constructed to assist clinicians in the treatment of osteosarcoma. However, the prognosis and treatment of patients with osteosarcoma remain unsatisfactory. Notably, super-enhancer (SE)-associated genes strongly promote the progression of osteosarcoma. In the present study, we constructed a novel effective prognostic model using SE-associated genes from osteosarcoma. Five SE-associated genes were initially screened through the least absolute shrinkage and selection operator (Lasso) penalized Cox regression, as well as univariate and multivariate Cox regression analyses. Meanwhile, a risk score model was constructed using the expression of these five genes. The excellent performance of the five-SE-associated-gene-based prognostic model was determined via time-dependent receiver operating characteristic (ROC) curves and Kaplan-Meier curves. Inferior outcome of overall survival (OS) was predicted in the high-risk group. A nomogram based on the polygenic risk score model was further established to validate the performance of the prognostic model. It showed that our prognostic model performed outstandingly in predicting 1-, 3-, and 5-year OS of patients with osteosarcoma. Meanwhile, these five genes also belonged to the hub genes associated with survival and necrosis of osteosarcoma according to the result of weighted gene co-expression network analysis based on the dataset of GSE39058. Therefore, we believe that the five-SE-associated-gene-based prognostic model established in this study can accurately predict the prognosis of patients with osteosarcoma and effectively assist clinicians in treating osteosarcoma in the future.

Identification of hub methylated-CpG sites and associated genes in oral squamous cell carcinoma

To improve personalized diagnosis and prognosis for oral squamous cell carcinoma (OSCC) by identification of hub methylated‐CpG sites and associated genes, weighted gene comethylation network analysis (WGCNA) was performed to examine and identify hub modules and CpG sites correlated with OSCC. Here, WGCNA modeling yielded blue and brown comethylation modules that were significantly associated with OSCC status. Following screening of the differentially expressed genes (DEGs) from gene expression microarrays and differentially methylated‐CpG sites (DCGs), integrated multiomics analysis of the DEGs, DCGs, and hub CpG sites from the modules was performed to investigate their correlations. Expression levels of 16 CpG sites‐associated genes were negatively correlated with methylation patterns of promoter. Moreover, Kaplan‐Meier survival analysis of the hub CpG sites and associated genes was carried out using 2 public databases, MethSurv and GEPIA. Only 5 genes, ACTA1, ACTN2, OSR1, SYNGR1, and ZNF677, had significant overall survival using GEPIA. Hypermethylated‐CpG sites ACTN2‐cg21376883 and OSR1‐cg06509239 were found to be associated with poor survival by MethSurv. Methylation status of specific site and expression levels of associated genes were determined using clinical samples by quantitative methylation‐specific PCR and real‐time PCR. Pearson's correlation analysis showed that methylation levels of cg06509239 and cg18335068 were negatively related to OSR1 and ZNF677 expression levels, respectively. Our classification schema using multiomics analysis represents a screening framework for identification of hub CpG sites and associated genes.