Role of very slightly deleterious mutations in molecular evolution and polymorphism

Several models of multiple slightly deleterious alleles are reviewed and theoretical consequences of slightly negative selection are discussed in conjunction with evolution and variation at the molecular level. Facts are organized which may be satisfactorily explained by the hypothesis of very slightly deleterious mutations. They are: (1) There appears to be an upper limit in heterozygosity for protein loci as measured by electrophoresis. (2) The excess of rare alleles is more pronounced in Drosophila than in man. (3) Correlation of heterozygosities at a locus among sibling species of the Drosophila willistoni group is too high compared to what is expected by the strict neutral theory, while it is not so among human races and between man and chimpanzee. (4) The rate of protein divergence is exceptionally high in Hawaiian Drosophila.     

Very slightly deleterious mutations and the molecular clock

Summary The model of very slightly deleterious mutations was examined from the standpoint of population genetics in relation to the molecular evolutionary clock. The distribution of selection coefficients of mutants (in terms of amino acid changes) with small effect is thought to be continuous around zero, with an average negative value. The variance of selection coefficients depends upon environmental diversity and hence on total population size of a species. By considering various examples of amino acid substitutions, the average and standard error of selection coefficients and the reciprocal of population size are assumed to have similar values. The model predicts negative correlation between evolutionary rate and population size. This effect is expected to be partially cancelled with the generation time effect of intrinsic mutation rate. Implications of this prediction on the molecular evolutionary clock are discussed.     

The McDonald-Kreitman test and slightly deleterious mutations

It is possible to estimate the proportion of substitutions that are due to adaptive evolution using the numbers of silent and nonsilent polymorphisms and substitutions in a McDonald and Kreitman-type analysis. Unfortunately, this estimate of adaptive evolution is biased downward by the segregation of slightly deleterious mutations. It has been suggested that 1 way to cope with the effects of these slightly deleterious mutations is to remove low-frequency polymorphisms from the analysis. We investigate the performance of this method theoretically. We show that although removing low-frequency polymorphisms does indeed reduce the bias in the estimate of adaptive evolution, the estimate is always downwardly biased, often to the extent that one would not be able to detect adaptive evolution, even if it existed. The method is reasonably satisfactory, only if the rate of adaptive evolution is high and the distribution of fitness effects for slightly deleterious mutations is very leptokurtic. Our analysis suggests that adaptive evolution could be quite prevalent in humans (>8%) and still not be detectable using current methodologies. Our analysis also suggests that the level of adaptive evolution has probably been underestimated, possibly substantially, in both bacteria and Drosophila.     

Quantifying the impact of protein tertiary structure on molecular evolution

To investigate the evolutionary impact of protein structure, the experimentally determined tertiary structure and the protein-coding DNA sequence were collected for each of 1,195 genes. These genes were studied via a model of sequence change that explicitly incorporates effects on evolutionary rates due to protein tertiary structure. In the model, these effects act via the solvent accessibility environments and pairwise amino acid interactions that are induced by tertiary structure. To compare the hypotheses that structure does and does not have a strong influence on evolution, Bayes factors were estimated for each of the 1,195 sequences. Most of the Bayes factors strongly support the hypothesis that protein structure affects protein evolution. Furthermore, both solvent accessibility and pairwise interactions among amino acids are inferred to have important roles in protein evolution. Our results also indicate that the strength of the relationship between tertiary structure and evolution has a weak but real correlation to the annotation information in the Gene Ontology database. Although their influences on rates of evolution vary among protein families, we find that the mean impacts of solvent accessibility and pairwise interactions are about the same.     

Alleviation of deleterious effects of protein mutation through inactivation of molecular chaperones

Molecular chaperones recognize and bind destabilized proteins. This can be especially important for proteins whose stability is reduced by mutations. We focused our study on a major chaperone system, RAC-Ssb, which assists folding of newly synthesized polypeptides in the yeast cytosol. A sensitive phenotypic assay, the red color of Ade2 mutants, was used to screen for variants with metabolic activity dependent on RAC-Ssb. None of the Ade2 mutants were found to exhibit lower metabolic activity after inactivation of RAC-Ssb. In order to explicitly test the relationship between protein instability and activity of chaperones, a series of temperature sensitive Ade2 mutants were tested in the presence or absence of RAC-Ssb. The growth of Ade2(ts) mutants at elevated temperatures was enhanced if chaperones were missing. Similar pattern was found for thermally sensitive mutants of several other genes. Because RAC-Ssb normally supports the folding of proteins, it appears paradoxical that catabolic activity of mutants is reduced when these chaperones are present. We suggest that under non-stressful conditions, molecular chaperones are tuned to support folding of native proteins, but not that of mutated ones.     

Evidence for abundant slightly deleterious polymorphisms in bacterial populations

The nearly neutral theory of molecular evolution predicts that slightly deleterious mutations subject to purifying selection are widespread in natural populations, particularly those of large effective population size. To test this hypothesis, the standardized difference between pairwise nucleotide difference and number of segregation sites (corrected for number of sequences) was estimated for 149 population data sets from 84 species of bacteria. This quantity (Tajima's D-statistic) was estimated separately for synonymous (D(syn)) and nonsynonymous (D(non)) polymorphisms. D(syn) was positive in 70% of data sets, and the overall median D(syn) (0.873) was positive. By contrast D(non) was negative in 68% of data sets, and the overall median D(non) (-0.656) was negative. The preponderance of negative values of D(non) is evidence that there are widespread rare nonsynonymous polymorphisms in the process of being eliminated by purifying selection, as predicted to occur in populations with large effective size by the nearly neutral theory. The major exceptions to this trend were seen among surface proteins, particularly those of bacteria parasitic on vertebrates, which included a number of cases of polymorphisms apparently maintained by balancing selection.     

Phylogenetic evidence for deleterious mutation load in RNA viruses and its contribution to viral evolution

Populations of RNA viruses are often characterized by abundant genetic variation. However, the relative fitness of these mutations is largely unknown, although this information is central to our understanding of viral emergence, immune evasion, and drug resistance. Here we develop a phylogenetic method, based on the distribution of nonsynonymous and synonymous changes, to assess the relative fitness of polymorphisms in the structural genes of 143 RNA viruses. This reveals that a substantial proportion of the amino acid variation observed in natural populations of RNA viruses comprises transient deleterious mutations that are later purged by purifying selection, potentially limiting virus adaptability. We also demonstrate, for the first time, the existence of a relationship between amino acid variability and the phylogenetic distribution of polymorphisms. From this relationship, we propose an empirical threshold for the maximum viable deleterious mutation load in RNA viruses.     

The influence of hitchhiking and deleterious mutation upon asexual mutation rates

The question of how natural selection affects asexual mutation rates has been considered since the 1930s, yet our understanding continues to deepen. The distribution of mutation rates observed in natural bacteria remains unexplained. It is well known that environmental constancy can favor minimal mutation rates. In contrast, environmental fluctuation (e.g., at period T) can create indirect selective pressure for stronger mutators: genes modifying mutation rate may "hitchhike" to greater frequency along with environmentally favored mutations they produce. This article extends a well-known model of Leigh to consider fitness genes with multiple mutable sites (call the number of such sites alpha). The phenotypic effect of such a gene is enabled if all sites are in a certain state and disabled otherwise. The effects of multiple deleterious loci are also included (call the number of such loci gamma). The analysis calculates the indirect selective effects experienced by a gene inducing various mutation rates for given values of alpha, gamma, and T. Finite-population simulations validate these results and let us examine the interaction of drift with hitchhiking selection. We close by commenting on the importance of other factors, such as spatiotemporal variation, and on the origin of variation in mutation rates.     

Accumulation of slightly deleterious mutations in the mitochondrial genome: A hallmark of animal domestication

The hypothesis that domestication leads to a relaxation of purifying selection on mitochondrial (mt) genomes was tested by comparative analysis of mt genes from dog, pig, chicken, and silkworm. The three vertebrate species showed mt genome phylogenies in which domestic and wild isolates were intermingled, whereas the domestic silkworm (Bombyx mori) formed a distinct cluster nested within its closest wild relative (Bombyx mandarina). In spite of these differences in phylogenetic pattern, significantly greater proportions of nonsynonymous SNPs than of synonymous SNPs were unique to the domestic populations of all four species. Likewise, in all four species, significantly greater proportions of RNA-encoding SNPs than of synonymous SNPs were unique to the domestic populations. Thus, domestic populations were characterized by an excess of unique polymorphisms in two categories generally subject to purifying selection: nonsynonymous sites and RNA-encoding sites. Many of these unique polymorphisms thus seem likely to be slightly deleterious; the latter hypothesis was supported by the generally lower gene diversities of polymorphisms unique to domestic populations in comparison to those of polymorphisms shared by domestic and wild populations.     

The temporal dynamics of slightly deleterious mutations in Escherichia coli and Shigella spp

The Shigella are recently emerged clones of Escherichia coli, which have independently adopted an intracellular pathogenic lifestyle. We examined the molecular evolutionary consequences of this niche specialization by comparing the normalized, directional frequency profiles of unique polymorphisms within 2,098 orthologues representing the intersection of five E. coli and four Shigella genomes. We note a surfeit of AT-enriching changes (GC-->AT), transversions, and nonsynonymous changes in the Shigella genomes. By examining these differences within a temporal framework, we conclude that our results are consistent with relaxed or inefficient selection in Shigella owing to a reduced effective population size. Alternative interpretations, and the interesting exception of Shigella sonnei, are discussed. Finally, this analysis lends support to the view that nucleotide composition typically does not lie at mutational equilibrium but that selection plays a role in maintaining a higher GC content than would result solely from mutation bias. This argument sheds light on the enrichment of adenine and thymine in the genomes of bacterial endosymbionts where purifying selection is very weak.     

Epistasis between deleterious mutations and the evolution of recombination

Epistasis and the evolution of recombination are closely intertwined: epistasis generates linkage disequilibria (i.e. statistical associations between alleles), whereas recombination breaks them up. The mutational deterministic hypothesis (MDH) states that high recombination rates are maintained because the breaking up of linkage disequilibria generated by negative epistasis enables more efficient purging of deleterious mutations. However, recent theoretical and experimental work challenges the MDH. Experimental evidence suggests that negative epistasis, required by the MDH, is relatively uncommon. On the theoretical side, population genetic models suggest that, compared with the combined effects of drift and selection, epistasis generates a negligible amount of linkage disequilibria. Here, we assess these criticisms and discuss to what extent they invalidate the MDH as an explanation for the evolution of recombination.     

Two optimal mutation rates in obligate pathogens subject to deleterious mutation

Pathogen species with high mutation rates are likely to accumulate deleterious mutations that reduce their reproductive potential within the host. By altering the within-host growth rate of the pathogen, the deleterious mutation load has the potential to affect epidemiological properties such as prevalence, mean pathogen load, and the mean duration of infections. Here, I examine an epidemiological model that allows for multiple segregating mutations that affect within-host replication efficiency. The model demonstrates a complex range of outcomes depending on pathogen mutation rate, including two distinct, widely separated mutation rates associated with high pathogen prevalence. The low mutation rate prevalence peak is associated with small amounts of genetic diversity within the pathogen population, relatively stable prevalence and infection dynamics, and genetic variation partitioned between hosts. The high mutation rate peak is characterized by considerable genetic diversity both within and between hosts, relatively frequent invasions by more virulent types, and is qualitatively similar to an RNA virus quasispecies. The two prevalence peaks are separated by a valley where natural selection favors evolution toward the optimal within-host state, which is associated with high virulence and relatively rapid host mortality. Both chronic and acute infections are examined using stochastic forward simulations.     

About a three states molecular evolution model

The Eigen model of molecular evolution has recently been reconsidered by taking into account the existence of non-reproducing sequences. New results obtained with this model are presented here: several time-scales characterize the population dynamics and the location of the error threshold, shifted towards larger values of the probability of mutation, depends on both the size of the population and time. Thus this model can be used to analyse extinction events. The distribution of individuals in sequence space beyond the threshold is also discussed.     

Genome-wide deleterious mutation favors dispersal and species integrity

Here I develop the idea that ubiquitous harmful genome-wide mutation with local differentiation favors dispersal, even though migration reduces average fitness. Historical contingency of the mutational process means that demes (sub-populations) differentiate from one another. Deleterious or lethal partially recessive mutations carried by migrants then do not encounter similar mutations in the recipient deme. Migrant offspring have higher fitness than offspring of residents, because migrant offspring are heterozygous rather than homozygous for harmful mutations. The advantage is inversely related to local inbreeding depression. Genome-wide deleterious mutation favors the evolution of dispersal, which in turn enhances the genetic integrity of the species.     

Quantifying the phylodynamic forces driving papillomavirus evolution

The associations between pathogens and their hosts are complex and can result from a variety of evolutionary processes including codivergence, lateral transfer, or duplication. Papillomaviruses (PVs) are double-stranded DNA viruses ubiquitously present in mammals and are a suitable target for rigorous statistical tests of potential virus-host codivergence. We analyze the evolutionary dynamics of PV diversification by comparing robust phylogenies of PVs and their respective hosts using different statistical approaches to assess topological and branch-length congruence. Mammalian PVs segregated into four diverse major clades that overlapped to varying degrees in terms of their mammalian host lineages. The hypothesis that PVs and hosts evolved independently was globally rejected (P = 0.0001), although only 90 of 207 virus-host associations (43%) were significant in individual tests. Virus-host codivergence accounted roughly for one-third of the evolutionary events required to reconcile PV-host evolutionary histories. When virus-host associations were analyzed locally within each of the four viral clades, numerous independent topological congruencies were identified that were incompatible with respect to the global trees. These results support an evolutionary scenario in which early PV radiation was followed by independent codivergence between viruses within each of the major clades and their hosts. Moreover, heterogeneous groups of closely related PVs infecting non-related hosts suggest several interspecies transmission events. Our results argue thus for the importance of alternative events in PV evolution, in contrast to the prevailing opinion that these viruses show a high degree of host specificity and codivergence.     

Local recombination and mutation effects on molecular evolution in Drosophila.

I studied the cause of the significant difference in the synonymous-substitution pattern found in the achaete-scute complex genes in two Drosophila lineages, higher codon bias in Drosophila yakuba, and lower bias in D. melanogaster. Besides these genes, the functionally unrelated yellow gene showed the same substitution pattern, suggesting a region-dependent phenomenon in the X-chromosome telomere. Because the numbers of A/T --> G/C substitutions were not significantly different from those of G/C --> A/T in the yellow noncoding regions of these species, a AT/GC mutational bias could not completely account for the synonymous-substitution biases. In contrast, we did find an approximately 14-fold difference in recombination rates in the X-chromosome telomere regions between the two species, suggesting that the reduction of recombination rates in this region resulted in the reduction of the efficacy of selection in D. melanogaster. In addition, the D. orena yellow showed a 5% increase in the G + C content at silent sites in the coding and noncoding regions since the divergence from D. erecta. This pattern was significantly different from those at the orena Adh and Amy loci. These results suggest that local changes in recombination rates and mutational pressures are contributing to the irregular synonymous-substitution patterns in Drosophila.     

Directional mutation pressure,mutator mutations,and dynamics of molecular evolution

Using a general form of the directional mutation theory, this paper analyzes the effect of mutations in mutator genes on the G + C content of DNA, the frequency of substitution mutations, and evolutionary changes (cumulative mutations) under various degrees of selective constraints. Directional mutation theory predicts that when the mutational bias between A/T and G/C nucleotide pairs is equilibrated with the base composition of a neutral set of DNA nucleotides, the mutation frequency per gene will be much lower than the frequency immediately after the mutator mutation takes place. This prediction explains the wide variation of the DNA G + C content among unicellular organisms and possibly also the wide intragenomic heterogeneity of third codon positions for the genes of multicellular eukaryotes. The present analyses lead to several predictions that are not consistent with a number of the frequently held assumptions in the field of molecular evolution, including belief in a constant rate of evolution, symmetric branching of phylogenetic trees, the generality of higher mutation frequency for neutral sets of nucleotides, the notion that mutator mutations are generally deleterious because of their high mutation rates, and teleological explanations of DNA base composition. Presented at the NATO Advanced Research Workshop onGenome Organization and Evolution, Spetsai, Greece, 16–22 September 1992     

Quantifying nonvertical inheritance in the evolution of Legionella pneumophila

The exchange of genetic material among bacterial strains and species is recognized as an important factor determining their evolutionary, population genetic, and epidemiological features. We present a detailed analysis of nonvertical inheritance in Legionella pneumophila, a human pathogen and facultative intracellular parasite of amoebas. We have analyzed the exchange of L. pneumophila genetic material with other bacteria at three different levels: population genetics, population genomics, and phylogenomics. At the population genetics level, we have analyzed 89 clinical and environmental isolates after sequencing six coding loci and three intergenic regions for a total of 3,923 bp. In the population genomics analysis, we have studied the roles of recombination and mutation in the common portion of the genome sequence of four L. pneumophila strains. In the phylogenomic analysis, we have studied the phylogenetic origin of 1,700 genes in the L. pneumophila pangenome. For this, we have considered 12 possible phylogenetic alternatives, derived from a reference tree obtained from 104 genes from 41 species, which have been tested under a rigorous statistical framework. The results obtained agree in assigning an important role to nonvertical inheritance in shaping the composition of the L. pneumophila genome and of the genetic variation in its populations. We have found a negative correlation between phylogenetic distance and likelihood of horizontal gene transfer. Phylogenetic proximity and increased chances resulting from sharing the ecological niche provided by the amoeba host have likely had a major influence on the rate of gene exchange in Legionella.     

Statistical properties of molecular tree construction methods under the neutral mutation model

Summary The statistical properties of three molecular tree construction methods—the unweighted pair-group arithmetic average clustering (UPG), Farris, and modified Farris methods—are examined under the neutral mutation model of evolution. The methods are compared for accuracy in construction of the topology and estimation of the branch lengths, using statistics of these two aspects. The distribution of the statistic concerning topological construction is shown to be as important as its mean and variance for the comparison.Of the three methods, the UPG method constructs the tree topology with the least variation. The modified Farris method, however, gives the best performance when the two aspects are considered simultaneously. It is also shown that a topology based on two genes is much more accurate than that based on one gene.There is a tendency to accept published molecular trees, but uncritical acceptance may lead one to spurious conclusions. It should always be kept in mind that a tree is a statistical result that is affected strongly by the stochastic error of nucleotide substitution and the error intrinsic to the tree construction method itself.