SUPPA2检测差异可变剪切

SUPPA2原理

假设我们已经了解了可变剪切类型:传送门，那么我们有许多检测可变剪切的软件，比方说rMAT，miso，SUPPA。而今天重点介绍下SUPPA的基本原理和用法

SUPPA2检测可变剪切是基于转录本定量来判断的，SUPPA2先是利用转录本定量的结果计算每一个样本的PSI(inclusion values)

A PSI value was given according to the ratio of the long form on total form present (short form and long form) to characterize inclusion of exon, differential splice-site choice, intron retention, etc

接下来SUPPA2考虑两个分布（或者说是变化量）：一个是每一个生物学重复之间的ΔPSI；另一个是不同条件（处理）之间的ΔPSI

首先考虑生物学重复之间PSI：

1. 其中a代表某一个基因的不同转录本
2. r代表不同生物学重复
3. Rc代表不同处理，如果是两个重复则可用1，2表示

那么每一个生物学重复之间的ΔPSI既是任意两个生物学重复某一个基因的不同PSI的差值，某可变剪切事件相同条件（处理）下各生物学重复的平均PSI在式子中用Erep表示

再来考虑不同条件（处理）之间的PSI：

1. 其中a代表某一个基因的不同转录本
2. r代表不同生物学重复
3. Rc代表不同处理，如果是两个重复则可用1，2表示

那么不同条件（处理）之间的ΔPSI既是任意两个条件（处理）中某一个基因的不同PSI的差值，某可变剪切事件在两个处理下总的平均PSI在式子中用Econd表示

我们可以计算每一个生物学重复之间的PSI与不同条件（处理）之间的平均PSI的差值，计算其最小值，Ei,rep代表第i个重复的PSI值

那么判断差异可变剪切的依据是组间的ΔPSI是否大于组内的ΔPSI

用法

该软件的参考文档：github地址

1.判断可变剪切事件

python3.4 suppa.py generateEvents -i  -o  -f ioe -e 
#-e为可变剪切类型
##SE: Skipping exon (SE)
##SS: Alternative 5' (A5) or 3' (A3) splice sites (generates both)
##MX: Mutually Exclusive (MX) exons
##RI: Retained intron (RI)
##FL: Alternative First (AF) and Last (AL) exons (generates both)

这一步输入GTF文件

input文件

得到文件如下：

ioi格式

2.计算PSI

python3.4 suppa.py psiPerIsoform -g  -e  -o 

如果是local AS 事件，输出文件为：

即每一个sample位置上的PSI值

3.计算差异PSI

python3.4 suppa.py diffSplice --method  --input  --psi   --tpm   --area <1000> --lower-bound <0.05> -gc -o 

对于local AS 事件：

1. Cond1_Cond2_dPSI: Event PSI difference (ΔPSI) between Cond1 and Cond2 (ΔPSI = PSI_2 - PSI_1).
2. Cond1_Cond2_pvalue: Significance of the difference of PSI between Cond1 and Cond2

4.计算平均PSI

input文件

python3.4 suppa.py clusterEvents --dpsi  --psivec  --sig-threshold <0.05> --eps <0.05> --min-pts <20> --groups <1-3,4-6> -o  

输出为：

参考：
《SUPPA2: fast, accurate, and uncertainty-aware differential splicing analysis across multiple conditions》