# 20191230对AS类型进行定义 前面通过Blast的方法对可变剪切在不同亚基因组间的保守性进行了分类，然后对每一类基因定义了一个概念，这个概念是紧紧围绕多倍化过程中，可变剪切在二倍中同源基因和四倍体同源基因的变化情况 ## 第一类比较保守的AS ***conservation AS pattern*** ![conservation AS pattern](https://43423.oss-cn-beijing.aliyuncs.com/img/20191230164457.png) 筛选条件 * 匹配长度大于200bp * 相似度大于90% * 覆盖度超过85% 代码 ```bash ## 在所有基因组中都保守的IR对应的bed awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_D.*/{print $0}' Dt_query_blast.txt |awk '$2~/Ghir_A/{a[$1][0]+=1}$2~/evm/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][0]>=1&&a[i][2]>=1&&a[i][1]>=1)print i}}' >converse/1_1 ## 都不发生IR事件的基因 cat 1_1.bed 2_2.bed 4_2.bed 1_3.bed 2_3.bed 5_1.bed 1_4.bed 3_2.bed 5_2.bed 2_1.bed 4_1.bed A2.bed At.bed D5.bed Dt.bed|awk '{print $NF}'|grep "Ghir_D" |awk -F "_" '$1~/Ghir/{print $1"_"$2}'|sort |uniq |xargs -I {} grep {} ../../../../GhDt_Gr_GhAt_Ga_end_noScaffold >tmp_Dt cat 1_1.bed 2_2.bed 4_2.bed 1_3.bed 2_3.bed 5_1.bed 1_4.bed 3_2.bed 5_2.bed 2_1.bed 4_1.bed A2.bed At.bed D5.bed Dt.bed|awk '{print $NF}'|grep "Ghir_A" |awk -F "_" '$1~/Ghir/{print $1"_"$2}'|sort |uniq |xargs -I {} grep {} ../../../../GhDt_Gr_GhAt_Ga_end_noScaffold >tmp_At cat 1_1.bed 2_2.bed 4_2.bed 1_3.bed 2_3.bed 5_1.bed 1_4.bed 3_2.bed 5_2.bed 2_1.bed 4_1.bed A2.bed At.bed D5.bed Dt.bed|awk '{print $NF}'|grep "Gor" |awk -F "_" '{print $1}'|sort |uniq |xargs -I {} grep {} ../../../../GhDt_Gr_GhAt_Ga_end_noScaffold >tmp_D5 cat 1_1.bed 2_2.bed 4_2.bed 1_3.bed 2_3.bed 5_1.bed 1_4.bed 3_2.bed 5_2.bed 2_1.bed 4_1.bed A2.bed At.bed D5.bed Dt.bed|awk '{print $NF}'|grep "evm" |awk -F "_" '{print $1}'|sort |uniq |xargs -I {} grep {} ../../../../GhDt_Gr_GhAt_Ga_end_noScaffold >tmp_A2 cat tmp* |sort |uniq |cat - ../../../../GhDt_Gr_GhAt_Ga_end_noScaffold |sort |uniq -u >../converseBed/1_2 ## 只在A2与At中发生 awk '$7>200&&$8>90&&$9>85&&$2!~/evm.*/{print $0}' A2_query_blast.txt |awk '$2~/Ghir_A/{a[$1][0]+=1}$2~/Ghir_D/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][1]==0&&a[i][0]>=1&&a[i][2]==0)print i}}' >converse/1_3 ## 只在D5与Dt中发生 awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_D.*/{print $0}' Dt_query_blast.txt |awk '$2~/Ghir_A/{a[$1][0]+=1}$2~/evm/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][0]==0&&a[i][2]==0&&a[i][1]>=1)print i}}' >converse/1_4 ``` 做成bed文件 ```bash ## 在所有基因组中都保守的IR对应的bed awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_D.*/{print $0}' Dt_query_blast.txt |awk '$2~/Ghir_A/{a[$1][0]+=1}$2~/evm/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][0]>=1&&a[i][2]>=1&&a[i][1]>=1)print i}}' |xargs -I {} grep {} Dt_query_blast.txt| awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_D.*/{print $1"\n"$2}'|awk -F "_" '$1~/Ghir/{print $3"_"$4,$5,$6,$8,$0}$1!~/Ghir/{print $2,$3,$4,$6,$0}' OFS="\t" |sort |uniq >converseBed/1_1.bed ## 只在A2与At中发生 awk '$7>200&&$8>90&&$9>85&&$2!~/evm.*/{print $0}' A2_query_blast.txt |awk '$2~/Ghir_A/{a[$1][0]+=1}$2~/Ghir_D/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][1]==0&&a[i][0]>=1&&a[i][2]==0)print i}}'|xargs -I {} grep {} A2_query_blast.txt |awk '$7>200&&$8>90&&$9>85&&$2!~/evm.*/{print $1"\n"$2}'|awk -F "_" '$1~/Ghir/{print $3"_"$4,$5,$6,$8,$0}$1!~/Ghir/{print $2,$3,$4,$6,$0}' OFS="\t" |sort |uniq >converseBed/1_3.bed ## 只在D5与Dt中发生 awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_D.*/{print $0}' Dt_query_blast.txt |awk '$2~/Ghir_A/{a[$1][0]+=1}$2~/evm/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][0]==0&&a[i][2]==0&&a[i][1]>=1)print i}}' |xargs -I {} grep {} Dt_query_blast.txt| awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_D.*/{print $1"\n"$2}'|awk -F "_" '$1~/Ghir/{print $3"_"$4,$5,$6,$8,$0}$1!~/Ghir/{print $2,$3,$4,$6,$0}' OFS="\t" |sort |uniq >converseBed/1_4.bed ``` ## 第二类多倍化后四倍体中AS发生丢失 ***AS loose*** * ***symmetry loose*** * ***At asymmetry loose*** * ***Dt asymmetry loose*** ![AS loose](https://43423.oss-cn-beijing.aliyuncs.com/img/20191230164455.png) 代码 ```bash ## Dt中AS 丢失 awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_A.*/{print $0}' At_query_blast.txt |awk '$2~/evm/{a[$1][0]+=1}$2~/Ghir_D/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][2]==0&&a[i][0]>=1&&a[i][1]>=1)print i}}' >converse/2_1 ## At中AS 丢失 awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_D.*/{print $0}' Dt_query_blast.txt |awk '$2~/Ghir_A/{a[$1][0]+=1}$2~/evm/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][0]==0&&a[i][2]>=1&&a[i][1]>=1)print i}}' >converse/2_2 ## At与Dt都丢失 awk '$7>200&&$8>90&&$9>85&&$2!~/evm.*/{print $0}' A2_query_blast.txt |awk '$2~/Ghir_A/{a[$1][0]+=1}$2~/Ghir_D/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][0]==0&&a[i][1]>=1&&a[i][2]==0)print i}}' >converse/2_3 ``` 做成bed文件 ```bash ## Dt中AS 丢失 awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_A.*/{print $0}' At_query_blast.txt |awk '$2~/evm/{a[$1][0]+=1}$2~/Ghir_D/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][2]==0&&a[i][0]>=1&&a[i][1]>=1)print i}}' |xargs -I {} grep {} At_query_blast.txt|awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_A.*/{print $1"\n"$2}'|awk -F "_" '$1~/Ghir/{print $3"_"$4,$5,$6,$8,$0}$1!~/Ghir/{print $2,$3,$4,$6,$0}' OFS="\t" |sort |uniq >converseBed/2_1.bed ## At中AS 丢失 awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_D.*/{print $0}' Dt_query_blast.txt |awk '$2~/Ghir_A/{a[$1][0]+=1}$2~/evm/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][0]==0&&a[i][2]>=1&&a[i][1]>=1)print i}}' |xargs -I {} grep {} Dt_query_blast.txt|awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_D.*/{print $1"\n"$2}'|awk -F "_" '$1~/Ghir/{print $3"_"$4,$5,$6,$8,$0}$1!~/Ghir/{print $2,$3,$4,$6,$0}' OFS="\t" |sort |uniq >converseBed/2_2.bed ## At与Dt都丢失 awk '$7>200&&$8>90&&$9>85&&$2!~/evm.*/{print $0}' A2_query_blast.txt |awk '$2~/Ghir_A/{a[$1][0]+=1}$2~/Ghir_D/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][0]==0&&a[i][1]>=1&&a[i][2]==0)print i}}'|xargs -I {} grep {} A2_query_blast.txt | awk '$7>200&&$8>90&&$9>85&&$2!~/evm.*/{print $1"\n"$2}'|awk -F "_" '$1~/Ghir/{print $3"_"$4,$5,$6,$8,$0}$1!~/Ghir/{print $2,$3,$4,$6,$0}' OFS="\t" |sort |uniq >converseBed/2_3.bed ``` ## 第三类多倍化后四倍体AS重新获得 ***AS gain*** * ***symmetry gain*** * ***At asymmetry gain*** * ***Dt asymmetry gain*** ![AS gain](https://43423.oss-cn-beijing.aliyuncs.com/img/20191230164452.png) 代码 :warning:单个基因组存在的事件得再converseBed目录下运行 ```bash ## 只有At 获得新的AS事件 cat 1_1.bed 1_3.bed 2_1.bed 3_2.bed 4_1.bed 4_2.bed 5_2.bed |awk -F "\t" '$5~/Ghir_A/{print ">"$5}'|cat - ../At_intron_junction.fasta|grep ">"|sort|uniq -u|sed 's/>//g' >../converse/3_1 ## At与Dt 获得同样的AS事件 awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_D.*/{print $0}' Dt_query_blast.txt |awk '$2~/Ghir_A/{a[$1][0]+=1}$2~/evm/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][1]==0&&a[i][2]==0&&a[i][0]>=1)print i}}' >converse/3_2 ## 只有Dt 获得新的AS事件 cat 1_1.bed 1_4.bed 2_2.bed 3_2.bed 4_1.bed 4_2.bed 5_1.bed|awk -F "\t" '$5~/Ghir_D/{print ">"$5}' |cat - ../Dt_intron_junction.fasta|grep ">"|sort|uniq -u|sed 's/>//g' >../converse/3_3 ``` 做成bed文件 ```bash ## At与Dt 获得同样的AS事件 awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_D.*/{print $0}' Dt_query_blast.txt |awk '$2~/Ghir_A/{a[$1][0]+=1}$2~/evm/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][1]==0&&a[i][2]==0&&a[i][0]>=1)print i}}'|xargs -I {} grep {} Dt_query_blast.txt|awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_D.*/{print $1"\n"$2}'|awk -F "_" '$1~/Ghir/{print $3"_"$4,$5,$6,$8,$0}$1!~/Ghir/{print $2,$3,$4,$6,$0}' OFS="\t" |sort |uniq >converseBed/3_2.bed ``` ## 第四类四倍体中两个亚基因组有趋同效应 ***convergence effect*** * ***convergence to gain*** * ***convergence to loose*** ![convergence effect](https://43423.oss-cn-beijing.aliyuncs.com/img/20191230164459.png) 代码 ```bash ## D5中不存在 awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_A.*/{print $0}' At_query_blast.txt |awk '$2~/evm/{a[$1][0]+=1}$2~/Ghir_D/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][2]>=1&&a[i][0]>=1&&a[i][1]==0)print i}}' >converse/4_1 ## A2中不存在 awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_D.*/{print $0}' Dt_query_blast.txt |awk '$2~/Ghir_A/{a[$1][0]+=1}$2~/evm/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][1]>=1&&a[i][2]==0&&a[i][0]>=1)print i}}' >converse/4_2 ## 只有A2中有 cat 1_1.bed 1_3.bed 2_1.bed 2_2.bed 2_3.bed 4_1.bed 5_1.bed|awk -F "\t" '$5~/evm/{print ">"$5}'|cat - ../A2_intron_junction.fasta|grep ">"|sort|uniq -u|sed 's/>//g' >../converse/4_3 ## 只有D5中有 cat 1_1.bed 1_4.bed 2_1.bed 2_2.bed 2_3.bed 4_2.bed 5_2.bed|awk -F "\t" '$5~/Gor/{print ">"$5}'|cat - ../D5_intron_junction.fasta|grep ">"|sort|uniq -u|sed 's/>//g' >../converse/4_4 ``` 做成bed文件 ```bash ## D5中不存在 awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_A.*/{print $0}' At_query_blast.txt |awk '$2~/evm/{a[$1][0]+=1}$2~/Ghir_D/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][2]>=1&&a[i][0]>=1&&a[i][1]==0)print i}}' |xargs -I {} grep {} At_query_blast.txt|awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_A.*/{print $1"\n"$2}'|awk -F "_" '$1~/Ghir/{print $3"_"$4,$5,$6,$8,$0}$1!~/Ghir/{print $2,$3,$4,$6,$0}' OFS="\t" |sort |uniq >converseBed/4_1.bed ## A2中不存在 awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_D.*/{print $0}' Dt_query_blast.txt |awk '$2~/Ghir_A/{a[$1][0]+=1}$2~/evm/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][1]>=1&&a[i][2]==0&&a[i][0]>=1)print i}}'|xargs -I {} grep {} Dt_query_blast.txt|awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_D.*/{print $1"\n"$2}'|awk -F "_" '$1~/Ghir/{print $3"_"$4,$5,$6,$8,$0}$1!~/Ghir/{print $2,$3,$4,$6,$0}' OFS="\t" |sort |uniq >converseBed/4_2.bed ``` ## 第五类 AS事件在不同亚基因组间保守 ```bash ## 只在A2与Dt中保守 awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_D.*/{print $0}' Dt_query_blast.txt |awk '$2~/Ghir_A/{a[$1][0]+=1}$2~/evm/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][0]==0&&a[i][2]>=1&&a[i][1]==0)print i}}' >converse/5_1 ## 只在D5与At中保守 awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_A.*/{print $0}' At_query_blast.txt |awk '$2~/evm/{a[$1][0]+=1}$2~/Ghir_D/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][2]==0&&a[i][1]>=1&&a[i][0]==0)print i}}' >converse/5_2 ``` 做成bed文件 ```bash ## 只在A2与Dt中保守 awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_D.*/{print $0}' Dt_query_blast.txt |awk '$2~/Ghir_A/{a[$1][0]+=1}$2~/evm/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][0]==0&&a[i][2]>=1&&a[i][1]==0)print i}}'|xargs -I {} grep {} Dt_query_blast.txt|awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_D.*/{print $1"\n"$2}'|awk -F "_" '$1~/Ghir/{print $3"_"$4,$5,$6,$8,$0}$1!~/Ghir/{print $2,$3,$4,$6,$0}' OFS="\t" |sort |uniq >converseBed/5_1.bed ## 只在D5与At中保守 awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_A.*/{print $0}' At_query_blast.txt |awk '$2~/evm/{a[$1][0]+=1}$2~/Ghir_D/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][2]==0&&a[i][1]>=1&&a[i][0]==0)print i}}'|xargs -I {} grep {} At_query_blast.txt|awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_A.*/{print $1"\n"$2}'|awk -F "_" '$1~/Ghir/{print $3"_"$4,$5,$6,$8,$0}$1!~/Ghir/{print $2,$3,$4,$6,$0}' OFS="\t" |sort |uniq >converseBed/5_2.bed ``` ## 求只在单个基因组中的剪切事件 在只某一个基因组中保守的类型，先统计与其他基因组中保守的类型的bed1文件，然后剔除掉这些就是只在某个基因组中单独存在的剪切事件了。 ```bash ## At和其他基因组中共有的类型 1_1 1_3 2_1 3_2 4_1 4_2 5_2 cat 1_1.bed 1_3.bed 2_1.bed 3_2.bed 4_1.bed 4_2.bed 5_2.bed |awk -F "\t" '$5~/Ghir_A/{print ">"$5}'|cat - ../At_intron_junction.fasta|grep ">"|sort|uniq -u|sed 's/>//g' |awk -F "_" '{print $3"_"$4,$5,$6,$8,$0}' OFS="\t" >At.bed ## Dt和其他基因组中共有的类型 cat 1_1.bed 1_4.bed 2_2.bed 3_2.bed 4_1.bed 4_2.bed 5_1.bed|awk -F "\t" '$5~/Ghir_D/{print ">"$5}' |cat - ../Dt_intron_junction.fasta|grep ">"|sort|uniq -u|sed 's/>//g' |awk -F "_" '{print $3"_"$4,$5,$6,$8,$0}' OFS="\t" >Dt.bed ## A2和其他基因组中共有的类型 cat 1_1.bed 1_3.bed 2_1.bed 2_2.bed 2_3.bed 4_1.bed 5_1.bed|awk -F "\t" '$5~/evm/{print ">"$5}'|cat - ../A2_intron_junction.fasta|grep ">"|sort|uniq -u|sed 's/>//g' |awk -F "_" '{print $2,$3,$4,$6,$0}' OFS="\t" >A2.bed ## D5和其他基因组共有的类型 cat 1_1.bed 1_4.bed 2_1.bed 2_2.bed 2_3.bed 4_2.bed 5_2.bed|awk -F "\t" '$5~/Gor/{print ">"$5}'|cat - ../D5_intron_junction.fasta|grep ">"|sort|uniq -u|sed 's/>//g' |awk -F "_" '{print $2,$3,$4,$6,$0}' OFS="\t" >D5.bed ``` ## A2与At 、D5与Dt间保守的剪切事件 ```bash ## A2与At中都存在的保守事件 awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_A.*/{print $0}' At_query_blast.txt |awk '$2~/evm/{a[$1][0]+=1}$2~/Ghir_D/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][0]>=1)print i}}' >converse/A2_At_converse.txt ## D5与Dt中都存在的保守事件 awk '$7>200&&$8>90&&$9>85&&$2!~/Ghir_D.*/{print $0}' Dt_query_blast.txt |awk '$2~/Ghir_A/{a[$1][0]+=1}$2~/evm/{a[$1][2]+=1}$2~/Gor/{a[$1][1]+=1}END{for(i in a){if(a[i][1]>=1)print i}}' >converse/D5_Dt_converse.txt ``` ## 统计基因数与事件数 总共有16类，对每一类都统计一下对应的事件的数目 ```bash ##第一大类 wc -l 1_1 awk -F "_" '$1~/Ghir/{print $1"_"$2}$1!~/Ghir/{print $1}' 1_1 |sort |uniq |wc -l wc -l 1_2 wc -l 1_3 awk -F "_" '$1~/Ghir/{print $1"_"$2}$1!~/Ghir/{print $1}' 1_3 |sort |uniq |wc -l wc -l 1_4 awk -F "_" '$1~/Ghir/{print $1"_"$2}$1!~/Ghir/{print $1}' 1_4 |sort |uniq |wc -l ##第二大类 wc -l 2_1 awk -F "_" '$1~/Ghir/{print $1"_"$2}$1!~/Ghir/{print $1}' 2_1 |sort |uniq |wc -l wc -l 2_2 awk -F "_" '$1~/Ghir/{print $1"_"$2}$1!~/Ghir/{print $1}' 2_2 |sort |uniq |wc -l wc -l 2_3 awk -F "_" '$1~/Ghir/{print $1"_"$2}$1!~/Ghir/{print $1}' 2_3|sort |uniq |wc -l ##第三大类 wc -l ../converseBed/At.bed cut -f5 ../converseBed/At.bed|awk -F "_" '$1~/Ghir/{print $1"_"$2}$1!~/Ghir/{print $1}' |sort |uniq |wc -l wc -l 3_2 awk -F "_" '$1~/Ghir/{print $1"_"$2}$1!~/Ghir/{print $1}' 3_2 |sort |uniq |wc -l wc -l ../converseBed/Dt.bed cut -f5 ../converseBed/Dt.bed |awk -F "_" '$1~/Ghir/{print $1"_"$2}$1!~/Ghir/{print $1}'|sort |uniq |wc -l ## 第四大类 wc -l 4_1 awk -F "_" '$1~/Ghir/{print $1"_"$2}$1!~/Ghir/{print $1}' 4_1 |sort |uniq |wc -l wc -l 4_2 awk -F "_" '$1~/Ghir/{print $1"_"$2}$1!~/Ghir/{print $1}' 4_2 |sort |uniq |wc -l wc -l ../converseBed/A2.bed cut -f5 ../converseBed/A2.bed |awk -F "_" '$1~/Ghir/{print $1"_"$2}$1!~/Ghir/{print $1}' |sort |uniq |wc -l wc -l ../converseBed/D5.bed cut -f5 ../converseBed/D5.bed |awk -F "_" '$1~/Ghir/{print $1"_"$2}$1!~/Ghir/{print $1}' |sort |uniq |wc -l ``` ## intronR | 事件 | 基因数目 | 事件数目 | | ---- | ---- | ----- | | 1\_1 | 840 | 1489 | | 1\_2 | 8071 | | | 1\_3 | 1122 | 1869 | | 1\_4 | 1438 | 2247 | | 2\_1 | 651 | 977 | | 2\_2 | 600 | 861 | | 2\_3 | 1043 | 1541 | | At | 3858 | 7683 | | 3\_2 | 1038 | 1571 | | Dt | 3924 | 8053 | | 4\_1 | 579 | 876 | | 4\_2 | 825 | 1245 | | A2 | 4087 | 8451 | | D5 | 5066 | 10217 | ## ExonS | 事件 | 基因数目 | 事件数目 | | ---- | ----- | ---- | | 1\_1 | 38 | 42 | | 1\_2 | 15390 | | | 1\_3 | 103 | 108 | | 1\_4 | 235 | 254 | | 2\_1 | 34 | 35 | | 2\_2 | 38 | 39 | | 2\_3 | 83 | 87 | | At | 712 | 800 | | 3\_2 | 66 | 68 | | Dt | 661 | 740 | | 4\_1 | 22 | 23 | | 4\_2 | 86 | 88 | | A2 | 656 | 759 | | D5 | 1406 | 1648 | ## AltA | 事件 | 基因数目 | 事件数目 | | ---- | ----- | ---- | | 1\_1 | 137 | 167 | | 1\_2 | 12517 | | | 1\_3 | 212 | 262 | | 1\_4 | 423 | 477 | | 2\_1 | 90 | 110 | | 2\_2 | 62 | 69 | | 2\_3 | 172 | 198 | | At | 1805 | 2312 | | 3\_2 | 367 | 421 | | Dt | 1777 | 2359 | | 4\_1 | 81 | 91 | | 4\_2 | 218 | 261 | | A2 | 1555 | 2126 | | D5 | 2127 | 2809 | ## AltD | 事件 | 基因数目 | 事件数目 | | ---- | ----- | ---- | | 1\_1 | 82 | 95 | | 1\_2 | 13275 | | | 1\_3 | 169 | 216 | | 1\_4 | 284 | 331 | | 2\_1 | 66 | 85 | | 2\_2 | 44 | 50 | | 2\_3 | 125 | 146 | | At | 1671 | 2122 | | 3\_2 | 316 | 364 | | Dt | 1714 | 2243 | | 4\_1 | 56 | 62 | | 4\_2 | 148 | 181 | | A2 | 1291 | 1688 | | D5 | 1550 | 1934 | ## 对同源基因间的剪切取交集 Dt中的可变剪切 ```bash 1_1 1_4 2_2 3_2 3_3 4_1 4_2 5_1 ``` At中的可变剪切 ```bash 1_1 1_3 2_1 3_1 3_2 4_1 4_2 5_2 ``` D5中的可变剪切 ```bash 1_1 1_4 2_1 2_2 2_3 4_2 4_4 5_2 ``` A2中的可变剪切 ```bash 1_1 1_3 2_1 2_2 2_3 4_1 4_3 5_1 ``` --- # Agent Instructions: Querying This Documentation If you need additional information that is not directly available in this page, you can query the documentation dynamically by asking a question. Perform an HTTP GET request on the current page URL with the `ask` query parameter: ``` GET https://zpliu.gitbook.io/booknote/ke-bian-jian-qie/shu-ju-chu-li-liu-cheng/20191230-dui-as-lei-xing-jin-hang-ding-yi.md?ask= ``` The question should be specific, self-contained, and written in natural language. The response will contain a direct answer to the question and relevant excerpts and sources from the documentation. Use this mechanism when the answer is not explicitly present in the current page, you need clarification or additional context, or you want to retrieve related documentation sections.