ncRNA, pseudogenes and lncRNA genes in highly variable features ( scRNA)
Is it okay to include those genes starts with "Gm" to variable features or should I remove them from variable features?
single-cell
• 1,690 views
•
link
updated
by
dsull
•
written
by
carolofharvest •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Correct way of doing enrichment analysis (ORA)
written by carolofharvest •Hello, Considering that MSigDB databases such as GO and Hallmark contain gene sets for different pathways, these gene sets may include genes that regulate the …
-
Is it reasonable to keep only protein-coding genes in highly variable features in scRNA-seq data?
written by carolofharvest •Hello, In my human cell line data I have a lot of psedo-genes and Seurat::FindVariableFeature() function also includes them into highly variable features. I noticed …
-
Removing features from RNA assay in seurat
written by Varun GuptaI have a seurat object which looks like this: An object of class Seurat 31061 features across 6471 samples within 2 assays Active assay: RNA …
-
Keep seeing same degs everytime after filtering and clustering (scRNA)
written by carolofharvest •Hi. I followed the standard Seurat workflow on my data and noticed a cluster with low UMI counts. When I performed differential expression analysis (DEG) …
-
Highest variable features in single cell data
written by carolofharvest •I used the 'FindVariableFeatures' function from the Seurat package to identify variable features, but some of the genes appearing in the results are only expressed …
-
Single cell
written by shiviagarwalpharma •HI Guys, I am analyzing a single cell data containing 22815 features (genes) across 1347 samples (cells) within 1 assay. I am following the Seurat …
-
topGO enrichment: Should the "Gene universe" also contain genes with no GO term?
written by YOUSEUFS •I have a set of interesting genes with GO terms, I want to find out which GO terms are enriched compared to a background set …
-
ENSEMBL annotation of genes start with "Gm" end with "Rik"
written by SMILEDears, When I did RNA-Seq and ChIP-Seq downstream analysis with GENE annotation of ENSEMBL. The final results are flooded with genes such as Gm's and …
-
Threshold of removing genes with low counts ?
written by hafiz.talhamalikWhat should be the threshold of removing genes with low counts in RNA-seq data ? i have removed all genes with sum of counts 0, …
-
How to deal with too-many-zeros genes from TCGA?
written by Xiaokang •I'm using the RNA-Seq data from TCGA of Level 3 with gene expression levels of normal and tumor patients to do classification. But found that, …
Which genome are you referring to?
Hi. I only have count matrix without empty droplets. I really dont know what they did in previous steps in cellranger etc. But it is a mouse data.
Single-cell RNA-seq is a lot of "exploratory" work; you can explore both approaches and see.
If a strong signal is being driven by a noncoding RNA, it might be something interesting to look into further. On the other hand, many ncRNAs are poorly annotated and/or have repetitive elements which causes reads to misalign to them.
If you're using the default cellranger genome, you can just use all those genes and proceed forward with your analysis like everyone else does, unless you encounter something funky downstream.
Hi,
I used those Gm42418 , Gm26917 and AY036118 genes in PercentageFeatureSet() function. Some cells have %60 of those genes. And also they are forming a single cluster in umap. I cant find a way to annotate them. In DEG results ın which I compared this cluster to other cluster, pct1 and pct2 are almost same .but this cluster expresses those genes highly.
What do you mean by "funky" ?
OK, good exploratory work -- yeah, you encountered the "funky" stuff I was talking about. Those are genes that have repetitive elements in them, so ribosomal RNAs present in your data will misalign to them.
I would remove them.
Hi, Would you remove genes or cells ?
Remove genes. After all, it’s those genes that are problematic.