Question: how to detect and get genes with outlier expression from RNA-seq
1
gravatar for rgescudero
9 months ago by
rgescudero30
Spain
rgescudero30 wrote:

I am interested in searching for genes having outlier overexpression from RNA-seq samples (hundreds). I have checked for outlier analysis tools, but most are focused on detecting outlier samples instead of outlier genes. Also, I would like to known which expression values (FPKM, RPKM, TPM, CPM, etc) are better to detect outlier expression. I would thank any help on tools, recommendations, etc.

ADD COMMENTlink modified 9 months ago by h.mon24k • written 9 months ago by rgescudero30

I don't completely understand what you want to do:

Do you want genes showing consistent higher or lower expression compared to the other genes - that is, for a particular hypothetical gene, all samples have higher or lower expression compared to the other genes. edit: strictly speaking, these wouldn't be outliers.

Or do you want genes that, in one or a few samples, have a higher or lower expression compared to the same gene from other samples?

ADD REPLYlink modified 9 months ago • written 9 months ago by h.mon24k

Second option: genes having higher outlier expression in a few samples compared with same genes in other samples. I'm interested in outlier overexpressed genes that may exist in some samples. It would be an unsupervised analysis because I do not have a priori information about which samples would have outlier expression of which genes.

Any help will be welcome

ADD REPLYlink written 9 months ago by rgescudero30
0
gravatar for h.mon
9 months ago by
h.mon24k
Brazil
h.mon24k wrote:

DESeq2 uses Cook's distance to identify count outiers, check this vignette. The DESeq() function automatically replaces outlier counts and refits the GLM, however the vignette shows how you can identify outlier counts. Keep in mind you need a good number of samples per group to have reliable detection of outliers - which probably won't be a problem in your case.

ADD COMMENTlink written 9 months ago by h.mon24k

Many thanks. I'll check it

ADD REPLYlink written 9 months ago by rgescudero30
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1933 users visited in the last hour