The origin of transcripts with putative frameshift in the de novo assembled transcriptome?
0
1
Entering edit mode
10.3 years ago
seta ★ 1.9k

Dear all,

I've performed a denovo transcriptome assembly of a non-model plants on trimmed and almost high quality illumina reads. I used Trapid tool to fast evaluation of assembled transcriptome, this tool reported high percentage of transcript with putative frameshift that is not desired. Could you please let me know why there is many transcript with putative frameshift in the transcriptome based on Trapid prediction, and how I can solve this problem? Thanks for your points and tips.

trapid next-gen Assembly frame-shift RNA-Seq • 2.4k views
ADD COMMENT
0
Entering edit mode

I suggest this might be an artifact, it is just difficult to exactly say why. One source might be the transcriptome assembly, which tool and parameters were used? Whatever you used, try a different assembly in addition. Another option is the predictions are bogus. Take some of those predicted frameshift containing transcripts and BlastX. Does the result show frameshifts as well: multiple HSPs in different frames for the best hits?

ADD REPLY
0
Entering edit mode

Thanks. I used CLC software for assembly of about 300 million PE reads. The high-quality assembly, in terms of N50 (1270bp), contig number and percentage of mapped back read (89%) was obtained in K-mer 64. However, Trapid tool detected many (about 60%) transcript with putative framshifts. I'm really concern about this results. Could you please share me other tools to detect putative framshifts?

ADD REPLY
0
Entering edit mode

I remember there was a similar post where we concluded that there was no evidence for a frame shift from blastx, correct? I think then we can regard this as resolved?

ADD REPLY
0
Entering edit mode

Hi Michael, yeah and you kindly suggested me check the frame shift using blastx and

  • Retain all best hits or all hits with score > some threshold
  • Retain all hits with more than one HSP
  • Retain all hits where at least one HSP has different frame from the others
  • All queries that pass these filters are frame shift candidates

You was supposed to share a script to evaluate the above parameters from blastx, could you please put it here? Thanks in advance

ADD REPLY
0
Entering edit mode
ADD REPLY

Login before adding your answer.

Traffic: 2859 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