Spatial transcriptomics reveals niche-specific enrichment and vulnerabilities of radial glial stem-like cells in malignant gliomas

Nanotechnology as a new strategy for the diagnosis and treatment of gliomas

Glioma is the most common malignant tumor of the central nervous system (CNS), and is characterized by high aggressiveness and a high recurrence rate. Currently, the main treatments for gliomas include surgical resection, temozolomide chemotherapy and radiotherapy. However, the prognosis of glioma patients after active standardized treatment is still poor, especially for glioblastoma (GBM); the median survival is still only 14.6 months, and the 5-year survival rate is only 4-5%. The current challenges in glioma treatment include difficulty in complete surgical resection, poor blood‒brain barrier (BBB) drug permeability, therapeutic resistance, and difficulty in tumor-specific targeting. In recent years, the rapid development of nanotechnology has provided new directions for diagnosing and treating gliomas. Nanoparticles (NPs) are characterized by excellent surface tunability, precise targeting, excellent biocompatibility, and high safety. In addition, NPs can be used for gene therapy, photodynamic therapy, and antiangiogenic therapy and can be combined with biomaterials for thermotherapy. In recent decades, breakthroughs in diagnosing and treating gliomas have been made with various functional NPs, and NPs are expected to become a new strategy for glioma diagnosis and treatment. In this paper, we review the main obstacles in the treatment of glioma and discuss the potential and challenges of the latest nanotechnology in the diagnosis and treatment of glioma.

The CEBPB+ glioblastoma subcluster specifically drives the formation of M2 tumor-associated macrophages to promote malignancy growth

Rationale: The heterogeneity of tumor cells within the glioblastoma (GBM) microenvironment presents a complex challenge in curbing GBM progression. Understanding the specific mechanisms of interaction between different GBM cell subclusters and non-tumor cells is crucial. Methods: In this study, we utilized a comprehensive approach integrating glioma single-cell and spatial transcriptomics. This allowed us to examine the molecular interactions and spatial localization within GBM, focusing on a specific tumor cell subcluster, GBM subcluster 6, and M2-type tumor-associated macrophages (M2 TAMs). Results: Our analysis revealed a significant correlation between a specific tumor cell subcluster, GBM cluster 6, and M2-type TAMs. Further in vitro and in vivo experiments demonstrated the specific regulatory role of the CEBPB transcriptional network in GBM subcluster 6, which governs its tumorigenicity, recruitment of M2 TAMs, and polarization. This regulation involves molecules such as MCP1 for macrophage recruitment and the SPP1-Integrin αvβ1-Akt signaling pathway for M2 polarization. Conclusion: Our findings not only deepen our understanding of the formation of M2 TAMs, particularly highlighting the differential roles played by heterogeneous cells within GBM in this process, but also provided new insights for effectively controlling the malignant progression of GBM.

Glioblastoma microenvironment-from biology to therapy

Glioblastoma (GBM) is the most aggressive primary brain cancer. These tumors exhibit high intertumoral and intratumoral heterogeneity in neoplastic and nonneoplastic compartments, low lymphocyte infiltration, and high abundance of myeloid subsets that together create a highly protumorigenic immunosuppressive microenvironment. Moreover, heterogeneous GBM cells infiltrate adjacent brain tissue, remodeling the neural microenvironment to foster tumor electrochemical coupling with neurons and metabolic coupling with nonneoplastic astrocytes, thereby driving growth. Here, we review heterogeneity in the GBM microenvironment and its role in low-to-high-grade glioma transition, concluding with a discussion of the challenges of therapeutically targeting the tumor microenvironment and outlining future research opportunities.

STimage-1K4M: A histopathology image-gene expression dataset for spatial transcriptomics

Recent advances in multi-modal algorithms have driven and been driven by the increasing availability of large image-text datasets, leading to significant strides in various fields, including computational pathology. However, in most existing medical image-text datasets, the text typically provides high-level summaries that may not sufficiently describe sub-tile regions within a large pathology image. For example, an image might cover an extensive tissue area containing cancerous and healthy regions, but the accompanying text might only specify that this image is a cancer slide, lacking the nuanced details needed for in-depth analysis. In this study, we introduce STimage-1K4M, a novel dataset designed to bridge this gap by providing genomic features for sub-tile images. STimage-1K4M contains 1,149 images derived from spatial transcriptomics data, which captures gene expression information at the level of individual spatial spots within a pathology image. Specifically, each image in the dataset is broken down into smaller sub-image tiles, with each tile paired with 15,000 - 30,000 dimensional gene expressions. With 4,293,195 pairs of sub-tile images and gene expressions, STimage-1K4M offers unprecedented granularity, paving the way for a wide range of advanced research in multi-modal data analysis an innovative applications in computational pathology, and beyond.

Systematic evaluation with practical guidelines for single-cell and spatially resolved transcriptomics data simulation under multiple scenarios

BACKGROUND

Single-cell RNA sequencing (scRNA-seq) and spatially resolved transcriptomics (SRT) have led to groundbreaking advancements in life sciences. To develop bioinformatics tools for scRNA-seq and SRT data and perform unbiased benchmarks, data simulation has been widely adopted by providing explicit ground truth and generating customized datasets. However, the performance of simulation methods under multiple scenarios has not been comprehensively assessed, making it challenging to choose suitable methods without practical guidelines.

RESULTS

We systematically evaluated 49 simulation methods developed for scRNA-seq and/or SRT data in terms of accuracy, functionality, scalability, and usability using 152 reference datasets derived from 24 platforms. SRTsim, scDesign3, ZINB-WaVE, and scDesign2 have the best accuracy performance across various platforms. Unexpectedly, some methods tailored to scRNA-seq data have potential compatibility for simulating SRT data. Lun, SPARSim, and scDesign3-tree outperform other methods under corresponding simulation scenarios. Phenopath, Lun, Simple, and MFA yield high scalability scores but they cannot generate realistic simulated data. Users should consider the trade-offs between method accuracy and scalability (or functionality) when making decisions. Additionally, execution errors are mainly caused by failed parameter estimations and appearance of missing or infinite values in calculations. We provide practical guidelines for method selection, a standard pipeline Simpipe ( https://github.com/duohongrui/simpipe ; https://doi.org/10.5281/zenodo.11178409 ), and an online tool Simsite ( https://www.ciblab.net/software/simshiny/ ) for data simulation.

CONCLUSIONS

No method performs best on all criteria, thus a good-yet-not-the-best method is recommended if it solves problems effectively and reasonably. Our comprehensive work provides crucial insights for developers on modeling gene expression data and fosters the simulation process for users.

Localization of protoporphyrin IX during glioma-resection surgery via paired stimulated Raman histology and fluorescence microscopy

The most widely used fluorophore in glioma-resection surgery, 5-aminolevulinic acid (5-ALA), is thought to cause the selective accumulation of fluorescent protoporphyrin IX (PpIX) in tumour cells. Here we show that the clinical detection of PpIX can be improved via a microscope that performs paired stimulated Raman histology and two-photon excitation fluorescence microscopy (TPEF). We validated the technique in fresh tumour specimens from 115 patients with high-grade gliomas across four medical institutions. We found a weak negative correlation between tissue cellularity and the fluorescence intensity of PpIX across all imaged specimens. Semi-supervised clustering of the TPEF images revealed five distinct patterns of PpIX fluorescence, and spatial transcriptomic analyses of the imaged tissue showed that myeloid cells predominate in areas where PpIX accumulates in the intracellular space. Further analysis of external spatially resolved metabolomics, transcriptomics and RNA-sequencing datasets from glioblastoma specimens confirmed that myeloid cells preferentially accumulate and metabolize PpIX. Our findings question 5-ALA-induced fluorescence in glioma cells and show how 5-ALA and TPEF imaging can provide a window into the immune microenvironment of gliomas.

Evidences of neurological injury caused by COVID-19 from glioma tissues and glioma organoids

INTRODUCTION

Despite the extensive neurological symptoms induced by COVID-19 and the identification of SARS-CoV-2 in post-mortem brain samples from COVID-19 patients months after death, the precise mechanisms of SARS-CoV-2 invasion into the central nervous system remain unclear due to the lack of research models.

METHODS

We collected glioma tissue samples from glioma patients who had a recent history of COVID-19 and examined the presence of the SARS-CoV-2 spike protein. Subsequently, spatial transcriptomic analyses were conducted on normal brain tissues, glioma tissues, and glioma tissues from glioma patients with recent COVID-19 history. Additionally, single-cell sequencing data from both glioma tissues and glioma organoids were collected and analyzed. Glioma organoids were utilized to evaluate the efficacy of potential COVID-19 blocking agents.

RESULTS

Glioma tissues from glioma patients with recent COVID-19 history exhibited the presence of the SARS-CoV-2 spike protein. Differences between glioma tissues from glioma patients who had a recent history of COVID-19 and healthy brain tissues primarily manifested in neuronal cells. Notably, neuronal cells within glioma tissues of COVID-19 history demonstrated heightened susceptibility to Alzheimer's disease, depression, and synaptic dysfunction, indicative of neuronal aberrations. Expressions of SARS-CoV-2 entry factors were confirmed in both glioma tissues and glioma organoids. Moreover, glioma organoids were susceptible to pseudo-SARS-CoV-2 infection and the infections could be partly blocked by the potential COVID-19 drugs.

CONCLUSIONS

Gliomas had inherent traits that render them susceptible to SARS-CoV-2 infection, leading to their representability of COVID-19 neurological symptoms. This established a biological foundation for the rationality and feasibility of utilization of glioma organoids as research and blocking drug testing model in SARS-CoV-2 infection within the central nervous system.

The Role of Mesenchymal Reprogramming in Malignant Clonal Evolution and Intra-Tumoral Heterogeneity in Glioblastoma

Glioblastoma (GBM) is the most common yet uniformly fatal adult brain cancer. Intra-tumoral molecular and cellular heterogeneities are major contributory factors to therapeutic refractoriness and futility in GBM. Molecular heterogeneity is represented through molecular subtype clusters whereby the proneural (PN) subtype is associated with significantly increased long-term survival compared to the highly resistant mesenchymal (MES) subtype. Furthermore, it is universally recognized that a small subset of GBM cells known as GBM stem cells (GSCs) serve as reservoirs for tumor recurrence and progression. The clonal evolution of GSC molecular subtypes in response to therapy drives intra-tumoral heterogeneity and remains a critical determinant of GBM outcomes. In particular, the intra-tumoral MES reprogramming of GSCs using current GBM therapies has emerged as a leading hypothesis for therapeutic refractoriness. Preventing the intra-tumoral divergent evolution of GBM toward the MES subtype via new treatments would dramatically improve long-term survival for GBM patients and have a significant impact on GBM outcomes. In this review, we examine the challenges of the role of MES reprogramming in the malignant clonal evolution of glioblastoma and provide future perspectives for addressing the unmet therapeutic need to overcome resistance in GBM.

Differential gene expression analysis of spatial transcriptomic experiments using spatial mixed models

Spatial transcriptomics (ST) assays represent a revolution in how the architecture of tissues is studied by allowing for the exploration of cells in their spatial context. A common element in the analysis is delineating tissue domains or "niches" followed by detecting differentially expressed genes to infer the biological identity of the tissue domains or cell types. However, many studies approach differential expression analysis by using statistical approaches often applied in the analysis of non-spatial scRNA data (e.g., two-sample t-tests, Wilcoxon's rank sum test), hence neglecting the spatial dependency observed in ST data. In this study, we show that applying linear mixed models with spatial correlation structures using spatial random effects effectively accounts for the spatial autocorrelation and reduces inflation of type-I error rate observed in non-spatial based differential expression testing. We also show that spatial linear models with an exponential correlation structure provide a better fit to the ST data as compared to non-spatial models, particularly for spatially resolved technologies that quantify expression at finer scales (i.e., single-cell resolution).

Identification of hypoxic macrophages in glioblastoma with therapeutic potential for vasculature normalization

Monocyte-derived tumor-associated macrophages (Mo-TAMs) intensively infiltrate diffuse gliomas with remarkable heterogeneity. Using single-cell transcriptomics, we chart a spatially resolved transcriptional landscape of Mo-TAMs across 51 patients with isocitrate dehydrogenase (IDH)-wild-type glioblastomas or IDH-mutant gliomas. We characterize a Mo-TAM subset that is localized to the peri-necrotic niche and skewed by hypoxic niche cues to acquire a hypoxia response signature. Hypoxia-TAM destabilizes endothelial adherens junctions by activating adrenomedullin paracrine signaling, thereby stimulating a hyperpermeable neovasculature that hampers drug delivery in glioblastoma xenografts. Accordingly, genetic ablation or pharmacological blockade of adrenomedullin produced by Hypoxia-TAM restores vascular integrity, improves intratumoral concentration of the anti-tumor agent dabrafenib, and achieves combinatorial therapeutic benefits. Increased proportion of Hypoxia-TAM or adrenomedullin expression is predictive of tumor vessel hyperpermeability and a worse prognosis of glioblastoma. Our findings highlight Mo-TAM diversity and spatial niche-steered Mo-TAM reprogramming in diffuse gliomas and indicate potential therapeutics targeting Hypoxia-TAM to normalize tumor vasculature.

Understanding the immunosuppressive microenvironment of glioma: mechanistic insights and clinical perspectives

Glioblastoma (GBM), the predominant and primary malignant intracranial tumor, poses a formidable challenge due to its immunosuppressive microenvironment, thereby confounding conventional therapeutic interventions. Despite the established treatment regimen comprising surgical intervention, radiotherapy, temozolomide administration, and the exploration of emerging modalities such as immunotherapy and integration of medicine and engineering technology therapy, the efficacy of these approaches remains constrained, resulting in suboptimal prognostic outcomes. In recent years, intensive scrutiny of the inhibitory and immunosuppressive milieu within GBM has underscored the significance of cellular constituents of the GBM microenvironment and their interactions with malignant cells and neurons. Novel immune and targeted therapy strategies have emerged, offering promising avenues for advancing GBM treatment. One pivotal mechanism orchestrating immunosuppression in GBM involves the aggregation of myeloid-derived suppressor cells (MDSCs), glioma-associated macrophage/microglia (GAM), and regulatory T cells (Tregs). Among these, MDSCs, though constituting a minority (4-8%) of CD45+ cells in GBM, play a central component in fostering immune evasion and propelling tumor progression, angiogenesis, invasion, and metastasis. MDSCs deploy intricate immunosuppressive mechanisms that adapt to the dynamic tumor microenvironment (TME). Understanding the interplay between GBM and MDSCs provides a compelling basis for therapeutic interventions. This review seeks to elucidate the immune regulatory mechanisms inherent in the GBM microenvironment, explore existing therapeutic targets, and consolidate recent insights into MDSC induction and their contribution to GBM immunosuppression. Additionally, the review comprehensively surveys ongoing clinical trials and potential treatment strategies, envisioning a future where targeting MDSCs could reshape the immune landscape of GBM. Through the synergistic integration of immunotherapy with other therapeutic modalities, this approach can establish a multidisciplinary, multi-target paradigm, ultimately improving the prognosis and quality of life in patients with GBM.

Cellular diversity through space and time: adding new dimensions to GBM therapeutic development

The current median survival for glioblastoma (GBM) patients is only about 16 months, with many patients succumbing to the disease in just a matter of months, making it the most common and aggressive primary brain cancer in adults. This poor outcome is, in part, due to the lack of new treatment options with only one FDA-approved treatment in the last decade. Advances in sequencing techniques and transcriptomic analyses have revealed a vast degree of heterogeneity in GBM, from inter-patient diversity to intra-tumoral cellular variability. These cutting-edge approaches are providing new molecular insights highlighting a critical role for the tumor microenvironment (TME) as a driver of cellular plasticity and phenotypic heterogeneity. With this expanded molecular toolbox, the influence of TME factors, including endogenous (e.g., oxygen and nutrient availability and interactions with non-malignant cells) and iatrogenically induced (e.g., post-therapeutic intervention) stimuli, on tumor cell states can be explored to a greater depth. There exists a critical need for interrogating the temporal and spatial aspects of patient tumors at a high, cell-level resolution to identify therapeutically targetable states, interactions and mechanisms. In this review, we discuss advancements in our understanding of spatiotemporal diversity in GBM with an emphasis on the influence of hypoxia and immune cell interactions on tumor cell heterogeneity. Additionally, we describe specific high-resolution spatially resolved methodologies and their potential to expand the impact of pre-clinical GBM studies. Finally, we highlight clinical attempts at targeting hypoxia- and immune-related mechanisms of malignancy and the potential therapeutic opportunities afforded by single-cell and spatial exploration of GBM patient specimens.

Opening the Black Box: Spatial Transcriptomics and the Relevance of Artificial Intelligence-Detected Prognostic Regions in High-Grade Serous Carcinoma

Image-based deep learning models are used to extract new information from standard hematoxylin and eosin pathology slides; however, biological interpretation of the features detected by artificial intelligence (AI) remains a challenge. High-grade serous carcinoma of the ovary (HGSC) is characterized by aggressive behavior and chemotherapy resistance, but also exhibits striking variability in outcome. Our understanding of this disease is limited, partly due to considerable tumor heterogeneity. We previously trained an AI model to identify HGSC tumor regions that are highly associated with outcome status but are indistinguishable by conventional morphologic methods. Here, we applied spatially resolved transcriptomics to further profile the AI-identified tumor regions in 16 patients (8 per outcome group) and identify molecular features related to disease outcome in patients who underwent primary debulking surgery and platinum-based chemotherapy. We examined formalin-fixed paraffin-embedded tissue from (1) regions identified by the AI model as highly associated with short or extended chemotherapy response, and (2) background tumor regions (not identified by the AI model as highly associated with outcome status) from the same tumors. We show that the transcriptomic profiles of AI-identified regions are more distinct than background regions from the same tumors, are superior in predicting outcome, and differ in several pathways including those associated with chemoresistance in HGSC. Further, we find that poor outcome and good outcome regions are enriched by different tumor subpopulations, suggesting distinctive interaction patterns. In summary, our work presents proof of concept that AI-guided spatial transcriptomic analysis improves recognition of biologic features relevant to patient outcomes.

Understanding glioblastoma at the single-cell level: Recent advances and future challenges

Glioblastoma, the most aggressive and prevalent form of primary brain tumor, is characterized by rapid growth, diffuse infiltration, and resistance to therapies. Intrinsic heterogeneity and cellular plasticity contribute to its rapid progression under therapy; therefore, there is a need to fully understand these tumors at a single-cell level. Over the past decade, single-cell transcriptomics has enabled the molecular characterization of individual cells within glioblastomas, providing previously unattainable insights into the genetic and molecular features that drive tumorigenesis, disease progression, and therapy resistance. However, despite advances in single-cell technologies, challenges such as high costs, complex data analysis and interpretation, and difficulties in translating findings into clinical practice persist. As single-cell technologies are developed further, more insights into the cellular and molecular heterogeneity of glioblastomas are expected, which will help guide the development of personalized and effective therapies, thereby improving prognosis and quality of life for patients.

Spatial transcriptomics reveals segregation of tumor cell states in glioblastoma and marked immunosuppression within the perinecrotic niche

Glioblastoma (GBM) remains an untreatable malignant tumor with poor patient outcomes, characterized by palisading necrosis and microvascular proliferation. While single-cell technology made it possible to characterize different lineage of glioma cells into neural progenitor-like (NPC-like), oligodendrocyte-progenitor-like (OPC-like), astrocyte-like (AC-like) and mesenchymal like (MES-like) states, it does not capture the spatial localization of these tumor cell states. Spatial transcriptomics empowers the study of the spatial organization of different cell types and tumor cell states and allows for the selection of regions of interest to investigate region-specific and cell-type-specific pathways. Here, we obtained paired 10x Chromium single-nuclei RNA-sequencing (snRNA-seq) and 10x Visium spatial transcriptomics data from three GBM patients to interrogate the GBM microenvironment. Integration of the snRNA-seq and spatial transcriptomics data reveals patterns of segregation of tumor cell states. For instance, OPC-like tumor and NPC-like tumor significantly segregate in two of the three samples. Our differentially expressed gene and pathway analyses uncovered significant pathways in functionally relevant niches. Specifically, perinecrotic regions were more immunosuppressive than the endogenous GBM microenvironment, and perivascular regions were more pro-inflammatory. Our gradient analysis suggests that OPC-like tumor cells tend to reside in areas closer to the tumor vasculature compared to tumor necrosis, which may reflect increased oxygen requirements for OPC-like cells. In summary, we characterized the localization of cell types and tumor cell states, the gene expression patterns, and pathways in different niches within the GBM microenvironment. Our results provide further evidence of the segregation of tumor cell states and highlight the immunosuppressive nature of the necrotic and perinecrotic niches in GBM.

Single-cell and spatial transcriptomics: Bridging current technologies with long-read sequencing

Single-cell technologies have transformed biomedical research over the last decade, opening up new possibilities for understanding cellular heterogeneity, both at the genomic and transcriptomic level. In addition, more recent developments of spatial transcriptomics technologies have made it possible to profile cells in their tissue context. In parallel, there have been substantial advances in sequencing technologies, and the third generation of methods are able to produce reads that are tens of kilobases long, with error rates matching the second generation short reads. Long reads technologies make it possible to better map large genome rearrangements and quantify isoform specific abundances. This further improves our ability to characterize functionally relevant heterogeneity. Here, we show how researchers have begun to combine single-cell, spatial transcriptomics, and long-read technologies, and how this is resulting in powerful new approaches to profiling both the genome and the transcriptome. We discuss the achievements so far, and we highlight remaining challenges and opportunities.

A longitudinal single-cell and spatial multiomic atlas of pediatric high-grade glioma

Pediatric high-grade glioma (pHGG) is an incurable central nervous system malignancy that is a leading cause of pediatric cancer death. While pHGG shares many similarities to adult glioma, it is increasingly recognized as a molecularly distinct, yet highly heterogeneous disease. In this study, we longitudinally profiled a molecularly diverse cohort of 16 pHGG patients before and after standard therapy through single-nucleus RNA and ATAC sequencing, whole-genome sequencing, and CODEX spatial proteomics to capture the evolution of the tumor microenvironment during progression following treatment. We found that the canonical neoplastic cell phenotypes of adult glioblastoma are insufficient to capture the range of tumor cell states in a pediatric cohort and observed differential tumor-myeloid interactions between malignant cell states. We identified key transcriptional regulators of pHGG cell states and did not observe the marked proneural to mesenchymal shift characteristic of adult glioblastoma. We showed that essential neuromodulators and the interferon response are upregulated post-therapy along with an increase in non-neoplastic oligodendrocytes. Through in vitro pharmacological perturbation, we demonstrated novel malignant cell-intrinsic targets. This multiomic atlas of longitudinal pHGG captures the key features of therapy response that support distinction from its adult counterpart and suggests therapeutic strategies which are targeted to pediatric gliomas.

Spatial transcriptomics reveals prognosis-associated cellular heterogeneity in the papillary thyroid carcinoma microenvironment

BACKGROUND

Papillary thyroid carcinoma (PTC) is the most common malignant endocrine tumour, and its incidence and prevalence are increasing considerably. Cellular heterogeneity in the tumour microenvironment is important for PTC prognosis. Spatial transcriptomics is a powerful technique for cellular heterogeneity study.

METHODS

In conjunction with a clinical pathologist identification method, spatial transcriptomics was employed to characterise the spatial location and RNA profiles of PTC-associated cells within the tissue sections. The spatial RNA-clinical signature genes for each cell type were extracted and applied to outlining the distribution regions of specific cells on the entire section. The cellular heterogeneity of each cell type was further revealed by ContourPlot analysis, monocle analysis, trajectory analysis, ligand-receptor analysis and Gene Ontology enrichment analysis.

RESULTS

The spatial distribution region of tumour cells, typical and atypical follicular cells (FCs and AFCs) and immune cells were accurately and comprehensively identified in all five PTC tissue sections. AFCs were identified as a transitional state between FCs and tumour cells, exhibiting a higher resemblance to the latter. Three tumour foci were shared among all patients out of the 13 observed. Notably, tumour foci No. 2 displayed elevated expression levels of genes associated with lower relapse-free survival in PTC patients. We discovered key ligand-receptor interactions, including LAMB3-ITGA2, FN1-ITGA3 and FN1-SDC4, involved in the transition of PTC cells from FCs to AFCs and eventually to tumour cells. High expression of these patterns correlated with reduced relapse-free survival. In the tumour immune microenvironment, reduced interaction between myeloid-derived TGFB1 and TGFBR1 in tumour focus No. 2 contributed to tumourigenesis and increased heterogeneity. The spatial RNA-clinical analysis method developed here revealed prognosis-associated cellular heterogeneity in the PTC microenvironment.

CONCLUSIONS

The occurrence of tumour foci No. 2 and three enhanced ligand-receptor interactions in the AFC area/tumour foci reduced the relapse-free survival of PTC patients, potentially leading to improved prognostic strategies and targeted therapies for PTC patients.

BANKSY unifies cell typing and tissue domain segmentation for scalable spatial omics data analysis

Spatial omics data are clustered to define both cell types and tissue domains. We present Building Aggregates with a Neighborhood Kernel and Spatial Yardstick (BANKSY), an algorithm that unifies these two spatial clustering problems by embedding cells in a product space of their own and the local neighborhood transcriptome, representing cell state and microenvironment, respectively. BANKSY's spatial feature augmentation strategy improved performance on both tasks when tested on diverse RNA (imaging, sequencing) and protein (imaging) datasets. BANKSY revealed unexpected niche-dependent cell states in the mouse brain and outperformed competing methods on domain segmentation and cell typing benchmarks. BANKSY can also be used for quality control of spatial transcriptomics data and for spatially aware batch effect correction. Importantly, it is substantially faster and more scalable than existing methods, enabling the processing of millions of cell datasets. In summary, BANKSY provides an accurate, biologically motivated, scalable and versatile framework for analyzing spatially resolved omics data.

Single-Cell Transcriptomics Reveals the Heterogeneity of the Immune Landscape of IDH-Wild-Type High-Grade Gliomas

Isocitrate dehydrogenase (IDH)-wild-type (WT) high-grade gliomas, especially glioblastomas, are highly aggressive and have an immunosuppressive tumor microenvironment. Although tumor-infiltrating immune cells are known to play a critical role in glioma genesis, their heterogeneity and intercellular interactions remain poorly understood. In this study, we constructed a single-cell transcriptome landscape of immune cells from tumor tissue and matching peripheral blood mononuclear cells (PBMC) from IDH-WT high-grade glioma patients. Our analysis identified two subsets of tumor-associated macrophages (TAM) in tumors with the highest protumorigenesis signatures, highlighting their potential role in glioma progression. We also investigated the T-cell trajectory and identified the aryl hydrocarbon receptor (AHR) as a regulator of T-cell dysfunction, providing a potential target for glioma immunotherapy. We further demonstrated that knockout of AHR decreased chimeric antigen receptor (CAR) T-cell exhaustion and improved CAR T-cell antitumor efficacy both in vitro and in vivo. Finally, we explored intercellular communication mediated by ligand-receptor interactions within the tumor microenvironment and PBMCs and revealed the unique cellular interactions present in the tumor microenvironment. Taken together, our study provides a comprehensive immune landscape of IDH-WT high-grade gliomas and offers potential drug targets for glioma immunotherapy.

Digital profiling of cancer transcriptomes from histology images with grouped vision attention

Cancer is a heterogeneous disease that demands precise molecular profiling for better understanding and management. Recently, deep learning has demonstrated potentials for cost-efficient prediction of molecular alterations from histology images. While transformer-based deep learning architectures have enabled significant progress in non-medical domains, their application to histology images remains limited due to small dataset sizes coupled with the explosion of trainable parameters. Here, we develop SEQUOIA, a transformer model to predict cancer transcriptomes from whole-slide histology images. To enable the full potential of transformers, we first pre-train the model using data from 1,802 normal tissues. Then, we fine-tune and evaluate the model in 4,331 tumor samples across nine cancer types. The prediction performance is assessed at individual gene levels and pathway levels through Pearson correlation analysis and root mean square error. The generalization capacity is validated across two independent cohorts comprising 1,305 tumors. In predicting the expression levels of 25,749 genes, the highest performance is observed in cancers from breast, kidney and lung, where SEQUOIA accurately predicts the expression of 11,069, 10,086 and 8,759 genes, respectively. The accurately predicted genes are associated with the regulation of inflammatory response, cell cycles and metabolisms. While the model is trained at the tissue level, we showcase its potential in predicting spatial gene expression patterns using spatial transcriptomics datasets. Leveraging the prediction performance, we develop a digital gene expression signature that predicts the risk of recurrence in breast cancer. SEQUOIA deciphers clinically relevant gene expression patterns from histology images, opening avenues for improved cancer management and personalized therapies.

SCAN: Spatiotemporal Cloud Atlas for Neural cells

The nervous system is one of the most complicated and enigmatic systems within the animal kingdom. Recently, the emergence and development of spatial transcriptomics (ST) and single-cell RNA sequencing (scRNA-seq) technologies have provided an unprecedented ability to systematically decipher the cellular heterogeneity and spatial locations of the nervous system from multiple unbiased aspects. However, efficiently integrating, presenting and analyzing massive multiomic data remains a huge challenge. Here, we manually collected and comprehensively analyzed high-quality scRNA-seq and ST data from the nervous system, covering 10 679 684 cells. In addition, multi-omic datasets from more than 900 species were included for extensive data mining from an evolutionary perspective. Furthermore, over 100 neurological diseases (e.g. Alzheimer's disease, Parkinson's disease, Down syndrome) were systematically analyzed for high-throughput screening of putative biomarkers. Differential expression patterns across developmental time points, cell types and ST spots were discerned and subsequently subjected to extensive interpretation. To provide researchers with efficient data exploration, we created a new database with interactive interfaces and integrated functions called the Spatiotemporal Cloud Atlas for Neural cells (SCAN), freely accessible at http://47.98.139.124:8799 or http://scanatlas.net. SCAN will benefit the neuroscience research community to better exploit the spatiotemporal atlas of the neural system and promote the development of diagnostic strategies for various neurological disorders.

NRP1 Induces Enhanced Stemness and Chemoresistance in Glioma Cells via YAP

Neuropilin-1 (NRP1), a transmembrane glycoprotein, plays an important role in the malignant progression of gliomas; however, its role in chemoresistance is not fully understood. In this study, we observed the effects of NRP1 on the stemness and chemoresistance of glioma cells and the mediating role of Yes-associated protein (YAP). We constructed NRP1 overexpressing LN-229 glioma cells. Cells were treated with recombinant NRP1 protein (rNRP1) and the YAP inhibitor Super-TDU when necessary. 3-(4,5-Dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) assay was used to detect the sensitivity of cells to temozolomide (TMZ). Sphere and clone formation assays were performed to detect the sphere- and clone-forming abilities of cells. Western blotting was performed to detect cellular CD133, CD44, p-LATS1, and p-YAP protein expression. Immunofluorescence and flow cytometry were used to detect the subcellular localization of YAP and apoptosis, respectively. We found that both NRP1 overexpression and rNRP1 treatment enhanced self-renewal, TMZ resistance, and CD133 and CD44 protein expression in LN-229 cells. NRP1 overexpression and rNRP1 treatment also induced LATS1 and YAP dephosphorylation and YAP nuclear translocation. Super-TDU inhibits NRP1 overexpression-induced enhanced self-renewal and TMZ resistance in LN-229 cells. Our study suggests that NRP1 induces increased stemness in glioma cells, resulting in chemoresistance, and that this effect is associated with YAP activation.

Cellular Components of the Tumor Environment in Gliomas-What Do We Know Today?

A generation ago, the molecular properties of tumor cells were the focus of scientific interest in oncology research. Since then, it has become increasingly apparent that the tumor environment (TEM), whose major components are non-neoplastic cell types, is also of utmost importance for our understanding of tumor growth, maintenance and resistance. In this review, we present the current knowledge concerning all cellular components within the TEM in gliomas, focusing on their molecular properties, expression patterns and influence on the biological behavior of gliomas. Insight into the TEM of gliomas has expanded considerably in recent years, including many aspects that previously received only marginal attention, such as the phenomenon of phagocytosis of glioma cells by macrophages and the role of the thyroid-stimulating hormone on glioma growth. We also discuss other topics such as the migration of lymphocytes into the tumor, phenotypic similarities between chemoresistant glioma cells and stem cells, and new clinical approaches with immunotherapies involving the cells of TEM.

AEP-cleaved DDX3X induces alternative RNA splicing events to mediate cancer cell adaptation in harsh microenvironments

Oxygen and nutrient deprivation are common features of solid tumors. Although abnormal alternative splicing (AS) has been found to be an important driving force in tumor pathogenesis and progression, the regulatory mechanisms of AS that underly the adaptation of cancer cells to harsh microenvironments remain unclear. Here, we found that hypoxia- and nutrient deprivation-induced asparagine endopeptidase (AEP) specifically cleaved DDX3X in a HIF1A-dependent manner. This cleavage yields truncated carboxyl-terminal DDX3X (tDDX3X-C), which translocates and aggregates in the nucleus. Unlike intact DDX3X, nuclear tDDX3X-C complexes with an array of splicing factors and induces AS events of many pre-mRNAs; for example, enhanced exon skipping (ES) in exon 2 of the classic tumor suppressor PRDM2 leads to a frameshift mutation of PRDM2. Intriguingly, the isoform ARRB1-Δexon 13 binds to glycolytic enzymes and regulates glycolysis. By utilizing in vitro assays, glioblastoma organoids, and animal models, we revealed that AEP/tDDX3X-C promoted tumor malignancy via these isoforms. More importantly, high AEP/tDDX3X-C/ARRB1-Δexon 13 in cancerous tissues was tightly associated with poor patient prognosis. Overall, our discovery of the effect of AEP-cleaved DDX3X switching on alternative RNA splicing events identifies a mechanism in which cancer cells adapt to oxygen and nutrient shortages and provides potential diagnostic and/or therapeutic targets.

Spatial transcriptomics in glioblastoma: is knowing the right zip code the key to the next therapeutic breakthrough?

Spatial transcriptomics, the technology of visualizing cellular gene expression landscape in a cells native tissue location, has emerged as a powerful tool that allows us to address scientific questions that were elusive just a few years ago. This technological advance is a decisive jump in the technological evolution that is revolutionizing studies of tissue structure and function in health and disease through the introduction of an entirely new dimension of data, spatial context. Perhaps the organ within the body that relies most on spatial organization is the brain. The central nervous system's complex microenvironmental and spatial architecture is tightly regulated during development, is maintained in health, and is detrimental when disturbed by pathologies. This inherent spatial complexity of the central nervous system makes it an exciting organ to study using spatial transcriptomics for pathologies primarily affecting the brain, of which Glioblastoma is one of the worst. Glioblastoma is a hyper-aggressive, incurable, neoplasm and has been hypothesized to not only integrate into the spatial architecture of the surrounding brain, but also possess an architecture of its own that might be actively remodeling the surrounding brain. In this review we will examine the current landscape of spatial transcriptomics in glioblastoma, outline novel findings emerging from the rising use of spatial transcriptomics, and discuss future directions and ultimate clinical/translational avenues.

Spatial transcriptomics in development and disease

The proper functioning of diverse biological systems depends on the spatial organization of their cells, a critical factor for biological processes like shaping intricate tissue functions and precisely determining cell fate. Nonetheless, conventional bulk or single-cell RNA sequencing methods were incapable of simultaneously capturing both gene expression profiles and the spatial locations of cells. Hence, a multitude of spatially resolved technologies have emerged, offering a novel dimension for investigating regional gene expression, spatial domains, and interactions between cells. Spatial transcriptomics (ST) is a method that maps gene expression in tissue while preserving spatial information. It can reveal cellular heterogeneity, spatial organization and functional interactions in complex biological systems. ST can also complement and integrate with other omics methods to provide a more comprehensive and holistic view of biological systems at multiple levels of resolution. Since the advent of ST, new methods offering higher throughput and resolution have become available, holding significant potential to expedite fresh insights into comprehending biological complexity. Consequently, a rapid increase in associated research has occurred, using these technologies to unravel the spatial complexity during developmental processes or disease conditions. In this review, we summarize the recent advancement of ST in historical, technical, and application contexts. We compare different types of ST methods based on their principles and workflows, and present the bioinformatics tools for analyzing and integrating ST data with other modalities. We also highlight the applications of ST in various domains of biomedical research, especially development and diseases. Finally, we discuss the current limitations and challenges in the field, and propose the future directions of ST.

Combining single-cell sequencing and spatial transcriptome sequencing to identify exosome-related features of glioblastoma and constructing a prognostic model to identify BARD1 as a potential therapeutic target for GBM patients

Background

Glioblastoma (GBM) is a malignant primary brain tumor. This study focused on exploring the exosome-related features of glioblastoma to better understand its cellular composition and molecular characteristics.

Methods

Single-cell RNA sequencing (scRNA-seq) and spatial transcriptome RNA sequencing (stRNA-seq) were used to analyze the heterogeneity of glioblastomas. After data integration, cell clustering, and annotation, five algorithms were used to calculate scores for exosome-related genes(ERGs). Cell trajectory analysis and intercellular communication analysis were performed to explore exosome-related communication patterns. Spatial transcriptome sequencing data were analyzed to validate the findings. To further utilize exosome-related features to aid in clinical decision-making, a prognostic model was constructed using GBM's bulk RNA-seq.

Results

Different cell subpopulations were observed in GBM, with Monocytes/macrophages and malignant cells in tumor samples showing higher exosome-related scores. After identifying differentially expressed ERGs in malignant cells, pseudotime analysis revealed the cellular status of malignant cells during development. Intercellular communication analysis highlighted signaling pathways and ligand-receptor interactions. Spatial transcriptome sequencing confirmed the high expression of exosome-related gene features in the tumor core region. A prognostic model based on six ERGs was shown to be predictive of overall survival and immunotherapy outcome in GBM patients. Finally, based on the results of scRNA-seq and prognostic modeling as well as a series of cell function experiments, BARD1 was identified as a novel target for the treatment of GBM.

Conclusion

This study provides a comprehensive understanding of the exosome-related features of GBM in both scRNA-seq and stRNA-seq, with malignant cells with higher exosome-related scores exhibiting stronger communication with Monocytes/macrophages. In terms of spatial data, highly scored malignant cells were also concentrated in the tumor core region. In bulk RNA-seq, patients with a high exosome-related index exhibited an immunosuppressive microenvironment, which was accompanied by a worse prognosis as well as immunotherapy outcomes. Prognostic models constructed using ERGs are expected to be independent prognostic indicators for GBM patients, with potential implications for personalized treatment strategies for GBM. Knockdown of BARD1 in GBM cell lines reduces the invasive and value-added capacity of tumor cells, and thus BARD1-positively expressing malignant cells are a risk factor for GBM patients.

Spatial cellular architecture predicts prognosis in glioblastoma

Intra-tumoral heterogeneity and cell-state plasticity are key drivers for the therapeutic resistance of glioblastoma. Here, we investigate the association between spatial cellular organization and glioblastoma prognosis. Leveraging single-cell RNA-seq and spatial transcriptomics data, we develop a deep learning model to predict transcriptional subtypes of glioblastoma cells from histology images. Employing this model, we phenotypically analyze 40 million tissue spots from 410 patients and identify consistent associations between tumor architecture and prognosis across two independent cohorts. Patients with poor prognosis exhibit higher proportions of tumor cells expressing a hypoxia-induced transcriptional program. Furthermore, a clustering pattern of astrocyte-like tumor cells is associated with worse prognosis, while dispersion and connection of the astrocytes with other transcriptional subtypes correlate with decreased risk. To validate these results, we develop a separate deep learning model that utilizes histology images to predict prognosis. Applying this model to spatial transcriptomics data reveal survival-associated regional gene expression programs. Overall, our study presents a scalable approach to unravel the transcriptional heterogeneity of glioblastoma and establishes a critical connection between spatial cellular architecture and clinical outcomes.

Dissecting the brain with spatially resolved multi-omics

Recent studies have highlighted spatially resolved multi-omics technologies, including spatial genomics, transcriptomics, proteomics, and metabolomics, as powerful tools to decipher the spatial heterogeneity of the brain. Here, we focus on two major approaches in spatial transcriptomics (next-generation sequencing-based technologies and image-based technologies), and mass spectrometry imaging technologies used in spatial proteomics and spatial metabolomics. Furthermore, we discuss their applications in neuroscience, including building the brain atlas, uncovering gene expression patterns of neurons for special behaviors, deciphering the molecular basis of neuronal communication, and providing a more comprehensive explanation of the molecular mechanisms underlying central nervous system disorders. However, further efforts are still needed toward the integrative application of multi-omics technologies, including the real-time spatial multi-omics analysis in living cells, the detailed gene profile in a whole-brain view, and the combination of functional verification.

Combined bulk RNA and single-cell RNA analyses reveal TXNL4A as a new biomarker for hepatocellular carcinoma

Introduction

Hepatocellular carcinoma (HCC) has a high mortality rate worldwide. The dysregulation of RNA splicing is a major event leading to the occurrence, progression, and drug resistance of cancer. Therefore, it is important to identify new biomarkers of HCC from the RNA splicing pathway.

Methods

We performed the differential expression and prognostic analyses of RNA splicing-related genes (RRGs) using The Cancer Genome Atlas-liver hepatocellular carcinoma (LIHC). The International Cancer Genome Consortium (ICGC)-LIHC dataset was used to construct and validate prognostic models, and the PubMed database was used to explore genes in the models to identify new markers. The screened genes were subjected to genomic analyses, including differential, prognostic, enrichment, and immunocorrelation analyses. Single-cell RNA (scRNA) data were used to further validate the immunogenetic relationship.

Results

Of 215 RRGs, we identified 75 differentially expressed prognosis-related genes, and a prognostic model incorporating thioredoxin like 4A (TXNL4A) was identified using least absolute shrinkage and selection operator regression analysis. ICGC-LIHC was used as a validation dataset to confirm the validity of the model. PubMed failed to retrieve HCC-related studies on TXNL4A. TXNL4A was highly expressed in most tumors and was associated with HCC survival. Chi-squared analyses indicated that TXNL4A expression positively correlated positively with the clinical features of HCC. Multivariate analyses revealed that high TXNL4A expression was an independent risk factor for HCC. Immunocorrelation and scRNA data analyses indicated that TXNL4A was correlated with CD8 T cell infiltration in HCC.

Conclusion

Therefore, we identified a prognostic and immune-related marker for HCC from the RNA splicing pathway.

Dissecting the tumor microenvironment of epigenetically driven gliomas: Opportunities for single-cell and spatial multiomics

Malignant gliomas are incurable brain neoplasms with dismal prognoses and near-universal fatality, with minimal therapeutic progress despite billions of dollars invested in research and clinical trials over the last 2 decades. Many glioma studies have utilized disparate histologic and genomic platforms to characterize the stunning genomic, transcriptomic, and immunologic heterogeneity found in gliomas. Single-cell and spatial omics technologies enable unprecedented characterization of heterogeneity in solid malignancies and provide a granular annotation of transcriptional, epigenetic, and microenvironmental states with limited resected tissue. Heterogeneity in gliomas may be defined, at the broadest levels, by tumors ostensibly driven by epigenetic alterations (IDH- and histone-mutant) versus non-epigenetic tumors (IDH-wild type). Epigenetically driven tumors are defined by remarkable transcriptional programs, immunologically distinct microenvironments, and incompletely understood topography (unique cellular neighborhoods and cell-cell interactions). Thus, these tumors are the ideal substrate for single-cell multiomic technologies to disentangle the complex intra-tumoral features, including differentiation trajectories, tumor-immune cell interactions, and chromatin dysregulation. The current review summarizes the applications of single-cell multiomics to existing datasets of epigenetically driven glioma. More importantly, we discuss future capabilities and applications of novel multiomic strategies to answer outstanding questions, enable the development of potent therapeutic strategies, and improve personalized diagnostics and treatment via digital pathology.