This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Transcript assemblers from RNA-seq reads - a request for tool recommendation

Hi all,

I would like to examine the complete transcripts sequence of specific genes on the RNA-seq. As input, I can use the RNA-seq FASTQ file. for example: Looking at all BRAF transcripts full transcript (all cDNA combinations) on the FASTQ file (could be pair end or not)

What will be the recommended way to do so? It will be great if I could avoid processing the whole sequence at some point of the processing pipeline and focus on the requested gene.

Many thanks, Eila

rna-seq reads assembler
• 1,936 views • link
written by eilalandkof • eilalandkof

2 answers

As far as I know, the only way to know which reads come from the gene BRAF is to align/pseudo-align your sequences. Of course you can do it on only the BRAF genome sequence, but I would advise against it. Why ? Because if a read comes from a transcript with local homology with BRAF, it might get aligned to BRAF if it is the only hit possible. In contrast, if you align your reads on the full genome/transcriptome, it will be aligned to BRAF only if this provides the best alignment compared to other possible genes. This kind of healthy competition between alignment is important to avoid spurious results.

• 0 views • link
written by Carlo Yague Carlo Yague

Hi Carlo, thank you for the quick answer. As far as I understand pseudo alignment is being used for genes / transcript expression levels calculation and does not “need” to read the full transcripts matching to do so. My experience with this kind of tools, such as Kallisto, and it does not report the full mRNA transcript sequence.

Are you familiar with and alignment tool that reports the full mRNA / cDNA transcript for the aligned transcript? I have added an image for what I am looking for.

Thanks, Eila enter image description here

• 0 views • link
written by eilalandkof • eilalandkof

You could try to align your reads to the genome (using hisat2 or start) and then use a genome-guided assembler to retrieve the mRNA sequences (for instance stringtie2). I would suggest aligning against the complete genome in the first step. You could try and filter the reads by position after that before the assembly, using samtools, to improve speed (and also only keep lightweight bam files for future use). If there is no genome avaiable I recomend assembling the whole transcriptome using trinity as refered above.

• 0 views • link
written by Rogerio Ribeiro • Rogerio Ribeiro

Along with the important caveats from Carlo Yague, you can check for BRAF transcripts de novo by assembling the whole transcriptome with Trinity https://github.com/trinityrnaseq/trinityrnaseq. I think this is available on Galaxy as well if you don't have local compute with big RAM (>64-128 GB likely necessary depending on your data).

• 0 views • link
written by colindaven colindaven

Thanks for the quick responses. I will try these recommendations. I have cloud computing available and using the human hg38 reference genome. So, should be fine to go with the non de-novo method

Eila

• 0 views • link
written by eilalandkof • eilalandkof

Log in to answer this question.

More posts like this