Mutational pressure as the main cause of molecular evolution and polymorphism

The effects of mutation and recombination rate heterogeneity on the inference of demography and the distribution of fitness effects

Disentangling the effects of demography and selection has remained a focal point of population genetic analysis. Knowledge about mutation and recombination is essential in this endeavour; however, despite clear evidence that both mutation and recombination rates vary across genomes, it is common practice to model both rates as fixed. In this study, we quantify how this unaccounted for rate heterogeneity may impact inference using common approaches for inferring selection (DFE-alpha, Grapes, and polyDFE) and/or demography (fastsimcoal2 and δaδi). We demonstrate that, if not properly modelled, this heterogeneity can increase uncertainty in the estimation of demographic and selective parameters and in some scenarios may result in mis-leading inference. These results highlight the importance of quantifying the fundamental evolutionary parameters of mutation and recombination prior to utilizing population genomic data to quantify the effects of genetic drift (i.e., as modulated by demographic history) and selection; or, at the least, that the effects of uncertainty in these parameters can and should be directly modelled in downstream inference.

An efficient and robust ABC approach to infer the rate and strength of adaptation

Inferring the effects of positive selection on genomes remains a critical step in characterizing the ultimate and proximate causes of adaptation across species, and quantifying positive selection remains a challenge due to the confounding effects of many other evolutionary processes. Robust and efficient approaches for adaptation inference could help characterize the rate and strength of adaptation in non-model species for which demographic history, mutational processes, and recombination patterns are not currently well-described. Here, we introduce an efficient and user-friendly extension of the McDonald-Kreitman test (ABC-MK) for quantifying long-term protein adaptation in specific lineages of interest. We characterize the performance of our approach with forward simulations and find that it is robust to many demographic perturbations and positive selection configurations, demonstrating its suitability for applications to non-model genomes. We apply ABC-MK to the human proteome and a set of known Virus Interacting Proteins (VIPs) to test the long-term adaptation in genes interacting with viruses. We find substantially stronger signatures of positive selection on RNA-VIPs than DNA-VIPs, suggesting that RNA viruses may be an important driver of human adaptation over deep evolutionary time scales.

Evolutionarily stable strategy analysis and its links to demography and genetics through invasion fitness

Evolutionarily stable strategy (ESS) analysis pioneered by Maynard Smith and Price took off in part because it often does not require explicit assumptions about the genetics and demography of a population in contrast to population genetic models. Though this simplicity is useful, it obscures the degree to which ESS analysis applies to populations with more realistic genetics and demography: for example, how does ESS analysis handle complexities such as kin selection, group selection and variable environments when phenotypes are affected by multiple genes? In this paper, I review the history of the ESS concept and show how early uncertainty about the method lead to important mathematical theory linking ESS analysis to general population genetic models. I use this theory to emphasize the link between ESS analysis and the concept of invasion fitness. I give examples of how invasion fitness can measure kin selection, group selection and the evolution of linked modifier genes in response to variable environments. The ESSs in these examples depend crucially on demographic and genetic parameters, which highlights how ESS analysis will continue to be an important tool in understanding evolutionary patterns as new models address the increasing abundance of genetic and long-term demographic data in natural populations. This article is part of the theme issue 'Half a century of evolutionary games: a synthesis of theory, application and future directions'.

Chance, Variation and Shared Ancestry: Population Genetics After the Synthesis

Chance has been a focus of attention ever since the beginning of population genetics, but neutrality has not, as natural selection once appeared to be the only worthwhile issue. Neutral change became a major source of interest during the neutralist-selectionist debate, 1970-1980. It retained interest beyond this period for two reasons that contributed to its becoming foundational for evolutionary reasoning. On the one hand, neutral evolution was the first mathematical prediction to emerge from Mendelian inheritance: until then evolution by natural selection was considered the alternative to the fixity of species; now it appears to be the alternative to continuous change. Second, neutral change generated a set of clear predictions on standing variation. These could be used as a reference for detecting more elusive alternative mechanisms of evolution including natural selection. In the wake of the transition from Mendelism to genomics, the combination of coalescent theory, DNA sequence variation, and numerical analysis made it possible to integrate contingent aspects of the history of species into a new null model, thus opening a new dimension in the concept of population that the Modern Synthesis formerly considered as a mere gene pool.

Population Genomics of Paramecium Species

Population-genomic analyses are essential to understanding factors shaping genomic variation and lineage-specific sequence constraints. The dearth of such analyses for unicellular eukaryotes prompted us to assess genomic variation in Paramecium, one of the most well-studied ciliate genera. The Paramecium aurelia complex consists of ∼15 morphologically indistinguishable species that diverged subsequent to two rounds of whole-genome duplications (WGDs, as long as 320 MYA) and possess extremely streamlined genomes. We examine patterns of both nuclear and mitochondrial polymorphism, by sequencing whole genomes of 10-13 worldwide isolates of each of three species belonging to the P. aurelia complex: P. tetraurelia, P. biaurelia, P. sexaurelia, as well as two outgroup species that do not share the WGDs: P. caudatum and P. multimicronucleatum. An apparent absence of global geographic population structure suggests continuous or recent dispersal of Paramecium over long distances. Intergenic regions are highly constrained relative to coding sequences, especially in P. caudatum and P. multimicronucleatum that have shorter intergenic distances. Sequence diversity and divergence are reduced up to ∼100-150 bp both upstream and downstream of genes, suggesting strong constraints imposed by the presence of densely packed regulatory modules. In addition, comparison of sequence variation at non-synonymous and synonymous sites suggests similar recent selective pressures on paralogs within and orthologs across the deeply diverging species. This study presents the first genome-wide population-genomic analysis in ciliates and provides a valuable resource for future studies in evolutionary and functional genetics in Paramecium.

New thoughts on an old riddle: What determines genetic diversity within and between species?

The question of what determines genetic diversity has long remained unsolved by the modern evolutionary theory (MET). However, it has not deterred researchers from producing interpretations of genetic diversity by using MET. We examine the two observations of genetic diversity made in the 1960s that contributed to the development of MET. The interpretations of these observations by MET are widely known to be inadequate. We review the recent progress of an alternative framework, the maximum genetic diversity (MGD) hypothesis, that uses axioms and natural selection to explain the vast majority of genetic diversity as being at equilibrium that is largely determined by organismal complexity. The MGD hypothesis absorbs the proven virtues of MET and considers its assumptions relevant only to a much more limited scope. This new synthesis has accounted for the overlooked phenomenon of progression towards higher complexity, and more importantly, been instrumental in directing productive research.

The distribution of deleterious genetic variation in human populations

Population genetic studies suggest that most amino-acid changing mutations are deleterious. Such mutations are of tremendous interest in human population genetics as they are important for the evolutionary process and may contribute risk to common disease. Genomic studies over the past 5 years have documented differences across populations in the number of heterozygous deleterious genotypes, number of homozygous derived deleterious genotypes, number of deleterious segregating sites and proportion of sites that are potentially deleterious. These differences have been attributed to population history affecting the ability of natural selection to remove deleterious variants from the population. However, recent studies have suggested that the genetic load is the same across populations and that the efficacy of natural selection has not differed across human populations. Here I show that these observations are not incompatible with each other and that the apparent differences are due to examining different features of the genetic data and differing definitions of terms.

Molecular evolutionary history of Sugarcane yellow leaf virus based on sequence analysis of RNA-dependent RNA polymerase and putative aphid transmission factor-coding genes

RNA-dependent RNA polymerase (RdRp) encoded by ORF2 and putative aphid transmission factor (PATF) encoded by ORF5 of Sugarcane yellow leaf virus (SCYLV) were detected in six sugarcane cultivars affected by yellow leaf using RT-PCR and real-time RT-PCR assays. Expression of both genes varied among infected plants, but overall expression of RdRp was higher than expression of PATF. Cultivar H87-4094 from Hawaii yielded the highest transcript levels of RdRp, whereas cultivar C1051-73 from Cuba exhibited the lowest levels. Sequence comparisons among 25 SCYLV isolates from various geographical locations revealed an amino acid similarity of 72.1-99.4 and 84.7-99.8 % for the RdRp and PATF genes, respectively. The 25 SCYLV isolates were separated into three (RdRp) and two (PATF) phylogenetic groups using the MEGA6 program that does not account for genetic recombination. However, the SCYLV genome contained potential recombination signals in the RdRp and PATF coding genes based on the GARD genetic algorithm. Use of this later program resulted in the reconstruction of phylogenies on the left as well as on the right sides of the putative recombination breaking points, and the 25 SCYLV isolates were distributed into three distinct phylogenetic groups based on either RdRp or PATF sequences. As a result, recombination reshuffled the affiliation of the accessions to the different clusters. Analysis of selection pressures exerted on RdRp and PATF encoded proteins revealed that ORF 2 and ORF 5 underwent predominantly purifying selection. However, a few sites were also under positive selection as assessed by various models such as FEL, IFEL, REL, FUBAR, MEME, GA-Branch, and PRIME.

Interaction between selection and biased gene conversion in mammalian protein-coding sequence evolution revealed by a phylogenetic covariance analysis

According to the nearly-neutral model, variation in long-term effective population size among species should result in correlated variation in the ratio of nonsynonymous over synonymous substitution rates (dN/dS). Previous empirical investigations in mammals have been consistent with this prediction, suggesting an important role for nearly-neutral effects on protein-coding sequence evolution. GC-biased gene conversion (gBGC), on the other hand, is increasingly recognized as a major evolutionary force shaping genome nucleotide composition. When sufficiently strong compared with random drift, gBGC may significantly interfere with a nearly-neutral regime and impact dN/dS in a complex manner. Here, we investigate the phylogenetic correlations between dN/dS, the equilibrium GC composition (GC*), and several life-history and karyotypic traits in placental mammals. We show that the equilibrium GC composition decreases with body mass and increases with the number of chromosomes, suggesting a modulation of the strength of biased gene conversion due to changes in effective population size and genome-wide recombination rate. The variation in dN/dS is complex and only partially fits the prediction of the nearly-neutral theory. However, specifically restricting estimation of the dN/dS ratio on GC-conservative transversions, which are immune from gBGC, results in correlations that are more compatible with a nearly-neutral interpretation. Our investigation indicates the presence of complex interactions between selection and biased gene conversion and suggests that further mechanistic development is warranted, to tease out mutation, selection, drift, and conversion.

Weak selection and protein evolution

The "nearly neutral" theory of molecular evolution proposes that many features of genomes arise from the interaction of three weak evolutionary forces: mutation, genetic drift, and natural selection acting at its limit of efficacy. Such forces generally have little impact on allele frequencies within populations from generation to generation but can have substantial effects on long-term evolution. The evolutionary dynamics of weakly selected mutations are highly sensitive to population size, and near neutrality was initially proposed as an adjustment to the neutral theory to account for general patterns in available protein and DNA variation data. Here, we review the motivation for the nearly neutral theory, discuss the structure of the model and its predictions, and evaluate current empirical support for interactions among weak evolutionary forces in protein evolution. Near neutrality may be a prevalent mode of evolution across a range of functional categories of mutations and taxa. However, multiple evolutionary mechanisms (including adaptive evolution, linked selection, changes in fitness-effect distributions, and weak selection) can often explain the same patterns of genome variation. Strong parameter sensitivity remains a limitation of the nearly neutral model, and we discuss concave fitness functions as a plausible underlying basis for weak selection.

Comparative analysis of encephalization in mammals reveals relaxed constraints on anthropoid primate and cetacean brain scaling

There is a well-established allometric relationship between brain and body mass in mammals. Deviation of relatively increased brain size from this pattern appears to coincide with enhanced cognitive abilities. To examine whether there is a phylogenetic structure to such episodes of changes in encephalization across mammals, we used phylogenetic techniques to analyse brain mass, body mass and encephalization quotient (EQ) among 630 extant mammalian species. Among all mammals, anthropoid primates and odontocete cetaceans have significantly greater variance in EQ, suggesting that evolutionary constraints that result in a strict correlation between brain and body mass have independently become relaxed. Moreover, ancestral state reconstructions of absolute brain mass, body mass and EQ revealed patterns of increase and decrease in EQ within anthropoid primates and cetaceans. We propose both neutral drift and selective factors may have played a role in the evolution of brain-body allometry.

The neutral theory of molecular evolution in the genomic era

The neutral theory of molecular evolution has been widely accepted and is the guiding principle for studying evolutionary genomics and the molecular basis of phenotypic evolution. Recent data on genomic evolution are generally consistent with the neutral theory. However, many recently published papers claim the detection of positive Darwinian selection via the use of new statistical methods. Examination of these methods has shown that their theoretical bases are not well established and often result in high rates of false-positive and false-negative results. When the deficiencies of these statistical methods are rectified, the results become largely consistent with the neutral theory. At present, genome-wide analyses of natural selection consist of collections of single-locus analyses. However, because phenotypic evolution is controlled by the interaction of many genes, the study of natural selection ought to take such interactions into account. Experimental studies of evolution will also be crucial.

The early phase of a bacterial insertion sequence infection

Bacterial insertion sequences are the simplest form of autonomous mobile DNA. It is unknown whether they need to have beneficial effects to infect and persist in bacterial populations, or whether horizontal gene transfer suffices for their persistence. We address this question by using branching process models to investigate the critical, early phase of an insertion sequence infection. We find that the probability of a successful infection is low and depends linearly on the difference between the rate of horizontal gene transfer and the fitness cost of the insertion sequences. Our models show that the median time to extinction of an insertion sequence that dies out is very short, while the median time for a successful infection to reach a modest population size is very long. We conclude that horizontal gene transfer is strong enough to allow the persistence of insertion sequences, although infection is an erratic and slow process.

Shifts in the intensity of purifying selection: an analysis of genome-wide polymorphism data from two closely related yeast species

How much does the intensity of purifying selection vary among populations and species? How uniform are the shifts in selective pressures across the genome? To address these questions, we took advantage of a recent, whole-genome polymorphism data set from two closely related species of yeast, Saccharomyces cerevisiae and S. paradoxus, paying close attention to the population structure within these species. We found that the average intensity of purifying selection on amino acid sites varies markedly among populations and between species. As expected in the presence of extensive weakly deleterious mutations, the effect of purifying selection is substantially weaker on single nucleotide polymorphisms (SNPs) segregating within populations than on SNPs fixed between population samples. Also in accordance with a Nearly Neutral model, the variation in the intensity of purifying selection across populations corresponds almost perfectly to simple measures of their effective size. As a first step toward understanding the processes generating these patterns, we sought to tease apart the relative importance of systematic, genome-wide changes in the efficacy of selection, such as those expected from demographic processes and of gene-specific changes, which may be expected after a shift in selective pressures. For that purpose, we developed a new model for the evolution of purifying selection between populations and inferred its parameters from the genome-wide data using a likelihood approach. We found that most, but not all changes seem to be explained by systematic shifts in the efficacy of selection. One population, the sake-derived strains of S. cerevisiae, however, also shows extensive gene-specific changes, plausibly associated with domestication. These findings have important implications for our understanding of purifying selection as well as for estimates of the rate of molecular adaptation in yeast and in other species.

Model of effectively neutral mutations in which selective constraint is incorporated

Based on the idea that selective neutrality is the limit when the selective disadvantage becomes indefinitely small, a model of neutral (and nearly neutral) mutations is proposed that assumes that the selection coefficient (s') against the mutant at various sites within a cistron (gene) follows a Gamma distribution; f(s') = alpha(beta)e(-alphas')s'(beta-1)/Gamma(beta), in which alpha = beta/ s' and s' is the mean selection coefficient against the mutants ( s' > 0; 1 [unk] beta > 0). The mutation rate for alleles whose selection coefficients s' lie in the range between 0 and 1/(2N(e)), in which N(e) is the effective population size, is termed the effectively neutral mutation rate (denoted by v(e)). Using the model of "infinite sites" in population genetics, formulas are derived giving the average heterozygosity ( h(e)) and evolutionary rate per generation (k(g)) in terms of mutant substitutions. It is shown that, with parameter values such as beta = 0.5 and s' = 0.001, the average heterozygosity increases much more slowly as N(e) increases, compared with the case in which a constant fraction of mutations are neutral. Furthermore, the rate of evolution per year (k(1)) becomes constant among various organisms, if the generation span (g) in years is inversely proportional to radicalN(e) among them and if the mutation rate per generation is constant. Also, it is shown that we have roughly k(g) = v(e). The situation becomes quite different if slightly advantageous mutations occur at a constant rate independent of environmental conditions. In this case, the evolutionary rate can become enormously higher in a species with a very large population size than in a species with a small population size, contrary to the observed pattern of evolution at the molecular level.

Protein variation in the marine bivalve Scrobicularia plana

Seventeen enzyme loci have been assayed for electrophoretically detectable variation in a population of the marine bivalve Scrobicularia plana. Mean heterozygosity is 0.120 +/- 0.033. In a comparison involving thirteen enzymes there is a significant correlation between heterozygosity in S. plana and Mytilus edulis and a suggestion of lower mean heterozygosity in S. plana. These findings are discussed in relation to current theories concerning the selective significance of protein variation.

Myoglobin polymorphism in the pollan (Osteichthyes: Coregoninae)

Three electrophoretically separable phenotypes of heart and lateral line muscle myoglobin were found in the Irish pollan (Coregonus pollan). This polymorphism appears to be under the control of two co-dominant alleles. The allele frequencies were found to be virtually identical in samples from two lakes which have probably been isolated since the close of the last Ice Age. A significant excess of heterozygotes was found in samples from both lakes. This myoglobin polymorphism appears to be balanced, maintained due to heterozygote superiority.

The frailty of adaptive hypotheses for the origins of organismal complexity

The vast majority of biologists engaged in evolutionary studies interpret virtually every aspect of biodiversity in adaptive terms. This narrow view of evolution has become untenable in light of recent observations from genomic sequencing and population-genetic theory. Numerous aspects of genomic architecture, gene structure, and developmental pathways are difficult to explain without invoking the nonadaptive forces of genetic drift and mutation. In addition, emergent biological features such as complexity, modularity, and evolvability, all of which are current targets of considerable speculation, may be nothing more than indirect by-products of processes operating at lower levels of organization. These issues are examined in the context of the view that the origins of many aspects of biological diversity, from gene-structural embellishments to novelties at the phenotypic level, have roots in nonadaptive processes, with the population-genetic environment imposing strong directionality on the paths that are open to evolutionary exploitation.

Selectionism and neutralism in molecular evolution

Charles Darwin proposed that evolution occurs primarily by natural selection, but this view has been controversial from the beginning. Two of the major opposing views have been mutationism and neutralism. Early molecular studies suggested that most amino acid substitutions in proteins are neutral or nearly neutral and the functional change of proteins occurs by a few key amino acid substitutions. This suggestion generated an intense controversy over selectionism and neutralism. This controversy is partially caused by Kimura's definition of neutrality, which was too strict (|2Ns|< or =1). If we define neutral mutations as the mutations that do not change the function of gene products appreciably, many controversies disappear because slightly deleterious and slightly advantageous mutations are engulfed by neutral mutations. The ratio of the rate of nonsynonymous nucleotide substitution to that of synonymous substitution is a useful quantity to study positive Darwinian selection operating at highly variable genetic loci, but it does not necessarily detect adaptively important codons. Previously, multigene families were thought to evolve following the model of concerted evolution, but new evidence indicates that most of them evolve by a birth-and-death process of duplicate genes. It is now clear that most phenotypic characters or genetic systems such as the adaptive immune system in vertebrates are controlled by the interaction of a number of multigene families, which are often evolutionarily related and are subject to birth-and-death evolution. Therefore, it is important to study the mechanisms of gene family interaction for understanding phenotypic evolution. Because gene duplication occurs more or less at random, phenotypic evolution contains some fortuitous elements, though the environmental factors also play an important role. The randomness of phenotypic evolution is qualitatively different from allele frequency changes by random genetic drift. However, there is some similarity between phenotypic and molecular evolution with respect to functional or environmental constraints and evolutionary rate. It appears that mutation (including gene duplication and other DNA changes) is the driving force of evolution at both the genic and the phenotypic levels.

The origins of eukaryotic gene structure

Most of the phenotypic diversity that we perceive in the natural world is directly attributable to the peculiar structure of the eukaryotic gene, which harbors numerous embellishments relative to the situation in prokaryotes. The most profound changes include introns that must be spliced out of precursor mRNAs, transcribed but untranslated leader and trailer sequences (untranslated regions), modular regulatory elements that drive patterns of gene expression, and expansive intergenic regions that harbor additional diffuse control mechanisms. Explaining the origins of these features is difficult because they each impose an intrinsic disadvantage by increasing the genic mutation rate to defective alleles. To address these issues, a general hypothesis for the emergence of eukaryotic gene structure is provided here. Extensive information on absolute population sizes, recombination rates, and mutation rates strongly supports the view that eukaryotes have reduced genetic effective population sizes relative to prokaryotes, with especially extreme reductions being the rule in multicellular lineages. The resultant increase in the power of random genetic drift appears to be sufficient to overwhelm the weak mutational disadvantages associated with most novel aspects of the eukaryotic gene, supporting the idea that most such changes are simple outcomes of semi-neutral processes rather than direct products of natural selection. However, by establishing an essentially permanent change in the population-genetic environment permissive to the genome-wide repatterning of gene structure, the eukaryotic condition also promoted a reliable resource from which natural selection could secondarily build novel forms of organismal complexity. Under this hypothesis, arguments based on molecular, cellular, and/or physiological constraints are insufficient to explain the disparities in gene, genomic, and phenotypic complexity between prokaryotes and eukaryotes.

Phylogenetic position of Phthiraptera (Insecta: Paraneoptera) and elevated rate of evolution in mitochondrial 12S and 16S rDNA

Phthiraptera (chewing and sucking lice) and Psocoptera (booklice and barklice) are closely related to each other and compose the monophyletic taxon Psocodea. However, there are two hypotheses regarding their phylogenetic relationship: (1) monophyletic Psocoptera is the sister group of Phthiraptera or (2) Psocoptera is paraphyletic, and Liposcelididae of Psocoptera is the sister group of Phthiraptera. Each hypothesis is supported morphologically and/or embryologically, and this problem has not yet been resolved. In the present study, the phylogenetic position of Phthiraptera was examined using mitochondrial 12S and 16S rDNA sequences, with three methods of phylogenetic analysis. Results of all analyses strongly supported the close relationship between Phthiraptera and Liposcelididae. Results of the present analyses also provided some insight into the elevated rate of evolution in mitochondrial DNA (mtDNA) in Phthiraptera. An elevated substitution rate of mtDNA appears to originate in the common ancestor of Phthiraptera and Liposcelididae, and directly corresponds to an increased G+C content. Therefore, the elevated substitution rate of mtDNA in Phthiraptera and Liposcelididae appears to be directional. A high diversity of 12S rDNA secondary structure was also observed in wide range of Phthiraptera and Liposcelididae, but these structures seem to have evolved independently in different clades.

Sequence variation of alcohol dehydrogenase (Adh) paralogs in cactophilic Drosophila

This study focuses on the population genetics of alcohol dehydrogenase (Adh) in cactophilic Drosophila. Drosophila mojavensis and D. arizonae utilize cactus hosts, and each host contains a characteristic mixture of alcohol compounds. In these Drosophila species there are two functional Adh loci, an adult form (Adh-2) and a larval and ovarian form (Adh-1). Overall, the greater level of variation segregating in D. arizonae than in D. mojavensis suggests a larger population size for D. arizonae. There are markedly different patterns of variation between the paralogs across both species. A 16-bp intron haplotype segregates in both species at Adh-2, apparently the product of an ancient gene conversion event between the paralogs, which suggests that there is selection for the maintenance of the intron structure possibly for the maintenance of pre-mRNA structure. We observe a pattern of variation consistent with adaptive protein evolution in the D. mojavensis lineage at Adh-1, suggesting that the cactus host shift that occurred in the divergence of D. mojavensis from D. arizonae had an effect on the evolution of the larval expressed paralog. Contrary to previous work we estimate a recent time for both the divergence of D. mojavensis and D. arizonae (2.4 +/- 0.7 MY) and the age of the gene duplication (3.95 +/- 0.45 MY).

Acceleration of genomic evolution caused by enhanced mutation rate in endocellular symbionts

Endosymbionts, which are widely observed in nature, have undergone reductive genome evolution because of their long-term intracellular lifestyle. Here we compared the complete genome sequences of two different endosymbionts, Buchnera and a protist mitochondrion, with their close relatives to study the evolutionary rates of functional genes in endosymbionts. The results indicate that the rate of amino acid substitution is two times higher in symbionts than in their relatives. This rate increase was observed uniformly among different functional classes of genes, although strong purifying selection may have counterbalanced the rate increase in a few cases. Our data suggest that, contrary to current views, neither the Muller's ratchet effect nor the slightly deleterious mutation theory sufficiently accounts for the elevated evolutionary rate. Rather, the elevated evolutionary rate appears to be mainly due to enhanced mutation rate, although the possibility of relaxation of purifying selection cannot be ruled out.

Neutralism and selectionism: the molecular clock

The neutrality theory predicts that the rate of molecular evolution will be constant over time, and thus that there is a molecular clock for timing evolutionary events. It has been observed that the variance of the rate of evolution is generally larger than expected according to the neutrality theory. Several modifications of the theory have been proposed to account for the 'overdispersion' of the molecular clock, by postulating effects attributed to generation-time, population size, slightly deleterious mutations, repair mechanisms, and the like. An extensive investigation of two proteins, glycerol-3-phosphate dehydrogenase (GPDH) and superoxide dismutase (SOD), manifests that none of these modifications can simultaneously account for the disparate patterns observed in both proteins. GPDH evolves very slowly in Drosophila species, but several times faster in mammals, other animals, plants, and fungi. SOD evolves very fast in Drosophila species and also in mammals, but much more slowly in other animals and still slower when plants and fungi are compared to one another, or to animals.

Extraordinary multilocus genetic organization in mole crickets, Gryllotalpidae

Allozymic diversity at 21 loci was analyzed in 470 individuals of three species of mole cricket superspecies, Gryllotalpa gryllotalpa (two new chromosomal species, G. tali and G. marismortui) and G. africana in Israel, which are distributed along a southward transect of increasing aridity. Two outstanding findings emerged in G. tali and G. marismortui: (1) genetic polymorphism was high but heterozygosity very low, indicating significant deviations from Hardy-Weinberg expectations; and (2) significant linkage disequilibria at an unprecedented level for outbreeders and remarkable intersite differences. The results may characterize subterranean gryllotalpids worldwide because a single sample of Neocurtilla hexadactyla from Tefé, Amazonia, shows the same features. Significant variation of heterozygote paucity among loci, combined with the biology of the species, rejects the simple explanation of inbreeding or any other single explanatory model. Likewise, direct selection against heterozygotes or specific multilocus associations can explain, but is not necessary nor likely to explain, the observed results in mole crickets. To explain these results, we developed a multiple-factor mathematical model combining niche viability selection, niche choice, and positive assortative mating. This model involves a special case of Wahlund effect and inbreeding. Simulations based on this model showed that a combination of these three mechanisms may produce the observed distribution of alleles, via selection on a few loci, to affect the entire genome organization.

MtDNA phylogeography and recent intra-island diversification among Canary Island Calathus beetles

The genus Calathus Bonelli comprises 24 species on the Canary Islands. Sequences of 927 and 687 bp of the mitochondrial cytochrome oxidase I and II genes, respectively, as well as the intervening tRNA leu gene in 21 of the 24 species, have identified three genetically divergent and unequivocally monophyletic groupings. A phylogeographic analysis is presented for the major monophyletic group comprising all the species of Gran Canaria, La Gomera, and El Hierro, and two Tenerifean species. A distance-based phylogenetic analysis and maximum parsimony analysis have clearly shown that this clade is composed of four distinct lineages. DNA sequence data suggest a recent origin for this clade and that lineages have not evolved at the same rate. Compared with diversification patterns observed in other Coleoptera on the Canary Islands, diversification has been recent relative to the time of colonization within the islands of Gran Canaria and La Gomera. Calathus diversification on La Gomera has been greater than on Gran Canaria. The influences of geological and ecological history are discussed in relation to Calathus diversification.

Deleterious mutations destabilize ribosomal RNA in endosymbiotic bacteria

In populations that are small and asexual, mutations with slight negative effects on fitness will drift to fixation more often than in large or sexual populations in which they will be eliminated by selection. If such mutations occur in substantial numbers, the combined effects of long-term asexuality and small population size may result in substantial accumulation of mildly deleterious substitutions. Prokaryotic endosymbionts of animals that are transmitted maternally for very long periods are effectively asexual and experience smaller effective population size than their free-living relatives. The contrast between such endosymbionts and related free-living bacteria allows us to test whether a population structure imposing frequent bottlenecks and asexuality does lead to an accumulation of slightly deleterious substitutions. Here we show that several independently derived insect endosymbionts, each with a long history of maternal transmission, have accumulated destabilizing base substitutions in the highly conserved 16S rRNA. Stabilities of Domain I of this subunit are 15-25% lower in endosymbionts than in closely related free-living bacteria. By mapping destabilizing substitutions onto a reconstructed phylogeny, we show that decreased ribosomal stability has evolved separately in each endosymbiont lineage. Our phylogenetic approach allows us to demonstrate statistical significance for this pattern: becoming endosymbiotic predictably results in decreased stability of rRNA secondary structure.

Conditions for positive and negative correlations between fitness and heterozygosity in equilibrium populations

The past decades have witnessed extensive efforts to correlate fitness traits with genomic heterozygosity. While positive correlations are revealed in most of the organisms studied, results of no/negative correlations are not uncommon. There has been little effort to reveal the genetic causes of these negative correlations. The positive correlations are regarded either as evidence for functional overdominance in large, randomly mating populations at equilibrium, or the results of populations at disequilibrium under dominance. More often, the positive correlations are viewed as a phenomenon of heterosis, so that it cannot possibly occur under within-locus additive allelic effects. Here we give exact genetic conditions that give rise to positive and negative correlations in populations at Hardy-Weinberg and linkage equilibria, thus offering a genetic explanation for the observed negative correlations. Our results demonstrate that the above interpretations concerning the positive correlations are not complete or even necessary. Such a positive correlation can result under dominance and potentially under additivity, even in populations where associated overdominance due to linked alleles at different loci is not significant. Additionally, negative correlations and heterosis can co-occur in a single population. Although our emphasis is on equilibrium populations and for biallelic genetic systems, the basic conclusions are generalized to non-equilibrium populations and for multi-allelic situations.

Marginal fitness contributions of nonessential genes in yeast

Analysis of the complete genome sequence of Saccharomyces cerevisiae confirms and extends earlier evidence that a majority of yeast genes are not essential, at least under laboratory conditions. Many fail to yield a discernible mutant phenotype even when disrupted. Genes not subject to natural selection would accumulate inactivating mutations, so these "cryptic" genes must have functions that are overlooked by the standard methods of yeast genetics. Two explanations seem possible: (i) They have important functions only in environments not yet duplicated in the laboratory and would have conditional phenotypes if tested appropriately. (ii) They make small, but significant, contributions to fitness even under routine growth conditions, but the effects are not large enough to be detected by conventional methods. We have tested the second "marginal benefit" hypothesis by measuring the fitnesses of a random collection of disruption mutants in direct competition with their wild-type progenitor. A substantial majority of mutant strains that lack obvious defects nevertheless are at a significant selective disadvantage just growing on rich medium under normal conditions. This result has important implications for efforts to understand the functions of novel genes revealed by sequencing projects.

Synonymous and nonsynonymous substitutions in mammalian genes and the nearly neutral theory

The nearly neutral theory of molecular evolution predicts larger generation-time effects for synonymous than for nonsynonymous substitutions. This prediction is tested using the sequences of 49 single-copy genes by calculating the average and variance of synonymous and nonsynonymous substitutions in mammalian star phylogenies (rodentia, artiodactyla, and primates). The average pattern of the 49 genes supports the prediction of the nearly neutral theory, with some notable exceptions. The nearly neutral theory also predicts that the variance of the evolutionary rate is larger than the value predicted by the completely neutral theory. This prediction is tested by examining the dispersion index (ratio of the variance to the mean), which is positively correlated with the average substitution number. After weighting by the lineage effects, this correlation almost disappears for nonsynonymous substitutions, but not quite so for synonymous substitutions. After weighting, the dispersion indices of both synonymous and nonsynonymous substitutions still exceed values expected under the simple Poisson process. The results indicate that both the systematic bias in evolutionary rate among the lineages and the episodic type of rate variation are contributing to the large variance. The former is more significant to synonymous substitutions than to nonsynonymous substitutions. Isochore evolution may be similar to synonymous substitutions. The rate and pattern found here are consistent with the nearly neutral theory, such that the relative contributions of drift and selection differ between the two types of substitutions. The results are also consistent with Gillespie's episodic selection theory.

The effect of background selection against deleterious mutations on weakly selected, linked variants

This paper analyses the effects of selection against deleterious alleles maintained by mutation ('background selection') on rates of evolution and levels of genetic diversity at weakly selected, completely linked, loci. General formulae are derived for the expected rates of gene substitution and genetic diversity, relative to the neutral case, as a function of selection and dominance coefficients at the loci in question, and of the frequency of gametes that are free of deleterious mutations with respect to the loci responsible for background selection. As in the neutral case, most effects of background selection can be predicted by considering the effective size of the population to be multiplied by the frequency of mutation-free gametes. Levels of genetic diversity can be sharply reduced by background selection, with the result that values for sites under selection approach those for neutral variants subject to the same regime of background selection. Rates of fixation of slightly deleterious mutations are increased by background selection, and rates of fixation of advantageous mutations are reduced. The properties of sex-linked and autosomal loci in random-mating populations are compared, and the effects of background selection on asexual and self-fertilizing populations are considered. The implications of these results for the interpretation of studies of molecular evolution and variation are discussed.