Positive selection of synonymous mutations in vesicular stomatitis virus

Prevailing evolutionary forces are typically deduced from the pattern of differences in synonymous and non-synonymous mutations, under the assumption of neutrality in the absence of amino acid change. We determined the complete sequence of ten vesicular stomatitis virus populations evolving under positive selection. A significant number of the mutations occurred independently in two or more strains, a process known as parallel evolution, and a substantial fraction of the parallel mutations were silent. Parallel evolution was also identified in non-coding regions. These results indicate that silent mutations can significantly contribute to adaptation in RNA viruses, and relative frequencies of synonymous and non-synonymous substitutions may not be useful to resolve their evolutionary history.     

Positive natural selection has driven the evolution of the Hsp70s in Diguetia spiders

Hsp70s are a ubiquitous family of highly conserved proteins. Hsp70s are chaperones and have important roles in both protein folding and thermotolerance. It has been widely assumed that Hsp70 sequence evolution is governed by the strong functional constraints imposed by its crucial cellular functions. In this study of cytosolic heat-inducible Hsp70s from three spider families, we have found clear evidence of positive natural selection altering Hsp70s in desert-dwelling and heat-loving Diguetidae spiders. These spiders are a small family restricted to deserts. They display heat-tolerant behaviours not seen in their closest relatives, the Pholcidae and Plectreuridae.     

Positive selection and subfunctionalization of duplicated CCT chaperonin subunits

To reach a functional and energetically stable conformation, many proteins need molecular helpers called chaperonins. Among the group II chaperonins, CCT proteins provide crucial machinery for the stabilization and proper folding of several proteins in the cytosol of eukaryotic cells through interactions that are subunit-specific and geometry-dependent. CCT proteins are made up of eight different subunits, all with similar sequences, positioned in a precise arrangement. Each subunit has been proposed to have a specialized function during the binding and folding of the CCT protein substrate. Here, we demonstrate that functional divergence occurred after several CCT duplication events due to the fixation of amino acid substitutions by positive selection. Sites critical for ATP binding and substrate binding were found to have undergone positive selection and functional divergence predominantly in subunits that bind tubulin but not actin. Furthermore, we show clear functional divergence between CCT subunits that bind the C-terminal domains of actin and tubulin and those that bind the N-terminal domains. Phylogenetic analyses could not resolve the deep relationships between most subunits, except for the groups alpha/beta/eta and delta/epsilon, suggesting several almost simultaneous ancient duplication events. Together, the results support the idea that, in contrast to homo-oligomeric chaperonins such as GroEL, the high divergence level between CCT subunits is the result of positive selection after each duplication event to provide a specialized role for each CCT subunit in the different steps of protein folding.     

Not so different after all: a comparison of methods for detecting amino acid sites under selection

We consider three approaches for estimating the rates of nonsynonymous and synonymous changes at each site in a sequence alignment in order to identify sites under positive or negative selection: (1) a suite of fast likelihood-based counting methods that employ either a single most likely ancestral reconstruction, weighting across all possible ancestral reconstructions, or sampling from ancestral reconstructions; (2) a random effects likelihood (REL) approach, which models variation in nonsynonymous and synonymous rates across sites according to a predefined distribution, with the selection pressure at an individual site inferred using an empirical Bayes approach; and (3) a fixed effects likelihood (FEL) method that directly estimates nonsynonymous and synonymous substitution rates at each site. All three methods incorporate flexible models of nucleotide substitution bias and variation in both nonsynonymous and synonymous substitution rates across sites, facilitating the comparison between the methods. We demonstrate that the results obtained using these approaches show broad agreement in levels of Type I and Type II error and in estimates of substitution rates. Counting methods are well suited for large alignments, for which there is high power to detect positive and negative selection, but appear to underestimate the substitution rate. A REL approach, which is more computationally intensive than counting methods, has higher power than counting methods to detect selection in data sets of intermediate size but may suffer from higher rates of false positives for small data sets. A FEL approach appears to capture the pattern of rate variation better than counting methods or random effects models, does not suffer from as many false positives as random effects models for data sets comprising few sequences, and can be efficiently parallelized. Our results suggest that previously reported differences between results obtained by counting methods and random effects models arise due to a combination of the conservative nature of counting-based methods, the failure of current random effects models to allow for variation in synonymous substitution rates, and the naive application of random effects models to extremely sparse data sets. We demonstrate our methods on sequence data from the human immunodeficiency virus type 1 env and pol genes and simulated alignments.     

Positive selectable marker genes for routine plant transformation

Summary Plant genetic transformation technologies rely upon the selection and recovery of transformed cells. Selectable marker genes used so far have been either antibiotic resistance genes or herbicide tolerance genes. There is a need to apply alternative principles of selection, as more transgenic traits have to be incorporated into a transgenic crop and because of concern that the use of conventional marker genes may pose a threat to humans and the environment. New classes of marker genes are now available, conferring metabolic advantage of the transgenic cells over the non-transformed cells. The new selection systems, as described in this review, are being used with success and superior performance over the traditional marker systems.     

Positive selection at reproductive ADAM genes with potential intercellular binding activity

Many genes with a role in reproduction, including those implicated in fertilization and spermatogenesis, have been shown to evolve at a faster rate relative to genes associated with other functions and tissues. These survey studies usually group a wide variety of genes with different characteristics and evolutionary histories as reproductive genes based on their site of expression or function. We have examined the molecular evolution of the ADAM (a disintegrin and metalloprotease) gene family, a structurally and functionally diverse group of genes expressed in reproductive and somatic tissue to test whether a variety of protein characteristics such as phylogenetic clusters, tissue of expression, and proteolytic and adhesive function can group fast evolving ADAM genes. We found that all genes were evolving under purifying selection (d(N)/d(S) < 1), although reproductive ADAMs, including those implicated in fertilization and spermatogenesis, evolved at the fastest rate. Genes with a role in binding to cell receptors in endogenous tissue appear to be evolving under purifying selection, regardless of the tissue of expression. In contrast, positive selection of codon sites in the disintegrin/cysteine-rich adhesion domains was detected exclusively in ADAMs 2 and 32, two genes expressed in the testis with a potential role in sperm-egg adhesion. Positive selection was detected in the transmembrane/cytosolic tail region of ADAM genes expressed in a variety of tissues.     

Positive selection of three chitinase genes of the family 18 of glycoside hydrolases in mammals

The digestive enzyme chitinase degrades chitin, and is found in a wide range of organisms, from prokaryotes to eukaryotes. Although mammals cannot synthesize or assimilate chitin, several proteins of the glycoside hydrolase (GH) chitinase family GH18, including some with enzymatic activity, have recently been identified from mammalian genomes. Consequently, there is growing interest in molecular evolution of this family of proteins. Here we report on the use of maximum likelihood methods to test for evidence of positive selection in three genes of the chitinase family GH18, all of which are found in mammals. These focal genes are CHIA, CHIT1 and CHI3L1, which encode the chitinase proteins acidic mammalian chitinase, chitotriosidase and cartilage protein 39, respectively. The results of our analyses indicate that each of these genes has undergone independent selective pressure in their evolution. Additionally, we have found evidence of a signature of positive natural selection, with most sites identified as being subject to adaptive evolution located in the catalytic domain. Our results suggest that positive selection on these genes stems from their function in digestion and/or immunity.     

Positive selection for indel substitutions in the rodent sperm protein catsper1

Catsper1 is a voltage-gated calcium channel located in the plasma membrane of the sperm tail and is necessary for sperm motility and fertility in mice. We here examine the evolutionary pattern of Catsper1 from nine species of the rodent subfamily Murinae of family Muridae. We show that the rate of insertion/deletion (indel) substitutions in exon 1 of the gene is 4-15 times that in introns or neutral genomic regions, suggesting the presence of strong positive selection that promotes fixations of indel mutations in exon 1. The number of indel polymorphisms within species appears higher than expected from interspecific comparisons, although there are too little data to provide a statistically significant conclusion. These results, together with an earlier report in primates, indicate that positive selection promoting length variation in Catsper1 may be widespread in mammals. A structural model of Catsper1 suggested the importance of the exon 1-encoded region in regulating channel inactivation, which may affect sperm mobility and sperm competition. Our findings provide a necessary foundation for future experimental investigations of Catsper1's function in sperm physiology and role in sperm competition using rodent models.     

Excess of amino acid substitutions relative to polymorphism between X-linked duplications in Drosophila melanogaster

We have obtained sequence polymorphism data from 13 genes belonging to 5 gene families in Drosophila melanogaster where the K(a)/K(s) between copies is greater than 1. Twelve of these 13 loci are X-linked. In general, there is evidence of purifying selection in all families, as inferred both from levels of silent and replacement variation and insertion/deletion variation, suggesting that the loci are likely functional. Shared polymorphisms indicative of gene conversion between paralogs are rare among the X-linked families, in contrast to available data from autosomal duplicates. McDonald-Kreitman tests between duplicates reveal an excess of amino-acid fixations between copies in the X-linked families, suggesting that the divergence between these loci was driven by positive selection. In contrast, available data from autosomal duplicates show a deficit of fixations, consistent with gene conversion being a strong homogenizing force.     

MtArt: a new model of amino acid replacement for Arthropoda

A statistical approach was applied to select those models that best fit each individual mitochondrial (mt) protein at different taxonomic levels of metazoans. The existing mitochondrial replacement matrices, MtREV and MtMam, were found to be the best-fit models for the mt-proteins of vertebrates, with the exception of Nd6, at different taxonomic levels. Remarkably, existing mitochondrial matrices generally failed to best-fit invertebrate mt-proteins. In an attempt to better model the evolution of invertebrate mt-proteins, a new replacement matrix, named MtArt, was constructed based on arthropod mt-proteomes. The new model was found to best fit almost all analyzed invertebrate mt-protein data sets. The observed pattern of model fit across the different data sets indicates that no single replacement matrix is able to describe the general evolutionary properties of mt-proteins but rather that taxonomical biases and/or the existence of different mt-genetic codes have great influence on which model is selected.     

The study of substrate specificity of a new peptide substrate of endothelin-converting enzyme

A new hexapeptide CMC-Ala-Gly-Gly-Trp-Phe-Arg was used as a substrate for assay of endothelin-converting (ECE; EC 3.4.24.71) and angiotensin-converting (ACE; EC 3.4.15.1) enzymes and of neutral endopeptidase (NEP; EC 3.4.24.11). The specific inhibitors lisinopril (for ACE) and thiorphan (for NEP) were used for discrimination between activities of these enzymes.     

Suppressors of RNAi from plant viruses are subject to episodic positive selection

Viral suppressors of RNAi (VSRs) are proteins that actively inhibit the antiviral RNA interference (RNAi) immune response, providing an immune evasion route for viruses. It has been hypothesized that VSRs are engaged in a molecular ‘arms race’ with RNAi pathway genes. Two lines of evidence support this. First, VSRs from plant viruses display high sequence diversity, and are frequently gained and lost over evolutionary time scales. Second, Drosophila antiviral RNAi genes show high rates of adaptive evolution. Here, we investigate whether VSRs diversify faster than other genes and, if so, whether this is a result of positive selection, as might be expected in an arms race. By analysis of 12 plant RNA viruses, we show that the relative rate of protein evolution is higher for VSRs than for other genes, but that this is not attributable to pervasive positive selection. We argue that, because evolutionary time scales are extremely different for viruses and eukaryotes, it is improbable that viral adaptation (as measured by the ratio of non-synonymous to synonymous change) will be dominated by one-to-one coevolution with eukaryotes. Instead, for plant virus VSRs, we find strong evidence of episodic selection—diversifying selection that acts on a subset of lineages—which might be attributable to frequent shifts between different host genotypes or species.     

A single amino acid change in the plant alternative oxidase alters the specificity of organic acid activation.

The alternative oxidase is a quinol oxidase of the respiratory chain of plants and some fungi and protists. Its activity is regulated by redox-sensitive disulphide bond formation between neighbouring subunits and direct interaction with certain alpha-ketoacids. To investigate these regulatory mechanisms, we undertook site-directed mutagenesis of soybean and Arabidopsis alternative oxidase cDNAs, and expressed them in tobacco plants and Escherichia coli, respectively. The homologous C99 and C127 residues of GmAOX3 and AtAOX1a, respectively, were changed to serine. In the plant system, this substitution prevented oxidative inactivation of alternative oxidase and rendered the protein insensitive to pyruvate activation, in agreement with the recent results from other laboratories [Rhoads et al. (1998) J. Biol. Chem. 273, 30750-30756; Vanlerberghe et al. (1998) Plant Cell 10, 1551-1560]. However, the mutated protein is instead activated specifically by succinate. Measurements of AtAOX1a activity in bacterial membranes lacking succinate dehydrogenase confirmed that the stimulation of the mutant protein's activity by succinate did not involve its metabolism. Examples of alternative oxidase proteins with the C to S substitution occur in nature and these oxidases are expected to be activated under most conditions in vivo, with implications for the efficiency of respiration in the tissues which express them.     

Insights from modeling protein evolution with context-dependent mutation and asymmetric amino acid selection

We develop an approximate maximum likelihood method to estimate flanking nucleotide context-dependent mutation rates and amino acid exchange-dependent selection in orthologous protein-coding sequences and use it to analyze genome-wide coding sequence alignments from mammals and yeast. Allowing context-dependent mutation provides a better fit to coding sequence data than simpler (context-independent or CpG "hotspot") models and significantly affects selection parameter estimates. Allowing asymmetric (nonreciprocal) selection on amino acid exchanges gives a better fit than simple dN/dS or symmetric selection models. Relative selection strength estimates from our models show good agreement with independent estimates derived from human disease-causing and engineered mutations. Selection strengths depend on local protein structure, showing expected biophysical trends in helical versus nonhelical regions and increased asymmetry on polar-hydrophobic exchanges with increased burial. The more stringent selection that has previously been observed for highly expressed proteins is primarily concentrated in buried regions, supporting the notion that such proteins are under stronger than average selection for stability. Our analyses indicate that a highly parameterized model of mutation and selection is computationally tractable and is a useful tool for exploring a variety of biological questions concerning protein and coding sequence evolution.     

Kinetics of synonymous codon change for an amino acid of arbitrary degeneracy

The kinetics of synonymous codon change and species divergence is described in a matrix formalism that is equally applicable to all levels of codon degeneracy and all levels of codon or nucleotide bias. Based on the formalism it is possible to calculate the sum of all the synonymous substitution rate constants from the observed sequence differences between two species. This sum, the relaxation rate, is equivalent to the LogDet transformation that has recently been proposed as a new measure of evolutionary distance (Lockhardt et al.Mol. Biol. Evol. 11(4): 605–612, 1994). The relationship between this measure and the average number of base changes per site (K) is discussed. The formalism is tested on some sets of simulated sequence divergence data.     

Positive selection on multiple antique allelic lineages of transferrin in the polyploid Carassius auratus

Transferrin polymorphism has been studied in the polyploid Carassius auratus by cloning and sequence analysis of cDNAs from its three subspecies C. auratus gibelio, C. auratus auratus, and C. auratus cuvieri. DNA polymorphism of extremely high extent was shown for the transferrin gene by the 248 segregation sites among coding region sequences of its alleles. The deduced amino acid sequences of the transferrin alleles showed variable theoretical physicochemical parameters, which might constitute molecular basis for their electrophoretic heterogeneity. Positive selection was inferred by the replacement/synonymous ratios larger than 1 in partial allelic lineages which was subsequently confirmed by likelihood simulation under neutral or selection models. Furthermore, the correspondent sites to these selected codons were collectively located at two planes in the crystallographic structure of rabbit transferrin, which suggested that the rapid evolution of C. auratus transferrin might correlate to its adaptation to variable environmental elements such as oxygen pressure. The minimal 26 recombination events were detected among coding sequences of C. auratus transferrin, with partial mosaic sequences and breakpoints identified by identity scanning and information site analyses. Phylogenetic analyses revealed multiple antique allelic lineages of transferrin, which was estimated to diverge fifteen to twenty MYA. All these features strongly suggested the role of balancing selection in long persistence of high transferrin polymorphism in C. auratus. Furthermore, owing to its particular evolutionary backgrounds, the silver crucian carp might possess a distinctive balancing selection mechanism.     

Positive Darwinian selection in gamete recognition proteins of Strongylocentrotus sea urchins

Gamete recognition proteins commonly experience positive Darwinian selection and evolve more rapidly than nonreproductive proteins, but the selective forces responsible for their adaptive diversification remain unclear. We examined the patterns of positive selection in the cognate interacting pair of proteins formed by sperm bindin and its egg receptor (EBR1) and in two regions of the sea urchin sperm receptor for egg jelly suREJ3 gene (exons 22 and 26) among four species of Strongylocentrotus sea urchins (S. purpuratus, S. droebachiensis, S. pallidus and S. franciscanus). The signatures of selection differed at each reproductive protein. A strong signal of positive selection was detected at bindin in all lineages even though the species compared had highly variable gamete traits and experience different intensities and forms of sexual selection and sexual conflict in nature. Weaker selection was observed at EBR1 but the small region studied precluded a clear understanding of the extent of sexual conflict between bindin and the EBR1 protein. At the suREJ3 locus, diversifying selection was observed in exon 22 but not exon 26, suggesting that these regions experience different selective pressures and evolutionary constraints. Positive selection was also detected within S. pallidus at suREJ‐22 because of the presence of 12 amino acid replacement mutations segregating at frequencies >0.10. Our results suggest that sexual conflict may be the predominant evolutionary mechanism driving the rapid diversification of reproductive proteins between, and polymorphism within, strongylocentrotid sea urchins.     

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.     

Evidence for positive selection on Drosophila melanogaster seminal fluid protease homologs

Proteins present in the seminal fluid of Drosophila melanogaster (accessory gland proteins Acps) contribute to female postmating behavioral changes, sperm storage, sperm competition, and immunity. Consequently, male-female coevolution and host-pathogen interactions are thought to underlie the rapid, adaptive evolution that characterizes several Acp-encoding genes. We propose that seminal fluid proteases are likely targets of selection due to their demonstrated or potential roles in between-sex interactions and immune processes. We use within- and between-species sequence data for 5 predicted protease-encoding Acp loci to test this hypothesis. Our polymorphism-based analyses find evidence for positive selection at 2 genes, both of which encode predicted serine protease homologs. One of these genes, CG6069, also shows evidence for consistent selection on a subset of codons over a deeper evolutionary time scale. The second gene, CG9997, was previously shown to be essential for normal sperm usage, suggesting that sexual selection may underlie its history of adaptation.     

A single amino acid change determines persistence of a chimeric Theiler's virus.

The DA strain of Theiler's virus persists in the central nervous system of mice and causes chronic inflammation and demyelination. On the other hand, the GDVII strain causes an acute encephalitis and does not persist in surviving animals. Series of recombinants between infectious cDNA clones of the genomes of DA and GDVII viruses have been constructed. The analysis of the phenotypes of the recombinant viruses has shown that determinants of persistence and demyelination are present in the capsid proteins of DA virus. Chimeric viruses constructed by the different research groups gave consistent results, with one exception. Chimeras GD1B-2A/DAFL3 and GD1B-2C/DAFL3, which contain part of capsid protein VP2, capsid proteins VP3 and VP1, and different portions of P2 of GDVII in a DA background, were able to persist and cause demyelination. Chimera R4, whose genetic map is identical to that of GD1B-2A/DAFL3, was not. After exchanging the viral chimeras between laboratories and verifying each other's observations, new chimeras were generated in order to explain this difference. Here we report that the discrepancy can be attributed to a single amino acid difference in the sequence of the capsid protein VP2 of the two parental DA strains. DAFL3 (University of Chicago) and the chimeras derived from it, GD1B-2A/DAFL3 and GD1B-2C/DAFL3, contain a Lys at position 141, while TMDA (Institut Pasteur) and R4, the chimera derived from it, contain an Asn in that position. This amino acid is located at the tip of the EF loop, on the rim of the depression spanning the twofold axis of the capsid. These results show that a single amino acid change can confer the ability to persist and demyelinate to a chimeric Theiler's virus, and they pinpoint a region of the viral capsid that is important for this phenotype.