What is your preferred way to map cDNA sequences to genome?
3
1
Entering edit mode
6.4 years ago
jon.brate ▴ 250

I have a set of cDNA-sequences representing mRNAs and I want to create a gtf/gff file of these. I have tried to use BLAST and BLAT, but I find it difficult to obtain a single reliabe hit that represents all the different exons. And especially with BLAT I get many small hits for each sequence. I also tried to map the sequences using TopHat, but for some reason it crashed, perhaps due to the long sequences.

Any suggestions on how to map such sequences to a genome with reliable information about the exon/intron junctions?

mapping alignment blat blast gff • 2.6k views
ADD COMMENT
7
Entering edit mode
6.4 years ago

Use gmap. It was build for that and outputs a gff file with each exon of every sequence in input. Very usefull !

ADD COMMENT
0
Entering edit mode

Thanks, I'll check it out. Does it also create an entry representing the entire gene in the gff-file?

ADD REPLY
0
Entering edit mode

Aligners do not really have the concept of a gene, only transcript.  If you want to treat each transcript as a gene, that can be done pretty easily by adding a tag to the resulting GFF annotation column in a post-alignment step.

ADD REPLY
2
Entering edit mode
6.4 years ago

I never had to solve this problem myself but I am guessing that a long read aligner  such as lastz executed with smaller (zero?) penalty for gap extension might work.

 

ADD COMMENT
0
Entering edit mode

Thanks! I never heard of that program. Will check it out.

ADD REPLY
0
Entering edit mode
6.4 years ago
Ram 34k

Have you tried blastx to compare against protein databases and then backtrack to the genes from the top protein hit? That could solve your problem with minimal human intervention, if you script the efetch to get gene from protein identifier. Exon junctions can then be an SQL query to UCSC Genome Browser.

ADD COMMENT
0
Entering edit mode

I didn't mention that, but in my case I have non-protein coding sequences from a non-model species.

ADD REPLY
0
Entering edit mode

Oh, I see. My solution is not useful then. Looks like there are better tools anyway!

ADD REPLY

Login before adding your answer.

Traffic: 2358 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6