This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Alternate Transcripts Paralogs

Hello,

I have a few doubts related to how the best way is to go about finding these:

  1. Finding unique regions when you compare various transcripts of a gene in one organism.
  2. Finding evidences for a particular transcript.
  3. Finding the % identity IN NCBI PROTEIN database for a particular sequence taken from uniprot. i.e i blast a protein sequence in uniprot database and get the paralogs for that along with percent identity. then i want to find the corresponding sequences for those paralogs in NCBI protein and find the %ID for my query sequence against those NCBI protein sequences. i tried ID mapping in uniprot and it gives me the corresponding IDs. but when i blast my query sequence against these IDs the %identity is way off from the %ID uniprot gives me.

any help/suggestions would be greatly appreciated.

gene uniprot ncbi
• 4,391 views • link
written by aj123 aj123

I am a bit confused. Do you mean you want to differentiate between alternative isoforms and paralogs?

• 0 views • link
written by Damian Kao Damian Kao

No, i just want to find orthologs and paralogs for a particular gene from a particular species (human). Apart from this, I also want to find alternate transcripts for a gene in human.

• 0 views • link
written by aj123 aj123

2 answers

  1. I assume that you want to identify isoform-specific exons for a gene, in that case, you need to overlap your annotations, and select the regions transcript-specific. You can do that with R:bioconductor and IRanges.
  2. Evidence of what? expression? protein? conservation?
  3. Are you using the same blast parameters? if the parameters are the same, the alignments must be similar.
• 1 views • link
written by JC JC

@JC Thanks so much for your reply.

1) is there a way to do this using just online bioinfo databases like NCBI, Ensembl etc with out using IRanges? also, how would i a) find the annotations for a transcript b)go aobut overlapping my annotations c)select the regions which are transcript specific? i mean how would i know a region is transcript specific?

2) evidence of expression

3)yes, all blast parameters used are standard ones. im not sure still how/why this is happening

• 0 views • link
written by aj123 aj123

1) and 2) depends on your model, human? mouse? 3) did you check if the alignments are the same? I'm thinking if you have 2 sequences, A with 100aa and B with 200aa, and you found a common alignment of 50aa, when you compare A->B the %ID is be 50% (based on A) but B->A is 25% (based on B)

• 0 views • link
written by JC JC

Im using Human. Im not sure how you mean w.r.t the second part. im blasting always in a standard format. i.e original sequence vs. subject sequences. Still a little confused w.r.t finding the alternate transcripts for a gene.

• 0 views • link
written by aj123 aj123

Default (standard) blast parameters differ between NCBI and UniProt blast sites. You have to be very careful about each setting.

• 0 views • link
written by Jerven • Jerven

@JC Thanks so much for your reply.

1) is there a way to do this using just online bioinfo databases like NCBI, Ensembl etc with out using IRanges? also, how would i a) find the annotations for a transcript b)go aobut overlapping my annotations c)select the regions which are transcript specific? i mean how would i know a region is transcript specific?

2) evidence of expression

3)yes, all blast parameters used are standard ones. im not sure still how/why this is happening

• 0 views • link
written by aj123 aj123

please add this as a comment to the answer then delete it - right now it is listed as an answer, but it is not really that

• 0 views • link
written by Istvan Albert Istvan Albert

oops! didnt know how to add it as a reply to a comment. have done it now. :)

• 0 views • link
written by aj123 aj123

Log in to answer this question.

More posts like this