Tree imbalance causes a bias in phylogenetic estimation of evolutionary timescales using heterochronous sequences

Phylogenetic estimation of evolutionary timescales has become routine in biology, forming the basis of a wide range of evolutionary and ecological studies. However, there are various sources of bias that can affect these estimates. We investigated whether tree imbalance, a property that is commonly observed in phylogenetic trees, can lead to reduced accuracy or precision of phylogenetic timescale estimates. We analysed simulated data sets with calibrations at internal nodes and at the tips, taking into consideration different calibration schemes and levels of tree imbalance. We also investigated the effect of tree imbalance on two empirical data sets: mitogenomes from primates and serial samples of the African swine fever virus. In analyses calibrated using dated, heterochronous tips, we found that tree imbalance had a detrimental impact on precision and produced a bias in which the overall timescale was underestimated. A pronounced effect was observed in analyses with shallow calibrations. The greatest decreases in accuracy usually occurred in the age estimates for medium and deep nodes of the tree. In contrast, analyses calibrated at internal nodes did not display a reduction in estimation accuracy or precision due to tree imbalance. Our results suggest that molecular‐clock analyses can be improved by increasing taxon sampling, with the specific aims of including deeper calibrations, breaking up long branches and reducing tree imbalance.     

Phylogenetic estimates of diversification rate are affected by molecular rate variation

Molecular phylogenies are increasingly being used to investigate the patterns and mechanisms of macroevolution. In particular, node heights in a phylogeny can be used to detect changes in rates of diversification over time. Such analyses rest on the assumption that node heights in a phylogeny represent the timing of diversification events, which in turn rests on the assumption that evolutionary time can be accurately predicted from DNA sequence divergence. But there are many influences on the rate of molecular evolution, which might also influence node heights in molecular phylogenies, and thus affect estimates of diversification rate. In particular, a growing number of studies have revealed an association between the net diversification rate estimated from phylogenies and the rate of molecular evolution. Such an association might, by influencing the relative position of node heights, systematically bias estimates of diversification time. We simulated the evolution of DNA sequences under several scenarios where rates of diversification and molecular evolution vary through time, including models where diversification and molecular evolutionary rates are linked. We show that commonly used methods, including metric‐based, likelihood and Bayesian approaches, can have a low power to identify changes in diversification rate when molecular substitution rates vary. Furthermore, the association between the rates of speciation and molecular evolution rate can cause the signature of a slowdown or speedup in speciation rates to be lost or misidentified. These results suggest that the multiple sources of variation in molecular evolutionary rates need to be considered when inferring macroevolutionary processes from phylogenies.     

Testing the molecular clock: molecular and paleontological estimates of divergence times in the Echinoidea (Echinodermata)

The phylogenetic relationships of 46 echinoids, with representatives from 13 of the 14 ordinal-level clades and about 70% of extant families commonly recognized, have been established from 3 genes (3,226 alignable bases) and 119 morphological characters. Morphological and molecular estimates are similar enough to be considered suboptimal estimates of one another, and the combined data provide a tree that, when calibrated against the fossil record, provides paleontological estimates of divergence times and completeness of their fossil record. The order of branching on the cladogram largely agrees with the stratigraphic order of first occurrences and implies that their fossil record is more than 85% complete at family level and at a resolution of 5-Myr time intervals. Molecular estimates of divergence times derived from applying both molecular clock and relaxed molecular clock models are concordant with estimates based on the fossil record in up to 70% of cases, with most concordant results obtained using Sanderson's semiparametric penalized likelihood method and a logarithmic-penalty function. There are 3 regions of the tree where molecular and fossil estimates of divergence time consistently disagree. Comparison with results obtained when molecular divergence dates are estimated from the combined (morphology + gene) tree suggests that errors in phylogenetic reconstruction explain only one of these. In another region the error most likely lies with the paleontological estimates because taxa in this region are demonstrated to have a very poor fossil record. In the third case, morphological and paleontological evidence is much stronger, and the topology for this part of the molecular tree differs from that derived from the combined data. Here the cause of the mismatch is unclear but could be methodological, arising from marked inequality of molecular rates. Overall, the level of agreement reached between these different data and methodological approaches leads us to believe that careful application of likelihood and Bayesian methods to molecular data provides realistic divergence time estimates in the majority of cases (almost 80% in this specific example), thus providing a remarkably well-calibrated phylogeny of a character-rich clade of ubiquitous marine benthic invertebrates.     

The low evolutionary rate of human T-cell lymphotropic virus type-1 confirmed by analysis of vertical transmission chains

The evolutionary rate of the human T-cell lymphotropic virus type-1 (HTLV-1) is considered to be very low, in strong contrast to the related human retrovirus HIV. However, current estimates of the HTLV-1 rate rely on the anthropological calibration of phylogenies using assumed dates of human migration events. To obtain an independent rate estimate, we analyzed two variable regions of the HTLV-1 genome (LTR and env) from eight infected families. Remarkable genetic stability was observed, as only two mutations in LTR (756 bp) and three mutations in env (522 bp) occurred within the 16 vertical transmission chains, including one ambiguous position in each region. The evolutionary rate in HTLV-1 was then calculated using a maximum-likelihood approach that used the highest and lowest possible times of HTLV-1 shared ancestry, given the known transmission histories. The rates for the LTR and env regions were 9.58 x 10(-8)-1.25 x 10(-5) and 7.84 x 10(-7) -2.33 x 10(-5)nucleotide substitutions per site per year, respectively. A more precise estimate was obtained for the combined LTR-env data set, which was 7.06 x 10(-7)-1.38 x 10(-5)substitutions per site per year. We also note an interesting correlation between the occurrence of mutations in HTLV-1 and the age of the individual infected.     

Variable evolutionary rates in the molecular evolution of mammalian growth hormones

In mammals pituitary growth hormone (GH) shows a slow basal rate of evolution (0.22 ± 0.03 × 10^–9 substitutions/amino acid site/year) which appears to have increased by at least 25–50-fold on two occasions, during the evolution of primates (to at least 10.8 ± 1.3 X 10^–9 substitutions/amino acid site/year) and artiodactyl ruminants (to at least 5.6 ± 1.3 X 10^–9 substitutions/amino acid site/year). That these rate increases are real, and not due to inadvertent comparison of nonorthologous genes, was established by showing that features of the GH gene sequences that are not expressed as mature hormone do not show corresponding changes in evolutionary rate. Thus, analysis of nonsynonymous substitutions in the coding sequence for the mature protein confirmed the rate increases seen in the primate and ruminant GHs, but analysis of nonsynonymous substitutions in the signal peptide sequence, synonymous substitutions in the coding sequence for signal peptide or mature protein, and 5 and 3 untranslated sequences showed no statistically significant changes in evolutionary rate. Evidence that the increases in evolutionary rate are probably due to positive selection is provided by the observation that in the cases of both ruminant and primate GHs the periods of rapid evolution were followed by a return to a slow rate similar to the basal rate seen in other mammalian GHs. Differences between the biological properties of GHs have been identified which may relate to these periods of rapid adaptive molecular evolution. On the basis of sequence data currently available (but excluding rodent GHs which show an intermediate rate, the basis of which is not clear) for most (90%) of evolutionary time mammalian GHs have been in the slow phase of evolution, with possibly most of the few amino acid substitutions that have occurred being neutral in nature. But most (80%) of the amino acid substitutions that have been introduced into GH during the course of mammalian evolution have been accepted during the rapid phases and were adaptive in nature.     

Avian influenza virus exhibits rapid evolutionary dynamics

Influenza A viruses from wild aquatic birds, their natural reservoir species, are thought to have reached a form of stasis, characterized by low rates of evolutionary change. We tested this hypothesis by estimating rates of nucleotide substitution in a diverse array of avian influenza viruses (AIV) and allowing for rate variation among lineages. The rates observed were extremely high, at >10(-3) substitutions per site, per year, with little difference among wild and domestic host species or viral subtypes and were similar to those seen in mammalian influenza A viruses. Influenza A virus therefore exhibits rapid evolutionary dynamics across its host range, consistent with a high background mutation rate and rapid replication. Using the same approach, we also estimated that the common ancestors of the hemagglutinin and neuraminidase sequences of AIV arose within the last 3,000 years, with most intrasubtype diversity emerging within the last 100 years and suggestive of a continual selective turnover.     

Primate phylogeny, evolutionary rate variations, and divergence times: a contribution from the nuclear gene IRBP

The first third (ca. 1200 bp) of exon 1 of the nuclear gene encoding the interstitial retinoid-binding protein (IRBP) has been sequenced for 12 representative primates belonging to Lemuriformes, Lorisiformes, Tarsiiformes, Platyrrhini, and Catarrhini, and combined with available data (13 other primates, 11 nonprimate placentals, and 2 marsupials). Phylogenetic analyses using maximum likelihood on nucleotides and amino acids robustly support the monophyly of primates, Strepsirrhini, Lemuriformes, Lorisiformes, Anthropoidea, Catarrhini, and Platyrrhini. It is interesting to note that 1) Tarsiidae grouped with Anthropoidea, and the support for this node depends on the molecular characters considered; 2) Cheirogaleidae grouped within Lemuriformes; and 3) Daubentonia was the sister group of all other Lemuriformes. Study of the IRBP evolutionary rate shows a high heterogeneity within placentals and also within primates. Maximum likelihood local molecular clocks were assigned to three clades displaying significantly contrasted evolutionary rates. Paenungulata were shown to evolve 2.5-3 times faster than Perissodactyla and Lemuriformes. Six independent calibration points were used to estimate splitting ages of the main primate clades, and their compatibility was evaluated. Divergence ages were obtained for the following crown groups: 13.8-14.2 MY for Lorisiformes, 26.5-27.2 MY for Lemuroidea, 39.6-40.7 MY for Lemuriformes, 45.4-46.7 MY for Strepsirrhini, and 56.7-58.4 MY for Haplorrhini. The incompatibility between some paleontological and molecular estimates may reflect the incompleteness of the placental fossil record, and/or indicate that the variable IRBP evolutionary rates are not fully accommodated by local molecular clocks.     

Neutral theory, phylogenies, and the relationship between phenotypic change and evolutionary rates

The neutral theory of molecular evolution predicts that rates of phenotypic change are largely independent from genotypic change. A recent study by Bromham et al. (2002) confirmed this expectation, finding no evidence for correlated phenotypic and molecular evolutionary rates in animals. We reevaluate this hypothesis, sampling at different taxonomic levels in plants and animals, using Bayesian inference to reconstruct phylogenetic trees and estimate rates of molecular evolution. We use independent contrasts in branch lengths to maximize the information extracted from each of the trees and nodal posterior probabilities to assess the influence of phylogenetic error. Our results indicate that in vascular plants between 2% and 11% of the variation in phenotypic rates of change can be explained by the rate of genotypic change. These results may be explained by the idea that processes that affect general evolutionary rates, such as body size, may also be expected to influence rates of morphological change.     

Complex evolution of orthologous and paralogous decarboxylase genes

The decarboxylases are involved in neurotransmitter synthesis in animals, and in pathways of secondary metabolism in plants. Different decarboxylase proteins are characterized for their different substrate specificities, but are encoded by homologous genes. We study, within a maximum-likelihood framework, the evolutionary relationships among dopa decarboxylase (Ddc), histidine decarboxylase (Hdc) and alpha-methyldopa hypersensitive (amd) in animals, and tryptophan decarboxylase (Wdc) and tyrosine decarboxylase (Ydc) in plants. The evolutionary rates are heterogeneous. There are differences between paralogous genes in the same lineages: 4.13 x 10(-10) nucleotide substitutions per site per year in mammalian Ddc vs. 1.95 in Hdc; between orthologous genes in different lineages, 7.62 in dipteran Ddc vs. 4.13 in mammalian Ddc; and very large temporal variations in some lineages, from 3.7 up to 54.9 in the Drosophila Ddc lineage. Our results are inconsistent with the molecular clock hypothesis.     

On inconsistency of the neighbor-joining, least squares, and minimum evolution estimation when substitution processes are incorrectly modeled

Using analytical methods, we show that under a variety of model misspecifications, Neighbor-Joining, minimum evolution, and least squares estimation procedures are statistically inconsistent. Failure to correctly account for differing rates-across-sites processes, failure to correctly model rate matrix parameters, and failure to adjust for parallel rates-across-sites changes (a rates-across-subtrees process) are all shown to lead to a "long branch attraction" form of inconsistency. In addition, failure to account for rates-across-sites processes is also shown to result in underestimation of evolutionary distances for a wide variety of substitution models, generalizing an earlier analytical result for the Jukes-Cantor model reported in Golding and a similar bias result for the GTR or REV model in Kelly and Rice (1996). Although standard rates-across-sites models can be employed in many of these cases to restore consistency, current models cannot account for other kinds of misspecification. We examine an idealized but biologically relevant case, where parallel changes in rates at sites across subtrees is shown to give rise to inconsistency. This changing rates-across-subtrees type model misspecification cannot be adjusted for with conventional methods or without carefully considering the rate variation in the larger tree. The results are presented for four-taxon trees, but the expectation is that they have implications for larger trees as well. To illustrate this, a simulated 42-taxon example is given in which the microsporidia, an enigmatic group of eukaryotes, are incorrectly placed at the archaebacteria-eukaryotes split because of incorrectly specified pairwise distances. The analytical nature of the results lend insight into the reasons that long branch attraction tends to be a common form of inconsistency and reasons that other forms of inconsistency like "long branches repel" can arise in some settings. In many of the cases of inconsistency presented, a particular incorrect topology is estimated with probability converging to one, the implication being that measures of uncertainty like bootstrap support will be unable to detect that there is a problem with the estimation. The focus is on distance methods, but previous simulation results suggest that the zones of inconsistency for distance methods contain the zones of inconsistency for maximum likelihood methods as well.     

Phylogeny,biogeography, and morphological ancestral character reconstruction in the Mediterranean genus Fumana (Cistaceae)

Fumana is a diverse genus of the Cistaceae family, consisting of 21 currently accepted species. In this study, nuclear (ITS) and plastid (matK, trnT‐L) molecular markers were used to reconstruct the phylogeny and to estimate divergence times, including 19 species of Fumana. Phylogenetic analyses (Bayesian Inference, Maximum Parsimony and Maximum Likelihood) confirmed the monophyly of Fumana and did not support the infrageneric divisions previously established. The results support four main clades that group species that differ in vegetative and reproductive characters. Given the impossibility to define morphological characters common to all species within the clades, our proposal is to reject infrageneric divisions. Molecular dating and ancestral area analyses provide evidence for a Miocene diversification of the genus in the north‐western Mediterranean. Ancestral state reconstructions revealed ancestral character states for some traits related to xeric and arid habitats, suggesting a preadaptation to the Mediterranean climate.     

Why do species vary in their rate of molecular evolution?

Despite hopes that the processes of molecular evolution would be simple, clock-like and essentially universal, variation in the rate of molecular evolution is manifest at all levels of biological organization. Furthermore, it has become clear that rate variation has a systematic component: rate of molecular evolution can vary consistently with species body size, population dynamics, lifestyle and location. This suggests that the rate of molecular evolution should be considered part of life-history variation between species, which must be taken into account when interpreting DNA sequence differences between lineages. Uncovering the causes and correlates of rate variation may allow the development of new biologically motivated models of molecular evolution that may improve bioinformatic and phylogenetic analyses.     

The polyphyletic origin of laboratory inbred mice and their rate of evolution

The analysis of various genetic components in the standard inbred laboratory strains of mice shows that several taxonomic units of the complex species Mus musculus have contributed to their genetic background. These laboratory lines are often taken as archetypes of the mammalian genome and since their genealogy is known for around 80 years they have been used to estimate mutation rates and various other evolutionary parameters. The knowledge of their origins is important if one is to avoid erroneous interpretations. Do they possess haplotypes that could have existed in natural populations of M. musculus domesticus or are they the fruit of artificial recombinations between divergent genomes? We discuss this in the light of various genetic systems.     

Mitochondrial evolutionary rate and speciation in termites: data on European Reticulitermes taxa (Isoptera, Rhinotermitidae)

Summary. The rate of mitochondrial DNA evolution and the speciation pattern in relation to glacial periods are tested in the European taxa of the eusocial genus Reticulitermes. The linearized tree obtained from cytochrome oxidase II sequences and a geological event calibration shows a substitution rate 100-fold higher than that usually applied for insect mitochondrial DNA. An accelerated rate of evolution has also been observed in social Vespidae (Hymenoptera); we therefore suggest the involvement of eusociality in mediating gene pool drift. The role of the last ice age in speciation pattern of Reticulitermes taxa is supported by molecular data, but a four refugia model better explains genetic diversity, phyletic relationships and present-day distribution of these termites.Received 30 March 2004; revised 29 July and 15 November 2004; accepted 1 December 2004.     

Extensive evolutionary rate variation in floral color determining genes in the genus Ipomoea

Flower color is an important adaptive trait in many plant species because it determines reproductive success through differential attractiveness to insect pollinators. The genus Ipomoea is a pan tropically distributed plant genus characterized by showy flowers that often differ in color among closely related species. Flower color is determined primarily by products of the anthocyanin biosynthesis pathway. To determine whether flower color evolution among members of the genus is driven by adaptive molecular evolution of the anthocyanin pathway genes, we analyzed data from 10 genes sequenced from 19 species of Ipomoea . Six protein-coding genes from the anthocyanin pathway were drawn for evolutionary analysis in addition to three genes from the unrelated leucine biosynthesis pathway and one MADS box regulatory gene for comparison. The analyses provided: (i) no convincing evidence for positive selection on anthocyanin pathway structural genes, or on the other sampled genes, despite shifts in flower color among species included in the sample; (ii) pathway position correlated weakly with estimates of the intensity of evolutionary constraint on the anthocyanin pathway enzyme coding genes; and (iii) there was substantial gene-specific heterogeneity in the rates of synonymous site evolution. Synonymous rate heterogeneity does not appear to be accounted for by codon bias or local contextual or compositional sequence differences, leading us to implicate heterogeneous rates of mutation among genes as the most probable cause of synonymous rate heterogeneity.     

Mode and tempo in the evolution of socio-political organization: reconciling 'Darwinian' and 'Spencerian' evolutionary approaches in anthropology

Traditional investigations of the evolution of human social and political institutions trace their ancestry back to nineteenth century social scientists such as Herbert Spencer, and have concentrated on the increase in socio-political complexity over time. More recent studies of cultural evolution have been explicitly informed by Darwinian evolutionary theory and focus on the transmission of cultural traits between individuals. These two approaches to investigating cultural change are often seen as incompatible. However, we argue that many of the defining features and assumptions of 'Spencerian' cultural evolutionary theory represent testable hypotheses that can and should be tackled within a broader 'Darwinian' framework. In this paper we apply phylogenetic comparative techniques to data from Austronesian-speaking societies of Island South-East Asia and the Pacific to test hypotheses about the mode and tempo of human socio-political evolution. We find support for three ideas often associated with Spencerian cultural evolutionary theory: (i) political organization has evolved through a regular sequence of forms, (ii) increases in hierarchical political complexity have been more common than decreases, and (iii) political organization has co-evolved with the wider presence of hereditary social stratification.     

戊型肝炎病毒遗传多样性和进化机制的研究进展

戊型肝炎病毒(hepatitis E virus, HEV)是一种引发全球急性病毒性肝炎的人兽共患病病原。HEV具有丰富的遗传多样性,不同基因型或基因亚型的流行与地理位置、宿主物种以及防控策略等密切相关。欧洲和美洲HEV流行株为HEV-3,包括3a-i亚型,而亚洲流行株含HEV-3和HEV-4;我国的流行毒株已从HEV-1进化到HEV-4。近年来,研究发现HEV进化的影响机制,包括同义密码子使用模式、氨基酸突变和基因重组等,其中氨基酸突变是病毒持续流行的主要驱动力。因此,本文就HEV的分类、全球流行特征、进化机制等进行综述,以期为戊型肝炎的有效防控以及疫苗开发提供参考。     

Extraordinary conservation of entire chromosomes in insects over long evolutionary periods

Comparison of the genomes of different Drosophila species has shown that six different chromosomes, the so‐called ‘‘Muller elements,” constitute the building blocks for all Drosophila species. Here, we confirm previous results suggesting that this conservation of the Muller elements extends far beyond Drosophila, to at least tephritid fruit flies, thought to have diverged from drosophilids 60–70 mYr ago. Less than 10 percent of genes differ in chromosome location between the two insect groups. Within chromosomes, however, the order is highly scrambled, as expected from the comparison between Drosophila species. The data also support the notion that the sex chromosomes of tephritid flies originated from an ancestor of the dot chromosome 4 of Drosophila. Overall, therefore, no new chromosome has been created for perhaps a billion generations over the two evolutionary lines. This stability at the chromosome level, which appears to extend to all Diptera including mosquitoes, is in stark contrast to other groups such as mammals, birds, fish and plants, in which chromosome numbers and organization vary enormously among species that have diverged over much fewer generations.