Question: Without genome reference, how to find exon-exon junctions given isoforms
0
gravatar for lwc628
4.7 years ago by
lwc628200
United States
lwc628200 wrote:

Is there any program that detects the exon-exon junctures when provided with isoforms without using the genome reference sequence?

 

rna-seq • 1.5k views
ADD COMMENTlink modified 4.7 years ago by spiderdijon20 • written 4.7 years ago by lwc628200

I'm not sure exactly what you're trying to do. You want to detect non-canonical exon-exon junctions correct? What do you mean when provided with isoforms - you want to use the isoform's exons as a basis for the spice-junctions? Why do you need to do this without a reference?
 

ADD REPLYlink written 4.7 years ago by spiderdijon20

Thanks for your comments. I generated a denovo transcriptome assembly using Trinity, and in it I have multiple isoforms for a gene model. Without reference, I don't know exon-exon junctions in those isoforms, which is why I asked this question.

ADD REPLYlink written 4.7 years ago by lwc628200

What organism are you working with?

ADD REPLYlink written 4.7 years ago by spiderdijon20
1
gravatar for spiderdijon
4.7 years ago by
spiderdijon20
Australia
spiderdijon20 wrote:

I'm going to assume your organism has a reference genome, otherwise this is a different question all together. From a scan of the Trinity documentation, it looks like it gives you the assembled fasta sequences of your genes and transcripts. So to find the exon-exon junctions, you're going to need to compare your transcript sequences to your reference organism's sequence. Try using something like BLAT to align your transcripts (https://genome.ucsc.edu/cgi-bin/hgBlat), which will give you the coordinates that each "exon" maps to. Depending on how many transcripts you need to do this for, you may need to run the command-line version of BLAT, then convert the output to something like BED format, assuming you need to visualise your transcript exon boundaries.

ADD COMMENTlink written 4.7 years ago by spiderdijon20
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1548 users visited in the last hour