Question

Which aligment tool is best for eukaryotic transcriptome?

1

Entering edit mode

6.3 years ago

sudha ▴ 20

Dear one, Which alignment tool is best to align eukaryotic (homo sapiens) transcriptome? Please comment your answer with reason why! Thanks in advance!

p.s. My sequences are paired end & human samples! My intention is to do differential gene expression study!

p.s. I have seen that there are many alignment tools are used and given in publications, so confused a lot! I am new to bioinformatics!

RNA-Seq alignment sequencing • 2.0k views

ADD COMMENT • link updated 6.3 years ago by GenoMax 141k • written 6.3 years ago by sudha ▴ 20

3

Entering edit mode

STAR is a good starter.

ADD REPLY • link 6.3 years ago by arta ▴ 670

score 2 · Answer 1 · 2018-01-16

2

Entering edit mode

6.3 years ago

Matteo Schiavinato ★ 3.6k

Alignment of RNASeq reads on a genomic reference

STAR (link) -> super-fast, gold standard for many people, awesome 2-pass mode (see manual)
HISAT2 (link) -> same as Tophat2, much, much faster, good precision
BWA (link) -> handles well also supplementary alignments (2048 bitwise flag) to then retrieve chimeras

all handle spliced alignment

Alignment of RNASeq reads on a transcriptomic reference

RSEM (link) -> accurate quantification of gene and isoform expression (probably your best choice, imho, for this)

ADD COMMENT • link 6.3 years ago by Matteo Schiavinato ★ 3.6k

2

Entering edit mode

personally don't agree with the "all handle spliced alignment" statement , they all can handle split or gapped alignment but that's a totally different thing than spliced alignment (== splice-site aware aligner) !

ADD REPLY • link 6.3 years ago by lieven.sterck 15k

1

Entering edit mode

You are right, that is what I meant, my bad.

ADD REPLY • link 6.3 years ago by Matteo Schiavinato ★ 3.6k

1

Entering edit mode

Thanks a lot for your answers, it made my day!

but (for alignment) RSEM uses bowtie2 only right ?

ADD REPLY • link 6.3 years ago by sudha ▴ 20

1

Entering edit mode

Yeah, or bowtie 1, if you specify it. But it works quite fine, imho. However, you can use HISAT2 on transcriptome references tailoring the output directly for cufflinks (--dta-cufflinks) if this helps you.

ADD REPLY • link 6.3 years ago by Matteo Schiavinato ★ 3.6k

score 2 · Answer 2 · 2018-01-16

2

Entering edit mode

6.3 years ago

Carlo Yague 8.7k

You should also consider pseudo-alignments tools, such as kallisto and salmon. They are extremely fast and lightweight, can quantify millions of reads under a few minutes on a laptop computer.

ADD COMMENT • link 6.3 years ago by Carlo Yague 8.7k

score 2 · Answer 3 · 2018-01-16

BBMap suite (includes scanning/trimming program bbduk.sh and many things other than the aligner, which is bbmap.sh). I will not enumerate its advantages, since they should be apparent, if/when you try it.

Many aligners will need ~ 30G of RAM for human genome. If you don't have that available then bwa may be the one to try. It minimally needs ~6 G RAM.