Harnessing computational spatial omics to explore the spatial biology intricacies

Spatial analyses revealed CXCL5 and SLC6A14 as the markers of microvascular invasion in intrahepatic cholangiocarcinoma

BACKGROUND

Microvascular invasion (MVI) is a critical prognostic factor in intrahepatic cholangiocarcinoma (ICC), strongly associated with postoperative recurrence. However, the phenotypic features and spatial organization of MVI remain inadequately understood.

METHODS

We performed a spatial transcriptomic analysis on 29,632 spots from six ICC samples, manually delineating MVI clusters using the cloupe software. Key biomarkers were identified and validated in an independent cohort of 135 ICC patients. Functional and survival analyses were conducted to assess clinical relevance, and cell-cell communication pathways were investigated.

RESULTS

MVI regions exhibited heightened proliferation, angiogenesis, and epithelial-mesenchymal transition, driven by increased expression of transcription factors SOX10, ZEB1, and SNAI2. CXCL5 and SLC6A14 were identified as potential MVI biomarkers and showed high expression in tumor-invasive areas. Serum CXCL5 demonstrated strong predictive power for vascular invasion (AUC = 0.92) and intrahepatic metastasis (AUC = 0.96). High expression of both CXCL5 and SLC6A14 was associated with the worst survival outcomes. MVI regions were enriched with immunosuppressive MRC1+ macrophages and exhibited elevated immune checkpoint expression, including HAVCR2 and TIGHT, indicative of immune resistance. Cell-cell communication analysis revealed CXCL5-CXCR2 and LGALS9-HAVCR2 as key ligand-receptor pairs contributing to the immunosuppressive microenvironment.

CONCLUSIONS

This study identifies CXCL5 and SLC6A14 as key biomarkers of MVI, highlighting their roles in tumor proliferation, immune resistance, and poor clinical outcomes. These findings provide valuable insights into the spatial organization of MVI and its contribution to ICC progression, offering potential therapeutic targets.

Deciphering spatial domains from spatially resolved transcriptomics through spatially regularized deep graph networks

BACKGROUND

Recent advancements in spatially resolved transcriptomics (SRT) have opened up unprecedented opportunities to explore gene expression patterns within spatial contexts. Deciphering spatial domains is a critical task in spatial transcriptomic data analysis, aiding in the elucidation of tissue structural heterogeneity and biological functions. However, existing spatial domain detection methods ignore the consistency of expression patterns and spatial arrangements between spots, as well as the severe gene dropout phenomenon present in SRT data, resulting in suboptimal performance in identifying tissue spatial heterogeneity.

RESULTS

In this paper, we introduce a novel framework, spatially regularized deep graph networks (SR-DGN), which integrates gene expression profiles with spatial information to learn spatially-consistent and informative spot representations. Specifically, SR-DGN employs graph attention networks (GAT) to adaptively aggregate gene expression information from neighboring spots, considering local expression patterns between spots. In addition, the spatial regularization constraint ensures the consistency of neighborhood relationships between physical and embedded spaces in an end-to-end manner. SR-DGN also employs cross-entropy (CE) loss to model gene expression states, effectively mitigating the impact of noisy gene dropouts.

CONCLUSIONS

Experimental results demonstrate that SR-DGN outperforms state-of-the-art methods in spatial domain identification across SRT data from different sequencing platforms. Moreover, SR-DGN is capable of recovering known microanatomical structures, yielding clearer low-dimensional visualizations and more accurate spatial trajectory inferences.

Integrated analysis of single-cell, spatial and bulk RNA-sequencing identifies a cell-death signature for predicting the outcomes of head and neck cancer

Background

Cell death plays an essential role in carcinogenesis, but its function in the recurrence and postoperative prognosis of head and neck cancer (HNC), which ranks as the 7th most common malignancy globally, remains unclear.

Methods

Data from five main subtypes of HNC related single-cell RNA sequencing (scRNA-seq) were recruited to establish a single-cell atlas, and the distribution of cell death models (CDMs) across different tissues as well as cell subtypes were analyzed. Bulk RNA-seq from the Cancer Genome Atlas Program (TCGA) dataset was subjected to a machine learning-based integrative procedure for constructing a consensus cell death-related signature risk score (CDRscore) model and validated by external data. The biofunctions including different expression analysis, immune cell infiltration, genomic mutations, enrichment analysis as well as cellchat analysis were compared between the high- and low- risk score groups categorized by this CDRscore model. Finally, samples from laryngeal squamous cell cancer (LSCC) were conducted by spatial transcriptomics (ST) to further validate the results of CDRscore model.

Results

T cells from HNC patients manifested the highest levels of cell death while HPV infection attenuates malignant cell death based on single-cell atlas. CDMs are positively correlated with the tumor-cell stemness, immune-related score and T cells are infiltrated. A CDRscore model was established based on the transcription of ten cell death prognostic genes (MRPL10, DDX19A, NDFIP1, PCMT1, HPRT1, SLC2A3, EFNB2, HK1, BTG3 and MAP2K7). It functions as an independent prognostic factor for overall survival in HNC and displays stable and powerful performance validated by GSE41613 and GSE65858 datasets. Patients in high CDRscore manifested worse overall survival, more active of epithelial mesenchymal transition, TGF-β-related pathways and hypoxia, higher transcription of T cell exhausted markers, and stronger TP53 mutation. ST from LSCC showed that spots with high-risk scores were colocalized with TGF-β and the proliferating malignant cells, additionally, the risk scores have a negative correlation with TCR signaling but positive association with LAG3 transcription.

Conclusion

The CDRscore model could be utilized as a powerful prognostic indicator for HNC.

Pianno: a probabilistic framework automating semantic annotation for spatial transcriptomics

Spatial transcriptomics has revolutionized the study of gene expression within tissues, while preserving spatial context. However, annotating spatial spots' biological identity remains a challenge. To tackle this, we introduce Pianno, a Bayesian framework automating structural semantics annotation based on marker genes. Comprehensive evaluations underscore Pianno's remarkable prowess in precisely annotating a wide array of spatial semantics, ranging from diverse anatomical structures to intricate tumor microenvironments, as well as in estimating cell type distributions, across data generated from various spatial transcriptomics platforms. Furthermore, Pianno, in conjunction with clustering approaches, uncovers a region- and species-specific excitatory neuron subtype in the deep layer 3 of the human neocortex, shedding light on cellular evolution in the human neocortex. Overall, Pianno equips researchers with a robust and efficient tool for annotating diverse biological structures, offering new perspectives on spatial transcriptomics data.

Applications of single‑cell omics and spatial transcriptomics technologies in gastric cancer (Review)

Gastric cancer (GC) is a prominent contributor to global cancer-related mortalities, and a deeper understanding of its molecular characteristics and tumor heterogeneity is required. Single-cell omics and spatial transcriptomics (ST) technologies have revolutionized cancer research by enabling the exploration of cellular heterogeneity and molecular landscapes at the single-cell level. In the present review, an overview of the advancements in single-cell omics and ST technologies and their applications in GC research is provided. Firstly, multiple single-cell omics and ST methods are discussed, highlighting their ability to offer unique insights into gene expression, genetic alterations, epigenomic modifications, protein expression patterns and cellular location in tissues. Furthermore, a summary is provided of key findings from previous research on single-cell omics and ST methods used in GC, which have provided valuable insights into genetic alterations, tumor diagnosis and prognosis, tumor microenvironment analysis, and treatment response. In summary, the application of single-cell omics and ST technologies has revealed the levels of cellular heterogeneity and the molecular characteristics of GC, and holds promise for improving diagnostics, personalized treatments and patient outcomes in GC.

Digital Spatial Profiling Identifies the Tumor Periphery as a Highly Active Biological Niche in Clear Cell Renal Cell Carcinoma

Clear cell renal cell carcinoma (ccRCC) is characterized by a high degree of intratumoral heterogeneity (ITH). Besides genomic ITH, there is considerable functional ITH, which encompasses spatial niches with distinct proliferative and signaling activities. The full extent of functional spatial heterogeneity in ccRCC is incompletely understood. In the present study, a total of 17 ccRCC tissue specimens from different sites (primary tumor, n = 11; local recurrence, n = 1; distant metastasis, n = 5) were analyzed using digital spatial profiling (DSP) of protein expression. A total of 128 regions of interest from the tumor periphery and tumor center were analyzed for the expression of 46 proteins, comprising three major signaling pathways as well as immune cell markers. Results were correlated to clinico-pathological variables. The differential expression of granzyme B was validated using conventional immunohistochemistry and was correlated to the cancer-specific patient survival. We found that a total of 37 proteins were differentially expressed between the tumor periphery and tumor center. Thirty-five of the proteins were upregulated in the tumor periphery compared to the center. These included proteins involved in cell proliferation, MAPK and PI3K/AKT signaling, apoptosis regulation, epithelial-to-mesenchymal transition, as well as immune cell markers. Among the most significantly upregulated proteins in the tumor periphery was granzyme B. Granzyme B upregulation in the tumor periphery correlated with a significantly reduced cancer-specific patient survival. In conclusion, this study highlights the unique cellular contexture of the tumor periphery in ccRCC. The correlation between granzyme B upregulation in the tumor periphery and patient survival suggests local selection pressure for aggressive tumor growth and disease progression. Our results underscore the potential of spatial biology for biomarker discovery in ccRCC and cancer in general.