SpaTalk Advanced Usage

# Custom LR pairs must have these columns:
# - ligand: ligand gene symbol
# - receptor: receptor gene symbol  
# - species: "Human" or "Mouse"

custom_lrpairs <- data.frame(
  ligand = c("CXCL12", "CCL2", "VEGFA"),
  receptor = c("CXCR4", "CCR2", "KDR"),
  species = "Human"
)

# Use in find_lr_path
obj <- find_lr_path(obj, lrpairs = custom_lrpairs, pathways = pathways)

# Custom pathways must have:
# - src: source gene
# - dest: destination gene (downstream)
# - type: interaction type (optional)

custom_pathways <- data.frame(
  src = c("CXCR4", "CXCR4", "STAT3"),
  dest = c("STAT3", "AKT1", "MYC"),
  type = c("activation", "activation", "transcription")
)

obj <- find_lr_path(obj, lrpairs = lrpairs, pathways = custom_pathways)

obj <- dec_celltype(
  object = obj,
  sc_data = sc_data,
  sc_celltype = sc_celltype,
  method = 1  # NNLM
)

# Requires: install.packages("spacexr")
obj <- dec_celltype(
  object = obj,
  sc_data = sc_data,
  sc_celltype = sc_celltype,
  method = 2  # RCTD
)

obj <- dec_celltype(
  object = obj,
  sc_data = sc_data,
  sc_celltype = sc_celltype,
  method = 3  # Seurat
)

# Requires: BiocManager::install("SPOTlight")
obj <- dec_celltype(
  object = obj,
  sc_data = sc_data,
  sc_celltype = sc_celltype,
  method = 4  # SPOTlight
)

# Requires: devtools::install_github("reneshbedre/deconvSeq")
obj <- dec_celltype(
  object = obj,
  sc_data = sc_data,
  sc_celltype = sc_celltype,
  method = 5  # deconvSeq
)

# Requires: conda environment with stereoscope
obj <- dec_celltype(
  object = obj,
  sc_data = sc_data,
  sc_celltype = sc_celltype,
  method = 6,  # stereoscope
  python_path = "/path/to/conda/envs/stereoscope/bin/python"
)

# Requires: conda environment with cell2location
obj <- dec_celltype(
  object = obj,
  sc_data = sc_data,
  sc_celltype = sc_celltype,
  method = 7,  # cell2location
  python_path = "/path/to/conda/envs/cell2location/bin/python"
)

library(doParallel)

# Detect available cores
n_cores <- parallel::detectCores() - 1
cat("Using", n_cores, "cores\n")

# Register parallel backend
registerDoParallel(cores = n_cores)

# Run with parallel enabled
obj <- dec_cci_all(
  object = obj,
  if_doParallel = TRUE
)

# Clean up
stopImplicitCluster()

# Process cell types in batches
celltypes <- unique(obj@meta$rawmeta$celltype)
batch_size <- 5

for(i in seq(1, length(celltypes), batch_size)) {
  batch <- celltypes[i:min(i + batch_size - 1, length(celltypes))]
  
  for(ct_sender in batch) {
    for(ct_receiver in batch) {
      obj <- dec_cci(obj, ct_sender, ct_receiver)
    }
  }
  
  gc()  # Force garbage collection
}

# Load Visium data (via Seurat)
library(Seurat)
visium <- Load10X_Spatial("/path/to/visium/")

# Extract data for SpaTalk
st_data <- GetAssayData(visium, slot = "counts")
st_meta <- data.frame(
  spot = colnames(visium),
  x = visium@images$slice1@coordinates$col,
  y = visium@images$slice1@coordinates$row
)

# Create SpaTalk object (spot-based)
obj <- createSpaTalk(
  st_data = st_data,
  st_meta = st_meta,
  species = "Human",
  if_st_is_sc = FALSE,
  spot_max_cell = 10  # Expected cells per spot
)

# Slide-seq typically has ~10 cells per bead
obj <- createSpaTalk(
  st_data = slideseq_counts,
  st_meta = slideseq_coords,
  species = "Mouse",
  if_st_is_sc = FALSE,
  spot_max_cell = 10
)

# Single-cell resolution - no deconvolution needed
obj <- createSpaTalk(
  st_data = starmap_counts,
  st_meta = starmap_coords,
  species = "Mouse",
  if_st_is_sc = TRUE,
  spot_max_cell = 1,
  celltype = cell_annotations  # Pre-annotated cell types
)

# Get significant LR pairs
lr_results <- obj@lrpair
sig_pairs <- lr_results[lr_results$lr_co_ratio_pvalue < 0.05, ]

# Export to CSV
write.csv(sig_pairs, "significant_lr_pairs.csv", row.names = FALSE)

# Get TF activity scores
tf_results <- obj@tf

# Filter by score threshold
active_tfs <- tf_results[tf_results$score > 0.1, ]

# Combine LR and TF results into network
network <- merge(
  obj@lrpair,
  obj@tf,
  by = c("celltype_sender", "celltype_receiver")
)

# Export for Cytoscape
write.csv(network, "cci_network.csv", row.names = FALSE)

# Check gene overlap
st_genes <- rownames(obj@data$rawdata)
lr_genes <- unique(c(lrpairs$ligand, lrpairs$receptor))
overlap <- intersect(st_genes, lr_genes)
cat("Overlapping genes:", length(overlap), "\n")

# Reduce data size
obj <- createSpaTalk(
  st_data = st_data[, sample(ncol(st_data), 5000)],  # Subsample spots
  ...
)

# Check expression thresholds
# Lower the expression cutoff
obj <- find_lr_path(
  obj, lrpairs, pathways,
  min_exp = 0.01  # Lower threshold
)