Covarion structure in plastid genome evolution: a new statistical test

Covarion models of molecular evolution allow the rate of evolution of a site to vary through time. There are few simple and effective tests for covarion evolution, and consequently, little is known about the presence of covarion processes in molecular evolution. We describe two new tests for covarion evolution and demonstrate with simulations that they perform well under a wide range of conditions. A survey of covarion evolution in sequenced plastid genomes found evidence of covarion drift in at least 26 out of 57 genes. Covarion evolution is most evident in first and second codon positions of the plastid genes, and there is no evidence of covarion evolution in third codon positions. Therefore, the significant covarion tests are likely due to changes in the selective constraints of amino acids. The frequency of covarion evolution within the plastid genome suggests that covarion processes of evolution were important in generating the observed patterns of sequence variation among plastid genomes.     

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.     

What can DNA Tell us About the Cambrian Explosion?

Molecular data is ideal for exploring deep evolutionary historybecause of its universality, stochasticity and abundance. Thesefeatures provide a means of exploring the evolutionary historyof all organisms (including those that do not tend to leavefossils), independently of morphological evolution, and withina statistical framework that allows testing of evolutionaryhypotheses. In particular, molecular data have an importantrole to play in examining hypotheses concerning the tempo andmode of evolution of animal body plans. Examples are given wheremolecular phylogenies have led to a re-examination of some fundamentalassumptions in metazoan evolution, such as the immutabilityof early developmental characters, and the evolvability of bauplancharacters. Molecular data is also providing a new and controversialtimescale for the evolution of animal phyla, pushing the majordivisions of the animal kingdom deep into the Precambrian. Therehave been many reasons to question the accuracy and precisionof molecular date estimates, such as the failure to accountfor lineage-specific rate variation and unreliable estimationof rates of molecular evolution. While these criticisms havebeen largely countered by recent studies, one problem has remaineda challenge: could temporal variation in the rate of molecularevolution, perhaps associated with "explosive" adaptive radiations,cause overestimation of diversification dates? Empirical evidencefor an effect of speciation rate, morphological evolution orecological diversification on rates of molecular evolution isexamined, and the potential for rate-variable methods for moleculardating are discussed.     

Statistical models of the overdispersed molecular clock

The most commonly used statistical model to describe the rate constancy of molecular evolution (molecular clock) is a simple Poisson process in which the variance of the number of amino acid or nucleotide substitutions in a particular gene should be equal to the mean and henceforth the dispersion index, the ratio of the variance to the mean, should be equal to one. Recent sequence data, however, have shown that the substitutional process in molecular evolution is often considerably overdispersed and have called into question the generality of using a simple Poisson process. Several efforts have been made to develop more realistic models of molecular evolution. In this paper, I will show that the spatial (site-specific) variation in the rate of molecular evolution is an improbable cause of the overdispersion and then review various statistical models which take the temporal variation into account. Although these models do not immediately specify what the mechanisms of molecular evolution might be, they do make qualitatively different predictions and give some insight into their inference. One way to distinguish them is suggested. In addition, effects of selected substitutions that presumably occur after a major change in a molecule are quasi-quantitatively examined. It is most likely that the overdispersion of molecular clock is due either to a major molecular reconfiguration (fluctuating neutral space) led by a series of subliminal neutral changes or to selected substitutions fine-tuning a molecule after a major molecular change. Although the latter possibility, of course, violates the simplest neutrality assumption, it would not impair the neutral theory as a whole.     

Selectionism and neutralism in molecular evolution

Charles Darwin proposed that evolution occurs primarily by natural selection, but this view has been controversial from the beginning. Two of the major opposing views have been mutationism and neutralism. Early molecular studies suggested that most amino acid substitutions in proteins are neutral or nearly neutral and the functional change of proteins occurs by a few key amino acid substitutions. This suggestion generated an intense controversy over selectionism and neutralism. This controversy is partially caused by Kimura's definition of neutrality, which was too strict (|2Ns|< or =1). If we define neutral mutations as the mutations that do not change the function of gene products appreciably, many controversies disappear because slightly deleterious and slightly advantageous mutations are engulfed by neutral mutations. The ratio of the rate of nonsynonymous nucleotide substitution to that of synonymous substitution is a useful quantity to study positive Darwinian selection operating at highly variable genetic loci, but it does not necessarily detect adaptively important codons. Previously, multigene families were thought to evolve following the model of concerted evolution, but new evidence indicates that most of them evolve by a birth-and-death process of duplicate genes. It is now clear that most phenotypic characters or genetic systems such as the adaptive immune system in vertebrates are controlled by the interaction of a number of multigene families, which are often evolutionarily related and are subject to birth-and-death evolution. Therefore, it is important to study the mechanisms of gene family interaction for understanding phenotypic evolution. Because gene duplication occurs more or less at random, phenotypic evolution contains some fortuitous elements, though the environmental factors also play an important role. The randomness of phenotypic evolution is qualitatively different from allele frequency changes by random genetic drift. However, there is some similarity between phenotypic and molecular evolution with respect to functional or environmental constraints and evolutionary rate. It appears that mutation (including gene duplication and other DNA changes) is the driving force of evolution at both the genic and the phenotypic levels.     

The first molecular phylogeny of Strepsiptera (Insecta) reveals an early burst of molecular evolution correlated with the transition to endoparasitism

A comprehensive model of evolution requires an understanding of the relationship between selection at the molecular and phenotypic level. We investigate this in Strepsiptera, an order of endoparasitic insects whose evolutionary biology is poorly studied. We present the first molecular phylogeny of Strepsiptera, and use this as a framework to investigate the association between parasitism and molecular evolution. We find evidence of a significant burst in the rate of molecular evolution in the early history of Strepsiptera. The evolution of morphological traits linked to parasitism is significantly correlated with the pattern in molecular rate. The correlated burst in genotypic-phenotypic evolution precedes the main phase of strepsipteran diversification, which is characterised by the return to a low and even molecular rate, and a period of relative morphological stability. These findings suggest that the transition to endoparasitism led to relaxation of selective constraint in the strepsipteran genome. Our results indicate that a parasitic lifestyle can affect the rate of molecular evolution, although other causal life-history traits correlated with parasitism may also play an important role.     

Detecting Site-Specific Biochemical Constraints Through Substitution Mapping

The neutral theory of molecular evolution states that most mutations are deleterious or neutral. It results that the evolutionary rate of a given position in an alignment is a function of the level of constraint acting on this position. Inferring evolutionary rates from a set of aligned sequences is hence a powerful method to detect functionally and/or structurally important positions in a protein. Some positions, however, may be constrained while having a high substitution rate, providing these substitutions do not affect the biochemical property under constraint. Here, I introduce a new evolutionary rate measure accounting for the evolution of specific biochemical properties (e.g., volume, polarity, and charge). I then present a new statistical method based on the comparison of two rate measures: a site is said to be constrained for property X if it shows an unexpectedly high conservation of X knowing its total evolutionary rate. Compared to single-rate methods, the two-rate method offers several advantages: it (i) allows assessment of the significance of the constraint, (ii) provides information on the type of constraint acting on each position, and (iii) detects positions that are not proposed by previous methods. I apply this method to a 200-sequence data set of triosephosphate isomerase and report significant cases of positions constrained for polarity, volume, or charge. The three-dimensional localization of these positions shows that they are of potential interest to the molecular evolutionist and to the biochemist.     

Origin and Evolution of a New Gene Descended from Alcohol Dehydrogenase in Drosophila

Drosophila alcohol dehydrogenase (Adh) is highly conserved in size, organization, and amino acid sequence. Adh-ψ was hypothesized to be a pseudogene derived from an Adh duplication in the repleta group of Drosophila; however, several results from molecular analyses of this gene conflict with currently held notions of molecular evolution. Perhaps the most difficult observations to reconcile with the pseudogene hypothesis are that the hypothetical replacement sites of Adh-ψ evolve only slightly more quickly than replacement sites of closely related, functional Adh genes, and that the replacement sites of the pseudogenes evolve considerably more slowly than neighboring silent sites. The data have been presented as a paradox that challenges our understanding of the mechanisms underlying DNA sequence divergence. Here I show that Adh-ψ is actually a new, functional gene recently descended from an Adh duplication. This descendant recruited ~60 new N-terminal amino acids, is considerably more basic than ADH, and is evolving at a faster rate than Adh. Furthermore, though the descendant is clearly functional, as inferred from molecular evolution and population genetic data, it retains no obvious ADH activity. This probably reflects functional divergence from its Adh ancestor.     

Standardized phylogenetic tree: a reference to discover functional evolution

Functional evolution is often driven by positive natural selection. Although it is thought to be rare in evolution at the molecular level, its effects may be observed as the accelerated evolutionary rates. Therefore one of the effective ways to identify functional evolution is to identify accelerated evolution. Many methods have been developed to test the statistical significance of the accelerated evolutionary rate by comparison with the appropriate reference rate. The rates of synonymous substitution are one of the most useful and popular references, especially for large-scale analyses. On the other hand, these rates are applicable only to a limited evolutionary time period because they saturate quickly--i.e., multiple substitutions happen frequently because of the lower functional constraint. The relative rate test is an alternative method. This technique has an advantage in terms of the saturation effect but is not sufficiently powerful when the evolutionary rate differs considerably among phylogenetic lineages. For the aim to provide a universal reference tree, we propose a method to construct a standardized tree which serves as the reference for accelerated evolutionary rate. The method is based upon multiple molecular phylogenies of single genes with the aim of providing higher reliability. The tree has averaged and normalized branch lengths with standard deviations for statistical neutrality limits. The standard deviation also suggests the reliability level of the branch order. The resulting tree serves as a reference tree for the reliability level of the branch order and the test of evolutionary rate acceleration even when some of the species lineages show an accelerated evolutionary rate for most of their genes due to bottlenecking and other effects.     

The Power of Relative Rates Tests Depends on the Data

One of the most useful features of molecular phylogenetic analyses is the potential for estimating dates of divergence of evolutionary lineages from the DNA of extant species. But lineage-specific variation in rate of molecular evolution complicates molecular dating, because a calibration rate estimated from one lineage may not be an accurate representation of the rate in other lineages. Many molecular dating studies use a ``clock test' to identify and exclude sequences that vary in rate between lineages. However, these clock tests should not be relied upon without a critical examination of their effectiveness at removing rate variable sequences from any given data set, particularly with regard to the sequence length and number of variable sites. As an illustration of this problem we present a power test of a frequently employed triplet relative rates test. We conclude that (1) relative rates tests are unlikely to detect moderate levels of lineage-specific rate variation (where one lineage has a rate of molecular evolution 1.5 to 4.0 times the other) for most commonly used sequences in molecular dating analyses, and (2) this lack of power is likely to result in substantial error in the estimation of dates of divergence. As an example, we show that the well-studied rate difference between murid rodents and great apes will not be detected for many of the sequences used to date the divergence between these two lineages and that this failure to detect rate variation is likely to result in consistent overestimation the date of the rodent–primate split. Received: 9 June 1999 / Accepted: 22 October 1999     

Congruence between comparative morphology and molecular phylogenies: Evolution of the male genital skeletal/muscular system in the subtribe Polyommatina (Lepidoptera,Lycaenidae)

The skeleton and musculature of male genitalia were studied in species of a model butterfly group (subtribe Polyommatina, Lycaenidae). In total, we analyzed 45 species of the tribe Polyommatini most of which were previously used in the molecular phylogenetic study (Talavera et al., 2013). The studied morphological characters were mapped on the molecular trees, which allowed us to reveal trends of morphological changes and to estimate the age of their origin. As a result, chronology of evolution of skeleton and musculature traits was established. It was shown that periods of slow morphological evolution alternated in the subtribe Polyommatina with those of a high rate of origin of new traits. For example, topography of the intravalvar muscles has not changed for 26 MY preserving their initial fan-shaped attachment. The evolution of intravalvar muscles started 10 MYA, proceeded slowly during the first 5 MY, and then accelerated during the last 5 MY resulting in the extensive splitting of the musculature in most monophyletic lineages. Mapping the morphological characters on the phylogeny demonstrated that the rates of skeleton and muscle evolution within the skeleton/musculature apparatus were different. In most cases the intravalvar musculature evolved much faster than the skeleton. The cladistic interpretation of states of morphological traits was found to be consistent with phylogenetic reconstructions based on analysis of multiple molecular markers. Moreover, morphological synapomorphies were found for the lineages Alpherakya + Glabroculus and Aricia + (Alpherakya + Glabroculus), which had low statistical support in molecular phylogenetic analysis. Additionally, in some cases molecular studies helped to reveal trends in the evolution of morphological traits. For example, the unpaired uncus and the compact juxta are not plesiomorphic for Cupidina as previously thought; instead, they were shown to have evolved secondarily within this subtribe.     

Rice yellow mottle virus, an RNA plant virus, evolves as rapidly as most RNA animal viruses

The rate of evolution of an RNA plant virus has never been estimated using temporally spaced sequence data, by contrast to the information available on an increasing range of animal viruses. Accordingly, the evolution rate of Rice yellow mottle virus (RYMV) was calculated from sequences of the coat protein gene of isolates collected from rice over a 40-year period in different parts of Africa. The evolution rate of RYMV was estimated by pairwise distance linear regression on five phylogeographically defined groups comprising a total of 135 isolates. It was further assessed from 253 isolates collected all over Africa by Bayesian coalescent methods under strict and relaxed molecular clock models and under constant size and skyline population genetic models. Consistent estimates of the evolution rate between 4 x 10(-4) and 8 x 10(-4) nucleotides (nt)/site/year were obtained whatever method and model were applied. The synonymous evolution rate was between 8 x 10(-4) and 11 x 10(-4) nt/site/year. The overall and synonymous evolution rates of RYMV were within the range of the rates of 50 RNA animal viruses, below the average but above the distribution median. Experimentally, in host change studies, substitutions accumulated at an even higher rate. The results show that an RNA plant virus such as RYMV evolves as rapidly as most RNA animal viruses. Knowledge of the molecular clock of plant viruses provides methods for testing a wide range of biological hypotheses.     

Can fast early rates reconcile molecular dates with the Cambrian explosion?

Molecular dates consistently place the divergence of major metazoan lineages in the Precambrian, leading to the suggestion that the 'Cambrian explosion' is an artefact of preservation which left earlier forms unrecorded in the fossil record. While criticisms of molecular analyses for failing to deal with variation in the rate of molecular evolution adequately have been countered by analyses which allow both site-to-site and lineage-specific rate variation, no analysis to date has allowed the rates to vary temporally. If the rates of molecular evolution were much higher early in the metazoan radiation, molecular dates could consistently overestimate the divergence times of lineages. Here, we use a new method which uses multiple calibration dates and an empirically determined range of possible substitution rates to place bounds on the basal date of divergence of lineages in order to ask whether faster rates of molecular evolution early in the metazoan radiation could possibly account for the discrepancy between molecular and palaeontological date estimates. We find that allowing basal (interphylum) lineages the fastest observed substitution rate brings the minimum possible divergence date (586 million years ago) to the Vendian period, just before the first multicellular animal fossils, but excludes divergence of the major metazoan lineages in a Cambrian explosion.     

Morphological innovation, diversification and invasion of a new adaptive zone

How ecological opportunity relates to diversification is a central question in evolutionary biology. However, there are few empirical examples of how ecological opportunity and morphological innovation open new adaptive zones, and promote diversification. We analyse data on diet, skull morphology and bite performance, and relate these traits to diversification rates throughout the evolutionary history of an ecologically diverse family of mammals (Chiroptera: Phyllostomidae). We found a significant increase in diversification rate driven by increased speciation at the most recent common ancestor of the predominantly frugivorous subfamily Stenodermatinae. The evolution of diet was associated with skull morphology, and morphology was tightly coupled with biting performance, linking phenotype to new niches through performance. Following the increase in speciation rate, the rate of morphological evolution slowed, while the rate of evolution in diet increased. This pattern suggests that morphology stabilized, and niches within the new adaptive zone of frugivory were filled rapidly, after the evolution of a new cranial phenotype that resulted in a certain level of mechanical efficiency. The tree-wide speciation rate increased non linearly with a more frugivorous diet, and was highest at measures of skull morphology associated with morphological extremes, including the most derived Stenodermatines. These results show that a novel stenodermatine skull phenotype played a central role in the evolution of frugivory and increasing speciation within phyllostomids.     

Does the VP1 gene of foot-and-mouth disease virus behave as a molecular clock?

Summary We have carried out a phylogenetic study of the evolution of the VP1 gene sequence from different serological types and subtypes of foot-and-mouth disease virus (FMDV). The maximum-likelihood method developed by Hasegawa and co-workers (Hasegawa et al. 1985) for the estimation of evolutionary parameters and branching dates has been used to decide between alternative models of evolution: constant versus variable rates. The results obtained indicate that a constant rate model, i.e., a molecular clock, seems to be the most plausible one. However, additional information suggests the possibility that the appearance of serotype CS has been accompanied by an episode of rapid evolution (Villaverde et al. 1991). We discuss the possibility that this evolution of RNA viruses was due to episodic positive Darwinian selection, which would have helped the new variant to escape the immunogenic pressure from the hosts. Offprint requests to: A. Moya     

Sociality and the rate of molecular evolution

The molecular clock does not tick at a uniform rate in all taxa but may be influenced by species characteristics. Eusocial species (those with reproductive division of labor) have been predicted to have faster rates of molecular evolution than their nonsocial relatives because of greatly reduced effective population size; if most individuals in a population are nonreproductive and only one or few queens produce all the offspring, then eusocial animals could have much lower effective population sizes than their solitary relatives, which should increase the rate of substitution of "nearly neutral" mutations. An earlier study reported faster rates in eusocial honeybees and vespid wasps but failed to correct for phylogenetic nonindependence or to distinguish between potential causes of rate variation. Because sociality has evolved independently in many different lineages, it is possible to conduct a more wide-ranging study to test the generality of the relationship. We have conducted a comparative analysis of 25 phylogenetically independent pairs of social lineages and their nonsocial relatives, including bees, wasps, ants, termites, shrimps, and mole rats, using a range of available DNA sequences (mitochondrial and nuclear DNA coding for proteins and RNAs, and nontranslated sequences). By including a wide range of social taxa, we were able to test whether there is a general influence of sociality on rates of molecular evolution and to test specific predictions of the hypothesis: (1) that social species have faster rates because they have reduced effective population sizes; (2) that mitochondrial genes would show a greater effect of sociality than nuclear genes; and (3) that rates of molecular evolution should be correlated with the degree of sociality. We find no consistent pattern in rates of molecular evolution between social and nonsocial lineages and no evidence that mitochondrial genes show faster rates in social taxa. However, we show that the most highly eusocial Hymenoptera do have faster rates than their nonsocial relatives. We also find that social parasites (that utilize the workers from related species to produce their own offspring) have faster rates than their social relatives, which is consistent with an effect of lower effective population size on rate of molecular evolution. Our results illustrate the importance of allowing for phylogenetic nonindependence when conducting investigations of determinants of variation in rate of molecular evolution.     

The genome as a life-history character: why rate of molecular evolution varies between mammal species

DNA sequences evolve at different rates in different species. This rate variation has been most closely examined in mammals, revealing a large number of characteristics that can shape the rate of molecular evolution. Many of these traits are part of the mammalian life-history continuum: species with small body size, rapid generation turnover, high fecundity and short lifespans tend to have faster rates of molecular evolution. In addition, rate of molecular evolution in mammals might be influenced by behaviour (such as mating system), ecological factors (such as range restriction) and evolutionary history (such as diversification rate). I discuss the evidence for these patterns of rate variation, and the possible explanations of these correlations. I also consider the impact of these systematic patterns of rate variation on the reliability of the molecular date estimates that have been used to suggest a Cretaceous radiation of modern mammals, before the final extinction of the dinosaurs.     

Temperature‐dependent mutational robustness can explain faster molecular evolution at warm temperatures,affecting speciation rate and global patterns of species diversity

Distribution of species across the Earth shows strong latitudinal and altitudinal gradients with the number of species decreasing with declining temperatures. While these patterns have been recognized for well over a century, the mechanisms generating and maintaining them have remained elusive. Here, we propose a mechanistic explanation for temperature‐dependent rates of molecular evolution that can influence speciation rates and global biodiversity gradients. Our hypothesis is based on the effects of temperature and temperature‐adaptation on stability of proteins and other catalytic biomolecules. First, due to the nature of physical forces between biomolecules and water, stability of biomolecules is maximal around + 20°C and decreases as temperature either decreases or increases. Second, organisms that have adapted to cold temperatures have evolved especially flexible (but unstable) proteins to facilitate catalytic reactions in cold, where molecular movements slow down. Both these effects should result in mutations being on average more detrimental at cold temperatures (i.e. lower mutational robustness in cold). At high temperatures, destabilizing water–biomolecule interactions, and the need to maintain structures that withstand heat denaturation, should decrease mutational robustness similarly. Decreased mutational robustness at extreme temperatures will slow down molecular evolution, as a larger fraction of new mutations will be removed by selection. Lower mutational robustness may also select for reduced mutation rates, further slowing down the rate of molecular evolution. As speciation requires the evolution of epistatic incompatibilities that prevent gene flow among incipient species, slow rate of molecular evolution at extreme temperatures will directly slow down the rate at which new species arise. The proposed mechanism can thus explain why molecular evolution is faster at warm temperatures, contributing to higher speciation rate and elevated species richness in environments characterized by stable and warm temperatures.     

Evolutionary assembly of the milkweed fauna: cytochrome oxidase I and the age of Tetraopes beetles

The insects that feed on the related plant families Apocynaceae and Asclepiadaceae (here collectively termed "milkweeds") comprise a "component community" of highly specialized, distinctive lineages of species that frequently sequester toxic cardiac glycosides from their host plants for defense against predators and are thus often aposematic, advertising their consequent unpalatability. Such sets of specialized lineages provide opportunities for comparative studies of the rate of adaptation, diversification, and habitat-related effects on molecular evolution. The cerambycid genus Tetraopes is the most diverse of the new world milkweed herbivores and the species are generally host specific, being restricted to single, different species of Asclepias, more often so than most other milkweed insects. Previous work revealed correspondence between the phylogeny of these beetles and that of their hosts. The present study provides analyses of near-complete DNA sequences for Tetraopes and relatives that are used to establish a molecular clock and temporal framework for Tetraopes evolution with their milkweed hosts.     

Emerging genomic and proteomic evidence on relationships among the animal, plant and fungal kingdoms

Sequence-based molecular phylogenies have provided new models of early eukaryotic evolution. This includes the widely accepted hypothesis that animals are related most closely to fungi, and that the two should be grouped together as the Opisthokonta. Although most published phylogenies have supported an opisthokont relationship, a number of genes contain a tree-building signal that clusters animal and green plant sequences, to the exclusion of fungi. The alternative tree-building signal is especially intriguing in light of emerging data from genomic and proteomic studies that indicate striking and potentially synapomorphic similarities between plants and animals. This paper reviews these new lines of evidence, which have yet to be incorporated into models of broad scale eukaryotic evolution.