Genetic Variation in a Heterogeneous Environment. I. Temporal Heterogeneity and the Absolute Dominance Model

The conditions for a stable polymorphism and the equilibrium gene frequency in an infinite population are compared when there is spatial or temporal environmental heterogeneity for the absolute dominance model. For temporal variation the conditions for stability are more restrictive and the equilibrium gene frequency is often at a low gene frequency. In a finite population, temporal environmental heterogeneity for the absolute dominance model was found to be quite ineffective in maintaining genetic variation and is often less effective than no selection at all. For comparison, the maximum maintenance for temporal variation is related to the overdominant model. In general, cyclic environmental variation was found to be more effective at maintaining genetic variation than where the environment varies stochastically. The importance of temporal environmental variation and the maintenance of genetic variation is discussed.     

Spatial heterogeneity in the strength of selection against deleterious alleles and the mutation load

According to current estimates of genomic deleterious mutation rates (which are often of the order 0.1-1) the mutation load (defined as a reduction in the average fitness of a population due to the presence of deleterious alleles) may be important in many populations. In this paper, I use multilocus simulations to explore the effect of spatial heterogeneity in the strength of selection against deleterious alleles on the mutation load (for example, it has been suggested that stressful environments may increase the strength of selection). These simulations show contrasted results: in some situations, spatial heterogeneity may greatly reduce the mutation load, due to the fact that migrants coming from demes under stronger selection carry relatively few deleterious alleles, and benefit from a strong advantage within demes under weaker selection (where individuals carry many more deleterious alleles); in other situations, however, deleterious alleles accumulate within demes under stronger selection, due to migration pressure from demes under weaker selection, leading to fitness erosion within those demes. This second situation is more frequent when the productivity of the different demes is proportional to their mean fitness. The effect of spatial heterogeneity is greatly reduced, however, when the response to environmental differences is inconsistent across loci.     

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.     

Population bottleneck and effective size in Bonamia ostreae-resistant populations of Ostrea edulis as inferred by microsatellite markers

Genetic variability at five microsatellite loci was analysed in three hatchery-propagated populations of the flat oyster, Ostrea edulis. These populations were part of a selection programme for resistance to the protozoan parasite Bonamia ostreae and were produced by mass spawns, without control of the genealogy. Evidence for population bottlenecks and inbreeding was sought. A reduction in the number of alleles, mainly due to the loss of rare alleles, was observed in all selected populations, relative to the natural population from which they were derived. Heterozygote excesses were observed in two populations, and were attributed to substructuring of the population into a small number of families. Pedigree reconstruction showed that these two populations were produced by at most two spawning events involving a limited number of parents. Most individuals within these populations are half or full-sib, as shown by relatedness coefficients. The occurrence of population bottlenecks was supported by estimates of effective number of breeders derived by three methods: temporal variance in allelic frequencies, heterozygote excess, and a new method based on reduction in the number of alleles. The estimates from the different methods were consistent. The evidence for bottleneck and small effective number of breeders are expected to lead to increasing inbreeding, and have important consequences for the future management of the three O. edulis selected populations.     

Comparison of selection methods for improvement of the population hybrid

Summary The objective of this study was to compare several selection procedures with respect to expected genetic gain in the population hybrid across a range of initial allelic frequencies, degrees of dominance, and environmental variances. The methods compared were intrapopulation recurrent selection using full-sib or S₁ families, full-sib and two half-sib reciprocal recurrent selection procedures, and convergent improvement applied to populations. Comparisons were made by calculating expected allelic frequency changes for each method. The optimal selection method for a given set of allelic frequencies and degree of dominance depended little on the environmental variance. Partly because of its short cycle, full-sib intrapopulation selection was the most effective method for the majority of allelic frequency combinations when the degree of dominance was small and an off-season nursery could be used to make recombinations. With larger values for the degree of dominance, S₁ and reciprocal full-sib methods became optimal, the former method especially when favorable alleles had a high frequency and the latter when populations were highly divergent. When off-season nursery use was restricted to making self-pollinations or was absent, S₁ selection was optimal for the majority of allelic frequency combinations. Convergent improvement was superior only for extremely divergent allelic frequencies and then only when the degree of dominance was less than 0.10. Half-sib reciprocal methods were never optimal, although the gain for the standard half-sib reciprocal procedure differed little from that of full-sib reciprocal selection when the degree of dominance was 0.75.     

Evolutionary Dynamics of Sporophytic Self-Incompatibility Alleles in Plants

The stationary frequency distribution and allelic dynamics in finite populations are analyzed through stochastic simulations in three models of single-locus, multi-allelic sporophytic self-incompatibility. The models differ in the dominance relationships among alleles. In one model, alleles act codominantly in both pollen and style (SSIcod), in the second, alleles form a dominance hierarchy in pollen and style (SSIdom). In the third model, alleles interact codominantly in the style and form a dominance hierarchy in the pollen (SSIdomcod). The SSIcod model behaves similarly to the model of gametophytic self-incompatibility, but the selection intensity is stronger. With dominance, dominant alleles invade the population more easily than recessive alleles and have a lower frequency at equilibrium. In the SSIdom model, recessive alleles have both a higher allele frequency and higher expected life span. In the SSIdomcod model, however, loss due to drift occurs more easily for pollen-recessive than for pollen-dominant alleles, and therefore, dominant alleles have a higher expected life span than the more recessive alleles. The process of allelic turnover in the SSIdomcod and SSIdom models is closely approximated by a random walk on a dominance ladder. Implications of the results for experimental studies of sporophytic self-incompatibility in natural populations are discussed.     

Molecular mechanisms of dominance evolution in Müllerian mimicry

Natural selection acting on dominance between adaptive alleles at polymorphic loci can be sufficiently strong for dominance to evolve. However, the molecular mechanisms underlying such evolution are generally unknown. Here, using Müllerian mimicry as a case‐study for adaptive morphological variation, we present a theoretical analysis of the invasion of dominance modifiers altering gene expression through different molecular mechanisms. Toxic species involved in Müllerian mimicry exhibit warning coloration, and converge morphologically with other toxic species of the local community, due to positive frequency‐dependent selection acting on these colorations. Polymorphism in warning coloration may be maintained by migration–selection balance with fine scale spatial heterogeneity. We modeled a dominance modifier locus altering the expression of the warning coloration locus, targeting one or several alleles, acting in cis or trans, and either enhancing or repressing expression. We confirmed that dominance could evolve when balanced polymorphism was maintained at the color locus. Dominance evolution could result from modifiers enhancing one allele specifically, irrespective of their linkage with the targeted locus. Nonspecific enhancers could also persist in populations, at frequencies tightly depending on their linkage with the targeted locus. Altogether, our results identify which mechanisms of expression alteration could lead to dominance evolution in polymorphic mimicry.     

COMPUTER PROGRAMS: nessi: a program for numerical estimations in sporophytic self-incompatibility genetic systems

nessi is a computer program generating predictions about allelic and genotypic frequencies at the S-locus in sporophytic self-incompatibility systems under finite and infinite populations. For any pattern of dominance relationships among self-incompatibility alleles, nessi computes deterministic equilibrium frequencies and estimates distributions in samples from finite populations of the number of alleles at equilibrium, allelic and genotypic frequencies at equilibrium and allelic and genotypic frequency changes in a single generation. These predictions can be used to rigorously test the impact of negative frequency-dependent selection on diversity patterns in natural populations.     

Evolution of mutation rates in bacteria

Evolutionary success of bacteria relies on the constant fine-tuning of their mutation rates, which optimizes their adaptability to constantly changing environmental conditions. When adaptation is limited by the mutation supply rate, under some conditions, natural selection favours increased mutation rates by acting on allelic variation of the genetic systems that control fidelity of DNA replication and repair. Mutator alleles are carried to high frequency through hitchhiking with the adaptive mutations they generate. However, when fitness gain no longer counterbalances the fitness loss due to continuous generation of deleterious mutations, natural selection favours reduction of mutation rates. Selection and counter-selection of high mutation rates depends on many factors: the number of mutations required for adaptation, the strength of mutator alleles, bacterial population size, competition with other strains, migration, and spatial and temporal environmental heterogeneity. Such modulations of mutation rates may also play a role in the evolution of antibiotic resistance.     

Distribution of Gene Frequency as a Test of the Theory of the Selective Neutrality of Polymorphisms

The variation in gene frequency among populations or between generations within a population is a result of breeding structure and selection. But breeding structure should affect all loci and alleles in the same way. If there is significant heterogeneity between loci in their apparent inbreeding coefficients F=s_p²/p(1-p), this heterogeneity may be taken as evidence for selection. We have given the statistical properties of F and shown how tests of heterogeneity can be made. Using data from human populations we have shown highly significant heterogeneity in F values for human polymorphic genes over the world, thus demonstrating that a significant fraction of human polymorphisms owe their current gene frequencies to the action of natural selection. We have also applied the method to temporal variation within a population for data on Dacus oleae and have found no significant evidence of selection.     

Linking genetic mechanisms of heterozygosity-fitness correlations to footprints of selection at single loci

Investigations of heterozygosity-fitness correlations (HFCs) are central to the understanding how genetic diversity is maintained in natural populations. Advanced genome-wide approaches will enrich the number of functional loci to be tested. We argue that a combined analysis of the genetic mechanisms of HFCs and selection signals at single loci will allow researchers to better understand the micro-evolutionary basis of HFCs. Different dominance relationships among the alleles at the locus can lead to positive, negative or null HFCs depending on the allele frequency distribution. These scenarios differ in the temporal stability of the HFCs and in the patterns of allele frequency changes over time. Here, we describe a simple theoretical framework that links the analyses of heterozygosity-fitness associations (ecological timescale) with tests for selection signals (evolutionary timescale). Different genomic footprints of selection can be expected for the different underlying genetic mechanisms of HFCs, and this information can be independently used for the classification of HFCs. We suggest that in addition to inbreeding and single-locus overdominant effects also loci under directional selection could play a significant role in the development of heterozygosity-fitness effects in large natural populations under recent or fluctuating ecological changes.     

Selection in autotetraploids

Summary Theoretical studies indicated that response to selection would always be greater in diploid than in autotetraploid populations when gene frequency was the same in both, and that situations in which little or no response to selection could be expected would be more frequent in autotetraploids. Interpretation of the coefficient of selection in terms of escape from infection in a program of selection for disease or insect resistance indicated that moderate levels of escape from infection can drastically reduce response to selection in some cases.The zygotic constitution of an autotetraploid population will change as it approaches a new random mating equilibrium once selection pressure is relaxed. The changes will result in no change in the population mean if the trait under selection exhibits no dominance, but the mean will decrease slightly if there is dominance.     

Selection, load and inbreeding depression in a large metapopulation

The subdivision of a species into local populations causes its response to selection to change, even if selection is uniform across space. Population structure increases the frequency of homozygotes and therefore makes selection on homozygous effects more effective. However, population subdivision can increase the probability of competition among relatives, which may reduce the efficacy of selection. As a result, the response to selection can be either increased or decreased in a subdivided population relative to an undivided one, depending on the dominance coefficient F(ST) and whether selection is hard or soft. Realistic levels of population structure tend to reduce the mean frequency of deleterious alleles. The mutation load tends to be decreased in a subdivided population for recessive alleles, as does the expected inbreeding depression. The magnitude of the effects of population subdivision tends to be greatest in species with hard selection rather than soft selection. Population structure can play an important role in determining the mean fitness of populations at equilibrium between mutation and selection.     

Kin selection and evolution of infectious disease resistance

Discoveries of mutations conferring resistance to infectious diseases have led to increased interest in the evolutionary dynamics of disease resistance. Several recent papers have estimated the historical strength of selection for mutations conferring disease resistance. These studies are based on simple population genetic models that do not take account of factors such as spatial and family structure. Such factors may have a substantial impact on the strength of natural selection through inclusive fitness effects. That is, people have a strong tendency to live with relatives and therefore have a high probability of transmitting infectious diseases to them. Thus, an allele that protects an individual against disease infection also protects that individual's family members. Because some of these family members are likely to also be carrying the allele, selection for that allele is magnified by family structure. In this paper, I use mathematical modeling techniques to explore the impact of such kin selection on the strength of selection for infectious disease resistance alleles. I show that if the resistance allele has the same proportional effect on both within- and between-family transmission, then the impact of kin selection is relatively minor. Selection coefficients are increased by 5-35%, with a greater benefit for weaker alleles. The reason is that an individual with a strong resistance allele does not need much protection from infection by family members and thus does not benefit much from their alleles. The effect of kin selection can be dramatic, however, if the resistance allele has a larger effect on between-family transmission than within-family transmission (which can occur if between-family infection rates are much smaller than within-family rates), increasing selection coefficients by as much as two- to threefold. These results show conditions when it is important to consider family structure in estimates of the strength of selection for infectious disease resistance alleles.     

Psorophora columbiae and Psorophora toltecum (Diptera: Culicidae) Colombian populations cannot be differentiated by isoenzymes

Two populations of the mosquito Psorophora columbiae from the central Andean area of Colombia and one population of Ps. toltecum from the Atlantic coast of Colombia were analyzed for 11 isoenzyme markers. Psorophora columbiae and Ps. toltecum are two of the main vectors of Venezuelan equine encephalitis. We found no conspicuous genetic differences between the two species. The relatively high gene flow levels among these populations indicate that these are not two different species or that there has been recent divergence between these taxa. In addition, no global differential selection among the loci was detected, although the alpha-GDH locus showed significantly less genetic heterogeneity than the remaining loci, which could mean that homogenizing natural selection acts at this locus. No isolation by distance was detected among the populations, and a spatial population analysis showed opposite spatial trends among the 31 alleles analyzed. Multiregression analyses showed that both expected heterozygosity and the average number of alleles per locus were totally determined by three variables: altitude, temperature and size of the human population at the locality. Individually, the expected heterozygosity is more related to these three variables than to the average number of alleles.     

Spatially and temporally varying selection on intrapopulation quantitative trait loci for a life history trade-off in Mimulus guttatus

Why do populations remain genetically variable despite strong continuous natural selection? Mutation reconstitutes variation eliminated by selection and genetic drift, but theoretical and experimental studies each suggest that mutation‐selection balance insufficient to explain extant genetic variation in most complex traits. The alternative hypothesis of balancing selection, wherein selection maintains genetic variation, is an aggregate of multiple mechanisms (spatial and temporal heterogeneity in selection, frequency‐dependent selection, antagonistic pleiotropy, etc.). Most of these mechanisms have been demonstrated for Mendelian traits, but there is little comparable data for loci affecting quantitative characters. Here, we report a 3‐year field study of selection on intrapopulation quantitative trait loci (QTL) of flower size, a highly polygenic trait in Mimulus guttatus. The QTL exhibit antagonistic pleiotropy: alleles that increase flower size, reduce viability, but increase fecundity. The magnitude and direction of selection fluctuates yearly and on a spatial scale of metres. This study provides direct evidence of balancing selection mechanisms on QTL of an ecologically relevant trait.     

Patterns of differentiation in a colour polymorphism and in neutral markers reveal rapid genetic changes in natural damselfly populations

The existence and mode of selection operating on heritable adaptive traits can be inferred by comparing population differentiation in neutral genetic variation between populations (often using F(ST) values) with the corresponding estimates for adaptive traits. Such comparisons indicate if selection acts in a diversifying way between populations, in which case differentiation in selected traits is expected to exceed differentiation in neutral markers [F(ST )(selected) > F(ST )(neutral)], or if negative frequency-dependent selection maintains genetic polymorphisms and pulls populations towards a common stable equilibrium [F(ST) (selected) < F(ST) (neutral)]. Here, we compared F(ST) values for putatively neutral data (obtained using amplified fragment length polymorphism) with estimates of differentiation in morph frequencies in the colour-polymorphic damselfly Ischnura elegans. We found that in the first year (2000), population differentiation in morph frequencies was significantly greater than differentiation in neutral loci, while in 2002 (only 2 years and 2 generations later), population differentiation in morph frequencies had decreased to a level significantly lower than differentiation in neutral loci. Genetic drift as an explanation for population differentiation in morph frequencies could thus be rejected in both years. These results indicate that the type and/or strength of selection on morph frequencies in this system can change substantially between years. We suggest that an approach to a common equilibrium morph frequency across all populations, driven by negative frequency-dependent selection, is the cause of these temporal changes. We conclude that inferences about selection obtained by comparing F(ST) values from neutral and adaptive genetic variation are most useful when spatial and temporal data are available from several populations and time points and when such information is combined with other ecological sources of data.     

Extinction time and age of an allele in a large finite population

For a single locus with two alleles we study the expected extinction and fixation times of the alleles under the influence of selection and genetic drift. Using a diffusion model we derive asymptotic approximations for these expected exit times for large populations. We consider the case where the fitness of the heterozygote is in between the fitnesses of the homozygotes (incomplete dominance). The asymptotic analysis not only yields new quantitative results but also reveals interesting features that remain hidden in the exact solution. Some of the outcomes are extensions of results known in the literature. The asymptotic approximations also apply to the expected first arrival time of an allele at a specified frequency and to the expected age of an allele.     

The spatial patterns of directional phenotypic selection

Local adaptation, adaptive population divergence and speciation are often expected to result from populations evolving in response to spatial variation in selection. Yet, we lack a comprehensive understanding of the major features that characterise the spatial patterns of selection, namely the extent of variation among populations in the strength and direction of selection. Here, we analyse a data set of spatially replicated studies of directional phenotypic selection from natural populations. The data set includes 60 studies, consisting of 3937 estimates of selection across an average of five populations. We performed meta‐analyses to explore features characterising spatial variation in directional selection. We found that selection tends to vary mainly in strength and less in direction among populations. Although differences in the direction of selection occur among populations they do so where selection is often weakest, which may limit the potential for ongoing adaptive population divergence. Overall, we also found that spatial variation in selection appears comparable to temporal (annual) variation in selection within populations; however, several deficiencies in available data currently complicate this comparison. We discuss future research needs to further advance our understanding of spatial variation in selection.     

Combined genetic analysis of partial blast resistance in an upland rice population and recurrent selection for line and hybrid values

The CNA-IRAT 5 upland rice population has been improved for 4 years by recurrent selection for blast resistance in Brazil. In order to predict the efficiency of recurrent selection in different test systems and to compare the relative advantage of hybrids versus pure line breeding, a combined genetic analysis of partial blast resistance in the CNA-IRAT 5 population was undertaken. A three-level hierarchical design in inbreeding and a factorial design were derived from the base population. Partial blast resistance of lines and hybrids was evaluated in the greenhouse and in the field by inoculation with one virulent blast isolate. The means and genetic variances of the hybrids and lines were estimated. Genetic advance by recurrent selection was predicted from estimates of variance components. The inheritance of partial blast resistance was mainly additive but non-additive effects were detected at both levels of means and variances. Mean heterosis ranged from 4%–8% for lesion size and lesion density to 10–12% for leaf and panicle resistance. High dominance or homozygous dominance variances relative to additive variance and negative covariance between additive and homozygous dominance effects were estimated. A low frequency of favourable alleles for partial resistance would explain the observed organisation of genetic variability in the base population. Recurrent selection will efficiently improve partial blast resistance of the CNA-IRAT 5 population. Genetic advance for line or hybrid values was expected to be higher testing doubled haploid lines than S1 lines, or than general combining ability. Two components of partial resistance assessed in the greenhouse, lesion size and lesion density, could be used as indirect selection criteria to improve field resistance. On the whole, hybrid breeding for partial blast resistance appeared to be slightly more advantageous than pure line breeding.