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
ADD COMMENT
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)

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

ADD COMMENT
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.

ADD REPLY
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,...

ADD REPLY
0
Entering edit mode

Thank you fracarb8 this is very helpful!

ADD REPLY

Login before adding your answer.

Traffic: 2204 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6