DNA DIVERGENCE AMONG HOMINOIDS

We have determined the degree of single-copy DNA divergence among the extant members of the Hominoidea employing the technique of DNA-DNA hybridization. The species studied include humans, two species of chimpanzees, gorillas, two subspecies of orangutans, and two species of gibbons; as an outgroup we have used a member of the Old World monkeys (Cercopithecidae), the baboon. Our methods are different from those previously used and allow us to control for two factors other than base-pair mismatch that can affect the thermal stability of DNA duplexes: the base composition and duplex length. In addition, we have studied more than one individual for most species and thus are able to assess the effect of intraspecific variation on phylogenetic conclusions. The results indicate that the closest extant relatives of humans are the chimpanzees. Gorillas are the next closest, followed by orangutans and gibbons. This result is strongly supported statistically, as there is virtually no overlap in measurements between different taxa. Our conclusions are in agreement with a growing amount of molecular evidence supporting this pattern of relatedness. The data behave as a reasonably good molecular clock, and we do not see an indication of slowdown in molecular evolution in the clade containing humans and African apes, contrary to what has been documented for protein-coding regions. Because of the clocklike nature of the results, we have estimated that the divergence of humans and chimpanzees occurred about 6-8 million years ago. Results from orangutans indicate that the Borneo and Sumatra populations are genetically distinct, about as different as the named species of chimpanzees.

Parameterizing sequence alignment with an explicit evolutionary model

Background

Inference of sequence homology is inherently an evolutionary question, dependent upon evolutionary divergence. However, the insertion and deletion penalties in the most widely used methods for inferring homology by sequence alignment, including BLAST and profile hidden Markov models (profile HMMs), are not based on any explicitly time-dependent evolutionary model. Using one fixed score system (BLOSUM62 with some gap open/extend costs, for example) corresponds to making an unrealistic assumption that all sequence relationships have diverged by the same time. Adoption of explicit time-dependent evolutionary models for scoring insertions and deletions in sequence alignments has been hindered by algorithmic complexity and technical difficulty.

Results

We identify and implement several probabilistic evolutionary models compatible with the affine-cost insertion/deletion model used in standard pairwise sequence alignment. Assuming an affine gap cost imposes important restrictions on the realism of the evolutionary models compatible with it, as single insertion events with geometrically distributed lengths do not result in geometrically distributed insert lengths at finite times. Nevertheless, we identify one evolutionary model compatible with symmetric pair HMMs that are the basis for Smith-Waterman pairwise alignment, and two evolutionary models compatible with standard profile-based alignment.

We test different aspects of the performance of these “optimized branch length” models, including alignment accuracy and homology coverage (discrimination of residues in a homologous region from nonhomologous flanking residues). We test on benchmarks of both global homologies (full length sequence homologs) and local homologies (homologous subsequences embedded in nonhomologous sequence).

Conclusions

Contrary to our expectations, we find that for global homologies a single long branch parameterization suffices both for distant and close homologous relationships. In contrast, we do see an advantage in using explicit evolutionary models for local homologies. Optimal branch parameterization reduces a known artifact called “homologous overextension”, in which local alignments erroneously extend through flanking nonhomologous residues.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0832-5) contains supplementary material, which is available to authorized users.

Analyses of mitochondrial amino acid sequence datasets support the proposal that specimens of Hypodontus macropi from three species of macropodid hosts represent distinct species

BACKGROUND

Hypodontus macropi is a common intestinal nematode of a range of kangaroos and wallabies (macropodid marsupials). Based on previous multilocus enzyme electrophoresis (MEE) and nuclear ribosomal DNA sequence data sets, H. macropi has been proposed to be complex of species. To test this proposal using independent molecular data, we sequenced the whole mitochondrial (mt) genomes of individuals of H. macropi from three different species of hosts (Macropus robustus robustus, Thylogale billardierii and Macropus [Wallabia] bicolor) as well as that of Macropicola ocydromi (a related nematode), and undertook a comparative analysis of the amino acid sequence datasets derived from these genomes.

RESULTS

The mt genomes sequenced by next-generation (454) technology from H. macropi from the three host species varied from 13,634 bp to 13,699 bp in size. Pairwise comparisons of the amino acid sequences predicted from these three mt genomes revealed differences of 5.8% to 18%. Phylogenetic analysis of the amino acid sequence data sets using Bayesian Inference (BI) showed that H. macropi from the three different host species formed distinct, well-supported clades. In addition, sliding window analysis of the mt genomes defined variable regions for future population genetic studies of H. macropi in different macropodid hosts and geographical regions around Australia.

CONCLUSIONS

The present analyses of inferred mt protein sequence datasets clearly supported the hypothesis that H. macropi from M. robustus robustus, M. bicolor and T. billardierii represent distinct species.

Efficient Methods for Estimating Amino Acid Replacement Rates

Replacement rate matrices describe the process of evolution at one position in a protein and are used in many applications where proteins are studied with an evolutionary perspective. Several general matrices have been suggested and have proved to be good approximations of the real process. However, there are data for which general matrices are inappropriate, for example, special protein families, certain lineages in the tree of life, or particular parts of proteins. Analysis of such data could benefit from adaption of a data-specific rate matrix. This paper suggests two new methods for estimating replacement rate matrices from independent pairwise protein sequence alignments and also carefully studies Müller-Vingron's resolvent method. Comprehensive tests on synthetic datasets show that both new methods perform better than the resolvent method in a variety of settings. The best method is furthermore demonstrated to be robust on small datasets as well as practical on very large datasets of real data. Neither short nor divergent sequence pairs have to be discarded, making the method economical with data. A generalization to multialignment data is suggested and used in a test on protein-domain family phylogenies, where it is shown that the method offers family-specific rate matrices that often have a significantly better likelihood than a general matrix.

Recreating ancestral proteins

Tracing the history of molecular changes using phylogenetic methods can provide powerful insights into how and why molecules work the way they do. It is now possible to recreate inferred ancestral proteins in the laboratory and study the function of these molecules. This provides a unique opportunity to study the paths and the mechanisms of functional change during molecular evolution. What insights have already emerged from such phylogenetic studies of protein evolution and function, what are the impediments to progress and what are the prospects for the future?

Phylogenetic information and experimental design in molecular systematics

Despite the widespread perception that evolutionary inference from molecular sequences is a statistical problem, there has been very little attention paid to questions of experimental design. Previous consideration of this topic has led to little more than an empirical folklore regarding the choice of suitable genes for analysis, and to dispute over the best choice of taxa for inclusion in data sets. I introduce what I believe are new methods that permit the quantification of phylogenetic information in a sequence alignment. The methods use likelihood calculations based on Markov-process models of nucleotide substitution allied with phylogenetic trees, and allow a general approach to optimal experimental design. Two examples are given, illustrating realistic problems in experimental design in molecular phylogenetics and suggesting more general conclusions about the choice of genomic regions, sequence lengths and taxa for evolutionary studies.

Probability distribution of molecular evolutionary trees: a new method of phylogenetic inference

A new method is presented for inferring evolutionary trees using nucleotide sequence data. The birth-death process is used as a model of speciation and extinction to specify the prior distribution of phylogenies and branching times. Nucleotide substitution is modeled by a continuous-time Markov process. Parameters of the branching model and the substitution model are estimated by maximum likelihood. The posterior probabilities of different phylogenies are calculated and the phylogeny with the highest posterior probability is chosen as the best estimate of the evolutionary relationship among species. We refer to this as the maximum posterior probability (MAP) tree. The posterior probability provides a natural measure of the reliability of the estimated phylogeny. Two example data sets are analyzed to infer the phylogenetic relationship of human, chimpanzee, gorilla, and orangutan. The best trees estimated by the new method are the same as those from the maximum likelihood analysis of separate topologies, but the posterior probabilities are quite different from the bootstrap proportions. The results of the method are found to be insensitive to changes in the rate parameter of the branching process.

Phylogenetic analysis using parsimony and likelihood methods

The assumptions underlying the maximum-parsimony (MP) method of phylogenetic tree reconstruction were intuitively examined by studying the way the method works. Computer simulations were performed to corroborate the intuitive examination. Parsimony appears to involve very stringent assumptions concerning the process of sequence evolution, such as constancy of substitution rates between nucleotides, constancy of rates across nucleotide sites, and equal branch lengths in the tree. For practical data analysis, the requirement of equal branch lengths means similar substitution rates among lineages (the existence of an approximate molecular clock), relatively long interior branches, and also few species in the data. However, a small amount of evolution is neither a necessary nor a sufficient requirement of the method. The difficulties involved in the application of current statistical estimation theory to tree reconstruction were discussed, and it was suggested that the approach proposed by Felsenstein (1981, J. Mol. Evol. 17: 368-376) for topology estimation, as well as its many variations and extensions, differs fundamentally from the maximum likelihood estimation of a conventional statistical parameter. Evidence was presented showing that the Felsenstein approach does not share the asymptotic efficiency of the maximum likelihood estimator of a statistical parameter. Computer simulations were performed to study the probability that MP recovers the true tree under a hierarchy of models of nucleotide substitution; its performance relative to the likelihood method was especially noted. The results appeared to support the intuitive examination of the assumptions underlying MP. When a simple model of nucleotide substitution was assumed to generate data, the probability that MP recovers the true topology could be as high as, or even higher than, that for the likelihood method. When the assumed model became more complex and realistic, e.g., when substitution rates were allowed to differ between nucleotides or across sites, the probability that MP recovers the true topology, and especially its performance relative to that of the likelihood method, generally deteriorates. As the complexity of the process of nucleotide substitution in real sequences is well recognized, the likelihood method appears preferable to parsimony. However, the development of a statistical methodology for the efficient estimation of the tree topology remains a difficult open problem.

Time and biosequences

In this paper we discuss and demonstrate the importance of several factors relative to the relationship between time and evolution of biosequences. In both quantitative and qualitative measurements of the genetic distances, the compositional constraints of the nucleotide sequences play a very important role. We demonstrate that when homologous sequences significantly differ in base composition we get erratic branching order and/or wrong evaluation of the evolutionary rates. We must consider that every gene may have a different evolutionary dynamic along its sequence, generally linked to its functional constraints; this too can seriously affect its clock-like behavior. We report some cases showing how these factors can affect the quantitative measurements of the genetic distances of biosequences.

Statistical tests of models of DNA substitution

Penny et al. have written that "The most fundamental criterion for a scientific method is that the data must, in principle, be able to reject the model. Hardly any [phylogenetic] tree-reconstruction methods meet this simple requirement." The ability to reject models is of such great importance because the results of all phylogenetic analyses depend on their underlying models--to have confidence in the inferences, it is necessary to have confidence in the models. In this paper, a test statistic suggested by Cox is employed to test the adequacy of some statistical models of DNA sequence evolution used in the phylogenetic inference method introduced by Felsenstein. Monte Carlo simulations are used to assess significance levels. The resulting statistical tests provide an objective and very general assessment of all the components of a DNA substitution model; more specific versions of the test are devised to test individual components of a model. In all cases, the new analyses have the additional advantage that values of phylogenetic parameters do not have to be assumed in order to perform the tests.

Heterogeneity in the substitution process of amino acid sites of proteins coded for by mitochondrial DNA

Several forms of maximum likelihood models are applied to aligned amino acid sequence data coded for in the mitochondrial DNA of six species (chicken, frog, human, bovine, mouse, and rat). These models range in form from relatively simple models of the type currently used for inferring phylogenetic tree structure to models more complex than those that have been used previously. No major discrepancies between the optimal trees inferred by any of these methods are found, but there are huge differences in adequacy of fit. A very significant finding is that the fit of any of these models is vastly improved by allowing a certain proportion of the amino acid sites to be invariant. An even more important, although disquieting, finding is that none of these models fits well, as judged by standard statistical criteria. The primary reason for this is that amino acid sites undergo substitution according to a process that is very heterogeneous. Because most phylogenetic inference is accomplished by choosing the optimal tree under the assumption that a homogeneous process is acting on the sites, the potential invalidity of some such conclusions is raised by this article's results. The seriousness of this problem depends upon the robustness of the phylogenetic inferential procedure to departures from the underlying model.

Different rates of substitution may produce different phylogenies of the eutherian mammals

In an attempt to resolve some points of branching order in the phylogeny of the eutherian mammals, a phylogenetic analysis of 26 nuclear and 6 mitochondrial genes was undertaken using a maximum likelihood method on a constant rate stochastic model of molecular evolution. Seventeen of the nuclear genes gave a primates/artiodactyls grouping highest support whereas three of the mitochondrial genes found a rodents/artiodactyls grouping to be best supported. The primates/rodents grouping was never the best supported. On the assumption that rodents are indeed an outgroup to primates and artiodactyls and that the latter taxa diverged 70 million years ago, an estimation was made, for each gene, of the time of divergence of the rodent lineage. In most cases such estimates were beyond the limits set by present interpretations of the paleontological record as were many estimates of the divergence time of mouse and rat. These results suggest that, although there is locus variation, the divergent position of the rodent lineage may be an artifact of an elevated rate of nucleotide substitution in this order.

Independent gene evolution in the potato actin gene family demonstrated by phylogenetic procedures for resolving gene conversions and the phylogeny of angiosperm actin genes

Nine different actin DNA sequences were isolated from the common potato, Solanum tuberosum, and the nucleotide sequence of five actin loci and of two allelic variants are presented. Unlike the wide variation in intron position among animal actin genes, the potato actin genes have three introns situated in the same positions as reported for all other angiosperm actin genes. Using a novel combination of analytical procedures (G-test and compatibility analysis), we could not find evidence of frequent large or small nonreciprocal exchanges of genetic material between the sequenced loci, although there were a few candidates. Resolution of such gene conversion events and the quantification of independence of gene evolution in multigene families is critical to the inference of phylogenetic relationships. Comparison with actin genes in other angiosperm species suggests that the actin multigene family can be divided into a number of subfamilies, evolved by descent rather than gene conversion, which are of possible functional origin, with one major subfamily diversification occurring before the divergence of monocots and dicots. The silent rate of nucleotide substitution was estimated to be similar to that suggested for a number of other plant nuclear genes, whereas the replacement rate was extremely slow, suggestive of selective constraints.

Evaluation of the maximum likelihood estimate of the evolutionary tree topologies from DNA sequence data, and the branching order in hominoidea

A maximum likelihood method for inferring evolutionary trees from DNA sequence data was developed by Felsenstein (1981). In evaluating the extent to which the maximum likelihood tree is a significantly better representation of the true tree, it is important to estimate the variance of the difference between log likelihood of different tree topologies. Bootstrap resampling can be used for this purpose (Hasegawa et al. 1988; Hasegawa and Kishino 1989), but it imposes a great computation burden. To overcome this difficulty, we developed a new method for estimating the variance by expressing it explicitly. The method was applied to DNA sequence data from primates in order to evaluate the maximum likelihood branching order among Hominoidea. It was shown that, although the orangutan is convincingly placed as an outgroup of a human and African apes clade, the branching order among human, chimpanzee, and gorilla cannot be determined confidently from the DNA sequence data presently available when the evolutionary rate constancy is not assumed.

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

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.

DNA hybridization evidence of hominoid phylogeny: results from an expanded data set

The living hominoids are human, the two species of chimpanzees, gorilla, orangutan, and nine species of gibbons. The cercopithecoids (Old World monkeys) are the sister group of the hominoids. A consensus about the phylogeny of the hominoids has been reached for the branching order of the gibbons (earliest) and the orangutan (next earliest), but the branching order among gorilla, chimpanzees, and human remains in contention. In 1984 we presented DNA-DNA hybridization data, based on 183 DNA hybrids, that we interpreted as evidence that the branching order, from oldest to most recent, was gibbons, orangutan, gorilla, chimpanzees, and human. In the present paper we report on an expanded data set totaling 514 DNA hybrids, which supports the branching order given above. The ranges for the datings of divergence nodes are Old World monkeys, 25-34 million years (Myr) ago; gibbons, 16.4-23 Myr ago; orangutan, 12.2-17 Myr ago; gorilla, 7.7-11 Myr ago; chimpanzees-human, 5.5-7.7 Myr ago. The possible effects of differences in age at first breeding are discussed, and some speculations about average genomic rates of evolution are presented.