Lineage-Specific Differences in the Amino Acid Substitution Process

In Darwinian evolution, mutations occur approximately at random in a gene, turned into amino acid mutations by the genetic code. Some mutations are fixed to become substitutions and some are eliminated from the population. Partitioning pairs of closely related species with complete genome sequences by average population size of each pair, we looked at the substitution matrices generated for these partitions and compared the substitution patterns between species. We estimated a population genetic model that relates the relative fixation probabilities of different types of mutations to the selective pressure and population size. Parameterizations of the average and distribution of selective pressures for different amino acid substitution types in different population size comparisons were generated with a Bayesian framework. We found that partitions in population size as well as in substitution type are required to explain the substitution data. Selection coefficients were found to decrease with increasingly radical amino acid substitution and with increasing effective population size.To further explore the role of underlying processes in amino acid substitution, we analyzed embryophyte (plant) gene families from TAED (The Adaptive Evolution Database), where solved structures for at least one member exist in the Protein Data Bank. Using PAML, we assigned branches to three categories: strong negative selection, moderate negative selection/neutrality, and positive diversifying selection. Focusing on the first and third categories, we identified sites changing along gene family lineages and observed the spatial patterns of substitution. Selective sweeps were expected to create primary sequence clustering under positive diversifying selection. Co-evolution through direct physical interaction was expected to cause tertiary structural clustering. Under both positive and negative selection, the substitution patterns were found to be nonrandom. Under positive diversifying selection, significant independent signals were found for primary and tertiary sequence clustering, suggesting roles for both selective sweeps and direct physical interaction. Under strong negative selection, the signals were not found to be independent. All together, a complex interplay of population genetic and protein thermodynamics forces is suggested.     

Hitchhiking both ways: effect of two interfering selective sweeps on linked neutral variation

The neutral polymorphism pattern in the vicinity of a selective sweep can be altered by both stochastic and deterministic factors. Here, we focus on the impact of another selective sweep in the region of influence of a first one. We study the signature left on neutral polymorphism by positive selection at two closely linked loci, when both beneficial mutations reach fixation. We show that, depending on the timing of selective sweeps and on their selection coefficients, the two hitchhiking effects can interfere with each other, leading to less reduction in heterozygosity than a single selective sweep of the same magnitude and more importantly to an excess of intermediate-frequency variants relative to neutrality under some parameter values. This pattern can be sustained and potentially alter the detection of positive selection, including by provoking spurious detection of balancing selection. In situations where positive selection is suspected a priori at several closely linked loci, the polymorphism pattern in the region may also be informative about their selective histories.     

Bayesian analysis suggests that most amino acid replacements in Drosophila are driven by positive selection

One of the principal goals of population genetics is to understand the processes by which genetic variation within species (polymorphism) becomes converted into genetic differences between species (divergence). In this transformation, selective neutrality, near neutrality, and positive selection may each play a role, differing from one gene to the next. Synonymous nucleotide sites are often used as a uniform standard of comparison across genes on the grounds that synonymous sites are subject to relatively weak selective constraints and so may, to a first approximation, be regarded as neutral. Synonymous sites are also interdigitated with nonsynonymous sites and so are affected equally by genomic context and demographic factors. Hence a comparison of levels of polymorphism and divergence between synonymous sites and amino acid replacement sites in a gene is potentially informative about the magnitude of selective forces associated with amino acid replacements. We have analyzed 56 genes in which polymorphism data from D. simulans are compared with divergence from a reference strain of D. melanogaster. The framework of the analysis is Bayesian and assumes that the distribution of selective effects (Malthusian fitnesses) is Gaussian with a mean that differs for each gene. In such a model, the average scaled selection intensity (gamma = N(e)s) of amino acid replacements eligible to become polymorphic or fixed is -7.31, and the standard deviation of selective effects within each locus is 6.79 (assuming homoscedasticity across loci). For newly arising mutations of this type that occur in autosomal or X-linked genes, the average proportion of beneficial mutations is 19.7%. Among the amino acid polymorphisms in the sample, the expected average proportion of beneficial mutations is 47.7%, and among amino acid replacements that become fixed the average proportion of beneficial mutations is 94.3%. The average scaled selection intensity of fixed mutations is +5.1. The presence of positive selection is pervasive with the single exception of kl-5, a Y-linked fertility gene. We find no evidence that a significant fraction of fixed amino acid replacements is neutral or nearly neutral or that positive selection drives amino acid replacements at only a subset of the loci. These results are model dependent and we discuss possible modifications of the model that might allow more neutral and nearly neutral amino acid replacements to be fixed.     

Evolution of the human immunodeficiency virus envelope gene is dominated by purifying selection

The evolution of the human immunodeficiency virus (HIV-1) during chronic infection involves the rapid, continuous turnover of genetic diversity. However, the role of natural selection, relative to random genetic drift, in governing this process is unclear. We tested a stochastic model of genetic drift using partial envelope sequences sampled longitudinally in 28 infected children. In each case the Bayesian posterior (empirical) distribution of coalescent genealogies was estimated using Markov chain Monte Carlo methods. Posterior predictive simulation was then used to generate a null distribution of genealogies assuming neutrality, with the null and empirical distributions compared using four genealogy-based summary statistics sensitive to nonneutral evolution. Because both null and empirical distributions were generated within a coalescent framework, we were able to explicitly account for the confounding influence of demography. From the distribution of corrected P-values across patients, we conclude that empirical genealogies are more asymmetric than expected if evolution is driven by mutation and genetic drift only, with an excess of low-frequency polymorphisms in the population. This indicates that although drift may still play an important role, natural selection has a strong influence on the evolution of HIV-1 envelope. A negative relationship between effective population size and substitution rate indicates that as the efficacy of selection increases, a smaller proportion of mutations approach fixation in the population. This suggests the presence of deleterious mutations. We therefore conclude that intrahost HIV-1 evolution in envelope is dominated by purifying selection against low-frequency deleterious mutations that do not reach fixation.     

Departure from neutrality at the mitochondrial NADH dehydrogenase subunit 2 gene in humans,but not in chimpanzees.

To test whether patterns of mitochondrial DNA (mtDNA) variation are consistent with a neutral model of molecular evolution, nucleotide sequences were determined for the 1041 bp of the NADH dehydrogenase subunit 2 (ND2) gene in 20 geographically diverse humans and 20 common chimpanzees. Contingency tests of neutrality were performed using four mutational categories for the ND2 molecule: synonymous and nonsynonymous mutations in the transmembrane regions, and synonymous and nonsynonymous mutations in the surface regions. The following three topological mutational categories were also used: intraspecific tips, intraspecific interiors, and interspecific fixed differences. The analyses reveal a significantly greater number of nonsynonymous polymorphisms within human transmembrane regions than expected based on interspecific comparisons, and they are inconsistent with a neutral equilibrium model. This pattern of excess nonsynonymous polymorphism is not seen within chimpanzees. Statistical tests of neutrality, such as TAJIMA''s D test, and the D and F tests proposed by FU and LI, indicate an excess of low frequency polymorphisms in the human data, but not in the chimpanzee data. This is consistent with recent directional selection, a population bottleneck or background selection of slightly deleterious mutations in human mtDNA samples. The analyses further support the idea that mitochondrial genome evolution is governed by selective forces that have the potential to affect its use as a "neutral" marker in evolutionary and population genetic studies.     

HIV evolution: CTL escape mutation and reversion after transmission

Within-patient HIV evolution reflects the strong selection pressure driving viral escape from cytotoxic T-lymphocyte (CTL) recognition. Whether this intrapatient accumulation of escape mutations translates into HIV evolution at the population level has not been evaluated. We studied over 300 patients drawn from the B- and C-clade epidemics, focusing on human leukocyte antigen (HLA) alleles HLA-B57 and HLA-B5801, which are associated with long-term HIV control and are therefore likely to exert strong selection pressure on the virus. The CTL response dominating acute infection in HLA-B57/5801-positive subjects drove positive selection of an escape mutation that reverted to wild-type after transmission to HLA-B57/5801-negative individuals. A second escape mutation within the epitope, by contrast, was maintained after transmission. These data show that the process of accumulation of escape mutations within HIV is not inevitable. Complex epitope- and residue-specific selection forces, including CTL-mediated positive selection pressure and virus-mediated purifying selection, operate in tandem to shape HIV evolution at the population level.     

Reversion of CTL escape-variant immunodeficiency viruses in vivo

Engendering cytotoxic T-lymphocyte (CTL) responses is likely to be an important goal of HIV vaccines. However, CTLs select for viral variants that escape immune detection. Maintenance of such escape variants in human populations could pose an obstacle to HIV vaccine development. We first observed that escape mutations in a heterogeneous simian immunodeficiency virus (SIV) isolate were lost upon passage to new animals. We therefore infected macaques with a cloned SIV bearing escape mutations in three immunodominant CTL epitopes, and followed viral evolution after infection. Here we show that each mutant epitope sequence continued to evolve in vivo, often re-establishing the original, CTL-susceptible sequence. We conclude that escape from CTL responses may exact a cost to viral fitness. In the absence of selective pressure upon transmission to new hosts, these original escape mutations can be lost. This suggests that some HIV CTL epitopes will be maintained in human populations.     

Selective sweeps in a 2-locus model for sex-ratio meiotic drive in Drosophila simulans

A way to identify loci subject to positive selection is to detect the signature of selective sweeps in given chromosomal regions. It is revealed by the departure of DNA polymorphism patterns from the neutral equilibrium predicted by coalescent theory. We surveyed DNA sequence variation in a region formerly identified as causing "sex-ratio" meiotic drive in Drosophila simulans. We found evidence that this system evolved by positive selection at 2 neighboring loci, which thus appear to be required simultaneously for meiotic drive to occur. The 2 regions are approximately 150-kb distant, corresponding to a genetic distance of 0.1 cM. The presumably large transmission advantage of chromosomes carrying meiotic drive alleles at both loci has not erased the individual signature of selection at each locus. This chromosome fragment combines a high level of linkage disequilibrium between the 2 critical regions with a high recombination rate. As a result, 2 characteristic traits of selective sweeps--the reduction of variation and the departure from selective neutrality in haplotype tests--show a bimodal pattern. Linkage disequilibrium level indicates that, in the natural population from Madagascar used in this study, the selective sweep may be as recent as 100 years.     

Microsatellite variation and evolution of human lactase persistence

The levels of haplotype diversity within the lineages defined by two single-nucleotide polymorphisms (SNPs) (–13910 C/T and –22018 G/A) associated with human lactase persistence were assessed with four fast-evolving microsatellite loci in 794 chromosomes from Portugal, Italy, Fulbe from Cameroon, São Tomé and Mozambique. Age estimates based on the intraallelic microsatellite variation indicate that the –13910*T allele, which is more tightly associated with lactase persistence, originated in Eurasia before the Neolithic and after the emergence of modern humans outside Africa. We detected significant departures from neutrality for the –13910*T variant in geographically and evolutionary distant populations from southern Europe (Portuguese and Italians) and Africa (Fulbe) by using a neutrality test based on the congruence between the frequency of the allele and the levels of intraallelic variability measured by the number of mutations in adjacent microsatellites. This result supports the role of selection in the evolution of lactase persistence, ruling out possible confounding effects from recombination suppression and population history. Reevaluation of the available evidence on variation of the –13910 and –22018 loci indicates that lactase persistence probably originated from different mutations in Europe and most of Africa, even if 13910*T is not the causal allele, suggesting that selective pressure could have promoted the convergent evolution of the trait. Our study shows that a limited number of microsatellite loci may provide sufficient resolution to reconstruct key aspects of the evolutionary history of lactase persistence, providing an alternative to approaches based on large numbers of SNPs.Electronic supplementary material Supplementary material is available for this article at     

Intra- and interspecific variation of the CCR5 gene in higher primates

We have evaluated the molecular evolution of the chemokine receptor CCR5 in primates. The chemokine receptor CCR5 serves as a major co-receptor for human immunodeficiency virus/simian immunodeficiency virus (HIV/SIV) infection. Knowledge of evolution of the CCR5 molecule and selection on the CCR5 gene may shed light on its functional role. The comparison of differences between intraspecific polymorphisms and interspecific fixed substitutions provides useful information regarding modes of selection during the course of evolution. There is marked polymorphism in the CCR5 gene sequence within different primate species, whereas sequence divergence between different species is small. By using contingency tests, we compared synonymous (SS) and nonsynonymous (NS) CCR5 mutations occurring within and between a broad range of primates. Our results demonstrate that CCR5 evolution did not follow expectations of strict neutrality at the level of the whole gene. The proportion of NS to SS at the intraspecific level was significantly higher than that observed at the interspecific level. These results suggest that most CCR5 NS polymorphisms are slightly deleterious. However, at domains more closely correlated with its known biological functions, there was no obvious evidence to support deviation from neutrality.     

Learning Natural Selection from the Site Frequency Spectrum

Genetic adaptation to external stimuli occurs through the combined action of mutation and selection. A central problem in genetics is to identify loci responsive to specific selective constraints. Many tests have been proposed to identify the genomic signatures of natural selection by quantifying the skew in the site frequency spectrum (SFS) under selection relative to neutrality. We build upon recent work that connects many of these tests under a common framework, by describing how selective sweeps affect the scaled SFS. We show that the specific skew depends on many attributes of the sweep, including the selection coefficient and the time under selection. Using supervised learning on extensive simulated data, we characterize the features of the scaled SFS that best separate different types of selective sweeps from neutrality. We develop a test, SFselect, that consistently outperforms many existing tests over a wide range of selective sweeps. We apply SFselect to polymorphism data from a laboratory evolution experiment of Drosophila melanogaster adapted to hypoxia and identify loci that strengthen the role of the Notch pathway in hypoxia tolerance, but were missed by previous approaches. We further apply our test to human data and identify regions that are in agreement with earlier studies, as well as many novel regions.     

Contrasting patterns of nonneutral evolution in proteins encoded in nuclear and mitochondrial genomes

We report that patterns of nonneutral DNA sequence evolution among published nuclear and mitochondrially encoded protein-coding loci differ significantly in animals. Whereas an apparent excess of amino acid polymorphism is seen in most (25/31) mitochondrial genes, this pattern is seen in fewer than half (15/36) of the nuclear data sets. This differentiation is even greater among data sets with significant departures from neutrality (14/15 vs. 1/6). Using forward simulations, we examined patterns of nonneutral evolution using parameters chosen to mimic the differences between mitochondrial and nuclear genetics (we varied recombination rate, population size, mutation rate, selective dominance, and intensity of germ line bottleneck). Patterns of evolution were correlated only with effective population size and strength of selection, and no single genetic factor explains the empirical contrast in patterns. We further report that in Arabidopsis thaliana, a highly self-fertilizing plant with effectively low recombination, five of six published nuclear data sets also exhibit an excess of amino acid polymorphism. We suggest that the contrast between nuclear and mitochondrial nonneutrality in animals stems from differences in rates of recombination in conjunction with a distribution of selective effects. If the majority of mutations segregating in populations are deleterious, high linkage may hinder the spread of the occasional beneficial mutation.     

Recurrent positive selection at bgcn, a key determinant of germ line differentiation, does not appear to be driven by simple coevolution with its partner protein bam

Surveys of nucleotide sequence polymorphism in Drosophila melanogaster and Drosophila simulans were performed at 2 interacting loci crucial for gametogenesis: bag-of-marbles (bam) and benign gonial cell neoplasm (bgcn). At the polymorphism level, both loci appear to be evolving under the expectations of the neutral theory. However, ratios of polymorphism and divergence for synonymous and nonsynonymous mutations depart significantly from neutral expectations for both loci consistent with a previous observation of positive selection at bam. The deviations suggest either an excess of synonymous polymorphisms or an excess of nonsynonymous fixations at both loci. Synonymous evolution appears to conform to neutrality at bam. At bgcn, there is evidence of positive selection affecting preferred synonymous mutations along the D. simulans lineage. However, there is also a significantly higher rate of nonsynonymous fixations at bgcn within D. simulans. Thus, the deviation from neutrality detected by the McDonald-Kreitman test at these 2 loci is likely due to the selective acceleration of nonsynonymous fixations. Differences in the pattern of amino acid fixations between these 2 interacting proteins suggest that the detected positive selection is not due to a simple model of coevolution.     

Robustness, evolvability, and neutrality

Biological systems, from macromolecules to whole organisms, are robust if they continue to function, survive, or reproduce when faced with mutations, environmental change, and internal noise. I focus here on biological systems that are robust to mutations and ask whether such systems are more or less evolvable, in the sense that they can acquire novel properties. The more robust a system is, the more mutations in it are neutral, that is, without phenotypic effect. I argue here that such neutral change--and thus robustness--can be a key to future evolutionary innovation, if one accepts that neutrality is not an essential feature of a mutation. That is, a once neutral mutation may cause phenotypic effects in a changed environment or genetic background. I argue that most, if not all, neutral mutations are of this sort, and that the essentialist notion of neutrality should be abandoned. This perspective reconciles two opposing views on the forces dominating organismal evolution, natural selection and random drift: neutral mutations occur and are especially abundant in robust systems, but they do not remain neutral indefinitely, and eventually become visible to natural selection, where some of them lead to evolutionary innovations.     

Simulated climate change provokes rapid genetic change in the Mediterranean shrub Fumana thymifolia

Rapid climate change will impose strong directional selection pressures on natural plant populations. Climate-linked genetic variation in natural populations indicates that an evolutionary response is possible. We investigated such a response by comparing individuals subjected to elevated drought and warming treatments with individuals establishing in an unmanipulated climate within the same population. We report that reduction in seedling establishment in response to climate manipulations is nonrandom and results from the selection pressure imposed by artificially warmed and droughted conditions. When compared against control samples, high single-locus genetic divergence occurred in drought and warming treatment samples, with genetic differentiation up to 37 times higher than background (mean neutral locus) genetic differentiation. These loci violate assumptions of selective neutrality, indicating the signature of natural selection by drought. Our results demonstrate that rapid evolution in response to climate change may be widespread in natural populations, based on genetic variation already present within the population.     

Theoretical Study of near Neutrality. I. Heterozygosity and Rate of Mutant Substitution

In order to clarify the nature of "near neutrality" in molecular evolution and polymorphism, extensive simulation studies were performed. Selection coefficients of new mutations are assumed to be small so that both random genetic drift and selection contribute to determining the behavior of mutants. The model also incorporates normally distributed spatial fluctuation of selection coefficients. If the system starts from "average neutrality," it will move to a better adapted state, and most new mutations will become "slightly deleterious." Monte Carlo simulations have indicated that such adaptation is attained, but that the rate of such "progress" is very low for weak selection. In general, the larger the population size, the more effective the selection becomes. Also, as selection becomes weaker, the behavior of the mutants approaches that of completely neutral genes. Thus, the weaker the selection, the smaller is the effect of population size on mutant dynamics. Increase of heterozygosity with population size is very pronounced for subdivided populations. The significance of these results is discussed in relation to various observed facts on molecular evolution and polymorphism, such as generation-time dependency and overdispersion of the molecular clock, or contrasting patterns of DNA and protein polymorphism among some closely related species.     

Positive and purifying selection influence the evolution of doublesex in the Anastrepha fraterculus species group

The gene doublesex (dsx) is considered to be under strong selective constraint along its evolutionary history because of its central role in somatic sex differentiation in insects. However, previous studies of dsx used global estimates of evolutionary rates to investigate its molecular evolution, which potentially miss signals of adaptive changes in generally conserved genes. In this work, we investigated the molecular evolution of dsx in the Anastrepha fraterculus species group (Diptera, Tephritidae), and test the hypothesis that this gene evolved solely by purifying selection using divergence-based and population-based methods. In the first approach, we compared sequences from Anastrepha and other Tephritidae with other Muscomorpha species, analyzed variation in nonsynonymous to synonymous rate ratios (dN/dS) in the Tephritidae, and investigated radical and conservative changes in amino acid physicochemical properties. We show a general selective constraint on dsx, but with signs of positive selection mainly in the common region. Such changes were localized in alpha-helices previously reported to be involved in dimer formation in the OD2 domain and near the C-terminal of the OD1 domain. In the population-based approach, we amplified a region of 540 bp that spanned almost all of the region common to both sexes from 32 different sites in Brazil. We investigated patterns of selection using neutrality tests based on the frequency spectrum and locations of synonymous and nonsynonymous mutations in a haplotype network. As in the divergence-based approach, these analyses showed that dsx has evolved under an overall selective constraint, but with some events of positive selection. In contrast to previous studies, our analyses indicate that even though dsx has indeed evolved as a conserved gene, the common region of dsx has also experienced bouts of positive selection, perhaps driven by sexual selection, during its evolution.     

The genealogy of a sequence subject to purifying selection at multiple sites

We investigate the effect of purifying selection at multiple sites on both the shape of the genealogy and the distribution of mutations on the tree. We find that the primary effect of purifying selection on a genealogy is to shift the distribution of mutations on the tree, whereas the shape of the tree remains largely unchanged. This result is relevant to the large number of coalescent estimation procedures, which generally assume neutrality for segregating polymorphisms--applying these estimators to evolutionarily constrained sequences could lead to a significant degree of bias. We also estimate the statistical power of several neutrality tests in detecting weak to moderate purifying selection and find that the power is quite good for some parameter combinations. This result contrasts with previous studies, which predicted low statistical power because of the minor effect that weak purifying selection has on the shape of a genealogy. Finally, we investigate the effect of Hill-Robertson interference among linked deleterious mutations on patterns of molecular variation. We find that dependence among selected loci can substantially reduce the efficacy of even fairly strong purifying selection.     

Large number of replacement polymorphisms in rapidly evolving genes of Drosophila. Implications for genome-wide surveys of DNA polymorphism

We present a survey of nucleotide polymorphism of three novel, rapidly evolving genes in populations of Drosophila melanogaster and D. simulans. Levels of silent polymorphism are comparable to other loci, but the number of replacement polymorphisms is higher than that in most other genes surveyed in D. melanogaster and D. simulans. Tests of neutrality fail to reject neutral evolution with one exception. This concerns a gene located in a region of high recombination rate in D. simulans and in a region of low recombination rate in D. melanogaster, due to an inversion. In the latter case it shows a very low number of polymorphisms, presumably due to selective sweeps in the region. Patterns of nucleotide polymorphism suggest that most substitutions are neutral or nearly neutral and that weak (positive and purifying) selection plays a significant role in the evolution of these genes. At all three loci, purifying selection of slightly deleterious replacement mutations appears to be more efficient in D. simulans than in D. melanogaster, presumably due to different effective population sizes. Our analysis suggests that current knowledge about genome-wide patterns of nucleotide polymorphism is far from complete with respect to the types and range of nucleotide substitutions and that further analysis of differences between local populations will be required to understand the forces more completely. We note that rapidly diverging and nearly neutrally evolving genes cannot be expected only in the genome of Drosophila, but are likely to occur in large numbers also in other organisms and that their function and evolution are little understood so far.