Question: Have Salmon output tabular files from Galaxy; need to import into R for DESeq2
0
gravatar for Dennis
11 months ago by
Dennis0
Dennis0 wrote:

Hi there,

I have Salmon output files from 9 samples - 3 biological replicates from 3 different treatment groups.

I need to perform DE analysis on these Salmon output files.

Firstly, from what I understand, it is uncouth to perform DE analysis on TPM values - I still don't understand why, since the values are normalized for library depth and transcript length, aren't these the best values to perform DE analysis on?

Secondly, is there a tutorial for importing Salmon files using tximport package into R that's written in layman English. I tried the Galaxy-based DEseq2 and picked the TPM option, but was informed by one of the moderators that the online tool has never worked properly with TPM values. So, I am left at the mercy of R, which I am not too friendly with.

Any help would be much appreciated.

Any recommendation on good DESeq2 tutorials in R would be much appreciated.

Cheers, Dennis

rna-seq tpm deseq2 salmon • 637 views
ADD COMMENTlink modified 11 months ago by WouterDeCoster34k • written 11 months ago by Dennis0
0
gravatar for WouterDeCoster
11 months ago by
Belgium
WouterDeCoster34k wrote:

Firstly, from what I understand, it is uncouth to perform DE analysis on TPM values - I still don't understand why, since the values are normalized for library depth and transcript length, aren't these the best values to perform DE analysis on?

TPM is a normalization, indeed, but it's a very simple one. DESeq2 is smarter, will do a better normalization, and therefore needs the raw data for properly estimating the dispersion. If you already "pre-normalize" a part of the information is lost (e.g. library size).

Secondly, is there a tutorial for importing Salmon files using tximport package into R that's written in layman English. I tried the Galaxy-based DEseq2 and picked the TPM option, but was informed by one of the moderators that the online tool has never worked properly with TPM values. So, I am left at the mercy of R, which I am not too friendly with.

I find this tutorial pretty easy to follow, and I really don't like R: Importing transcript abundance datasets with tximport.
This bioconductor workflow RNA-seq workflow: gene-level exploratory analysis and differential expression is also a good starting point for DESeq2.

ADD COMMENTlink written 11 months ago by WouterDeCoster34k
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: 956 users visited in the last hour