Elevated rates of nonsynonymous substitution in island birds

Slightly deleterious mutations are expected to fix at relatively higher rates in small populations than in large populations. Support for this prediction of the nearly-neutral theory of molecular evolution comes from many cases in which lineages inferred to differ in long-term average population size have different rates of nonsynonymous substitution. However, in most of these cases, the lineages differ in many other ways as well, leaving open the possibility that some factor other than population size might have caused the difference in substitution rates. We compared synonymous and nonsynonymous substitutions in the mitochondrial cyt b and ND2 genes of nine closely related island and mainland lineages of ducks and doves. We assumed that island taxa had smaller average population sizes than those of their mainland sister taxa for most of the time since they were established. In all nine cases, more nonsynonymous substitutions occurred on the island branch, but synonymous substitutions showed no significant bias. As in previous comparisons of this kind, the lineages with smaller populations might differ in other respects that tend to increase rates of nonsynonymous substitution, but here such differences are expected to be slight owing to the relatively recent origins of the island taxa. An examination of changes to apparently "preferred" and "unpreferred" synonymous codons revealed no consistent difference between island and mainland lineages.     

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.     

Rates of molecular evolution in nuclear genes of east Mediterranean scorpions

Scorpions of the genus Mesobuthus represent a useful terrestrial model system for studying molecular evolution. They are distributed on several Aegean islands and the adjacent mainland, they are believed to have low rates of dispersal, and evolutionary divergence dates of taxa are available based on biogeographic events that separated islands from each other and the mainland. Here, we present data on polymorphism and synonymous (Ks) and non-synonymous (Ka) substitution rates for nine nuclear protein-coding genes of two east Mediterranean scorpion species, Mesobuthus gibbosus and M. cyprius (Buthidae). Levels of polymorphism tend to be lower in populations from islands (mean nucleotide diversity pi = 0.0071 +/- 0.0028) than in mainland populations (mean pi = 0.0201 +/- 0.0085). By using linear regression of genetic divergence versus isolation time, we estimate Ks to be 3.17 +/- 1.54 per (site x 10(9) years), and Ka to be 0.39 +/- 0.94 per (site x 10(9) years). These estimates for both Ks and Ka are considerably lower than for many other invertebrates, such as Drosophila, and may be attributed to scorpions' mammal-like generation times (approximately 2 years) and low metabolic rates. Phylogenetic analysis using maximum likelihood revealed a phylogeny that is congruent with that expected based on biogeographic events and in which divergences at synonymous sites are proportional to the dates that the taxa are believed to have split. Tests of equality of branch lengths for the Cyprus and Crete lineages revealed that Ks-estimates are about the same in both lineages, as expected from the biogeographic events that separated the islands, but Ka was increased in the Cyprus lineage compared to the Cretan lineage.     

Effective population size and the rate and pattern of nucleotide substitutions

Both the overall rate of nucleotide substitution and the relative proportions of synonymous and non-synonymous substitutions are predicted to vary between species that differ in effective population size (N_e). Our understanding of the genetic processes underlying these lineage-specific differences in molecular evolution is still developing. Empirical analyses indicate that variation in substitution rates and patterns caused by differences in N_e is often substantial, however, and must be accounted for in analyses of molecular evolution.     

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.     

百合目叶绿体rbcL基因核苷酸替代速率检测

分析了百合目主要类群叶绿体中编码核酮糖1,5二磷酸羧化氧化酶大亚基rbcL基因的42条序列,使用RRTree相对速率检测方法,详细研究rbcL基因在百合目7科间同义替代速率和非同义替代速率的变化.相对速率检测显示:百合目内秋水仙科(Colchicaceae)的同义替代速率和非同义替代速率均最快,金梅草科(Campynemat-aceae)同义替代速率最慢,百合科(Liliaceae)的非同义替代速率最慢,但在百合目各科间,无论同义替代速率还是非同义替代速率差异均不显著.     

Rates of nucleotide substitution in Cornaceae (Cornales)-Pattern of variation and underlying causal factors

Identifying causes of genetic divergence is a central goal in evolutionary biology. Although rates of nucleotide substitution vary among taxa and among genes, the causes of this variation tend to be poorly understood. In the present study, we examined the rate and pattern of molecular evolution for five DNA regions over a phylogeny of Cornus, the single genus of Cornaceae. To identify evolutionary mechanisms underlying the molecular variation, we employed Bayesian methods to estimate divergence times and to infer how absolute rates of synonymous and nonsynonymous substitutions and their ratios change over time. We found that the rates vary among genes, lineages, and through time, and differences in mutation rates, selection type and intensity, and possibly genetic drift all contributed to the variation of substitution rates observed among the major lineages of Cornus. We applied independent contrast analysis to explore whether speciation rates are linked to rates of molecular evolution. The results showed no relationships for individual genes, but suggested a possible localized link between species richness and rate of nonsynonymous nucleotide substitution for the combined cpDNA regions. Furthermore, we detected a positive correlation between rates of molecular evolution and morphological change in Cornus. This was particularly pronounced in the dwarf dogwood lineage, in which genome-wide acceleration in both molecular and morphological evolution has likely occurred.     

Synonymous and nonsynonymous rate variation in nuclear genes of mammals

A maximum likelihood approach was used to estimate the synonymous and nonsynonymous substitution rates in 48 nuclear genes from primates, artiodactyls, and rodents. A codon-substitution model was assumed, which accounts for the genetic code structure, transition/transversion bias, and base frequency biases at codon positions. Likelihood ratio tests were applied to test the constancy of nonsynonymous to synonymous rate ratios among branches (evolutionary lineages). It is found that at 22 of the 48 nuclear loci examined, the nonsynonymous/synonymous rate ratio varies significantly across branches of the tree. The result provides strong evidence against a strictly neutral model of molecular evolution. Our likelihood estimates of synonymous and nonsynonymous rates differ considerably from previous results obtained from approximate pairwise sequence comparisons. The differences between the methods are explored by detailed analyses of data from several genes. Transition/transversion rate bias and codon frequency biases are found to have significant effects on the estimation of synonymous and nonsynonymous rates, and approximate methods do not adequately account for those factors. The likelihood approach is preferable, even for pairwise sequence comparison, because more-realistic models about the mutation and substitution processes can be incorporated in the analysis. Received: 17 May 1997 / Accepted: 28 September 1997     

Synonymous and nonsynonymous substitutions in mammalian genes and the nearly neutral theory

The nearly neutral theory of molecular evolution predicts larger generation-time effects for synonymous than for nonsynonymous substitutions. This prediction is tested using the sequences of 49 single-copy genes by calculating the average and variance of synonymous and nonsynonymous substitutions in mammalian star phylogenies (rodentia, artiodactyla, and primates). The average pattern of the 49 genes supports the prediction of the nearly neutral theory, with some notable exceptions.The nearly neutral theory also predicts that the variance of the evolutionary rate is larger than the value predicted by the completely neutral theory. This prediction is tested by examining the dispersion index (ratio of the variance to the mean), which is positively correlated with the average substitution number. After weighting by the lineage effects, this correlation almost disappears for nonsynonymous substitutions, but not quite so for synonymous substitutions. After weighting, the dispersion indices of both synonymous and nonsynonymous substitutions still exceed values expected under the simple Poisson process. The results indicate that both the systematic bias in evolutionary rate among the lineages and the episodic type of rate variation are contributing to the large variance. The former is more significant to synonymous substitutions than to nonsynonymous substitutions. Isochore evolution may be similar to synonymous substitutions. The rate and pattern found here are consistent with the nearly neutral theory, such that the relative contributions of drift and selection differ between the two types of substitutions. The results are also consistent with Gillespie's episodic selection theory.     

Function-specific accelerations in rates of sequence evolution suggest predictable epistatic responses to reduced effective population size

Changes in effective population size impinge on patterns of molecular evolution. Notably, slightly deleterious mutations are more likely to drift to fixation in smaller populations, which should typically also lead to an overall acceleration in the rates of evolution. This prediction has been validated empirically for several endosymbiont and island taxa. Here, we first show that rate accelerations are also evident in bacterial pathogens whose recent shifts in virulence make them prime candidates for reduced effective population size: Bacillus anthracis, Bordetella parapertussis, Mycobacterium leprae, Salmonella enterica typhi, Shigella spp., and Yersinia pestis. Using closely related genomes to analyze substitution rate dynamics across six phylogenetically independent bacterial clades, we demonstrate that relative rates of coding sequence evolution are biased according to gene functional category. Notably, genes that buffer against slightly deleterious mutations, such as chaperones, experience stronger rate accelerations than other functional classes at both nonsynonymous and synonymous sites. Although theory predicts altered evolutionary dynamics for buffer loci in the face of accumulating deleterious mutations, to observe even stronger rate accelerations is surprising. We suggest that buffer loci experience elevated substitution rates because the accumulation of deleterious mutations in the remainder of the genome favors compensatory substitutions in trans. Critically, the hyper-acceleration is evident across phylogenetically independent clades, supporting the hypothesis that reductions in effective population size predictably induce epistatic responses in genes that buffer against slightly deleterious mutations.     

Sequence evolution in bacterial endosymbionts having extreme base compositions.

A major limitation on ability to reconstruct bacterial evolution is the lack of dated ancestors that might be used to evaluate and calibrate molecular clocks. Vertically transmitted symbionts that have cospeciated with animal hosts offer a firm basis for calibrating sequence evolution in bacteria, since fossils of the hosts can be used to date divergence events. Sequences for a functionally diverse set of genes have been obtained for bacterial endosymbionts (Buchnera) from two pairs of aphid host species, each pair diverging 50-70 MYA. Using these dates and estimated numbers of Buchnera generations per year, we calculated rates of base substitution for neutral and selected sites of protein-coding genes and overall rates for rRNA genes. Buchnera shows homogeneity among loci with regard to synonymous rate. The Buchnera synonymous rate is about twice that for low-codon-bias genes of Escherichia coli-Salmonella typhimurium on an absolute timescale, and fourfold higher on a generational timescale. Nonsynonymous substitutions show a greater rate disparity in favor of Buchnera, a result consistent with a genomewide decrease in selection efficiency in Buchnera. Ratios of synonymous to nonsynonymous substitutions differ for the two pairs of Buchnera, indicating that selection efficiency varies among lineages. Like numerous other intracellular bacteria, such as Rickettsia and Wolbachia, Buchnera has accumulated amino acids with codons rich in A or T. Phylogenetic reconstruction of amino acid replacements indicates that replacements yielding increased A + T predominated early in the evolution of Buchnera, with the trend slowing or stopping during the last 50 Myr. This suggests that base composition in Buchnera has approached a limit enforced by selective constraint acting on protein function.     

Apparent Selection Intensity for the Cytochrome Oxidase Subunit I Gene Varies with Mode of Reproduction in Echinoderms

When most amino acid substitutions in protein-coding genes are slightly deleterious rather than selectively neutral, life history differences can potentially modify the effective population size or the selective regime, resulting in altered ratios of non-synonymous to synonymous substitutions among taxa. We studied substitution patterns for the mitochondrial cytochrome oxidase subunit I (COI) gene in a sea star genus (Leptasterias spp.) with an obligate brood-protecting mode of reproduction and small-scale population genetic subdivision, and compared the results to available COI sequences in nine other genera of echinoderms with pelagic larvae: three sea stars, five sea urchins and one brittle star. We predicted that this life history difference would be associated with differences in the ratio of non-synonymous (dN) to synonymous (dS) substitution rates. Leptasterias had a significantly greater dN/dS ratio (both between species and within species), a significantly smaller transition/transversion rate ratio, and a significantly lower average nucleotide diversity within species, than did the non-brooding genera. Other explanations for the results, such as altered mutation rates or selective sweeps, were not supported by the data analysis. These findings highlight the potential influence of reproductive traits and other life history factors on patterns of nucleotide substitution within and between species.     

Chromosomal rearrangements are associated with higher rates of molecular evolution in mammals

Evolutionary rates are not uniformly distributed across the genome. Knowledge about the biological causes of this observation is still incomplete, but its exploration has provided valuable insight into the genomical, historical and demographical variables that influence rates of genetic divergence. Recent studies suggest a possible association between chromosomal rearrangements and regions of greater divergence, but evidence is limited and contradictory. Here, we test the hypothesis of a relationship between chromosomal rearrangements and higher rates of molecular evolution by studying the genomic distribution of divergence between 12,000 human-mouse orthologous genes. Our results clearly show that genes located in genomic regions that have been highly rearranged between the two species present higher rates of synonymous (0.7686 vs. 0.7076) and non-synonymous substitution (0.1014 vs. 0.0871), and that synonymous substitution rates are higher in genes close to the breakpoints of individual rearrangements. The many potential causes of such striking are discussed, particularly in the light of speciation models suggesting that chromosomal rearrangements may have contributed to some of the speciation processes along the human and mouse lineages. Still, there are other possible causes and further research is needed to properly explore them.     

Life history influences rates of climatic niche evolution in flowering plants

Across angiosperms, variable rates of molecular substitution are linked with life-history attributes associated with woody and herbaceous growth forms. As the number of generations per unit time is correlated with molecular substitution rates, it is expected that rates of phenotypic evolution would also be influenced by differences in generation times. Here, we make the first broad-scale comparison of growth-form-dependent rates of niche evolution. We examined the climatic niches of species on large time-calibrated phylogenies of five angiosperm clades and found that woody lineages have accumulated fewer changes per million years in climatic niche space than related herbaceous lineages. Also, climate space explored by woody lineages is consistently smaller than sister lineages composed mainly of herbaceous taxa. This pattern is probably linked to differences in the rate of climatic niche evolution. These results have implications for niche conservatism; in particular, the role of niche conservatism in the distribution of plant biodiversity. The consistent differences in the rate of climatic niche evolution also emphasize the need to incorporate models of phenotypic evolution that allow for rate heterogeneity when examining large datasets.     

Synonymous rates at the RpII215 gene of Drosophila: variation among species and across the coding region

The region encompassing the RpII215 gene that encodes the largest component of the RNA polymerase II complex (1889 amino acids) has been sequenced in Drosophila subobscura, D. madeirensis, D. guanche, and D. pseudoobscura. Nonsynonymous divergence estimates (Ka) indicate that this gene has a very low rate of amino acid replacements. Given its low Ka and constitutive expression, synonymous substitution rates are, however, unexpectedly high. Sequence comparisons have allowed the molecular clock hypothesis to be tested. D. guanche is an insular species and it is therefore expected to have a reduced effective size relative to D. subobscura. The significantly higher rate of synonymous substitutions detected in the D. guanche lineage could be explained if synonymous mutations behave as nearly neutral. Significant departure from the molecular clock hypothesis for synonymous and nonsynonymous substitutions was detected when comparing the D. subobscura, D. pseudoobscura, and D. melanogaster lineages. Codon bias and synonymous divergence between D. subobscura and D. melanogaster were negatively correlated across the RpII215 coding region, which indicates that selection coefficients for synonymous mutations vary across the gene. The C-terminal domain (CTD) of the RpII215 protein is structurally and functionally differentiated from the rest of the protein. Synonymous substitution rates were significantly different in both regions, which strongly indicates that synonymous mutations in the CTD and in the non-CTD regions are under detectably different selection coefficients.     

Effects of metabolic rate on protein evolution

Since the modern evolutionary synthesis was first proposed early in the twentieth century, attention has focused on assessing the relative contribution of mutation versus natural selection on protein evolution. Here we test a model that yields general quantitative predictions on rates of protein evolution by combining principles of individual energetics with Kimura's neutral theory. The model successfully predicts much of the heterogeneity in rates of protein evolution for diverse eukaryotes (i.e. fishes, amphibians, reptiles, birds, mammals) from different thermal environments. Data also show that the ratio of non-synonymous to synonymous nucleotide substitution is independent of body size, and thus presumably of effective population size. These findings indicate that rates of protein evolution are largely controlled by mutation rates, which in turn are strongly influenced by individual metabolic rate.     

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.     

A comprehensive phylogeny of Neurospora reveals a link between reproductive mode and molecular evolution in fungi

The filamentous ascomycete genus Neurospora encompasses taxa with a wide range of reproductive modes. Sexual reproduction in this genus can be divided into three major modes; heterothallism (self-incompatibility), homothallism (self-compatibility) and pseudohomothallism (partial self-compatibility). In addition to the sexual pathway, most of the heterothallic taxa propagate with morphologically distinct, vegetative dissemination propagules (macroconidia), while this feature is undetected in the majority of the homothallic taxa. In this study, we used sequence information of seven nuclear gene loci from 43 taxa (295 of the possible 301 locus-by-taxon combinations) to create a phylogeny of Neurospora. The results suggest that transitions in reproductive mode have occurred at multiple times within this group of fungi. Although a homothallic ancestor would imply fewer switches in reproductive mode, we argue that the ancestor of Neurospora was likely heterothallic and that homothallism has evolved independently at least six times in the evolutionary history of the genus. Furthermore, the two pseudohomothallic taxa of Neurospora (N. tetrasperma and N. tetraspora) represent two independent origins of pseudohomothallism. Likelihood ratio tests of substitution rates among branches in the phylogeny indicate that reproductive mode is an important factor driving genome evolution in Neurospora. First, an increased level of non-synonymous/synonymous substitutions in branches delineating homothallic taxa was found, suggesting a reduced efficiency of purifying selection in these taxa. Furthermore, elevated nucleotide substitution rates were found in heterothallic, conidia-producing, lineages as compared to the homothallic non-conidiating lineages. The latter finding is likely due to the presence of conidia, i.e., a higher rate of mitotic divisions inducing mutations, and/or that the homothallic taxa have evolved a lower mutation rate to avoid genomic degeneration.     

Correlates of substitution rate variation in mammalian protein-coding sequences

Background  

Rates of molecular evolution in different lineages can vary widely, and some of this variation might be predictable from aspects of species' biology. Investigating such predictable rate variation can help us to understand the causes of molecular evolution, and could also help to improve molecular dating methods. Here we present a comprehensive study of the life history correlates of substitution rate variation across the mammals, comparing results for mitochondrial and nuclear loci, and for synonymous and non-synonymous sites. We use phylogenetic comparative methods, refined to take into account the special nature of substitution rate data. Particular attention is paid to the widespread correlations between the components of mammalian life history, which can complicate the interpretation of results.     

Statistical methods for detecting molecular adaptation

The past few years have seen the development of powerful statistical methods for detecting adaptive molecular evolution. These methods compare synonymous and nonsynonymous substitution rates in protein-coding genes, and regard a nonsynonymous rate elevated above the synonymous rate as evidence for darwinian selection. Numerous cases of molecular adaptation are being identified in various systems from viruses to humans. Although previous analyses averaging rates over sites and time have little power, recent methods designed to detect positive selection at individual sites and lineages have been successful. Here, we summarize recent statistical methods for detecting molecular adaptation, and discuss their limitations and possible improvements.