PHYML User's guide (commandline interface)


Overview
PHYML is a software implementing a new method for building phylogenies
from DNA and protein sequences using maximum likelihood. Data sets can
be analyzed under several models of evolution (JC69, K80, F81, F84,
HKY85, TN93 and GTR for nucleotides and Dayhoff, JTT, mtREV, WAG,
DCMut, RtREV, CpREV, VT, Blosum62 and MtMam for amino acids). A discretegamma model (Yang, 1994) is implemented to accommodate rate variation among sites. Invariable sites can also be taken into account. PHYML has been compared to several other softwares using extensive simulations. The results indicate that its topological accuracy is at least as high as that of fastDNAml, while being much faster.


The commandline interfaceDownload the binary files ; you can execute PHYML by typing "./phyml" followd by a list of parameters. Type 13 parameters for DNA sequences :./phyml sequences file data type sequence format nb data sets nb bootstrapped data sets substitution model ts/tv ratio prop. invariable sites nb categories gamma parameter starting tree optimise topology optimise branch lengths and rate parameters Example : ./phyml seqs1 0 i 2 0 HKY 4.0 e 1 1.0 BIONJ y y Type 12 parameters for aminoacids sequences : ./phyml sequences file data type sequence format nb data sets nb bootstrapped data sets substitution model prop. invariable sites nb categories gamma parameter starting tree optimise topology optimise branch lengths and rate parameters Example : ./phyml seqs2 1 i 1 0 JTT 0.0 4 1.0 BIONJ n nFor complete details type './phyml h' or see the 'commandline' specific comments in the 'Options' section below. PHYML enables to analyze one or several data sets in conjunction with one or several starting trees. PHYML produces several results files : Here are the possible uses of PHYML : Standard analysis under a given substitution model, PHYML then returns the inferred tree. Moreover, a special option allows to perform nonparametric bootstrapp analysis on the original data set. PHYML then returns the bootstrap tree with branch lengths and bootstrap values, using standard NEWICK format (an option gives the pseudo trees in a *_boot_trees.txt file). Several standard analysis start from the same intial tree with different data sets, without the bootstrap option. The results are given in the order of the data sets. This can be used to process multiple genes in a supertree approach. Several standard analysis of the same data set using different starting tree situations, without the bootstrap option. All results are given in the order of the trees. Moreover, the most likely tree is provided in the *_best_stat.txt and *_best_tree.txt files. This should be used to avoid being trapped into local optima and then obtain better trees. Fast parsimony methods can be used to obtain a set of starting trees. Several standard runs, where each data set is analysed with the corresponding starting tree, without the bootstrap option. The results are given in the order of the data sets. This can be used when comparing the likelihood of various trees regarding different data sets. 

OptionsThe input sequence file is a standard PHYLIP file of aligned DNA or aminoacids sequences. It should look like this in interleaved format : 5 60 Tax1 CCATCTCACGGTCGGTACGATACACCTGCTTTTGGCAG Tax2 CCATCTCACGGTCAGTAAGATACACCTGCTTTTGGCGG Tax3 CCATCTCCCGCTCAGTAAGATACCCCTGCTGTTGGCGG Tax4 TCATCTCATGGTCAATAAGATACTCCTGCTTTTGGCGG Tax5 CCATCTCACGGTCGGTAAGATACACCTGCTTTTGGCGG GAAATGGTCAATATTACAAGGT GAAATGGTCAACATTAAAAGAT GAAATCGTCAATATTAAAAGGT GAAATGGTCAATCTTAAAAGGT GAAATGGTCAATATTAAAAGGTThe same data set in sequential format: 5 60 Tax1 CCATCTCACGGTCGGTACGATACACCTGCTTTTGGCAGGAAATGGTCAATATTACAAGGT Tax2 CCATCTCACGGTCAGTAAGATACACCTGCTTTTGGCGGGAAATGGTCAACATTAAAAGAT Tax3 CCATCTCCCGCTCAGTAAGATACCCCTGCTGTTGGCGGGAAATCGTCAATATTAAAAGGT Tax4 TCATCTCATGGTCAATAAGATACTCCTGCTTTTGGCGGGAAATGGTCAATCTTAAAAGGT Tax5 CCATCTCACGGTCGGTAAGATACACCTGCTTTTGGCGGGAAATGGTCAATATTAAAAGGT On the first line is the number of taxa, a space, then the number of characters for each taxon. The maximum number of characters in species name MUST not exceed 50. Blanks within the species name are NOT allowed. However, blanks (one or more) MUST appear at the end of each species name. In a sequence, three special characters '.', '', and '?' may be used: a dot '.' means the same character as in the first sequence, a dash '' means an alignment gap and a question mark '?' means an undetermined nucleotide. Sites at which one or more sequences involve '' are NOT excluded from the analysis. Therefore, gaps are treated as unknown character (like '?') on the grounds that ''we don't know what would be there if something were there'' (J. Felsenstein, PHYLIP documentation). Finally, standard ambiguity characters for nucleotides are accepted (Table 1).
commandline : type the sequence filename with its path from the current directory This indicates if the sequence file contains DNA or aminoacids. The default choice is to analyze DNA sequences. commandline : the data type is specified by the number of parameters, 9 for DNA, 8 for aminoacids The input sequences can be either in interleaved (default) or sequential format, see "Sequences" above. commandline : type 'i' or 's' 5 60 Tax1 CCATCTCACGGTCGGTACGATACACCTGCTTTTGGCAGGAAATGGTCAATATTACAAGGT Tax2 CCATCTCACGGTCAGTAAGATACACCTGCTTTTGGCGGGAAATGGTCAACATTAAAAGAT Tax3 CCATCTCCCGCTCAGTAAGATACCCCTGCTGTTGGCGGGAAATCGTCAATATTAAAAGGT Tax4 TCATCTCATGGTCAATAAGATACTCCTGCTTTTGGCGGGAAATGGTCAATCTTAAAAGGT Tax5 CCATCTCACGGTCGGTAAGATACACCTGCTTTTGGCGGGAAATGGTCAATATTAAAAGGT 5 60 Tax1 CCATCTCACGGTCGGTACGATACACCTGCTTTTGGCAGGAAATGGTCAATATTACAAGGT Tax2 CCATCTCACGGTCAGTAAGATACACCTGCTTTTGGCGGGAAATGGTCAACATTAAAAGAT Tax3 CCATCTCCCGCTCAGTAAGATACCCCTGCTGTTGGCGGGAAATCGTCAATATTAAAAGGT Tax4 TCATCTCATGGTCAATAAGATACTCCTGCTTTTGGCGGGAAATGGTCAATCTTAAAAGGT Tax5 CCATCTCACGGTCGGTAAGATACACCTGCTTTTGGCGGGAAATGGTCAATATTAAAAGGT 5 60 Tax1 CCATCTCACGGTCGGTACGATACACCTGCTTTTGGCAGGAAATGGTCAATATTACAAGGT Tax2 CCATCTCACGGTCAGTAAGATACACCTGCTTTTGGCGGGAAATGGTCAACATTAAAAGAT Tax3 CCATCTCCCGCTCAGTAAGATACCCCTGCTGTTGGCGGGAAATCGTCAATATTAAAAGGT Tax4 TCATCTCATGGTCAATAAGATACTCCTGCTTTTGGCGGGAAATGGTCAATCTTAAAAGGT Tax5 CCATCTCACGGTCGGTAAGATACACCTGCTTTTGGCGGGAAATGGTCAATATTAAAAGGT commandline : type the value of this parameter When there is only one data sets you can ask PHYML to generate bootstrapped pseudo data sets from this original data set. PHYML returns the bootstrap tree with branchlengths and bootstrap values, using standard NEWICK format. The pseudo trees are given in a *_boot_trees.txt file. commandline : type the value of this parameter A nucleotide or aminoacid substitution model. For DNA sequences, the default choice is HKY85 (Hasegawa et al., 1985). This model is analogous to K80 (Kimura, 1980), but allows for different base frequencies. The other models are JC69 (Jukes and Cantor, 1969), K80 (Kimura, 1980), F81 (Felsenstein, 1981), F84 (Felsenstein, 1989), TN93 (Tamura and Nei, 1993) and GTR (e.g., Lanave et al. 1984, Tavaré 1986, Rodriguez et al. 1990). The rate matrices of these models are given in Swofford et al. (1996). For AminoAcid sequences, the default choice is JTT (Jones, Taylor and Thornton, 1992). The other models are Dayhoff (Dayhoff et al., 1978), mtREV (as implemented in Yang's PAML), WAG (Whelan and Goldman, 2001), DCMut (Kosiol and Goldman, 2005), RtREV (Dimmic et al.), CpREV (Adachi et al., 2000) VT (Muller and Vingron, 2000), Blosum62 (Henikoff anf Henikoff, 1992) and MtMam (Cao, 1998). commandline : type the name of the model With DNA sequences, it is possible to set the transition/transversion ratio, except for the JC69 and F81 models, or to estimate its value by maximizing the likelihood of the phylogeny. The later makes the program slower. The default value is 4.0. The definition of the transition/transversion ratio is the same as in PAML (Yang, 1994). In PHYLIP, the ''transition/transversion rate ratio'' is used instead. 4.0 in PHYML roughly corresponds to 2.0 in PHYLIP. commandline : type the value of this parameter or type 'e' to estimate it The default is to consider that the data set does not contain invariable sites (0.0). However, this proportion can be set to any value in the 0.01.0 range. This parameter can also be estimated by maximizing the likelihood of the phylogeny. The later makes the program slower. commandline : type the value of this parameter or type 'e' to estimate it The default is having all the sites evolving at the same rate, hence having one substitution rate category. A discretegamma distribution can be used to account for variable substitution rates among sites, in which case the number of categories that defines this distribution is supplied by the user. The higher this number, the better is the goodnessoffit regarding the continuous distribution. The default is to use four categories, in this case the likelihood of the phylogeny at one site is averaged over four conditional likelihoods corresponding to four rates and the computation of the likelihood is four times slower than with a unique rate. Number of categories less than four or higher than eight are not recommended. In the first case, the discrete distribution is a poor approximation of the continuous one. In the second case, the computational burden becomes high and an higher number of categories is not likely to enhance the accuracy of phylogeny estimation. commandline : type the value of this parameter The shape of a gamma distribution is defined by this numerical parameter. The higher its value, the lower the variation of substitution rates among sites (this option is used when having more than 1 substitution rate category). The default value is 1.0. It corresponds to a moderate variation. Values less than say 0.7 correspond to high variations. Values between 0.7 and 1.5 corresponds to moderate variations. Higher values correspond to low variations. This value can be fixed by the user. It can also be estimated by maximizing the likelihood of the phylogeny. commandline : type the value of this parameter or type 'e' to estimate it Used as the starting tree(s) to be refined by the maximum likelihood algorithm. The default is to use a BIONJ distancebased tree. It is also possible to supply one or several trees in NEWICK format, one per line in the file, which must be written in the standard parenthesis representation (NEWICK format) ; the branch lengths must be given, and the tree(s) must be unrooted. Labels on branches (such as bootstrap proportions) are supported. Therefore, a tree with four taxa named A, B, C, and D with a bootstrap value equals to 90 on its internal branch, should look like this: (A:0.02,B:0.004,(C:0.1,D:0.04)90:0.05); If you give several trees and analyse several data sets the two numbers must match. commandline : type the tree filename with its path from the current directory or type 'BIONJ' You can optimise the starting tree(s) in three ways :  You can optimise the topology, the branch lengths and rate parameters (transition/transversion ratio, proportion of invariant sites, gamma distribution parameter), commandline : type 'y' and 'y'  You can keep the topology and optimise the branch lengths and rate parameters (it is not possible to optimise the tree topology and keep the branch lengths and rate parameters), commandline : type 'n' and 'y'  You can ask for no optimisation, PHYML just computes the likelihood of the starting tree(s). commandline : type 'n' and 'n' 

References 
