Hello,
I took a 10x matrix from a collaborator and created a Seurat object. The issue I am having which I only realised when attempting to visualise my favourite genes, is that the original matrix has gene names in the format of "gene name - ensembl"
example :
Dcn-ENSMUSG00000019929.16, Inmt-ENSMUSG00000003477.5, Mfap4-ENSMUSG00000042436.12
So I can not search for my genes of interest unless I know what the corresponding ENSMUS following it is. Is there I way I can remove at this stage everything after (and including) the "-" so that I am left only with the gene name?
So convert "Dcn-ENSMUSG00000019929.16" to simply "Dcn".
Thank you :)
•
written
by
V
4 answers
Huh, you are quite right that my comment on jomo's answer doesn't work. Too used to SingleCellExperiment objects, where that does work. Anyway, there are two workarounds. The annoying one is to extract the counts from assay data and remake the Seurat object:
# Get counts
count.data <- GetAssayData(object = pbmc_small[["RNA"]], slot = "counts")
count.data <- as.matrix(count.data)
# Rename rows
rownames(count.data) <- unlist(map(strsplit(rownames(count.data),"-ENSMUS",fixed=T),1))
# Generate new Seurat object.
new.obj <- CreateSeuratObject(
count.data,
project = "SeuratProject",
assay = "RNA",
min.cells = 0,
min.features = 0,
names.field = 1,
names.delim = "_",
meta.data = NULL
)
Obviously not ideal. You'll have to redo the normalization. Other option is to convert to a SingleCellExperiment object, rename the rownames as above, then convert back to a Seurat object. This at least retains the normalization and any dimensionality reductions, metadata, clustering, etc.
I am afraid not. See: https://github.com/satijalab/seurat/issues/2617
Should just be:
rownames(SeuratObject)=unlist(map(strsplit(rownames(SeuratObject),"-ENSMUS",fixed=T),1))
Ensembl IDs are always unique, gene names are not, so making the row names not unique could cause problems. 10x usually makes a matrix such that row names are Enesmbl IDs, but software like Loupe and Seurat can work in gene names with that data. Your matrix really isn't made like that?
This should work.
RNA <- obj@assays$RNA
RNA@counts@Dimnames[[1]] = NEW_NAMES
RNA@data@Dimnames[[1]] = NEW_NAMES
rownames(RNA@meta.features) = NEW_NAMES
obj@assays$RNA <- RNA
Log in to answer this question.
-
How to use gsub inside mapply to assign names to Seurat objects
written by kousi31 •
-
CellRanger output more cells than specified using --force-cells? Why?
written by rohitsatyam102
-
How to load count matrix of scRNAseq data from cel-seq2 (96 well plates) protocol into a Seurat obj…
written by ashkan
-
how to subset Seurat object by gene names
written by leranwangcs •
-
Change location of 'fragments' in multiome Seurat file
written by V
-
Running scVelo on Seurat
written by V
-
Seurat FeaturePlot across stims
written by V
-
Seurat scRNAseq analysis: what is the mito gene list
written by wenbinm •
-
CellCycleScoring in Seurat with custom cell-cycle gene list
written by dtatarak •
-
Pick out genes based on protein function
written by V
Simply grep it while ensuring not to get false partial matches:
I've tried this and it only returns the row number of where my gene is :(
Hi, I am running into a similar problem in Seurat: My genes have the following format ENSEMBL-gene-biotype (e.g. ENSG00000000003_TSPAN6_ProteinCoding) so I run into the same problem whenever I need a reference (cell cycle scoring, label transfer,...). I would like to remove everything before and after the gene name (including the two "_"). And ideally I would like to store this original format so that I can later change it again (to know the biotype).
Thanks so much for your help!
Best Julia