How to label metadata of Seurat object based gene read count

I have a scRNA-seq dataset that was mapped to both the mouse and mCMV genomes. All the genes from mCMV begin with "392XXXX" and there are 161 of them

For each cell in the dataset, I want to add a metadata identifier if at least one read was detected from at least one of these genes.

Any tips greatly appreciated!

thanks

seurat scrna-seq

• 938 views • link

updated 4.4 years ago by rpolicastro

• written 4.4 years ago by cook.675 •

Untested but should work in theory. If you get an error you may need to convert the counts from a sparse to dense matrix with as.matrix.

library("stringr")
library("Seurat")

mCMV_data <- GetAssayData(seu, assay="RNA", slot="counts")
mCMV_data <- mCMV_data[str_starts(rownames(mCMV_data), "392"), ]

seu$mCMVExpr <- rowSums(mCMV_data > 0) > 1

• 0 views • link

written 4.4 years ago by rpolicastro

0 answers

No answers yet.

More posts like this

Creating a signature matrix using CibersortX
written 2.0 years ago by Aspire

(1) It is possible to create a signature matrix with CibersortX from single-cell data, or alternatively from RNA-seq of sorted cell populations. Is there an …
featureCount Error "No paired-end reads were detected in paired-end read library"
written 2.7 years ago by cook.675 •

I created a combined mm39 and MHV-A59 (Viral) reference genome and aligned my paired end reads using STAR with the following input commands: STAR --runMode …
Constructing GRNs from multimodal data?
written 5.0 years ago by kgoss •

Hello all, I have a multimodal single cell dataset consisting of scRNA-seq and scATAC-seq. I'm wanting to construct gene regulatory networks, and was looking into …
correlation of WGBS dataset and RNAseq dataset
written 5.3 years ago by frylhc •

Hi, everyone. I am analyzing the correlation of WGBS dataset and RNAseq dataset. Generally, more than one methylated CpGs could be detected in the transcriptional …
customized gene not show in the output from cellranger count even if I have added to the reference …
written 5.3 years ago by leranwangcs •

Hi, I have 3 genes that I need to add to the cellranger mouse reference genome database. I followed the "Add a Marker Gene to …
Barnyard analysis advice
written 6.0 years ago by cf556 •

I have a scRNA dataset with both human and mouse cells. Due to the upstream experiment, there is no way to physically separate the cells …
How to extract a specific gene sequence for a specific cell from scRNA-seq fastq file
written 6.7 years ago by cook.675 •

I have a scRNA seq experiment where I would like to look at the raw sequencing data for one specific gene. We have 3' ends, …
DEG data per-cell
written 6.8 years ago by cook.675 •

I have 2 scRNA-seq data sets which come from different treatment conditions (Drug vs. Control) I have been using Seurat and running the data through …
QC filtering by Count depth
written 6.8 years ago by cook.675 •

Im working with a scRNA data set and trying out some different filtering parameters other than using the Seurat vignettes. Specifically I am working at …
Summing data from 3 different parts of an R object at once
written 6.8 years ago by cook.675 •

I have 3 replicate samples from a scRNA seq experiment; call them V1, V2, V3. I merged them into one large object "Object". When I …