SEARCH · PubMed Health

Results for “branch length estimation”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2Linked to original sources

PubMed2000

An updated and comprehensive rRNA phylogeny of (crown) eukaryotes based on rate-calibrated evolutionary distances.

Recent experience with molecular phylogeny has shown that all molecular markers have strengths and weaknesses. Nonetheless, despite several notable discrepancies with phylogenies obtained from protein data, the merits of the small subunit ribosomal RNA (SSU rRNA) as a molecular phylogenetic marker remain indisputable. Over the last 10 to 15 years a massive SSU rRNA database has been gathered, including more then 3000 complete sequences from eukaryotes. This creates a huge computational challenge, which is exacerbated by phenomena such as extensive rate variation among sites in the molecule. A few years ago, a fast phylogenetic method was developed that takes into account among-site rate variation in the estimation of evolutionary distances. This "substitution rate calibration" (SRC) method not only corrects for a major source of artifacts in phylogeny reconstruction but, because it is based on a distance approach, allows comprehensive trees including thousands of sequences to be constructed in a reasonable amount of time. In this study, a nucleotide variability map and a phylogenetic tree were constructed, using the SRC method, based on all available (January 2000) complete SSU rRNA sequences (2551) for species belonging to the so-called eukaryotic crown. The resulting phylogeny constitutes the most complete description of overall eukaryote diversity and relationships to date. Furthermore, branch lengths estimated with the SRC method better reflect the huge differences in evolutionary rates among and within eukaryotic lineages. The ribosomal RNA tree is compared with a recent protein phylogeny obtained from concatenated actin, alpha-tubulin, beta-tubulin, and elongation factor 1-alpha amino acid sequences. A consensus phylogeny of the eukaryotic crown based on currently available molecular data is discussed, as well as specific problems encountered in analyzing sequences when large differences in substitution rate are present, either between different sequences (rate variation among lineages) or between different positions within the same sequence (among-site rate variation).

Base Sequence↗

PubMed1995

Improved dating of the human/chimpanzee separation in the mitochondrial DNA tree: heterogeneity among amino acid sites.

The internal branch lengths estimated by distance methods such as neighbor-joining are shown to be biased to be short when the evolutionary rate differs among sites. The variable-invariable model for site heterogeneity fits the amino acid sequence data encoded by the mitochondrial DNA from Hominoidea remarkably well. By assuming the orangutan separation to be 13 or 16 Myr old, a maximum-likelihood analysis estimates a young date of 3.6 +/- 0.6 or 4.4 +/- 0.7 Myr (+/- 1 SE) for the human/chimpanzee separation, and these estimates turn out to be robust against differences in the assumed model for amino acid substitutions. Although some uncertainties still exist in our estimates, this analysis suggests that humans separated from chimpanzees some 4-5 Myr ago.

Amino Acid Sequence↗

PubMed1981

Simple method for constructing phylogenetic trees from distance matrices.

A simple method is proposed for constructing phylogenetic trees from distance matrices. The procedure for constructing tree topologies is similar to that of the unweighted pair-group method (UPG method) but makes corrections for unequal rates of evolution among lineages. The procedure for estimating branch lengths is the same as that of the Fitch and Margoliash method (F-M method) except that it allows no negative branch lengths. The performance of the present procedure for the construction of tree topologies is compared with that of the UPG method, the F-M method, Farris' method, and the modified Farris method by using Tateno's simulation outputs for nucleotide sequence divergence and his results for the performances of the latter four methods [Tateno, Y. (1978) Dissertation (Univ. Texas, Houston, TX). In this limited comparison, the present method performs considerably better than the UPG method and the F-M method and about equally well as the last two methods. The present method appears to be preferable to the UPG method for analysis of data from populations that have not differentiated much. Indeed, an application of the present method to gene frequency data from some Amerindian populations gives a tree topology far more reasonable than that obtained by the UPG method.

Biological Evolution↗

PubMed2000

War of the Iguanas: conflicting molecular and morphological phylogenies and long-branch attraction in iguanid lizards.

Recent studies based on different types of data (i.e., morphology, molecules) have found strongly conflicting phylogenies for the genera of iguanid lizards but have been unable to explain the basis for this incongruence. We reanalyze published data from morphology and from the mitochondrial ND4, cytochrome b, 12S, and 16S genes to explore the sources of incongruence and resolve these conflicts. Much of the incongruence centers on the genus Cyclura, which is the sister taxon of Iguana, according to parsimony analyses of the morphology and the ribosomal genes, but is the sister taxon of all other Iguanini, according to the protein-coding genes. Maximum likelihood analyses show that there has been an increase in the rate of nucleotide substitution in Cyclura in the two protein-coding genes (ND4 and cytochrome b), although this increase is not as clear when parsimony is used to estimate branch lengths. Parametric simulations suggest that Cyclura may be misplaced by the protein-coding genes as a result of long-branch attraction; even when Cyclura and Iguana are sister taxa in a simulated phylogeny, Cyclura is still placed as the basal member of the Iguanini by parsimony analysis in 55% of the replicates. A similar long-branch attraction problem may also exist in the morphological data with regard to the placement of Sauromalus with the Galápagos iguanas (Amblyrhynchus and Conolophus). The results have many implications for the analysis of diverse data sets, the impact of long branches on parsimony and likelihood methods, and the use of certain protein-coding genes in phylogeny reconstruction.

Animals↗

PubMed2000

Integrating ambiguously aligned regions of DNA sequences in phylogenetic analyses without violating positional homology.

Phylogenetic analyses of non-protein-coding nucleotide sequences such as ribosomal RNA genes, internal transcribed spacers, and introns are often impeded by regions of the alignments that are ambiguously aligned. These regions are characterized by the presence of gaps and their uncertain positions, no matter which optimization criteria are used. This problem is particularly acute in large-scale phylogenetic studies and when aligning highly diverged sequences. Accommodating these regions, where positional homology is likely to be violated, in phylogenetic analyses has been dealt with very differently by molecular systematists and evolutionists, ranging from the total exclusion of these regions to the inclusion of every position regardless of ambiguity in the alignment. We present a new method that allows the inclusion of ambiguously aligned regions without violating homology. In this three-step procedure, first homologous regions of the alignment containing ambiguously aligned sequences are delimited. Second, each ambiguously aligned region is unequivocally coded as a new character, replacing its respective ambiguous region. Third, each of the coded characters is subjected to a specific step matrix to account for the differential number of changes (summing substitutions and indels) needed to transform one sequence to another. The optimal number of steps included in the step matrix is the one derived from the pairwise alignment with the greatest similarity and the least number of steps. In addition to potentially enhancing phylogenetic resolution and support, by integrating previously nonaccessible characters without violating positional homology, this new approach can improve branch length estimations when using parsimony.

Alligators and Crocodiles↗

PubMed1996

A Hidden Markov Model approach to variation among sites in rate of evolution.

The method of Hidden Markov Models is used to allow for unequal and unknown evolutionary rates at different sites in molecular sequences. Rates of evolution at different sites are assumed to be drawn from a set of possible rates, with a finite number of possibilities. The overall likelihood of phylogeny is calculated as a sum of terms, each term being the probability of the data given a particular assignment of rates to sites, times the prior probability of that particular combination of rates. The probabilities of different rate combinations are specified by a stationary Markov chain that assigns rate categories to sites. While there will be a very large number of possible ways of assigning rates to sites, a simple recursive algorithm allows the contributions to the likelihood from all possible combinations of rates to be summed, in a time proportional to the number of different rates at a single site. Thus with three rates, the effort involved is no greater than three times that for a single rate. This "Hidden Markov Model" method allows for rates to differ between sites and for correlations between the rates of neighboring sites. By summing over all possibilities it does not require us to know the rates at individual sites. However, it does not allow for correlation of rates at nonadjacent sites, nor does it allow for a continuous distribution of rates over sites. It is shown how to use the Newton-Raphson method to estimate branch lengths of a phylogeny and to infer from a phylogeny what assignment of rates to sites has the largest posterior probability. An example is given using beta-hemoglobin DNA sequences in eight mammal species; the regions of high and low evolutionary rates are inferred and also the average length of patches of similar rates.

Animals↗

PubMed1987

Molecular evolution of mammalian class I alcohol dehydrogenase.

Phylogenetic relationship and the rates of evolution of mammalian alcohol dehydrogenases (ADHs) have been studied by using the amino acid sequences from the human (ADH alpha, ADH beta, and ADH gamma), rat, mouse, and horse (ADH E and ADH S). With the maize ADH1 and ADH2 used as references, the patterns of the amino acid replacements in the beta-sheets, alpha-helices, and random coils in each of the catalytic and coenzyme-binding domains were analyzed separately. The phylogenetic trees based on the different sets of amino acid substitutions consistently showed that (1) multiple ADHs in human and horse have arisen after mammalian radiation, (2) the common ancestor of human ADHs alpha and beta diverged from the ancestor of ADH gamma first and the former two ADHs diverged from each other more recently, and (3) the human ADHs are more closely related to the rodent ADHs than to the horse ADHs. Furthermore, the estimated branch lengths showed that the rodent ADHs are evolving faster than the other ADHs. This difference in evolutionary rate between the two groups of organisms is explainable either in terms of the difference in the number of cell generations per year or in terms of reduction of functional constraints.

Alcohol Dehydrogenase↗

PubMed2006

Parametric alignment of Drosophila genomes.

The classic algorithms of Needleman-Wunsch and Smith-Waterman find a maximum a posteriori probability alignment for a pair hidden Markov model (PHMM). To process large genomes that have undergone complex genome rearrangements, almost all existing whole genome alignment methods apply fast heuristics to divide genomes into small pieces that are suitable for Needleman-Wunsch alignment. In these alignment methods, it is standard practice to fix the parameters and to produce a single alignment for subsequent analysis by biologists. As the number of alignment programs applied on a whole genome scale continues to increase, so does the disagreement in their results. The alignments produced by different programs vary greatly, especially in non-coding regions of eukaryotic genomes where the biologically correct alignment is hard to find. Parametric alignment is one possible remedy. This methodology resolves the issue of robustness to changes in parameters by finding all optimal alignments for all possible parameters in a PHMM. Our main result is the construction of a whole genome parametric alignment of Drosophila melanogaster and Drosophila pseudoobscura. This alignment draws on existing heuristics for dividing whole genomes into small pieces for alignment, and it relies on advances we have made in computing convex polytopes that allow us to parametrically align non-coding regions using biologically realistic models. We demonstrate the utility of our parametric alignment for biological inference by showing that cis-regulatory elements are more conserved between Drosophila melanogaster and Drosophila pseudoobscura than previously thought. We also show how whole genome parametric alignment can be used to quantitatively assess the dependence of branch length estimates on alignment parameters.

Algorithms↗

PubMed1983

Accuracy of estimated phylogenetic trees from molecular data. II. Gene frequency data.

The accuracies and efficiencies of three different methods of making phylogenetic trees from gene frequency data were examined by using computer simulation. The methods examined are UPGMA, Farris' (1972) method, and Tateno et al.'s (1982) modified Farris method. In the computer simulation eight species (or populations) were assumed to evolve according to a given model tree, and the evolutionary changes of allele frequencies were followed by using the infinite-allele model. At the end of the simulated evolution five genetic distance measures (Nei's standard and minimum distances, Rogers' distance, Cavalli-Sforza's f theta, and the modified Cavalli-Sforza distance) were computed for all pairs of species, and the distance matrix obtained for each distance measure was used for reconstructing a phylogenetic tree. The phylogenetic tree obtained was then compared with the model tree. The results obtained indicate that in all tree-making methods examined the accuracies of both the topology and branch lengths of a reconstructed tree (rooted tree) are very low when the number of loci used is less than 20 but gradually increase with increasing number of loci. When the expected number of gene substitutions (M) for the shortest branch is 0.1 or more per locus and 30 or more loci are used, the topological error as measured by the distortion index (dT) is not great, but the probability of obtaining the correct topology (P) is less than 0.5 even with 60 loci. When M is as small as 0.004, P is substantially lower. In obtaining a good topology (small dT and high P) UPGMA and the modified Farris method generally show a better performance than the Farris method. The poor performance of the Farris method is observed even when Rogers' distance which obeys the triangle inequality is used. The main reason for this seems to be that the Farris method often gives overestimates of branch lengths. For estimating the expected branch lengths of the true tree UPGMA shows the best performance. For this purpose Nei's standard distance gives a better result than the others because of its linear relationship with the number of gene substitutions. Rogers' or Cavalli-Sforza's distance gives a phylogenetic tree in which the parts near the root are condensed and the other parts are elongated. It is recommended that more than 30 loci, including both polymorphic and monomorphic loci, be used for making phylogenetic trees. The conclusions from this study seem to apply also to data on nucleotide differences obtained by the restriction enzyme techniques.

Algorithms↗

PubMed1992

Interpretation of time constant and electrotonic length estimates in multicylinder or branched neuronal structures.

1. We have investigated the theoretical and practical problems associated with the interpretation of time constants and the estimation of electrotonic length with equivalent cylinder formulas for neurons best represented as multiple cylinders or branched structures. Two analytic methods were used to compute the time constants and coefficients of passive voltage transients (and time constants of current transients under voltage clamp). One method, suitable for simple geometries, involves analytic solutions to boundary value problems. The other, suitable for neurons of any geometric complexity, is an algebraic approach based on compartmental models. Neither of these methods requires the simulation of transients. 2. We computed the time constants and coefficients of voltage transients for several hypothetical neurons and also for a spinal motoneuron whose morphology was characterized from serial reconstructions. These time constants and coefficients were used to generate voltage transients. Then exponential peeling, nonlinear regression, and transform methods were applied to these transients to test how well these procedures estimate the underlying time constants and coefficients. 3. For a serially reconstructed motoneuron with 732 compartments, we found that the theoretical and peeled tau 0 values were nearly equal, but the theoretical tau 1 was much larger than the peeled tau 1. The theoretical tau 1 could not be peeled because it was associated with a coefficient, C1, that had a very small value. In fact, there were 156 time constants between 1.0 and 6.0 ms, most of which had very small coefficients; none had a coefficient larger than 2% of the signal. The peeled value of tau 1 (called tau 1 peel) can be viewed as some sort of a weighted average of the time constants having the largest coefficients. 4. We studied simple hypothetical neurons to determine what interpretation could be applied to the multitude of theoretical time constants. We found that after tau 0, there was a group of time constants associated with eigenfunctions that were odd (or approximately odd) functions with respect to the soma. These time constants could be interpreted as "equalizing" time constants along particular paths between different pairs of dendritic terminals in the neuron. After this group of time constants, there was one that we call tau even because it was associated with an eigen-function that was approximately even with respect to the soma. This tau even could be interpreted as an equalizing time constant for charge equalization between proximal membrane (soma and proximal dendrites) and distal membrane (including all distal dendrites).4=

Animals↗

PubMed1992

Statistical method for estimating the standard errors of branch lengths in a phylogenetic tree reconstructed without assuming equal rates of nucleotide substitution among different lineages.

A statistical method is developed for estimating the standard errors of branch lengths in a phylogenetic tree reconstructed without assuming equal rates of nucleotide substitution among different lineages. This method can be easily used for testing whether the length of an interior branch in a reconstructed tree is positive, i.e., whether the topology of the tree is correct. Computer simulations indicate that this method is appropriate for a statistical test. As an example, this method is applied to phylogenetic trees reconstructed for the four hominoid species: human, chimpanzee, gorilla, and orangutan. The results obtained show that the present method provides a powerful statistical test.

Animals↗

PubMed1994

Estimating errors and confidence intervals for branch lengths in phylogenetic trees by a bootstrap approach.

A method, based on the bootstrap procedure, is proposed for the estimation of branch-length errors and confidence intervals in a phylogenetic tree for which equal rates of substitution among lineages do not necessarily hold. The method can be used to test whether an estimated internodal distance is significantly greater than zero. In the application of the method, any estimator of genetic distances, as well as any tree reconstruction procedure (based on distance matrices), can be used. Also the method is not limited by the number of species involved in the phylogenetic tree. An example of the application of the method in the reconstruction of the phylogenetic tree for the four hominoid species--human, chimpanzee, gorilla, and orangutan--is shown.

Animals↗

PubMed1985

High-performance liquid chromatography of in vitro synthesized poly(ADP-ribose) on ion-exchange columns, separation of oligomers of varying chain length and estimation of apparent branching.

Separated macromolecular fractions of in vitro synthesized poly(ADP-ribose) by liver nuclei were subjected to ion-exchange chromatography in a programmed high-performance liquid chromatographic elution system. The effects of ionic strength, pH and temperature on the separation of poly(ADP-ribose) chains were determined. Short chain oligomers (up to n = 11) were fractionated into individual components by baseline separation. Each fraction was analyzed for chain length. Trace amounts of Ado(P)Rib(P)Rib(P) found in phosphodiesterase digests were taken as indication of apparent branching. In phosphodiesterase digests of the shorter oligomers, besides traces of the above component, two other digestion products were also observed, presumably representing oligomer termini, one terminal fragment being dominant in short oligomers. Medium and long chain oligomers were partly resolved to individual components, and especially the long oligomers exhibited marked temperature dependent elution patterns. Apparent branching increased with increasing chain length up to about 3% for n = 44 and components presumably indicating termini diminished to mere traces. The adenine spectra of all fractions identified individual components.

Chemical Phenomena↗

PubMed1989

A statistical test of phylogenies estimated from sequence data.

A simple approach to testing the significance of the branching order, estimated from protein or DNA sequence data, of three taxa is proposed. The branching order is inferred by the transformed-distance method, under the assumption that one or two outgroups are available, and the branch lengths are estimated by the least-squares method. The inferred branching order is considered significant if the estimated internodal distance is significantly greater than zero. To test this, a formula for the variance of the internodal distance has been developed. The statistical test proposed has been checked by computer simulation. The same test also applies to the case of four taxa with no outgroup, if one considers an unrooted tree. Formulas for the variances of internodal distances have also been developed for the case of five taxa. Conditions are given under which it is more efficient to add the sequence of a fifth taxon than to do 25% more nucleotide sequencing in each of the original four. A method is presented for combining analyses of disparate data to get a single P value. Finally, the test, applied to the human-chimpanzee-gorilla problem, shows that the issue is not yet resolved.

Animals↗

PubMed2001

Empirical and hierarchical Bayesian estimation of ancestral states.

Several methods have been proposed to infer the states at the ancestral nodes on a phylogeny. These methods assume a specific tree and set of branch lengths when estimating the ancestral character state. Inferences of the ancestral states, then, are conditioned on the tree and branch lengths being true. We develop a hierarchical Bayes method for inferring the ancestral states on a tree. The method integrates over uncertainty in the tree, branch lengths, and substitution model parameters by using Markov chain Monte Carlo. We compare the hierarchical Bayes inferences of ancestral states with inferences of ancestral states made under the assumption that a specific tree is correct. We find that the methods are correlated, but that accommodating uncertainty in parameters of the phylogenetic model can make inferences of ancestral states even more uncertain than they would be in an empirical Bayes analysis.

Animals↗

PubMed2004

The importance of proper model assumption in bayesian phylogenetics.

We studied the importance of proper model assumption in the context of Bayesian phylogenetics by examining >5,000 Bayesian analyses and six nested models of nucleotide substitution. Model misspecification can strongly bias bipartition posterior probability estimates. These biases were most pronounced when rate heterogeneity was ignored. The type of bias seen at a particular bipartition appeared to be strongly influenced by the lengths of the branches surrounding that bipartition. In the Felsenstein zone, posterior probability estimates of bipartitions were biased when the assumed model was underparameterized but were unbiased when the assumed model was overparameterized. For the inverse Felsenstein zone, however, both underparameterization and overparameterization led to biased bipartition posterior probabilities, although the bias caused by overparameterization was less pronounced and disappeared with increased sequence length. Model parameter estimates were also affected by model misspecification. Underparameterization caused a bias in some parameter estimates, such as branch lengths and the gamma shape parameter, whereas overparameterization caused a decrease in the precision of some parameter estimates. We caution researchers to assure that the most appropriate model is assumed by employing both a priori model choice methods and a posteriori model adequacy tests.

Bayes Theorem↗

PubMed2002

Variation in specific needle area of old-growth Douglas-fir in relation to needle age, within-crown position and epicormic shoot production.

Variation in specific needle area (SNA; cm2 projected fresh needle area g-1 oven-dried needle weight) was investigated in relation to needle age, within-crown position and epicormic shoot production in 450-year-old Douglas-fir (Pseudotsuga menziesii Mirb. (Franco) var. menziesii) trees. Specific needle area decreased with increasing needle age. The magnitude and rate of change in SNA with needle age were greatest for lower-crown branches, and decreased toward the middle- and upper-crown branches. For all branches, there was no difference between regular and epicormic shoots in the relationship between SNA and needle age. Specific needle area decreased with increasing distance from branch base, and this relationship was significant for the majority of needle age classes of the upper- and middle-crown branches. In the lowercrown branches, SNA did not vary with distance from branch base for the majority of needle age classes. For all branches, there was no difference between regular and epicormic shoots in the relationship between SNA and distance from branch base for the majority of needle age classes. These results indicate that renewal of foliage by epicormic shoot production maintains needle quality. Branch SNA increased linearly with decreasing height in the crown at a mean rate of 0.951 +/- 0.110 cm2 g-1 per vertical meter. Total needle area of branches was estimated from total needle dry weight taking into account within-branch variation in SNA. Analyses of allometric relationships between branch size and foliage amount (needle area and needle dry weight) showed that branch length was a better predictor of foliage amount than branch diameter for old Douglas-fir trees. Total needle dry weight and needle area of the sample trees, estimated from branch length and branch height and taking into account vertical within-crown variation in branch SNA, ranged from 42.4 to 154.2 kg and from 246.2 to 816.0 m2 per tree, respectively.

Plant Leaves↗

PubMed2000

Complexity of the simplest phylogenetic estimation problem.

The maximum-likelihood (ML) solution to a simple phylogenetic estimation problem is obtained analytically The problem is estimation of the rooted tree for three species using binary characters with a symmetrical rate of substitution under the molecular clock. ML estimates of branch lengths and log-likelihood scores are obtained analytically for each of the three rooted binary trees. Estimation of the tree topology is equivalent to partitioning the sample space (space of possible data outcomes) into subspaces, within each of which one of the three binary trees is the ML tree. Distance-based least squares and parsimony-like methods produce essentially the same estimate of the tree topology, although differences exist among methods even under this simple model. This seems to be the simplest case, but has many of the conceptual and statistical complexities involved in phylogeny estimation. The solution to this real phylogeny estimation problem will be useful for studying the problem of significance evaluation.

Likelihood Functions↗