Recent advances in spatial transcriptomics technologies have enabled comprehensive characterization of gene expression patterns in the context of tissue microenvironment. To elucidate spatial gene expression variation, we present SpaGCN, a graph convolutional network approach that integrates gene expression, spatial location and histology in spatial transcriptomics data analysis. Through graph convolution, SpaGCN aggregates gene expression of each spot from its neighboring spots, which enables the identification of spatial domains with coherent expression and histology. The subsequent domain guided differential expression analysis then detects genes with enriched expression patterns in the identified domains.Analyzing five spatially resolved transcriptomics datasets using SpaGCN, we show it can detect genes with much more enriched spatial expression patterns than existing methods. Furthermore, genes detected by SpaGCN are transferrable and can be utilized to study spatial variation of gene expression in other datasets. SpaGCN is computationally fast, making it a desirable tool for spatial transcriptomics studies.

Recent advances in spatial transcriptomics technologies have enabled comprehensive characterization of gene expression patterns in the context of tissue microenvironment. To elucidate spatial gene expression variation, we present SpaGCN, a graph convolutional network approach that integrates gene expression, spatial location and histology in spatial transcriptomics data analysis. Through graph convolution, SpaGCN aggregates gene expression of each spot from its neighboring spots, which enables the identification of spatial domains with coherent expression and histology. The subsequent domain guided differential expression analysis then detects genes with enriched expression patterns in the identified domains. Analyzing five spatially resolved transcriptomics datasets using SpaGCN, we show it can detect genes with much more enriched spatial expression patterns than existing methods. Furthermore, genes detected by SpaGCN are transferrable and can be utilized to study spatial variation of gene expression in other datasets. SpaGCN is computationally fast, making it a desirable tool for spatial transcriptomics studies.

Graphical abstract

Kidneys have one of the most complex three-dimensional cellular organizations in the body, but the spatial molecular principles of kidney health and disease are poorly understood. Here we generate high-quality single cell (sc), single nuclear (sn), spatial (sp) RNA expression and sn open chromatin datasets for 73 samples, capturing half a million cells from healthy, diabetic, and hypertensive diseased human kidneys. Combining the sn/sc and sp RNA information, we identify > 100 cell types and states and successfully map them back to their spatial locations. Computational deconvolution of spRNA-seq identifies glomerular/vascular, tubular, immune, and fibrotic spatial microenvironments (FMEs). Although injured proximal tubule cells appear to be the nidus of fibrosis, we reveal the complex, heterogenous cellular and spatial organization of human FMEs, including the highly intricate and organized immune environment. We demonstrate the clinical utility of the FME spatial gene signature for the classification of a large number of human kidneys for disease severity and prognosis. We provide a comprehensive spatially-resolved molecular roadmap for the human kidney and the fibrotic process and demonstrate the clinical utility of spatial transcriptomics.

No abstract