Overdispersion of the molecular clock: temporal variation of gene-specific substitution rates in Drosophila

Simple models of molecular evolution assume that sequences evolve by a Poisson process in which nucleotide or amino acid substitutions occur as rare independent events. In these models, the expected ratio of the variance to the mean of substitution counts equals 1, and substitution processes with a ratio greater than 1 are called overdispersed. Comparing the genomes of 10 closely related species of Drosophila, we extend earlier evidence for overdispersion in amino acid replacements as well as in four-fold synonymous substitutions. The observed deviation from the Poisson expectation can be described as a linear function of the rate at which substitutions occur on a phylogeny, which implies that deviations from the Poisson expectation arise from gene-specific temporal variation in substitution rates. Amino acid sequences show greater temporal variation in substitution rates than do four-fold synonymous sequences. Our findings provide a general phenomenological framework for understanding overdispersion in the molecular clock. Also, the presence of substantial variation in gene-specific substitution rates has broad implications for work in phylogeny reconstruction and evolutionary rate estimation.     

Neighboring-nucleotide effects on the rates of germ-line single-base-pair substitution in human genes.

The spectrum of single-base-pair substitutions logged in The Human Gene Mutation Database (HGMD), comprising 7,271 different lesions in the coding regions of 547 different human genes, was analyzed for nearest-neighbor effects on relative mutation rates. Owing to its retrospective nature, HGMD allows mutation rates to be estimated only in relative terms. Therefore, a novel methodology was devised in order to obtain these estimates in iterative fashion, correcting, at the same time, for the confounding effects of differential codon usage and for the fact that different types of amino acid replacement come to clinical attention with different probabilities. Over and above the hypermutability of CpG dinucleotides, reflected in transition rates five times the base mutation rate, only a subtle and locally confined influence of the surrounding DNA sequence on relative single-base-pair substitution rates was observed, which extended no farther than 2 bp from the substitution site. A disparity between the two DNA strands was evidenced by the fact that, when substitution rates were estimated conditional on the 5' and 3' flanking nucleotides, a significant rate difference emerged for 10 of 96 possible pairs of complementary substitutional events. Mutational bias, favoring substitutions toward flanking bases, a phenomenon reminiscent of misalignment mutagenesis, was apparent and exhibited both directionality and reading-frame sensitivity. No specific preponderance of repeat-sequence motifs was observed in the vicinity of nucleotide substitutions, but a moderate correlation between the relative mutability and thermodynamic stability of DNA triplets emerged, suggesting either inefficient DNA replication in regions of high stability or the transient stabilization of misaligned intermediates.     

Partitioning the variation in mammalian substitution rates

We have used analysis of variance to partition the variation in synonymous and amino acid substitution rates between three effects (gene, lineage, and a gene-by-lineage interaction) in mammalian nuclear and mitochondrial genes. We find that gene effects are stronger for amino acid substitution rates than for synonymous substitution rates and that lineage effects are stronger for synonymous substitution rates than for amino acid substitution rates. Gene-by-lineage interactions, equivalent to overdispersion corrected for lineage effects, are found in amino acid substitutions but not in synonymous substitutions. The variance in the ratio of amino acid and synonymous substitution rates is dominated by gene effects, but there is also a significant gene-by-lineage interaction.     

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.     

Correlation of Functional Domains and Rates of Nucleotide Substitution in Cytochromeb

Distinguishing noise from signal presents a problem when DNA sequences are used for phylogeny reconstruction. Multiple substitutions at sites are a primary cause of noise and this is compounded by variation in substitution rates among sites. For protein-coding genes, one method used to determine if data are noisy is to assess levels of saturation of substitutions by codon position. However, this procedure may not be a fine enough filter for assessing noise. Variation in substitution rates may also be caused by constraints on change imposed by the function of the protein product. Using a structural model of the cytochromebprotein as a template, I divided cytbsequence data for species within the avian family Falconidae (falcons and caracaras) into three functional domains. Saturation of substitutions of sequences within these regions was assessed graphically. This qualitative determination of saturation was then used to differentially weight phylogenetic analysis, resulting in an hypothesis congruent with existing cladistic analyses and traditional morphology. These results demonstrate that saturation of substitutions is correlated with functional regions of cytochromeband that using this information improves phylogenetic inference.     

Rapidly evolving mitochondrial genome and directional selection in mitochondrial genes in the parasitic wasp nasonia (hymenoptera: pteromalidae)

We sequenced the nearly complete mtDNA of 3 species of parasitic wasps, Nasonia vitripennis (2 strains), Nasonia giraulti, and Nasonia longicornis, including all 13 protein-coding genes and the 2 rRNAs, and found unusual patterns of mitochondrial evolution. The Nasonia mtDNA has a unique gene order compared with other insect mtDNAs due to multiple rearrangements. The mtDNAs of these wasps also show nucleotide substitution rates over 30 times faster than nuclear protein-coding genes, indicating among the highest substitution rates found in animal mitochondria (normally <10 times faster). A McDonald and Kreitman test shows that the between-species frequency of fixed replacement sites relative to silent sites is significantly higher compared with within-species polymorphisms in 2 mitochondrial genes of Nasonia, atp6 and atp8, indicating directional selection. Consistent with this interpretation, the Ka/Ks (nonsynonymous/synonymous substitution rates) ratios are higher between species than within species. In contrast, cox1 shows a signature of purifying selection for amino acid sequence conservation, although rates of amino acid substitutions are still higher than for comparable insects. The mitochondrial-encoded polypeptides atp6 and atp8 both occur in F0F1ATP synthase of the electron transport chain. Because malfunction in this fundamental protein severely affects fitness, we suggest that the accelerated accumulation of replacements is due to beneficial mutations necessary to compensate mild-deleterious mutations fixed by random genetic drift or Wolbachia sweeps in the fast evolving mitochondria of Nasonia. We further propose that relatively high rates of amino acid substitution in some mitochondrial genes can be driven by a "Compensation-Draft Feedback"; increased fixation of mildly deleterious mutations results in selection for compensatory mutations, which lead to fixation of additional deleterious mutations in nonrecombining mitochondrial genomes, thus accelerating the process of amino acid substitutions.     

Molecular evolutionary rates of oncogenes

Summary Using nine sets of viral and cellular oncogenes, the rates of nucleotide substitutions were computed by using Gojobori and Yokoyama's (1985) method. The results obtained confirmed our previous conclusion that the rates of nucleotide substitution for the viral oncogenes are about a million times higher than those for their cellular counterparts. For cellular oncogenes and most viral oncogenes, however, the rate of synonymous substitution is higher than that of nonsynonymous substitution. Moreover, the pattern of nucleotide substitutions for viral oncogenes is more similar to that for functional genes (such as cellular oncongenes) than for pseudogenes. This implies that nucleotide substitutions in viral oncogenes may be functionally constrained. Thus, our observation supports that nucleotide substitutions for the oncogenes in those DNA and RNA genomes are consistent with Kimura's neutral theory of molecular evolution (Kimura 1968, 1983).     

Molecular Evolution of the Plant R Regulatory Gene Family

Anthocyanin pigmentation patterns in different plant species are controlled in part by members of the myc-like R regulatory gene family. We have examined the molecular evolution of this gene family in seven plant species. Three regions of the R protein show sequence conservation between monocot and dicot R genes. These regions encode the basic helix-loop-helix domain, as well as conserved N-terminal and C-terminal domains; mean replacement rates for these conserved regions are 1.02 X 10(-9) nonsynonymous nucleotide substitutions per site per year. More than one-half of the protein, however, is diverging rapidly, with nonsynonymous substitution rates of 4.08 X 10(-9) substitutions per site per year. Detailed analysis of R homologs within the grasses (Poaceae) confirm that these variable regions are indeed evolving faster than the flanking conserved domains. Both nucleotide substitutions and small insertion/deletions contribute to the diversification of the variable regions within these regulatory genes. These results demonstrate that large tracts of sequence in these regulatory loci are evolving at a fairly rapid rate.     

Phylogeny, rate variation, and genome size evolution of Pelargonium (Geraniaceae)

The phylogeny of 58 Pelargonium species was estimated using five plastid markers (rbcL, matK, ndhF, rpoC1, trnL-F) and one mitochondrial gene (nad5). The results confirmed the monophyly of three major clades and four subclades within Pelargonium but also indicate the need to revise some sectional classifications. This phylogeny was used to examine karyotype evolution in the genus: plotting chromosome sizes, numbers and 2C-values indicates that genome size is significantly correlated with chromosome size but not number. Accelerated rates of nucleotide substitution have been previously detected in both plastid and mitochondrial genes in Pelargonium, but sparse taxon sampling did not enable identification of the phylogenetic distribution of these elevated rates. Using the multigene phylogeny as a constraint, we investigated lineage- and locus-specific heterogeneity of substitution rates in Pelargonium for an expanded number of taxa and demonstrated that both plastid and mitochondrial genes have had accelerated substitution rates but with markedly disparate patterns. In the plastid, the exons of rpoC1 have significantly accelerated substitution rates compared to its intron and the acceleration was mainly due to nonsynonymous substitutions. In contrast, the mitochondrial gene, nad5, experienced substantial acceleration of synonymous substitution rates in three internal branches of Pelargonium, but this acceleration ceased in all terminal branches. Several lineages also have dN/dS ratios significantly greater than one for rpoC1, indicating that positive selection is acting on this gene, whereas the accelerated synonymous substitutions in the mitochondrial gene are the result of elevated mutation rates.     

Rapid evolution of the sex-determining gene,transformer: structural diversity and rate heterogeneity among sibling species of Drosophila

While developmentally regulated genes are generally conserved, transformer (tra), a key locus involved in the regulation of sexual differentiation, is highly diverged between species of Drosophila. With an aim to understand its divergence between sibling species, we investigated tra sequence variation among members of the Drosophila melanogaster species complex, D. melanogaster, D. simulans, D. mauritiana, and D. sechellia. In this species group, tra divergence is rapid yet clocklike and exhibits large differences in protein size. D. melanogaster contains a 13-amino acid tandem duplication, whereas D. sechellia possesses a 72-amino acid tandem duplication representing a 30% increase in total amino acid residues. We also found evidence of a nonrandom distribution of replacement substitutions and heterogeneity in substitution rates using clustering statistics and a codon substitution model. We show that tra's rapid divergence in this species complex is the result of generally lower selective constraints around regions that encode arginine-serine (RS) domains and a significantly higher rate of substitutions around the insertion site of D. sechellia's large duplication. The proximity of rapidly diverged regions to sites of nucleotide insertion suggests that higher local rates of mutation may provide a causal mechanism for TRA's rapid divergence in this subgroup. A comparison of tra orthologs across the genus Drosophila suggest that TRA maintains an assortment of RS domains for proper sex determining function while much of the protein evolves relatively unconstrained.     

The effect of tandem substitutions on the correlation between synonymous and nonsynonymous rates in rodents.

Nonsynonymous substitutions in DNA cause amino acid substitutions while synonymous substitutions in DNA leave amino acids unchanged. The cause of the correlation between the substitution rates at nonsynonymous (K(A)) and synonymous (K(S)) sites in mammals is a contentious issue, and one that impacts on many aspects of molecular evolution. Here we use a large set of orthologous mammalian genes to investigate the causes of the K(A)-K(S) correlation in rodents. The strength of the K(A)-K(S) correlation exceeds the neutral theory expectation when substitution rates are estimated using algorithmic methods, but not when substitution rates are estimated by maximum likelihood. Irrespective of this methodological uncertainty the strength of the K(A)-K(S) correlation appears mostly due to tandem substitutions, an excess of which is generated by substitutional nonindependence. Doublet mutations cannot explain the excess of tandem synonymous-nonsynonymous substitutions, and substitution patterns indicate that selection on silent sites is the likely cause. We find no evidence for selection on codon usage. The nature of the relationship between synonymous divergence and base composition is unclear because we find a significant correlation if we use maximum-likelihood methods but not if we use algorithmic methods. Finally, we find that K(S) is reduced at the start of genes, which suggests that selection for RNA structure may affect silent sites in mammalian protein-coding genes.     

Models of amino acid substitution and applications to mitochondrial protein evolution

Models of amino acid substitution were developed and compared using maximum likelihood. Two kinds of models are considered. "Empirical" models do not explicitly consider factors that shape protein evolution, but attempt to summarize the substitution pattern from large quantities of real data. "Mechanistic" models are formulated at the codon level and separate mutational biases at the nucleotide level from selective constraints at the amino acid level. They account for features of sequence evolution, such as transition-transversion bias and base or codon frequency biases, and make use of physicochemical distances between amino acids to specify nonsynonymous substitution rates. A general approach is presented that transforms a Markov model of codon substitution into a model of amino acid replacement. Protein sequences from the entire mitochondrial genomes of 20 mammalian species were analyzed using different models. The mechanistic models were found to fit the data better than empirical models derived from large databases. Both the mutational distance between amino acids (determined by the genetic code and mutational biases such as the transition-transversion bias) and the physicochemical distance are found to have strong effects on amino acid substitution rates. A significant proportion of amino acid substitutions appeared to have involved more than one codon position, indicating that nucleotide substitutions at neighboring sites may be correlated. Rates of amino acid substitution were found to be highly variable among sites.      

Molecular considerations in the evolution of bacterial genes

Summary Synonymous and nonsynonymous substitution rates at the loci encoding glyceraldehyde-3-phosphate dehydrogenase (gap) and outer membrane protein 3A (ompA) were examined in 12 species of enteric bacteria. By examining homologous sequences in species of varying degrees of relatedness and of known phylogenetic relationships, we analyzed the patterns of synonymous and nonsynonymous substitutions within and among these genes. Although both loci accumulate synonymous substitutions at reduced rates due to codon usage bias, portions of thegap andompA reading frames show significant deviation in synonymous substitution rates not attributable to local codon bias. A paucity of synonymous substitutions in portions of theompA gene may reflect selection for a novel mRNA secondary structure. In addition, these studies allow comparisons of homologous protein-coding sequences (gap) in plants, animals, and bacteria, revealing differences in evolutionary constraints on this glycolytic enzyme in these lineages.     

Codon Substitution in Evolution and the "Saturation" of Synonymous Changes

A mathematical model for codon substitution is presented, taking into account unequal mutation rates among different nucleotides and purifying selection. This model is constructed by using a 61 X 61 transition probability matrix for the 61 nonterminating codons. Under this model, a computer simulation is conducted to study the numbers of silent (synonymous) and amino acid-altering (nonsynonymous) nucleotide substitutions when the underlying mutation rates among the four kinds of nucleotides are not equal. It is assumed that the substitution rates are constant over evolutionary time, the codon frequencies being in equilibrium, and, thus, the numbers of synonymous and nonsynonymous substitutions both increase linearly with evolutionary time. It is shown that, when the mutation rates are not equal, the estimate of synonymous substitutions obtained by F. Perler, A. Efstratiadis, P. Lomedico, W. Gilbert, R. Kolodner and J. Dodgson's "Percent Corrected Divergence" method increases nonlinearly, although the true number of synonymous substitutions increases linearly. It is, therefore, possible that the "saturation" of synonymous substitutions observed by Perler et al. is due to the inefficiency of their method to detect all synonymous substitutions.     

Relative rates of synonymous substitutions in the mitochondrial, chloroplast and nuclear genomes of seed plants

Previous studies have estimated that, in angiosperms, the synonymous substitution rate of chloroplast genes is three times higher than that of mitochondrial genes and that of nuclear genes is twelve times higher than that of mitochondrial genes. Here we used 12 genes in 27 seed plant species to investigate whether these relative rates of substitutions are common to diverse seed plant groups. We find that the overall relative rate of synonymous substitutions of mitochondrial, chloroplast and nuclear genes of all seed plants is 1:3:10, that these ratios are 1:2:4 in gymnosperms but 1:3:16 in angiosperms and that they go up to 1:3:20 in basal angiosperms. Our results show that the mitochondrial, chloroplast and nuclear genomes of seed plant groups have different synonymous substitutions rates, that these rates are different in different seed plant groups and that gymnosperms have smaller ratios than angiosperms.     

Empirical relationship between the number of nucleotide substitutions and interspecific identity of amino acid sequences in some proteins

There are three different methods of estimating the number of nucleotide substitutions between a pair of species from amino acid sequence data, i.e. the Poisson correction method, random evolutionary hit method, and counting the actual but minimum number of nucleotide substitutions. In this paper the relationships among the estimates obtained by these methods are studied empirically. The results obtained indicate that there is a high correlation among these estimates and in practice any of the three methods may be used for constructing evolutionary trees or relating nucleotide substitutions to evolutionary time. The effects of varying rates of nucleotide substition among different sites on the Poisson correction and random evolutionary hit methods are also studied mathematically. It is shown that these two methods are quite insensitive to the variation of the rate of nucleotide substitution.     

Uniformity in the nonsynonymous substitution rates of embryonic β-globin genes of several vertebrate species

Summary The nucleotide substitution rate in structural portions of the embryonic β-globin genes of placental mammals is lower than that for the adult β-globin genes. This difference occurs entirely within the class of substitutions that result in nonsynonymous (replacement) differences between these genes, and therefore represents a constraint on the structure of the mammalian embryonic β-globin proteins relative to the adult proteins (Shapiro et al. 1983; Hardison 1984). A similar effect has also been observed in marsupial mammals (Koop and Goodman 1988). In an effort to determine whether the observed rates are evidence of a uniform degree of selective constraint on the embryonic β-globin genes, analyses were performed that compared replacement substitution rates. The analyses reveal that embryonic β-globin genes appear to have been fixing replacement substitutions at nearly the same average rate not only in placental and marsupial mammals but in avian and amphibian species as well. In contrast, the adult β-globin genes from these organisms appear to have a more variable rate of replacement substitution with an especially low rate for birds. In the chicken (Gallus gallus), the adult β-globin gene replacement substitution rate appears to be lower than the embryonic replacement substitution rate.     

A large variation in the rates of synonymous substitution for RNA viruses and its relationship to a diversity of viral infection and transmission modes

RNA viruses successfully adapt to various environments by repeatedly producing new mutants, often through generating a number of nucleotide substitutions. To estimate the degree of variation in mutation rates of RNA viruses and to understand the source of such variation, we studied the synonymous substitution rate because synonymous substitution is exempt from functional constraints at the protein level, and its rate reflects the mutation rate to a great extent. We estimated the synonymous substitution rates for a total of 49 different species of RNA viruses, and we found that the rates had tremendous variation by 5 orders of magnitude (from 1.3 x 10(-7) to 6.2 x 10(-2) /synonymous site/year). Comparing the synonymous substitution rates with the replication frequencies and replication error rates for the RNA viruses, we found that the main source of the rate variation was differences in the replication frequency because the rates of replication error were roughly constant over different RNA viruses. Moreover, we examined a relationship between viral life strategies and synonymous substitution rates to understand which viral life strategies affect replication frequencies. The results show that the variation of synonymous substitution rates has been influenced most by either the difference in the infection modes or the differences in the transmission modes. In conclusion, the variation of mutation rates for RNA viruses is caused by different replication frequencies, which are affected strongly by the infection and transmission modes.     

Nucleic acid composition,codon usage,and the rate of synonymous substitution in protein-coding genes

Summary Based on the rates of synonymous substitution in 42 protein-codin gene pairs from rat and human, a correlation is shown to exist between the frequency of the nucleotides in all positions of the codon and the synonymous substitution rate. The correlation coefficients were positive for A and T and negative for C and G. This means that AT-rich genes accumulate more synonymous substitutions than GC-rich genes. Biased patterns of mutation could not account for this phenomenon. Thus, the variation in synonymous substitution rates and the resulting unequal codon usage must be the consequence of selection against A and T in synonymous positions. Most of the varition in rates of synonymous substitution can be explained by the nucleotide composition in synonymous positions. Codon-anticodon interactions, dinucleotide frequencies, and contextual factors influence neither the rates of synonymous substitution nor codon usage. Interestingly, the nucleotide in the second position of codons (always a nonsynonymous position) was found to affect the rate of synonymous substitution. This finding links the rate of nonsynonymous substitution with the synonymous rate. Consequently, highly conservative proteins are expected to be encoded by genes that evolve slowly in terms of synonymous substitutions, and are consequently highly biased in their codon usage.