Comparative analysis of molecular sequence data is essential for reconstructing the evolutionary histories of species and inferring the nature and extent of selective forces shaping the evolution of genes and species. Here, we announce the release of Molecular Evolutionary Genetics Analysis version 5 (MEGA5), which is a user-friendly software for mining online databases, building sequence alignments and phylogenetic trees, and using methods of evolutionary bioinformatics in basic biology, biomedicine, and evolution. The newest addition in MEGA5 is a collection of maximum likelihood (ML) analyses for inferring evolutionary trees, selecting best-fit substitution models (nucleotide or amino acid), inferring ancestral states and sequences (along with probabilities), and estimating evolutionary rates site-by-site. In computer simulation analyses, ML tree inference algorithms in MEGA5 compared favorably with other software packages in terms of computational efficiency and the accuracy of the estimates of phylogenetic trees, substitution parameters, and rate variation among sites. The MEGA user interface has now been enhanced to be activity driven to make it easier for the use of both beginners and experienced scientists. This version of MEGA is intended for the Windows platform, and it has been configured for effective use on Mac OS X and Linux desktops. It is available free of charge from http://www.megasoftware.net.

/pdf/mega5-molecular-evolutionary-genetics-analysis-using-maximum-2fhrdnt1zl.pdf

MEGA5: Molecular Evolutionary Genetics Analysis using Maximum Likelihood, Evolutionary Distance, and Maximum Parsimony Methods

Summary: The program MODELTEST uses log likelihood scores to establish the model of DNA evolution that best fits the data. Availability: The MODELTEST package, including the source code and some documentation is available at http://bioag.byu.edu/zoology/crandall―lab/modeltest.html. Contact: dp47@email.byu.edu.

/pdf/modeltest-testing-the-model-of-dna-substitution-2dmippvvyz.pdf

MODELTEST: testing the model of DNA substitution.

The increase in the number of large data sets and the complexity of current probabilistic sequence evolution models necessitates fast and reliable phylogeny reconstruction methods. We describe a new approach, based on the maximum- likelihood principle, which clearly satisfies these requirements. The core of this method is a simple hill-climbing algorithm that adjusts tree topology and branch lengths simultaneously. This algorithm starts from an initial tree built by a fast distance-based method and modifies this tree to improve its likelihood at each iteration. Due to this simultaneous adjustment of the topology and branch lengths, only a few iterations are sufficient to reach an optimum. We used extensive and realistic computer simulations to show that the topological accuracy of this new method is at least as high as that of the existing maximum-likelihood programs and much higher than the performance of distance-based and parsimony approaches. The reduction of computing time is dramatic in comparison with other maximum-likelihood packages, while the likelihood maximization ability tends to be higher. For example, only 12 min were required on a standard personal computer to analyze a data set consisting of 500 rbcL sequences with 1,428 base pairs from plant plastids, thus reaching a speed of the same order as some popular distance-based and parsimony algorithms. This new method is implemented in the PHYML program, which is freely available on our web page: http://www.lirmm.fr/w3ifa/MAAS/. (Algorithm; computer simulations; maximum likelihood; phylogeny; rbcL; RDPII project.) The size of homologous sequence data sets has in- creased dramatically in recent years, and many of these data sets now involve several hundreds of taxa. More- over, current probabilistic sequence evolution models (Swofford et al., 1996 ; Page and Holmes, 1998 ), notably those including rate variation among sites (Uzzell and Corbin, 1971 ; Jin and Nei, 1990 ; Yang, 1996 ), require an increasing number of calculations. Therefore, the speed of phylogeny reconstruction methods is becoming a sig- nificant requirement and good compromises between speed and accuracy must be found. The maximum likelihood (ML) approach is especially accurate for building molecular phylogenies. Felsenstein (1981) brought this framework to nucleotide-based phy- logenetic inference, and it was later also applied to amino acid sequences (Kishino et al., 1990). Several vari- ants were proposed, most notably the Bayesian meth- ods (Rannala and Yang 1996; and see below), and the discrete Fourier analysis of Hendy et al. (1994), for ex- ample. Numerous computer studies (Huelsenbeck and Hillis, 1993; Kuhner and Felsenstein, 1994; Huelsenbeck, 1995; Rosenberg and Kumar, 2001; Ranwez and Gascuel, 2002) have shown that ML programs can recover the cor- rect tree from simulated data sets more frequently than other methods can. Another important advantage of the ML approach is the ability to compare different trees and evolutionary models within a statistical framework (see Whelan et al., 2001, for a review). However, like all optimality criterion-based phylogenetic reconstruction approaches, ML is hampered by computational difficul- ties, making it impossible to obtain the optimal tree with certainty from even moderate data sets (Swofford et al., 1996). Therefore, all practical methods rely on heuristics that obtain near-optimal trees in reasonable computing time. Moreover, the computation problem is especially difficult with ML, because the tree likelihood not only depends on the tree topology but also on numerical pa- rameters, including branch lengths. Even computing the optimal values of these parameters on a single tree is not an easy task, particularly because of possible local optima (Chor et al., 2000). The usual heuristic method, implemented in the pop- ular PHYLIP (Felsenstein, 1993 ) and PAUP ∗ (Swofford, 1999 ) packages, is based on hill climbing. It combines stepwise insertion of taxa in a growing tree and topolog- ical rearrangement. For each possible insertion position and rearrangement, the branch lengths of the resulting tree are optimized and the tree likelihood is computed. When the rearrangement improves the current tree or when the position insertion is the best among all pos- sible positions, the corresponding tree becomes the new current tree. Simple rearrangements are used during tree growing, namely "nearest neighbor interchanges" (see below), while more intense rearrangements can be used once all taxa have been inserted. The procedure stops when no rearrangement improves the current best tree. Despite significant decreases in computing times, no- tably in fastDNAml (Olsen et al., 1994 ), this heuristic becomes impracticable with several hundreds of taxa. This is mainly due to the two-level strategy, which sepa- rates branch lengths and tree topology optimization. In- deed, most calculations are done to optimize the branch lengths and evaluate the likelihood of trees that are finally rejected. New methods have thus been proposed. Strimmer and von Haeseler (1996) and others have assembled four- taxon (quartet) trees inferred by ML, in order to recon- struct a complete tree. However, the results of this ap- proach have not been very satisfactory to date (Ranwez and Gascuel, 2001 ). Ota and Li (2000, 2001) described

A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood.

The evolutionary analysis of molecular sequence variation is a statistical enterprise. This is reflected in the increased use of probabilistic models for phylogenetic inference, multiple sequence alignment, and molecular population genetics. Here we present BEAST: a fast, flexible software architecture for Bayesian analysis of molecular sequences related by an evolutionary tree. A large number of popular stochastic models of sequence evolution are provided and tree-based models suitable for both within- and between-species sequence data are implemented. BEAST version 1.4.6 consists of 81000 lines of Java source code, 779 classes and 81 packages. It provides models for DNA and protein sequence evolution, highly parametric coalescent analysis, relaxed clock phylogenetics, non-contemporaneous sequence data, statistical alignment and a wide range of options for prior distributions. BEAST source code is object-oriented, modular in design and freely available at http://beast-mcmc.googlecode.com/
 under the GNU LGPL license. BEAST is a powerful and flexible evolutionary analysis package for molecular sequence variation. It also provides a resource for the further development of new models and statistical methods of evolutionary analysis.

/pdf/beast-bayesian-evolutionary-analysis-by-sampling-trees-2t1t8jmg8x.pdf

BEAST: Bayesian evolutionary analysis by sampling trees

PAML, currently in version 4, is a package of programs for phylogenetic analyses of DNA and protein sequences using maximum likelihood (ML). The programs may be used to compare and test phylogenetic trees, but their main strengths lie in the rich repertoire of evolutionary models implemented, which can be used to estimate parameters in models of sequence evolution and to test interesting biological hypotheses. Uses of the programs include estimation of synonymous and nonsynonymous rates (d(N) and d(S)) between two protein-coding DNA sequences, inference of positive Darwinian selection through phylogenetic comparison of protein-coding genes, reconstruction of ancestral genes and proteins for molecular restoration studies of extinct life forms, combined analysis of heterogeneous data sets from multiple gene loci, and estimation of species divergence times incorporating uncertainties in fossil calibrations. This note discusses some of the major applications of the package, which includes example data sets to demonstrate their use. The package is written in ANSI C, and runs under Windows, Mac OSX, and UNIX systems. It is available at -- (http://abacus.gene.ucl.ac.uk/software/paml.html).

/pdf/paml-4-phylogenetic-analysis-by-maximum-likelihood-2gg19rsgos.pdf

PAML 4: Phylogenetic Analysis by Maximum Likelihood

A new statistical method for estimating divergence dates of species from DNA sequence data by a molecular clock approach is developed. This method takes into account effectively the information contained in a set of DNA sequence data. The molecular clock of mitochondrial DNA (mtDNA) was calibrated by setting the date of divergence between primates and ungulates at the Cretaceous-Tertiary boundary (65 million years ago), when the extinction of dinosaurs occurred. A generalized least-squares method was applied in fitting a model to mtDNA sequence data, and the clock gave dates of 92.3 +/- 11.7, 13.3 +/- 1.5, 10.9 +/- 1.2, 3.7 +/- 0.6, and 2.7 +/- 0.6 million years ago (where the second of each pair of numbers is the standard deviation) for the separation of mouse, gibbon, orangutan, gorilla, and chimpanzee, respectively, from the line leading to humans. Although there is some uncertainty in the clock, this dating may pose a problem for the widely believed hypothesis that the pipedal creature Australopithecus afarensis, which lived some 3.7 million years ago at Laetoli in Tanzania and at Hadar in Ethiopia, was ancestral to man and evolved after the human-ape splitting. Another likelier possibility is that mtDNA was transferred through hybridization between a proto-human and a proto-chimpanzee after the former had developed bipedalism.

Dating of the human-ape splitting by a molecular clock of mitochondrial DNA.

Phylogenetic trees among eukaryotic kingdoms were inferred for large- and small-subunit rRNAs by using a maximum-likelihood method developed by Felsenstein. Although Felsenstein's method assumes equal evolutionary rates for transitions and transversions, this is apparently not the case for these data. Therefore, only transversiontype substitutions were taken into account. The molecules used were large-subunit rRNAs fromXenopus laevis (Animalia), rice (Plantae),Saccharomyces cerevisiae (Fungi),Dictyostelium discoideum (Protista), andPhysarum polycephalum (Protista); and small-subunit rRNAs from maize (Plantae),S. cerevisiae, X. laevis, rat (Animalia), andD. discoideum. Only conservative regions of the nucleotide sequences were considered for this study. In the maximum-likelihood trees for both large- and small-subunit rRNAs, Animalia and Fungi were the most closely related eukaryotic kingdoms, and Plantae is the next most closely related kingdom, although other branching orders among Plantae, Animalia, and Fungi were not excluded by this work. These three eukaryotic kingdoms apparently shared a common ancestor after the divergence of the two species of Protista,D. discoideum andP. polycephalum. These two species of Protista do not form a clade, andP. polycephalum diverged first andD. discoideum second from the line leading to the common ancestor of Plantae, Animalia, and Fungi. The sequence data indicate that a drastic change occurred in the nucleotide sequences of rRNAs during the evolutionary separation between prokaryote and eukaryote.

Phylogenetic relationships among eukaryotic kingdoms inferred from ribosomal RNA sequences.

Divergence dates among primates were estimated by molecular clock analysis of DNA sequence data. A molecular clock of eta-globin pseudogene was calibrated by setting the date of divergence between Catarrhini and Platyrrhini at 38 million years (Myr) ago. The clock gave dates of 25.3 +/- 2.4, 11.9 +/- 1.7, 5.9 +/- 1.2, and 4.9 +/- 1.2 Myr ago ( +/- refers to standard error) for the separation of rhesus monkey, orangutan, gorilla, and chimpanzee, respectively, from the line leading to humans. In placing confidence intervals of the estimates in a robust way, a bootstrap method was used. The 95% confidence intervals are 20.5-29.5, 9.0-14.8, 4.1-7.8, and 3.1-7.0 Myr ago for the separation of rhesus monkey, orangutan, gorilla, and chimpanzee, respectively. By a molecular clock dating of the Prosimii-Anthropoidea splitting, it was suggested that the evolutionary rate of the eta-globin gene was high early in primate evolution and subsequently decreased in the line of Anthropoidea. And, by a relative rate test using bootstrap sampling, the possibility of further decrease of the rate (more than 10%) in the line of Hominoidea compared with that of Cercopithecoidea was suggested. Therefore, the above dating of the splittings within Hominoidea may be biased slightly toward younger dates. On the other hand, mitochondrial DNA (mtDNA) seems to have evolved in mammals with a more uniform rate than the eta-globin gene. The ratio of the dates of orangutan splitting to chimpanzee splitting is larger for the mtDNA clock than that for the eta-globin clock, suggesting the possibilities of mtDNA introgression among the early hominids and the early African apes, and/or of mtDNA polymorphism within the common ancestral species of orangutan and the African apes that obscures the date of the true species separation of orangutans.

Man's place in Hominoidea as inferred from molecular clocks of DNA.

Estimation of branching dates among primates by molecular clocks of nuclear DNA which slowed down in Hominoidea

The approximate constancy of the rate of change is one feature of molecular evolution. It is consistent with the neutral theory1 and serves as a clock to date a divergence time of organisms.2>,3) Since Sarich and Wilson,4~ many researchers have estimated the divergence time between man and the African apes using the molecular clock (for review, see ref. 5). In spite of diverse materials and methods used, their results are similar and show a recent divergence of 4-& million years (Myr) ago. This dating does not contradict with any fossil evidence, if Ramapithecus is placed on orangutan line.5~ '0> The molecular clock, however, is not yet free of any contradiction with hominoid fossil records. The clock estimate of the divergence time between orangutan and man is between 8 Myr4~ and 10±3 Myr ago,5~ though recent finding indicates that orangutan and man diverged at a minimum of 13 Myr, probably more than 14.5 Myr. 7 ) The clockk dated the split of gibbons from man between 10 Myr4~ to 12±3 Myr ago.5~ But the fossil evidence indicates that Epipliopithecius and Micropithecus are 16 and 20 Myr old, respectively, and are ancestors of gibbons.8~ Since the previous estimates of genetic distances are indirect, we performed a direct comparison among mitochondrial. DNA (mtDNA) sequences. Our data are the L-strand mtDNA sequences of 896 nucleotides from five species of hominoids,9~ from bovine,10~ and from mouse.'1 This segment contains genes for three tRNAs and parts of two proteins. In animal mtDNA, it is known that transition predominates over transversion,9~ and that the substitution at the third colon position is much more frequent than other nucleotide positions. The numbers of transition and transversion type differences between species are listed in Table I separately for the third colon positions (class 1 sites) and for the rest (class 2 sites). Table I indicates that mouse (splitting date : t1 Myr ago), bovine

Taka-aki Yano

Papers

Dating of the human-ape splitting by a molecular clock of mitochondrial DNA.

Phylogenetic relationships among eukaryotic kingdoms inferred from ribosomal RNA sequences.

Man's place in Hominoidea as inferred from molecular clocks of DNA.

Estimation of branching dates among primates by molecular clocks of nuclear DNA which slowed down in Hominoidea

A New Molecular Clock of Mitochondrial DNA and the Evolution of Hominoids