determination of alternative splicing events
0
0
Entering edit mode
6.5 years ago

hi all,

for my project I need to determine and classify the alternate splicing (AS) events. I have used ASTALAVASTA software for this purpose that has successfully determined the AS events with code for the gtf annotation input file generated by stringtie. the problem is that it is not providing me any information about the classification of these events like whetehr they are exon skipping, intron retension, alternative doner or alterntive acceptors. it provides codes only for these 5 classes but my results contain large variety of other types (mixed, have transcription start/end or novel for which it provide no info). I have tried to group the one that have mixed anatomy in these 5 classes but it is pretty confusing. can anyone help me how could i classify these events???

or may suggest any other good tool for this purpose plz.

kind help would be highly appreciated.

my output file of ASTALAVASTA looks like this:

  chr1  Undefined   as_event    581052  581358  .   +   .   transcript_id "STRG.5.2,STRG.5.1"; gene_id "chr1:581052-583794W"; flanks "581052[,581358-"; structure "1^,2^"; splice_chain "581257^,581276^"; sources "Undefined,Undefined"; NMorMRNA "null"; degree "4"; dimension "2_2"; 
chr1    Undefined   as_event    581052  582101  .   +   .   transcript_id "STRG.5.3,STRG.5.2"; gene_id "chr1:581052-583794W"; flanks "581052[,582101]"; structure "0,1^2-"; splice_chain ",581257^581358-"; sources "Undefined,Undefined"; NMorMRNA "null"; degree "4"; dimension "2_2"; 
chr1    Undefined   as_event    581052  583794  .   +   .   transcript_id "STRG.5.1,STRG.5.3"; gene_id "chr1:581052-583794W"; flanks "581052[,null"; structure "1^2-4^5-6],3]"; splice_chain "581276^581358-582109^583479-583794],582101]"; sources "Undefined,Undefined"; NMorMRNA "null"; degree "4"; dimension "2_3"; 
chr1    Undefined   as_event    950006  957972  .   +   .   transcript_id "STRG.34.1,STRG.35.1"; gene_id "chr1:950006-957972W"; flanks "null,null"; structure "1[2^7-8],3[4^5-6]"; splice_chain "950006[950090^957762-957972],950774[950968^951662-952237]"; sources "Undefined,Undefined"; NMorMRNA "null"; degree "4"; dimension "2_2"; 
chr1    Undefined   as_event    974385  976220  .   +   .   transcript_id "STRG.40.1,STRG.40.2"; gene_id "chr1:973650-976220W"; flanks "974385[,976220]"; structure "1^2-4^6-,3^5-"; splice_chain "975147^975385-975564^975721-,975477^975634-"; sources "Undefined,Undefined"; NMorMRNA "null"; degree "4"; dimension "2_2"; 
chr1    Undefined   as_event    973650  976220  .   +   .   transcript_id "STRG.40.3,STRG.40.1"; gene_id "chr1:973650-976220W"; flanks "null,null"; structure "1[7],2[3^4-5^6-8]"; splice_chain "973650[975869],974385[975147^975385-975564^975721-976220]"; sources "Undefined,Undefined"; NMorMRNA "null"; degree "4"; dimension "2_3"; 
chr1    Undefined   as_event    973650  976220  .   +   .   transcript_id "STRG.40.3,STRG.40.2"; gene_id "chr1:973650-976220W"; flanks "null,null"; structure "1[5],2[3^4-6]"; splice_chain "973650[975869],974385[975477^975634-976220]"; sources "Undefined,Undefined"; NMorMRNA "null"; degree "4"; dimension "2_3"; 
chr1    Undefined   as_event    1297130 1298623 .   +   .   transcript_id "STRG.53.1,STRG.53.2"; gene_id "chr1:1297130-1298623W"; flanks "1297130[,null"; structure "1^3-4],2]"; splice_chain "1298245^1298432-1298623],1298254]"; sources "Undefined,Undefined"; NMorMRNA "null"; degree "4"; dimension "2_2"; 
chr1    Undefined   as_event    1646433 1647884 .   +   .   transcript_id "STRG.70.1,STRG.70.2"; gene_id "chr1:1646433-1647884W"; flanks "null,1647884]"; structure "1[2^4-,3["; splice_chain "1646433[1646513^1646837-,1646754["; sources "Undefined,Undefined"; NMorMRNA "null"; degree "4"; dimension "2_2"; 
chr1    Undefined   as_event    1869897 1870747 .   +   .   transcript_id "STRG.79.1/STRG.79.3,STRG.79.2"; gene_id "chr1:1868804-1873243W"; flanks "1869897^,1870747^"; structure "1-,2-"; splice_chain "1870664-,1870669-"; sources "Undefined,Undefined"; NMorMRNA "null"; degree "4"; dimension "2_2"; 
chr1    Undefined   as_event    1871874 1873243 .   +   .   transcript_id "STRG.79.1/STRG.79.2,STRG.79.3"; gene_id "chr1:1868804-1873243W"; flanks "1871874-,null"; structure "1^3-4],2]"; splice_chain "1872150^1873095-1873243],1872197]"; sources "Undefined,Undefined"; NMorMRNA "null"; degree "4"; dimension "2_2"; 
chr1    Undefined   as_event    1886878 1887691 .   +   .   transcript_id "STRG.81.1,STRG.81.2"; gene_id "chr1:1884860-1891139W"; flanks "1886878^,1887691^"; structure "1-,2-"; splice_chain "1887614-,1887619-"; sources "Undefined,Undefined"; NMorMRNA "null"; degree "4"; dimension "2_2"; 
chr1    Undefined   as_event    2119498 2121140 .   +   .   transcript_id "STRG.98.1,STRG.98.2"; gene_id "chr1:2118059-2121140W"; flanks "2119498-,null"; structure "1^3-4^5-6],2]"; splice_chain "2119561^2119968-2120288^2120951-2121140],2119570]"; sources "Undefined,Undefined"; NMorMRNA "null"; degree "4"; dimension "2_2"; 
chr1    Undefined   as_event    2381238 2381766 .   +   .   transcript_id "STRG.120.2,STRG.120.1"; gene_id "chr1:2380953-2383475W"; flanks "2381238^,2381766-"; structure "0,1-2^"; splice_chain ",2381490-2381595^"; sources "Undefined,Undefined"; NMorMRNA "null"; degree "4"; dimension "2_2";
next-gen • 1.0k views
ADD COMMENT

Login before adding your answer.

Traffic: 2197 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