Seurat24式太極拳之4左右摟膝拗步---整合多個空間轉錄組切片

四、左右摟膝拗步
1、左摟膝拗步
轉體擺臂 擺臂收腳 上步屈肘 弓步摟推
2、右摟膝拗步
后坐撇腳 擺臂收腳 上步屈肘 弓步摟推
3、左摟膝拗步
后坐撇腳 擺臂收腳 上步屈肘 弓步摟推

在Seurat中處理多個切片

小鼠大腦的該數據集包含與大腦另一半相對應的另一個切片。在這里,我們將其讀入并執(zhí)行相同的初始歸一化。

brain2 <- LoadData("stxBrain", type = "posterior1")
brain2 <- SCTransform(brain2, assay = "Spatial", verbose = FALSE)

為了在同一個Seurat對象中使用多個切片,我們提供了該merge功能。

brain.merge <- merge(brain, brain2)

然后,這使得聯(lián)合維數減少并在基礎RNA表達數據上聚類。

DefaultAssay(brain.merge) <- "SCT"
VariableFeatures(brain.merge) <- c(VariableFeatures(brain), VariableFeatures(brain2))
brain.merge <- RunPCA(brain.merge, verbose = FALSE)
brain.merge <- FindNeighbors(brain.merge, dims = 1:30)
brain.merge <- FindClusters(brain.merge, verbose = FALSE)
brain.merge <- RunUMAP(brain.merge, dims = 1:30)

最后,可以在單個UMAP圖中共同可視化數據。SpatialDimPlot()SpatialFeaturePlot()`默認將所有切片繪制為列,將分組/功能繪制為行。

DimPlot(brain.merge, reduction = "umap", group.by = c("ident", "orig.ident"))
image
SpatialDimPlot(brain.merge)
image
SpatialFeaturePlot(brain.merge, features = c("Hpca", "Plp1"))
image.png

致謝

我們要感謝Nigel Delaney和Stephen Williams對新的空間Seurat代碼的有益反饋和貢獻

?著作權歸作者所有,轉載或內容合作請聯(lián)系作者
【社區(qū)內容提示】社區(qū)部分內容疑似由AI輔助生成,瀏覽時請結合常識與多方信息審慎甄別。
平臺聲明:文章內容(如有圖片或視頻亦包括在內)由作者上傳并發(fā)布,文章內容僅代表作者本人觀點,簡書系信息發(fā)布平臺,僅提供信息存儲服務。

相關閱讀更多精彩內容

友情鏈接更多精彩內容