Nuclear–mitochondrial epistasis: a gene's eye view of genomic conflict

We use population genetic models to investigate the cooperative and conflicting synergistic fitness effects between genes from the nucleus and the mitochondrion. By varying fitness parameters, we examine the scope for conflict relative to cooperation among genomes and the utility of the “gene's eye view” analytical approach, which is based on the marginal average fitness of specific alleles. Because sexual conflict can maintain polymorphism of mitochondrial haplotypes, we can explore two types of evolutionary conflict (genomic and sexual) with one epistatic model. We find that the nuclear genetic architecture (autosomal, X‐linked, or Z‐linked) and the mating system change the regions of parameter space corresponding to the evolution by sexual and genomic conflict. For all models, regardless of conflict or cooperation, we find that population mean fitness increases monotonically as evolution proceeds. Moreover, we find that the process of gene frequency change with positive, synergistic fitnesses is self‐accelerating, as the success of an allele in one genome or in one sex increases the frequency of the interacting allele upon which its success depends. This results in runaway evolutionary dynamics caused by the positive intergenomic associations generated by selection. An inbreeding mating system tends to further accelerate these runaway dynamics because it maintains favorable host–symbiont or male–female gene combinations. In contrast, where conflict predominates, the success of an allele in one genome or in one sex diminishes the frequency of the corresponding allele in the other, resulting in considerably slower evolutionary dynamics. The rate of change of mean fitness is also much faster with positive, synergistic fitnesses and much slower where conflict is predominant. Consequently, selection rapidly fixes cooperative gene combinations, while leaving behind a slowing evolving residue of conflicting gene combinations at mutation–selection balance. We discuss how an emphasis on marginal fitness averages may obscure the interdependence of allelic fitness across genomes, making the evolutionary trajectories appear independent of one another when they are not.     

The Alternative Fitness Sets Which Preserve Allele Trajectories: A General Treatment

A general solution is presented of the problem of specifying all alternative, generally frequency-dependent, (absolute) fitness sets which give rise to the same allele frequency changes and population dynamics as a given fitness set. The one- and two-locus cases are analyzed in detail and the method is then extended to the n-locus case. It is shown that if biological constraints can be used to specify the mean fitness of the population and the relative fitnesses of the heterozygotes, then the allele frequency trajectories determine a unique fitness set.     

Genetic Divergence under Uniform Selection. II. Different Responses to Selection for Knockdown Resistance to Ethanol among DROSOPHILA MELANOGASTER Populations and Their Replicate Lines

We have tested the hypothesis that genetic differences among conspecific populations may result in diverse responses to selection, using natural populations of Drosophila melanogaster. Selection for ethanol tolerance in a tube measuring knockdown resistance was imposed on five West Coast populations. In 24 generations the selected lines increased their mean knockdown times, on average, by a factor of 2.7. An initially weak latitudinal cline was steepened by selection. The two southernmost populations showed the same increases in the selected character, but differed consistently in their correlated responses in characters related to ethanol tolerance. This result indicates that the populations responded to selection by different genetic changes. Selection decreased female body weight and increased resistance to acetone, suggesting components of the response unrelated to ethanol metabolism. The Adhs allele was favored by selection in all populations at the onset, but increased in frequency only in the selected lines of the southernmost population. There was a correlation between latitude and Adh frequency changes, suggesting that fitnesses of the Adh alleles were dependent on the genetic background. Genetic background also had a large effect on the loss of fitness due to selection. Genetic drift between replicate lines caused more variation in selection response than initial genetic differences between populations. This result demonstrates the importance of genetic drift in divergence among natural populations undergoing uniform selection, since the effective population sizes approached those of small natural populations. Drift caused greater divergence between selected replicates than control replicates. Implications of this result for the genetic model of selection response are discussed.     

Models of Frequency-Dependent Selection with Mutation from Parental Alleles

Frequency-dependent selection (FDS) remains a common heuristic explanation for the maintenance of genetic variation in natural populations. The pairwise-interaction model (PIM) is a well-studied general model of frequency-dependent selection, which assumes that a genotype’s fitness is a function of within-population intergenotypic interactions. Previous theoretical work indicated that this type of model is able to sustain large numbers of alleles at a single locus when it incorporates recurrent mutation. These studies, however, have ignored the impact of the distribution of fitness effects of new mutations on the dynamics and end results of polymorphism construction. We suggest that a natural way to model mutation would be to assume mutant fitness is related to the fitness of the parental allele, i.e., the existing allele from which the mutant arose. Here we examine the numbers and distributions of fitnesses and alleles produced by construction under the PIM with mutation from parental alleles and the impacts on such measures due to different methods of generating mutant fitnesses. We find that, in comparison with previous results, generating mutants from existing alleles lowers the average number of alleles likely to be observed in a system subject to FDS, but produces polymorphisms that are highly stable and have realistic allele-frequency distributions.     

Dynamical behavior for population genetics models of differential and difference equations with nonmonotone fitnesses

This paper studies the dynamical behavior of classical 2-dimensional models of continuously and discretely reproducing diploid populations with two alleles at one locus. The phase variables are allele frequency and population density. The genotype fitnesses are not assumed to be monotonically decreasing functions of density. Hence the mean fitness curve is more complicated than in the monotonic case. If genotype fitnesses are only density dependent, results concerning equilibrium stability are obtained similar to those for the monotonic case, and periodic solutions are precluded in the differential equation model. An example with one-hump genotype fitnesses is presented and analyzed.Research supported by funds provided by the USDA-Forest Service, Southeastern Forest Experiment Station, Pioneering (Population Genetics of Forest Trees) Research Unit, Raleigh, North Carolina     

Evolutionary principles for general frequency-dependent two-phenotype models in sexual populations

The evolutionary dynamics in general two-sex two-phenotype frequency-dependent selection models are studied with respect to underlying multi-allele one-locus genetic systems. Two classes of equilibria come into play: genotypic equilibria, with equilibrium allelic frequencies independent of the phenotype, and phenotypic equilibria, which are characterized by equal mean phenotypic fitnesses. The exact conditions for genotypic equilibria to exist and be stable and for phenotypic equilibria to exist and be evolutionarily attractive are examined. Using adequate definitions of mean fitnesses in general contexts of frequency-dependent selection in dioecious populations, we show that two phenotypes, when they can coexist in the population, tend to balance their fitnesses as far as is allowed by the genetic system as more alleles responsible for phenotype determination are introduced into the population.     

Stochastic selection in large and small populations

We describe two models of stochastic variation in selection intensity. In both models the arithmetic mean fitness of all genotypes is equal; in both models the geometric mean fitness of the heterozygous genotype is greater than that of both homozygous genotypes. In one model the correlation between the fitnesses of the homozygous genotypes is +1; in the other it is −1. We show that the expected time to absorption of an allele in a finite population is significantly retarded for all initial gene frequencies in the former model. The expected time to absorption of an allele in the latter model is retarded only at extreme initial gene frequencies; at intermediate initial gene frequencies the expected time to absorption is accelerated. We conclude that the criterion for polymorphism based on the geometric mean of the heterozygote being greater than that of both homozygotes provides only limited information about the fate of gene frequency.     

The generation and maintenance of genetic variation by frequency-dependent selection: constructing polymorphisms under the pairwise interaction model

Frequency-dependent selection remains the most commonly invoked heuristic explanation for the maintenance of genetic variation. For polymorphism to exist, new alleles must be both generated and maintained in the population. Here we use a construction approach to model frequency-dependent selection with mutation under the pairwise interaction model. The pairwise interaction model is a general model of frequency-dependent selection at the genotypic level. We find that frequency-dependent selection is able to generate a large number of alleles at a single locus. The construction process generates multiallelic polymorphisms with a wide range of allele-frequency distributions and genotypic fitness relationships. Levels of polymorphism and mean fitness are uncoupled, so constructed polymorphisms remain permanently invasible to new mutants; thus the model never settles down to an equilibrium state. Analysis of constructed fitness sets reveals signatures of heterozygote advantage and positive frequency dependence.     

THE EVOLUTION OF FITNESS

In every generation, the mean fitness of populations increases because of natural selection and decreases because of mutations and changes in the environment. The magnitudes of these effects can be measured in two ways: either directly, by comparing the fitnesses of selected and unselected populations, or indirectly, by measuring the additive variance of fitness and making use of the fundamental theorem of natural selection. The available data suggest that the amount by which natural selection increases mean fitness each generation (or degradation decreases mean fitness) will usually be between 0.1% and 30%; more tentatively, it is suggested that values will typically fall between 1% and 10%. These values can be used to set an upper limit of 5%–10% on the genetic advantage of mate choice.     

Models of general frequency-dependent selection and mating-interaction effects and the analysis of selection patterns in Drosophila inversion polymorphisms

We investigate mechanisms of balancing selection by extending two deterministic models of selection in a one-locus two-allele genetic system to allow for frequency-dependent fitnesses. Specifically we extend models of constant selection to allow for general frequency-dependent fitness functions for sex-dependent viabilities and multiplicative fertilities, while non-multiplicative mating-dependent components remain constant. We compute protected polymorphism conditions that take the form of harmonic means involving both the frequency- and the mating-dependent parameters. This allows for a direct comparison of the equilibrium properties of the frequency-dependent models with those of the constant models and for an analysis of equilibrium of the general model of constant fertility. We then apply the theory to analyze the maintenance of inversion polymorphisms in Drosophila subobscura and D. pseudoobscura, for which data on empirical fitness component estimates are available in the literature. Regression on fitness estimates obtained at different starting frequencies enables us to implement explicit fitness functions in the models and therefore to perform complete studies of equilibrium and stability for particular sets of data. The results point to frequency dependence of fitness components as the main mechanism responsible for the maintenance of the inversion polymorphisms considered, particularly in relation to heterosis, although we also discuss the contribution of other selection mechanisms.     

Frequency-dependent selection and the maintenance of genetic variation: exploring the parameter space of the multiallelic pairwise interaction model

When individuals' fitnesses depend on the genetic composition of the population in which they are found, selection is then frequency dependent. Frequency-dependent selection (FDS) is often invoked as a heuristic explanation for the maintenance of large numbers of alleles at a locus. The pairwise interaction model is a general model of FDS via intraspecific competition at the genotypic level. Here we use a parameter-space approach to investigate the full potential for the maintenance of multiallelic equilibria under the pairwise interaction model. We find that FDS maintains full polymorphism more often than classic constant-selection models and produces more skewed equilibrium allele frequencies. Fitness sets with some degree of rare advantage maintained full polymorphism most often, but a wide variety of nonobvious fitness patterns were also found to have positive potential for polymorphism. An example is put forth suggesting possible explanations for multiallelic polymorphisms maintained despite positive FDS on individual alleles.     

Convergence to constant equilibrium for a density-dependent selection model with diffusion

We consider the classical single locus two alleles selection model with diffusion where the fitnesses of the genotypes are density dependent. Using a theorem of Peter Brown, we show that in a bounded domain with homogeneous Neumann boundary conditions, the allele frequency and population density converge to a constant equilibrium lying on the zero population mean fitness curve. The results agree with the case without diffusion obtained by Selgrade and Namkoong. Frequency and density dependent selection is also considered.Research partially supported by NSF grant DMS-8601585     

Frequency-dependent viability selection (a single-locus, multi-phenotype model)

Frequency-dependent natural selection models are examined where the viability of an individual in the diploid population is determined by its phenotype and the frequency of other phenotypes present. The equilibria of the multi-phenotypic system are characterized through local mean fitness functions. It is shown that stability can best be analyzed by combining the principles of maximization of population mean fitness with the evolutionary stability conditions that apply when phenotypic fitnesses relative to the genetic constraints are equal.     

Evolution of recombination due to random drift

In finite populations subject to selection, genetic drift generates negative linkage disequilibrium, on average, even if selection acts independently (i.e., multiplicatively) upon all loci. Negative disequilibrium reduces the variance in fitness and hence, by Fisher's (1930) fundamental theorem, slows the rate of increase in mean fitness. Modifiers that increase recombination eliminate the negative disequilibria that impede selection and consequently increase in frequency by "hitchhiking." Thus, stochastic fluctuations in linkage disequilibrium in finite populations favor the evolution of increased rates of recombination, even in the absence of epistatic interactions among loci and even when disequilibrium is initially absent. The method developed within this article allows us to quantify the strength of selection acting on a modifier allele that increases recombination in a finite population. The analysis indicates that stochastically generated linkage disequilibria do select for increased recombination, a result that is confirmed by Monte Carlo simulations. Selection for a modifier that increases recombination is highest when linkage among loci is tight, when beneficial alleles rise from low to high frequency, and when the population size is small.     

Density-dependent selection migration model with non-monotone fitness functions

This paper studies the classical single locus, diallelic selection model with diffusion for a continuously reproducing population. The phase variables are population density and allele frequency (or allele density). The genotype fitness depend only on population density but include one-hump functions of the density variable. With mild assumptions on genotype fitnesses, we study the geometry of the nullclines and the asymptotic behavior of solutions of the selection model without diffusion. For the diffusion model with zero Neumann boundary conditions, we use this geometric information to show that if the initial data satisfy certain conditions then the corresponding solution to the reaction-diffusion equation converges to the spatially constant stable equilibrium which is closest to the initial data.Research partially supported by NSF grant DMS-8920597Research supported by funds provided by the USDA-Forest Service, Southeastern Forest Experiment Station, Pioneering (Population Genetics of Forest Trees) Research Unit, Raleigh, North Carolina     

Frequency dependence and cooperation: theory and a test with bacteria

Hamilton's inclusive fitness theory provides a leading explanation for the problem of cooperation. A general result from inclusive fitness theory is that, except under restrictive conditions, cooperation should not be subject to frequency-dependent selection. However, several recent studies in microbial systems have demonstrated that the relative fitness of cheaters, which do not cooperate, is greater when cheaters are rarer. Here we demonstrate theoretically that such frequency-dependent selection can occur in microbes when there is (1) sufficient population structuring or (2) an association between the level of cooperation and total population growth. We test prediction (2) and its underlying assumption, using the pathogenic bacterium Pseudomonas aeruginosa, by competing strains that produce iron-scavenging siderophore molecules (cooperators) with nonproducers (cheaters) at various ratios, under conditions that minimize population structuring. We found that both the relative fitness of cheaters and the productivity of the mixed culture were significantly negatively related to initial cheater frequency. Furthermore, when the period of population growth was experimentally shortened, the strength of frequency dependence was reduced. More generally, we argue that frequency-dependent selection on cooperative traits may be more common in microbes than in metazoans because strong selection, structuring, and cooperation-dependent growth will be more common in microbial populations.     

Adaptive Topography in Fertility-Viability Selection Models: The Haplodiploid Case

A linear combination of partial changes of mean fitnesses from one generation to the next one is shown to be approximately equal to the additive genetic variance in fitness after enough generations and away from equilibrium in random mating haplodiploid populations under arbitrary weak frequency-dependent selection on sex-differentiated viability of individuals and sex-differentiated fertility of matings controlled at a single multiallelic locus. The result can be applied to X-linked locus models in diploid populations. The result is used to deduce approximate adaptive topographies far frequency-independent selection models in the cases of nonsex-differentiated fertilities and multiplicative sex-differentiated fertilities and for kin selection models in family-structured populations under the assumptions of single insemination and multiple insemination of females. Multiple insemination creates frequency-dependent selection regimes.     

"Runaway" social evolution: reinforcing selection for inbreeding and altruism.

Kin selection theory predicts that altruistic behaviors, those that decrease the fitness of the individual performing the behavior but increase the fitness of the recipient, can increase in frequency if the individuals interacting are closely related. Several studies have shown that inbreeding therefore generally increases the effectiveness of kin selection when fitnesses are linear, additive functions of the number of altruists in the family, although with extreme forms of altruism, inbreeding can actually retard the evolution of altruism. These models assume that a constant proportion of the population mates at random and a constant proportion practices some form of inbreeding. In order to investigate the effect of inbreeding on the evolution of altruistic behavior when the mating structure is allowed to evolve, we examined a two-locus model by computer simulation of a diploid case and illustrated the important qualitative features by mathematical analysis of a haploid case. One locus determines an individual's propensity to perform altruistic social behavior and the second locus determines the probability that an individual will mate within its sibship. We assumed positive selection for altruism and no direct selection at the inbreeding locus. We observed that the altruistic allele and the inbreeding allele become positively associated, even when the initial conditions of the model assume independence between these loci. This linkage disequilibrium becomes established, because the altruistic allele increases more rapidly in the inbreeding segment of the population. This association subsequently results in indirect selection on the inbreeding locus. However, the dynamics of this model go beyond a simple "hitch-hiking" effect, because high levels of altruism lead to increased inbreeding, and high degrees of inbreeding accelerate the rate of change of the altruistic allele in the entire population. Thus, the dynamics of this model are similar to those of "runaway" sexual selection, with gene frequency change at the two loci interactively causing rapid evolutionary change.     

A method to infer positive selection from marker dynamics in an asexual population

MOTIVATION: The observation of positive selection acting on a mutant indicates that the corresponding mutation has some form of functional relevance. Determining the fitness effects of mutations thus has relevance to many interesting biological questions. One means of identifying beneficial mutations in an asexual population is to observe changes in the frequency of marked subsets of the population. We here describe a method to estimate the establishment times and fitnesses of beneficial mutations from neutral marker frequency data. RESULTS: The method accurately reproduces complex marker frequency trajectories. In simulations for which positive selection is close to 5% per generation, we obtain correlations upwards of 0.91 between correct and inferred haplotype establishment times. Where mutation selection coefficients are exponentially distributed, the inferred distribution of haplotype fitnesses is close to being correct. Applied to data from a bacterial evolution experiment, our method reproduces an observed correlation between evolvability and initial fitness defect.     

Selection, Generalized Transmission and the Evolution of Modifier Genes. I. the Reduction Principle

Modifier gene models are used to explore the evolution of features of organisms, such as the genetic system, that are not directly involved in the determination of fitness. Recent work has shown that a general "reduction principle" holds in models of selectively neutral modifiers of recombination, mutation, and migration. Here we present a framework for models of modifier genes that shows these reduction results to be part of a more general theory, for which recombination and mutation are special cases. The deterministic forces that affect the genetic composition of a population can be partitioned into two categories: selection and transmission. Selection includes differential viabilities, fertilities, and mating success. Imperfect transmission occurs as a result of such phenomena as recombination, mutation and migration, meiosis, gene conversion, and meiotic drive. Selectively neutral modifier genes affect transmission, and a neutral modifier gene can evolve only by generating association with selected genes whose transmission it affects. We show that, in randomly mating populations at equilibrium, imperfect transmission of selected genes allows a variance in their marginal fitnesses to be maintained. This variance in the marginal fitnesses of selected genes is what drives the evolution of neutral modifier genes. Populations with a variance in marginal fitnesses at equilibrium are always subject to invasion by modifier genes that bring about perfect transmission of the selected genes. It is also found, within certain constraints, that for modifier genes producing what we call "linear variation" in the transmission processes, a new modifier allele can invade a population at equilibrium if it reduces the level of imperfect transmission acting on the selected genes, and will be expelled if it increases the level of imperfect transmission. Moreover, the strength of the induced selection on the modifier gene is shown to range up to the order of the departure of the genetic system from perfect transmission.