|
1
|
1 <tool id="plant_tribes_kaks_analysis" name="Perform orthologous or paralogous ks analyses" version="0.7.0">
|
|
|
2 <description>of coding sequences and amino acid sequences</description>
|
|
0
|
3 <requirements>
|
|
|
4 <requirement type="package" version="0.7">plant_tribes_kaks_analysis</requirement>
|
|
|
5 </requirements>
|
|
|
6 <stdio>
|
|
|
7 <!-- Anything other than zero is an error -->
|
|
|
8 <exit_code range="1:" />
|
|
|
9 <exit_code range=":-1" />
|
|
|
10 <!-- In case the return code has not been set propery check stderr too -->
|
|
|
11 <regex match="Error:" />
|
|
|
12 <regex match="Exception:" />
|
|
|
13 </stdio>
|
|
|
14 <command>
|
|
|
15 <![CDATA[
|
|
|
16 #if str($options_type.options_type_selector) == 'advanced':
|
|
|
17 #set codeml_ctl_file_cond = $options_type.codeml_ctl_file_cond
|
|
|
18 #set codeml_ctl_file_select = $codeml_ctl_file_cond.codeml_ctl_file_select
|
|
|
19 #set fit_components_cond = $options_type.fit_components_cond
|
|
|
20 #set fit_components = $fit_components_cond.fit_components
|
|
|
21 #end if
|
|
|
22 KaKsAnalysis
|
|
|
23 --config_dir '$scaffold.fields.path'
|
|
|
24 --num_threads \${GALAXY_SLOTS:-4}
|
|
|
25 --coding_sequences_species_1 '$coding_sequences_species_1'
|
|
|
26 --proteins_species_1 '$proteins_species_1'
|
|
|
27 --comparison $comparison
|
|
3
|
28 #if str($comparison_cond.comparison) == 'orthologs':
|
|
0
|
29 --coding_sequences_species_2 '$comparison_cond.coding_sequences_species_2'
|
|
|
30 --proteins_species_2 '$comparison_cond.proteins_species_2'
|
|
|
31 #end if
|
|
|
32 #if str($options_type.options_type_selector) == 'advanced':
|
|
|
33 --min_coverage $min_coverage
|
|
|
34 --recalibration_rate $recalibration_rate
|
|
|
35 #if str($codeml_ctl_file_select) == 'yes':
|
|
|
36 --codeml_ctl_file '$codeml_ctl_file_cond.codeml_ctl_file'
|
|
|
37 # No else block needed here because the default codeml_ctl config
|
|
|
38 # will be used if the --codeml_ctl_file flag is missing.
|
|
|
39 #end if
|
|
|
40 #if str($fit_components) == 'yes':
|
|
|
41 --num_of_components $fit_components_cond.num_of_components
|
|
|
42 --min_ks $fit_components_cond.min_ks
|
|
|
43 --max_ks $fit_components_cond.max_ks
|
|
|
44 #end if
|
|
|
45 #end if
|
|
|
46 ]]>
|
|
|
47 </command>
|
|
|
48 <inputs>
|
|
2
|
49 <param name="scaffold" type="select" label="Orthogroups or gene families proteins scaffold">
|
|
|
50 <options from_data_table="plant_tribes_scaffolds" />
|
|
|
51 <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."/>
|
|
|
52 </param>
|
|
0
|
53 <param name="coding_sequences_species_1" format="fasta" type="data" label="Coding sequences (CDS) fasta file for the species" />
|
|
3
|
54 <param name="proteins_species_1" format="fasta" type="data" label="Aamino acids (proteins) sequences fasta file for the species" />
|
|
0
|
55 <conditional name="comparison_cond">
|
|
|
56 <param name="comparison" type="select" label="Select method for pairwise sequence comparison to determine homolgous pairs" help="Cross species comparison requires selection of inputs for second species">
|
|
|
57 <option value="paralogs" selected="true">Self species comparison</option>
|
|
|
58 <option value="orthologs">Cross species comparison</option>
|
|
|
59 </param>
|
|
|
60 <when value="paralogs" />
|
|
|
61 <when value="orthologs">
|
|
|
62 <param name="coding_sequences_species_2" format="fasta" type="data" label="Coding sequences (CDS) fasta file for the second species" />
|
|
3
|
63 <param name="proteins_species_2" format="fasta" type="data" label="Aamino acids (proteins) sequences fasta file for the second species" />
|
|
0
|
64 </when>
|
|
|
65 </conditional>
|
|
|
66 <conditional name="options_type">
|
|
|
67 <param name="options_type_selector" type="select" label="Options Configuration">
|
|
|
68 <option value="basic" selected="true">Basic</option>
|
|
|
69 <option value="advanced">Advanced</option>
|
|
|
70 </param>
|
|
|
71 <when value="basic" />
|
|
|
72 <when value="advanced">
|
|
|
73 <param name="min_coverage" type="float" value="0.5" min="0.3" max="1.0" label="Minimum sequence pairwise coverage length between homologous pairs" />
|
|
|
74 <param name="recalibration_rate" type="float" value="0.0" min="0.0" label="Predetermined evolutionary rate for recalibrating synonymous subsitutions (ks) of species" />
|
|
|
75 <conditional name="codeml_ctl_file_cond">
|
|
|
76 <param name="codeml_ctl_file_select" type="select" label="Select PAML codeml control file?" help="Used for ML analysis of protein-coding DNA sequences using codon substitution models, select No to use the default control file">
|
|
|
77 <option value="no" selected="true">No</option>
|
|
|
78 <option value="yes">Yes</option>
|
|
|
79 </param>
|
|
|
80 <when value="no" />
|
|
|
81 <when value="yes">
|
|
|
82 <param name="codeml_ctl_file" format="txt" type="data" label="PAML codeml control file" />
|
|
|
83 </when>
|
|
|
84 </conditional>
|
|
|
85 <conditional name="fit_components_cond">
|
|
|
86 <param name="fit_components" type="select" label="Fit a mixture model of multivariate normal components to synonymous (ks) distribution?" help="Used to identify significant duplication events in a genome">
|
|
|
87 <option value="no" selected="true">No</option>
|
|
|
88 <option value="yes">Yes</option>
|
|
|
89 </param>
|
|
|
90 <when value="no" />
|
|
|
91 <when value="yes">
|
|
|
92 <param name="num_of_components" type="integer" value="0" min="0" label="Number components to fit to synonymous subsitutions (ks) distribution" />
|
|
|
93 <param name="min_ks" type="float" value="0.0" min="0.0" label="Lower limit of synonymous subsitutions (ks)" help="Reduces background noise from young paralogous pairs due to normal gene births and deaths in a genome" />
|
|
|
94 <param name="max_ks" type="float" value="0.0" min="0.0" label="Upper limit of synonymous subsitutions (ks)" help="Excludes likey ancient paralogous pairs" />
|
|
|
95 </when>
|
|
|
96 </conditional>
|
|
|
97 </when>
|
|
|
98 </conditional>
|
|
|
99 </inputs>
|
|
|
100 <outputs>
|
|
|
101 <collection name="ks_files" type="list">
|
|
|
102 <discover_datasets pattern="__name__" directory="kaksAnalysis_dir" visible="false" ext="fasta" />
|
|
|
103 </collection>
|
|
|
104 </outputs>
|
|
|
105 <tests>
|
|
|
106 <test>
|
|
|
107 </test>
|
|
|
108 </tests>
|
|
|
109 <help>
|
|
|
110 This tool is one of the PlantTribes' collection of automated modular analysis pipelines that utilize objective classifications of
|
|
1
|
111 complete protein sequences from sequenced plant genomes to perform comparative evolutionary studies. This tool performs orthologous
|
|
|
112 or paralogous ks analyses of coding sequences and amino acid sequences.
|
|
0
|
113
|
|
|
114 -----
|
|
|
115
|
|
|
116 **Options**
|
|
|
117
|
|
|
118 * **Required options**
|
|
|
119
|
|
|
120 - **Select gene family clusters** - Sequences classified into gene family clusters, optionally including corresponding coding sequences.
|
|
|
121 - **Orthogroups or gene families proteins scaffold** - PlantTribes scaffolds data.
|
|
|
122 - **Protein clustering method** - One of GFam (domain architecture based clustering), OrthoFinder (broadly defined clusters) or OrthoMCL (narrowly defined clusters).
|
|
|
123
|
|
|
124 * **Multiple sequence alignments options**
|
|
|
125
|
|
|
126 - **Select method for multiple sequence alignments** - Method used for setting multiple sequence alignments.
|
|
|
127 - **Input sequences include corresponding coding sequences?** - Selecting 'Yes' for this option requires that the selected input data format is 'ptorthocs'.
|
|
|
128 - **Construct orthogroup multiple codon alignments?** - Construct orthogroup multiple codon alignments.
|
|
|
129 - **Sequence type used in the phylogenetic inference** - Sequence type (dna or amino acid) used in the phylogenetic inference.
|
|
|
130 - **Use corresponding coding sequences?** - Selecting 'Yes' for this option requires that the selected input data format is 'ptorthocs' or this tool will produce an error.
|
|
|
131
|
|
|
132 * **Phylogenetic trees options**
|
|
|
133
|
|
|
134 - **Phylogenetic trees inference method** - Phylogenetic trees inference method.
|
|
|
135 - **Select rooting order configuration for rooting trees??** - If 'No' is selected, trees will be rooted using the most distant taxon present in the orthogroup.
|
|
|
136 - **Number of replicates for rapid bootstrap analysis and search for the best-scoring ML tree** - Number of replicates for rapid bootstrap analysis and search for the best-scoring ML tree.
|
|
|
137 - **Maximum number of sequences in orthogroup alignments** - Maximum number of sequences in orthogroup alignments.
|
|
|
138 - **Minimum number of sequences in orthogroup alignments** - Minimum number of sequences in orthogroup alignments.
|
|
|
139
|
|
|
140 * **MSA quality control options**
|
|
|
141
|
|
|
142 - **Remove sequences with gaps of** - Removes gappy sequences in alignments (i.e., 0.5 removes sequences with 50% gaps).
|
|
|
143 - **Select process used for gap trimming** - Either nucleotide based trimming or alignments are trimed using using trimAl's ML heuristic trimming approach.
|
|
|
144 - **Remove sites in alignments with gaps of** - If the process used for gap trimming is nucleotide based, this is the gap value used when removing gappy sites in alignments (i.e., 0.1 removes sites with 90% gaps).
|
|
|
145
|
|
|
146 </help>
|
|
|
147 <citations>
|
|
|
148 <citation type="bibtex">
|
|
|
149 @unpublished{None,
|
|
|
150 author = {Eric Wafula},
|
|
|
151 title = {None},
|
|
|
152 year = {None},
|
|
|
153 url = {https://github.com/dePamphilis/PlantTribes}
|
|
|
154 }
|
|
|
155 </citation>
|
|
|
156 <citation type="doi">10.1186/1471-2105-10-421</citation>
|
|
|
157 <citation type="doi">10.1093/molbev/msm088</citation>
|
|
|
158 <citation type="doi">10.18637/jss.v004.i02</citation>
|
|
|
159 </citations>
|
|
|
160 </tool>
|
