TopHat: Difference between revisions

From Cheaha
Jump to navigation Jump to search
Line 125: Line 125:




[[Category: NGS]][[Category:Bioinformatics]]
[[Category: NGS]][[Category:Software]][[Category:Bioinformatics]]

Revision as of 19:06, 18 June 2013

TopHat is a fast splice junction mapper for RNA-Seq reads. It aligns RNA-Seq reads to mammalian-sized genomes using the ultra high-throughput short read aligner Bowtie, and then analyzes the mapping results to identify splice junctions between exons.

TopHat is a collaborative effort between the University of Maryland Center for Bioinformatics and Computational Biology and the University of California, Berkeley Departments of Mathematics and Molecular and Cell Biology.

The latest stable release on Cheaha is 1.2.0 as part of the Galaxy module. TopHat can run through the UAB Galaxy interface or directly on the cluster by submitting jobs to SGE.

Using TopHat

TopHat is a free software that runs on Linux and Mac OS X.

TopHat on your Desktop

TopHat can be downloaded and installed on your desktop from http://tophat.cbcb.umd.edu/

Installation instructions for TopHat are available at: http://tophat.cbcb.umd.edu/tutorial.html

TopHat on Cheaha

TopHat is pre-installed on the Cheaha research computing system. This allows users to run TopHat directly on the cluster without any need to install software.

There are two methods of using TopHat on Cheaha: 1. Direct use on the cluster using SGE submit scripts. 2. TopHat use through the UAB Galaxy interface.

Direct Use through SGE submit scripts

These instructions provide an example of how to create and submit a TopHat job on Cheaha.

First, create the working directory for the job- Replace 'USERNAME' with the account associated username on Cheaha. You can create any directory to run your job. It is recommended that the job directory be on the scratch (i.e. lustre filesystem) instead of the user home directory.

$ mkdir -p /lustre/scratch/USERNAME/jobs/tophat 
$ cd /lustre/scratch/USERNAME/jobs/tophat

Next, Copy all the files required for TopHat to the working directory.

TopHat requires at least the reference index database and the sequence reads files. If the Genome index database is not already built, you build it using the following steps. Load the appropriate modules file to load TopHat and Bowtie.

$ module load module load galaxy/galaxy
$ bowtie-build /lustre/scratch/USERNAME/jobs/tophat/example.fa example_bowtie_idx

This builds the Bowtie index for the genome 'example'. You can also download the pre-built index file from other external sources.

Next, create a job submit script as shown below called 'tophatSubmit', make sure to edit the following parameters:

* s_rt to an appropriate soft wall time limit
* h_rt to the maximum wall time for your job
* -N - job name
* -M - user email
* -pe smp numberOfProcessors  (-pe smp 8 - run the code in parallel on 8  processors on the same node of  Cheaha. 
   Also, -pe smp 12 is the maximum available number of processors)
* -l vf to the maximum memory needed for each task
* In the following example the genome is 'example', the reads are read1,read2, output will be written to testTophat/ folder. 
  The -r parameter is the Mean Inner Distance between Mate Pairs for pair-ended library. 
  This parameter can be left out for a single-ended library.
#!/bin/bash
#$ -S /bin/bash
#
# Execute in the current working directory
#$ -cwd
#
# Job runtime (10 hours)
#$ -l h_rt=10:00:00,s_rt=9:55:00
#$ -j y
#
# Job Name and email
#$ -N tophat
#$ -M username@uab.edu
#
#$ -pe smp 8
# Load the appropriate module(s)
module load galaxy/galaxy
#
#$ -V
# Amount of Memory needed (RAM) 3G
#$ -l vf=3G

 tophat -p  8 -r 20 -o testTophat/ example read1 read2

Submit the script to the scheduler with

$ qsub tophatSubmit

The output will be

Your job 8078636 ("tophat") has been submitted

You can check the status of the jobs using the 'qstat' command

qstat -u $USER

job-ID  prior   name       user         state submit/start at     queue                          slots ja-task-ID 
-----------------------------------------------------------------------------------------------------------------
8078636 0.52111 tophat     tanthony     r     08/03/2011 11:42:41 sipsey.q@sipsey-compute-2-15.l       8        

The job output can be found in the output directory specified in the submit script. You can view the output using samtools.

$ cd testTophat
$ samtools view accepted_hits.bam

Using the UAB Galaxy interface

  1. Login to the UAB Galaxy interface at: http://galaxy.uabgrid.uab.edu/
  2. Upload the sequences file(s) using the 'Get Data' tool. (data should be available in the History pane)
  3. Click 'NGS:RNA Analyasis' > 'Tophat' under the NGS Toolbox
  4. Select reference genome and RNA sequence file
  5. Execute

Galaxy support

In order to facilitate interaction among UAB Galaxy users, share experience, and provide peer-support we have established a galaxy-users group. To join this group and participate in email discussions please subscribe to the galaxy-user group. On-line archives of these discussions are available here. Please note, the email discussions are a public forum. You are advised to only post information you are authorized to share and comfortable with being public.

TopHat Tutorials