This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to cluster transcripts assembled using TRINITY into unigene after blast?

I have assembled a transcriptome using TRINITY with 100,000 transcripts. I have perform blastx against NR database (-outfmt 5). I want to remove redundancy in the assembled transcript before proceed to further processing.

1) How can I cluster them together into unigenes and remove redundancy?

2) Or how to select transcipts with longest sequences if they return same hit?

3) Any software or program suggestion for doing this?

4) Is it necessary for removing redundancy?

Thank you very much.

rna-seq assembly gene
• 7,532 views • link
updated by lakhujanivijay lakhujanivijay • written by sbchua.1990 • sbchua.1990

1 answer

You will get transcripts after assembling short reads. Depending on alternative splicing,one gene can have more than one transcript.

Clustering or grouping puts all these similar sequences together and help to make a set of transcript for one gene.

On the other hand assembly process can create some transcripts that are not real (sequences/transcripts with more than 95% identity to another sequence) and clustering helps identifying them. Therefore, the flow goes like:

Assembly -> Transcripts -> clustering transcripts -> unigenes -> further downstream processing.

You can use CD-HIT for this purpose.

• 0 views • link
written by lakhujanivijay lakhujanivijay

Log in to answer this question.

More posts like this