Selection and Biased Gene Conversion in a Multigene Family: Consequences of Interallelic Bias and Threshold Selection

In a previous paper, I investigated the interactions in a gene family of additive selection and biased gene conversion in a finite population when conversion events are rare. Here I extend my "weak-conversion limit" model by allowing biased interallelic conversion (conversion between alleles at the same locus) of arbitrary frequency and various threshold selection schemes for rare interlocus conversion events. I suggest that it is not unreasonable for gene families to experience threshold fitness functions, and show that certain types of thresholds can greatly constrain the rate at which advantageous alleles are fixed as compared to other fitness schemes, such as additive selection. It is also shown that the double sampling process operating on a gene family in a finite population (sampling over the number of genes in the gene family and over the number of individuals in the population) can have interesting consequences. For selectively neutral alleles that experience interallelic bias, the probability of fixation at each single locus may be essentially neutral, but the cumulative effects on the entire gene family of small departures from neutrality can be significant, especially if the gene family is large. Thus, in some situations, gene families can respond to directional forces that are weak in comparison to drift at single loci.     

Neutral mutation as the source of genetic variation in life history traits

The mechanism underlying the maintenance of adaptive genetic variation is a long-standing question in evolutionary genetics. There are two concepts (mutation-selection balance and balancing selection) which are based on the phenotypic differences between alleles. Mutation - selection balance and balancing selection cannot properly explain the process of gene substitution, i.e. the molecular evolution of quantitative trait loci affecting fitness. I assume that such loci have non-essential functions (small effects on fitness), and that they have the potential to evolve into new functions and acquire new adaptations. Here I show that a high amount of neutral polymorphism at these loci can exist in real populations. Consistent with this, I propose a hypothesis for the maintenance of genetic variation in life history traits which can be efficient for the fixation of alleles with very small selective advantage. The hypothesis is based on neutral polymorphism at quantitative trait loci and both neutral and adaptive gene substitutions. The model of neutral - adaptive conversion (NAC) assumes that neutral alleles are not neutral indefinitely, and that in specific and very rare situations phenotypic (relative fitness) differences between them can appear. In this paper I focus on NAC due to phenotypic plasticity of neutral alleles. The important evolutionary consequence of NAC could be the increased adaptive potential of a population. Loci responsible for adaptation should be fast evolving genes with minimally discernible phenotypic effects, and the recent discovery of genes with such characteristics implicates them as suitable candidates for loci involved in adaptation.     

Interchromosomal Biased Gene Conversion, Mutation and Selection in a Multigene Family

A mathematical model of the effects of interchromosomal biased gene conversion, mutation and natural selection on a multigene family is developed and analyzed. The model assumes two allelic states at each of n loci. The effects of genetic drift are ignored. The model is developed under the assumption of no recombination, but the analysis shows that, at equilibrium, there is no linkage disequilibrium, which implies that the conclusions are valid for arbitrary recombination among loci. At equilibrium, the balance between mutation, gene conversion and selection depends on the ratio of the mutation rates to the quantity [s + g(2α - 1)/ n], where s is the increment or decrement in relative fitness with each additional copy of one of the alleles, g is the conversion rate, and α is a measure of the bias in favor of one of the alleles. When this quantity is large relative to the mutation rates, the allele that has the net advantage, combining the effects of selection and conversion, will be nearly fixed in the multigene family. A comparison of these results with those from a comparable model of intrachromosomal biased conversion shows that biased interchromosomal conversion leads to approximately the same equilibrium copy number as does intrachromosomal conversion of the same strength. Interchromosomal conversion is much more effective in causing the substitution of one allele by another. The relative frequencies of interchromosomal and intrachromosomal conversion is indicated by the extent of the linkage disequilibrium among the loci in a multigene family.     

SHORT INDELS ARE SUBJECT TO INSERTION‐BIASED GENE CONVERSION

Recombination between homologous loci is accompanied by formation of heteroduplexes. Repairing mismatches in heteroduplexes often leads to single nucleotide substitutions in a process known as gene conversion. Gene conversion was shown to be GC‐biased in different organisms; that is, a W(A or T)→S(G or C) substitution is more likely in this process than a S→W substitution. Here, we show that the insertion/deletion ratio for short noncoding indels that reach fixation between species is positively correlated with the recombination rate in Drosophila melanogaster, Homo sapiens, and Saccharomyces cerevisiae. This correlation is both due to an increase of the fixation rate of insertions and decrease of the fixation rate of deletions in regions of high recombination. Whole‐genome data on indel polymorphism and divergence in D. melanogaster rule out mutation biases and selection as the cause of this trend, pointing to insertion‐biased gene conversion as the most likely explanation. The bias toward insertions is the strongest for single‐nucleotide indels, and decreases with indel length. In regions of high recombination rate this bias leads to an up to ～5‐fold excess of fixed short insertions over deletions, and substantially affects the evolution of DNA segments.     

Mating system and recombination affect molecular evolution in four Triticeae species.

Mating systems and recombination are thought to have a deep impact on the organization and evolution of genomes. Because of the decline in effective population size and the interference between linked loci, the efficacy of selection is expected to be reduced in regions with low recombination rates and in the whole genome of self-fertilizing species. At the molecular level, relaxed selection is expected to result in changes in the rate of protein evolution and the pattern of codon bias. It is increasingly recognized that recombination also affects non-selective processes such as the biased gene conversion towards GC alleles (bGC). Like selection, this kind of meiotic drive in favour of GC over AT alleles is expected to be reduced in weakly recombining regions and genomes. Here, we investigated the effect of mating system and recombination on molecular evolution in four Triticeae species: two outcrossers (Secale cereale and Aegilops speltoides) and two selfers (Triticum urartu and Triticum monococcum). We found that GC content, possibly driven by bGC, is affected by mating system and recombination as theoretically predicted. Selection efficacy, however, is only weakly affected by mating system and recombination. We investigated the possible reasons for this discrepancy. A surprising one is that, in outcrossing lineages, selection efficacy could be reduced because of high substitution rates in favour of GC alleles. Outcrossers, but not selfers, would thus suffer from a 'GC-induced' genetic load. This result sheds new light on the evolution of mating systems.     

Substitution rates in Drosophila nuclear genes: implications for translational selection

The relationships between synonymous and nonsynonymous substitution rates and between synonymous rate and codon usage bias are important to our understanding of the roles of mutation and selection in the evolution of Drosophila genes. Previous studies used approximate estimation methods that ignore codon bias. In this study we reexamine those relationships using maximum-likelihood methods to estimate substitution rates, which accommodate the transition/transversion rate bias and codon usage bias. We compiled a sample of homologous DNA sequences at 83 nuclear loci from Drosophila melanogaster and at least one other species of Drosophila. Our analysis was consistent with previous studies in finding that synonymous rates were positively correlated with nonsynonymous rates. Our analysis differed from previous studies, however, in that synonymous rates were unrelated to codon bias. We therefore conducted a simulation study to investigate the differences between approaches. The results suggested that failure to properly account for multiple substitutions at the same site and for biased codon usage by approximate methods can lead to an artifactual correlation between synonymous rate and codon bias. Implications of the results for translational selection are discussed.     

Adaptation or biased gene conversion? Extending the null hypothesis of molecular evolution

The analysis of evolutionary rates is a popular approach to characterizing the effect of natural selection at the molecular level. Sequences contributing to species adaptation are expected to evolve faster than nonfunctional sequences because favourable mutations have a higher fixation probability than neutral ones. Such an accelerated rate of evolution might be due to factors other than natural selection, in particular GC-biased gene conversion. This is true of neutral sequences, but also of constrained sequences, which can be illustrated using the mouse Fxy gene. Several criteria can discriminate between the natural selection and biased gene conversion models. These criteria suggest that the recently reported human accelerated regions are most likely the result of biased gene conversion. We argue that these regions, far from contributing to human adaptation, might represent the Achilles' heel of our genome.     

The rate of synonymous substitution in enterobacterial genes is inversely related to codon usage bias

Genes sequences from Escherichia coli, Salmonella typhimurium, and other members of the Enterobacteriaceae show a negative correlation between the degree of synonymous-codon usage bias and the rate of nucleotide substitution at synonymous sites. In particular, very highly expressed genes have very biased codon usage and accumulate synonymous substitutions very slowly. In contrast, there is little correlation between the degree of codon bias and the rate of protein evolution. It is concluded that both the rate of synonymous substitution and the degree of codon usage bias largely reflect the intensity of selection at the translational level. Because of the high variability among genes in rates of synonymous substitution, separate molecular clocks of synonymous substitution might be required for different genes.      

Effect of balancing selection on spatial genetic structure within populations: theoretical investigations on the self-incompatibility locus and empirical studies in Arabidopsis halleri

The effect of selection on patterns of genetic structure within and between populations may be studied by contrasting observed patterns at the genes targeted by selection with those of unlinked neutral marker loci. Local directional selection on target genes will produce stronger population genetic structure than at neutral loci, whereas the reverse is expected for balancing selection. However, theoretical predictions on the intensity of this signal under precise models of balancing selection are still lacking. Using negative frequency-dependent selection acting on self-incompatibility systems in plants as a model of balancing selection, we investigated the effect of such selection on patterns of spatial genetic structure within a continuous population. Using numerical simulations, we tested the effect of the type of self-incompatibility system, the number of alleles at the self-incompatibility locus and the dominance interactions among them, the extent of gene dispersal, and the immigration rate on spatial genetic structure at the selected locus and at unlinked neutral loci. We confirm that frequency-dependent selection is expected to reduce the extent of spatial genetic structure as compared to neutral loci, particularly in situations with low number of alleles at the self-incompatibility locus, high frequency of codominant interactions among alleles, restricted gene dispersal and restricted immigration from outside populations. Hence the signature of selection on spatial genetic structure is expected to vary across species and populations, and we show that empirical data from the literature as well as data reported here on three natural populations of the herb Arabidopsis halleri confirm these theoretical results.     

Effect of parasitic sex-ratio distorters on host gene frequencies in a mainland-island context

It was previously argued that infection by parasitic sex-ratio distorters can enhance both random genetic drift and genetic influx from outside the population. However, these two enhancement effects have been studied independently. Here, we study the equilibrium frequencies of alleles (neutral and selected) in a mainland-island scenario where both genetic drift and genetic influx are enhanced due to infection by a cytoplasmic feminizing element. Interestingly, our model reveals that at neutral loci, the two effects almost exactly cancel each other out, such that infection has only a very minor effect on the equilibrium frequency distributions of alleles. At selected loci, in contrast, the two effects are unbalanced and infection has conspicuous effects. Despite the cryptic effects of infection at neutral loci, we demonstrate that temporally spaced data can be used to evaluate the effect of infection on genetic drift and that on gene flow separately.     

Effect of strong directional selection on weakly selected mutations at linked sites: implication for synonymous codon usage

The fixation of weakly selected mutations can be greatly influenced by strong directional selection at linked loci. Here, I investigate a two-locus model in which weakly selected, reversible mutations occur at one locus and recurrent strong directional selection occurs at the other locus. This model is analogous to selection on codon usage at synonymous sites linked to nonsynonymous sites under strong directional selection. Two approximations obtained here describe the expected frequency of the weakly selected preferred alleles at equilibrium. These approximations, as well as simulation results, show that the level of codon bias declines with an increasing rate of substitution at the strongly selected locus, as expected from the well-understood theory that selection at one locus reduces the efficacy of selection at linked loci. These solutions are used to examine whether the negative correlation between codon bias and nonsynonymous substitution rates recently observed in Drosophila can be explained by this hitchhiking effect. It is shown that this observation can be reasonably well accounted for if a large fraction of the nonsynonymous substitutions on genes in the data set are driven by strong directional selection.     

The anomalous effects of biased mutation

A model is presented in which alleles at a number of loci combine to influence the value of a quantitative trait that is subject to stabilizing selection. Mutations can occur to alleles at the loci under consideration. Some of these mutations will tend to increase the value of the trait, while others will tend to decrease it. In contrast to most previous models, we allow the mean effect of mutations to be nonzero. This means that, on average, mutations can have a bias, such that they tend to either increase or decrease the value of the trait. We find, unsurprisingly, that biased mutation moves the equilibrium mean value of the quantitative trait in the direction of the bias. What is more surprising is the behavior of the deviation of the equilibrium mean value of the trait from its optimal value. This has a nonmonotonic dependence on the degree of bias, so that increasing the degree of bias can actually bring the mean phenotype closer to the optimal phenotype. Furthermore, there is a definite maximum to the extent to which biased mutation can cause a difference between the mean phenotype and the optimum. For plausible parameter values, this maximum-possible difference is small. Typically, quantitative-genetics models assume an unconstrained model of mutation, where the expected difference in effect between a parental allele and a mutant allele is independent of the current state of the parental allele. Our results show that models of this sort can easily lead to biologically implausible consequences when mutations are biased. In particular, unconstrained mutation typically leads to a continual increase or decrease in the mean allelic effects at all trait-controlling loci. Thus at each of these loci, the mean allelic effect eventually becomes extreme. This suggests that some of the models of mutation most commonly used in quantitative genetics should be modified so as to introduce genetic constraints.     

Multilocus patterns of polymorphism and selection across the X chromosome of Caenorhabditis remanei

Natural selection and neutral processes such as demography, mutation, and gene conversion all contribute to patterns of polymorphism within genomes. Identifying the relative importance of these varied components in evolution provides the principal challenge for population genetics. To address this issue in the nematode Caenorhabditis remanei, I sampled nucleotide polymorphism at 40 loci across the X chromosome. The site-frequency spectrum for these loci provides no evidence for population size change, and one locus presents a candidate for linkage to a target of balancing selection. Selection for codon usage bias leads to the non-neutrality of synonymous sites, and despite its weak magnitude of effect (N(e)s approximately 0.1), is responsible for profound patterns of diversity and divergence in the C. remanei genome. Although gene conversion is evident for many loci, biased gene conversion is not identified as a significant evolutionary process in this sample. No consistent association is observed between synonymous-site diversity and linkage-disequilibrium-based estimators of the population recombination parameter, despite theoretical predictions about background selection or widespread genetic hitchhiking, but genetic map-based estimates of recombination are needed to rigorously test for a diversity-recombination relationship. Coalescent simulations also illustrate how a spurious correlation between diversity and linkage-disequilibrium-based estimators of recombination can occur, due in part to the presence of unbiased gene conversion. These results illustrate the influence that subtle natural selection can exert on polymorphism and divergence, in the form of codon usage bias, and demonstrate the potential of C. remanei for detecting natural selection from genomic scans of polymorphism.     

Sojourn times and substitution rate at overdominant and linked neutral loci

The sojourn times until fixation of an overdominant allele were investigated based on the diffusion equation. Furthermore, the rate of accumulation of mutations, or the substitution rate, was predicted from the mean extinction time of a common overdominant allele. The substitution rate calculated theoretically agreed well with that determined by computer simulation. Overdominant selection enhances the polymorphism at linked loci, while its effect on the sojourn times and the substitution rate at linked loci has not been studied yet. To solve these problems, a model that assumed two linked loci, each with infinite alleles, was examined by computer simulation. A decrease in the recombination rate between two loci markedly changed the distribution of sojourn times of a neutral allele. Although overdominant selection obviously increased the sojourn times and the polymorphism at a linked locus, the rate of nucleotide substitution at the neutral locus was not influenced significantly even if complete linkage was assumed. These results suggest that, in regions containing overdominant genes, linked neutral loci will exhibit elevated levels of polymorphism, but their rate of molecular evolution remains that predicted by neutral theory.     

Evolution of MHC class IIB in the genome of wild and ornamental guppies, Poecilia reticulata

This is the first study to quantify genomic sequence variation of the major histocompatibility complex (MHC) in wild and ornamental guppies, Poecilia reticulata. We sequenced 196-219 bp of exon 2 MHC class IIB (DAB) in 56 wild Trinidadian guppies and 14 ornamental strain guppies. Each of two natural populations possessed high allelic richness (15-16 alleles), whereas only three or fewer DAB alleles were amplified from ornamental guppies. The disparity in allelic richness between wild and ornamental fish cannot be fully explained by fixation of alleles by inbreeding, nor by the presence of non-amplified sequences (ie null alleles). Rather, we suggest that the same allele is fixed at duplicated MHC DAB loci owing to gene conversion. Alternatively, the number of loci in the ornamental strains has contracted during >100 generations in captivity, a hypothesis consistent with the accordion model of MHC evolution. We furthermore analysed the substitution patterns by making pairwise comparisons of sequence variation at the putative peptide binding region (PBR). The rate of non-synonymous substitutions (dN) only marginally exceeded synonymous substitutions (dS) in PBR codons. Highly diverged sequences showed no evidence for diversifying selection, possibly because synonymous substitutions have accumulated since their divergence. Also, the substitution pattern of similar alleles did not show evidence for diversifying selection, plausibly because advantageous non-synonymous substitutions have not yet accumulated. Intermediately diverged sequences showed the highest relative rate of non-synonymous substitutions, with dN/dS>14 in some pairwise comparisons. Consequently, a curvilinear relationship was observed between the dN/dS ratio and the level of sequence divergence.     

Evidence of selection on silent site base composition in mammals: potential implications for the evolution of isochores and junk DNA.

It has been suggested that mutation bias is the major determinant of base composition bias at synonymous, intron, and flanking DNA sites in mammals. Here I test this hypothesis using population genetic data from the major histocompatibility genes of several mammalian species. The results of two tests are inconsistent with the mutation hypothesis in coding, noncoding, CpG-island, and non-CpG-island DNA, but are consistent with selection or biased gene conversion. It is argued that biased gene conversion is unlikely to affect silent site base composition in mammals. The results therefore suggest that selection is acting upon silent site G + C content. This may have broad implications, since silent site base composition reflects large-scale variation in G + C content along mammalian chromosomes. The results therefore suggest that selection may be acting upon the base composition of isochores and large sections of junk DNA.     

Can mutation or fixation biases explain the allele frequency distribution of human single nucleotide polymorphisms (SNPs)?

One of the most abiding controversies in evolutionary biology concerns the role of neutral processes in molecular evolution. A main focus of the debate has been the evolution of isochores, the strong and systematic variation of base composition in mammalian genomes. One set of hypotheses argue that regions of similar GC are owing to localised mutational biases coupled with neutral evolution. The alternatives point to either selection or biased gene conversion as mechanisms to preferentially remove A or T bases, favouring G and C instead. Using a novel method, we compare models including such fixation biases to models based on mutation bias alone, under the assumption that non-coding, non-repetitive human DNA is at compositional equilibrium. While failing to fully explain the allele frequency distributions of recent single nucleotide polymorphism data, we show that the data are best fitted if the mutation bias is assumed to be constant across the genome, while fixation bias varies with GC content. We also attempt to estimate the strength of fixation bias, which increases linearly with increasing GC. Our approximation suggests that this force exists within the necessary parameter range: it is not so weak as to be drowned by random drift, but not so strong as to lead to exclusive use of G and C alone. Together these results demonstrate that mutation bias fails to explain the evolution of isochores, and suggest that either selection or biased gene conversion are involved.     

The Correlation Between Recombination Rate and Dinucleotide Bias in Drosophila melanogaster

Revealing how recombination affects genomic sequence is of great significance to our understanding of genome evolution. The present paper focuses on the correlation between recombination rate and dinucleotide bias in Drosophila melanogaster genome. Our results show that the overall dinucleotide bias is positively correlated with recombination rate for genomic sequences including untranslated regions, introns, intergenic regions, and coding sequences. The correlation patterns of individual dinucleotide biases with recombination rate are presented. Possible mechanisms of interaction between recombination and dinucleotide bias are discussed. Our data indicate that there may be a genome-wide universal mechanism acting between recombination rate and dinucleotide bias, which is likely to be neighbor-dependent biased gene conversion.     

HYBRID ZONES WITH DOBZHANSKY-TYPE EPISTATIC SELECTION

Dobzhansky's model of epistatic selection assumes that viable genotypes form “clusters” in genotype space so that populations can evolve from one state to a reproductively isolated state following a “ridge” of well-fit genotypes without crossing any deep adaptive valleys. Recently, the importance of Dobzhansky-type models in evolutionary studies has been reemphasized by Gavrilets (1997a) and Gavrilets and Gravner (1997) who argue that the existence of “ridges” of well-fit genotypes connecting reproductively isolated genotypes is actually a general property of multidimensional adaptive landscapes. Using rigorous techniques and numerical simulations, I analyze clines in the frequencies of selected and neutral alleles maintained by a balance of migration and Dobzhansky-type epistatic selection acting on two diallelic loci. I show that Dobzhansky-type epistatic selection can build up a very strong barrier to neutral gene flow. I describe properties of clines that are indicative of Dobzhansky-type selection.