The raw ST reads were aligned to the same reference package as for the snRNA-seq dataset using SpaceRanger (version 1.2.2, 10x Genomics) software to generate spot barcode-to-gene feature matrix for downstream analysis. Individual ST object was created by Load10X_Spatial() function with Seurat (v3.2.2) package (146 (link)) and normalized with SCTransform() function (147 (link)). The merge() function was used to aggregate 16 slices into one ST object (Visi) containing images, and the Visi@assays$SCT@counts was pulled to create a separate ST object (VisiDot) by CreateSeuratObject() function to get aggregated spots without images. A processing pipeline similar to that of snRNA-seq was employed for the VisiDot object; specifically, Harmony() was applied over “IL01_uniqueID” variable to integrate 16 samples. To match the terminal MRI images, the orientation of the Visium images was corrected by transforming the image array stored at Visi@images$slice1@image, and the spot information stored at Visi@images$slice1$@coordinates was swapped accordingly (See GitHub post for detail). The microenvironment (ME) cluster annotated in the VisiDot object was transferred back to the Visi object for spatial visualization.
Free full text: Click here