From UABgrid Documentation
(added category bioinformatics)
m (→TopHat Tutorials)
|Line 125:||Line 125:|
Revision as of 13:02, 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.
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 firstname.lastname@example.org # #$ -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 email@example.com 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
- Login to the UAB Galaxy interface at: http://galaxy.uabgrid.uab.edu/
- Upload the sequences file(s) using the 'Get Data' tool. (data should be available in the History pane)
- Click 'NGS:RNA Analyasis' > 'Tophat' under the NGS Toolbox
- Select reference genome and RNA sequence file
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 can be found at: http://tophat.cbcb.umd.edu/manual.html
- TopHat Tutorials: http://icb.med.cornell.edu/wiki/index.php/Tophat_tutorial