|
0
|
1 <tool id="antismash" name="Secondary Metabolites" version="2.0.2.0">
|
|
4
|
2 <description>and Antibiotics Analysis (antiSMASH)</description>
|
|
|
3 <requirements>
|
|
|
4 <requirement type="package" version="3.0">hmmer</requirement>
|
|
|
5 <requirement type="package" version="2.2.28">blast+</requirement>
|
|
|
6 <requirement type="package" version="3.8.31">muscle</requirement>
|
|
17
|
7 <requirement type="package" version="1.4.0-post-1">straight.plugin</requirement>
|
|
|
8 <requirement type="package" version="0.1.2">helperlibs</requirement>
|
|
4
|
9 <requirement type="package" version="2.0.2">antismash</requirement>
|
|
|
10 </requirements>
|
|
|
11 <command>
|
|
15
|
12 #import os, glob
|
|
|
13 #set $outputfolder = os.path.join($html.files_path, 'antismash')
|
|
|
14 ln -s $infile #echo 'input_tempfile.' + $infile.ext#;
|
|
|
15 mkdir -p $outputfolder;
|
|
4
|
16 run_antismash.py
|
|
|
17 --cpus 4
|
|
13
|
18 --enable $types
|
|
15
|
19 --input-type 'nucl'
|
|
4
|
20 $smcogs
|
|
|
21 $clusterblast
|
|
|
22 $subclusterblast
|
|
|
23 $inclusive
|
|
|
24 $full_hmmer
|
|
|
25 $full_blast
|
|
0
|
26
|
|
15
|
27 #if str($pfam_database) != "None":
|
|
|
28 --pfamdir $pfam_database.fields.path
|
|
|
29 #end if
|
|
|
30
|
|
|
31 --disable-embl
|
|
|
32 --outputfolder $outputfolder
|
|
|
33
|
|
|
34
|
|
|
35 #echo 'input_tempfile.' + $infile.ext#
|
|
0
|
36
|
|
9
|
37 ## leave out the start and end features, it can be easily replaced with Galaxy tools
|
|
|
38 ##--from START Start analysis at nucleotide specified
|
|
15
|
39 ##--to END
|
|
|
40
|
|
|
41 ;
|
|
|
42 ## shuffling files to create the correct outputs for Galaxy
|
|
|
43
|
|
|
44 ## html output
|
|
16
|
45 #if 'html' in str($outputs).split(','):
|
|
15
|
46 cp #echo os.path.join($outputfolder, 'index.html')# $html;
|
|
|
47 #end if
|
|
|
48
|
|
|
49 ## gene clusters
|
|
16
|
50 #if 'geneclusterprots_tabular' in str($outputs).split(','):
|
|
15
|
51 cp #echo os.path.join($outputfolder, 'geneclusters.txt')# $geneclusterprots;
|
|
|
52 #end if
|
|
|
53
|
|
|
54 ##SVG images
|
|
16
|
55 #if 'archive_svgs' in str($outputs).split(','):
|
|
15
|
56 cd #echo os.path.join($outputfolder, 'svg')#;
|
|
16
|
57 tar cfz $archive_svgs *_all.svg genecluster*;
|
|
15
|
58 #end if
|
|
|
59
|
|
|
60 ##all files in a archive
|
|
16
|
61 #if 'archive' in str($outputs).split(','):
|
|
15
|
62 cd $outputfolder;
|
|
|
63 tar cf $archive *.zip;
|
|
|
64 #end if
|
|
|
65
|
|
|
66 ## genbank
|
|
16
|
67 #if 'gb' in str($outputs).split(','):
|
|
15
|
68 #for $gb_file in glob.glob( os.path.join($outputfolder, '*.gbk') ):
|
|
|
69 cat $gb_file >> $genbank;
|
|
|
70 #end for
|
|
|
71 #end if
|
|
0
|
72
|
|
|
73 </command>
|
|
|
74 <inputs>
|
|
|
75 <param name="infile" type="data" format="gb,embl" label="Nucleotide sequence file in GenBank or EMBL format"/>
|
|
|
76
|
|
11
|
77 <param name="clusterblast" type="boolean" label="BLAST identified clusters against known clusters" truevalue="--clusterblast" falsevalue="" checked="True" />
|
|
|
78 <param name="smcogs" type="boolean" label="analysis of secondary metabolism gene families (smCOGs)"
|
|
|
79 falsevalue="" truevalue="--smcogs" checked="True" />
|
|
|
80
|
|
|
81 <param name="full_blast" type="boolean" label="Run a whole-genome BLAST analysis" truevalue="--full-blast" falsevalue="" checked="False" />
|
|
|
82 <param name="subclusterblast" type="boolean" label="Subcluster Blast analysis" truevalue="--subclusterblast" falsevalue="" checked="false" />
|
|
|
83 <param name="full_hmmer" type="boolean" label="Run a whole-genome Pfam analysis" truevalue="--full-hmmer" falsevalue="" checked="false" />
|
|
0
|
84
|
|
|
85 <param name="inclusive" type="boolean" label="Use inclusive algorithm for cluster detection" truevalue="--inclusive" falsevalue="" checked="false" />
|
|
|
86
|
|
15
|
87 <param name="pfam_database" type="select" optional="true" label="Pfam database" help="Pfam Covariance models">
|
|
4
|
88 <options from_file="antismash.loc">
|
|
|
89 <column name="value" index="0"/>
|
|
|
90 <column name="name" index="1"/>
|
|
|
91 <column name="path" index="2"/>
|
|
|
92 </options>
|
|
|
93 </param>
|
|
|
94
|
|
0
|
95 <param name="types" type="select" display="checkboxes" multiple="true" label="Gene cluster types to search">
|
|
9
|
96 <option value="t1pks" selected="True">type I polyketide synthases</option>
|
|
|
97 <option value="t2pks" selected="True">type II polyketide synthases</option>
|
|
|
98 <option value="t3pks" selected="True">type III polyketide synthases</option>
|
|
|
99 <option value="t4pks" selected="True">type IV polyketide synthases</option>
|
|
|
100 <option value="transatpks" selected="True">trans-AT PKS</option>
|
|
|
101 <option value="nrps" selected="True">nonribosomal peptide synthetases</option>
|
|
|
102 <option value="terpene" selected="True">terpene synthases</option>
|
|
|
103 <option value="lantipeptide" selected="True">lantipeptides</option>
|
|
|
104 <option value="bacteriocin" selected="True">bacteriocins</option>
|
|
|
105 <option value="blactam" selected="True">beta-lactams</option>
|
|
|
106 <option value="amglyccycl" selected="True">aminoglycosides / aminocyclitols</option>
|
|
|
107 <option value="aminocoumarin" selected="True">aminocoumarins</option>
|
|
|
108 <option value="siderophore" selected="True">siderophores</option>
|
|
|
109 <option value="ectoine" selected="True">ectoines</option>
|
|
|
110 <option value="butyrolactone" selected="True">butyrolactones</option>
|
|
|
111 <option value="indole" selected="True">indoles</option>
|
|
|
112 <option value="nucleoside" selected="True">nucleosides</option>
|
|
|
113 <option value="phosphoglycolipid" selected="True">phosphoglycolipids</option>
|
|
|
114 <option value="oligosaccharide" selected="True">oligosaccharides</option>
|
|
|
115 <option value="furan" selected="True">furans</option>
|
|
|
116 <option value="hserlactone" selected="True">hserlactones</option>
|
|
|
117 <option value="thiopeptide" selected="True">thiopeptides</option>
|
|
|
118 <option value="phenazine" selected="True">phenazines</option>
|
|
|
119 <option value="phosphonate" selected="True">phosphonates</option>
|
|
15
|
120 <option value="other" selected="True">others</option>
|
|
|
121 </param>
|
|
|
122
|
|
|
123 <param name="outputs" type="select" multiple="true" label="Output formats">
|
|
|
124 <option value="geneclusterprots_fasta" selected="True">Gene cluster proteins (FASTA)</option>
|
|
|
125 <option value="geneclusterprots_tabular">Gene cluster proteins (Tabular)</option>
|
|
|
126 <option value="html" selected="True">Interactive summary page (html)</option>
|
|
|
127 <option value="archive_svgs">All clusters as image (compressed)</option>
|
|
|
128 <option value="archive">All files compressed</option>
|
|
|
129 <option value="gb">Annotated genome (GenBank)</option>
|
|
0
|
130 </param>
|
|
|
131
|
|
|
132 </inputs>
|
|
|
133 <outputs>
|
|
15
|
134 <data format="fasta" name="geneclusterprots" label="${tool.name} on ${on_string} (Gen Cluster Proteins)">
|
|
|
135 <filter>'geneclusterprots_fasta' in outputs</filter>
|
|
|
136 </data>
|
|
|
137 <data format="fasta" name="geneclusterprots" label="${tool.name} on ${on_string} (Gen Cluster Proteins)">
|
|
|
138 <filter>'geneclusterprots_tabular' in outputs</filter>
|
|
|
139 </data>
|
|
|
140 <data format="gzipped" name="archive" label="${tool.name} on ${on_string} (all files compressed)">
|
|
|
141 <filter>'archive' in outputs</filter>
|
|
|
142 </data>
|
|
|
143 <data format="gzipped" name="archive_svgs" label="${tool.name} on ${on_string} (SVG images)">
|
|
|
144 <filter>'archive_svgs' in outputs</filter>
|
|
|
145 </data>
|
|
|
146 <data format="html" name="html" label="${tool.name} on ${on_string} (html report)">
|
|
|
147 <filter>'html' in outputs</filter>
|
|
|
148 </data>
|
|
|
149 <data name="genbank" format="genbank" label="${tool.name} on ${on_string} (genbank)">
|
|
|
150 <filter>'gb' in outputs</filter>
|
|
0
|
151 </data>
|
|
|
152 </outputs>
|
|
|
153 <help>
|
|
|
154
|
|
|
155 .. class:: infomark
|
|
|
156
|
|
|
157 **What it does**
|
|
|
158
|
|
|
159 antiSMASH allows the rapid genome-wide identification, annotation and analysis of secondary metabolite biosynthesis gene clusters in bacterial and fungal genomes.
|
|
|
160 It integrates and cross-links with a large number of in silico secondary metabolite analysis tools that have been published earlier.
|
|
|
161
|
|
|
162
|
|
|
163 **Input**
|
|
|
164
|
|
11
|
165 The ideal input for antiSMASH is an annotated nucleotide file in Genbank format or EMBL format. If no annotation is available,
|
|
|
166 we recommend running your sequence through an annotation pipeline like RAST are one included in Galaxy.
|
|
|
167
|
|
|
168
|
|
|
169 There are several optional analyses that may or may not be run on your sequence.
|
|
|
170 Highly recommended is the Gene Cluster Blast Comparative Analysis, which runs BlastP using each amino acid sequence from a detected gene cluster as a
|
|
|
171 query on a large database of predicted protein sequences from secondary metabolite biosynthetic gene clusters, and pools the results to identify
|
|
|
172 the gene clusters that are most homologous to the gene cluster that was detected in your query nucleotide sequence.
|
|
|
173
|
|
|
174
|
|
|
175 Also available is the analysis of secondary metabolism gene families (smCOGs).
|
|
|
176 This analysis attempts to allocate each gene in the detected gene clusters to a secondary metabolism-specific gene
|
|
|
177 family using profile hidden Markov models specific for the conserved sequence region characteristic of this family.
|
|
|
178 Additionally, a phylogenetic tree is constructed of each gene together with the (max. 100) sequences of the smCOG seed alignment.
|
|
|
179
|
|
|
180
|
|
|
181 For the most thorough genome analysis, we provide genome-wide PFAM HMM analysis of all genes in the genome through modules of the CLUSEAN pipeline.
|
|
|
182 Of course, some regions important to secondary metabolism may have been missed in the gene cluster identification stage
|
|
|
183 (e.g. because they represent the biosynthetic pathway of a yet unknown secondary metabolite).
|
|
|
184 Therefore, when genome-wide PFAM HMM analysis is selected, the PFAM frequencies are also used to find all genome regions in which PFAM domains typical for secondary metabolism are overrepresented.
|
|
0
|
185
|
|
|
186
|
|
|
187 **References**
|
|
|
188
|
|
|
189 Marnix H. Medema, Kai Blin, Peter Cimermancic, Victor de Jager, Piotr Zakrzewski, Michael A. Fischbach, Tilmann Weber,
|
|
|
190 Rainer Breitling and Eriko Takano (2011). antiSMASH: Rapid identification, annotation and analysis of secondary metabolite biosynthesis gene clusters. Nucleic Acids Research, doi: 10.1093/nar/gkr466.
|
|
|
191
|
|
|
192 http://antismash.secondarymetabolites.org/help.html
|
|
|
193
|
|
|
194 </help>
|
|
|
195 </tool>
|
