-
Notifications
You must be signed in to change notification settings - Fork 47
/
Readme.txt
executable file
·150 lines (100 loc) · 7.89 KB
/
Readme.txt
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
SURPI
v1.0.7
April 2014
Note: For the most up to date version of the SURPI source code, go to this website: http://chiulab.ucsf.edu/surpi.
SURPI has been tested on Ubuntu 12.04. It will likely function properly on other Linux distributions, but this has not been tested.
Hardware Requirements
SURPI requires a machine with high RAM in order to run efficiently. This is mainly due to SNAP, which gains its speed by loading the reference databases completely into RAM. We’ve run SURPI successfully on machines with 60.5GB RAM. SURPI will use all cores on a machine by default, though the number of cores used can be adjusted within the config file. Much of SURPI is parallelized, so it benefits from using as many cores as possible.
The steps to install SURPI on a machine are as follows:
1. Install all software dependencies
2. Decompress SURPI.tar.gz, and place all files into a directory included in your $PATH
3. Create the databases
4. Customize certain SURPI files as shown below
5. Run SURPI
1. Install Software Dependencies
The below software must be installed for SURPI to function properly.
• fastQValidator
http://genome.sph.umich.edu/wiki/FastQValidator
• Minimo (v1.6)
http://sourceforge.net/projects/amos/files/amos/3.1.0/
• Abyss (v1.3.5)
http://www.bcgsc.ca/platform/bioinfo/software/abyss
• RAPSearch (v2.12)
http://omics.informatics.indiana.edu/mg/RAPSearch2/
• seqtk (v 1.0r31)
https://github.com/lh3/seqtk
• SNAP (v0.15)
http://snap.cs.berkeley.edu
• gt (v1.5.1)
http://genometools.org/index.html
• fastq
https://github.com/brentp/bio-playground/tree/master/reads-utils
• fqextract
https://gist.github.com/drio/1168330
• cutadapt (v1.2.1)
https://code.google.com/p/cutadapt/
• prinseq-lite.pl
http://prinseq.sourceforge.net
• dropcache
http://stackoverflow.com/questions/13646925/allowing-a-non-root-user-to-drop-cache
2. Decompress SURPI
Use tar to decompress SURPI.tar.gz. Something like the following should work:
tar xvfz SURPI.tar.gz
3. Create Databases
The content and creation of the SNAP databases is documented in the paper in the Reference Databases section, which is duplicated below:
A 3.1 gigabase (Gb) human nucleotide database (human DB) was constructed from a combination of human genomic DNA (GRCh37 / hg19), rRNA (RefSeq), mRNA (RefSeq), and mitochondrial RNA (RefSeq) sequences in NCBI as of March of 2012. The bacterial nucleotide, viral nucleotide, and viral protein databases used by SURPI in fast mode (bacterial DB, viral nucleotide DB, and viral protein DB, respectively) were also constructed from sequences in NCBI as of March of 2012. The 3 Gb bacterial DB was constructed from all bacterial RefSeq entries and consisted of 348,922 unique accessioned sequences, each with a minimum length of 100 bp. The 1.4 Gb viral nucleotide DB included 1,193,607 entries and was constructed by searching for all viral sequences in the 42 Gb National Center for Biotechnology Information (NCBI) nt collection using the query term “viridae[Organism]” in BioPython. The viral protein DB was similarly constructed by extracting viral sequences from the NCBI nr DB collection. Index tables for SNAP (v0.15) were generated with an empirically determined default seed size of 20 for the human DB and viral nucleotide DB, and seed size of 16 for the bacterial DB. Index tables for RAPSearch (v2.09) were generated from the viral protein DB using default parameters.
To generate the National Center for Biotechnology Information (NCBI) nucleotide (nt) collection (NCBI nt DB) used by SURPI in comprehensive mode, the complete 42 Gb nucleotide collection (nt) was downloaded from NCBI in January of 2013. This collection consists of a comprehensive archive of sequences from multiple sources, including GenBank, European Molecular Biology Laboratory (EMBL), DNA Data Bank of Japan (DDBJ), and Protein Data Bank (PDB), and is the richest collection of annotated microbial sequence data publicly available. As SNAP uses 32-bit offsets in the reference genome during hashing, the aligner restricts the size of the reference genome to an absolute maximum of 2^32 bases, or ~4.2 Gb. Thus, the 42 Gb NCBI nt collection was first split into 29 sub-databases, each approximately 1.5 Gb in size. Each sub-database was then indexed separately by SNAP at default parameters with a seed size of 20. This generated 29 SNAP indexed databases, each approximately 27 GB in size, with the aggregate of all 29 databases referred to as the NCBI nt DB.
SNAP Databases
• Human DB
Comprehensive Mode
• NCBI nr DB
• Viral protein DB
• NCBI nt DB
Fast Mode
• Viral nt DB
• Bacterial DB
Taxonomy Databases
These databases can be created using the shell script titled: create_taxonomy_db.sh. Depending on your internet connection speed, and the speed of your system, this script may take several hours to complete creation of the database.
• gi_taxid_prot.db
• gi_taxid_nucl.db
• names_nodes_scientific.db
To use this script, execute the create_taxonomy_db.sh program. This script will download the necessary data from NCBI, and generate the above 3 databases.
4. SURPI file customization
Below are some notes on files that may need to be modified to run SURPI.
• cutadapt_quality.csh - specify location of /tmp folder
cutadapt_quality.csh defaults to using /tmp for temporary file storage. If using a system with limited space in this location, change the location to a directory with more storage space available.
• taxonomy_lookup_embedded.pl
Set database_directory to the location of the taxonomy databases created below.
• tweet.pl
SURPI has the ability to send out notifications via Twitter at various stages within the pipeline. If this feature is desired, you will need to set up a Twitter application within your account for this purpose. See https://dev.twitter.com/apps for more details.
Once an application has been set up, fill in the below parameters to the tweet.pl program.
consumer_key
consumer_secret
oauth_token
oauth_token_secret
• perl modules to install
Net::Twitter::Lite::WithAPIv1_1
Net::OAuth
5. Run SURPI
To run SURPI, execute the following in a directory containing your FASTQ input file.
1. This command will create the necessary config file to run SURPI:
SURPI.sh -z <INPUTFILE>
After typing the above line, a config file and a “go” file will be created. The config file will contain default values for many parameters - these parameters may need to be modified depending on your environment. The config file has descriptions of the options allowed by SURPI.
2. Once the config file has been customized, the SURPI pipeline can be initiated by typing in the name of the go file that was created. Below is an example (boldfaced text is inputted by the user):
sfederman@tribble:/data/inputfile/test$ ls -laF
total 750212
drwxrwxr-x 2 sfederman sfederman 4096 Jan 20 16:45 ./
drwxrwxr-x 11 sfederman sfederman 61440 Jan 20 16:45 ../
-rw-rw-r-- 1 sfederman sfederman 768143660 Jan 20 16:45 inputfile.fastq
sfederman@tribble:/data/inputfile/test$ SURPI.sh -z inputfile.fastq
inputfile.config generated. Please edit it to contain the proper parameters for your analysis.
go_ inputfile generated. Initiate the pipeline by running this program. (./go_inputfile)
sfederman@tribble:/data/inputfile/test$ ls -laF
total 750220
drwxrwxr-x 2 sfederman sfederman 4096 Jan 20 16:47 ./
drwxrwxr-x 11 sfederman sfederman 61440 Jan 20 16:45 ../
-rw-rw-r-- 1 sfederman sfederman 1976 Jan 20 16:47 inputfile.config
-rw-rw-r-- 1 sfederman sfederman 768143660 Jan 20 16:45 inputfile.fastq
-rwxrwxr-x 1 sfederman sfederman 84 Jan 20 16:47 go_inputfile*
sfederman@tribble:/data/inputfile/test$ ./go_inputfile &
Progression of the pipeline can be followed by monitoring the log file (titled inputfile.SURPI.log, in the above example). We have also find it useful to monitor the status of the pipeline with the program htop.