|
9
|
1 <tool id="phyml-sms" name="PhyML-SMS" version="1.3.1">
|
|
|
2 <description>Maximum likelihood-based inference of phylogenetic trees with Smart Model Selection</description>
|
|
12
|
3 <macros>
|
|
|
4 <import>macro_seqtypedetect.xml</import>
|
|
|
5 </macros>
|
|
9
|
6 <requirements>
|
|
|
7 <requirement type="package" version="1.3">phyml-sms</requirement>
|
|
|
8 </requirements>
|
|
|
9 <command><![CDATA[
|
|
12
|
10
|
|
|
11
|
|
|
12
|
|
|
13 ## Automatic sequence detection
|
|
|
14 ## read an info file to choose which option set
|
|
|
15 #set $info = open( str($input_info) ).read()
|
|
|
16
|
|
9
|
17 sms.sh
|
|
|
18 -i $input
|
|
|
19 -o .
|
|
|
20 -p $output_models
|
|
|
21 -t
|
|
12
|
22
|
|
|
23 #if $sequence.seqtype == "cfg"
|
|
|
24 ## Automatic sequence detection
|
|
|
25 ## read an info file to choose which option set
|
|
|
26 #set $info = open( str($input_info) ).read()
|
|
|
27
|
|
|
28 #if 'dna' in $info:
|
|
|
29 -d nt
|
|
|
30 #else if 'protein' in $info :
|
|
|
31 -d aa
|
|
|
32 #else
|
|
|
33 -d $sequence.seqtype
|
|
9
|
34
|
|
|
35 -c $stat_crit
|
|
|
36
|
|
|
37 #if $inpuTree.inputtree == "true" :
|
|
|
38 -u $inpuTree.userInpuTree
|
|
|
39 #end if
|
|
|
40 -s $move
|
|
|
41
|
|
|
42 #if $support_condition.support == "sh":
|
|
|
43 -b -4
|
|
|
44 #else if $support_condition.support == "aBayes":
|
|
|
45 -b -5
|
|
|
46 #else if $support_condition.support == "no":
|
|
|
47 -b 0
|
|
|
48 #else if $support_condition.support == "boot":
|
|
|
49 -b $support_condition.boot_number
|
|
|
50 #end if
|
|
|
51
|
|
|
52 #if $randstart.value != 0 and $move.value == "SPR" :
|
|
|
53 -r $randstart
|
|
|
54 #end if
|
|
|
55 > $output_stdout
|
|
|
56 ;
|
|
|
57 mv *_phyml_tree.txt $output_tree;
|
|
|
58 mv *_phyml_stats.txt $output_stats;
|
|
|
59 ]]>
|
|
|
60 </command>
|
|
|
61 <inputs>
|
|
|
62 <param format="phylip" name="input" type="data" label="Alignment file" help="phylip format"/>
|
|
|
63 <conditional name="sequence">
|
|
|
64 <param name="seqtype" type="select" label="Data type">
|
|
|
65 <option value="nt">Nucleic acids</option>
|
|
|
66 <option value="aa">Amino acids</option>
|
|
12
|
67 <when value="cfg">
|
|
|
68 <param name="input_info" type="data" format="txt" label="info" help="Precompute file containning sequence description (dna or protein)" />
|
|
|
69 </when>
|
|
9
|
70 </param>
|
|
|
71 <when value="nt" />
|
|
|
72 <when value="aa" />
|
|
12
|
73 </conditional>
|
|
9
|
74 <param name="stat_crit" type="select" label="Statistical criterion to select the model">
|
|
|
75 <option value="aic">AIC</option>
|
|
|
76 <option value="bic">BIC</option>
|
|
|
77 </param>
|
|
|
78 <param name="move" type="select" label="Tree topology search" display="radio">
|
|
|
79 <option value="NNI">NNI (Nearest Neighbor Interchange)</option>
|
|
|
80 <option value="SPR">SPR (Subtree Pruning and Regraphing)</option>
|
|
|
81 </param>
|
|
|
82 <conditional name="support_condition">
|
|
|
83 <param type="select" name="support" label="Branch support" help="Use aLRT or aBayes to save computing time">
|
|
|
84 <option value="sh">SH-like aLRT</option>
|
|
|
85 <option value="aBayes">aBayes</option>
|
|
|
86 <option value="boot">Bootstrap</option>
|
|
|
87 <option value="no">No branch support</option>
|
|
|
88 </param>
|
|
|
89 <when value="sh"/>
|
|
|
90 <when value="aBayes"/>
|
|
|
91 <when value="boot">
|
|
|
92 <param type="integer" name="boot_number" min="1" value="100" label="Number of bootstrap replicates" help="Must be a positive integer"/>
|
|
|
93 </when>
|
|
|
94 <when value="no"/>
|
|
|
95 </conditional>
|
|
|
96 <conditional name="inpuTree">
|
|
|
97 <param name="inputtree" type="boolean" truevalue="true" falsevalue="false" checked="False" label="Use input starting tree" />
|
|
|
98 <when value="true">
|
|
|
99 <param name="userInpuTree" type="data" label="Tree file" help="newick format"/>
|
|
|
100 </when>
|
|
|
101 <when value="false"/>
|
|
|
102 </conditional>
|
|
|
103 <param type="integer" name="randstart" value="0" min="0" max="10" label="Number of random starting trees" />
|
|
|
104 </inputs>
|
|
|
105 <outputs>
|
|
|
106 <data format="nhx" name="output_tree" label="PhyML Newick tree"/>
|
|
|
107 <data format="txt" name="output_stats" label="PhyML Statistics"/>
|
|
|
108 <data format="txt" name="output_stdout" label="SMS Best Model"/>
|
|
|
109 <data format="txt" name="output_models" label="SMS compare models"/>
|
|
|
110 </outputs>
|
|
|
111 <tests>
|
|
|
112 <test>
|
|
|
113 <param name="input" value="phylip" />
|
|
|
114 <param name="seqtype" value="nt" />
|
|
|
115 <output name="output_models" file="sms.csv" />
|
|
|
116
|
|
|
117 </test>
|
|
|
118 </tests>
|
|
|
119 <help><![CDATA[
|
|
|
120
|
|
|
121 .. class:: infomark
|
|
|
122
|
|
|
123 This script runs SMS to select the substitution model which best fits the input data.
|
|
|
124 It may also run PhyML with the selected model.
|
|
12
|
125 SMS options :
|
|
9
|
126 -h = Help
|
|
|
127 -i = **Mandatory** Input alignment file in PHYLIP format
|
|
|
128 -d = **Mandatory** Data type : 'aa' or 'nt'
|
|
|
129 -o = Output directory
|
|
|
130 -p = Output CSV filename
|
|
|
131 -c = Statistical criterion to select the model : 'AIC' (default) or 'BIC'
|
|
|
132 -u = Input starting tree (Newick format)
|
|
|
133 -t = Use this option if you want SMS to infer a tree with PhyML using the SMS selected model
|
|
12
|
134 PhyML options :
|
|
9
|
135 -s = Type of tree improvement : 'NNI (default)' or 'SPR'
|
|
|
136 -r = Number of random starting trees : 0 (default)
|
|
|
137 -b = Branch Support : >0 for bootstraps, -4 for aLRT, 0 (default)
|
|
|
138
|
|
|
139 **PhyML 20120412**
|
|
|
140
|
|
|
141 -----
|
|
|
142
|
|
|
143
|
|
|
144 ===========
|
|
|
145 Overview:
|
|
|
146 ===========
|
|
|
147
|
|
|
148 PhyML is a phylogeny software based on the maximum-likelihood principle. Early PhyML versions used a fast algorithm to perform Nearest Neighbor Interchanges (NNIs), in order to improve a reasonable starting tree topology. Since the original publication (Guindon and Gascuel 2003), PhyML has been widely used due to its simplicity and a fair accuracy/speed compromise. In the mean time research around PhyML has continued.
|
|
|
149
|
|
|
150 We designed an efficient algorithm to search the tree space using Subtree Pruning and Regrafting (SPR) topological moves (Hordijk and Gascuel 2005), and proposed a fast branch test based on an approximate likelihood ratio test (Anisimova and Gascuel 2006). However, these novelties were not included in the official version of PhyML, and we found that improvements were still needed in order to make them effective in some practical cases. PhyML 3.0 achieves this task.
|
|
|
151
|
|
|
152 It implements new algorithms to search the space of tree topologies with user-defined intensity. A non-parametric, Shimodaira-Hasegawa-like branch test is also available. The program provides a number of new evolutionary models and its interface was entirely re-designed. We tested PhyML 3.0 on a large collection of real data sets to ensure that the new version is stable, ready-to-use and still reasonably fast and accurate.
|
|
|
153
|
|
|
154 -----
|
|
|
155
|
|
|
156 For further informations, please visite the PhyML_ and SMS_ website.
|
|
|
157
|
|
|
158
|
|
|
159 .. _PhyML: http://www.atgc-montpellier.fr/phyml/
|
|
|
160 .. _SMS: http://www.atgc-montpellier.fr/phyml-sms/
|
|
|
161
|
|
|
162
|
|
|
163 Models References:
|
|
|
164
|
|
|
165 - Dayhoff : Dayhoff, M., Schwartz, R. & Orcutt, B.
|
|
|
166 A model of evolutionary change in proteins.
|
|
|
167 In Dayhoff, M. (ed.) Atlas of Protein Sequence and Structure, vol. 5, 345–352 (National Biomedical Research Foundation, Washington, D. C., 1978)
|
|
|
168
|
|
|
169 - JTT : Jones, D., Taylor, W. & Thornton, J.
|
|
|
170 The rapid generation of mutation data matrices from protein sequences.
|
|
|
171 Computer Applications in the Biosciences (CABIOS) 8, 275–282 (1992).
|
|
|
172
|
|
|
173 - Blosum62 : Henikoff, S. & Henikoff, J.
|
|
|
174 Amino acid substitution matrices from protein blocks.
|
|
|
175 Proceedings of the National Academy of Sciences of the United States of America (PNAS) 89, 10915–10919 (1992).
|
|
|
176
|
|
|
177 - MtREV : Adachi, J. & Hasegawa, M.
|
|
|
178 MOLPHY version 2.3. programs for molecular phylogenetics based on maximum likelihood.
|
|
|
179 In Ishiguro, M. et al. (eds.) Computer Science Monographs, 28 (The Institute of Statistical Mathematics, Tokyo,1996).
|
|
|
180
|
|
|
181 - MtMam : Cao, Y., A. Janke , P. J. Waddell, M. Westerman, O. Takenaka, S. Murata, N. Okada, S. Paabo, and M. Hasegawa
|
|
|
182 Conflict among individual mitochondrial proteins in resolving the phylogeny of eutherian orders.
|
|
|
183 Journal of Molecular Evolution 47, 307–322 (1998).
|
|
|
184
|
|
|
185 - CpREV : Adachi, J., P., W., Martin, W. & Hasegawa, M.
|
|
|
186 Plastid genome phylogeny and a model of amino acid substitution for proteins encoded by chloroplast DNA.
|
|
|
187 Journal of Molecular Evolution 50, 348–358 (2000).
|
|
|
188
|
|
|
189 - VT : Muller, T. & Vingron, M.
|
|
|
190 Modeling amino acid replacement.
|
|
|
191 Journal of Computational Biology 7, 761–776 (2000).
|
|
|
192
|
|
|
193 - WAG : Whelan, S. & Goldman, N.
|
|
|
194 A general empirical model of protein evolution derived from multiple protein families using a maximum-likelihood approach.
|
|
|
195 Molecular Biology and Evolution 18, 691–699 (2001).
|
|
|
196
|
|
|
197 - RtREV : Dimmic, M., Rest, J., Mindell, D. & Goldstein, D.
|
|
|
198 rtREV : an amino acid substitution matrix for inference of retrovirus and reverse transcriptase phylogeny.
|
|
|
199 Journal of Molecular Evolution 55, 65–73 (2002).
|
|
|
200
|
|
|
201 - DCMut : Kosiol, C. & Goldman, N.
|
|
|
202 Different versions of the Dayhoff rate matrix.
|
|
|
203 Molecular Biology and Evolution 22, 193–199 (2004).
|
|
|
204
|
|
|
205 - MtArt : Abascal F, Posada D, Zardoya R.
|
|
|
206 MtArt: a new model of amino acid replacement for Arthropoda.
|
|
|
207 Mol Biol Evol. 2007 Jan;24(1):1-5. Epub 2006 Oct 16.
|
|
|
208
|
|
|
209 - HIVb - HIVw : Nickle DC, Heath L, Jensen MA, Gilbert PB, Mullins JI, Kosakovsky Pond SL.
|
|
|
210 HIV-Specific Probabilistic Models of Protein Evolution.
|
|
|
211 PLoS ONE. 2007 Jun 6;2:e503.
|
|
|
212
|
|
|
213 - LG : Le, S. & Gascuel, O.
|
|
|
214 An improved general amino-acid replacement matrix.
|
|
|
215 Mol. Biol. Evol. 25, 1307–1320 (2008)
|
|
|
216
|
|
|
217 - MtZOA : Rota-Stabelli O, Yang Z, Telford MJ
|
|
|
218 MtZoa: A general mitochondrial amino acid substitutions model for animal evolutionary studies
|
|
|
219 Molecular Phylogenetics and Evolution 52 (2009) 268–272
|
|
|
220
|
|
|
221 - FLU : Cuong Cao Dang, Quang Si Le2, Olivier Gascuel and Vinh Sy Le
|
|
|
222 FLU, an amino acid substitution model for influenza proteins.
|
|
|
223 BMC Evolutionary Biology 2010 Apr 12;10:99
|
|
|
224
|
|
|
225 - AB :
|
|
|
226 Alexander Mirsky,Linda Kazandjian and Maria Anisimova
|
|
|
227 Antibody-specific Model of Amino Acid Substitution for Immunological Inferences from Alignments of Antibody Sequences
|
|
|
228 Mol Biol Evol (2014) doi: 10.1093/molbev/msu340
|
|
|
229 ]]>
|
|
|
230 </help>
|
|
|
231 <citations>
|
|
|
232 <citation type="doi">10.1093/sysbio/syq010</citation>
|
|
|
233 <citation type="doi">10.1093/nar/gki352</citation>
|
|
|
234 <citation type="doi">10.1080/10635150390235520</citation>
|
|
|
235 </citations>
|
|
|
236 </tool>
|
