Maximum likelihood phylogenetic estimation from DNA sequences with variable rates over sites: Approximate methods

Abstract

Abstract. Two approximate methods are proposedfor maximum likelihood phylogenetic estimation, which allow variable rates of substitution across nucleotide sites. Three data sets with quite different characteristics were analyzed to examine empirically the performance of these methods. The first, called the "discrete gamma model," uses several categories of rates to approximate the gamma distribution, with equal probability for each category. The mean of each category is used to represent all the rates falling in the category. The performance of this method is found to be quite good, and four such categories appear to be sufficient to produce both an optimum, or near-optimum fit by the model to the data, and also an acceptable approximation to the continuous distribution. The second method, called "fixed-rates model," classifies sites into several classes according to their rates predicted assuming the star tree. Sites in different classes are then assumed to be evolving at these fixed rates when other tree topologies are evaluated. Analyses of the data sets suggest that this method can produce reasonable results, but it seems to share some properties of a least-squares pairwise comparison; for example, interior branch lengths in nonbest trees are often found to be zero. The computational requirements of the two methods are comparable to that of Felsenstein's (1981, J Mol Evol 17:368-376) model, which assumes a single rate for all the sites.

Genomic adaptation to polyphagy and insecticides in a major East Asian noctuid pest

“…When the alignment showed highly conservative and non-conservative regions, only the conservative regions were retained for further analysis. Model selection was conducted by MEGA version 6 75 and the LG+ Gamma+ I mode [76][77][78] . The maximum likelihood tree was inferred by RaxML version 8 79 using the LG+ Gamma+ I model.…”

Estimating Maximum Likelihood Phylogenies with PhyML

“…PhyML v3.0 also implements Ziheng Yang's discrete gamma model (27) to describe the variability of substitution rates across nucleotide or amino acids positions. Users can specify the number of substitution rate categories (Substitution Model sub-menu option 'C' or '-c' argument from the command line) and choose to estimate the gamma shape parameter from the data or fix its value a priori (option 'A' or '-a').…”

Age and biogeography of major clades in sturgeons and paddlefishes (Pisces: Acipenseriformes)

“…This parametric approach relaxes the assumption of a strict Table 2 Twenty-four mitochondrial cytochrome b gene sequences used for dating the divergence time of acipenseriform species retrieved from GenBank for the second dating analysis a Nine distribution areas were coded from A to I according to the analysis of Bemis and Kynard (1997) (Yang, 1994) along with their variance-covariance matrix. Parameters for the F84 + G model were estimated using the BASEML program in PAML (Yang, 2000).…”