Analysis of gene co‑expression network reveals prognostic significance of CNFN in patients with head and neck cancer

Identification of exosome-related gene signature as a promising diagnostic and therapeutic tool for breast cancer

Background

Exosomes are promising tools for the development of new diagnostic and therapeutic approaches. Exosomes possess the ability to activate signaling pathways that contribute to the remodeling of the tumor microenvironment, angiogenesis, and the regulation of immune responses. We aimed to develop a prognostic score based on exosomes derived from breast cancer.

Materials and methods

Training was conducted on the TCGA-BRCA dataset, while validation was conducted on GSE20685, GSE5764, GSE7904, and GSE29431. A total of 121 genes related to exosomes were retrieved from the ExoBCD database. The Cox proportional hazards model is used to develop risk score model. The GSVA package was utilized to analyze single-sample gene sets and identify exosome signatures, while the WGCNA package was utilized to identify gene modules associated with clinical outcomes. The clusterProfiler and GSVA R packages facilitated gene set enrichment and variation analyses. Furthermore, CIBERSORT quantified immune infiltration, and a correlation between gene expression and drug sensitivity was assessed using the TIDE algorithm.

Results

An exosome-related prognostic score was established using the following selected genes: ABCC9, PIGR, CXCL13, DOK7, CD24, and IVL. Various immune cells that promote cancer immune evasion were associated with a high-risk prognostic score, which was an independent predictor of outcome. High-risk and low-risk groups exhibited significantly different infiltration abundances (p < 0.05). By conducting a sensitivity comparison, we found that patients with high-risk scores exhibited more favorable responses to immunotherapy than those with low-risk scores.

Conclusion

The exosome-related gene signature exhibits outstanding performance in predicting the prognosis and cancer status of patients with breast cancer and guiding immunotherapy.

Identification and validation of novel prognostic signatures based on m5C methylation patterns and tumor EMT profiles in head and neck squamous cell carcinoma

The role of 5-methylcytosine (m5C) in tumor initiation and progression has been increasingly recognized. However, the precise association between the regulation of m5C and the progression, metastasis, and prognosis of head and neck squamous cell carcinoma (HNSCC) has not yet been fully explored. Data from 545 HNSCC patients obtained from The Cancer Genome Atlas (TCGA) database were analyzed. Unsupervised cluster analysis was conducted using the expression levels of m5C regulatory genes. Additionally, gene set variation analysis (GSVA), single-sample gene set enrichment analysis (ssGSEA), and Cox regression analysis were utilized. Quantitative reverse transcription polymerase chain reaction (RT-qPCR), colony formation assay, transwell experiments and western blots were performed in the HNSCC cell line UM-SCC-17B to assess the expression and functional role of one of the novel signatures, CNFN. Significant expression differences were found in m5C regulatory genes between tumor and normal tissues in HNSCC. Two distinct m5C modification patterns, characterized by substantial prognostic differences, were identified. Cluster-2, which exhibited a strong association with epithelial-mesenchymal transition (EMT), was found to be associated with a poorer prognosis. Based on the m5C clusters and EMT status, differentially expressed genes (DEGs) were identified. Using DEGs, an 8-gene signature (CAMK2N1, WNT7A, F2RL1, AREG, DEFB1, CNFN, TGFBI, and CAV1) was established to develop a prognostic model. The performance of this signature was validated in both the training and external validation datasets, demonstrating its promising efficacy. Furthermore, additional investigations using RT-qPCR on clinical specimens and experimental assays in cell lines provided compelling evidence suggesting that CNFN, one of the genes in the signature, could play a role in HNSCC progression and metastasis through the EMT pathway. This study highlighted the role of m5C in HNSCC progression and metastasis. The relationship between m5C and EMT has been elucidated for the first time. A robust prognostic model was developed for accurately predicting HNSCC patients' survival outcomes. Potential molecular mechanisms underlying these associations have been illuminated through this research.

FOXA1 repression drives lineage plasticity and immune heterogeneity in bladder cancers with squamous differentiation

Cancers arising from the bladder urothelium often exhibit lineage plasticity with regions of urothelial carcinoma adjacent to or admixed with regions of divergent histomorphology, most commonly squamous differentiation. To define the biologic basis for and clinical significance of this morphologic heterogeneity, here we perform integrated genomic analyses of mixed histology bladder cancers with separable regions of urothelial and squamous differentiation. We find that squamous differentiation is a marker of intratumoral genomic and immunologic heterogeneity in patients with bladder cancer and a biomarker of intrinsic immunotherapy resistance. Phylogenetic analysis confirms that in all cases the urothelial and squamous regions are derived from a common shared precursor. Despite the presence of marked genomic heterogeneity between co-existent urothelial and squamous differentiated regions, no recurrent genomic alteration exclusive to the urothelial or squamous morphologies is identified. Rather, lineage plasticity in bladder cancers with squamous differentiation is associated with loss of expression of FOXA1, GATA3, and PPARG, transcription factors critical for maintenance of urothelial cell identity. Of clinical significance, lineage plasticity and PD-L1 expression is coordinately dysregulated via FOXA1, with patients exhibiting morphologic heterogeneity pre-treatment significantly less likely to respond to immune checkpoint inhibitors.

Suprabasin: Role in human cancers and other diseases

Suprabasin (SBSN), a gene with unknown function located in q13 region of chromosome 19, was first found to be expressed in the basal layer of the stratified epithelium in mouse and human tissues and was thought to be a potential precursor of keratinized capsules. However, in recent years, significant progress has been made in the study of SBSN in a variety of human diseases. One common theme appears to be the effect of SBSN on tumor progression, such as invasion, metastasis and resistance. However, the function and mechanism of action of SBSN is still elusive. In this study, we reviewed the literature on SBSN in the PubMed database to identify the basic characteristics, biological functions, and roles of SBSN in cancer and other diseases. In particular, we focused on the potential mechanisms of SBSN activity, to improve our understanding of the complex function of this protein and provide a theoretical basis for further research on the role of SBSN in cancer and other diseases.

Transcriptome research identifies four hub genes related to primary myelofibrosis: a holistic research by weighted gene co-expression network analysis

Objectives: This study aimed to identify specific diagnostic as well as predictive targets of primary myelofibrosis (PMF).

Methods: The gene expression profiles of GSE26049 were obtained from Gene Expression Omnibus (GEO) dataset, WGCNA was constructed to identify the most related module of PMF. Subsequently, Gene Ontology (GO), Kyoto Encyclopedia Genes and Genomes (KEGG), Gene Set Enrichment Analysis (GSEA) and Protein-Protein interaction (PPI) network were conducted to fully understand the detailed information of the interested green module. Machine learning, Principal component analysis (PCA), and expression pattern analysis including immunohistochemistry and immunofluorescence of genes and proteins were performed to validate the reliability of these hub genes.

Results: Green module was strongly correlated with PMF disease after WGCNA analysis. 20 genes in green module were identified as hub genes responsible for the progression of PMF. GO, KEGG revealed that these hub genes were primarily enriched in erythrocyte differentiation, transcription factor binding, hemoglobin complex, transcription factor complex and cell cycle, etc. Among them, EPB42, CALR, SLC4A1 and MPL had the most correlations with PMF. Machine learning, Principal component analysis (PCA), and expression pattern analysis proved the results in this study.

Conclusions: EPB42, CALR, SLC4A1 and MPL were significantly highly expressed in PMF samples. These four genes may be considered as candidate prognostic biomarkers and potential therapeutic targets for early stage of PMF. The effects are worth expected whether in the diagnosis at early stage or as therapeutic target.

Paradoxical Roles of Desmosomal Components in Head and Neck Cancer

Desmosomes are intercellular adhesion complexes involved in various aspects of epithelial pathophysiology, including tissue homeostasis, morphogenesis, and disease development. Recent studies have reported that the abnormal expression of various desmosomal components correlates with tumor progression and poor survival. In addition, desmosomes have been shown to act as a signaling platform to regulate the proliferation, invasion, migration, morphogenesis, and apoptosis of cancer cells. The occurrence and progression of head and neck cancer (HNC) is accompanied by abnormal expression of desmosomal components and loss of desmosome structure. However, the role of desmosomal components in the progression of HNC remains controversial. This review aims to provide an overview of recent developments showing the paradoxical roles of desmosomal components in tumor suppression and promotion. It offers valuable insights for HNC diagnosis and therapeutics development.

A ceRNA-associated risk model predicts the poor prognosis for head and neck squamous cell carcinoma patients

Head and neck squamous cell carcinoma (HNSCC) is one of the most malignant cancers with poor prognosis worldwide. Emerging evidence indicates that competing endogenous RNAs (ceRNAs) are involved in various diseases, however, the regulatory mechanisms of ceRNAs underlying HNSCC remain unclear. In this study, we retrieved differentially expressed long non-coding RNAs (DElncRNAs), messenger RNAs (DEmRNAs) and microRANs (DEmiRNAs) from The Cancer Genome Atlas database and constructed a ceRNA-based risk model in HNSCC by integrated bioinformatics approaches. Functional enrichment analyses showed that DEmRNAs might be involved in extracellular matrix related biological processes, and protein–protein interaction network further selected out prognostic genes, including MYL1 and ACTN2. Importantly, co-expressed RNAs identified by weighted co-expression gene network analysis constructed the ceRNA networks. Moreover, AC114730.3, AC136375.3, LAT and RYR3 were highly correlated to overall survival of HNSCC by Kaplan–Meier method and univariate Cox regression analysis, which were subsequently implemented multivariate Cox regression analysis to build the risk model. Our study provides a deeper understanding of ceRNAs on the regulatory mechanisms, which will facilitate the expansion of the roles on the ceRNAs in the tumorigenesis, development and treatment of HNSCC.

KNSTRN promotes tumorigenesis and gemcitabine resistance by activating AKT in bladder cancer

KNSTRN is a component of the mitotic spindle, which was rarely investigated in tumorigenesis. AKT plays an essential role in tumorigenesis by modulating the phosphorylation of various substrates. The activation of AKT is regulated by PTEN and PIP₃. Here, we prove KNSTRN is positively correlated with malignancy of bladder cancer and KNSTRN activates AKT phosphorylation at Thr308 and Ser473. More importantly, our study reveals that both KNSTRN and PTEN interact with PH domain of AKT at cell membrane. The amount of KNSTRN interacted with AKT is negatively related to PTEN. Furthermore, PIP₃ pull-down assay proves that KNSTRN promoted AKT movement to PIP₃. These data suggest KNSTRN may activate AKT phosphorylation by promoting AKT movement to PIP₃ and alleviating PTEN suppression. Based on the activation of AKT phosphorylation, our study demonstrates that KNSTRN promotes bladder cancer metastasis and gemcitabine resistance in vitro and in vivo. Meanwhile, the effect of KNSTRN on tumorigenesis and gemcitabine resistance could be restored by AKT specific inhibitor MK2206 or AKT overexpression. In conclusion, we identify an oncogene KNSTRN that promotes tumorigenesis and gemcitabine resistance by activating AKT phosphorylation and may serve as a therapeutic target in bladder cancer.

Suprabasin-A Review

Among the ~22,000 human genes, very few remain that have unknown functions. One such example is suprabasin (SBSN). Originally described as a component of the cornified envelope, the function of stratified epithelia-expressed SBSN is unknown. Both the lack of knowledge about the gene role under physiological conditions and the emerging link of SBSN to various human diseases, including cancer, attract research interest. The association of SBSN expression with poor prognosis of patients suffering from oesophageal carcinoma, glioblastoma multiforme, and myelodysplastic syndromes suggests that SBSN may play a role in human tumourigenesis. Three SBSN isoforms code for the secreted proteins with putative function as signalling molecules, yet with poorly described effects. In this first review about SBSN, we summarised the current knowledge accumulated since its original description, and we discuss the potential mechanisms and roles of SBSN in both physiology and pathology.

WFDC2 suppresses prostate cancer metastasis by modulating EGFR signaling inactivation

WAP four-disulfide core domain 2 (WFDC2) is a small secretory protein that has been widely studied in ovarian cancer. It has been proven that WFDC2 promotes proliferation and metastasis in ovarian cancer, and serves as a diagnostic biomarker. However, the specific function of WFDC2 in prostate cancer has not been reported. Here, we first screened the diagnostic marker and favorable prognostic factor WFDC2 in prostate cancer by bioinformatics. WFDC2 expression was negatively correlated with Gleason score and metastasis in prostate cancer. Then, we revealed that overexpression of WFDC2, and addition of recombinant protein HE4 can significantly inhibit prostate cancer metastasis in vivo and in vitro. By co-immunoprecipitation and co-localization assays, we proved that WFDC2 binds to the extracellular domain of epidermal growth factor receptor (EGFR). Immunoblot showed that WFDC2 overexpression and recombinant protein HE4 addition inactivated the EGFR/AKT/GSK3B/Snail signaling pathway, and then restrained the progression of epithelial-mesenchymal transition. In conclusion, our study identified that the tumor suppressor WFDC2 can suppress prostate cancer metastasis by inactivating EGFR signaling.

Identification of Hub Genes as Biomarkers Correlated with the Proliferation and Prognosis in Lung Cancer: A Weighted Gene Co-Expression Network Analysis

Lung cancer is one of the most malignant tumors in the world. Early diagnosis and treatment of lung cancer are vitally important to reduce the mortality of lung cancer patients. In the present study, we attempt to identify the candidate biomarkers for lung cancer by weighted gene co-expression network analysis (WGCNA). Gene expression profile of GSE30219 was downloaded from the gene expression omnibus (GEO) database. The differentially expressed genes (DEGs) were analyzed by the limma package, and the co-expression modules of genes were built by WGCNA. UALCAN was used to analyze the relative expression of normal group and tumor subgroups based on tumor individual cancer stages. Survival analysis for the hub genes was performed by Kaplan–Meier plotter analysis with the TCGA database. A total of 2176 genes (745 upregulated and 1431 downregulated genes) were obtained from the GSE30219 database. Seven gene co-expression modules were conducted by WGCNA and the blue module might be inferred as the most crucial module in the pathogenesis of lung cancer. In the pathway enrichment analysis of KEGG, the candidate genes were enriched in the “DNA replication,” “Cell cycle,” and “P53 signaling pathway” pathways. Among these, the cell cycle pathway was the most significant pathway in the blue module with four hub genes CCNB1, CCNE2, MCM7, and PCNA which were selected in our study. Kaplan–Meier plotter analysis indicated that the high expressions of four hub genes were correlated with a worse overall survival (OS) and advanced tumors. qRT-PCR showed that mRNA expression levels of MCM7 (p = 0.038) and CCNE2 (0.003) were significantly higher in patients with the TNM stage. In summary, the high expression of the MCM7 and CCNE2 were significantly related with advanced tumors and worse OS in lung cancer. Thus, the MCM7 and CCNE2 genes can be good indicators for cellular proliferation and prognosis in lung cancer.

Identification and validation of key modules and hub genes associated with the pathological stage of oral squamous cell carcinoma by weighted gene co-expression network analysis

Background

Oral squamous cell carcinoma (OSCC) is a major lethal malignant cancer of the head and neck region, yet its molecular mechanisms of tumourigenesis are still unclear.

Patients and methods

We performed weighted gene co-expression network analysis (WGCNA) on RNA-sequencing data with clinical information obtained from The Cancer Genome Atlas (TCGA) database. The relationship between co-expression modules and clinical traits was investigated by Pearson correlation analysis. Furthermore, the prognostic value and expression level of the hub genes of these modules were validated based on data from the TCGA database and other independent datasets from the Gene Expression Omnibus (GEO) database and the Human Protein Atlas database. The significant modules and hub genes were also assessed by functional analysis and gene set enrichment analysis (GSEA).

Results

We found that the turquoise module was strongly correlated with pathologic T stage and significantly enriched in critical functions and pathways related to tumourigenesis. PPP1R12B, CFD, CRYAB, FAM189A2 and ANGPTL1 were identified and statistically validated as hub genes in the turquoise module and were closely implicated in the prognosis of OSCC. GSEA indicated that five hub genes were significantly involved in many well-known cancer-related biological functions and signaling pathways.

Conclusion

In brief, we systematically discovered a co-expressed turquoise module and five hub genes associated with the pathologic T stage for the first time, which provided further insight that WGCNA may reveal the molecular regulatory mechanism involved in the carcinogenesis and progression of OSCC. In addition, the five hub genes may be considered candidate prognostic biomarkers and potential therapeutic targets for the precise early diagnosis, clinical treatment and prognosis of OSCC in the future.

An outcome model for human bladder cancer: A comprehensive study based on weighted gene co-expression network analysis

The precision evaluation of prognosis is crucial for clinical treatment decision of bladder cancer (BCa). Therefore, establishing an effective prognostic model for BCa has significant clinical implications. We performed WGCNA and DEG screening to initially identify the candidate genes. The candidate genes were applied to construct a LASSO Cox regression analysis model. The effectiveness and accuracy of the prognostic model were tested by internal/external validation and pan‐cancer validation and time‐dependent ROC. Additionally, a nomogram based on the parameter selected from univariate and multivariate cox regression analysis was constructed. Eight genes were eventually screened out as progression‐related differentially expressed candidates in BCa. LASSO Cox regression analysis identified 3 genes to build up the outcome model in E‐MTAB‐4321 and the outcome model had good performance in predicting patient progress free survival of BCa patients in discovery and test set. Subsequently, another three datasets also have a good predictive value for BCa patients' OS and DFS. Time‐dependent ROC indicated an ideal predictive accuracy of the outcome model. Meanwhile, the nomogram showed a good performance and clinical utility. In addition, the prognostic model also exhibits good performance in pan‐cancer patients. Our outcome model was the first prognosis model for human bladder cancer progression prediction via integrative bioinformatics analysis, which may aid in clinical decision‐making.

A Chemical Probe for Tudor Domain Protein Spindlin1 to Investigate Chromatin Function

Modifications of histone tails, including lysine/arginine methylation, provide the basis of a "chromatin or histone code". Proteins that contain "reader" domains can bind to these modifications and form specific effector complexes, which ultimately mediate chromatin function. The spindlin1 (SPIN1) protein contains three Tudor methyllysine/arginine reader domains and was identified as a putative oncogene and transcriptional coactivator. Here we report a SPIN1 chemical probe inhibitor with low nanomolar in vitro activity, exquisite selectivity on a panel of methyl reader and writer proteins, and with submicromolar cellular activity. X-ray crystallography showed that this Tudor domain chemical probe simultaneously engages Tudor domains 1 and 2 via a bidentate binding mode. Small molecule inhibition and siRNA knockdown of SPIN1, as well as chemoproteomic studies, identified genes which are transcriptionally regulated by SPIN1 in squamous cell carcinoma and suggest that SPIN1 may have a role in cancer related inflammation and/or cancer metastasis.