exon bed file

[ArashDepp]

Hi,
I have two questions regarding the exon bed file. Would appreciate your time and help in this.

1. Should I provide non-overlapping disjoint exons for a gene OR all the exons from gencode gtf file would work (which often overlap at their boundaries)?

2. In case I provide only the non-overlapping intervals for a gene, will free PSI infer other possible isoforms having exons with different ends than the provided bed files?

Best.

[JY-Zhou]

Sorry for my delayed reply, and appreciate your interest in FreePSI. :)

For the first question, FreePSI can only deal with non-overlapping disjoint exons for a gene. If there are overlapped exons, I suggest preprocessing the annotation files into non-overlapping exons.
This preprocessing can always work, e.g. if A and B are overlapped, it can be represented as A, A∩B, B.
This script can help you to do the preprocess: scripts/mergeToFlatGene.py.

For the second question, the answer is no. FreePSI only estimates PSI value for the provided non-overlapping exons. However, if you have some potential regions that might be new exons, you can directly add these regions to the annotation file, and check the estimated PSI values for these regions. (Note that FreePSI doesn't care about the composition of isoforms.)

Hope this would help, and really sorry for my late reply.

[ArashDepp]

Hi. Got it. Thanks for helping me with this. I will post in case I have any more concerns. Best regards.

…

On Sun, Mar 13, 2022 at 10:07 AM Arrows ***@***.***> wrote: Sorry for my delayed reply, and appreciate your interest in FreePSI. :) For the first question, FreePSI can only deal with non-overlapping disjoint exons for a gene. If there are overlapped exons, I suggest preprocessing the annotation files into non-overlapping exons. This preprocessing can always work, e.g. if A and B are overlapped, it can be represented as A, A∩B, B. This script can help you to do the preprocess: scripts/mergeToFlatGene.py <https://github.com/JY-Zhou/FreePSI/blob/82830c5df94d33fe9e0df4054905e4c3bb47b381/scripts/mergeToFlatGene.py> . For the second question, the answer is no. FreePSI only estimates PSI value for the provided non-overlapping exons. However, if you have some potential regions that might be new exons, you can directly add these regions to the annotation file, and check the estimated PSI values for these regions. (Note that FreePSI doesn't care about the composition of isoforms.) Hope this would help, and really sorry for my late reply. — Reply to this email directly, view it on GitHub <#7 (comment)>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AKI6R77NRFC3YGGGI35IGITU7XZBNANCNFSM5QHAYNXA> . Triage notifications on the go with GitHub Mobile for iOS <https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675> or Android <https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub>. You are receiving this because you authored the thread.Message ID: ***@***.***>