{"id":1427,"date":"2014-04-09T16:51:33","date_gmt":"2014-04-09T16:51:33","guid":{"rendered":"http:\/\/ri.itservices.manchester.ac.uk\/csf-apps\/?page_id=1427"},"modified":"2015-07-09T09:49:08","modified_gmt":"2015-07-09T09:49:08","slug":"qiime","status":"publish","type":"page","link":"https:\/\/ri.itservices.manchester.ac.uk\/csf-apps\/software\/applications\/qiime\/","title":{"rendered":"QIIME"},"content":{"rendered":"

Overview<\/h2>\n

QIIME<\/a> (Quantitative Insights Into Microbial Ecology) is an open source software package for comparison and analysis of microbial communities, primarily based on high-throughput amplicon sequencing data (such as SSU rRNA) generated on a variety of platforms, but also supporting analysis of other types of data (such as shotgun metagenomic data).<\/p>\n

Versions 1.9.1 and 1.8.0 are installed on the CSF.<\/p>\n

QIMME is often used with USEARCH or vsearch. These pieces of software are installed separately on the CSF. If you wish to use them together you need to load QIIME (see below) and USEARCH<\/a> or<\/strong> vsearch<\/a>.<\/p>\n

Restrictions on use<\/h2>\n

None. Please read the citation instructions<\/a> for publishing work that has used Qiime.<\/p>\n

Set up procedure<\/h2>\n

To access the software you must first load the modulefile for the version you require:<\/p>\n

\r\nmodule load apps\/gcc\/qiime\/1.9.1\r\n<\/pre>\n
This will automatically load the apps\/binapps\/anaconda\/2.3.0<\/code> modulefile providing Anaconda Python.<\/p>\n
OR<\/p>\n
\r\nmodule load apps\/gcc\/qiime\/1.8.0\r\n<\/pre>\n
This will automatically load the apps\/binapps\/anaconda\/1.8.0<\/code> modulefile providing Anaconda Python.<\/p>\n
Running the application<\/h2>\n
Please do not run Qiime on the login node. Jobs should be submitted to the compute nodes via batch.<\/p>\n
Serial batch job submission<\/h3>\n
Make sure you have the modulefile loaded and that you have a created a directory for the job to run in and that that directory contains all the files and other directories needed by your job. Then in the man directory for the job create a batch submission script, for example:<\/p>\n
\r\n#!\/bin\/bash\r\n#$ -S \/bin\/bash\r\n#$ -cwd             # Job will run from the current directory\r\n#$ -V               # Job will inherit current environment settings\r\n\r\npick_open_reference_otus.py -i \/scratch\/$USER\/data\/combined_seqs.fna -m usearch \\\r\n-r \/scratch\/$USER\/Genes\/abc.fasta -o \/scratch\/$USER\/uclust_picked_otus \\\r\n-p \/scratch\/$USER\/qiime_parameters.txt --suppress_align_and_tree\r\n<\/pre>\n
where pick_open_reference_otus.py<\/em> is the Qiime tool you wish to use. To see a list of those available run the following on the login node:<\/p>\n
\r\nls $QIIME_HOME\/bin\r\n<\/pre>\n
Now submit the job with:<\/p>\n
qsub jobscript<\/pre>\n
Where jobscript<\/em> is the name of the submission script.<\/p>\n
Parallel batch job submission – 2 to 16 cores <\/h3>\n
As above, but we need to add a line to the submission script to ask for more than one core (in this example 4) and then we need to tell the software this as well – in the case of the component we are using this is done with -a -O $NSLOTS<\/code>:<\/p>\n
\r\n#!\/bin\/bash\r\n#$ -S \/bin\/bash\r\n#$ -cwd             # Job will run from the current directory\r\n#$ -V               # Job will inherit current environment settings\r\n#$ -pe smp.pe 4     # 4 core parallel job request\r\n\r\npick_open_reference_otus.py -a -O $NSLOTS -i \/scratch\/$USER\/data\/combined_seqs.fna -m usearch61 \\\r\n-r \/scratch\/$USER\/Genes\/abc.fasta -o \/scratch\/$USER\/uclust_picked_otus \\\r\n-p \/scratch\/$USER\/qiime_parameters.txt --suppress_align_and_tree\r\n<\/pre>\n
Note 1:<\/strong> Using $NSLOTS ensures that the number of cores you have requested (on the #$ -pe<\/code> line) and the number of cores the software actually tries to use match. If it doesn’t it can cause nodes to be overloaded (the software is using more cores than you reserved) or resources to be wasted (the software is using fewer cores than you reserved) which can then have negative impacts on the CSF service. DO NOT edit your config or jobscript to use more than 16 cores as it will overload nodes.<\/p>\n
Note 2:<\/strong> Not all tools in qiime can be run in parallel. Some may require options as in the example above, others may not. Some are named parallel_tool<\/em>.py<\/code> (where tool<\/em> is the qiime tool you are using). Please consult the qiime documentation for further details.<\/p>\n
Local Config File<\/h2>\n
The Qiime configuration options can be displayed by running print_qiime_config.py<\/code>. If you wish to modify any options you should create a local configuration file as follows:<\/p>\n
\r\ncp $QIIME_CONFIG_FP ~\/.qiime_config\r\n<\/pre>\n
and then edit your local ~\/.qiime_config<\/code> file.<\/p>\n
DO NOT make changes to this file that if you are not confident about how they will affect your job – errors in this config can cause issues in the batch system affecting not only your jobs, but those of other users. <\/p>\n
Further info<\/h2>\n