I am trying to compute PSI values for es splicing events detected by diffSplice and have noticed that some splicing events belong to more than one gene. For example:
This is the list of isoforms that contribute to a particular es event:
m54284U_211218235859/153224706/ccsENSG00000156875.13,m54284U_211226025210/100141500/ccsENSG00000283761.1,m54284U_211116_013100/148899678/ccs_ENSG00000283761.1,m54284U_211226_025210/166922717/ccs_ENSG00000156875.13,m54284U_211203_202148/119604216/ccs_ENSG00000156875.13,m54284U_211218_235859/121045159/ccs_ENSG00000156875.13,m54284U_211218_235859/168560587/ccs_ENSG00000156875.13,m54284U_211226_025210/133956376/ccs_ENSG00000156875.13,m54284U_200726_233028/93325783/ccs_ENSG00000283761.1,m54284U_200720_151958/73663327/ccs_ENSG00000156875.13,m54284U_211116_013100/106430925/ccs_ENSG00000156875.13,ENST00000370152.7_ENSG00000156875.13,ENST00000640357.1_ENSG00000283761.1,m54284U_211228_151634/81396504/ccs_ENSG00000156875.13,m54284U_200726_233028/166724650/ccs_ENSG00000283761.1,m54284U_211218_235859/118948776/ccs_ENSG00000283761.1,ENST00000639037.1_ENSG00000283761.1
As I've highlighted in bold, there are isoforms from 2 different genes that contribute to this particular event. This behavior seemed strange and counter-intuitive to me, as I would expect each splicing event to belong to a single gene.
Do you have any insight on whether this is a bug or something else? Do you have any ideas for if / when / how often we should see this kind of behavior? I think it would greatly affect the interpretations being made from the data. For example, let's say we compute the PSI value for this particular event. It will tell us how often the exon is included, but it's an aggregate across two different genes and does not tell us the relative inclusion of the event in one gene or the other.
Thank you very much for your help and time, it is much appreciated!
I am trying to compute PSI values for es splicing events detected by diffSplice and have noticed that some splicing events belong to more than one gene. For example:
This is the list of isoforms that contribute to a particular es event: m54284U_211218235859/153224706/ccsENSG00000156875.13,m54284U_211226025210/100141500/ccsENSG00000283761.1,m54284U_211116_013100/148899678/ccs_ENSG00000283761.1,m54284U_211226_025210/166922717/ccs_ENSG00000156875.13,m54284U_211203_202148/119604216/ccs_ENSG00000156875.13,m54284U_211218_235859/121045159/ccs_ENSG00000156875.13,m54284U_211218_235859/168560587/ccs_ENSG00000156875.13,m54284U_211226_025210/133956376/ccs_ENSG00000156875.13,m54284U_200726_233028/93325783/ccs_ENSG00000283761.1,m54284U_200720_151958/73663327/ccs_ENSG00000156875.13,m54284U_211116_013100/106430925/ccs_ENSG00000156875.13,ENST00000370152.7_ENSG00000156875.13,ENST00000640357.1_ENSG00000283761.1,m54284U_211228_151634/81396504/ccs_ENSG00000156875.13,m54284U_200726_233028/166724650/ccs_ENSG00000283761.1,m54284U_211218_235859/118948776/ccs_ENSG00000283761.1,ENST00000639037.1_ENSG00000283761.1
As I've highlighted in bold, there are isoforms from 2 different genes that contribute to this particular event. This behavior seemed strange and counter-intuitive to me, as I would expect each splicing event to belong to a single gene.
Do you have any insight on whether this is a bug or something else? Do you have any ideas for if / when / how often we should see this kind of behavior? I think it would greatly affect the interpretations being made from the data. For example, let's say we compute the PSI value for this particular event. It will tell us how often the exon is included, but it's an aggregate across two different genes and does not tell us the relative inclusion of the event in one gene or the other.
Thank you very much for your help and time, it is much appreciated!