This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Forcing Seurat to include certain genes during integration

Is there a way of forcing Seurat to include a certain set of genes during integration? I think I saw somewhere that this was possible but after looking through the documentation for the different functions used during integration, and also at various questions on this site, I haven't been able to find anything.

integration seurat
• 2,545 views • link
written by Aaron • Aaron

1 answer

You can pass the list of genes to FindIntegrationAnchors

# find common genes for the integration 
features <- SelectIntegrationFeatures(your_list_of_seurat_objects)

# add your genes
genes2add <- c("CD8A","CD14",...)
features <- c(genes2add, features)
features <- unique(features)

# Integrate the datasets using the new anchors
anchors <- FindIntegrationAnchors(object.list = your_list_of_seurat_objects , anchor.features = features)
exp_Integrated <- IntegrateData(anchorset = anchors)

You can then DefaultAssay(exp_Integrated) <- "integrated", and continue with your analysis (PCA, UMAP, clustering,...)

• 0 views • link
written by fracarb8 fracarb8

This is interesting/cool. fracarb8 I am assuming you have done this/do this (include a certain set of genes during integration). What effect have you seen from doing this? I kind of have an ambiguous notion of what could happen.

• 0 views • link
written by Pratik Pratik

Actually, I have never done it. I have removed some genes (e.g. genes associated to contamination), and modify the values of VariableFeatures for training/testing classification models, but I never touched them during integration. The effects of the alteration can vary greatly, mainly because there are many factors involved: the biological system, the number of genes added/removed, the experiment design, the questions you are trying to answer,...

• 0 views • link
written by fracarb8 fracarb8

Log in to answer this question.

More posts like this