SpotClean is a computational method to adjust for spot swapping in spatial transcriptomics data. Recent spatial transcriptomics experiments utilize slides containing thousands of spots with spot-specific barcodes that bind mRNA. Ideally, unique molecular identifiers at a spot measure spot-specific expression, but this is often not the case due to bleed from nearby spots, an artifact we refer to as spot swapping. SpotClean is able to estimate the contamination rate in observed data and decontamin
中文简介正在准备中。当前显示原文说明。
实测评分 未实测
查看安装方法安装前可以查看交给 AI 代理的指令内容。
Load raw 10x Visium Space Ranger output, perform decontamination, visualize results, and convert to a Seurat object.
library(SpotClean)
library(S4Vectors)
# 1. Load raw data and slide information
data(mbrain_raw)
spatial_dir <- system.file(file.path("extdata", "V1_Adult_Mouse_Brain_spatial"), package = "SpotClean")
mbrain_slide_info <- read10xSlide(
tissue_csv_file = file.path(spatial_dir, "tissue_positions_list.csv"),
tissue_img_file = file.path(spatial_dir, "tissue_lowres_image.png"),
scale_factor_file = file.path(spatial_dir, "scalefactors_json.json")
)
# 2. Create the slide object
slide_obj <- createSlide(count_mat = mbrain_raw, slide_info = mbrain_slide_info)
# 3. Visualize raw data
visualizeSlide(slide_obj = slide_obj)
visualizeHeatmap(slide_obj, "Mbp")
# 4. Decontaminate raw data
decont_obj <- spotclean(slide_obj, maxit = 10, candidate_radius = 20)
# 5. Visualize decontaminated gene and contamination rate
visualizeHeatmap(decont_obj, "Mbp")
visualizeHeatmap(
decont_obj,
metadata(decont_obj)$contamination_rate,
logged = FALSE,
legend_title = "contamination rate",
legend_range = c(0, 1)
)
# 6. Convert to Seurat object
seurat_obj <- convertToSeurat(decont_obj, image_dir = spatial_dir)
Input: Raw count matrix and 10x Visium slide metadata. Output: A decontaminated Seurat object ready for downstream analysis.
Decontaminate spatial transcriptomics data stored in a SpatialExperiment object.
library(SpotClean)
library(SpatialExperiment)
# Assuming slide_obj is a SpatialExperiment loaded with read10xVisium(..., data = "raw")
# slide_obj <- read10xVisium(samples = "/path/to/spaceranger/output/", data = "raw")
# decont_obj <- spotclean(slide_obj)
Input: A SpatialExperiment object containing raw spatial transcriptomics data. Output: A decontaminated SpatialExperiment object.
spotclean().arcScore().visualizeLabel() and visualizeHeatmap().convertToSeurat().spotclean() relies on background spots to estimate contamination.spotclean() reassigns UMIs and changes estimated sequencing depth, which may conflict with standard depth normalization.read10xRaw()).read10xSlide()).SpatialExperiment object constructed with data = "raw".createSlide().createSlide().30): Maximum number of iterations for the decontamination optimization in spotclean().spotclean().convertToSeurat().TRUE): Logical flag in visualizeHeatmap() to apply log-scaling to the heatmap values.spotclean(), specifying data = "raw" if loading via read10xVisium().visualizeHeatmap() before and after running spotclean() to verify the reduction of background noise.spotclean() on slides with nearly 100% tissue coverage leads to unreliable contamination estimates; ensure at least 25% of spots are background.read10xRaw() or read10xVisium(..., data = "raw").Seurat: For general spatial transcriptomics downstream analysis and clustering.scater: For single-cell and spatial QC and visualization.scran: For variance modeling and normalization of single-cell/spatial data.还没有评价。欢迎分享你的使用感受。
无需登录即可发表。欢迎分享使用感受、提问，或说明在哪些环境下无法运行。
还没有评论。来写第一条吧。
The package contains functions for exploratory oligonucleotide array analysis. The dependence on tkWidgets only concerns few convenience functions. 'affy' is fully functional without it.
中文简介正在准备中。当前显示原文说明。
This tool uses the label-tree function from HyPhy to annotate a phylogenetic tree. It allows users to select a subset of leaves using either a regular expression or a list of sequence names, and then apply a specified label to these selecte
中文简介正在准备中。当前显示原文说明。
This package provides a client for the Bioconductor AnnotationHub web resource. The AnnotationHub web resource provides a central location where genomic files (e.g., VCF, bed, wig) and other resources from standard locations (e.g., UCSC, En
中文简介正在准备中。当前显示原文说明。
apeglm provides Bayesian shrinkage estimators for effect sizes for a variety of GLM models, using approximation of the posterior for individual coefficients.
中文简介正在准备中。当前显示原文说明。
ASURAT is a software for single-cell data analysis. Using ASURAT, one can simultaneously perform unsupervised clustering and biological interpretation in terms of cell type, disease, biological process, and signaling pathway activity. Inputting a single-cell RNA-seq data and knowledge-based databases, such as Cell Ontology, Gene Ontology, KEGG, etc., ASURAT transforms gene expression tables into original multivariate tables, termed sign-by-sample matrices (SSMs).
中文简介正在准备中。当前显示原文说明。
AUCell allows to identify cells with active gene sets (e.g. signatures, gene modules...) in single-cell RNA-seq data. AUCell uses the "Area Under the Curve" (AUC) to calculate whether a critical subset of the input gene set is enriched with
中文简介正在准备中。当前显示原文说明。