Molecular differentiation and phylogeny of entomopathogenic nematodes (rhabditida: heterorhabditidae) based on ND4 gene sequences of Mitochondrial DNA.

We determined partial ND4 gene sequences of mitochondrial DNA from 15 heterorhabditid nematode isolates, representing 5 species collected from different regions of the world, by using polymerase chain reaction (PCR) and direct-sequencing of PCR products. Aligned nucleotide as well as amino acid sequences were used to differentiate nematode species by comparing sequence divergence and to infer phylogeny of the nematodes by using maximum parsimony and likelihood methods. Robustness of our phylogenetic trees was checked by bootstrap tests. The 15 nematode isolates can be divided into 7 haplotypes based on DNA sequences. On a larger scale, the sequence divergence revealed 4 distinct groups corresponding to 4 described species. No sequence divergence was detected from 5 isolates of Heterorhabditis bacteriophora or between Heterorhabditis marelatus to Heterorhabditis hepialius. Our sequence data yielded phylogenetic trees with identical topologies when different tree-building methods were used. Most relationships were also confirmed by using amino acid sequences in maximum parsimony analysis. Our molecular phylogeny of Heterorhabditis species support an existing taxonomy that is based largely on morphology and the sequence divergence of the ND4 gene permits species identification.     

Evolution of genes and taxa: a primer

The rapidly growing fields of molecular evolution and systematics have much to offer to molecular biology, but like any field have their own repertoire of terms and concepts. Homology, for example, is a central theme in evolutionary biology whose definition is complex and often controversial. Homology extends to multigene families, where the distinction between orthology and paralogy is key. Nucleotide sequence alignment is also a homology issue, and is a key stage in any evolutionary analysis of sequence data. Models based on our understanding of the processes of nucleotide substitution are used both in the estimation of the number of evolutionary changes between aligned sequences and in phylogeny reconstruction from sequence data. The three common methods of phylogeny reconstruction – parsimony, distance and maximum likelihood – differ in their use of these models. All three face similar problems in finding optimal – and reliable – solutions among the vast number of possible trees. Moreover, even optimal trees for a given gene may not reflect the relationships of the organisms from which the gene was sampled. Knowledge of how genes evolve and at what rate is critical for understanding gene function across species or within gene families. The Neutral Theory of Molecular Evolution serves as the null model of molecular evolution and plays a central role in data analysis. Three areas in which the Neutral Theory plays a vital role are: interpreting ratios of nonsynonymous to synonymous nucleotide substitutions, assessing the reliability of molecular clocks, and providing a foundation for molecular population genetics.     

A New Method of Inference of Ancestral Nucleotide and Amino Acid Sequences

A statistical method was developed for reconstructing the nucleotide or amino acid sequences of extinct ancestors, given the phylogeny and sequences of the extant species. A model of nucleotide or amino acid substitution was employed to analyze data of the present-day sequences, and maximum likelihood estimates of parameters such as branch lengths were used to compare the posterior probabilities of assignments of character states (nucleotides or amino acids) to interior nodes of the tree; the assignment having the highest probability was the best reconstruction at the site. The lysozyme c sequences of six mammals were analyzed by using the likelihood and parsimony methods. The new likelihood-based method was found to be superior to the parsimony method. The probability that the amino acids for all interior nodes at a site reconstructed by the new method are correct was calculated to be 0.91, 0.86, and 0.73 for all, variable, and parsimony-informative sites, respectively, whereas the corresponding probabilities for the parsimony method were 0.84, 0.76, and 0.51, respectively. The probability that an amino acid in an ancestral sequence is correctly reconstructed by the likelihood analysis ranged from 91.3 to 98.7% for the four ancestral sequences.     

A hominoid-specific nuclear insertion of the mitochondrial D-loop: implications for reconstructing ancestral mitochondrial sequences

A nuclear integration of a mitochondrial control region sequence on human chromosome 9 has been isolated. PCR analyses with primers specific for the respective insertion-flanking nuclear regions showed that the insertion took place on the lineage leading to Hominoidea (gibbon, orangutan, gorilla, chimpanzee, and human) after the Old World monkey-Hominoidea split. The sequences of the control region integrations were determined for humans, chimpanzees, gorillas, orangutans, and siamangs. These sequences were then used to construct phylogenetic trees with different methods, relating them with several hominoid, Old Work monkey, and New World monkey mitochondrial control region sequences. Applying maximum-likelihood, neighbor-joining, and parsimony algorithms, the insertion clade was attached to the branch leading to the hominoid mitochondrial sequences as expected from the PCR-determined presence/absence of this integration. An unexpected long branch leading to the internal node that connects all insertion sequences was observed for the different phylogeny reconstruction procedures. This finding is not totally compatible with the lower evolutionary rate in the nucleus than in the mitochondrial compartment. We determined the unambiguous substitutions on the branch leading to the most recent common ancestor (MRCA) of the mitochondrial inserts according to the parsimony criterium. We propose that they are unlikely to have been caused by damage of the transposing nucleic acid and that they are probably due to a change in the evolutionary mode after the transposition.      

Log Transformation Improves Dating of Phylogenies

Phylogenetic trees inferred from sequence data often have branch lengths measured in the expected number of substitutions and therefore, do not have divergence times estimated. These trees give an incomplete view of evolutionary histories since many applications of phylogenies require time trees. Many methods have been developed to convert the inferred branch lengths from substitution unit to time unit using calibration points, but none is universally accepted as they are challenged in both scalability and accuracy under complex models. Here, we introduce a new method that formulates dating as a nonconvex optimization problem where the variance of log-transformed rate multipliers is minimized across the tree. On simulated and real data, we show that our method, wLogDate, is often more accurate than alternatives and is more robust to various model assumptions.     

Phylogeny and divergence of the pinnipeds (Carnivora: Mammalia) assessed using a multigene dataset

Background  

Phylogenetic comparative methods are often improved by complete phylogenies with meaningful branch lengths (e.g., divergence dates). This study presents a dated molecular supertree for all 34 world pinniped species derived from a weighted matrix representation with parsimony (MRP) supertree analysis of 50 gene trees, each determined under a maximum likelihood (ML) framework. Divergence times were determined by mapping the same sequence data (plus two additional genes) on to the supertree topology and calibrating the ML branch lengths against a range of fossil calibrations. We assessed the sensitivity of our supertree topology in two ways: 1) a second supertree with all mtDNA genes combined into a single source tree, and 2) likelihood-based supermatrix analyses. Divergence dates were also calculated using a Bayesian relaxed molecular clock with rate autocorrelation to test the sensitivity of our supertree results further.     

Molecular evolution and phylogeny of the atpB-rbcL spacer of chloroplast DNA in the true mosses.

The nucleotide variation of a noncoding region between the atpB and rbcL genes of the chloroplast genome was used to estimate the phylogeny of 11 species of true mosses (subclass Bryidae). The A+T rich (82.6%) spacer sequence is conserved with 48% of bases showing no variation between the ingroup and outgroup. Rooted at liverworts, Marchantia and Bazzania, the monophyly of true mosses was supported cladistically and statistically. A nonparametric Wilcoxon Signed-Ranks test Ts statistic for testing the taxonomic congruence showed no significant differences between gene trees and organism trees as well as between parsimony trees and neighbor-joining trees. The reconstructed phylogeny based on the atpB-rbcL spacer sequences indicated the validity of the division of acrocarpous and pleurocarpous mosses. The size of the chloroplast spacer in mosses fits into an evolutionary trend of increasing spacer length from liverworts through ferns to seed plants. According to the relative rate tests, the hypothesis of a molecular clock was supported in all species except for Thuidium, which evolved relatively fast. The evolutionary rate of the chloroplast DNA spacer in mosses was estimated to be (1.12 +/- 0.019) x 10(-10) nucleotides per site per year, which is close to the nonsynonymous substitution rates of the rbcL gene in the vascular plants. The constrained molecular evolution (total nucleotide substitutions, K approximately 0.0248) of the chloroplast DNA spacer is consistent with the slow evolution in morphological traits of mosses. Based on the calibrated evolutionary rate, the time of the divergence of true mosses was estimated to have been as early as 220 million years ago.     

The Pattern of Mammalian Evolution and the Relative Rate of Molecular Evolution

The rates of nucleotide substitution at four genes in four orders of eutherian mammals are compared in relative rate tests using marsupial orthologs for reference. There is no evidence of systematic variation in evolutionary rate among the orders. The sequences are used to reconstruct the phylogeny of the orders using maximum likelihood, parsimony and compatibility methods. A branching order of rodent then ungulate then primate and lagomorph is overwhelmingly indicated. The nodes of the nucleotide based cladograms are widely separated in relation to the total lengths of the branches. The assumption of a star phylogeny that underlies Kimura's test for molecular evolutionary rate variation is shown to be invalid for eutherian mammals. Excess variance in nucleotide or amino acid differences between mammalian orders, above that predicted by neutral theory is explained better by variation in divergence time than by variation in evolutionary rate.     

Cytochrome b phylogeny of North American hares and jackrabbits (Lepus, lagomorpha) and the effects of saturation in outgroup taxa

Jackrabbits and hares, members of the genus Lepus, comprise over half of the species within the family Leporidae (Lagomorpha). Despite their ecological importance, potential economic impact, and worldwide distribution, the evolution of hares and jackrabbits has been poorly studied. We provide an initial phylogenetic framework for jackrabbits and hares so that explicit hypotheses about their evolution can be developed and tested. To this end, we have collected DNA sequence data from a 702-bp region of the mitochondrial cytochrome b gene and reconstructed the evolutionary history (via parsimony, neighbor joining, and maximum likelihood) of 11 species of Lepus, focusing on North American taxa. Due to problems of saturation, induced by multiple substitutions, at synonymous coding positions between the ingroup taxa and the outgroups (Oryctolagus and Sylvilagus), both rooted and unrooted trees were examined. Variation in tree topologies generated by different reconstruction methods was observed in analyses including the outgroups, but not in the analyses of unrooted ingroup networks. Apparently, substitutional saturation hindered the analyses when outgroups were considered. The trees based on the cytochrome b data indicate that the taxonomic status of some species needs to be reassessed and that species of Lepus within North America do not form a monophyletic entity.     

A MOLECULAR PHYLOGENY OF THE HETEROKONT ALGAE BASED ON ANALYSES OF CHLOROPLAST-ENCODED rbcL SEQUENCE DATA1

Nearly complete ribulose-1,5-bisphosphate carboxylase/ oxygenase (rbcL)sequences from 27 taxa of heterokont algae were determined and combined with rbcL sequences obtained from GenBank for four other heterokont algae and three red algae. The phylogeny of the morphologically diverse haterokont algae was inferred from an unambiguously aligned data matrix using the red algae as the root, Significantly higher levels of mutational saturation in third codon positions were found when plotting the pair-wise substitutions with and without corrections for multiple substitutions at the same site for first and second codon positions only and for third positions only. In light of this observation, third codon positions were excluded from phylogenetic analyses. Both weighted-parsimony and maximum-likelihood analyses supported with high bootstrap values the monophyly of the nine currently recognized classes of heterokont algae. The Eustigmatophyceae were the most basal group, and the Dictyochophyceae branched off as the second most basal group. The branching pattern for the other classes was well supported in terms of bootstrap values in the weightedparsimony analysis but was weakly supported in the maximum-likelihood analysis (<50%). In the parsimony analysis, the diatoms formed a sister group to the branch containing the Chrysophyceae and Synurophyceae. This clade, charactetized by siliceous structures (frustules, cysts, scales), was the sister group to the Pelagophyceae/Sarcinochrysidales and Phaeo-/Xantho-/ Raphidophyceae clades. In the latter clade, the raphido-phytes were sister to the Phaeophyceae and Xanthophyceae. A relative rate test revealed that the rbcL gene in the Chrysophyceae and Synurophyceae has experienced a significantly different rate of substitutions compared to other classes of heterokont algae. The branch lengths in the maximum-likelihood reconstruction suggest that these two classes have evolved at an accelerated rate. Six major carotenoids were analyzed cladistically to study the usefulness of carotenoid pigmentation as a class-level character in the heterokont algae. In addition, each carotenoid was mapped onto both the rbcL tree and a consensus tree derived from nuclear-encoded small-subunit ribosomal DNA (SSU rDNA) sequences. Carotenoid pigmentation does not provide unambiguous phylogenetic information, whether analyzed cladistically by itself or when mapped onto phylogenetic trees based upon molecular sequence data.     

Rate and mode differences between nuclear and mitochondrial small-subunit rRNA genes in mushrooms.

Sequences from homologous regions of the nuclear and mitochondrial small-subunit rRNA genes from 10 members of the mushroom order Boletales were used to construct evolutionary trees and to compare the rates and modes of evolution. Trees constructed independently for each gene by parsimony and tested by bootstrap analysis have identical topologies in all statistically significant branches. Examination of base substitutions revealed that the nuclear gene is biased toward C-T transitions and that the distribution of transversions in the mitochondrial gene is strongly effected by an A-T bias. When only homologous regions of the two genes were compared, base substitutions per nucleotide were roughly 16-fold greater in the mitochondrial gene. The difference in the frequency of length mutations was at least as great but was impossible to estimate accurately because of their absence in the nuclear gene. Maximum likelihood was used to show that base-substitution rates vary dramatically among the branches. A significant part of the rate inconstancy was caused by an accelerated nuclear rate in one branch and a retarded mitochondrial rate in a different branch. A second part of the rate variability involved a consistent inconstancy: short branches exhibit ratios of mitochondrial to nuclear divergences of less than 1, while longer branches had ratios of approximately 4:1-8:1. This pattern suggests a systematic error in the branch length calculation. The error may be related to the simplicity of the divergence estimates, which assumes that all base positions have an equal probability of change.     

A comparison of cytochrome b and DNA hybridization data bearing on the phylogeny of swallows (Aves: Hirundinidae)

Cytochrome b sequence data from 17 species representing 16 genera of swallows (Aves: Hirundinidae) were compared with DNA-DNA hybridization data from the same species in a taxonomic congruence assessment of swallow phylogeny. In the process, subsets (partitions) of the cytochrome b sequence data were examined in light of the DNA hybridization distances to assess their potential phylogenetic informativeness. When the sequence data were weighted-with or without reference to the DNA hybridization data-they produced parsimony and maximum likelihood (but not distance) trees that were largely congruent with the DNA hybridization tree. To this extent, the cytochrome b data supported many of the phylogenetic conclusions based on the DNA hybridization tree and vice versa. However, the cytochrome b data produced largely unresolved trees when branch robustness was tested by bootstrapping and other methods. This poor resolution appeared to be caused by a lack of hierarchical structure in the cytochrome b distances, which were confined to a narrow range (between 10-13%), compressed by saturation, and noisy. Partition analysis by codon sites and protein domains yielded typical avian cytochrome b patterns, except for idiosyncrasies attributable to the genetic divergence level of swallows in comparison to other groups of birds whose cytochrome b sequences have been analyzed.     

The likelihood node density effect and consequences for evolutionary studies of molecular rates

Many molecular phylogenies show longer root-to-tip path lengths in species-rich groups, encouraging hypotheses linking cladogenesis with accelerated molecular evolution. However, the pattern can also be caused by an artifact called the node density effect (NDE): this effect occurs when the method used to reconstruct a tree underestimates multiple hits that would have been revealed by extra nodes, leading to longer root-to-tip path lengths in clades with more terminal taxa. Here we use a twofold approach to demonstrate that maximum likelihood and Bayesian methods also suffer from the NDE known to affect parsimony. First, simulations deliberately mismatching the simulation and reconstruction models show that the greater the model disparity, the greater the gap between actual and reconstructed tree lengths, and the greater the NDE. Second, taxon sampling manipulation with empirical data shows that NDE can still be present when using optimized models: across 12 datasets, 70 out of 109 sister path comparisons showed significant evidence of NDE. Unless the model fairly accurately reconstructs the real tree length-and given the complexity of real sequence evolution this may be uncommon -- it will consistently produce a node density artifact. At commonly encountered divergence levels, a 10% underestimation of tree length results in > or = 80% of simulated phylogenies showing a positive NDE. Bayesian trees have a slight but consistently stronger effect. This pervasive methodological artifact increases apparent rate heterogeneity, and can compromise investigations of factors influencing molecular evolutionary rate that use path lengths in topologically asymmetric trees.     

Determining evolutionary distances from highly diverged nucleic acid sequences: Operator metrics

Summary Operator metrics are explicity designed to measure evolutionary distances from nucleic acid sequences when substitution rates differ greatly among the organisms being compared, or when substitutions have been extensive. Unlike lengths calculated by the distance matrix and parsimony methods, in which substitutions in one branch of a tree can alter the measured length of another branch, lengths determined by operator metrics are not affected by substitutions outside the branch.In the method, lengths (operator metrics) corresponding to each of the branches of an unrooted tree are calculated. The metric length of a branch reconstructs the number of (transversion) differences between sequences at a tip and a node (or between nodes) of a tree. The theory is general and is fundamentally independent of differences in substitution rates among the organisms being compared. Mathematically, the independence has been obtained becuase the metrics are eigen vectors of fundamental equations which describe the evolution of all unrooted trees.Even under conditions when both the distance matrix method or a simple parsimony length method are show to indicate lengths than are an order of magnitude too large or too small, the operator metrics are accurate. Examples, using data calculated with evolutionary rates and branchings designed to confuse the measurement of branch lengths and to camouflage the topology of the true tree, demonstrate the validity of operator metrics. The method is robust. Operator metric distances are easy to calculated, can be extended to any number of taxa, and provide a statistical estimate of their variances.The utility of the method is demonstrated by using it to analyze the origins and evolutionary of chloroplasts, mitochondria, and eubacteria.     

Phylogeny Derived from Coding Retroviral Genome Organization

When divergence between viral species is large, the analysis and comparison of nucleotide or protein sequences are dependent on mutation biases and multiple substitutions per site leading, among other things, to the underestimation of branch lengths in phylogenetic trees. To avoid the problem of multiply substituted sites, a method not directly based on the nucleic or protein sequences has been applied to retroviruses. It consisted of asking questions about genome structure or organization, and gene function, the series of answers creating coded sequences analyzed by phylogenic software. This method recovered the principal retroviral groups such as the lentiviruses and spumaviruses and highlighted questions and answers characteristic of each group of retroviruses. In general, there was reasonable concordance between the coded genome methodology and that based on conventional phylogeny of the integrase protein sequence, indicating that integrase was fixing mutations slowly enough to marginalize the problem of multiple substitutions at sites. To a first approximation, this suggests that the acquisition of novel genetic features generally parallels the fixation of amino acid substitutions. Received: 18 May 2001 / Accepted: 7 September 2001     

The relative contribution of band number to phylogenetic accuracy in AFLP data sets

We examined the effect of increasing the number of sampled amplified fragment length polymorphism (AFLP) bands to reconstruct an accurate and well-supported AFLP-based phylogeny. In silico AFLP was performed using simulated DNA sequences evolving along balanced and unbalanced model trees with recent, uniform and ancient radiations and average branch lengths (from the most internal node to the tip) ranging from 0.02 to 0.05 substitutions per site. Trees were estimated by minimum evolution (ME) and maximum parsimony (MP) methods from both DNA sequences and virtual AFLP fingerprints. The comparison of the true tree with the estimated AFLP trees suggests that moderate numbers of AFLP bands are necessary to recover the correct topology with high bootstrap support values (i.e. >70%). Fewer numbers of bands are necessary for shorter tree lengths and for balanced than for unbalanced tree topologies. However, branch length estimation was rather unreliable and did not improve substantially after a certain number of bands were sampled. These results hold for different levels of genome coverage and number of taxa analysed. In silico AFLP using bacterial genomic DNA sequences recovered a well-supported tree topology that mirrored an empirical phylogeny based on a set of 31 orthologous gene sequences when as few as 263 AFLP bands were scored. These results suggest that AFLPs may be an efficient alternative to traditional DNA sequencing for accurate topology reconstruction of shallow trees when not very short ancestral branches exist.     

A rate-independent technique for analysis of nucleic acid sequences: evolutionary parsimony

The method of evolutionary parsimony--or operator invariants--is a technique of nucleic acid sequence analysis related to parsimony analysis and explicitly designed for determining evolutionary relationships among four distantly related taxa. The method is independent of substitution rates because it is derived from consideration of the group properties of substitution operators rather than from an analysis of the probabilities of substitution in branches of a tree. In both parsimony and evolutionary parsimony, three patterns of nucleotide substitution are associated one-to-one with the three topologically linked trees for four taxa. In evolutionary parsimony, the three quantities are operator invariants. These invariants are the remnants of substitutions that have occurred in the interior branch of the tree and are analogous to the substitutions assigned to the central branch by parsimony. The two invariants associated with the incorrect trees must equal zero (statistically), whereas only the correct tree can have a nonzero invariant. The chi 2-test is used to ascertain the nonzero invariant and the statistically favored tree. Examples, obtained using data calculated with evolutionary rates and branchings designed to camouflage the true tree, show that the method accurately predicts the tree, even when substitution rates differ greatly in neighboring peripheral branches (conditions under which parsimony will consistently fail). As the number of substitutions in peripheral branches becomes fewer, the parsimony and the evolutionary-parsimony solutions converge. The method is robust and easy to use.      

Molecular phylogeny and evidence for an adaptive radiation of geophagine cichlids from South America (Perciformes: Labroidei)

Nucleotide sequences from the mitochondrial ND4 gene and the nuclear RAG2 gene were used to derive the most extensive molecular phylogeny to date for the Neotropical cichlid subfamily Geophaginae. Previous hypotheses of relationships were tested in light of these new data and a synthesis of all existing molecular information was provided. Novel phylogenetic findings included support for : (1) a 'Big Clade' containing the genera Geophagus sensu lato, Gymnogeophagus, Mikrogeophagus, Biotodoma, Crenicara, and Dicrossus; (2) a clade including the genera Satanoperca, Apistogramma, Apistogrammoides, and Taeniacara; and (3) corroboration for Kullander's clade Acarichthyini. ND4 demonstrated saturation effects at the third code position and lineage-specific rate heterogeneity, both of which influenced phylogeny reconstruction when only equal weighted parsimony was employed. Both branch lengths and internal branch tests revealed extremely short basal nodes that add support to the idea that geophagine cichlids have experienced an adaptive radiation sensu Schluter that involved ecomorphological specializations and life history diversification.     

Early evolutionary relationships among known life forms inferred from elongation factor EF-2/EF-G sequences: Phylogenetic coherence and structure of the archaeal domain

Summary Phylogenies were inferred from both the gene and the protein sequences of the translational elongation factor termed EF-2 (for Archaea and Eukarya) and EF-G (for Bacteria). All treeing methods used (distance-matrix, maximum likelihood, and parsimony), including evolutionary parsimony, support the archaeal tree and disprove the eocyte tree (i.e., the polyphyly and paraphyly of the Archaea). Distance-matrix trees derived from both the amino acid and the DNA sequence alignments (first and second codon positions) showed the Archaea to be a monophyletia-holophyletic grouping whose deepest bifurcation divides a Sulfolobus branch from a branch comprising Methanococcus, Halobacterium, and Thermoplasma. Bootstrapped distance-matrix treeing confirmed the monophyly-holophyly of Archaea in 100% of the samples and supported the bifurcation of Archaea into a Sulfolobus branch and a methanogen-halophile branch in 97% of the samples. Similar phylogenies were inferred by maximum likelihood and by maximum (protein and DNA) parsimony. DNA parsimony trees essentially identical to those inferred from first and second codon positions were derived from alternative DNA data sets comprising either the first or the second position of each codon. Bootstrapped DNA parsimony supported the monophyly-holophyly of Archaea in 100% of the bootstrap samples and confirmed the division of Archaea into a Sulfolobus branch and a methanogen-halophile branch in 93% of the bootstrap samples. Distance-matrix and maximum likelihood treeing under the constraint that branch lengths must be consistent with a molecular clock placed the root of the universal tree between the Bacteria and the bifurcation of Archaea and Eukarya. The results support the division of Archaea into the kingdoms Crenarchaeota (corresponding to the Sulfolobus branch and Euryarchaeota). This division was not confirmed by evolutionary parsimony, which identified Halobacterium rather than Sulfolobus as the deepest offspring within the Archaea.Offprint requests to: P. Cammarano     

Evolutionary change of restriction cleavage sites and phylogenetic inference for man and apes

A mathematical theory for the evolutionary change of restriction endonuclease cleavage sites is developed, and the probabilities of various types of restriction-site changes are evaluated. A computer simulation is also conducted to study properties of the evolutionary change of restriction sites. These studies indicate that parsimony methods of constructing phylogenetic trees often make erroneous inferences about evolutionary changes of restriction sites unless the number of nucleotide substitutions per site is less than 0.01 for all branches of the tree. This introduces a systematic error in estimating the number of mutational changes for each branch and, consequently, in constructing phylogenetic trees. Therefore, parsimony methods should be used only in cases where nucleotide sequences are closely related. Reexamination of Ferris et al.'s data on restriction-site differences of mitochondrial DNAs does not support Templeton's conclusions regarding the phylogenetic tree for man and apes and the molecular clock hypothesis. Templeton's claim that Nei and Li's method of estimating the number of nucleotide substitutions per site is seriously affected by parallel losses and loss-gains of restriction sites is also unsupported.