Polymorphism with selection and genotype-dependent migration

For a single autosomal locus with multiple alleles both an island and a multiple-niche model with discrete nonoverlapping generations are formulated for the maintenance of genetic variability. Both models incorporate viability selection in an arbitrary way and allow for genotypic differences in the pertinent migration structure. Random drift is ignored, and mating is at random. A global analysis is given for the island model in the neutral case. For a subdivided population, conditions are derived for the existence of a protected polymorphism, and the model is examined in some special two-niche cases. Of particular consideration is the loss of neutral alleles due solely to population regulation and genotype-dependent migration, and the possible existence of equilibrium clines without selection.M. M. was supported by USPHS Pre-doctoral training grant No. GM 7197 to the University of Chicago; this work represents part of the author's Doctoral dissertation.     

Migration-selection polymorphism in dioecious populations

Summary Conditions are derived for a protected polymorphism in a dioecious population subdivided into an arbitrary number of demes which exchange migrants. Generations are discrete and nonoverlapping; mutation and random drift are neglected. The analysis is restricted to a diallelic autosomal locus. In contrast to the monoecious case, the protection criteria depend on the order of migration and selection; they become identical for adult and juvenile migration if both the male and female backward migration matrices are symmetric, or the migration or selection patterns in the two sexes are the same. The protection conditions are presented explicitly for the Levene model. A recessive allele is protected in a panmictic dioecious population if the unweighted average of the recessive-to-dominant fitness ratios in the two sexes exceeds unity.Supported by the National Science Foundation (Grant No. DEB77-21494)     

Population trajectories for single locus additive fecundity selection and related selection models

A selection model which comprises models of additive fecundities as well as models of viability, fecundity, or differential mating selection acting only in one sex, is investigated for an autosomal gene locus in a population reproducing in nonoverlapping generations. The recurrence equations and basic properties of the genotypic population trajectories and equilibrium points are formulated for the multiallelic case. For the diallelic case, the trajectory development is discussed in more detail, and it is proven that every population trajectory converges to a Hardy-Weinberg equilibrium point.     

The effect on neutral gene flow of selection at a linked locus

The effect on gene flow at a neutral locus of a selective cline at a linked locus is investigated. A diffusion approximation for a two-locus island model is derived in which only one locus is subject to selection. The moments of the stationary distribution are obtained and compared to the corresponding moments from a one-locus, neutral island model. This comparison yields an effective migration rate. The effective migration rate is always less than the actual migration rate, but this effect is seen to be small for weak selection and loose linkage in the case of adult migration. The importance of selection at linked loci to the question of genetic differentiation in a subdivided population is discussed.     

The diffusion model for migration and selection in a plant population

 The diffusion approximation is derived for migration and selection at a multiallelic locus in a partially selfing plant population subdivided into a lattice of colonies. Generations are discrete and nonoverlapping; both pollen and seeds disperse. In the diffusion limit, the genotypic frequencies at each point are those determined at equilibrium by the local rate of selfing and allelic frequencies. If the drift and diffusion coefficients are taken as the appropriate linear combination of the corresponding coefficients for pollen and seeds, then the migration terms in the partial differential equation for the allelic frequencies have the standard form for a monoecious animal population. The selection term describes selection on the local genotypic frequencies. The boundary conditions and the unidimensional transition conditions for a geographical barrier and for coincident discontinuities in the carrying capacity and migration rate have the standard form. In the diallelic case, reparametrization renders the entire theory of clines and of the wave of advance of favorable alleles directly applicable to plant populations. Received 30 August 1995; received in revised form 23 February 1996     

A diffusion model for geographically structured populations

Summary A diffusion model is derived for the evolution of a diploid monoecious population under the influence of migration, mutation, selection, and random genetic drift. The population occupies an unbounded linear habitat; migration is independent of genotype, symmetric, and homogeneous. The treatment is restricted to a single diallelic locus without dominance. With the customary diffusion hypotheses for migration and the assumption that the mutation rates, selection coefficient, variance of the migrational displacement, and reciprocal of the population density are all small and of the same order of magnitude, a boundary value problem is deduced for the mean gene frequency and the covariance between the gene frequencies at any two points in the habitat. Supported by the National Science Foundation (Grant No. DEB77-21494).     

Geographical Variation in a Quantitative Character

A model for the evolution of the local averages of a quantitative character under migration, selection, and random genetic drift in a subdivided population is formulated and investigated. Generations are discrete and nonoverlapping; the monoecious, diploid population mates at random in each deme. All three evolutionary forces are weak, but the migration pattern and the local population numbers are otherwise arbitrary. The character is determined by purely additive gene action and a stochastically independent environment; its distribution is Gaussian with a constant variance; and it is under Gaussian stabilizing selection with the same parameters in every deme. Linkage disequilibrium is neglected. Most of the results concern the covariances of the local averages. For a finite number of demes, explicit formulas are derived for (i) the asymptotic rate and pattern of convergence to equilibrium, (ii) the variance of a suitably weighted average of the local averages, and (iii) the equilibrium covariances when selection and random drift are much weaker than migration. Essentially complete analyses of equilibrium and convergence are presented for random outbreeding and site homing, the Levene and island models, the circular habitat and the unbounded linear stepping-stone model in the diffusion approximation, and the exact unbounded stepping-stone model in one and two dimensions.     

The Maintenance of Single-Locus Polymorphism. I. Numerical Studies of a Viability Selection Model

The ability of viability selection to maintain single-locus polymorphism is investigated with two models in which the population is bombarded with a series of mutations with random fitnesses. In the first model, the population is allowed to reach equilibrium before mutation resumes; in the second the iterations and mutation occur simultaneously. Monte Carlo simulations of these models show that viability selection is easily able to maintain stable 6- or 7-allele polymorphisms and that monomorphisms and diallelic polymorphisms are uncommon. The question of how monomorphisms arise is also discussed.     

The island model with stochastic migration

The island model with stochastically variable migration rate and immigrant gene frequency is investigated. It is supposed that the migration rate and the immigrant gene frequency are independent of each other in each generation, and each of them is independently and identically distributed in every generation. The treatment is confined to a single diallelic locus without mutation. If the diploid population is infinite, selection is absent and the immigrant gene frequency is fixed, then the gene frequency on the island converges to the immigrant frequency, and the logarithm of the absolute value of its deviation from it is asymptotically normally distributed. Assuming only neutrality, the evolution of the exact mean and variance of the gene frequency are derived for an island with finite population. Selection is included in the diffusion approximation: if all evolutionary forces have comparable roles, the gene frequency will be normally distributed at all times. All results in the paper are completely explicit.     

The robustness of neutral models of geographical variation

The approximation of diploid migration by gametic dispersion is studied. The monoecious, diploid population is subdivided into panmictic colonies that exchange migrants. Generations are discrete and nonoverlapping; the analysis is restricted to a single locus in the absence of selection; every allele mutates to a new allele at the same rate u. Diploid-migration models without self-fertilization and with selfing at the “random” rate (equal to the reciprocal of the deme size in each deme) are investigated; in the gametic-dispersion models, selfing occurs at the random rate. It is shown for the unbounded stepping-stone model in one and two dimensions, the circular stepping-stone model, and the island model that the probabilitities of identity in state at equilibrium for diploid migration are close to those for gametic dispersion if the mutation rate is small or the deme size is large. Explicit error bounds are presented in all the above cases. It is also proved that if the number of demes is finite and the migration matrix is arbitrary but time independent and ergodic, then in the strong-migration approximation the equilibrium and the ultimate rate and pattern of convergence of both diploid-dispersion models are close to the corresponding gametic-dispersion formulae. For the strong-migration approximation at equilibrium, migration must dominate both mutation and random drift; for the convergence results, it suffices that migration dominate random drift. All the results apply to a dioecious population if the migration pattern and mutation rate are sex independent.     

Multinomial-Sampling Models for Random Genetic Drift

Three different derivations of models with multinomial sampling of genotypes in a finite population are presented. The three derivations correspond to the operation of random drift through population regulation, conditioning on the total number of progeny, and culling, respectively. Generations are discrete and nonoverlapping; the diploid population mates at random. Each derivation applies to a single multiallelic locus in a monoecious or dioecious population; in the latter case, the locus may be autosomal or X-linked. Mutation and viability selection are arbitrary; there are no fertility differences. In a monoecious population, the model yields the Wright-Fisher model (i.e., multinomial sampling of genes) if and only if the viabilities are multiplicative. In a dioecious population, the analogous reduction does not occur even for pure random drift. Thus, multinomial sampling of genotypes generally does not lead to multinomial sampling of genes. Although the Wright-Fisher model probably lacks a sound biological basis and may be inaccurate for small populations, it is usually (perhaps always) a good approximation for genotypic multinomial sampling in large populations.     

Dynamic Regimes in a Model of Single-Locus Density-Dependent Selection

A model of density-dependent selection in a Mendelian single-locus population was analyzed in the case where the fitnesses of genotypic forms are exponential functions of the population size. Analytical and numerical studies of the model were performed for a diallelic locus, and parametric regions were established for different dynamic behaviors of the model. The diallelic model of density-dependent selection was generalized to a multiallelic locus; the results of its analysis are described.     

Dynamic regimes in a model of single-locus density-dependent selection

A model of density-dependent selection in a Mendelian single-locus population was analyzed in the case where the fitnesses of genotypic forms are exponential functions of the population size. Analytical and numerical studies of the model were performed for a diallelic locus, and parametric regions were established for different dynamic behaviors of the model. The diallelic model of density-dependent selection was generalized to a multiallelic locus; the results of its analysis are described.     

The diffusion model for migration and selection in a dioecious population

The diffusion approximation is derived for migration and selection at a multiallelic locus in a dioecious population subdivided into a lattice of panmictic colonies. Generations are discrete and nonoverlapping; autosomal and X-linked loci are analyzed. The relation between juvenile and adult subpopulation numbers is very general and includes both soft and hard selection; the zygotic sex ratio is the same in every colony. All the results hold for both adult and juvenile migration. If ploidy-weighted average selection, drift, and diffusion coefficients are used, then the ploidy-weighted average allelic frequencies satisfy the corresponding partial differential equation for a monoecious population. The boundary conditions and the unidimensional transition conditions for coincident discontinuities in the carrying capacity and migration rate extend identically. The previous unidimensional formulation and analysis of symmetric, nearest-neighbor migration of a monoecious population across a geographical barrier is generalized to symmetric migration of arbitrary finite range, and the transition conditions are shown to hold for a dioecious population. Thus, the entire theory of clines and of the wave of advance of favorable alleles is applicable to dioecious populations.This work was supported by National Science Foundation grant BSR-9006285     

Further properties of Gavrilets' one-locus two-allele model of maternal selection

I derive several properties of the model proposed by Gavrilets for maternal selection at a single diallelic locus. Most notably, (i) stable oscillations of genotype frequencies (i.e., cycling) can occur and (ii) in the special case in which maternal effects and standard viability selection act multiplicatively, maternal selection effectively acts on maternally derived alleles only.     

Mutation-Selection Balance under Genomic Imprinting at an Autosomal Locus

I model the effect of genomic imprinting on the equilibrium allele frequencies at an autosomal diallelic locus subject to viability selection and mutation. The population size is assumed to be very large; male and female mutation rates may be unequal. Different models examine cases of the inactivation of one gene (with both complete and partial penetrance) and of differential expression of genes according to the parent of origin. In the simplest cases the frequency of the deleterious allele is approximately twice that of a dominant nonimprinting mutant, but considerably less than that of a recessive nonimprinting mutant. Under imprinting, selection and unequal mutation rates interact: other things being equal, male-biased mutation leads to lower mutant frequencies under maternal imprinting and higher frequencies under paternal imprinting. I also model cases where just one allele is imprintable (and the other not). These models allow us to predict the frequency of a failure to imprint in a normally imprinting system, as well as the frequency of imprinting at a standard nonimprinting locus.     

Conditions for the Existence of Clines

A very general partial differential equation in space and time satisfied by the gene frequency in a monoecious population distributed continuously over an arbitrary habitat is derived. The treatment is restricted to a single diallelic locus in the absence of mutation and random drift, and it is supposed that time is continuous, births and deaths occur at random, and migration is independent of genotype. With the further assumptions that migration is isotropic and homogeneous, the population density is constant and uniform (as permitted by the population regulation mechanism included in the formulation), and Hardy-Weinberg proportions obtain locally, this partial differential equation reduces to the simplest multidimensional generalization of the classical Fisher-Haldane cline model. The efficacy of migration and selection in maintaining genetic variability at equilibrium in this model is investigated by deducing conditions for the existence of clines under various circumstances. The effects of the degree of dominance, a neutral belt between the regions where a particular allele is advantageous and deleterious, finiteness of the habitat, and habitat dimensionality are evaluated. Provided at least one of the alleles is favored only in a finite region, excluding the special case in which its total effective selective coefficient is zero, if conditions for supporting a cline are too unfavorable because migration is too strong, selection is too weak, or both, a cline cannot exist at all. Thus, unless there is overdominance, the population must be monomorphic. It is possible for a cline which can barely exist under the prevailing ecological circumstances to show a large amount of variation in gene frequency.     

Effective size and F-statistics of subdivided populations for sex-linked loci.

For a population subdivided into an arbitrary number (s) of subpopulations, each consisting of different numbers of separate sexes, with arbitrary distributions of family size and variable migration rates by males (dm) and females (df), the recurrence equations for inbreeding coefficient and coancestry between individuals within and among subpopulations for a sex-linked locus are derived and the corresponding expressions for asymptotic effective size are obtained by solving the recurrence equations. The usual assumptions are made which are stable population size and structure, discrete generations, the island migration model, and without mutation and selection. The results show that population structure has an important effect on the inbreeding coefficients in any generation, asymptotic effective size, and F-statistics. Gene exchange among subpopulations inhibits inbreeding in initial generations but increases inbreeding in later generations. The larger the migration rate, the greater the final inbreeding coefficients and the smaller the effective size. Thus if the inbreeding coefficient is to be restricted to a specific value within a given number of generations, the appropriate population structure (the values of s, dm, and df) can be obtained by using the recurrence equations. It is shown that the greater the extent of subdivision (large s, small dm and df), the larger the effective size. For a given subdivided population, the effective size for a sex-linked locus may be larger or smaller than that for an autosomal locus, depending on the sex ratio, variance and covariance of family size, and the extend of subdivision. For the special case of a single unsubdivided population, our recurrence equations for inbreeding coefficient and coancestry and formulas for effective size reduce to the simple expressions derived by previous authors.     

On the opportunity for polymorphism with sex-linkage or haplodiploidy

This paper addresses the assertion that X-linked and haplodiploid genetic systems are inherently limited with respect to the potential for selectively maintained genetic polymorphisms. Using a variation of Haldane and Jayakar's (1964) parameterization of selection on an X-linked locus, analytical expressions are derived for the proportion of the total parameter space (P) in which stable diallelic polymorphism is attained. P is a function of the ratio of selection coefficients (r) associated with homozygous and hemizygous genotypes, and the intensity of selection (s). Analytical expressions for the opportunity for polymorphism at an autosomal locus (P(a)) are also derived for comparison to the X-linked case. P and P(a) are maximal and equal if the ratios of selection coefficients are -1 and selection is intense. Otherwise, P is slightly less than P(a), but the difference between autosomal and sex-linked loci is less than the range of values of P obtained over the range of r. Several arguments are presented suggesting that polymorphism arising from differential selection in the sexes (r < 0) is probabilistically and biologically feasible.     

Evolution under multiallelic migration-selection models

The loss of a specified allele and the convergence of the gene frequencies at a single multiallelic locus under the joint action of migration and viability selection are investigated. The monoecious, diploid population is subdivided into finitely many panmictic colonies that exchange adult migrants independently of genotype. Sufficient conditions are established for global fixation and for global loss of a particular allele. When migration is either sufficiently weak or sufficiently strong relative to selection, the equilibria are described, convergence of the gene frequencies is demonstrated, and sufficient conditions for the increase of a suitably defined mean fitness are offered. If the selection pattern is the same in every colony and such that in a panmictic population there is a globally asymptotically stable, internal (i.e., completely polymorphic) equilibrium point, then under certain weak assumptions on migration, the gene frequencies in the subdivided population converge globally to that equilibrium point. Thus, in this case, the ultimate state of the population is unaffected by geographical structure.