24
|
1 <tool id="plant_tribes_assembly_post_processor" name="Postprocess de novo assembly transcripts" version="0.4.0">
|
5
|
2 <description>into putative coding sequences</description>
|
0
|
3 <requirements>
|
14
|
4 <requirement type="package" version="0.4">plant_tribes_assembly_post_processor</requirement>
|
0
|
5 </requirements>
|
|
6 <stdio>
|
|
7 <!-- Anything other than zero is an error -->
|
|
8 <exit_code range="1:" />
|
6
|
9 <exit_code range=":-1" />
|
|
10 <!-- In case the return code has not been set propery check stderr too -->
|
0
|
11 <regex match="Error:" />
|
|
12 <regex match="Exception:" />
|
|
13 </stdio>
|
|
14 <command>
|
|
15 <![CDATA[
|
15
|
16 #import os
|
0
|
17 AssemblyPostProcesser
|
15
|
18 --transcripts '$input'
|
0
|
19 --prediction_method $prediction_method_cond.prediction_method
|
|
20 #if str($prediction_method_cond.prediction_method) == 'estscan':
|
15
|
21 --score_matrices '$score_matrices'
|
0
|
22 #end if
|
|
23 #if str($options_type.options_type_selector) == 'advanced':
|
|
24 #if str($options_type.target_gene_family_assembly_cond.target_gene_family_assembly) == 'yes':
|
24
|
25 --config_dir '$options_type.target_gene_family_assembly_cond.scaffold.fields.path'
|
15
|
26 --gene_family_search '$options_type.target_gene_family_assembly_cond.orthogroups'
|
21
|
27 --scaffold '$options_type.target_gene_family_assembly_cond.scaffold.fields.path'
|
15
|
28 --method '$options_type.target_gene_family_assembly_cond.method'
|
0
|
29 #end if
|
|
30 --gap_trimming $options_type.gap_trimming
|
|
31 #if str($options_type.strand_specific) == 'yes':
|
|
32 --strand_specific
|
|
33 #end if
|
|
34 #if str($options_type.dereplicate) == 'yes':
|
|
35 --dereplicate
|
|
36 #end if
|
|
37 --min_length $options_type.min_length
|
|
38 #end if
|
|
39 --num_threads \${GALAXY_SLOTS:-4}
|
25
|
40 >/dev/null
|
0
|
41 ]]>
|
|
42 </command>
|
|
43 <inputs>
|
|
44 <param name="input" format="fasta" type="data" label="De novo transcriptome assembly fasta file"/>
|
|
45 <conditional name="prediction_method_cond">
|
|
46 <param name="prediction_method" type="select" label="Prediction method for coding regions">
|
|
47 <option value="transdecoder" selected="true">TransDecoder</option>
|
|
48 <option value="estscan">ESTScan</option>
|
|
49 </param>
|
|
50 <when value="transdecoder" />
|
|
51 <when value="estscan">
|
|
52 <param name="score_matrices" format="smat" type="data" label="Scores matrices"/>
|
|
53 </when>
|
|
54 </conditional>
|
|
55 <conditional name="options_type">
|
|
56 <param name="options_type_selector" type="select" label="Options Configuration">
|
|
57 <option value="basic" selected="true">Basic</option>
|
|
58 <option value="advanced">Advanced</option>
|
|
59 </param>
|
|
60 <when value="basic" />
|
|
61 <when value="advanced">
|
|
62 <conditional name="target_gene_family_assembly_cond">
|
|
63 <param name="target_gene_family_assembly" type="select" label="Target gene family assembly?">
|
|
64 <option value="no" selected="true">No</option>
|
|
65 <option value="yes">Yes</option>
|
|
66 </param>
|
|
67 <when value="no" />
|
|
68 <when value="yes">
|
|
69 <param name="orthogroups" format="tabular" type="data" label="List of orthogroup identifiers for target gene families to assemble"/>
|
|
70 <param name="scaffold" type="select" label="Orthogroups or gene families proteins scaffold">
|
|
71 <options from_data_table="plant_tribes_scaffolds" />
|
10
|
72 <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
|
73 </param>
|
|
74 <param name="method" type="select" label="Protein clustering method">
|
|
75 <option value="gfam" selected="true">GFam</option>
|
|
76 <option value="orthofinder">OrthoFinder</option>
|
|
77 <option value="orthomcl">OrthoMCL</option>
|
|
78 </param>
|
|
79 </when>
|
|
80 </conditional>
|
3
|
81 <param name="gap_trimming" type="float" value="0.1" min="0" max="1.0" label="Remove sites in alignments with gaps of" help="0.1 removes sites with 90% gaps"/>
|
0
|
82 <param name="strand_specific" type="select" label="De novo transcriptome assembly was performed with strand-specific library?">
|
|
83 <option value="no" selected="true">No</option>
|
|
84 <option value="yes">Yes</option>
|
|
85 </param>
|
|
86 <param name="dereplicate" type="select" label="Remove repeated sequences in predicted coding regions?">
|
|
87 <option value="no" selected="true">No</option>
|
|
88 <option value="yes">Yes</option>
|
|
89 </param>
|
|
90 <param name="min_length" type="integer" value="200" label="Minimum sequence length of predicted coding regions"/>
|
|
91 </when>
|
|
92 </conditional>
|
|
93 </inputs>
|
|
94 <outputs>
|
|
95 <collection name="transcripts" type="list">
|
1
|
96 <discover_datasets pattern="__name__" directory="assemblyPostProcessing_dir" visible="true" ext="fasta" />
|
0
|
97 </collection>
|
|
98 </outputs>
|
|
99 <tests>
|
|
100 <test>
|
|
101 <param name="input" value="assembly.fasta" ftype="fasta" />
|
|
102 <param name="prediction_method" value="transdecoder" />
|
6
|
103 <output_collection name="transcripts" type="list">
|
15
|
104 <element name="transcripts.cds" file="transcripts.cds" ftype="fasta"/>
|
|
105 <element name="transcripts.cleaned.cds" file="transcripts.cleaned.cds" ftype="fasta"/>
|
16
|
106 <element name="transcripts.cleaned.pep" file="transcripts.cleaned.pep" ftype="fasta"/>
|
15
|
107 <element name="transcripts.pep" file="transcripts.pep" ftype="fasta"/>
|
6
|
108 </output_collection>
|
|
109 </test>
|
|
110 <test>
|
|
111 <param name="input" value="assembly.fasta" ftype="fasta" />
|
|
112 <param name="prediction_method" value="transdecoder" />
|
0
|
113 <param name="options_type_selector" value="advanced" />
|
7
|
114 <param name="gap_trimming" value="0.1" />
|
0
|
115 <param name="dereplicate" value="yes" />
|
7
|
116 <param name="min_length" value="200" />
|
0
|
117 <output_collection name="transcripts" type="list">
|
16
|
118 <element name="transcripts.cds" file="transcripts.cds" ftype="fasta"/>
|
|
119 <element name="transcripts.cleaned.cds" file="transcripts.cleaned.cds" ftype="fasta"/>
|
|
120 <element name="transcripts.cleaned.nr.cds" file="transcripts.cleaned.nr.cds" ftype="fasta"/>
|
|
121 <element name="transcripts.cleaned.nr.pep" file="transcripts.cleaned.nr.pep" ftype="fasta"/>
|
|
122 <element name="transcripts.cleaned.pep" file="transcripts.cleaned.pep" ftype="fasta"/>
|
|
123 <element name="transcripts.pep" file="transcripts.pep" ftype="fasta"/>
|
0
|
124 </output_collection>
|
|
125 </test>
|
|
126 </tests>
|
|
127 <help>
|
|
128 This tool is one of the PlantTribes' collection of automated modular analysis pipelines that utilize objective classifications of
|
5
|
129 complete protein sequences from sequenced plant genomes to perform comparative evolutionary studies. It postprocesses de novo assembly
|
|
130 transcripts into putative coding sequences and their corresponding amino acid translations, locally assembling targeted gene families.
|
0
|
131
|
|
132 -----
|
|
133
|
|
134 **Options**
|
|
135
|
|
136 * **Prediction method for coding regions** - The prediction method for coding regions; one of ESTScan or TransDecoder.
|
|
137 * **Scores matrices** - Scores matrices that reflect the codons preferences in the studied organisms.
|
5
|
138 * **Target gene family assembly?** - Select 'Yes' to target gene families to assemble.
|
|
139 * **List of orthogroup identifiers for target gene families to assemble** - History item with a list of orthogroup identifiers for target gene families to assemble.
|
23
|
140 * **Orthogroups or gene families proteins scaffold** - PlantTribes scaffolds data installed into Galaxy by the PlantTribes Scaffolds Download Data Manager tool.
|
5
|
141 * **Protein clustering method** - One of GFam (domain architecture based clustering), OrthoFinder (broadly defined clusters) or OrthoMCL (narrowly defined clusters).
|
|
142 * **Remove sites in alignments with gaps of** - Removes gappy sites in alignments (i.e. 0.1 removes sites with 90% gaps): [0.0 to 1.0]
|
|
143 * **De novo transcriptome assembly was performed with strand-specific library?** - Select 'Yes' if de novo transcriptome assembly was performed with strand-specific library.
|
0
|
144 * **Remove repeated sequences in predicted coding regions?** - Select yes to remove repeated sequences in predicted coding regions.
|
|
145 * **Minimum sequence length of predicted coding regions** - Set the minimum sequence length of predicted coding regions.
|
|
146
|
|
147 </help>
|
|
148 <citations>
|
24
|
149 <citation type="bibtex">
|
|
150 @unpublished{None,
|
|
151 author = {Eric Wafula},
|
|
152 title = {None},
|
|
153 year = {None},
|
|
154 url = {https://github.com/dePamphilis/PlantTribes}
|
|
155 }</citation>
|
|
156 <citation type="bibtex">
|
|
157 @published{Proc Int Conf Intell Syst Mol Biol,
|
|
158 author = {Iseli C, Jongeneel CV, Bucher P.},
|
|
159 title = {ESTScan: a program for detecting, evaluating, and reconstructing potential coding regions in EST sequences.},
|
|
160 year = {1999},
|
|
161 url = {http://estscan.sourceforge.net/}
|
|
162 }</citation>
|
|
163 <citation type="doi">10.1038/nprot.2013.084</citation>
|
|
164 <citation type="doi">10.1109/tcbb.2013.68</citation>
|
|
165 <citation type="bibtex">
|
|
166 @unpublished{None,
|
|
167 author = {None},
|
|
168 title = {HMMER 3.1+ hmmscan search sequence(s) against a profile database},
|
|
169 year = {2013},
|
|
170 url = {http://hmmer.org/}
|
|
171 }</citation>
|
|
172 <citation type="bibtex">
|
|
173 @published{Genome Research,
|
|
174 author = {Huang X, Madan A},
|
|
175 title = {CAP3: A DNA sequence assembly program},
|
|
176 year = {1999},
|
|
177 url = {http://seq.cs.iastate.edu/cap3.html}
|
|
178 }</citation>
|
0
|
179 </citations>
|
|
180 </tool>
|