This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Aligning Short Reads To Homologous Sequences?

I'm working on a de novo transcriptome project and the genome has not been sequenced for our plant. In general, trinity is working pretty well, but we are seeing some spurious fusion genes (N terminus of gene A fused to C terminus of gene B). For the questions we are asking, having the correct full-length protein sequence is important.

I was wondering if we could get around this problem by first mapping short reads to homologous protein sequences from Arabidopsis and then feeding each of those subcollections into a contig builder? Is there any software that acts like blastx but works well on short Illumina reads?

Thanks

transcriptome short aligner
• 4,031 views • link
updated by rob234king rob234king • written by Goldbear • Goldbear

Hi, length of the reads and the divergence between your species and Arabidopsis will be very useful to know first.

• 0 views • link
written by Haibao Tang Haibao Tang

Hi, my reads are about 90bp. unfortunately the divergence is pretty far. Our strategy may be to fully sequence one representative species from our clade and then use that as a reference for the other species. Thanks!

• 0 views • link
written by Goldbear • Goldbear

2 answers

There is a blastx-based pipeline called STM (a Genome Research paper) which uses de novo contigs and blastx search results to proteins from a closely related species to improve de novo assembly. Maybe it's also not hard to put together a similar pipeline using perl, in which blastx outputs can be dissected more closely. However, I think choosing the right reference taxon should be a very important factor.

• 1 views • link
written by Vitis Vitis

Thanks! This was exactly what I was looking for.

• 0 views • link
written by Goldbear • Goldbear

There is a setting on trinity to attempt to address fusion transcripts, try again but use --jaccard_clip

• 0 views • link
written by rob234king rob234king

Log in to answer this question.

More posts like this