Genomic Instability Profiles at the Single Cell Level in Mouse Colorectal Cancers of Defined Genotypes

Integrating bulk and single-cell RNA sequencing data reveals epithelial-mesenchymal transition molecular subtype and signature to predict prognosis, immunotherapy efficacy, and drug candidates in low-grade gliomas

Objective: Epithelial-mesenchymal transition (EMT) is a tightly regulated and dynamic process occurring in both embryonic development and tumor progression. Our study aimed to comprehensively explore the molecular subtypes, immune landscape, and prognostic signature based on EMT-related genes in low-grade gliomas (LGG) in order to facilitate treatment decision-making and drug discovery. Methods: We curated EMT-related genes and performed molecular subtyping with consensus clustering algorithm to determine EMT expression patterns in LGG. The infiltration level of diverse immune cell subsets was evaluated by implementing the single-sample gene set enrichment analysis (ssGSEA) and ESTIMATE algorithms. The distinctions in clinical characteristics, mutation landscape, and immune tumor microenvironment (TME) among the subtypes were subjected to further investigation. Gene Set Variation Analysis (GSVA) was performed to explore the biological pathways that were involved in subtypes. The chemo drug sensitivity and immunotherapy of subtypes were estimated through GDSC database and NTP algorithm. To detect EMT subtype-related prognostic gene modules, the analysis of weighted gene co-expression network (WGCNA) was performed. The LASSO algorithm was utilized to construct a prognostic risk model, and its efficacy was verified through an independent CGGA dataset. Finally, the expression of the hub genes from the prognostic model was evaluated through the single-cell dataset and in-vitro experiment. Results: The TCGA-LGG dataset revealed the creation of two molecular subtypes that presented different prognoses, clinical implications, TME, mutation landscapes, chemotherapy, and immunotherapy. A three-gene signature (SLC39A1, CTSA and CLIC1) based on EMT expression pattern were established through WGCNA analysis. Low-risk patients showed a positive outlook, increased immune cell presence, and higher expression of immune checkpoint proteins. In addition, several promising drugs, including birinapant, fluvastatin, clofarabine, dasatinib, tanespimycin, TAK-733, GDC-0152, AZD8330, trametinib and ingenol-mebutate had great potential to the treatment of high risk patients. Finally, CTSA and CLIC1 were highly expressed in monocyte cell through single-cell RNA sequencing analysis. Conclusion: Our research revealed non-negligible role of EMT in the TME diversity and complexity of LGG. A prognostic signature may contribute to the personalized treatment and prognostic determination.

Prognostic prediction and immunotherapy response analysis of the fatty acid metabolism-related genes in clear cell renal cell carcinoma

Background

Clear cell renal cell carcinoma (ccRCC) is a common urinary cancer. Although diagnostic and therapeutic approaches for ccRCC have been improved, the survival outcomes of patients with advanced ccRCC remain unsatisfactory. Fatty acid metabolism (FAM) has been increasingly recognized as a critical modulator of cancer development. However, the significance of the FAM in ccRCC remains unclear. Herein, we explored the function of a FAM-related risk score in the stratification and prediction of treatment responses in patients with ccRCC.

Methods

First, we applied an unsupervised clustering method to categorize patients from The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) datasets into subtypes and retrieved FAM-related genes from the MSigDB database. We discern differentially expressed genes (DEGs) among different subtypes. Then, we applied univariate Cox regression analysis followed by least absolute shrinkage and selection operator (LASSO) linear regression based on DEGs expression to establish a FAM-related risk score for ccRCC.

Results

We stratified the three ccRCC subtypes based on FAM-related genes with distinct overall survival (OS), clinical features, immune infiltration patterns, and treatment sensitivities. We screened nine genes from the FAM-related DEGs in the three subtypes to establish a risk prediction model for ccRCC. Nine FAM-related genes were differentially expressed in the ccRCC cell line ACHN compared to the normal kidney cell line HK2. High-risk patients had worse OS, higher genomic heterogeneity, a more complex tumor microenvironment (TME), and elevated expression of immune checkpoints. This phenomenon was validated in the ICGC cohort.

Conclusion

We constructed a FAM-related risk score that predicts the prognosis and therapeutic response of ccRCC. The close association between FAM and ccRCC progression lays a foundation for further exploring FAM-related functions in ccRCC.

Comprehensive analysis of cellular senescence-related genes in the prognosis, tumor microenvironment, and immunotherapy/chemotherapy of clear cell renal cell carcinoma

Background

The transcriptome public database and advances in biological discoveries contributed to significant progresses in identifying the drivers of cancer progression. Cellular senescence (CS) is considered as a leading factor resulting in cancer development. The purpose of this study was to explore the significance of CS-related genes in the molecular classification and survival outcome of clear cell renal cell carcinoma (ccRCC).

Methods

CS-related genes were obtained from the CellAge database, and patients from TCGA-KIRC dataset and ICGC dataset were clustered by ConsesusClusterPlus. The characteristics of overall survival (OS), genomic variation, and tumor microenvironment (TME) of each cluster were analyzed. Least Absolute Shrinkage and Selection Operator (LASSO) Cox regression analysis was conducted to develop a CS-related risk model to score ccRCC patients and assess the risk scores in predicting patients’ response to immunotherapy and chemotherapy. A nomogram based on the risk model was established to improve the risk stratification of patients.

Results

CcRCC was divided into three molecular subtypes based on CS-related genes. The three molecular phenotypes showed different OS and clinical manifestations, mutation patterns, and TME states. Five genes were obtained from nine differentially expressed CS-related genes in the three molecular subtypes to develop a risk model. Patients with ccRCC were divided into high- and low-risk subgroups. The former showed an unfavorable OS, with a significantly higher genomic variation rate, TME score, and numerous immune checkpoint expressions when compared to the low-risk subgroup. Risk score reflected the response of patients to axitinib, bortezomib, sorafenib, sunitinib, and temsirolimus.

Conclusions

In general, CS-related genes divided ccRCC into three molecular subtypes with distinct OS, mutation patterns, and TME states. The risk model based on the five CS-related genes can predict the prognosis and therapeutic outcome of ccRCC patients, providing a theoretical basis for further study on the molecular mechanism of CS-related ccRCC.

Primary and metastatic peritoneal surface malignancies

Peritoneal surface malignancies comprise a heterogeneous group of primary tumours, including peritoneal mesothelioma, and peritoneal metastases of other tumours, including ovarian, gastric, colorectal, appendicular or pancreatic cancers. The pathophysiology of peritoneal malignancy is complex and not fully understood. The two main hypotheses are the transformation of mesothelial cells (peritoneal primary tumour) and shedding of cells from a primary tumour with implantation of cells in the peritoneal cavity (peritoneal metastasis). Diagnosis is challenging and often requires modern imaging and interventional techniques, including surgical exploration. In the past decade, new treatments and multimodal strategies helped to improve patient survival and quality of life and the premise that peritoneal malignancies are fatal diseases has been dismissed as management strategies, including complete cytoreductive surgery embedded in perioperative systemic chemotherapy, can provide cure in selected patients. Furthermore, intraperitoneal chemotherapy has become an important part of combination treatments. Improving locoregional treatment delivery to enhance penetration to tumour nodules and reduce systemic uptake is one of the most active research areas. The current main challenges involve not only offering the best treatment option and developing intraperitoneal therapies that are equivalent to current systemic therapies but also defining the optimal treatment sequence according to primary tumour, disease extent and patient preferences. New imaging modalities, less invasive surgery, nanomedicines and targeted therapies are the basis for a new era of intraperitoneal therapy and are beginning to show encouraging outcomes.

Identification of a Somatic Mutation-Derived Long Non-Coding RNA Signatures of Genomic Instability in Renal Cell Carcinoma

Background

Renal cell carcinoma (RCC) is a malignant tumor with high morbidity and mortality. It is characterized by a large number of somatic mutations and genomic instability. Long non-coding RNAs (lncRNAs) are widely involved in the expression of genomic instability in renal cell carcinoma. But no studies have identified the genome instability-related lncRNAs (GInLncRNAs) and their clinical significances in RCC.

Methods

Clinical data, gene expression data and mutation data of 943 RCC patients were downloaded from The Cancer Genome Atlas (TCGA) database. Based on the mutation data and lncRNA expression data, GInLncRNAs were screened out. Co-expression analysis, Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis were conducted to explore their potential functions and related signaling pathways. A prognosis model was further constructed based on genome instability-related lncRNAs signature (GInLncSig). And the efficiency of the model was verified by receiver operating characteristic (ROC) curve. The relationships between the model and clinical information, prognosis, mutation number and gene expression were analyzed using correlation prognostic analysis. Finally, the prognostic model was verified in clinical stratification according to TCGA dataset.

Results

A total of 45 GInLncRNAs were screened out. Functional analysis showed that the functional genes of these GInLncRNAs were mainly enriched in chromosome and nucleoplasmic components, DNA binding in molecular function, transcription and complex anabolism in biological processes. Univariate and Multivariate Cox analyses further screened out 11 GInLncSig to construct a prognostic model (AL031123.1, AC114803.1, AC103563.7, AL031710.1, LINC00460, AC156455.1, AC015977.2, 'PRDM16-dt', AL139351.1, AL035661.1 and LINC01606), and the coefficient of each GInLncSig in the model was calculated. The area under the curve (AUC) value of the ROC curve was 0.770. Independent analysis of the model showed that the GInLncSig model was significantly correlated with the RCC patients' overall survival. Furthermore, the GInLncSig model still had prognostic value in different subgroups of RCC patients.

Conclusion

Our study preliminarily explored the relationship between genomic instability, lncRNA and clinical characteristics of RCC patients, and constructed a GInLncSig model consisted of 11 GInLncSig to predict the prognosis of patients with RCC. At the same time, our study provided theoretical support for the exploration of the formation and development of RCC.