Hello,
I am performing RNA-Seq analysis on 6 samples (3 infected & 3 control). I aligned my samples to GENCODE's release 41 primary assembly and quantified my fragments using RSubread::featurecounts(). So now I have a counts matrix and have TMM normalized my counts in EdgeR. I removed genes that are lowly expressed and would like to remove Ribosomal genes.
I am thinking that I could just remove rows in my matrix that correspond to the 62 genes classified as Ribosomal at HGNC's website : https://www.genenames.org/data/genegroup/#!/group/1379. Is this an effective method of removing ribosomal rna counts from my analysis?
Thank You, Sara
1 answer
In my experience, this is a most reliable way of detecting rRNA sequences:
Log in to answer this question.
-
Comparing P-Values & FDR adjusted p-values between RNA-Seq experiments?
written by Saran •
-
Flitering Genes in RNA-Seq causes more significant FDR Adjusted P-Val DEGs?
written by Saran •
-
is normalize RNA seq raw counts used correctly?
written by kng •
-
edgeR normalization
written by rheab1230 •
-
limma voom -- all genes after filtering are differentially expressed
written by wiscoyogi •
-
Remove HB genes and recalculate read depth
written by Cece •
-
Cleaning the expression matrix for the purpose of deconvolution of RNA-Seq data
written by jane.merlevede •
-
Analysis of GEO dataset normalized by FPKM
written by GSAENZDEPIP •
-
how to obtain a normalised expression matrix from read counts and scaling factors
written by user31888
-
Effective length in normalisation for reads counts
written by Xiaokang •
You'll want to remove reads mapping to rRNA before normalizing with EdgeR.
Yes, that is true, I will re-run normalization when I figure out how to remove rRNA genes.