SSUsearch is pipeline for identify SSU rRNA gene and use them for diversity analysis. THe pipeline requires HMMER3.1, mothur, RDP mcclust, and python numpy, pandas, scipy, matplotlib, and screed package. The pipeline has two key features:
-
unsupervised (OTU based) community analysis with shotgun data
-
Scalibility: 5 Gb peak mem and 5 CPU hours on about 40 Gb data.
The manuscript has been accepted in AEM (doi: 10.1128/AEM.02772-15). It is also available on research gate.
The tutorials are written in ipython notebook. The easiest way to run it is using amazon EC2 instances with ami ami-7c82af16 and add security groups https. Add more storage according to your data size (the default is 20 Gb). Here is tutorial on how to setup EC2 instances. Notebooks could be accessed through https using browser (chrome or firefox NOT safari). Briefly, connect into your machine by using "https://" plus your machine name or IP address, and accept the “broken certificate” message. Password to access https is openscience. There are some introduction here.
If you want to run the notebooks on your computer, ipython (v4) with notebook dependencies are needed.
If you already have Python:
pip install "ipython[notebook]"
If you are new to Python, go to http://continuum.io/downloads and copy the link for Linux installer (pay attention to 64 or 32 bit version based on your cpu type). This tutorial use 64 bit version. At the time of this writing the 64 bit installer is here.
wget -c https://3230d63b5fc54e62148e-c95ac804525aac4b6dba79b00b39d1d3.ssl.cf1.rackcdn.com/Anaconda-2.2.0-Linux-x86_64.sh
bash Anaconda-2.2.0-Linux-x86_64.sh
Accept the license. Choose an installation directory (default is fine). Add to default path at the end of installation.
Then
source ~/.bashrc
Ipython is ready to go :)
Download SSUsearch:
cd ~/Desktop
git clone https://github.com/jiarong/SSUsearch
Go to notebook directory:
cd SSUsearch/notebooks-pc-linux
Start ipython notebook:
ipython notebook
In your default web browser, you will see a page showing several *.ipynb links.
First open overview.ipynb, and follow the instructions there. You will need to change the data directory settings for your own data.
If you prefer to run the tutorial in terminal, please go to http://microbial-ecology-protocols.readthedocs.org/en/latest/SSUsearch/overview.html
There is also a Makefile implemention shown below. It is great for automation but not as easy to read and modify. Plus I am not interested in maintaining api or wrapers. I keep it here mainly for some early users already using it. Thus the above notebooks are recommended to run the pipeline for new users.
This pipeline requires: HMMER3.1, mothur, RDP mcclust, and python pandas, scipy, matplotlib, and screed package. Following steps should work for linux machines.
git clone https://github.com/jiarong/SSUsearch.git
cd SSUsearch
mkdir -p external_tools; cd external_tools
wget http://selab.janelia.org/software/hmmer3/3.1b1/hmmer-3.1b1-linux-intel-x86_64.tar.gz -O hmmer-3.1b1-linux-intel-x86_64.tar.gz
tar -xzvf hmmer-3.1b1-linux-intel-x86_64.tar.gz
HMMSEARCH_BIN=$(readlink -f hmmer-3.1b1-linux-intel-x86_64/binaries/hmmsearch)
wget http://www.mothur.org/w/images/8/88/Mothur.cen_64.zip -O mothur.zip
unzip mothur.zip
MOTHUR_BIN=$(readlink -f mothur/mothur)
# RDP McClust
# Updated version can be found github: https://github.com/rdpstaff/RDPTools.git
# Here an older version is used
wget https://s3.amazonaws.com/ssusearchdb/Clustering.tgz
tar -xzvf Clustering.tgz
MCCLUST_JAR=$(readlink -f Clustering/dist/Clustering.jar)
# Python packages
# Use virtualenv (assuming virtualenv installed). virtualenv installation guide: https://virtualenv.pypa.io/en/latest/virtualenv.html#installation
virtualenv ssusearch_pyenv
source ssusearch_pyenv/bin/activate
pip install numpy pandas scipy screed matplotlib
# Or use [anaconda](http://continuum.io/downloads) as described above
# test if tools are properly installed
cd ..
make -f Makefile tool_check Hmmsearch=$HMMSEARCH_BIN Mothur=$MOTHUR_BIN Mcclust_jar=$MCCLUST_JAR
In the above, variables are set in command line arguments. It is better to modify the value of these variables to absolute path of binaries in Makefile, so there is no need to put them in argument in future. If tools are installed system wide, just use hmmsearch and mothur for the values of these three variables.
Several databases and hidden markov models are needed:
curl -O https://s3.amazonaws.com/ssusearchdb/SSUsearch_db.tgz
tar -xzvf SSUsearch_db.tgz
Parameters related to databases are Ali_template, Copy_db, Gene_db_cc, Gene_db, Gene_model_org, Gene_tax_cc, Gene_tax, and Hmm. First part of each database files its variable in Makefile. e.g. Ali_template.silva_lsu.fasta is for Ali_template. Set values of these parameters to absolute paths of their database in Makefile. Again, we can also set variables as commandline arguments. The default value should be OK for this tutorial if you follow the above steps.
Other than wgs data files, a design file is also needed for diversity analysis in mothur. Some description can be found in mothur wiki. Use the first column as sample tag and second as treatment. Sample tag is the first part of wgs data file names. e.g. For M1.fa.gz, the tag is M1.
# download a test dataset and its design file
wget https://s3.amazonaws.com/ssusearchdb/test.tgz
tar -xzvf test.tgz
# check wgs data files
ls test/data/*.fa
# check design file
cat test/SS.design
# run the pipeline
make -f Makefile Seqfiles="test/data/1c.fa test/data/1d.fa test/data/2c.fa test/data/2d.fa" Design=test/SS.design Script_dir=./scripts SSUsearch_db=./SSUsearch_db Method=ssusearch_no_qc Hmmsearch=$HMMSEARCH_BIN Mothur=$MOTHUR_BIN Mcclust_jar=$MCCLUST_JAR
Taxonomy results are here:
ls test/data/*.fa.ssu.out/*.taxonomy
By default, V4 will be used for OTU based analysis. The diversity analysis can be found at:
# no copy correction
ls diversity.ssu
# copy correction
ls diversity_cc.ssu
This pipeline includes PCoA, beta-diversity indces, weighted UNIFRAC and AMOVA. Familiarity with mothur are needed to find and plot the specific results. With SS.shared, SS.names, SS.groups, SS.design, most analysis in mothur can be done.
Makefile includes many variables for steps in pipeline. To see a full list of variables in Makefile:
make -f Makefile help
The soil datasets used in the paper are pair end merged long reads (longest is ~300 bp) and the default Start, End, and Len_cutoff are set for those datasets. For datasets have 100bp reads, Start=577, End=657, Len_cutoff=75 is recommended. Rule of thumb is to pick a region with more reads with larger overlap. Details are in "Testing the effect of target region size and variable region on clustering" of the paper.