WGLab / PhenoSV

PhenoSV: Interpretable phenotype-aware model for the prioritization of genes affected by structural variants.
MIT License
16 stars 3 forks source link

Insertion marked as "Exonic" but start position is in an Intron #2

Closed MariosGvr closed 8 months ago

MariosGvr commented 8 months ago

Hi, thank you for the incredible tool!

While reviewing my results table, I noticed that certain insertions marked as exonic actually start in an intronic position. Considering both the start and end positions, these insertions potentially "cover" an exonic region.

However, in my understanding, an insertion that begins in an intron should be classified as intronic, as it "pushes" the exon further away rather than spanning it. Is there a specific notion behind which PhenoSV marks those insertions as exonic?

I'm using the terminal version of the PhenoSV tool for annotating multiple SVs, and the web version gives me the same results.

Example:

chrom   start       end     length  svtype      Elements    Type
chr14   39064016    39067194    3178    insertion   SEC23A      Exonic

Start position on an intronic region: image

Range of SV (Start position marked with the blue line, SV spans on the right): image

Thank you!

kaichop commented 8 months ago

If it is just an insertion at a location, then the start and end needs to be identical to each other. Here you have a 3.2kb genomic region in the file, and this is not an insertion. Clearly this region covers exons.

On Tue, Jan 23, 2024 at 2:03 PM MariosGvr @.***> wrote:

Hi, thank you for the incredible tool!

While reviewing my results table, I noticed that certain insertions marked as exonic actually start in an intronic position. Considering both the start and end positions, these insertions potentially "cover" an exonic region.

However, in my understanding, an insertion that begins in an intron should be classified as intronic, as it "pushes" the exon further away rather than spanning it. Is there a specific notion behind which PhenoSV marks those insertions as exonic?

I'm using the terminal version of the PhenoSV tool for annotating multiple SVs, and the web version gives me the same results.

Example:

chrom start end length svtype Elements Type chr14 39064016 39067194 3178 insertion SEC23A Exonic

Start position on an intronic region: image.png (view on web) https://github.com/WGLab/PhenoSV/assets/101292511/30fc9cb7-2910-45d7-a6ae-e64ab66512ef

Range of SV (Start position marked with the blue line, SV spans on the right): image.png (view on web) https://github.com/WGLab/PhenoSV/assets/101292511/8bd95878-0f32-467e-a78b-bda93ffe2037

Thank you!

— Reply to this email directly, view it on GitHub https://github.com/WGLab/PhenoSV/issues/2, or unsubscribe https://github.com/notifications/unsubscribe-auth/ABNG3OGZAA6ASB4BA4W5KYDYQACRZAVCNFSM6AAAAABCHQS7EGVHI2DSMVQWIX3LMV43ASLTON2WKOZSGA4TMNZZGEZTGOI . You are receiving this because you are subscribed to this thread.Message ID: @.***>

MariosGvr commented 8 months ago

Thank you very much for your input! I was able to work around the issue.