This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Why does the transcriptome from squid sp U. aligns only 1% with squid sp. E genome?

I downloaded the transcriptome from a published paper where they sequenced the transcriptome from different organs of 2 different species of squid (squid U and squid E). Our lab recently sequenced the whole genome of Squid E. When I aligned the transcriptome of squid E from the previously published paper to its own genome, it aligns 80% (I used bowtie2 and hisat). However, the transcriptome of squid U processed the exact same way as squid E aligns only 1% or even less to the genome. The RNA was isolated and sequenced the exact same way in both the species and I processed them both the exact same way. Does anyone have theories why this might be the case?

This is the paper from which I used their transcriptomes: https://www.pnas.org/content/111/44/E4736.short

rna-seq
• 911 views • link
written by nidhi.vijayan13 • nidhi.vijayan13

When I aligned the transcriptome

Did you align the sequencing reads, or the asembled transcriptome?

(I am assuming you tried to map the assembled transcriptome to a reference genome.) If squid U isn't very similar to squid E, bowtie2 and hisat aren't good tools for the job. You may try minimap2, gmap or spaln, they will do a better job for somewhat divergent species.

• 0 views • link
written by h.mon h.mon

0 answers

No answers yet.

Log in to answer this question.

More posts like this