1
$\begingroup$

I would like to ask why in RNA seq analysis (alignment step) we use sometimes reference genome instead of reference transcriptome? thank you!

$\endgroup$
7
$\begingroup$

Pseudoaligmnent to a reference transcriptome with something like Kallisto is perfectly valid, assuming you are in an organism whose reference transcriptome is very well documented.

But if you want to find novel transcripts or novel splicing junctions, you need to align to genome.

$\endgroup$
11
  • $\begingroup$ thank you! Could I use STAR or a similar program (not Kallisto or Salmon etc) with reference transcriptome? $\endgroup$
    – marilu
    Aug 5 '20 at 20:57
  • 2
    $\begingroup$ You can make a STAR reference out of a transcriptome, but why would you do that? It wasn't designed to be run like that. $\endgroup$
    – swbarnes2
    Aug 5 '20 at 21:07
  • $\begingroup$ if i am only interested in differential expression, could i use star with reference transcriptome to decrease computational effort? $\endgroup$
    – marilu
    Aug 5 '20 at 21:24
  • 2
    $\begingroup$ As a rule, do not use programs for things they were not designed for unless you are absolutely sure you understand what you are doing and why. If you want something faster than STAR, use Kallisto. $\endgroup$
    – swbarnes2
    Aug 5 '20 at 21:34
  • $\begingroup$ I actually think you could run STAR on a transcriptome. In their manuscript, you can run alignment using a GTF file. $\endgroup$
    – Phoenix Mu
    Aug 5 '20 at 21:36

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.