|
133
|
1 <tool id="plant_tribes_gene_family_classifier" name="GeneFamilyClassifier" version="@WRAPPER_VERSION@.0">
|
|
134
|
2 <description>classifies gene sequences into pre-computed orthologous gene family clusters</description>
|
|
131
|
3 <macros>
|
|
|
4 <import>macros.xml</import>
|
|
|
5 </macros>
|
|
|
6 <expand macro="requirements_gene_family_classifier" />
|
|
140
|
7 <command detect_errors="exit_code">
|
|
0
|
8 <![CDATA[
|
|
21
|
9 #if str($options_type.options_type_selector) == 'advanced':
|
|
51
|
10 #set specify_super_orthogroups_cond = $options_type.specify_super_orthogroups_cond
|
|
|
11 #set specify_super_orthogroups = $specify_super_orthogroups_cond.specify_super_orthogroups
|
|
17
|
12 #set create_orthogroup_cond = $options_type.create_orthogroup_cond
|
|
|
13 #set create_orthogroup = $create_orthogroup_cond.create_orthogroup
|
|
44
|
14 #set specify_single_copy_cond = $options_type.specify_single_copy_cond
|
|
46
|
15 #set specify_single_copy = $specify_single_copy_cond.specify_single_copy
|
|
51
|
16 #if str($specify_super_orthogroups) == 'yes':
|
|
|
17 #set specify_super_orthos = True
|
|
|
18 #set super_orthogroups = $specify_super_orthogroups_cond.super_orthogroups
|
|
|
19 #else:
|
|
|
20 #set specify_super_orthos = False
|
|
|
21 #end if
|
|
21
|
22 #if str($create_orthogroup) == 'yes':
|
|
44
|
23 #set create_ortho_sequences = True
|
|
21
|
24 #set create_corresponding_coding_sequences_cond = $create_orthogroup_cond.create_corresponding_coding_sequences_cond
|
|
|
25 #if str($create_corresponding_coding_sequences_cond.create_corresponding_coding_sequences) == 'yes':
|
|
|
26 #set create_corresponding_coding_sequences = True
|
|
|
27 #else:
|
|
|
28 #set create_corresponding_coding_sequences = False
|
|
44
|
29 #end if
|
|
21
|
30 #else:
|
|
|
31 #set create_ortho_sequences = False
|
|
23
|
32 #set create_corresponding_coding_sequences = False
|
|
21
|
33 #end if
|
|
52
|
34 #if str($specify_single_copy) == 'yes':
|
|
|
35 #set single_copy_orthogroup = True
|
|
|
36 #set single_copy_cond = $specify_single_copy_cond.single_copy_cond
|
|
|
37 #set single_copy = $single_copy_cond.single_copy
|
|
|
38 #else:
|
|
|
39 #set single_copy_orthogroup = False
|
|
|
40 #end if
|
|
13
|
41 #else:
|
|
47
|
42 #set single_copy_orthogroup = False
|
|
13
|
43 #set create_ortho_sequences = False
|
|
23
|
44 #set create_corresponding_coding_sequences = False
|
|
4
|
45 #end if
|
|
20
|
46
|
|
98
|
47 python $__tool_directory__/gene_family_classifier.py
|
|
100
|
48 --input '$input'
|
|
67
|
49 --scaffold '$scaffold.fields.path'
|
|
0
|
50 --method $method
|
|
13
|
51 --classifier $save_hmmscan_log_cond.classifier
|
|
69
|
52 --config_dir '$scaffold.fields.path'
|
|
30
|
53 --num_threads \${GALAXY_SLOTS:-4}
|
|
20
|
54
|
|
0
|
55 #if str($options_type.options_type_selector) == 'advanced':
|
|
51
|
56 #if specify_super_orthos:
|
|
|
57 --super_orthogroups $super_orthogroups
|
|
|
58 #end if
|
|
47
|
59 #if $single_copy_orthogroup:
|
|
68
|
60 #if str($single_copy) == 'custom':
|
|
70
|
61 #set single_copy_custom_cond = $single_copy_cond.single_copy_custom_cond
|
|
|
62 #set single_copy_custom = $single_copy_custom_cond.single_copy_custom
|
|
|
63 #if str($single_copy_custom) == 'no':
|
|
72
|
64 --single_copy_custom 'default'
|
|
70
|
65 #else:
|
|
|
66 --single_copy_custom '$single_copy_custom_cond.single_copy_custom_config'
|
|
|
67 #end if
|
|
68
|
68 #else:
|
|
133
|
69 #if str($single_copy_cond.single_copy_taxa):
|
|
126
|
70 --single_copy_taxa $single_copy_cond.single_copy_taxa
|
|
|
71 #end if
|
|
133
|
72 #if str($single_copy_cond.taxa_present):
|
|
126
|
73 --taxa_present $single_copy_cond.taxa_present
|
|
|
74 #end if
|
|
68
|
75 #end if
|
|
0
|
76 #end if
|
|
89
|
77 #if $create_ortho_sequences:
|
|
99
|
78 --orthogroup_fasta "true"
|
|
18
|
79 #if $create_corresponding_coding_sequences:
|
|
35
|
80 --coding_sequences '$create_corresponding_coding_sequences_cond.coding_sequences'
|
|
17
|
81 #end if
|
|
0
|
82 #end if
|
|
|
83 #end if
|
|
20
|
84
|
|
97
|
85 #if str($save_hmmscan_log_cond.classifier) == 'hmmscan' or str($save_hmmscan_log_cond.classifier) == 'both' and str($save_hmmscan_log_cond.save_hmmscan_log) == 'yes':
|
|
|
86 --save_hmmscan_log
|
|
|
87 --hmmscan_log '$hmmscan_log'
|
|
12
|
88 #end if
|
|
58
|
89 #if $create_ortho_sequences:
|
|
62
|
90 #if $create_corresponding_coding_sequences:
|
|
97
|
91 --output_ptorthocs '$output_ptorthocs'
|
|
102
|
92 --output_ptorthocs_dir '$output_ptorthocs.files_path'
|
|
62
|
93 #else:
|
|
97
|
94 --output_ptortho '$output_ptortho'
|
|
102
|
95 --output_ptortho_dir '$output_ptortho.files_path'
|
|
18
|
96 #end if
|
|
4
|
97 #end if
|
|
47
|
98 #if $single_copy_orthogroup:
|
|
58
|
99 #if $create_ortho_sequences:
|
|
97
|
100 --output_ptsco '$output_ptsco'
|
|
102
|
101 --output_ptsco_dir '$output_ptsco.files_path'
|
|
48
|
102 #end if
|
|
44
|
103 #end if
|
|
0
|
104 ]]>
|
|
|
105 </command>
|
|
|
106 <inputs>
|
|
121
|
107 <param name="input" format="fasta" type="data" label="Proteins fasta file"/>
|
|
|
108 <param name="scaffold" type="select" label="Gene family scaffold">
|
|
0
|
109 <options from_data_table="plant_tribes_scaffolds" />
|
|
39
|
110 <validator type="no_options" message="No PlantTribes scaffolds are available. Use the PlantTribes Scaffolds Download Data Manager tool in Galaxy to install and populate the PlantTribes scaffolds data table."/>
|
|
0
|
111 </param>
|
|
|
112 <param name="method" type="select" label="Protein clustering method">
|
|
|
113 <option value="gfam" selected="true">GFam</option>
|
|
|
114 <option value="orthofinder">OrthoFinder</option>
|
|
|
115 <option value="orthomcl">OrthoMCL</option>
|
|
|
116 </param>
|
|
13
|
117 <conditional name="save_hmmscan_log_cond">
|
|
121
|
118 <param name="classifier" type="select" label="Protein classifier">
|
|
13
|
119 <option value="blastp" selected="true">blastp</option>
|
|
130
|
120 <option value="hmmscan">hmmscan</option>
|
|
|
121 <option value="both">Both blastp and hmmscan</option>
|
|
13
|
122 </param>
|
|
|
123 <when value="blastp" />
|
|
|
124 <when value="hmmscan">
|
|
121
|
125 <param name="save_hmmscan_log" type="select" label="Save hmmscan log?">
|
|
13
|
126 <option value="no" selected="true">No</option>
|
|
|
127 <option value="yes">Yes</option>
|
|
|
128 </param>
|
|
|
129 </when>
|
|
|
130 <when value="both">
|
|
121
|
131 <param name="save_hmmscan_log" type="select" label="Save hmmscan log?">
|
|
13
|
132 <option value="no" selected="true">No</option>
|
|
|
133 <option value="yes">Yes</option>
|
|
|
134 </param>
|
|
|
135 </when>
|
|
|
136 </conditional>
|
|
0
|
137 <conditional name="options_type">
|
|
121
|
138 <param name="options_type_selector" type="select" label="Options configuration">
|
|
0
|
139 <option value="basic" selected="true">Basic</option>
|
|
|
140 <option value="advanced">Advanced</option>
|
|
|
141 </param>
|
|
|
142 <when value="basic" />
|
|
|
143 <when value="advanced">
|
|
51
|
144 <conditional name="specify_super_orthogroups_cond">
|
|
121
|
145 <param name="specify_super_orthogroups" type="select" label="Super orthogroups configuration">
|
|
51
|
146 <option value="no" selected="true">No</option>
|
|
|
147 <option value="yes">Yes</option>
|
|
|
148 </param>
|
|
|
149 <when value="no"/>
|
|
|
150 <when value="yes">
|
|
121
|
151 <param name="super_orthogroups" type="select" label="Clustering distance measure">
|
|
130
|
152 <option value="min_evalue" selected="true">minimum e-value</option>
|
|
|
153 <option value="avg_evalue">average e-value</option>
|
|
51
|
154 </param>
|
|
|
155 </when>
|
|
|
156 </conditional>
|
|
44
|
157 <conditional name="specify_single_copy_cond">
|
|
121
|
158 <param name="specify_single_copy" type="select" label="Single copy orthogroups configuration">
|
|
44
|
159 <option value="no" selected="true">No</option>
|
|
|
160 <option value="yes">Yes</option>
|
|
0
|
161 </param>
|
|
44
|
162 <when value="no"/>
|
|
|
163 <when value="yes">
|
|
|
164 <conditional name="single_copy_cond">
|
|
121
|
165 <param name="single_copy" type="select" label="Selection criterion">
|
|
130
|
166 <option value="taxa" selected="true">Global selection</option>
|
|
|
167 <option value="custom">Custom selection</option>
|
|
44
|
168 </param>
|
|
70
|
169 <when value="custom">
|
|
|
170 <conditional name="single_copy_custom_cond">
|
|
121
|
171 <param name="single_copy_custom" type="select" label="Custom selection configuration">
|
|
70
|
172 <option value="no" selected="true">No</option>
|
|
|
173 <option value="yes">Yes</option>
|
|
|
174 </param>
|
|
|
175 <when value="no"/>
|
|
|
176 <when value="yes">
|
|
121
|
177 <param name="single_copy_custom_config" format="txt" type="data" label="Custom selection file"/>
|
|
70
|
178 </when>
|
|
|
179 </conditional>
|
|
|
180 </when>
|
|
44
|
181 <when value="taxa">
|
|
133
|
182 <param name="single_copy_taxa" type="integer" optional="true" min="0" label="Minimum single copy taxa"/>
|
|
|
183 <param name="taxa_present" type="integer" optional="true" min="0" label="Minimum taxa present"/>
|
|
44
|
184 </when>
|
|
|
185 </conditional>
|
|
0
|
186 </when>
|
|
|
187 </conditional>
|
|
|
188 <conditional name="create_orthogroup_cond">
|
|
121
|
189 <param name="create_orthogroup" type="select" label="Orthogroups fasta configuration">
|
|
0
|
190 <option value="no" selected="true">No</option>
|
|
|
191 <option value="yes">Yes</option>
|
|
|
192 </param>
|
|
|
193 <when value="no" />
|
|
|
194 <when value="yes">
|
|
17
|
195 <conditional name="create_corresponding_coding_sequences_cond">
|
|
121
|
196 <param name="create_corresponding_coding_sequences" type="select" label="Orthogroups coding sequences">
|
|
17
|
197 <option value="no" selected="true">No</option>
|
|
|
198 <option value="yes">Yes</option>
|
|
|
199 </param>
|
|
|
200 <when value="no" />
|
|
|
201 <when value="yes">
|
|
121
|
202 <param name="coding_sequences" format="fasta" type="data" label="Coding sequences fasta file"/>
|
|
17
|
203 </when>
|
|
|
204 </conditional>
|
|
0
|
205 </when>
|
|
|
206 </conditional>
|
|
|
207 </when>
|
|
|
208 </conditional>
|
|
|
209 </inputs>
|
|
|
210 <outputs>
|
|
15
|
211 <data name="hmmscan_log" format="txt" label="Protein classification hmmscan.log on ${on_string}">
|
|
|
212 <filter>save_hmmscan_log_cond['classifier'] in ['hmmscan', 'both'] and save_hmmscan_log_cond['save_hmmscan_log'] == 'yes'</filter>
|
|
|
213 </data>
|
|
63
|
214 <data name="output_ptortho" format="ptortho" label="Gene family clusters on ${on_string}">
|
|
33
|
215 <filter>options_type['options_type_selector'] == 'advanced' and options_type['create_orthogroup_cond']['create_orthogroup'] == 'yes' and options_type['create_orthogroup_cond']['create_corresponding_coding_sequences_cond']['create_corresponding_coding_sequences'] == 'no'</filter>
|
|
31
|
216 </data>
|
|
63
|
217 <data name="output_ptorthocs" format="ptorthocs" label="Gene family clusters and corresponding coding sequences on ${on_string}">
|
|
31
|
218 <filter>options_type['options_type_selector'] == 'advanced' and options_type['create_orthogroup_cond']['create_orthogroup'] == 'yes' and options_type['create_orthogroup_cond']['create_corresponding_coding_sequences_cond']['create_corresponding_coding_sequences'] == 'yes'</filter>
|
|
9
|
219 </data>
|
|
63
|
220 <data name="output_ptsco" format="tabular" label="Single copy orthogroups on ${on_string}">
|
|
53
|
221 <filter>options_type['options_type_selector'] == 'advanced' and options_type['create_orthogroup_cond']['create_orthogroup'] == 'yes' and options_type['specify_single_copy_cond']['specify_single_copy'] == 'yes'</filter>
|
|
63
|
222 <change_format>
|
|
|
223 <when input="options_type.create_orthogroup_cond.create_corresponding_coding_sequences_cond.create_corresponding_coding_sequences" value="no" format="ptortho" />
|
|
|
224 <when input="options_type.create_orthogroup_cond.create_corresponding_coding_sequences_cond.create_corresponding_coding_sequences" value="yes" format="ptorthocs" />
|
|
|
225 </change_format>
|
|
52
|
226 </data>
|
|
30
|
227 <collection name="orthos" type="list">
|
|
12
|
228 <discover_datasets pattern="__name__" directory="geneFamilyClassification_dir" visible="false" ext="tabular" />
|
|
0
|
229 </collection>
|
|
|
230 </outputs>
|
|
|
231 <tests>
|
|
132
|
232 <!-- Need to figure out how to test using scaffolds data.
|
|
0
|
233 <test>
|
|
30
|
234 <param name="input" value="transcripts.cleaned.nr.pep" ftype="fasta" />
|
|
132
|
235 <param name="scaffold" value="plant_tribes/scaffolds/22Gv1.1"/>
|
|
30
|
236 <param name="method" value="orthomcl"/>
|
|
|
237 <param name="classifier" value="blastp"/>
|
|
0
|
238 <output_collection name="orthos" type="list">
|
|
132
|
239 <element name="proteins.blastp.22Gv1.1" file="proteins.blastp.22Gv1.1" ftype="tabular" compare="contains"/>
|
|
35
|
240 <element name="proteins.blastp.22Gv1.1.bestOrthos" file="proteins.blastp.22Gv1.1.bestOrthos" ftype="tabular"/>
|
|
|
241 <element name="proteins.blastp.22Gv1.1.bestOrthos.summary" file="proteins.blastp.22Gv1.1.bestOrthos.summary" ftype="tabular"/>
|
|
0
|
242 </output_collection>
|
|
|
243 </test>
|
|
132
|
244 -->
|
|
0
|
245 </tests>
|
|
|
246 <help>
|
|
121
|
247 This tool is one of the PlantTribes collection of automated modular analysis pipelines for comparative and evolutionary
|
|
|
248 analyses of genome-scale gene families and transcriptomes. This tool classifies gene coding sequences either produced by
|
|
|
249 the AssemblyPostProcessor tool or from an external source into pre-computed orthologous gene family clusters (orthogroups)
|
|
|
250 of a PlantTribes scaffold. Classified sequences are then assigned with the corresponding orthogroups’ metadata that includes
|
|
|
251 gene counts of backbone taxa, super clusters (super orthogoups) at multiple stringencies, and functional annotations from
|
|
134
|
252 sources such as Gene Ontology (GO), InterPro protein domains, TAIR, UniProtKB/TrEMBL, and UniProtKB/Swiss-Prot. Additionally,
|
|
|
253 sequences belonging to single/low-copy gene families that are mainly utilized in species tree inference can be determined.
|
|
0
|
254
|
|
|
255 -----
|
|
|
256
|
|
121
|
257 **Required options**
|
|
|
258
|
|
|
259 * **Proteins fasta file** - proteins fasta file either produced by the AssemblyPostProcessor tool or an external source selected from your history.
|
|
|
260 * **Gene family scaffold** - one of the PlantTribes gene family scaffolds [2-4] installed into Galaxy by the PlantTribes Scaffold Data Manager tool.
|
|
|
261 * **Protein clustering method** - gene family scaffold protein clustering method as described in the AssemblyPostProcessor tool.
|
|
130
|
262 * **Protein classifier** - classifier to assign protein sequences into a specified scaffold orthogroups. PlantTribes implements three classification approaches; blastp (faster)[5], hmmscan (slower but more sensitive assignment of divergent homologs)[6], and both blastp and hmmscan (disagreements resolved in favor of hmmscan; more exhaustive).
|
|
121
|
263
|
|
|
264 **Other options**
|
|
|
265
|
|
130
|
266 * **Super orthogroups configuration** - select ‘Yes’ to enable super orthogroups configuration options. Super orthogroups[7] are constructed through a second iteration of MCL clustering to connect distant, but potentially related orthogroup clusters.
|
|
121
|
267
|
|
130
|
268 * **Clustering distance measure** - distance measure used in merging orthogroup clusters into super orthogroup clusters. PlantTribes pre-computed super orthogroups are based on the minimum and average blastp e-value between all pairs of scaffold orthogroups used as the input matrix for MCL clustering algorithm[8].
|
|
121
|
269
|
|
|
270 * **Single copy orthogroups configuration** - select ‘Yes’ to enable single/low-copy orthogroups selection configuration options.
|
|
0
|
271
|
|
130
|
272 * **Selection criterion** - single/low-copy orthogroups selection criterion. PlantTribes provides custom and global selection criteria for selecting user-defined single/low-copy scaffold orthogoups.
|
|
|
273
|
|
133
|
274 * **Global selection configuration** - the upper limit values of the following two parameters vary depending on the selected gene family scaffold, and the tool will produce an error if the value exceeds the number of species in the circumscribed scaffold.
|
|
130
|
275
|
|
|
276 * **Minimum single copy taxa** - minimum number of taxa with single copy genes in the orthogroup.
|
|
|
277 * **Minimum taxa present** - minimum number of taxa present in the orthogroup.
|
|
121
|
278
|
|
135
|
279 * **Custom selection configuration** - select ‘Yes’ to enable selection of a single copy configuration file. Scaffold configuration templates (.singleCopy.config) of how to customize single/low-copy orthogroups selection can be found in the scaffold data installed into Galaxy via the PlantTribes Scaffolds Download Data Manager tool, and also available at the PlantTribes GitHub `repository`_. Single/low-copy settings shown in these templates are used as defaults if ‘No’ is selected.
|
|
|
280
|
|
121
|
281 * **Custom selection file** - select a single/low-copy customized configuration file from your history.
|
|
|
282
|
|
|
283 * **Orthogroups fasta configuration** - select ‘Yes’ to create proteins orthogroups fasta files for the classified sequences.
|
|
|
284
|
|
138
|
285 * **Orthogroups coding sequences** - select ‘Yes’ to create corresponding coding sequences orthogroup fasta files for the classified protein sequences. Requires coding sequences fasta file corresponding to the proteins fasta file to be selected from your history.
|
|
121
|
286
|
|
|
287 * **Coding sequences fasta file** - select coding sequences fasta file corresponding to the proteins fasta file from your history.
|
|
|
288
|
|
137
|
289 .. _repository: https://github.com/dePamphilis/PlantTribes/tree/master/config
|
|
|
290
|
|
0
|
291 </help>
|
|
|
292 <citations>
|
|
131
|
293 <expand macro="citation1" />
|
|
134
|
294 <citation type="bibtex">
|
|
|
295 @article{Sasidharan2012,
|
|
|
296 journal = {Nucleic Acids Research},
|
|
|
297 author = {2. Sasidharan R, Nepusz T, Swarbreck D, Huala E, Paccanaro A},
|
|
|
298 title = {GFam: a platform for automatic annotation of gene families},
|
|
|
299 year = {2012},
|
|
|
300 pages = {gks631},}
|
|
|
301 </citation>
|
|
|
302 <citation type="bibtex">
|
|
|
303 @article{Li2003,
|
|
|
304 journal = {Genome Research}
|
|
|
305 author = {3. Li L, Stoeckert CJ, Roos DS},
|
|
|
306 title = {OrthoMCL: identification of ortholog groups for eukaryotic genomes},
|
|
|
307 year = {2003},
|
|
|
308 volume = {13},
|
|
|
309 number = {9},
|
|
|
310 pages = {2178-2189},}
|
|
|
311 </citation>
|
|
|
312 <citation type="bibtex">
|
|
|
313 @article{Emms2015,
|
|
|
314 journal = {Genome Biology}
|
|
|
315 author = {4. Emms DM, Kelly S},
|
|
|
316 title = {OrthoFinder: solving fundamental biases in whole genome comparisons dramatically improves orthogroup inference accuracy},
|
|
|
317 year = {2015},
|
|
|
318 volume = {16},
|
|
|
319 number = {1},
|
|
|
320 pages = {157},}
|
|
|
321 </citation>
|
|
121
|
322 <citation type="bibtex">
|
|
|
323 @article{Altschul1990,
|
|
|
324 journal = {Journal of molecular biology}
|
|
|
325 author = {5. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ},
|
|
|
326 title = {Basic local alignment search tool},
|
|
|
327 year = {1990},
|
|
|
328 volume = {215},
|
|
|
329 number = {3},
|
|
|
330 pages = {403-410},}
|
|
|
331 </citation>
|
|
|
332 <citation type="bibtex">
|
|
|
333 @article{Eddy2009,
|
|
|
334 journal = {Genome Inform},
|
|
|
335 author = {6. Eddy SR},
|
|
|
336 title = {A new generation of homology search tools based on probabilistic inference},
|
|
|
337 year = {2009},
|
|
|
338 volume = {23},
|
|
|
339 number = {1},
|
|
|
340 pages = {205-211},}
|
|
|
341 </citation>
|
|
|
342 <citation type="bibtex">
|
|
130
|
343 @article{Wall2008,
|
|
|
344 journal = {Nucleic Acids Research},
|
|
|
345 author = {7. Wall PK, Leebens-Mack J, Muller KF, Field D, Altman NS},
|
|
|
346 title = {PlantTribes: a gene and gene family resource for comparative genomics in plants},
|
|
|
347 year = {2008},
|
|
|
348 volume = {36},
|
|
|
349 number = {suppl 1},
|
|
|
350 pages = {D970-D976},}
|
|
|
351 </citation>
|
|
|
352 <citation type="bibtex">
|
|
121
|
353 @article{Enright2002,
|
|
|
354 journal = {Nucleic acids research},
|
|
130
|
355 author = {8. Enright AJ, Van Dongen S, Ouzounis CA},
|
|
121
|
356 title = {n efficient algorithm for large-scale detection of protein families},
|
|
|
357 year = {2002},
|
|
|
358 volume = {30},
|
|
|
359 number = {7},
|
|
|
360 pages = {1575-1584},}
|
|
|
361 </citation>
|
|
0
|
362 </citations>
|
|
|
363 </tool>
|
