Search

A phylogenetic tree can reflect the evolutionary relationships between species or gene families, and they play a critical role in modern biological research. In this review, we summarize common methods for constructing phylogenetic trees, including distance methods, maximum parsimony, maximum likelihood, Bayesian inference, and tree-integration methods (supermatrix and supertree). Here we discuss the advantages, shortcomings, and applications of each method and offer relevant codes to construct phylogenetic trees from molecular data using packages and algorithms in R. This review aims to provide comprehensive guidance and reference for researchers seeking to construct phylogenetic trees while also promoting further development and innovation in this field. By offering a clear and concise overview of the different methods available, we hope to enable researchers to select the most appropriate approach for their specific research questions and datasets. [ABSTRACT FROM AUTHOR]

A phylogenetic tree can reflect the evolutionary relationships between species or gene families, and they play a critical role in modern biological research. In this review, we summarize common methods for constructing phylogenetic trees, including distance methods, maximum parsimony, maximum likelihood, Bayesian inference, and tree-integration methods (supermatrix and supertree). Here we discuss the advantages, shortcomings, and applications of each method and offer relevant codes to construct phylogenetic trees from molecular data using packages and algorithms in R. This review aims to provide comprehensive guidance and reference for researchers seeking to construct phylogenetic trees while also promoting further development and innovation in this field. By offering a clear and concise overview of the different methods available, we hope to enable researchers to select the most appropriate approach for their specific research questions and datasets.

Maximum parsimony and maximum likelihood are two contrasting approaches for reconstructing phylogenetic trees from sequence and character data. We establish analytic links between these methods (extending connections reported earlier) under the simple Poisson model of substitutions in two settings. First, we show that if the underlying state space is sufficiently large, then the maximum likelihood estimate phylogenetic tree is always a maximum parsimony tree for the data. Second, we show that a sufficiently dense sampling of sequences ensures that the most parsimonious likelihood tree is always a maximum parsimony tree. [Copyright &y& Elsevier]

Chou-Fasman parameters, measuring preferences of each amino acid for different conformational regions in proteins, were used to obtain an amino acid difference index of conformational parameter distance (CPD) values. CPD values were found to be significantly lower for amino acid exchanges representing in the genetic code transitions of purines, G⇆A than for exchanges representing either transitions of pyrimidines, C⇆U, or transversions of purines and pyrimidines. Inasmuch as the distribution of CPD values in these non G⇆A exchanges resembles that obtained for amino acid pairs with double or triple base differences in their underlying codons, we conclude that the genetic code was not particularly designed to minimize effects of mutation on protein conformation. That natural selection minimizes these changes, however, was shown by tabulating results obtained by the maximum parsimony method for eight protein genealogies with a total occurrence of 4574 base substitutions. At the beginning position of the codons G⇆A transitions were in very great excess over other base substitutions, and, conversely, C⇆U transitions were deficient. At the middle position of the codons only fast evolving proteins showed an excess of G⇆A transitions, as though selection mainly preserved conformation in these proteins while weeding out mutations affecting chemical properties of functional sites in slow evolving proteins. In both fast and slow evolving proteins the net direction of transitions and transversions was found to be from G beginning codons to non-G beginning codons resulting in more commonly occurring amino acids, especially alanine with its generalized conformational properties, being replaced at suitable sites by amino acids with more specialized conformational and chemical properties. Historical circumstances pertaining to the origin of the genetic code and the nature of primordial proteins could account for such directional changes leading to increases in the functional density of proteins. In order to further explore the course of protein evolution, a modified parsimony algorithm was developed for constructing protein genealogies on the basis of minimum CPD length. The algorithm's ability to judge with finer discrimination that in protein evolution certain pathways of amino acid substitution should occur more readily than others was considered a potential advantage over strict maximum parsimony. In developing this CPD algorithm, the path of minimum CPD length through intermediate amino acids allowed by the genetic code for each pair of amino acids was determined. It was found that amino acid exchanges representing two base changes have a considerably lower average CPD value per base substitution than the amino acid exchanges representing single base changes. Amino acid exchanges representing three base changes have yet a further marked reduction in CPD per base change. This shows how extreme constraining effects of stabilizing selection can be circumvented, for by way of intermediate amino acids almost any amino acid can ultimately be substituted for another without damage to an evolving protein's conformation during the process. [ABSTRACT FROM AUTHOR]

Goodman et al's (1975) claim of accelerated evolution in the early stages of globin evolution is based on an erroneous assignment of the time of divergence of vertebrate myoglobin and hemoglobin. When this is corrected, there is no basis for their claim. The data are much more consistent with the nearly constant rate expected on the neutral mutation-random drift hypothesis than with the uneven rates expected if most amino acid changes were caused by substitution of favorable mutants through Darwinian selection. In addition, the majority of the codons determined by their maximum parsimony method have turned out to be wrong when compared to the actual nucleotide sequences of rabbit α and human β hemoglobins determined by direct sequencing. [ABSTRACT FROM AUTHOR]

Kimura mistook ambiguous maximum parsimony codons for wrong codons. The maximum parsimony method performed well as judged by the two classes of serine codons (which can not be connected by silent mutations) on comparing the parsimony codons for serines in human, rabbit, and mouse α hemoglobin chains to actual codons determined by nucleotide sequencing. In genealogical reconstructions involving 247 eucaryotic globins, the maximum parsimony distances separating the contemporary sequences show that Kimura's Poisson and Dayhoff's PAM estimates of rate of globin evolution miss most of the superimposed replacements and are therefore seriously in error. Nor is Kimura's constant rate assumption and his belief in a single origin of myoglobin supported. Lamprey myoglobin appears to be most like lamprey hemoglobin, while gnathostome myoglobin seems closest to gnathostome hemoglobin. It was found that the three types of gnathostome globins (Mb, α Hb, β Hb) evolved between the shark-boney vertebrate and bird-mammal ancestors at a much faster rate than from the latter ancestor to the present. The data indicate that rates were exceedingly fast during the origin of these globin chains because a high proportion of substitutions were adaptive. It was concluded that wherever strong stabilizing selection acts on a protein, somewhere in the past positive Darwinian selection must have spread the amino acid substitutions now being preserved. [ABSTRACT FROM AUTHOR]