STAR: ultrafast universal RNA-seq aligner

[JonB]

I have problems generating a new genome.
My genome is more than 300.000 contigs and ~300MB in size.

Feb 18 11:25:11 ..... Started STAR run
Feb 18 11:25:11 ... Starting to generate Genome files
/var/spool/slurmd/job1306271/slurm_script: line 9: 12129 Killed

[NGSfan]

That is a great feature, Alex. I am looking forward to it.

I currently convert cufflinks output into read counts for DESeq.

You read in the length and coverage columns from the isoforms.fpkm_tracking file as well as the read length of your reads and then caclulate:

my $reads = $length * $coverage / $readlength;

I have compared the values to htseq-count and found them to be nearly the same.

The advantage over htseq is you can multi-thread cufflinks.

[iris_aurelia]

Hi Alex,

Yes it is a pity that samtools doesn't support multiple threads. We've tried Novosort and that seems to work fine, however you need a license in order to run it in a multithreaded way.

HT-seq is what I usually use to get the readcounts. Are you gonna make this step multithreaded as well? It would indeed be nice if you could just run a sample and get the whole bunch of information all in once.

[alexdobin]

@volks:
Good example from @iris_aurelia.
If you have multiple files you wish to map in one run, they should be separated by commas, while paired-end mates are separated by space:
--readFilesIn Read1a.gz,Read1b.gz Read2a.gz,Read2b.gz

'zcat' should be on your path, or you can use absolute path such as '/bin/zcat'


@iris_aurelia:
Presently I am working on outputting sorted BAM, since 'samtools sort' is a bottleneck for many users (because samtools is not threaded). Then I will code "read counts per gene" (similar to HT-seq counts). Then you will be able to map the reads and get the counts in one run. Hopefully, it will be ready in 1-2 months.

[iris_aurelia]

I've used it as follows, which works fine:

STAR_2.3.0e.Linux_x86_64/STAR --genomeDir /dir/STAR/Genome --readFilesIn /dir/*.fastq.gz --readFilesCommand zcat --outFileNamePrefix SampleA --runThreadN 20

I hope the implementation of the readcounts will follow soon, since I'm curious to compare this tool to the tools I'm using at this moment.

[volks]

i also like STAR a lot. so much simpler and faster than other tools around.

could you give an example how to use gzipped input files (--readFilesCommand)? i cannot get it to work.

thanks!

[NGSfan]

Yes, I am also happy with the results from STAR. It is incredibly fast and easy to use. I will test the shared memory option again now that I have my own workstation machine to play with.

Compatibility with Cufflinks has also been great. Works with HT-seq as well.

Thank you!!

[aldeluis]

Your program is great. It's solving my alignments in a fraction of time respect TopHat. Thank you!!

[alexdobin]

Dear All,

I would like to formally introduce on SEQanswers our RNA-seq mapper STAR.
Its advantages include:

• Very high mapping speed:
  on a modest 12-core cluster STAR maps 400 Million pairs per hour for human 2x100 Illumina reads (>50 times faster than TopHat).
• Accurate alignment of contiguous and spliced reads:
  in our tests on real and simulated data STAR showed better sensitivity and precision than TopHat.
• Detection of polyA-tails, non-canonical splices and chimeric (fusion) junctions.
• Mapping reads of any length:
  STAR can efficiently map reads of any length generated by current or emerging sequencing platforms, starting from ~15 bases (small RNA) and up to full length transcripts several kilobases long.
• Thorough testing on large ENCODE datasets:
  STAR was used to map 64 Billion reads of long RNA-seq and 16 Billion reads of short RNA-seq, and will be used to map RNA-seq data in the next ENCODE phase.

STAR requires ~30GB of RAM for mapping to the human genome (could be reduced to 16GB in the "sparse" mode with some speed loss).

More information can be found in out recent paper.
If you decide to try it out, please download one of the latest STAR releases.
I will be happy to answer any questions via SEQanswers, STAR discussion forum, or by e-mail:[email protected]

Cheers
Alex