Positive selection in the evolution of cancer

We hypothesize that forms of antagonistic coevolution have forged strong links between positive selection at the molecular level and increased cancer risk. By this hypothesis, evolutionary conflict between males and females, mothers and foetuses, hosts and parasites, and other parties with divergent fitness interests has led to rapid evolution of genetic systems involved in control over fertilization and cellular resources. The genes involved in such systems promote cancer risk as a secondary effect of their roles in antagonistic coevolution, which generates evolutionary disequilibrium and maladaptation. Evidence from two sources: (1) studies on specific genes, including SPANX cancer/testis antigen genes, several Y-linked genes, the pem homebox gene, centromeric histone genes, the breast cancer gene BRCA1, the angiogenesis gene ANG, cadherin genes, cytochrome P450 genes, and viral oncogenes; and (2) large-scale database studies of selection on different functional categories of genes, supports our hypothesis. These results have important implications for understanding the evolutionary underpinnings of cancer and the dynamics of antagonistically-coevolving molecular systems.     

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 selection for elevated gene expression noise in yeast

It is well known that the expression noise is lessened by natural selection for genes that are important for cell growth or are sensitive to dosage. In theory, expression noise can also be elevated by natural selection when noisy gene expression is advantageous. Here we analyze yeast genome‐wide gene expression noise data and show that plasma‐membrane transporters show significantly elevated expression noise after controlling all confounding factors. We propose a model that explains why and under what conditions elevated expression noise may be beneficial and subject to positive selection. Our model predicts and the simulation confirms that, under certain conditions, expression noise also increases the evolvability of gene expression by promoting the fixation of favorable expression level‐altering mutations. Indeed, yeast genes with higher noise show greater between‐strain and between‐species divergences in expression, even when all confounding factors are excluded. Together, our theoretical model and empirical results suggest that, for yeast genes such as plasma‐membrane transporters, elevated expression noise is advantageous, is subject to positive selection, and is a facilitator of adaptive gene expression evolution.     

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.     

Natural selection promotes divergence of transferrin among salmonid species

Transferrin is an iron-binding protein that plays an important role in iron metabolism and resistance to bacterial infection in a variety of organisms. A comparison of transferrin coding sequences from four salmonid species shows that the rate of evolution at nonsynonymous sites is significantly higher than the rate at synonymous sites, suggesting that positive natural selection for new alleles has played an important role in the evolution of transferrin in some salmon species. We hypothesize that the selective agent driving rapid divergence is interactions between host transferrin and the iron-scavenging proteins of pathogenic bacteria.     

Positive selection of IL‐33 in adaptive immunity of domestic Chinese goats

The identification of the candidate genes that play key role in phenotypic variation in livestock populations can provide new information about evolution and positive selection. IL‐33 (71954) (Interleukin) gene is associated with the increased nematode resistance in small ruminants; however, the role of IL‐33 for the genetic control of different diseases in Chinese goat breeds is poorly described in scientific literature. Therefore, the current investigation was performed for the better understanding of the molecular evolution and the positive selection of single‐nucleotide polymorphism in IL‐33 gene. Fixation Index (F_ST)‐based method was used for the outlier loci determination and found that IL‐33 was present in outlier area with the provisional combined allocation of mean heterozygosity and F_ST. Positively selected IL‐33 gene was significantly, that is, p(Simul F_ST < sample F_ST = 0.98*) present in corresponding positive selection area. Hence, our study provided novel information about the nucleotide variations in IL‐33 gene and found to be nonsynonymous which may helpful for the genetic control of diseases by enhancing the immune system in local Chinese goat breeds as well as in other analyzed vertebrate species.     

Modulating the proteinase inhibitory profile of a plant cystatin by single mutations at positively selected amino acid sites

Cysteine proteinase inhibitors of the cystatin superfamily have several important functions in plants, including the inhibition of exogenous cysteine proteinases during herbivory or infection. Here we used a maximum-likelihood approach to assess whether plant cystatins, like other proteins implicated in host-pest interactions, have been subject to positive selection during the course of their evolution. Several amino acid sites were identified as being positively selected in cystatins from either Poaceae (monocots) and Solanaceae (dicots). These hypervariable sites were located at strategic positions on the protein: on each side of the conserved glycine residues in the N-terminal trunk, within the first and second inhibitory loops entering the active site of target enzymes, and surrounding the larfav motif, a sequence of unknown function conserved among plant cystatins. Supporting the assumption that positively selected, hypervariable sites are indicative of amino acid sites implicated in functional diversity, mutants of the 8th cystatin unit of tomato multicystatin including alternative residues at positively selected sites in the N-terminal trunk exhibited highly variable affinities for the cysteine proteases papain, cathepsin B and cathepsin H. Overall, these observations support the hypothesis that plant cystatins have been under selective pressure to evolve in response to predatory challenges by herbivorous enemies. They also indicate the potential of site-directed mutagenesis at positively selected sites for the generation of cystatins with improved binding properties.     

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 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 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 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.     

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.     

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.     

Positive selection drives adaptive diversification of the 4‐coumarate: CoA ligase (4CL) gene in angiosperms

Lignin and flavonoids play a vital role in the adaption of plants to a terrestrial environment. 4‐Coumarate: coenzyme A ligase (4CL) is a key enzyme of general phenylpropanoid metabolism which provides the precursors for both lignin and flavonoids biosynthesis. However, very little is known about how such essential enzymatic functions evolve and diversify. Here, we analyze 4CL sequence variation patterns in a phylogenetic framework to further identify the evolutionary forces that lead to functional divergence. The results reveal that lignin‐biosynthetic 4CLs are under positive selection. The majority of the positively selected sites are located in the substrate‐binding pocket and the catalytic center, indicating that nonsynonymous substitutions might contribute to the functional evolution of 4CLs for lignin biosynthesis. The evolution of 4CLs involved in flavonoid biosynthesis is constrained by purifying selection and maintains the ancestral role of the protein in response to biotic and abiotic factors. Overall, our results demonstrate that protein sequence evolution via positive selection is an important evolutionary force driving adaptive diversification in 4CL proteins in angiosperms. This diversification is associated with adaption to a terrestrial environment.