Genic mutation rates in mammals: local similarity,chromosomal heterogeneity,and X-versus-autosome disparity

The reduction of mutation rates on the mammalian X chromosome relative to autosomes is most often explained in the literature as evidence of male-driven evolution. This hypothesis attributes lowered mutation rates on the X chromosome to the fact that this chromosome spends less time in the germline of males than in the germline of females. In contrast to this majority view, two articles argued that the patterns of mutation rates across chromosomes are inconsistent with male-driven evolution. One article reported a 40% reduction in synonymous substitution rates (Ks) for X-linked genes relative to autosomes in the mouse-rat lineage. The authors argued that this reduction is too dramatic to be explained by male-driven evolution and concluded that selection has systematically reduced mutation rate on the X chromosome to a level optimal for this male-hemizygous chromosome. More recently, a second article found that chromosomal mutation rates in both the human-mouse and mouse-rat lineages were so heterogeneous that the X chromosome was not an outlier. Here again, the authors argued that this is at odds with male-driven evolution and suggested that selection has modulated chromosomal mutation rates to locally optimal levels, thus extending the argument of the first mentioned article to include autosomes. Here, we reexamine these conclusions using mouse-rat and human-mouse coding-region data. We find a more modest reduction of Ks on the X chromosome, but our results contradict the finding that the X chromosome is not distinct from autosomes. Multiple statistical tests show that Ks rates on the X chromosome differ systematically from the autosomes in both lineages. We conclude that the moderate reduction of mutation rate on the X chromosome of both lineages is consistent with male-driven evolution; however, the large variance in mutation rates across chromosomes suggests that mutation rates are affected by additional factors besides male-driven evolution. Investigation of mutation rates by synteny reveals that synteny blocks, rather than entire chromosomes, might represent the unit of mutation rate variation.     

Maximum-likelihood phylogenetic analysis under a covarion-like model

Here, a model allowing covarion-like evolution of DNA sequences is introduced. In contrast to standard representation of the distribution of evolutionary rates, this model allows the site-specific rate to vary between lineages. This is achieved by adding as few as two parameters to the widely used among-site rate variation model, namely, (1) the proportion of sites undergoing rate changes and (2) the rate of rate change. This model is implemented in the likelihood framework, allowing parameter estimation, comparison of models, and tree reconstruction. An application to ribosomal RNA sequences suggests that covarions (i.e., site-specific rate changes) play an important role in the evolution of these molecules. Neglecting them results in a severe underestimate of the variance of rates across sites. It has, however, little influence on the estimation of ancestral G+C contents obtained from a nonhomogeneous model, or on the resulting inferences about the evolution of thermophyly. This theoretical effort should be useful for the study of protein adaptation, which presumably proceeds in a typical covarion-like manner.     

Exploring the correlations between sequence evolution rate and phenotypic divergence across the Mammalian tree provides insights into adaptive evolution

Sequence evolution behaves in a relatively consistent manner, leading to one of the fundamental paradigms in biology, the existence of a ??molecular clock??. The molecular clock can be distilled to the concept of accumulation of substitutions, through time yielding a stable rate from which we can estimate lineage divergence. Over the last 50?years, evolutionary biologists have obtained an in-depth understanding of this clock??s nuances. It has been fine-tuned by taking into account the vast heterogeneity in rates across lineages and genes, leading to ??relaxed?? molecular clock methods for timetree reconstruction. Sequence rate varies with life history traits including body size, generation time and metabolic rate, and we review recent studies on this topic. However, few studies have explicitly examined correlates between molecular evolution and morphological evolution. The patterns observed across diverse lineages suggest that rates of molecular and morphological evolution are largely decoupled. We discuss how identifying the molecular mechanisms behind rapid functional radiations are central to understanding evolution. The vast functional divergence within mammalian lineages that have relatively ??slow?? sequence evolution refutes the hypotheses that pulses in diversification yielding major phenotypic change are the result of steady accumulation of substitutions. Patterns rather suggest phenotypic divergence is likely caused by regulatory alterations mediated through mechanisms such as insertions/deletions in functional regions. These can rapidly arise and sweep to fixation faster than predicted from a lineage??s sequence neutral substitution rate, enabling species to leapfrog between phenotypic ??islands??. We suggest research directions that could illuminate mechanisms behind the functional diversity we see today.     

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.     

Universal Pacemaker of Genome Evolution

A fundamental observation of comparative genomics is that the distribution of evolution rates across the complete sets of orthologous genes in pairs of related genomes remains virtually unchanged throughout the evolution of life, from bacteria to mammals. The most straightforward explanation for the conservation of this distribution appears to be that the relative evolution rates of all genes remain nearly constant, or in other words, that evolutionary rates of different genes are strongly correlated within each evolving genome. This correlation could be explained by a model that we denoted Universal PaceMaker (UPM) of genome evolution. The UPM model posits that the rate of evolution changes synchronously across genome-wide sets of genes in all evolving lineages. Alternatively, however, the correlation between the evolutionary rates of genes could be a simple consequence of molecular clock (MC). We sought to differentiate between the MC and UPM models by fitting thousands of phylogenetic trees for bacterial and archaeal genes to supertrees that reflect the dominant trend of vertical descent in the evolution of archaea and bacteria and that were constrained according to the two models. The goodness of fit for the UPM model was better than the fit for the MC model, with overwhelming statistical significance, although similarly to the MC, the UPM is strongly overdispersed. Thus, the results of this analysis reveal a universal, genome-wide pacemaker of evolution that could have been in operation throughout the history of life.     

Effects of models of rate evolution on estimation of divergence dates with special reference to the metazoan 18S ribosomal RNA phylogeny

The molecular clock, i.e., constancy of the rate of evolution over time, is commonly assumed in estimating divergence dates. However, this assumption is often violated and has drastic effects on date estimation. Recently, a number of attempts have been made to relax the clock assumption. One approach is to use maximum likelihood, which assigns rates to branches and allows the estimation of both rates and times. An alternative is the Bayes approach, which models the change of the rate over time. A number of models of rate change have been proposed. We have extended and evaluated models of rate evolution, i.e., the lognormal and its recent variant, along with the gamma, the exponential, and the Ornstein-Uhlenbeck processes. These models were first applied to a small hominoid data set, where an empirical Bayes approach was used to estimate the hyperparameters that measure the amount of rate variation. Estimation of divergence times was sensitive to these hyperparameters, especially when the assumed model is close to the clock assumption. The rate and date estimates varied little from model to model, although the posterior Bayes factor indicated the Ornstein-Uhlenbeck process outperformed the other models. To demonstrate the importance of allowing for rate change across lineages, this general approach was used to analyze a larger data set consisting of the 18S ribosomal RNA gene of 39 metazoan species. We obtained date estimates consistent with paleontological records, the deepest split within the group being about 560 million years ago. Estimates of the rates were in accordance with the Cambrian explosion hypothesis and suggested some more recent lineage-specific bursts of evolution.     

Unequal evolutionary rates between annual and perennial lineages of checker mallows (Sidalcea, Malvaceae): evidence from 18S-26S rDNA internal and external transcribed spacers

Heterogeneous DNA substitution rates were found in the 18S-26S nuclear ribosomal DNA internal transcribed spacer (ITS) and external transcribed spacer (ETS) regions of Sidalcea (Malvaceae), a putatively young genus of annuals and perennials. The majority of comparisons revealed that the annual species had significantly higher molecular evolutionary rates than the perennials, whereas rates were consistently homogenous between obligate annual species. These findings led us to conclude that generation time or possibly another biological factor distinguishing annuals and perennials has influenced rates of molecular evolution in SIDALCEA: The congruence of relative-rate test results across both spacer regions reinforced the association between life history and rate of rDNA evolution across lineages of checker mallows. Evolutionary rate variation within perennials mainly involved three basally divergent lineages. The faster rate in one lineage, Sidalcea stipularis, compared with other perennials may be the result of genetic drift in the only known, small, population. The other two basally divergent lineages had slower evolutionary rates compared with the remaining perennials; possible explanations for these differences include rate-reducing effects of a suffrutescent (rather than herbaceous) habit and seed dormancy.     

Reconciling extreme branch length differences: decoupling time and rate through the evolutionary history of filmy ferns

The rate of molecular evolution is not constant across the Tree of Life. Characterizing rate discrepancies and evaluating the relative roles of time and rate along branches through the past are both critical to a full understanding of evolutionary history. In this study, we explore the interactions of time and rate in filmy ferns (Hymenophyllaceae), a lineage with extreme branch length differences between the two major clades. We test for the presence of significant rate discrepancies within and between these clades, and we separate time and rate across the filmy fern phylogeny to simultaneously yield an evolutionary time scale of filmy fern diversification and reconstructions of ancestral rates of molecular evolution. Our results indicate that the branch length disparity observed between the major lineages of filmy ferns is indeed due to a significant difference in molecular evolutionary rate. The estimation of divergence times reveals that the timing of crown group diversification was not concurrent for the two lineages, and the reconstruction of ancestral rates of molecular evolution points to a substantial rate deceleration in one of the clades. Further analysis suggests that this may be due to a genome-wide deceleration in the rate of nucleotide substitution.     

A compound poisson process for relaxing the molecular clock

The molecular clock hypothesis remains an important conceptual and analytical tool in evolutionary biology despite the repeated observation that the clock hypothesis does not perfectly explain observed DNA sequence variation. We introduce a parametric model that relaxes the molecular clock by allowing rates to vary across lineages according to a compound Poisson process. Events of substitution rate change are placed onto a phylogenetic tree according to a Poisson process. When an event of substitution rate change occurs, the current rate of substitution is modified by a gamma-distributed random variable. Parameters of the model can be estimated using Bayesian inference. We use Markov chain Monte Carlo integration to evaluate the posterior probability distribution because the posterior probability involves high dimensional integrals and summations. Specifically, we use the Metropolis-Hastings-Green algorithm with 11 different move types to evaluate the posterior distribution. We demonstrate the method by analyzing a complete mtDNA sequence data set from 23 mammals. The model presented here has several potential advantages over other models that have been proposed to relax the clock because it is parametric and does not assume that rates change only at speciation events. This model should prove useful for estimating divergence times when substitution rates vary across lineages.     

Angiosperm divergence times: the effect of genes, codon positions, and time constraints

An understanding of the evolution of modern terrestrial ecosystems requires an understanding of the dynamics associated with angiosperm evolution, including the timing of their origin and diversification into their extraordinary present-day diversity. Molecular estimates of angiosperm age have varied widely, and many substantially predate the Early Cretaceous fossil appearance of the group. In this study, the effect of different genes, codon positions, and chronological constraints on node ages are examined on divergence time estimates across seed plants, with a special focus on angiosperms. Penalized likelihood was used to estimate divergence times on a phylogenetic hypothesis for seed plants derived from Bayesian analysis, with branch lengths estimated with maximum likelihood. The plastid genes atpB, psaA, psbB, and rbcL were used individually and in combination, using first and second, third, and the three codon positions, including and excluding age constraints on 20 nodes derived from a critical examination of the land-plant fossil record. The optimal level of rate smoothing according to each unconstrained and constrained dataset was obtained with penalized likelihood. Tests for a molecular clock revealed significantly unclocklike rates in all datasets. Addition of fossil constraints resulted in even greater departures from constancy. Consistently with significant deviations from a clock, estimated optimal smoothing values were low, but a strict correlation between rate heterogeneity and optimal smoothing value was not found. Age estimates for nodes across the phylogeny varied, sometimes substantially, with gene and codon position. Nevertheless, estimates based on the four concatenated genes are very similar to the mean of the four individual gene estimates. For any given node, unconstrained age estimates are more variable than constrained estimates and are frequently younger than well-substantiated fossil members of the clade. Constrained estimates of ages of clades are older than unconstrained estimates and oldest fossil representatives, sometimes substantially so. Angiosperm age estimates decreased as rate smoothing increased. Whereas the range of unconstrained angiosperm age estimates spans the fossil age of the clade, the range of constrained estimates is narrower (and older) than the earliest angiosperm fossils. Results unambiguously indicate the relevance of constraints in reducing the variability of ages derived from different partitions of the data and diminishing the effect of the smoothing parameter. Constrained optimizations of divergence times and substitution rates across the phylogeny suggest appreciably different evolutionary dynamics for angiosperms and for gymnosperms. Whereas the gymnosperm crown group originated shortly after the origin of seed plants, a long time elapsed before the origin of crown group angiosperms. Although absolute age estimates of angiosperms and angiosperm clades are older than their earliest fossils, the estimated pace of phylogenetic diversification largely agrees with the rapid appearance of angiosperm lineages in stratigraphic sequences.     

Calibrating the avian molecular clock

Molecular clocks are widely used to date phylogenetic events, yet evidence supporting the rate constancy of molecular clocks through time and across taxonomic lineages is weak. Here, we present 90 candidate avian clock calibrations obtained from fossils and biogeographical events. Cross-validation techniques were used to identify and discard 16 inconsistent calibration points. Molecular evolution occurred in an approximately clock-like manner through time for the remaining 74 calibrations of the mitochondrial gene, cytochrome b . A molecular rate of approximately 2.1% (± 0.1%, 95% confidence interval) was maintained over a 12-million-year interval and across most of 12 taxonomic orders. Minor but significant variance in rates occurred across lineages but was not explained by differences in generation time, body size or latitudinal distribution as previously suggested.     

Population size and molecular evolution on islands

The nearly neutral theory predicts that the rate and pattern of molecular evolution will be influenced by effective population size (Ne), because in small populations more slightly deleterious mutations are expected to drift to fixation. This important prediction has not been widely empirically tested, largely because of the difficulty of comparing rates of molecular evolution in sufficient numbers of independent lineages which differ only in Ne. Island endemic species provide an ideal test of the effect of Ne on molecular evolution because species restricted to islands frequently have smaller Ne than closely related mainland species, and island endemics have arisen from mainland lineages many times in a wide range of taxa. We collated a dataset of 70 phylogenetically independent comparisons between island and mainland taxa, including vertebrates, invertebrates and plants, from 19 different island groups. The rate of molecular evolution in these lineages was estimated by maximum likelihood using two measures: overall substitution rate and the ratio of non-synonymous to synonymous substitution rates. We show that island lineages have significantly higher ratios of non-synonymous to synonymous substitution rates than mainland lineages, as predicted by the nearly neutral theory, although overall substitution rates do not differ significantly.     

Developmental life history is associated with variation in rates of climatic niche evolution in a salamander adaptive radiation*

Rates of climatic niche evolution vary widely across the tree of life and are strongly associated with rates of diversification among clades. However, why the climatic niche evolves more rapidly in some clades than others remains unclear. Variation in life history traits often plays a key role in determining the environmental conditions under which species can survive, and therefore, could impact the rate at which lineages can expand in available climatic niche space. Here, we explore the relationships among life-history variation, climatic niche breadth, and rates of climatic niche evolution. We reconstruct a phylogeny for the genus Desmognathus, an adaptive radiation of salamanders distributed across eastern North America, based on nuclear and mitochondrial genes. Using this phylogeny, we estimate rates of climatic niche evolution for species with long, short, and no aquatic larval stage. Rates of climatic niche evolution are unrelated to the mean climatic niche breadth of species with different life histories. Instead, we find that the evolution of a short larval period promotes greater exploration of climatic space, leading to increased rates of climatic niche evolution across species having this trait. We propose that morphological and physiological differences associated with variation in larval stage length underlie the heterogeneous ability of lineages to explore climatic niche space. Rapid rates of climatic niche evolution among species with short larval periods were an important dimension of the clade's adaptive radiation and likely contributed to the rapid rate of lineage accumulation following the evolution of an aquatic life history in this clade. Our results show how variation in a key life-history trait can constrain or promote divergence of the climatic niche, leading to variation in rates of climatic niche evolution among species.     

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.      

Temporal Patterns of Diversification across Global Cichlid Biodiversity (Acanthomorpha: Cichlidae)

The contrasting distribution of species diversity across the major lineages of cichlids makes them an ideal group for investigating macroevolutionary processes. In this study, we investigate whether different rates of diversification may explain the disparity in species richness across cichlid lineages globally. We present the most taxonomically robust time-calibrated hypothesis of cichlid evolutionary relationships to date. We then utilize this temporal framework to investigate whether both species-rich and depauperate lineages are associated with rapid shifts in diversification rates and if exceptional species richness can be explained by clade age alone. A single significant rapid rate shift increase is detected within the evolutionary history of the African subfamily Pseudocrenilabrinae, which includes the haplochromins of the East African Great Lakes. Several lineages from the subfamilies Pseudocrenilabrinae (Australotilapiini, Oreochromini) and Cichlinae (Heroini) exhibit exceptional species richness given their clade age, a net rate of diversification, and relative rates of extinction, indicating that clade age alone is not a sufficient explanation for their increased diversity. Our results indicate that the Neotropical Cichlinae includes lineages that have not experienced a significant rapid burst in diversification when compared to certain African lineages (rift lake). Neotropical cichlids have remained comparatively understudied with regard to macroevolutionary patterns relative to African lineages, and our results indicate that of Neotropical lineages, the tribe Heroini may have an elevated rate of diversification in contrast to other Neotropical cichlids. These findings provide insight into our understanding of the diversification patterns across taxonomically disparate lineages in this diverse clade of freshwater fishes and one of the most species-rich families of vertebrates.     

Mitochondrial rate variation among lineages of passerine birds

The order Passeriformes comprises the majority of extant avian species. Analyses of molecular data have provided important insights into the evolution of this diverse order. However, molecular estimates of the evolutionary and demographic timescales of passerine species have been hindered by a lack of reliable calibrations. This has led to a reliance on the application of standard substitution rates to mitochondrial DNA data, particularly rates estimated from analyses of the gene encoding cytochrome b (CYTB). To investigate patterns of rate variation across passerine lineages, we used a Bayesian phylogenetic approach to analyse the protein‐coding genes of 183 mitochondrial genomes. We found that the most commonly used mitochondrial marker, CYTB, has low variation in rates across passerine lineages. This lends support to its widespread use as a molecular clock in birds. However, we also found that the patterns of among‐lineage rate variation in CYTB are only weakly related to the evolutionary rate of the mitochondrial genome as a whole. Our analyses confirmed the presence of mutational saturation at third codon positions across the protein‐coding genes of the mitochondrial genome, reinforcing the view that these sites should be excluded in studies of deep passerine relationships. The results of our analyses have provided information that will be useful for molecular‐clock studies of passerine evolution.     

Contrasting the efficacy of selection on the X and autosomes in Drosophila

To investigate the relative efficacy of both positive and purifying natural selection on the X chromosome and the autosomes in Drosophila, we compared rates and patterns of molecular evolution between these chromosome sets using the newly available alignments of orthologous genes from 12 species. Parameters that may influence the relative X versus autosomal substitution rates include the relative effective population sizes, the male and female germline mutation rates, the distribution of allelic effects on fitness, and the degree of dominance of novel mutations. Our analysis reveals that codon usage bias is consistently greater for X-linked genes, suggesting that purifying selection consistently has greater efficacy on the X chromosome than on the autosomes across the Drosophila phylogeny. However, our results are less consistent with respect to the efficacy of positive selection, with only some lineages showing a higher substitution rate on the X chromosome. This suggests that either the distribution of selective effects of mutations or other relevant parameters are sufficiently variable across species to tip the balance in different ways in individual lineages. These data suggest that rates of substitution are not solely governed by adaptive evolution. This genome-wide analysis provides a clear picture that the efficacy of selection varies intragenomically and that this effect is markedly more consistent across the phylogeny in the case of purifying selection. Our results also suggest that simple models that predict systematic differences in rates of evolution between the X and the autosomes can only be made to be compatible with these Drosophila data if the relevant population genetic parameters that drive substitution rates differ among species and chromosomal contexts.     

Artifactual phylogenies caused by correlated distribution of substitution rates among sites and lineages: the good, the bad, and the ugly

Despite the advances in understanding molecular evolution, current phylogenetic methods barely take account of a fraction of the complexity of evolution. We are chiefly constrained by our incomplete knowledge of molecular evolutionary processes and the limits of computational power. These limitations lead to the establishment of either biologically simplistic models that rarely account for a fraction of the complexity involved or overfitting models that add little resolution to the problem. Such oversimplified models may lead us to assign high confidence to an incorrect tree (inconsistency). Rate-across-site (RAS) models are commonly used evolutionary models in phylogenetic studies. These account for heterogeneity in the evolutionary rates among sites but do not account for changing within-site rates across lineages (heterotachy). If heterotachy is common, using RAS models may lead to systematic errors in tree inference. In this work we show possible misleading effects in tree inference when the assumption of constant within-site rates across lineages is violated using maximum likelihood. Using a simulation study, we explore the ways in which gamma stationary models can lead to wrong topology or to deceptive bootstrap support values when the within-site rates change across lineages. More precisely, we show that different degrees of heterotachy mislead phylogenetic inference when the model assumed is stationary. Finally, we propose a geometry-based approach to visualize and to test for the possible existence of bias due to heterotachy.     

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.