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.     

Role of gene duplication in evolution

It is now known that many multigene and supergene families exist in eukaryote genomes: multigene families with uniform copy members like genes for ribosomal RNA, those with variable members like immunoglobulin genes, and supergene families such as those for various growth factor and hormone receptors. Many such examples indicate that gene duplication and subsequent differentiation are extremely important for organismal evolution. In particular, gene duplication could well have been the primary mechanism for the evolution of complexity in higher organisms. Population genetic models for the origin of gene families with diverse functions are presented, in which natural selection favors those genomes with more useful mutants in duplicated genes. Since any gene has a certain probability of degenerating by mutation, success versus failure in acquiring a new gene by duplication may be expressed as the ratio of probabilities of spreading of useful versus detrimental mutations in redundant gene copies. Also examined are the effects of gene duplication on evolution by compensatory advantageous mutations. Results of the analyses show that both natural selection and random drift are important for the origin of gene families. In addition, interaction between molecular mechanisms such as unequal crossing-over and gene conversion, and selection or drift is found to have a large effect on evolution by gene duplication.     

Pleiotropic Models of Polygenic Variation, Stabilizing Selection, and Epistasis

We show that in polymorphic populations many polygenic traits pleiotropically related to fitness are expected to be under apparent ``stabilizing selection' independently of the real selection acting on the population. This occurs, for example, if the genetic system is at a stable polymorphic equilibrium determined by selection and the nonadditive contributions of the loci to the trait value either are absent, or are random and independent of those to fitness. Stabilizing selection is also observed if the polygenic system is at an equilibrium determined by a balance between selection and mutation (or migration) when both additive and nonadditive contributions of the loci to the trait value are random and independent of those to fitness. We also compare different viability models that can maintain genetic variability at many loci with respect to their ability to account for the strong stabilizing selection on an additive trait. Let V(m) be the genetic variance supplied by mutation (or migration) each generation, V(g) be the genotypic variance maintained in the population, and n be the number of the loci influencing fitness. We demonstrate that in mutation (migration)-selection balance models the strength of apparent stabilizing selection is order V(m)/V(g). In the overdominant model and in the symmetric viability model the strength of apparent stabilizing selection is approximately 1/(2n) that of total selection on the whole phenotype. We show that a selection system that involves pairwise additive by additive epistasis in maintaining variability can lead to a lower genetic load and genetic variance in fitness (approximately 1/(2n) times) than an equivalent selection system that involves overdominance. We show that, in the epistatic model, the apparent stabilizing selection on an additive trait can be as strong as the total selection on the whole phenotype.     

The evolutionary rate of duplicated genes under concerted evolution

The effect of directional selection on the fixation process of a single mutation that spreads in a multigene family by gene conversion is investigated. A simple two-locus model with two alleles, A and a, is first considered in a random-mating diploid population with size N. There are four haplotypes, AA, Aa, aA, and aa, and selection works on the number of alleles A in a diplod (i = 0, 1, 2, 3, 4). Because gene conversion is allowed between the two loci, when the mutation rate is very low, either AA or aa will fix in the population eventually. We consider a situation where a single mutant, A, arises in one locus when a is fixed in both loci. Then, we derive the fixation probability analytically, and the fixation time is investigated by simulations. It is found that gene conversion has an effect to increase the "effective" population size, so that weak selection works more efficiently in a multigene family. With these results, we discuss the effect of gene conversion on the rate of molecular evolution in a multigene family undergoing concerted evolution. We also argue about the applicability of the theoretical results to models of multigene families with more than two loci.     

A model of evolution for accumulating genetic information

By taking into account recent knowledge of multigene families and other repetitive DNA sequences, a model of evolution by gene duplication for accumulating genetic information is studied. Genetic information is defined as the sum of distinct functions that the gene family can perform. A coefficient, "genetic diversity" is defined and used in this study, that is highly correlated with genetic information. Initially, a multigene family with a few gene copies is assumed, and natural selection starts to work on this gene family to increase genetic diversity contained in the gene family. As an important mechanism, unequal crossing-over is incorporated. Together with mutation, it is responsible for supplying genetic variability among individuals for selection to work. A specific model, in which individuals with less genetic diversity are selectively disadvantageous, has been studied in detail. Through approximate theoretical analysis and extensive Monte Carlo studies, it has been shown that the system is an extremely efficient way to accumulate genetic information. For attaining one gene, the genetic load is much smaller under this model than under the traditional model of natural selection. The model may be applied to the process of origin of multigene families with diverse copy members such as those of immunoglobulin or cytochrome P450. In general, the process of creating new genes by duplication might be somewhere between the present and the traditional models.     

Simulation Study of a Multigene Family, with Special Reference to the Evolution of Compensatory Advantageous Mutations

We investigate the evolution of a multigene family incorporating the forces of drift, mutation, gene conversion, unequal crossing over and selection. The use of simulation studies is required due to the complexity of the model. Selection is modeled in two modes: positive selection as a function of the number of different beneficial alleles and negative selection against deleterious alleles. We assume that gene conversion is unbiased, and that all mutations are initially deleterious. Compensation between mutants creates beneficial and neutral alleles, and allowances are made for compensatory mutations either within or between the members of a multigene family. We find that gene conversion can enhance the rate of acquisition of compensatory advantageous mutations when genes are redundant.     

FURTHER SIMULATION STUDIES ON EVOLUTION BY GENE DUPLICATION

In order to understand the origin of multigene families, Monte Carlo simulations were performed to see how a genetic system evolves under unequal crossing-over, mutation, random genetic drift and natural selection, starting from a single gene copy. Both haploid and diploid models were examined. Beneficial, neutral, and detrimental mutations were incorporated, and “positive” selection favors those chromosomes (haploid) or individuals (diploid) with more beneficial mutations than others. The same model for haploids was previously investigated with special reference to the evolution of gene organization, and the ratio of the numbers of beneficial genes to pseudogenes was found to be a rough indicator of the relative strengths of positive and negative (against deleterious alleles) natural selection (Ohta, 1987b). In the present paper, the evolution of gene organization and of sequence divergence among genes in the multigene family is examined. It is shown that positive selection accelerates the accumulation of arrays containing different beneficial mutations, but that total divergence including both neutral and beneficial mutations is not very sensitive to positive selection, under this model. The proportion of beneficial mutations in the total mutations accumulated is a better indicator of positive selection than is the total divergence. It is pointed out that various observed examples in which amino-acid substitutions are accelerated, as compared with synonymous substitutions in duplicated genes (Li, 1985), may reflect the effect of selection similar to the present scheme. The diploid model is shown to be more efficient for accumulating beneficial mutations in duplicated genes than the haploid one, and the relevance of this finding to the advantage of sexual reproduction is discussed.     

Mutation-selection balance: ancestry,load, and maximum principle

We analyze the equilibrium behavior of deterministic haploid mutation-selection models. To this end, both the forward and the time-reversed evolution processes are considered. The stationary state of the latter is called the ancestral distribution, which turns out as a key for the study of mutation-selection balance. We find that the ancestral genotype frequencies determine the sensitivity of the equilibrium mean fitness to changes in the corresponding fitness values and discuss implications for the evolution of mutational robustness. We further show that the difference between the ancestral and the population mean fitness, termed mutational loss, provides a measure for the sensitivity of the equilibrium mean fitness to changes in the mutation rate. The interrelation of the loss and the mutation load is discussed. For a class of models in which the number of mutations in an individual is taken as the trait value, and fitness is a function of the trait, we use the ancestor formulation to derive a simple maximum principle, from which the mean and variance of fitness and the trait may be derived; the results are exact for a number of limiting cases, and otherwise yield approximations which are accurate for a wide range of parameters. These results are applied to threshold phenomena caused by the interplay of selection and mutation (known as error thresholds). They lead to a clarification of concepts, as well as criteria for the existence of error thresholds.     

Theoretical population genetics of repeated genes forming a multigene family

The evolution of repeated genes forming a multigene family in a finite population is studied with special reference to the probability of gene identity, i.e., the identity probability of two gene units chosen from the gene family. This quantity is called clonality and is defined as the sum of squares of the frequencies of gene lineages in the family. The multigene family is undergoing continuous unequal somatic crossing over, ordinary interchromosomal crossing over, mutation and random frequency drift. Two measures of clonality are used: clonality within one chromosome and that between two different chromosomes. The equilibrium properties of the means, the variances and the covariance of the two measures of clonality are investigated by using the diffusion equation method under the assumption of constant number of gene units in the multigene family. Some models of natural selection based on clonality are considered. The possible significance of the variance and covariance of clonality among the chromosomes on the adaptive differentiation of gene families such as those producing antibodies is discussed.     

The coalescent with selection on copy number variants

We develop a coalescent-based simulation tool to generate patterns of single nucleotide polymorphisms (SNPs) in a wide region encompassing both the original and duplicated genes. Selection on the new duplicated copy and interlocus gene conversion between the two copies are incorporated. This simulation enables us to explore how selection on duplicated copies affects the pattern of SNPs. The fixation of an advantageous duplicated copy causes a strong reduction in polymorphism not only in the duplicated copy but also in its flanking regions, which is a typical signature of a selective sweep by positive selection. After fixation, polymorphism gradually increases by accumulating neutral mutations and eventually reaches the equilibrium value if there is no gene conversion. When gene conversion is active, the number of SNPs in the duplicated copy quickly increases by transferring SNPs from the original copy; therefore, the time when we can recognize the signature of selection is decreased. Because this effect of gene conversion is restricted only to the duplicated region, more power to detect selection is expected if a flanking region to the duplicated copy is used.     

A Model for Damage Load and Its Implications for the Evolution of Bacterial Aging

Deleterious mutations appearing in a population increase in frequency until stopped by natural selection. The ensuing equilibrium creates a stable frequency of deleterious mutations or the mutational load. Here I develop the comparable concept of a damage load, which is caused by harmful non-heritable changes to the phenotype. A damage load also ensues when the increase of damage is opposed by selection. The presence of a damage load favors the evolution of asymmetrical transmission of damage by a mother to her daughters. The asymmetry is beneficial because it increases fitness variance, but it also leads to aging or senescence. A mathematical model based on microbes reveals that a cell lineage dividing symmetrically is immortal if lifetime damage rates do not exceed a threshold. The evolution of asymmetry allows the lineage to persist above the threshold, but the lineage becomes mortal. In microbes with low genomic mutation rates, it is likely that the damage load is much greater than the mutational load. In metazoans with higher genomic mutation rates, the damage and the mutational load could be of the same magnitude. A fit of the model to experimental data shows that Escherichia coli cells experience a damage rate that is below the threshold and are immortal under the conditions examined. The model estimates the asymmetry level of E. coli to be low but sufficient for persisting at higher damage rates. The model also predicts that increasing asymmetry results in diminishing fitness returns, which may explain why the bacterium has not evolved higher asymmetry.     

Recombination,selection, and the evolution of tandem gene arrays

Multigene families—immunity genes or sensory receptors, for instance—are often subject to diversifying selection. Allelic diversity may be favored not only through balancing or frequency-dependent selection at individual loci but also by associating different alleles in multicopy gene families. Using a combination of analytical calculations and simulations, we explored a population genetic model of epistatic selection and unequal recombination, where a trade-off exists between the benefit of allelic diversity and the cost of copy abundance. Starting from the neutral case, where we showed that gene copy number is Gamma distributed at equilibrium, we derived also the mean and shape of the limiting distribution under selection. Considering a more general model, which includes variable population size and population substructure, we explored by simulations mean fitness and some summary statistics of the copy number distribution. We determined the relative effects of selection, recombination, and demographic parameters in maintaining allelic diversity and shaping the mean fitness of a population. One way to control the variance of copy number is by lowering the rate of unequal recombination. Indeed, when encoding recombination by a rate modifier locus, we observe exactly this prediction. Finally, we analyzed the empirical copy number distribution of 3 genes in human and estimated recombination and selection parameters of our model.     

Evolution of Bacterial Transformation: Is Sex with Dead Cells Ever Better than No Sex at All?

Computer simulations of bacterial transformation are used to show that, under a wide range of biologically reasonable assumptions, transforming populations undergoing deleterious mutation and selection have a higher mean fitness at equilibrium than asexual populations. The source of transforming DNA, the amount of DNA taken up by each transforming cell, and the relationship between number of mutations and cell viability (the fitness function) are important factors. When the DNA source is living cells, transformation resembles meiotic sex. When the DNA source is cells killed by selection against mutations, transformation increases the average number of mutations per genome but can nevertheless increase the mean fitness of the population at equilibrium. In a model of regulated transformation, in which the most fit cells of a transforming population do not transform, transforming populations are always fitter at equilibrium than asexual populations. These results show that transformation can reduce mutation load.     

Colony-level selection in the social insects: Single locus additive and nonadditive models

Genetic models of colony-level selection applicable to diploids (termites) and haplodiploids (social Hymenoptera) are analysed. In the Additive model colony fitnesses are just the arithmetic average of the contribution of the worker genotypes. In the Nonadditive model the fitness of the heterogenotypic colonies (those comprised of more than one worker genotype) may be altered due to interaction between the different worker genotypes. This is modelled by multiplying the additive fitness by the variable, e_i. With additive selection the same equilibrium gene frequency occurs in diploids and in haplodiploids with both once and twice mated queens. In haplodiploids if selection is nonadditive and strong, up to three polymorphic equilibria can exist; however, only a maximum of two are possible with weak selection. Multiple mating by queens increases the number of equilibria possible. Worker-produced males alter the conditions for the existence of a polymorphic equilibrium, and shift the male and female equilibrium gene frequencies.     

The balance between pleiotropic mutation and selection,when alleles have discrete effects

The theory of pleiotropic mutation and selection is investigated and developed for a large population of asexual organisms. Members of the population are subject to stabilising selection on Omega phenotypic characters, which each independently affect fitness. Pleiotropy is incorporated into the model by allowing each mutation to simultaneously affect all characters. To expose differences with continuous-allele models, the characters are taken to originate from discrete-effect alleles and thus have discrete genotypic effects. Each character can take the values nxDelta where n=0,+/-1,+/-2, em leader, and the splitting in character effects, Delta, is a parameter of the model. When the distribution of mutant effects is normally distributed around the parental value, and Delta is large, a "stepwise" model of mutation arises, where only adjacent trait effects are accessible from a single mutation. The present work is primarily concerned with the opposite limit, where Delta is small and many different trait effects are accessible from a single mutation.In contrast to what has been established for continuous-effect models, discrete-effect models do not yield a singular equilibrium distribution of genotypic effects for any value of Omega. Instead, for different values of Omega, the equilibrium frequencies of trait values have very different dependencies on Delta. For Omega=1 and 2, decreasing Delta broadens the width of the frequency distribution and hence increases the equilibrium level of polymorphism. For all sufficiently large values of Omega, however, decreasing Delta decreases the width of the frequency distribution and the equilibrium level of polymorphism. The connection with continuous trait models follows when the limit Delta-->0 is considered, and a singular probability density of trait values is obtained for all sufficiently large Omega.     

King''s Formula for the Mutation Load with Epistasis

A formula by J. L. King gives the equilibrium mutation load as L = 2 sigma ui(1 - qi)/z - x) in which ui is the mutation rate to deleterious alleles at the ith locus, qi is the frequency of mutant alleles at this locus, x is the mean number of such mutant genes per individual before selection, z is the mean number in individuals eliminated by selection, and the summation is over all relevant loci. We show that this rule is inaccurate for intense selection and that a correct formula is L = 2 sigma ui(1 - qi) w/(z - x) = 2U w/(z - x) = 2U/(z - x + 2U) in which U is the mean number of new mutations per haploid genome in the population and w is the mean relative fitness before selection. If w/(z - x) less than 1/2, the mutation load is less than the Haldane value (U less than or equal to L less than or equal to 2U) and can be considerably less. In a diploid asexual population, however, with independent occurrence of mutations, L = 1 - e-2U regardless of the mode of selection.     

The Impact of Message Mutation on the Fitness of a Genetic Code

The Standard Genetic Code is organized such that similar codons encode similar amino acids. One explanation suggested that the Standard Code is the result of natural selection to reduce the fitness ``load' that derives from the mutation and mistranslation of protein-coding genes. We review the arguments against the mutational load-minimizing hypothesis and argue that they need to be reassessed. We review recent analyses of the organization of the Standard Code and conclude that under cautious interpretation they support the mutational load-minimizing hypothesis. We then present a deterministic asexual model with which we study the mode of selection for load minimization. In this model, individual fitness is determined by a protein phenotype resulting from the translation of a mutable set of protein-coding genes. We show that an equilibrium fitness may be associated with a population with the same genetic code and that genetic codes that assign similar codons to similar amino acids have a higher fitness. We also show that the number of mutant codons in each individual at equilibrium, which determines the strength of selection for load minimization, reflects a long-term evolutionary balance between mutations in messages and selection on proteins, rather than the number of mutations that occur in a single generation, as has been assumed by previous authors. We thereby establish that selection for mutational load minimization acts at the level of an individual in a single generation. We conclude with comments on the shortcomings and advantages of load minimization over other hypotheses for the origin of the Standard Code. Received: 4 April 2001 / Accepted: 22 October 2001     

Social selection in human populations: differential modification of the fitness by the sex of an affected parent

Social selection model of two alleles at a locus has been constructed where the fitness of offspring is modified differentially by the sex of an affected parent. It has been shown that the equilibrium frequency of deleterious alleles depends on the loss of fitness of an individual due to the trait (gamma), due to an affected mother (beta 1) and affected father (beta 2), and the probability that the heterozygote develops the trait (h). The equilibrium gene frequency is given approximately by alpha/hs for 0 less than h less than or equal to 1 and square root alpha/s for h = 0, where s = gamma + (beta 1 + beta 2) (1 - gamma)/2 for multiplicative models. The same result holds for the additive model if the value of gamma is small.     

Evolution of fitnesses and allele frequencies in a population with spatially heterogeneous selection pressures

The level of gene flow considerably influences the outcome of evolutionary processes in structured populations with spatial heterogeneity in selection pressures; low levels of gene flow may allow local adaptation whereas high levels of gene flow may oppose this process thus preventing the stable maintenance of polymorphism. Indeed, proportions of fitness space that successfully maintain polymorphism are substantially larger in spatially heterogenous populations with lower to moderate levels of gene flow when compared to single-deme models. Nevertheless, the effect of spatial heterogeneity on the evolutionary construction of polymorphism is less clear. We have investigated the levels of polymorphism resulting from a simple two-deme construction model, which incorporates recurrent mutation as well as selection. We further compared fitness properties, stability of equilibria, and frequency distribution patterns emerging from the construction approach and compared these to the static fitness-space approach. The construction model either promotes or constrains the level of polymorphisms, depending on the levels of gene flow. Comparison of the fitness properties resulting from both approaches shows that they maintain variation in different parts of fitness space. The part of fitness space resulting from construction is more stable than that implied by the ahistoric fitness-space approach. Finally, the equilibrium allele-frequency distribution patterns vary substantially with different levels of gene flow, underlining the importance of correctly sampling spatial structure if these patterns are to be used to estimate population-genetic processes.     

Molecular Evolution in a Multisite Nearly Neutral Mutation Model

A simple nearly neutral mutation model of protein evolution was studied using computer simulation assuming a constant population size. In this model, a gene consists of a finite number of codons and there is no recombination within a gene. Each codon has two replacement and one silent sites. The fitness of a gene was determined multiplicatively by amino acids specified by codons (the independent multicodon model). Nucleotide diversity at replacement sites decreases as selection becomes stronger. A reduction of nucleotide diversity at silent sites also occurs as selection intensifies but the magnitude of the reduction is not a monotone function of the intensity of selection. The dispersion index is close to one. The average value of Tajima's and Fu and Li's statistics are negative and their absolute values increases as selection intensifies. However, their powers of detecting selection under the present model were not high unless the number of sites is large or mutation rate is high. The MK test was shown to detect intermediate selection fairly well. For comparison, the house-of-cards model was also investigated and its behavior was shown to be more sensitive to changes of population size than that of the independent multicodon model. The relevance of the present model for explaining protein evolution was discussed comparing its prediction and recent DNA data. Received: 24 May 1999 / Accepted: 17 August 1999