Maximum-likelihood analysis of molecular adaptation in abalone sperm lysin reveals variable selective pressures among lineages and sites

Maximum-likelihood models of codon substitution were used to analyze sperm lysin genes of 25 abalone (HALIOTIS:) species to identify lineages and amino acid sites under diversifying selection. The models used the nonsynonymous/synonymous rate ratio (omega = d(N)/d(S)) as an indicator of selective pressure and allowed the ratio to vary among lineages or sites. Likelihood ratio tests suggested significant variation in selective pressure among lineages. The variable selective pressure provided an explanation for the previous observation that the omega ratio is >1 in comparisons of closely related species and <1 in comparisons of distantly related species. Computer simulations demonstrated that saturation of nonsynonymous substitutions and constraint on lysin structure were unlikely to account for the observed pattern. Lineages linking closely related sympatric species appeared to be under diversifying selection, while lineages separating distantly related species from different geographic locations were associated with low evolutionary rates. The selective pressure indicated by the omega ratio was found to vary greatly among amino acid sites in lysin. Sites under potential diversifying selection were identified. Ancestral lysins were inferred to trace the route of evolution at individual sites and to provide lysin sequences for future laboratory studies.     

Looking for darwin in genomic sequences--validity and success of statistical methods

The use of codon substitution models to compare synonymous and nonsynonymous substitution rates is a widely used approach to detecting positive Darwinian selection affecting protein evolution. However, in several recent papers, Hughes and colleagues claim that codon-based likelihood-ratio tests (LRTs) are logically flawed as they lack prior hypotheses and fail to accommodate random fluctuations in synonymous and nonsynonymous substitutions Friedman and Hughes (2007) also used site-based LRTs to analyze 605 gene families consisting of human and mouse paralogues. They found that the outcome of the tests was largely determined by irrelevant factors such as the GC content at the third codon positions and the synonymous rate d(S), but not by the nonsynonymous rate d(N) or the d(N)/d(S) ratio, factors that should be related to selection. Here, we reanalyze those data. Contra Friedman and Hughes, we found that the test results are related to sequence length and the average d(N)/d(S) ratio. We examine the criticisms of Hughes and suggest that they are based on misunderstandings of the codon models and on statistical errors. Our analyses suggest that codon-based tests are useful tools for comparative analysis of genomic data sets.     

Positive and negative selection in the DAZ gene family

Because a microdeletion containing the DAZ gene is the most frequently observed deletion in infertile men, the DAZ gene was considered a strong candidate for the azoospermia factor. A recent evolutionary analysis, however, suggested that DAZ was free from functional constraints and consequently played little or no role in human spermatogenesis. The major evidence for this surprising conclusion is that the nonsynonymous substitution rate is similar to the synonymous rate and to the rate in introns. In this study, we reexamined the evolution of the DAZ gene family by using maximum-likelihood methods, which accommodate variable selective pressures among sites or among branches. The results suggest that DAZ is not free from functional constraints. Most amino acids in DAZ are under strong selective constraint, while a few sites are under diversifying selection with nonsynonymous/ synonymous rate ratios (d(N)/d(S)) well above 1. As a result, the average d(N)/d(S) ratio over sites is not a sensible measure of selective pressure on the protein. Lineage-specific analysis indicated that human members of this gene family were evolving by positive Darwinian selection, although the evidence was not strong.     

Adaptive evolution of the IgA hinge region in primates

IgA is a major component that prevents the penetration of pathogenic bacteria into mucosal surfaces. The IgA antibody is cleaved at the IgA hinge region with high specificity by IgA-specific proteases produced by several pathogenic bacteria. We conducted a genomic sequence analysis of the IgA genes of a wide spectrum of primates, including the first intron and second exon, which consist of the hinge region and the CH2 domain, to find evidence of positive selection. Because the hinge region is quite small, we combined the largest collection of sequences that could be clearly aligned and evaluated the total numbers of synonymous and nonsynonymous substitutions on the phylogenetic tree. The nonsynonymous to synonymous substitution ratio (d(N)/d(S) test) showed that hominoids, Old World monkeys, and New World monkeys have d(N)/d(S) ratios of 5.4, 6.3, and 4.2, respectively. Fisher's exact probability tests showed statistical significance for the Old World monkey. Because the substitution rates of the flanking sequences are more or less similar to the synonymous rates of the hinge region, these high values of d(N)/d(S) should be the result of positive selection at the hinge region. Combining the high sequence variability in each population and the highly accelerated nonsynonymous substitution rates in the hinge region, we conclude that this unusual IgA evolution is a molecular evidence of adaptive evolution possibly caused by the host-parasite relationship.     

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.     

Molecular evolution of the COX7A gene family in primates.

COX VIIa is one of 10 nuclear-encoded subunits of the COX holoenzyme, and one of three that have isoforms with tissue-specific differences in expression. Analysis of nucleotide substitution rates revealed an accelerated rate of nonsynonymous substitutions relative to that of synonymous substitutions for the heart isoform gene (COX7AH) in six primate lineages. Rate accelerations have been noted for four other COX-related genes in this time period, suggesting that the COX holoenzyme has experienced an episode of adaptive evolution. A third member of the gene family, COX7AR, has recently been described. Although its function is currently unknown, low nonsynonymous substitution/synonymous substitution (N/S) ratios in mammalian evolution suggest that COX7AR is of functional importance. When the COX7A isoforms were divided into domains, examination of nucleotide substitution rates suggested that mitochondrial targeting residues experienced an accelerated nonsynonymous substitution rate in the period following gene duplication. In contrast, paralogous comparisons of the targeting residues of each isoform show they have been relatively conserved in mammalian evolution. This pattern is consistent with the evolution of tissue-specific function.     

Molecular evolution of FLORICAULA/LEAFY orthologs in the Andropogoneae (Poaceae)

Members of the grass family (Poaceae) exhibit a broad range of inflorescence structures and other morphologies, making the grasses an interesting model system for studying the evolution of development. Here we present an analysis of the molecular evolution of FLORICAULA/LEAFY-like genes, which are important developmental regulatory loci known to affect inflorescence development in a wide range of flowering plant species. We have focused on sequences from the Andropogoneae, a tribe within the grass family that includes maize (Zea mays ssp. mays) and Sorghum (Sorghum bicolor). The FLORICAULA/LEAFY gene phylogeny we generated largely agrees with previously published phylogenies for the Andropogoneae using other nuclear genes but is unique in that it includes both members of one of the many duplicate gene sets present in maize. The placement of these sequences in the phylogeny suggests that the duplication of the maize FLORICAULA/LEAFY orthologs, zfl1 and zfl2, is a consequence of a proposed tetraploidy event that occurred in the common ancestor of Zea and a closely related genus, Tripsacum. Our data are consistent with the hypothesis that the transcribed regions of the FLORICAULA/LEAFY-like genes in the Andropogoneae are functionally constrained at both nonsynonymous and synonymous sites and show no evidence of directional selection. We also examined conservation of short noncoding sequences in the first intron, which may play a role in gene regulation. Finally, we investigated the genetic diversity of one of the two maize FLORICAULA/LEAFY orthologs, zfl2, in maize and its wild ancestor, teosinte (Z. mays ssp. parviglumis), and found no evidence for selection pressure resulting from maize domestication within the zfl2-coding region.     

Molecular evolution of nuclear genes in Cupressacea, a group of conifer trees

We surveyed the molecular evolutionary characteristics of 11 nuclear genes from 10 conifer trees belonging to the Taxodioideae, the Cupressoideae, and the Sequoioideae. Comparisons of substitution rates among the lineages indicated that the synonymous substitution rates of the Cupressoideae lineage were higher than those of the Taxodioideae. This result parallels the pattern previously found in plastid genes. Likelihood-ratio tests showed that the nonsynonymous-synonymous rate ratio did not change significantly among lineages. In addition, after adjustments for lineage effects, the dispersion indices of synonymous and nonsynonymous substitutions were considerably reduced, and the latter was close to 1. These results indicated that the acceleration of evolutionary rates in the Cupressoideae lineage occurred in both the nuclear and plastid genomes, and that generally, this lineage effect affected synonymous and nonsynonymous substitutions similarly. We also investigated the relationship of synonymous substitution rates with the nonsynonymous substitution rate, base composition, and codon bias in each lineage. Synonymous substitution rates were positively correlated with nonsynonymous substitution rates and GC content at third codon positions, but synonymous substitution rates were not correlated with codon bias. Finally, we tested the possibility of positive selection at the protein level, using maximum likelihood models, assuming heterogeneous nonsynonymous-synonymous rate ratios among codon (amino acid) sites. Although we did not detect strong evidence of positively selected codon sites, the analysis suggested that significant variation in nonsynonymous-synonymous rate ratio exists among the sites. The most likely sites for action of positive selection were found in the ferredoxin gene, which is an important component of the apparatus for photosynthesis.     

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.     

Codon-substitution models for detecting molecular adaptation at individual sites along specific lineages

The nonsynonymous (amino acid-altering) to synonymous (silent) substitution rate ratio (omega = d(N)/d(S)) provides a measure of natural selection at the protein level, with omega = 1, >1, and <1, indicating neutral evolution, purifying selection, and positive selection, respectively. Previous studies that used this measure to detect positive selection have often taken an approach of pairwise comparison, estimating substitution rates by averaging over all sites in the protein. As most amino acids in a functional protein are under structural and functional constraints and adaptive evolution probably affects only a few sites at a few time points, this approach of averaging rates over sites and over time has little power. Previously, we developed codon-based substitution models that allow the omega ratio to vary either among lineages or among sites. In this paper we extend previous models to allow the omega ratio to vary both among sites and among lineages and implement the new models in the likelihood framework. These models may be useful for identifying positive selection along prespecified lineages that affects only a few sites in the protein. We apply those branch-site models as well as previous branch- and site-specific models to three data sets: the lysozyme genes from primates, the tumor suppressor BRCA1 genes from primates, and the phytochrome (PHY) gene family in angiosperms. Positive selection is detected in the lysozyme and BRCA genes by both the new and the old models. However, only the new models detected positive selection acting on lineages after gene duplication in the PHY gene family. Additional tests on several data sets suggest that the new models may be useful in detecting positive selection after gene duplication in gene family evolution.     

Likelihood ratio tests for detecting positive selection and application to primate lysozyme evolution

An excess of nonsynonymous substitutions over synonymous ones is an important indicator of positive selection at the molecular level. A lineage that underwent Darwinian selection may have a nonsynonymous/synonymous rate ratio (dN/dS) that is different from those of other lineages or greater than one. In this paper, several codon-based likelihood models that allow for variable dN/dS ratios among lineages were developed. They were then used to construct likelihood ratio tests to examine whether the dN/dS ratio is variable among evolutionary lineages, whether the ratio for a few lineages of interest is different from the background ratio for other lineages in the phylogeny, and whether the dN/dS ratio for the lineages of interest is greater than one. The tests were applied to the lysozyme genes of 24 primate species. The dN/dS ratios were found to differ significantly among lineages, indicating that the evolution of primate lysozymes is episodic, which is incompatible with the neutral theory. Maximum- likelihood estimates of parameters suggested that about nine nonsynonymous and zero synonymous nucleotide substitutions occurred in the lineage leading to hominoids, and the dN/dS ratio for that lineage is significantly greater than one. The corresponding estimates for the lineage ancestral to colobine monkeys were nine and one, and the dN/dS ratio for the lineage is not significantly greater than one, although it is significantly higher than the background ratio. The likelihood analysis thus confirmed most, but not all, conclusions Messier and Stewart reached using reconstructed ancestral sequences to estimate synonymous and nonsynonymous rates for different lineages.      

Estimating a nucleotide substitution rate for maize from polymorphism at a major domestication locus

To estimate a rate for single nucleotide substitutions for maize (Zea mays ssp. mays), we have taken advantage of data from genetic and archaeological studies of the domestication of maize from its wild ancestor, teosinte (Z. mays ssp. parviglumis). Genetic studies have shown that the teosinte branched1 (tb1) gene was a major target of human selection during maize domestication, and sequence diversity in the intergenic region 5' to the tb1-coding sequence is extraordinarily low. We show that polymorphism in this region is consistent with new mutation following fixation for a small number of tb1 haplotypes during domestication. Archeological studies suggest that maize was domesticated approximately 6,250-10,000 years ago and subsequently the size of the maize population is thought to have expanded rapidly. Using the observed number of mutations within the region of selection at tb1, the approximate age of maize domestication, and approximations for the maize genealogy, we have derived estimates for the nucleotide substitution rate for the tb1 intergenic region. Using two approaches, one of which is a coalescent approach, we obtain rate estimates of approximately 2.9 x 10(-8) and 3.3 x 10(-8) substitutions per site per year. We also show that the pattern of polymorphism in the tb1 intergenic region appears to have been strongly affected by the mutagenic effect of DNA methylation. Excluding target sites of symmetric DNA methylation (CG and CNG sites) from analysis, the mutation rate estimates are reduced by approximately 50%-60%, while the rates for CG and CNG sites are nearly an order of magnitude higher. We use rate estimates from the tb1 region to estimate the timing of expansion of transposable elements in the maize genome and suggest that this expansion occurred primarily within the last million years.     

Molecular phylogenies and evolution of crt genes in algae

The carotenoids constitute the most widespread class of pigments in nature. Most previous work has concentrated on the identification and characterization of their chemical physical properties and bioavailability. In recent years, significant amounts of research have been conducted in an attempt to analyze the genes and the molecular regulation of the genes involved in the biosynthesis of carotenoids. However, it is important not to lose sight of the early evolution of carotenoid biosynthesis. One of the major obstacles in understanding the evolution of the respective enzymes and their patterns of selection is a lack of a well-supported phylogenic analysis. In the present research, a major long-term objective was to provide a clearer picture of the evolutionary history of genes, together with an evaluation of the patterns of selection in algae. These phylogenies will be important in studies characterizing the evolution of algae. The gene sequences of the enzymes involved in the major steps of the carotenoid biosynthetic pathway in algae (cyanobacteria, rhofophyta, chlorophyta) have been analyzed. Phylogenetic relationships among protein-coding DNA sequences were reconstructed by neighbor-joining (NJ) analysis for the respective carotenoid biosynthetic pathway genes (crt) in algae. The analysis also contains an estimation of the rate of nonsynonymous nucleotide substitutions per nonsynonymous site (d(N)), synonymous nucleotide substitution per synonymous site (d(S)), and the ratio of nonsynonmous (d(N)/d(S)) for the test of selection patterns. The phylogenetic trees show that the taxa of some genera have a closer evolutionary relationship with other genera in some gene sequences, which suggests a common ancient origin and that lateral gene transfer has occurred among unrelated genera. The d(N) values of crt genes in the early pathway are relatively low, while those of the following steps are slightly higher, while the d(N) values of crt genes in chlorophyta are higher than those in cyanobacteria. Most of the d(N)/d(S) values exceed 1. The phylogenetic analysis revealed that lateral gene transfer may have taken place across algal genomes and the d(N) values suggest that most of the early crt genes are well conserved compared to the later crt genes. Furthermore, d(N) values also revealed that the crt genes of chlorophyta are more evolutionary than cyanobacteria. The amino acids' changes are mostly adaptive evolution under the influence of positive diversity selection.     

Overestimation of nonsynonymous/synonymous rate ratio by reverse-translation of aligned amino acid sequences

In the analysis of protein-coding nucleotide sequences, the ratio of the number of nonsynonymous substitutions to that of synonymous substitutions (d(N)/d(S)) is used as an indicator for the direction and magnitude of natural selection operating at the amino acid sequence level. The d(S) and d(N) values are estimated based on the comparison of homologous codons, which are often identified by converting (reverse-translating) aligned amino acid sequences into codon sequences. In this method, however, homologous codons may be mis-identified when frame-shifts occurred or amino acid sequences were mis-aligned, which may lead to overestimation of the d(N)/d(S) ratio. Here the effect of reverse-translating aligned amino acid sequences on the estimation of d(N)/d(S) ratio was examined through a large-scale analysis of protein-coding nucleotide sequences from vertebrate species. Apparently, 1-9% of codon sites that were identified as homologous with reverse-translation contained non-homologous codons, where the d(N)/d(S) ratio was unduly high. By correcting the d(N)/d(S) ratio for these codon sites, it was inferred that the ratio was 5-43% overestimated with reverse-translation. These results suggest that caution should be exerted in the study of natural selection using the d(N)/d(S) ratio by reverse-translating aligned amino acid sequences.     

Sexual selection and the adaptive evolution of mammalian ejaculate proteins

An elevated rate of substitution characterizes the molecular evolution of reproductive proteins from a wide range of taxa. Although the selective pressures explaining this rapid evolution are yet to be resolved, recent evidence implicates sexual selection as a potentially important explanatory factor. To investigate this hypothesis, we sought evidence of a high rate of adaptive gene evolution linked to postcopulatory sexual selection in muroid rodents, a model vertebrate group displaying a broad range of mating systems. Specifically, we sequenced 7 genes from diverse rodents that are expressed in the testes, prostate, or seminal vesicles, products of which have the potential to act in sperm competition. We inferred positive Darwinian selection in these genes by estimation of the ratio of nonsynonymous (d(N), amino acid changing) to synonymous (d(S), amino acid retaining) substitution rates (omega = d(N)/d(S)). Next, we tested whether variation in this ratio among lineages could be attributed to interspecific variation in mating systems, as inferred from the variation in these rodents' relative testis sizes (RTS). Four of the 7 genes examined (Prm1, Sva, Acrv1, and Svs2, but not Svp2, Msmb, or Spink3) exhibit unambiguous evidence of positive selection. One of these, the seminal vesicle-derived protein Svs2, also shows some evidence for a concentration of positive selection in those lineages in which sperm competition is common. However, this was not a general trend among all the rodent genes we examined. Using the same methods, we then reanalyzed previously published data on 2 primate genes, SEMG1 and SEMG2. Although SEMG2 also shows evidence of positive selection concentrated in lineages subject to high levels of sperm competition, no such trend was found for SEMG1. Overall, despite a high rate of positive selection being a feature of many ejaculate proteins, these results indicate that the action of sexual selection potentially responsible for elevated evolutionary rates may be difficult to detect on a gene-by-gene basis. Although the extreme diversity of reproductive phenotypes exhibited in nature attests to the power of sexual selection, the extent to which this force predominates in driving the rapid molecular evolution of reproductive genes therefore remains to be determined.     

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 population genetics of the gene encoding the human fertilization protein zonadhesin reveals rapid adaptive evolution

A hallmark of positive selection (adaptive evolution) in protein-coding regions is a d(N)/d(S) ratio >1, where d(N) is the number of nonsynonymous substitutions/nonsynonymous sites and d(S) is the number of synonymous substitutions/synonymous sites. Zonadhesin is a male reproductive protein localized on the sperm head, comprising many domains known to be involved in cell-cell interaction or cell adhesion. Previous studies have shown that VWD domains (homologous to the D domains of the von Willebrand factor) are involved directly in binding to the female zona pellucida (ZP) in a species-specific manner. In this study, we sequenced 47 coding exons in 12 primate species and, by using maximum-likelihood methods to determine sites under positive selection, we show that VWD2, membrane/A5 antigen mu receptor, and mucin-like domains in zonadhesin are rapidly evolving and, thus, may be involved in binding to the ZP in a species-specific manner in primates. In addition, polymorphism data from 48 human individuals revealed significant polymorphism-to-divergence heterogeneity and a significant departure from equilibrium-neutral expectations in the frequency spectrum, suggesting balancing selection and positive selection occurring in zonadhesin (ZAN) within human populations. Finally, we observe adaptive evolution in haplotypes segregating for a frameshift mutation that was previously thought to indicate that ZAN was a potential pseudogene.     

Decreased diversity but increased substitution rate in host mtDNA as a consequence of Wolbachia endosymbiont infection

A substantial fraction of insects and other terrestrial arthropods are infected with parasitic, maternally transmitted endosymbiotic bacteria that manipulate host reproduction. In addition to imposing direct selection on the host to resist these effects, endosymbionts may also have indirect effects on the evolution of the mtDNA with which they are cotransmitted. Patterns of mtDNA diversity and evolution were examined in Drosophila recens, which is infected with the endosymbiont Wolbachia, and its uninfected sister species D. subquinaria. The level of mitochondrial, but not nuclear, DNA diversity is much lower in D. recens than in D. subquinaria, consistent with the hypothesized diversity-purging effects of an evolutionarily recent Wolbachia sweep. The d(N)/d(S) ratio in mtDNA is significantly greater in D. recens, suggesting that Muller's ratchet has brought about an increased rate of substitution of slightly deleterious mutations. The data also reveal elevated rates of synonymous substitutions in D. recens, suggesting that these sites may experience weak selection. These findings show that maternally transmitted endosymbionts can severely depress levels of mtDNA diversity within an infected host species, while accelerating the rate of divergence among mtDNA lineages in different species.