Computational solutions for spatial transcriptomics

Transcriptome level expression data connected to the spatial organization of the cells and molecules would allow a comprehensive understanding of how gene expression is connected to the structure and function in the biological systems. The spatial transcriptomics platforms may soon provide such information. However, the current platforms still lack spatial resolution, capture only a fraction of the transcriptome heterogeneity, or lack the throughput for large scale studies. The strengths and weaknesses in current ST platforms and computational solutions need to be taken into account when planning spatial transcriptomics studies. The basis of the computational ST analysis is the solutions developed for single-cell RNA-sequencing data, with advancements taking into account the spatial connectedness of the transcriptomes. The scRNA-seq tools are modified for spatial transcriptomics or new solutions like deep learning-based joint analysis of expression, spatial, and image data are developed to extract biological information in the spatially resolved transcriptomes. The computational ST analysis can reveal remarkable biological insights into spatial patterns of gene expression, cell signaling, and cell type variations in connection with cell type-specific signaling and organization in complex tissues. This review covers the topics that help choosing the platform and computational solutions for spatial transcriptomics research. We focus on the currently available ST methods and platforms and their strengths and limitations. Of the computational solutions, we provide an overview of the analysis steps and tools used in the ST data analysis. The compatibility with the data types and the tools provided by the current ST analysis frameworks are summarized.

Identifies microtubule-binding protein CSPP1 as a novel cancer biomarker associated with ferroptosis and tumor microenvironment

Centrosome and spindle pole-associated protein (CSPP1) is a centrosome and microtubule-binding protein that plays a role in cell cycle-dependent cytoskeleton organization and cilia formation. Previous studies have suggested that CSPP1 plays a role in tumorigenesis; however, no pan-cancer analysis has been performed. This study systematically investigates the expression of CSPP1 and its potential clinical outcomes associated with diagnosis, prognosis, and therapy. CSPP1 is widely present in tissues and cells and its aberrant expression serves as a diagnostic biomarker for cancer. CSPP1 dysregulation is driven by multi-dimensional mechanisms involving genetic alterations, DNA methylation, and miRNAs. Phosphorylation of CSPP1 at specific sites may play a role in tumorigenesis. In addition, CSPP1 correlates with clinical features and outcomes in multiple cancers. Take brain low-grade gliomas (LGG) with a poor prognosis as an example, functional enrichment analysis implies that CSPP1 may play a role in ferroptosis and tumor microenvironment (TME), including regulating epithelial-mesenchymal transition, stromal response, and immune response. Further analysis confirms that CSPP1 dysregulates ferroptosis in LGG and other cancers, making it possible for ferroptosis-based drugs to be used in the treatment of these cancers. Importantly, CSPP1-associated tumors are infiltrated in different TMEs, rendering immune checkpoint blockade therapy beneficial for these cancer patients. Our study is the first to demonstrate that CSPP1 is a potential diagnostic and prognostic biomarker associated with ferroptosis and TME, providing a new target for drug therapy and immunotherapy in specific cancers.

Discovery of tumor immune infiltration-related snoRNAs for predicting tumor immune microenvironment status and prognosis in lung adenocarcinoma

Lung adenocarcinoma (LUAD) has a high mortality rate and is difficult to diagnose and treat in its early stage. Previous studies have demonstrated that small nucleolar RNAs (snoRNAs) play a critical role in tumor immune infiltration and the development of a variety of solid tumors. However, there have been no studies on the correlation between tumor-infiltrating immune-related snoRNAs (TIISRs) and LUAD. In this study, we filtered six immune-related snoRNAs based on the tissue specificity index (TSI) and expression profile of all snoRNAs between all LUAD cell lines from the Cancer Cell Line Encyclopedia and 21 types of immune cells from the Gene Expression Omnibus database. Further, we performed real-time quantitative polymerase chain reaction (RT-qPCR) to validate the expression status of these snoRNAs on peripheral blood mononuclear cells (PBMCs) and lung cancer cell lines. Next, we developed a TIISR signature based on the expression profiles of snoRNAs from 479 LUAD patients filtered by the random survival forest algorithm. We then analyzed the value of this TIISR signature (TIISR risk score) for assessing tumor immune infiltration, immune checkpoint inhibitor (ICI) treatment response, and the prognosis of LUAD between groups with high and low TIISR risk score. Further, we found that the TIISR risk score groups showed significant differences in biological characteristics and that the risk score could be used to assess the level of tumor immune cell infiltration, thereby predicting prognosis and responsiveness to immunotherapy in LUAD patients.

The molecular feature of macrophages in tumor immune microenvironment of glioma patients

Background

Gliomas are one of the most common types of primary tumors in central nervous system. Previous studies have found that macrophages actively participate in tumor growth.

Methods

Weighted gene co-expression network analysis was used to identify meaningful macrophage-related gene genes for clustering. Pamr, SVM, and neural network were applied for validating clustering results. Somatic mutation and methylation were used for defining the features of identified clusters. Differentially expressed genes (DEGs) between the stratified groups after performing elastic regression and principal component analyses were used for the construction of MScores. The expression of macrophage-specific genes were evaluated in tumor microenvironment based on single cell sequencing analysis. A total of 2365 samples from 15 glioma datasets and 5842 pan-cancer samples were used for external validation of MScore.

Results

Macrophages were identified to be negatively associated with the survival of glioma patients. Twenty-six macrophage-specific DEGs obtained by elastic regression and PCA were highly expressed in macrophages at single-cell level. The prognostic value of MScores in glioma was validated by the active proinflammatory and metabolic profile of infiltrating microenvironment and response to immunotherapies of samples with this signature. MScores managed to stratify patient survival probabilities in 15 external glioma datasets and pan-cancer datasets, which predicted worse survival outcome. Sequencing data and immunohistochemistry of Xiangya glioma cohort confirmed the prognostic value of MScores. A prognostic model based on MScores demonstrated high accuracy rate.

Conclusion

Our findings strongly support a modulatory role of macrophages, especially M2 macrophages in glioma progression and warrants further experimental studies.