The efficiency of purifying selection in Mammals vs. Drosophila for metabolic genes

The nearly neutral theory of molecular evolution states that the efficiency of natural selection depends on the effective population size. By using a wide range of multispecies data on nucleotide polymorphism, we have tried to ascertain whether there are any differences in the level of selective constraints of metabolic process genes between Mammals and Drosophila species. The results are consistent with a higher selective constraint in Drosophila than in Mammals, according to the expected under the nearly neutral model: purifying selection seems to be more efficient in species with a larger effective population size.     

Mating systems and the efficacy of selection at the molecular level

Mating systems are thought to play a key role in molecular evolution through their effects on effective population size (N(e)) and effective recombination rate. Because of reduced N(e), selection in self-fertilizing species is supposed to be less efficient, allowing fixation of weakly deleterious alleles or lowering adaptation, which may jeopardize their long-term evolution. Relaxed selection pressures in selfers should be detectable at the molecular level through the analyses of the ratio of nonsynonymous and synonymous divergence, D(n)/D(s), or the ratio of nonsynonymous and synonymous polymorphism, pi(n)/pi(s). On the other hand, selfing reveals recessive alleles to selection (homozygosity effect), which may counterbalance the reduction in N(e). Through population genetics models, this study investigates which process may prevail in natural populations and which conditions are necessary to detect evidence for relaxed selection signature at the molecular level in selfers. Under a wide range of plausible population and mutation parameters, relaxed selection against deleterious mutations should be detectable, but the differences between the two mating systems can be weak. At equilibrium, differences between outcrossers and selfers should be more pronounced using divergence measures (D(n)/D(s) ratio) than using polymorphism data (pi(n)/pi(s) ratio). The difference in adaptive substitution rates between outcrossers and selfers is much less predictable because it critically depends on the dominance levels of new advantageous mutations, which are poorly known. Different ways of testing these predictions are suggested, and implications of these results for the evolution of self-fertilizing species are also discussed.     

Population genomics of eusocial insects: the costs of a vertebrate‐like effective population size

The evolution of reproductive division of labour and social life in social insects has lead to the emergence of several life‐history traits and adaptations typical of larger organisms: social insect colonies can reach masses of several kilograms, they start reproducing only when they are several years old, and can live for decades. These features and the monopolization of reproduction by only one or few individuals in a colony should affect molecular evolution by reducing the effective population size. We tested this prediction by analysing genome‐wide patterns of coding sequence polymorphism and divergence in eusocial vs. noneusocial insects based on newly generated RNA‐seq data. We report very low amounts of genetic polymorphism and an elevated ratio of nonsynonymous to synonymous changes – a marker of the effective population size – in four distinct species of eusocial insects, which were more similar to vertebrates than to solitary insects regarding molecular evolutionary processes. Moreover, the ratio of nonsynonymous to synonymous substitutions was positively correlated with the level of social complexity across ant species. These results are fully consistent with the hypothesis of a reduced effective population size and an increased genetic load in eusocial insects, indicating that the evolution of social life has important consequences at both the genomic and population levels.     

Comparative genomics and the study of evolution by natural selection

Genomics profoundly affects most areas of biology, including ecology and evolutionary biology. By examining genome sequences from multiple species, comparative genomics offers new insight into genome evolution and the way natural selection moulds DNA sequence evolution. Functional divergence, as manifested in the accumulation of nonsynonymous substitutions in protein-coding genes, differs among lineages in a manner seemingly related to population size. For example, the ratio of nonsynonymous to synonymous substitution (d_N/d_S) is higher in apes than in rodents, compatible with Ohta's nearly neutral theory of molecular evolution, which suggests that the fixation of slightly deleterious mutations contributes to protein evolution at an extent negatively correlated with effective population size. While this supports the idea that functional evolution is not necessarily adaptive, comparative genomics is uncovering a role for positive Darwinian selection in 10–40% of all genes in different lineages, estimates that are likely to increase when the addition of more genomes gives increased power. Again, population size seems to matter also in this context, with a higher proportion of fixed amino acid changes representing advantageous mutations in large populations. Genes that are particularly prone to be driven by positive selection include those involved with reproduction, immune response, sensory perception and apoptosis. Genetic innovations are also frequently obtained by the gain or loss of complete gene sequences. Moreover, it is increasingly realized, from comparative genomics, that purifying selection conserves much more than just the protein-coding part of the genome, and this points at an important role for regulatory elements in trait evolution. Finally, genome sequencing using outbred or multiple individuals has provided a wealth of polymorphism data that gives information on population history, demography and marker evolution.     

Studying Patterns of Recent Evolution at Synonymous Sites and Intronic Sites in Drosophila melanogaster

Most previous studies of the evolution of codon usage bias (CUB) and intronic GC content (iGC) in Drosophila melanogaster were based on between-species comparisons, reflecting long-term evolutionary events. However, a complete picture of the evolution of CUB and iGC cannot be drawn without knowledge of their more recent evolutionary history. Here, we used a polymorphism dataset collected from Zimbabwe to study patterns of the recent evolution of CUB and iGC. Analyzing coding and intronic data jointly with a model which can simultaneously estimate selection, mutational, and demographic parameters, we have found that: (1) natural selection is probably acting on synonymous codons; (2) a constant population size model seems to be sufficient to explain most of the observed synonymous polymorphism patterns; (3) GC is favored over AT in introns. In agreement with the long-term evolutionary patterns, ongoing selection acting on X-linked synonymous codons is stronger than that acting on autosomal codons. The selective differences between preferred and unpreferred codons tend to be greater than the differences between GC and AT in introns, suggesting that natural selection, not just biased gene conversion, may have influenced the evolution of CUB. Interestingly, evidence for non-equilibrium evolution comes exclusively from the intronic data. However, three different models, an equilibrium model with two classes of selected sites and two non-equilibrium models with changes in either population size or mutational parameters, fit the intronic data equally well. These results show that using inadequate selection (or demographic) models can result in incorrect estimates of demographic (or selection) parameters.     

Inferring Weak Selection from Patterns of Polymorphism and Divergence at ``silent' Sites in Drosophila DNA

Patterns of codon usage and ``silent'''' DNA divergence suggest that natural selection discriminates among synonymous codons in Drosophila. ``Preferred'''' codons are consistently found in higher frequencies within their synonymous families in Drosophila melanogaster genes. This suggests a simple model of silent DNA evolution where natural selection favors mutations from unpreferred to preferred codons (preferred changes). Changes in the opposite direction, from preferred to unpreferred synonymous codons (unpreferred changes), are selected against. Here, selection on synonymous DNA mutations is investigated by comparing the evolutionary dynamics of these two categories of silent DNA changes. Sequences from outgroups are used to determine the direction of synonymous DNA changes within and between D. melanogaster and Drosophila simulans for five genes. Population genetics theory shows that differences in the fitness effect of mutations can be inferred from the comparison of ratios of polymorphism to divergence. Unpreferred changes show a significantly higher ratio of polymorphism to divergence than preferred changes in the D. simulans lineage, confirming the action of selection at silent sites. An excess of unpreferred fixations in 28 genes suggests a relaxation of selection on synonymous mutations in D. melanogaster. Estimates of selection coefficients for synonymous mutations (3.6 <|N(e)s| < 1.3) in D. simulans are consistent with the reduced efficacy of natural selection (|N(e)s| < 1) in the three- to sixfold smaller effective population size of D. melanogaster. Synonymous DNA changes appear to be a prevalent class of weakly selected mutations in Drosophila.     

The neutral theory is dead. Long live the neutral theory

The neutral theory of molecular evolution has been instrumental in organizing our thinking about the nature of evolutionary forces shaping variation at the DNA level. More importantly, it has provided empiricists with a strong set of testable predictions and hence, a useful null hypothesis against which to test for the presence of selection. Evidence indicates that the neutral theory cannot explain key features of protein evolution nor patterns of biased codon usage in certain species. Whereas we now have a reasonable model of selection acting on synonymous changes in Drosophila, protein evolution remains poorly understood. Despite limitations in the applicability of the neutral theory, it is likely to remain an integral part of the quest to understand molecular evolution.     

Evidence for widespread positive and purifying selection across the European rabbit (Oryctolagus cuniculus) genome

The nearly neutral theory of molecular evolution predicts that the efficacy of both positive and purifying selection is a function of the long-term effective population size (N(e)) of a species. Under this theory, the efficacy of natural selection should increase with N(e). Here, we tested this simple prediction by surveying ~1.5 to 1.8 Mb of protein coding sequence in the two subspecies of the European rabbit (Oryctolagus cuniculus algirus and O. c. cuniculus), a mammal species characterized by high levels of nucleotide diversity and N(e) estimates for each subspecies on the order of 1 × 10(6). When the segregation of slightly deleterious mutations and demographic effects were taken into account, we inferred that >60% of amino acid substitutions on the autosomes were driven to fixation by positive selection. Moreover, we inferred that a small fraction of new amino acid mutations (<4%) are effectively neutral (defined as 0 < N(e)s < 1) and that this fraction was negatively correlated with a gene's expression level. Consistent with models of recurrent adaptive evolution, we detected a negative correlation between levels of synonymous site polymorphism and the rate of protein evolution, although the correlation was weak and nonsignificant. No systematic X chromosome-autosome difference was found in the efficacy of selection. For example, the proportion of adaptive substitutions was significantly higher on the X chromosome compared with the autosomes in O. c. algirus but not in O. c. cuniculus. Our findings support widespread positive and purifying selection in rabbits and add to a growing list of examples suggesting that differences in N(e) among taxa play a substantial role in determining rates and patterns of protein evolution.     

Effective population size and the rate and pattern of nucleotide substitutions

Both the overall rate of nucleotide substitution and the relative proportions of synonymous and non-synonymous substitutions are predicted to vary between species that differ in effective population size (N_e). Our understanding of the genetic processes underlying these lineage-specific differences in molecular evolution is still developing. Empirical analyses indicate that variation in substitution rates and patterns caused by differences in N_e is often substantial, however, and must be accounted for in analyses of molecular evolution.     

Adaptive Protein Evolution in Animals and the Effective Population Size Hypothesis

The rate at which genomes adapt to environmental changes and the prevalence of adaptive processes in molecular evolution are two controversial issues in current evolutionary genetics. Previous attempts to quantify the genome-wide rate of adaptation through amino-acid substitution have revealed a surprising diversity of patterns, with some species (e.g. Drosophila) experiencing a very high adaptive rate, while other (e.g. humans) are dominated by nearly-neutral processes. It has been suggested that this discrepancy reflects between-species differences in effective population size. Published studies, however, were mainly focused on model organisms, and relied on disparate data sets and methodologies, so that an overview of the prevalence of adaptive protein evolution in nature is currently lacking. Here we extend existing estimators of the amino-acid adaptive rate by explicitly modelling the effect of favourable mutations on non-synonymous polymorphism patterns, and we apply these methods to a newly-built, homogeneous data set of 44 non-model animal species pairs. Data analysis uncovers a major contribution of adaptive evolution to the amino-acid substitution process across all major metazoan phyla—with the notable exception of humans and primates. The proportion of adaptive amino-acid substitution is found to be positively correlated to species effective population size. This relationship, however, appears to be primarily driven by a decreased rate of nearly-neutral amino-acid substitution because of more efficient purifying selection in large populations. Our results reveal that adaptive processes dominate the evolution of proteins in most animal species, but do not corroborate the hypothesis that adaptive substitutions accumulate at a faster rate in large populations. Implications regarding the factors influencing the rate of adaptive evolution and positive selection detection in humans vs. other organisms are discussed.     

The Genomic Selfing Syndrome Accompanies the Evolutionary Breakdown of Heterostyly

The evolutionary transition from outcrossing to selfing can have important genomic consequences. Decreased effective population size and the reduced efficacy of selection are predicted to play an important role in the molecular evolution of the genomes of selfing species. We investigated evidence for molecular signatures of the genomic selfing syndrome using 66 species of Primula including distylous (outcrossing) and derived homostylous (selfing) taxa. We complemented our comparative analysis with a microevolutionary study of P. chungensis, which is polymorphic for mating system and consists of both distylous and homostylous populations. We generated chloroplast and nuclear genomic data sets for distylous, homostylous, and distylous–homostylous species and identified patterns of nonsynonymous to synonymous divergence (d_N/d_S) and polymorphism (π_N/π_S) in species or lineages with contrasting mating systems. Our analysis of coding sequence divergence and polymorphism detected strongly reduced genetic diversity and heterozygosity, decreased efficacy of purifying selection, purging of large-effect deleterious mutations, and lower rates of adaptive evolution in samples from homostylous compared with distylous populations, consistent with theoretical expectations of the genomic selfing syndrome. Our results demonstrate that self-fertilization is a major driver of molecular evolutionary processes with genomic signatures of selfing evident in both old and relatively young homostylous populations.     

Patterns of protein evolution in Tetrahymena thermophila: implications for estimates of effective population size

High levels of synonymous substitutions among alleles of the surface antigen SerH led to the hypothesis that Tetrahymena thermophila has a tremendously large effective population size, one that is greater than estimated for many prokaryotes (Lynch, M., and J. S. Conery. 2003. Science 302:1401-1404.). Here we show that SerH is unusual as there are substantially lower levels of synonymous variation at five additional loci (four nuclear and one mitochondrial) characterized from T. thermophila populations. Hence, the effective population size of T. thermophila, a model single-celled eukaryote, is lower and more consistent with estimates from other microbial eukaryotes. Moreover, reanalysis of SerH polymorphism data indicates that this protein evolves through a combination of vertical transmission of alleles and concerted evolution of repeat units within alleles. SerH may be under balancing selection due to a mechanism analogous to the maintenance of antigenic variation in vertebrate immune systems. Finally, the dual nature of ciliate genomes and particularly the amitotic divisions of processed macronuclear genomes may make it difficult to estimate accurately effective population size from synonymous polymorphisms. This is because selection and drift operate on processed chromosomes in macronuclei, where assortment of alleles, disruption of linkage groups, and recombination can alter the genetic landscape relative to more canonical eukaryotic genomes.     

Codon usage bias: causative factors,quantification methods and genome‐wide patterns: with emphasis on insect genomes

Codon usage bias refers to the phenomenon where specific codons are used more often than other synonymous codons during translation of genes, the extent of which varies within and among species. Molecular evolutionary investigations suggest that codon bias is manifested as a result of balance between mutational and translational selection of such genes and that this phenomenon is widespread across species and may contribute to genome evolution in a significant manner. With the advent of whole‐genome sequencing of numerous species, both prokaryotes and eukaryotes, genome‐wide patterns of codon bias are emerging in different organisms. Various factors such as expression level, GC content, recombination rates, RNA stability, codon position, gene length and others (including environmental stress and population size) can influence codon usage bias within and among species. Moreover, there has been a continuous quest towards developing new concepts and tools to measure the extent of codon usage bias of genes. In this review, we outline the fundamental concepts of evolution of the genetic code, discuss various factors that may influence biased usage of synonymous codons and then outline different principles and methods of measurement of codon usage bias. Finally, we discuss selected studies performed using whole‐genome sequences of different insect species to show how codon bias patterns vary within and among genomes. We conclude with generalized remarks on specific emerging aspects of codon bias studies and highlight the recent explosion of genome‐sequencing efforts on arthropods (such as twelve Drosophila species, species of ants, honeybee, Nasonia and Anopheles mosquitoes as well as the recent launch of a genome‐sequencing project involving 5000 insects and other arthropods) that may help us to understand better the evolution of codon bias and its biological significance.     

The global impact of Wolbachia on mitochondrial diversity and evolution

The spread of maternally inherited microorganisms, such as Wolbachia bacteria, can induce indirect selective sweeps on host mitochondria, to which they are linked within the cytoplasm. The resulting reduction in effective population size might lead to smaller mitochondrial diversity and reduced efficiency of natural selection. While documented in several host species, it is currently unclear if such a scenario is common enough to globally impact the diversity and evolution of mitochondria in Wolbachia‐infected lineages. Here, we address this question using a mapping of Wolbachia acquisition/extinction events on a large mitochondrial DNA tree, including over 1000 species. Our analyses indicate that on a large phylogenetic scale, other sources of variation, such as mutation rates, tend to hide the effects of Wolbachia. However, paired comparisons between closely related infected and uninfected taxa reveal that Wolbachia is associated with a twofold reduction in silent mitochondrial polymorphism, and a 13% increase in nonsynonymous substitution rates. These findings validate the conjecture that the widespread distribution of Wolbachia infections throughout arthropods impacts the effective population size of mitochondria. These effects might in part explain the disconnection between genetic diversity and demographic population size in mitochondria, and also fuel red‐queen‐like cytonuclear co‐evolution through the fixation of deleterious mitochondrial alleles.     

Multilocus analyses reveal little evidence for lineage‐wide adaptive evolution within major clades of soft pines (Pinus subgenus Strobus)

Estimates from molecular data for the fraction of new nonsynonymous mutations that are adaptive vary strongly across plant species. Much of this variation is due to differences in life history strategies as they influence the effective population size (N_e). Ample variation for these estimates, however, remains even when comparisons are made across species with similar values of N_e. An open question thus remains as to why the large disparity for estimates of adaptive evolution exists among plant species. Here, we have estimated the distribution of deleterious fitness effects (DFE) and the fraction of adaptive nonsynonymous substitutions (α) for 11 species of soft pines (subgenus Strobus) using DNA sequence data from 167 orthologous nuclear gene fragments. Most newly arising nonsynonymous mutations were inferred to be so strongly deleterious that they would rarely become fixed. Little evidence for long‐term adaptive evolution was detected, as all 11 estimates for α were not significantly different from zero. Nucleotide diversity at synonymous sites, moreover, was strongly correlated with attributes of the DFE across species, thus illustrating a strong consistency with the expectations from the Nearly Neutral Theory of molecular evolution. Application of these patterns to genome‐wide expectations for these species, however, was difficult as the loci chosen for the analysis were a biased set of conserved loci, which greatly influenced the estimates of the DFE and α. This implies that genome‐wide parameter estimates will need truly genome‐wide data, so that many of the existing patterns documented previously for forest trees (e.g. little evidence for signature of selection) may need revision.     

Gene expression and protein length influence codon usage and rates of sequence evolution in Populus tremula

Codon bias is generally thought to be determined by a balance between mutation, genetic drift, and natural selection on translational efficiency. However, natural selection on codon usage is considered to be a weak evolutionary force and selection on codon usage is expected to be strongest in species with large effective population sizes. In this paper, I study associations between codon usage, gene expression, and molecular evolution at synonymous and nonsynonymous sites in the long-lived, woody perennial plant Populus tremula (Salicaceae). Using expression data for 558 genes derived from expressed sequence tags (EST) libraries from 19 different tissues and developmental stages, I study how gene expression levels within single tissues as well as across tissues affect codon usage and rates sequence evolution at synonymous and nonsynonymous sites. I show that gene expression have direct effects on both codon usage and the level of selective constraint of proteins in P. tremula, although in different ways. Codon usage genes is primarily determined by how highly expressed a genes is, whereas rates of sequence evolution are primarily determined by how widely expressed genes are. In addition to the effects of gene expression, protein length appear to be an important factor influencing virtually all aspects of molecular evolution in P. tremula.     

FOUNDER EFFECTS AND THE RATE OF MITOCHONDRIAL DNA EVOLUTION IN HAWAIIAN DROSOPHILA

The nearly-neutral-mutation theory predicts that populations with small effective sizes will undergo more rapid molecular evolution than populations with very large effective sizes. In particular, Ohta (1976) predicted that populations of Hawaiian Drosophila that are known to have small population sizes and to experience repeated population bottlenecks due to founder events should show a more rapid rate of molecular evolution relative to the rate of molecular evolution of species with large population sizes such as the continental Drosophila. In this paper we test this prediction by comparing the rate of molecular evolution in two closely related lineages of Hawaiian Drosophila that have experienced very different evolutionary histories. Both lineages belong to the planitibia subgroup of Hawaiian Drosophila. The beta lineage (which includes D. silvestris, D. planitibia, D. differens, and D. hemipeza) has undergone repeated founder events, as evidenced by their geographic distribution and behavioral biology. On the other hand, evidence on geographic distribution and behavior indicates that the alpha lineage (which includes D. melanocephala, D. cyrtaloma, and D. neoperkinsi) has arisen from large ancestral populations without founder effects. The mitochondrial DNA data reveal that, within a lineage, the rate of molecular evolution is rather uniform, while all comparisons between the two lineages show that the rate of molecular evolution in the beta lineage is three times that of the alpha lineage. This analysis strongly supports the predictions made by Ohta.     

Impact of mating systems on patterns of sequence polymorphism in flowering plants

A fundamental challenge in population genetics and molecular evolution is to understand the forces shaping the patterns of genetic diversity within and among species. Among them, mating systems are thought to have important influences on molecular diversity and genome evolution. Selfing is expected to reduce effective population size, Ne, and effective recombination rates, directly leading to reduced polymorphism and increased linkage disequilibrium compared with outcrossing. Increased isolation between populations also results directly from selfing or indirectly from evolutionary changes, such as small flowers and low pollen output, leading to greater differentiation of molecular markers than under outcrossing. The lower effective recombination rate increases the likelihood of hitch-hiking, further reducing within-deme diversity of selfers and thus increasing their genetic differentiation. There are also indirect effects on molecular evolutionary processes. Low Ne reduces the efficacy of selection; in selfers, selection should thus be less efficient in removing deleterious mutations. The rarity of heterozygous sites in selfers leads to infrequent action of biased conversion towards GC, which tends to increase sequences' GC content in the most highly recombining genome regions of outcrossers. To test these predictions in plants, we used a newly developed sequence polymorphism database to investigate the effects of mating system differences on sequence polymorphism and genome evolution in a wide set of plant species. We also took into account other life-history traits, including life form (whether annual or perennial herbs, and woody perennial) and the modes of pollination and seed dispersal, which are known to affect enzyme and DNA marker polymorphism. We show that among various life-history traits, mating systems have the greatest influence on patterns of polymorphism.     

The Nonrandom Location of Synonymous Codons Suggests That Reading Frame-Independent Forces Have Patterned Codon Preferences

Biased codon usage is common in eukaryotic and prokaryotic genes. Evidence from Escherichia, Saccharomyces, and Drosophila indicates that it favors translational efficiency and accuracy. However, to date no functional advantages have been identified in the codon–anticodon interactions involving the most frequently used (preferred) codons. Here we present evidence that forces not related to the individual codon–anticodon interaction may be involved in determining which synonymous codons are preferred or avoided. We show that the ``off-frame' trinucleotide motif preferences inferrable from Drosophila coding regions are often in the same direction as Drosophila's ``in-frame' codon preferences, i.e., its codon usage. The off-frame preferences were inferred from the nonrandomness of the location of confamilial synonymous codons along coding regions—a pattern often described as a context dependence of nucleotide choice at synonymous positions or as codon-pair bias. We relied on randomizations of the location of confamilial codons that do not alter, and cannot be influenced by, the encoded amino acid sequences, codon usage, or base composition of the genes examined. The statistically significant congruency of in-frame and off-frame trinucleotide preferences suggests that the same kind of reading-frame-independent force(s) may also influence synonymous codon choice. These forces may have produced biases in codon usage that then led to the evolution of the translational advantages of these motifs as preferred codons. Under this scenario, tRNA pool size differences between preferred and nonpreferred codons initially were evolved to track the default overrepresentation of codons with preferred motifs. The motif preference hypothesis can explain the structuring of codon preferences and the similarities in the codon usages of distantly related organisms. Received: 10 November 1998 / Accepted: 23 February 1999