Simulating and detecting autocorrelation of molecular evolutionary rates among lineages

Evolutionary timescales can be estimated from genetic data using phylogenetic methods based on the molecular clock. To account for molecular rate variation among lineages, a number of relaxed‐clock models have been developed. Some of these models assume that rates vary among lineages in an autocorrelated manner, so that closely related species share similar rates. In contrast, uncorrelated relaxed clocks allow all of the branch‐specific rates to be drawn from a single distribution, without assuming any correlation between rates along neighbouring branches. There is uncertainty about which of these two classes of relaxed‐clock models are more appropriate for biological data. We present an R package, NELSI, that allows the evolution of DNA sequences to be simulated according to a range of clock models. Using data generated by this package, we assessed the ability of two Bayesian phylogenetic methods to distinguish among different relaxed‐clock models and to quantify rate variation among lineages. The results of our analyses show that rate autocorrelation is typically difficult to detect, even when there is complete taxon sampling. This provides a potential explanation for past failures to detect rate autocorrelation in a range of data sets.     

Molecular evolution of the ent-kaurenoic acid oxidase gene in Oryzeae

We surveyed the substitution patterns in the ent-kaurenoic acid oxidase (KAO) gene in 11 species of Oryzeae with an outgroup in the Ehrhartoidaea. The synonymous and non-synonymous substitution rates showed a high positive correlation with each other, but were negatively correlated with codon usage bias and GC content at third codon positions. The substitution rate was heterogenous among lineages. Likelihood-ratio tests showed that the non-synonymous/synonymous rate ratio changed significantly among lineages. Site-specific models provided no evidence for positive selection of particular amino acid sites in any codon of the KAO gene. This finding suggested that the significant rate heterogeneity among some lineages may have been caused by variability in the relaxation of the selective constraint among lineages or by neutral processes.     

Estimating the rate of evolution of the rate of molecular evolution

A simple model for the evolution of the rate of molecular evolution is presented. With a Bayesian approach, this model can serve as the basis for estimating dates of important evolutionary events even in the absence of the assumption of constant rates among evolutionary lineages. The method can be used in conjunction with any of the widely used models for nucleotide substitution or amino acid replacement. It is illustrated by analyzing a data set of rbcL protein sequences.      

Rates of molecular evolution in bacteria are relatively constant despite spore dormancy

Rates of molecular evolution are known to vary considerably among lineages, partially due to differences in life-history traits such as generation time. The generation-time effect has been well documented in some eukaryotes, but its prevalence in prokaryotes is unknown. "Because many species of Firmicute bacteria spend long periods of time as metabolically dormant spores, which could result in fewer DNA substitutions per unit time, they present an excellent system for testing predictions of the molecular clock hypothesis." To test whether spore-forming bacteria evolve more slowly than their non-spore-forming relatives, I used phylogenetic methods to determine if there were differences in rates of amino acid substitution between spore-forming and non-spore-forming lineages of Firmicute bacteria. Although rates of evolution do vary among lineages, I find no evidence for an effect of spore-formation on evolutionary rate and, furthermore, evolutionary rates are similar to those calculated for enteric bacteria. These results support the notion that variation in generation time does not affect evolutionary rates in bacterial lineages.     

Evolutionary distance estimation under heterogeneous substitution pattern among lineages

Most of the sophisticated methods to estimate evolutionary divergence between DNA sequences assume that the two sequences have evolved with the same pattern of nucleotide substitution after their divergence from their most recent common ancestor (homogeneity assumption). If this assumption is violated, the evolutionary distance estimated will be biased, which may result in biased estimates of divergence times and substitution rates, and may lead to erroneous branching patterns in the inferred phylogenies. Here we present a simple modification for existing distance estimation methods to relax the assumption of the substitution pattern homogeneity among lineages when analyzing DNA and protein sequences. Results from computer simulations and empirical data analyses for human and mouse genes are presented to demonstrate that the proposed modification reduces the estimation bias considerably and that the modified method performs much better than the LogDet methods, which do not require the homogeneity assumption in estimating the number of substitutions per site. We also discuss the relationship of the substitution and mutation rate estimates when the substitution pattern is not the same in the lineages leading to the two sequences compared.     

Maximum-likelihood models for combined analyses of multiple sequence data

Models of nucleotide substitution were constructed for combined analyses of heterogeneous sequence data (such as those of multiple genes) from the same set of species. The models account for different aspects of the heterogeneity in the evolutionary process of different genes, such as differences in nucleotide frequencies, in substitution rate bias (for example, the transition/transversion rate bias), and in the extent of rate variation across sites. Model parameters were estimated by maximum likelihood and the likelihood ratio test was used to test hypotheses concerning sequence evolution, such as rate constancy among lineages (the assumption of a molecular clock) and proportionality of branch lengths for different genes. The example data from a segment of the mitochondrial genome of six hominoid species (human, common and pygmy chimpanzees, gorilla, orangutan, and siamang) were analyzed. Nucleotides at the three codon positions in the protein-coding regions and from the tRNA-coding regions were considered heterogeneous data sets. Statistical tests showed that the amount of evolution in the sequence data reflected in the estimated branch lengths can be explained by the codon-position effect and lineage effect of substitution rates. The assumption of a molecular clock could not be rejected when the data were analyzed separately or when the rate variation among sites was ignored. However, significant differences in substitution rate among lineages were found when the data sets were combined and when the rate variation among sites was accounted for in the models. Under the assumption that the orangutan and African apes diverged 13 million years ago, the combined analysis of the sequence data estimated the times for the human-chimpanzee separation and for the separation of the gorilla as 4.3 and 6.8 million years ago, respectively.     

Rate heterogeneity, ancestral character state reconstruction, and the evolution of limb morphology in Lerista (Scincidae, Squamata)

Rates of phenotypic evolution derive from numerous interrelated processes acting at varying spatial and temporal scales and frequently differ substantially among lineages. Although current models employed in reconstructing ancestral character states permit independent rates for distinct types of transition (forward and reverse transitions and transitions between different states), these rates are typically assumed to be identical for all branches in a phylogeny. In this paper, I present a general model of character evolution enabling rate heterogeneity among branches. This model is employed in assessing the extent to which the assumption of uniform transition rates affects reconstructions of ancestral limb morphology in the scincid lizard clade Lerista and, accordingly, the potential for rate variability to mislead inferences of evolutionary patterns. Permitting rate variation among branches significantly improves model fit for both the manus and the pes. A constrained model in which the rate of digit acquisition is assumed to be effectively zero is strongly supported in each case; when compared with a model assuming unconstrained transition rates, this model provides a substantially better fit for the manus and a nearly identical fit for the pes. Ancestral states reconstructed assuming the constrained model imply patterns of limb evolution differing significantly from those implied by reconstructions for uniform-rate models, particularly for the pes; whereas ancestral states for the uniform-rate models consistently entail the reacquisition of pedal digits, those for the model incorporating among-lineage rate heterogeneity imply repeated, unreversed digit loss. These results indicate that the assumption of identical transition rates for all branches in a phylogeny may be inappropriate in modeling the evolution of phenotypic traits and emphasize the need for careful evaluation of phylogenetic tests of Dollo's law.     

Why mammalian lineages respond differently to sexual selection: metabolic rate constrains the evolution of sperm size

The hypothesis that sperm competition should favour increases in sperm size, because it results in faster swimming speeds, has received support from studies on many taxa, but remains contentious for mammals. We suggest that this may be because mammalian lineages respond differently to sexual selection, owing to major differences in body size, which are associated with differences in mass-specific metabolic rate. Recent evidence suggests that cellular metabolic rate also scales with body size, so that small mammals have cells that process energy and resources from the environment at a faster rate. We develop the 'metabolic rate constraint hypothesis' which proposes that low mass-specific metabolic rate among large mammals may limit their ability to respond to sexual selection by increasing sperm size, while this constraint does not exist among small mammals. Here we show that among rodents, which have high mass-specific metabolic rates, sperm size increases under sperm competition, reaching the longest sperm sizes found in eutherian mammals. By contrast, mammalian lineages with large body sizes have small sperm, and while metabolic rate (corrected for body size) influences sperm size, sperm competition levels do not. When all eutherian mammals are analysed jointly, our results suggest that as mass-specific metabolic rate increases, so does maximum sperm size. In addition, species with low mass-specific metabolic rates produce uniformly small sperm, while species with high mass-specific metabolic rates produce a wide range of sperm sizes. These findings support the hypothesis that mass-specific metabolic rates determine the budget available for sperm production: at high levels, sperm size increases in response to sexual selection, while low levels constrain the ability to respond to sexual selection by increasing sperm size. Thus, adaptive and costly traits, such as sperm size, may only evolve under sexual selection when metabolic rate does not constrain cellular budgets.     

Rates of single-copy DNA evolution in herons

DNA-DNA hybridization was used to discover the extent of single-copy DNA similarity among 13 species of herons and one ibis. Genetic distances among taxa were summarized as Tm values in a folded matrix. From this matrix, trees with the same branching pattern were constructed by least squares under one of two assumptions: (1) that sister branches are equal in length and (2) that sister branches are not necessarily equal in length. The residual sums of squares of these trees were compared by F-test to see whether the branches of the tree built under assumption (2) fit the matrix data significantly better than those of the tree built under assumption (1). By this method the existence of different rates of DNA evolution in different heron lineages was established. Bittern single-copy DNA has evolved at a rate approximately 25% faster, and boat-billed heron (Cochearius) and rufescent tiger heron (Tigrisoma lineatum) DNA has evolved approximately 19% slower, than that of day and night herons. It appears that the differences in rates of DNA evolution may increase proportionally with genealogical distance.      

Substitution rates in hepatitis delta virus

Substitution rates were estimated for the coding and noncoding regions of the hepatitis delta virus (HDV). The estimated rates of synonymous substitution in HDV were lower than the rates of substitution at nonsynonymous sites and in the noncoding region. HDV has lower synonymous substitution rates than the hepatitis C virus, though both are RNA viruses. The relatively low rate of synonymous substitution in HDV may be due to a strong preference of G and C nucleotides at third codon positions. Variation in substitution rate among HDV lineages may be correlated with the clinical development of the HDV-induced hepatitis. The phylogenetic tree inferred for 24 HDV strains reveals similarities between lineages isolated from the same geographic region. Correspondence to: W.-H. Li     

Consistent variation in amino-acid substitution rate, despite uniformity of mutation rate: protein evolution in mammals is not neutral

Variation in mutation rate, attributed to differences in both generation time and in metabolic rate, has been invoked under the neutral theory of molecular evolution to account for differences in substitution rate among mammalian lineages. We show that substitution rates at fourfold-degenerate sites and at sites in noncoding regions do not vary between the primate and rodent lineages, implying mutation- rate uniformity. In contrast, the substitution rates at nondegenerate sites vary both within and between lineages. This difference in substitution-rate pattern between the two types of site is incompatible with neutral theory but may result from substitutions occurring by fixation of slightly deleterious mutations. Variation in the rate of protein evolution among mammalian lineages appears to be due more to differences in population fixation rates than to biochemical or physiological differences affecting mutation rates.      

Estimating absolute rates of synonymous and nonsynonymous nucleotide substitution in order to characterize natural selection and date species divergences

The rate of molecular evolution can vary among lineages. Sources of this variation have differential effects on synonymous and nonsynonymous substitution rates. Changes in effective population size or patterns of natural selection will mainly alter nonsynonymous substitution rates. Changes in generation length or mutation rates are likely to have an impact on both synonymous and nonsynonymous substitution rates. By comparing changes in synonymous and nonsynonymous rates, the relative contributions of the driving forces of evolution can be better characterized. Here, we introduce a procedure for estimating the chronological rates of synonymous and nonsynonymous substitutions on the branches of an evolutionary tree. Because the widely used ratio of nonsynonymous and synonymous rates is not designed to detect simultaneous increases or simultaneous decreases in synonymous and nonsynonymous rates, the estimation of these rates rather than their ratio can improve characterization of the evolutionary process. With our Bayesian approach, we analyze cytochrome oxidase subunit I evolution in primates and infer that nonsynonymous rates have a greater tendency to change over time than do synonymous rates. Our analysis of these data also suggests that rates have been positively correlated.     

Rates of genome evolution and branching order from whole genome analysis

Accurate estimation of any phylogeny is important as a framework for evolutionary analysis of form and function at all levels of organization from sequence to whole organism. Using alignments of nonrepetitive components of opossum, human, mouse, rat, and dog genomes we evaluated two alternative tree topologies for eutherian evolution. We show with very high confidence that there is a basal split between rodents (as represented by the mouse and rat) and a branch joining primates (as represented by humans) and carnivores (as represented by dogs), consistent with some but not the most widely accepted mammalian phylogenies. The result was robust to substitution model choice with equivalent inference returned from a spectrum of models ranging from a general time reversible model, a model that treated nucleotides as either purines and pyrimidines, and variants of these that incorporated rate heterogeneity among sites. By determining this particular branching order we are able to show that the rate of molecular evolution is almost identical in rodent and carnivore lineages and that sequences evolve approximately 11%-14% faster in these lineages than in the primate lineage. In addition by applying the chicken as outgroup the analyses suggested that the rate of evolution in all eutherian lineages is approximately 30% slower than in the opossum lineage. This pattern of relative rates is inconsistent with the hypothesis that generation time is an important determinant of substitution rates and, by implication, mutation rates. Possible factors causing rate differences between the lineages include differences in DNA repair and replication enzymology, and shifts in nucleotide pools. Our analysis demonstrates the importance of using multiple sequences from across the genome to estimate phylogeny and relative evolutionary rate in order to reduce the influence of distorting local effects evident even in relatively long sequences.     

Aging asexual lineages and the evolutionary maintenance of sex

Finite populations of asexual and highly selfing species suffer from a reduced efficacy of selection. Such populations are thought to decline in fitness over time due to accumulating slightly deleterious mutations or failing to adapt to changing conditions. These within‐population processes that lead nonrecombining species to extinction may help maintain sex and outcrossing through species level selection. Although inefficient selection is proposed to elevate extinction rates over time, previous models of species selection for sex assumed constant diversification rates. For sex to persist, classic models require that asexual species diversify at rates lower than sexual species; the validity of this requirement is questionable, both conceptually and empirically. We extend past models by allowing asexual lineages to decline in diversification rates as they age, that is nonrecombining lineages “senesce” in diversification rates. At equilibrium, senescing diversification rates maintain sex even when asexual lineages, at young ages, diversify faster than their sexual progenitors. In such cases, the age distribution of asexual lineages contains a peak at intermediate values rather than showing the exponential decline predicted by the classic model. Coexistence requires only that the average rate of diversification in asexuals be lower than that of sexuals.     

The molecular clock ticks regularly in muroid rodents and hamsters

Summary Extensive DNA sequence data are used to compare the rates of nucleotide substitution in the mouse, rat, and hamster lineages. A relative rate test using hamster sequences as references shows that the rates of synonymous and nonsynonymous substitution in the mouse and rat lineages are nearly equal and a test using human sequences as references shows that the rates in the mouse, rat, and hamster lineages are also nearly equal. Under the assumptions that the guinea pig lineage and the myomorph (mouse, rat, and hamster) lineage diverged 70–100 million years (Myr) ago and that the rate of nucleotide substitution has been constant in all these lineages since their divergence, the date of the mouse-rat split is estimated to be between 20 and 29 Myr ago, which is considerably older than the date ( 12 Myr) suggested by available rodent fossils and considerably younger than the date ( 35 Myr) suggested by Wilson and colleagues. The murid-hamster split is estimated to be 1.6 times older than the mouse-rat split.     

Ancient DNA enables timing of the pleistocene origin and holocene expansion of two adélie penguin lineages in antarctica

The timing of divergent events in history is one of the central goals of contemporary evolutionary biology. Such studies are however dependent on accurate evolutionary rates. Recent developments in ancient DNA analysis enable the estimation of more accurate evolutionary rates and therefore more accurate timing of divergence events. Consequently, this leads to a better understanding of changes in populations through time. We use an evolutionary rate calculated from ancient DNA of Adélie penguins (Pygoscelis adeliae) to time divergent events in their history. We report the presence of two distinct and highly variable mitochondrial DNA lineages and track changes in these lineages through space and time. When the ancient DNA and the phylogenetic rates are used to estimate the time of origin of the lineages, two very different estimates resulted. In addition, these same rates provide very different estimates of the time of expansion of these lineages. We suggest that the rate calculated from ancient DNA is more consistent with the glacial history of Antarctica and requires fewer assumptions than does a narrative based on the phylogenetic rate. Finally, we suggest that our study indicates an important new role for ancient DNA studies in the timing of divergent events in history.     

EXTINCTION RATES SHOULD NOT BE ESTIMATED FROM MOLECULAR PHYLOGENIES

Molecular phylogenies contain information about the tempo and mode of species diversification through time. Because extinction leaves a characteristic signature in the shape of molecular phylogenetic trees, many studies have used data from extant taxa only to infer extinction rates. This is a promising approach for the large number of taxa for which extinction rates cannot be estimated from the fossil record. Here, I explore the consequences of violating a common assumption made by studies of extinction from phylogenetic data. I show that when diversification rates vary among lineages, simple estimators based on the birth–death process are unable to recover true extinction rates. This is problematic for phylogenetic trees with complete taxon sampling as well as for the simpler case of clades with known age and species richness. Given the ubiquity of variation in diversification rates among lineages and clades, these results suggest that extinction rates should not be estimated in the absence of fossil data.     

Patterns and relative rates of nucleotide and insertion/deletion evolution at six chloroplast intergenic regions in new world species of the Lecythidaceae

Insertions and deletions (indels) in chloroplast noncoding regions are common genetic markers to estimate population structure and gene flow, although relatively little is known about indel evolution among recently diverged lineages such as within plant families. Because indel events tend to occur nonrandomly along DNA sequences, recurrent mutations may generate homoplasy for indel haplotypes. This is a potential problem for population studies, because indel haplotypes may be shared among populations after recurrent mutation as well as gene flow. Furthermore, indel haplotypes may differ in fitness and therefore be subject to natural selection detectable as rate heterogeneity among lineages. Such selection could contribute to the spatial patterning of cpDNA haplotypes, greatly complicating the interpretation of cpDNA population structure. This study examined both nucleotide and indel cpDNA variation and divergence at six noncoding regions (psbB-psbH, atpB-rbcL, trnL-trnH, rpl20-5'rps12, trnS-trnG, and trnH-psbA) in 16 individuals from eight species in the Lecythidaceae and a Sapotaceae outgroup. We described patterns of cpDNA changes, assessed the level of indel homoplasy, and tested for rate heterogeneity among lineages and regions. Although regression analysis of branch lengths suggested some degree of indel homoplasy among the most divergent lineages, there was little evidence for indel homoplasy within the Lecythidaceae. Likelihood ratio tests applied to the entire phylogenetic tree revealed a consistent pattern rejecting a molecular clock. Tajima's 1D and 2D tests revealed two taxa with consistent rate heterogeneity, one showing relatively more and one relatively fewer changes than other taxa. In general, nucleotide changes showed more evidence of rate heterogeneity than did indel changes. The rate of evolution was highly variable among the six cpDNA regions examined, with the trnS-trnG and trnH-psbA regions showing as much as 10% and 15% divergence within the Lecythidaceae. Deviations from rate homogeneity in the two taxa were constant across cpDNA regions, consistent with lineage-specific rates of evolution rather than cpDNA region-specific natural selection. There is no evidence that indels are more likely than nucleotide changes to experience homoplasy within the Lecythidaceae. These results support a neutral interpretation of cpDNA indel and nucleotide variation in population studies within species such as Corythophora alta.     

Rates of gene rearrangement and nucleotide substitution are correlated in the mitochondrial genomes of insects

A number of studies indicated that lineages of animals with high rates of mitochondrial (mt) gene rearrangement might have high rates of mt nucleotide substitution. We chose the hemipteroid assemblage and the Insecta to test the idea that rates of mt gene rearrangement and mt nucleotide substitution are correlated. For this purpose, we sequenced the mt genome of a lepidopsocid from the Psocoptera, the only order of hemipteroid insects for which an entire mtDNA sequence is not available. The mt genome of this lepidopsocid is circular, 16,924 bp long, and contains 37 genes and a putative control region; seven tRNA genes and a protein-coding gene in this genome have changed positions relative to the ancestral arrangement of mt genes of insects. We then compared the relative rates of nucleotide substitution among species from each of the four orders of hemipteroid insects and among the 20 insects whose mt genomes have been sequenced entirely. All comparisons among the hemipteroid insects showed that species with higher rates of gene rearrangement also had significantly higher rates of nucleotide substitution statistically than did species with lower rates of gene rearrangement. In comparisons among the 20 insects, where the mt genomes of the two species differed by more than five breakpoints, the more rearranged species always had a significantly higher rate of nucleotide substitution than the less rearranged species. However, in comparisons where the mt genomes of two species differed by five or less breakpoints, the more rearranged species did not always have a significantly higher rate of nucleotide substitution than the less rearranged species. We tested the statistical significance of the correlation between the rates of mt gene rearrangement and mt nucleotide substitution with nine pairs of insects that were phylogenetically independent from one another. We found that the correlation was positive and statistically significant (R2 = 0.73, P = 0.01; Rs = 0.67, P < 0.05). We propose that increased rates of nucleotide substitution may lead to increased rates of gene rearrangement in the mt genomes of insects.     

Mutational meltdown or controlled chain reaction: The dynamics of rapid plastome evolution in the hyperdiversity of Poaceae

The study of genomic structural evolution associated with accelerated evolutionary rates that result in avoidance of meltdown and increase biodiversity is becoming ever more possible as the number of available plastomes increases. To more comprehensively analyze rate heterogeneity among monocots and within Poaceae, we sequenced plastomes from four Poaceae species, combined them with publicly available data from ~200 plastomes, and conducted comparative analyses to quantify the pattern of rate heterogeneity between different lineages, functional groups, and periods of evolutionary time. We compared structural differences across the Poaceae to quantify how changes in plastome size correspond to different genomic subunits and the evolution of IR–SC junction boundaries. The substitution rates among ancestral Poaceae were inferred to be exceptionally rapid compared to other monocots but slowed after divergence into extant lineages, which could not be sufficiently explained by positive selection. As such, rapid rates in the ancestral lineage leading to Poaceae might be more closely linked to large-scale structural changes like the loss of ycf1 and ycf2. The total increase in plastome size across Poaceae was positively correlated with the total length of intergenic spacers, tandem repeats, and dispersed repeats as well as large single copy, and inverted repeats (IRs). The continuous evolution of IR–SC junction boundaries was asynchronous with sizes of total genome and subunits across Poaceae. Future work is needed to better understand what factors in ancestral Poaceae evolved to harness such rapid rates of plastome evolution, avoid a mutational meltdown, and escape the stagnation of strong purifying selection as well as if these factors could be utilized to synthetically control rates.