Can I use abundance from Tximport to compare the expression level of transcript in same condition?
1
0
Entering edit mode
10 weeks ago
aimanbarki ▴ 20

I wanted to compare different transcript of the same gene to get an idea which one is expressed more than other.

For this I have 3 samples. I did the alignment using Salmon and then I tximport to pull out abundances.

Can I take a average of 3 samples and make a bar graph to see the transcript expression of a gene?

enter image description here

Salmon Tximport • 487 views
ADD COMMENT
1
Entering edit mode

Can I

Yes. Should you? Probably not. There must be better, statistically more robust ways of quantifying based on biological replicates.

ADD REPLY
0
Entering edit mode

Can you please mention any of it. I could not just find any good example from literature as I am only interested in knowing transcript of one gene in one condition (normal)

ADD REPLY
0
Entering edit mode

Like ntsopoul said, explore DESeq2.

ADD REPLY
0
Entering edit mode

DESeq2 is for Differential gene analysis. I also tried DEXSeq (DEU) and DRIMSeq DTU) but my issue is the following:

  1. I have only one condition
  2. I only wanted to know the relative expression of transcript of one gene in normal sample.
ADD REPLY
1
Entering edit mode

I have only one condition

There is no "differential" then. RNAseq is a relative measure, so comparing single-condition replicates using trivial mathematical operations will be inherently flawed. Fold change metrics make sense, but any normalization or aggregation you do will only serve one purpose: look good in a bar plot; it won't be robust data you can base an investigation on. I deal with data following your exact premise (RNA-seq on 100+ baseline samples, no other condition) and it's difficult to justify any kind of arithmetic we do to say sample X has n times the expression for gene G1 compared to sample Y. We can only say that sample X seems to express more of gene G1 than sample Y.

ADD REPLY
0
Entering edit mode

Thanks, I am not doing any differential. But, Can we say that 6th transcript is relatively more abundance in our sample than other transcripts of that gene?

ADD REPLY
1
Entering edit mode

Oh, you can just use TPM to do that. Within sample comparison is a completely different (and much easier) game.

ADD REPLY
0
Entering edit mode

So can I average 6 replicates??

ADD REPLY
1
Entering edit mode

No. I think you can compare average TPMs as long as you explicitly mention that anywhere you report metrics but don't average counts.

ADD REPLY
0
Entering edit mode

Rob Please, Can you or people from your group answer . As I am using one of your paper as my reference.

ADD REPLY
0
Entering edit mode
10 weeks ago
ntsopoul ▴ 50

No, you should not do that since the count matrix is not normalized. Also, the statistics for getting a pValue requires consideration of all genes.

If you have reached the step where you can retrieve the count data via tximport you should use one of the many differential gene expression analysis tools. I recommend deseq2, since the vignette describes exactly how to proceed from tximport onwards

https://bioconductor.org/packages/release/bioc/html/DESeq2.html

Good luck!

ADD COMMENT
0
Entering edit mode

Salmon is giving me TPM at transcript level. So, What I understand from that TPM is normalized for sequencing depth and Transcript length. I am using Tximport to generate a matrix to get those abundance/TPM. According to following link enter image description here

ADD REPLY

Login before adding your answer.

Traffic: 1604 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6