Question: Alternative splicing: explanation about delta psi values
gravatar for user31888
2.0 years ago by
United States
user3188870 wrote:

I am a bit confused.

I compared differentially expressed spliced transcripts between a sample 1 and a sample 2 using MISO algorithm.

The final (filtered) result table notably mentions the following info (this is an example for a transcript with 3 isoforms):

. delta PSI: -0.67,0.19,-0.87

. isoforms: isoform1,isoform2,isoform3 (isoform 1 is the inclusive form (with all exons); isoforms 2 and 3 are spliced)

According to MISO paper:

. PSI is the fraction of the inclusive transcript in a sample (i.e. not spliced if I understand well)

. delta PSI = PSI(sample1) - PSI(sample2)

Thus, for example:

. delta PSI(isoform1) = -0.67 means there are more inclusive transcript in sample 2 than sample 1,

. delta PSI(isoform2) = 0.19 means there are more inclusive transcript in sample 1 than sample 2, but how is it related to isoform2 which is known to be spliced? Does MISO report only the spliced isoforms in its output?

alternative splicing miso psi • 2.0k views
ADD COMMENTlink modified 2.0 years ago by kristoffer.vittingseerup2.4k • written 2.0 years ago by user3188870

Can I ask how you filtered your output when your sample has more than two isoforms? I'm aware that the documentation allows you to filter using the filter_events script, but thats only for two isoform events. You can filter on bayes-factor on >2 isoform events but this still outputs an note/error for me: "filter_events: filtering MISO pairwise comparison output.

Note: This utility is only works on MISO output for two-isoform event annotations.

Error: is only defined for MISO output on two-isoform alternative events. Found a non-two isoform event: ENSMUSG00000032469"

ADD REPLYlink written 10 months ago by aarin.jones0
gravatar for kristoffer.vittingseerup
2.0 years ago by
European Union
kristoffer.vittingseerup2.4k wrote:

I must admit I'm also a bit confused. Note sure how to interpret PSI values for isoforms (instead of an exon) since it have to be different than Isoform Fractions (as they are not summing to zero).

You could take a look at IsoformSwitchAnalyzeR to see whether that does what you want.

ADD COMMENTlink written 2.0 years ago by kristoffer.vittingseerup2.4k

I agree. It does not make any sense for isoform-centric analysis. I will definetely try your program. Does it require technical/biological replicates?

ADD REPLYlink modified 2.0 years ago • written 2.0 years ago by user3188870

Yes (like everything else where you want to do statistics). It can however be run without but then the results are not trustworthy.

ADD REPLYlink written 10 months ago by kristoffer.vittingseerup2.4k
Please log in to add an answer.


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