Forcing Seurat to include certain genes during integration
1
1
Entering edit mode
8 months ago
Aaron ▴ 10

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 • 595 views
4
Entering edit mode
8 months ago
fracarb8 ▴ 830

You can pass the list of genes to FindIntegrationAnchors

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

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,...)

1
Entering edit mode

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.

2
Entering edit mode

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
Entering edit mode

Thank you fracarb8 this is very helpful!