Equilibria and dynamics of selection at a diallelic autosomal locus in the Nagylaki-Crow continuous model of a monoecious population

Let birth rates and death rates be constant, birth rates positive, fertilities additive, and each birth rate not larger than twice any other birth rate. Global convergence to equilibria is proved for the model in the title. There is at most one polymorphic equilibrium or there are a continuum of equilibria. The phase portraits are given. If there is a polymorphic equilibrium, then the largest negatively invariant set in the state space is a continuous curve connecting the two fixation equilibria. This curve coincides with the Hardy-Weinberg manifold exactly when the death rate is additive. Disregarding extinction, the polymorphic equilibria are the same for the continuous model as for the corresponding discrete model exactly when the death rate is additive.     

Equilibria and dynamics of selection at a diallelic autosomal locus in the Nagylaki-Crow continuous model of a monoecious population with random mating

Global convergence to equilibria is proved for the model in the title under the assumption that death rates (and fertilities) are constant and fertilities are additive and positive. The phase portraits are determined. There is at most one polymorphic equilibrium or there are a continuum of equilibria.     

Evolution under Fertility and Viability Selection

Evolution at a single multiallelic locus under arbitrary weak selection on both fertility and viability is investigated. Discrete, nonoverlapping generations are posited for autosomal and X-linked loci in dioecious populations, but monoecious populations are studied in both discrete and continuous time. Mating is random. The results hold after several generations have elapsed. With an error of order s [i.e., O(s)], where s represents the selection intensity, the population evolves in Hardy-Weinberg proportions. Provided the change per generation of the fertilities and viabilities due to their explicit time dependence (if any) is O(s2), the rate of change of the deviation from Hardy-Weinberg proportions is O(s2). If the change per generation of the viabilities and genotypic fertilities is smaller than second order [i.e., o(s2)], then to O(s2) the rate of change of the mean fitness is equal to the genic variance. The mean fitness is the product of the mean fertility and the mean viability; in dioecious populations, the latter is the unweighted geometric mean of the mean viabilities of the two sexes. Hence, as long as there is significant gene frequency change, the mean fitness increases. If it is the fertilities of matings that change slowly [at rate o(s2)], the above conclusions apply to a modified mean fitness, defined as the product of the mean viability and the square root of the mean fertility.     

The multiple-locus symmetric fertility model

Selection due to variation in the fecundity among matings of genotypes with respect to many loci each with two alleles is studied. The fitness of a mating depends only on the genotypic distinction between homozygote and heterozygote at each locus in the two individuals, and differences among loci are allowed. This symmetric fertility model is therefore a generalization of the multiple-locus symmetric viability model. The phenomena seen in the two-locus symmetric fertility model generalize—e.g., the possibility of joint stability of equilibria with linkage equilibrium and with linkage disequilibrium, and the existence of different types of totally polymorphic equilibria with the gametic proportions in linkage equilibrium. The central equilibrium with genotypic frequencies in Hardy-Weinberg proportions and gametic frequencies in Robbins proportions exists for all symmetric fertility models. For some symmetric fertility regimes additional equilibria exist with gametic frequencies in linkage equilibrium and with genotypic frequencies in Hardy-Weinberg proportions at all except one locus. These equilibria may exist in the dioecious symmetric viability model, and then they will be locally stable. For free recombination the stable equilibria show linkage equilibrium, but several of these with different numbers of polymorphic loci may be stable simultaneously.     

Deviations of Genotypic Structures from Hardy-Weinberg Proportions under Random Mating and Differential Selection between the Sexes

Population genetic models, such as differential viability selection between the sexes and differential multiplicative fecundity contributions of the sexes, are considered for a single multiallelic locus. These selection models usually produce deviations of the zygotic genotype frequencies from Hardy-Weinberg proportions. The deviations are investigated (with special emphasis put on equilibrium states) to quantify the effect of selective asymmetry in the two sexes. For many selection regimes, the present results demonstrate a strong affinity of zygotic genotype frequencies for Hardy-Weinberg proportions after two generations, at the latest. It is shown that the deviations of genotypic equilibria from the corresponding Hardy-Weinberg proportions can be expressed and estimated by means of selection components of only that sex with the lower selection intensity. This corresponds to the well-known fact that viability selection acting in only one sex yields Hardy-Weinberg equilibria.     

The Evolution of One- and Two-Locus Systems. II

Weak selection in a monoecious population is studied in two multiallelic panmictic models. In the first, a single locus is considered with continuous time and age-independent fertilities and mortalities. If the fertilities of the various matings and the genotypic mortalities may be expressed with an error at most of the second order in s (i.e., O(s ²)), where s is the intensity of selection, as sums of terms corresponding to the different genotypes and alleles, respectively, then after several generations the deviations from Hardy-Weinberg proportions are of O(s²). In the second model, two loci are treated with discrete nonoverlapping generations. It is shown that if the epistatic parameters are of O(s²), then after several generations the linkage disequilibria are reduced to O(s²). Assuming only weak selection, it is proved that in both models, after several generations, the total change in mean fitness is generally positive. It is likely that the exclusion of the initial period is usually unnecessary in natural populations. Exceptions are discussed.     

Hardy-Weinberg equilibria in random mating populations

The structure of multiloci random mating populations is examined. Sufficient conditions for the existence of stable local Hardy-Weinberg equilibria for n loci and an arbitrary number of alleles per locus, are then derived for specified situations under the assumption of multiplicative gene action between loci. It is shown that a stable Hardy-Weinberg equilibrium can not be a local maximum of the mean fitness function with multiplicative gene action between loci. The stability of Hardy-Weinberg type border points and the condition for the increase of newly introduced genes are topics on which some n-loci results are also obtained for an arbitrary number of alleles per locus in systems that allow Hardy-Weinberg equilibria.     

A slow selection analysis of two locus, two allele traits.

A deterministic, continuous time model for the dynamics of two locus, two allele Mendelian traits in a large randomly mating diploid population is derived. The model allows for frequency and time dependent birth and death rates. It is analyzed under the assumption that the selective forces acting in the population are small. Slow selection approximations to the system's solution are then constructed. Two particular cases are studied. First, when linkage between loci is tight, the population is shown to rapidly approach Hardy-Weinberg proportions, which then may vary on a (slow) time scale determined by differential fitness. In the case of constant birth and death rates, a measure of the population's fitness is shown to increase on the slow time scale after an initial rapid adjustment. The second case considered is for loose linkage; a population near linkage equilibrium is studied. It is shown that the epistatic parameters cancel and that the results agree with the tight linkage case to leading order. The linkage disequlibrium is described in both cases.     

Comparisons of positive assortative mating and sexual selection models

A series of theoretical models of positive assortative mating and sexual selection are contrasted. It is established that for a dominant trait partial positive assortative mating generally implies some fixation, whereas sexual selection exhibits a unique globally stable polymorphism exhibiting Hardy-Weinberg proportions. The effects of monogamy against polygamy do not qualitatively alter the equilibrium outcomes, although the rate of evolutionary change is generally slowed with monogamy vis-à-vis polygamy. For sexual selection the influence of timing of random mating as against preferential mating causes no change in the equilibrium states, although the rates of convergence can be slowed if sexual selection occurs late in the breeding season. Under assortative mating the timing can alter the equilibrium outcomes. The amount of heterozygosity is always deficient in cases of assortative mating, but always exhibits Hardy-Weinberg ratios under a sexual selection mechanism. This suggests that observations consistent with Hardy-Weinberg equilibrium states cannot preclude ipso facto certain forms of selection forces, including mating patterns and some natural selection structures.     

The multilocus Hardy-Weinberg law

Work on the genetic algebra of multilocus genetic systems is reviewed, with particular emphasis on aspects of Hardy-Weinberg theory, including existence and stability of equilibria, global convergence and disequilibrium functions. It is pointed out that the non-uniqueness of the disequilibrium functions does not necessarily invalidate proof of the multilocus Hardy-Weinberg law.     

Continuous selective models with mutation and migration

The continuous selective model formulated previously for a single locus with multiple alleles in a monoecious population is extended to include mutation and migration. Somatic and germ line genotypic frequencies are distinguished, and the alternative hypotheses of constant mutation rates and age-independent mutation frequencies are analyzed in detail for arbitrary selection and mating schemes. With any mating pattern, if there is no selection, the equilibrium allelic frequencies are shown to be unaffected by the generalizations introduced in this paper. If, in addition, mating is at random, the equilibrium genotypic frequencies are proved to be in Hardy-Weinberg proportions. For both models, the nature of the approach to equilibrium is discussed. Migration is treated in the island model.     

General analysis of frequency-dependent sexual selection at a multi-allelic locus

A general model is analyzed in which arbitrarily frequency-dependent selection acts on one sex of a diploid population with several alleles at one locus, as a result of viability or mating-success differences. The existence of boundary and polymorphic equilibria is examined, and conditions for local stability, internal and external, are obtained. The status of Hardy-Weinberg approximations in studying stability and approach to equilibria is also considered. The general principles are then applied to two specific models: one where genotypes fall into two phenotypic classes; and one with a hierarchy of dominance where viability and sexual selection are opposed. In the latter case it is found that, of all the equilibria present, there is one and only one which could possibly be stable: the existence of a unique globally stable equilibrium might then be inferred.     

A multi-locus continuous-time selection model

Summary A continuous time selection model is formulated for a diploid monoecious population with multiple alleles at each of an arbitrary number of loci, incorporating differential fertility and mortality as well as arbitrary mating and age structure. The model is simplified in the case of age-independence and for the case of a stable age distribution. The age-independent model is examined in detail for the special case of multiple alleles at each of two loci. This model is analyzed under the assumptions of random mating and additive fertilities, with close attention given to the behavior of the system with respect to Hardy-Weinberg proportions and linkage equilibrium.M. M. was supported by a U.S. Public Health Service training grant (Grant No. GM780).     

Statistical Inference for Hardy-Weinberg Proportions in the Presence of Missing Genotype Information

In genetic association studies, tests for Hardy-Weinberg proportions are often employed as a quality control checking procedure. Missing genotypes are typically discarded prior to testing. In this paper we show that inference for Hardy-Weinberg proportions can be biased when missing values are discarded. We propose to use multiple imputation of missing values in order to improve inference for Hardy-Weinberg proportions. For imputation we employ a multinomial logit model that uses information from allele intensities and/or neighbouring markers. Analysis of an empirical data set of single nucleotide polymorphisms possibly related to colon cancer reveals that missing genotypes are not missing completely at random. Deviation from Hardy-Weinberg proportions is mostly due to a lack of heterozygotes. Inbreeding coefficients estimated by multiple imputation of the missings are typically lowered with respect to inbreeding coefficients estimated by discarding the missings. Accounting for missings by multiple imputation qualitatively changed the results of 10 to 17% of the statistical tests performed. Estimates of inbreeding coefficients obtained by multiple imputation showed high correlation with estimates obtained by single imputation using an external reference panel. Our conclusion is that imputation of missing data leads to improved statistical inference for Hardy-Weinberg proportions.     

Host mating system and the prevalence of disease in a plant population

A modified susceptible-infected-recovered (SIR) host-pathogen model is used to determine the influence of plant mating system on the outcome of a host-pathogen interaction. Unlike previous models describing how interactions between mating system and pathogen infection affect individual fitness, this model considers the potential consequences of varying mating systems on the prevalence of resistance alleles and disease within the population. If a single allele for disease resistance is sufficient to confer complete resistance in an individual and if both homozygote and heterozygote resistant individuals have the same mean birth and death rates, then, for any parameter set, the selfing rate does not affect the proportions of resistant, susceptible or infected individuals at equilibrium. If homozygote and heterozygote individual birth rates differ, however, the mating system can make a difference in these proportions. In that case, depending on other parameters, increased selfing can either increase or decrease the rate of infection in the population. Results from this model also predict higher frequencies of resistance alleles in predominantly selfing compared to predominantly outcrossing populations for most model conditions. In populations that have higher selfing rates, the resistance alleles are concentrated in homozygotes, whereas in more outcrossing populations, there are more resistant heterozygotes.     

Selection and mutation at a diallelic X-linked locus

Equilibria and convergence of gene frequencies are studied in the case of a diallelic X-linked locus under the influence of selection and mutation. The model used is that of an infinite diploid population with nonoverlapping discrete generations and random mating. It is proved that if the mutation rates and fitnesses are constant and the mutation rates are less than one-third, then global convergence of gene frequencies to equilibria occurs. The phase portraits of the dynamical system describing the change of allelic frequencies from one generation to the next are determined. Convergence of gene frequencies is monotone from a certain generation on if every other generation is skipped. In the case without mutation, our proof of this monotone convergence simplifies G. Palm's original proof [37].     

Pseudo-Random Mating Populations. in Celebration of the 80th Anniversary of the Hardy-Weinberg Law

That random mating leads to Hardy-Weinberg distribution of genotypes is well known. This report is to show that, if the deviations from random mating are of a certain pattern, the offspring generation will also be in Hardy-Weinberg proportions. This brings out the fact that random mating is a sufficient condition, not a necessary one, for the attainment of the Hardy-Weinberg proportions. Such nonrandom-mating populations are tentatively said to be pseudo-random mating. Pseudo-random-mating populations exist for both autosomal and sex-linked systems with two or multiple alleles. This report covers the basic case of a two-allele autosomal locus in detail, but the possible extension to two loci and cytonuclear systems have also been mentioned in discussion.     

Evolutionary game dynamics with non-uniform interaction rates

The classical setting of evolutionary game theory, the replicator equation, assumes uniform interaction rates. The rate at which individuals meet and interact is independent of their strategies. Here we extend this framework by allowing the interaction rates to depend on the strategies. This extension leads to non-linear fitness functions. We show that a strict Nash equilibrium remains uninvadable for non-uniform interaction rates, but the conditions for evolutionary stability need to be modified. We analyze all games between two strategies. If the two strategies coexist or exclude each other, then the evolutionary dynamics do not change qualitatively, only the location of the equilibrium point changes. If, however, one strategy dominates the other in the classical setting, then the introduction of non-uniform interaction rates can lead to a pair of interior equilibria. For the Prisoner's Dilemma, non-uniform interaction rates allow the coexistence between cooperators and defectors. For the snowdrift game, non-uniform interaction rates change the equilibrium frequency of cooperators.     

Genetics of cat hemoglobins: A quantitative polymorphism

Blood samples from domestic cats, representing the parents (26) and offspring (91) from 33 matings, were analyzed to determine the proportions of hemoglobins A and B. Evidence is presented that the proportions of the two hemoglobins are genetically determined and that cats may be divided into three groups on the basis of this characteristic. Family data and Hardy-Weinberg analysis support the hypothesis that the three groups represent the possible combinations of an allelic gene pair. In an attempt to explain the observed phenotypic differences, possible points of action of this gene pair are discussed, including effects on rates of synthesis or epigenetic modification of globin chains.     

The Maintenance of Genetic Variability in Two-Locus Models of Stabilizing Selection

The maintenance of genetic variability at two diallelic loci under stabilizing selection is investigated. Generations are discrete and nonoverlapping; mating is random; mutation and random genetic drift are absent; selection operates only through viability differences. The determination of the genotypic values is purely additive. The fitness function has its optimum at the value of the double heterozygote and decreases monotonically and symmetrically from its optimum, but is otherwise arbitrary. The resulting fitness scheme is identical to the symmetric viability model. Linkage disequilibrium is neglected, but the results are otherwise exact. Explicit formulas are found for all the equilibria, and explicit conditions are derived fro their existence and stability. A complete classification of the six possible global convergence patterns is presented. In addition to the symmetric equilibrium (with gene frequency 1/2 at both loci), a pair of unsymmetric equilibria may exist; the latter are usually, but not always, unstable. If the ratio of the effect of the major locus to that of the minor one exceeds a critical value, both loci will be stably polymorphic. If selection is weak at the minor locus, the more rapidly the fitness function decreases near the optimum, the lower is this critical value; for rapidly decreasing fitness functions, the critical value is close to one. If the fitness function is smooth at the optimum, then a stable polymorphism exists at both loci only if selection is strong at the major locus.