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.     

Morphological and genetic characteristics of genetic tumor in Nicotiana glauca

It is well understood that genetic tumors develop in certain interspecific Nicotiana hybrids. Nicotiana species are divided into “plus” and “minus” groups and crosses between “plus” and “minus” species give rise to tumorous hybrids. However, it has been proposed that parents and hybrids derived from crosses among members within the same group do not produce tumors. In this study, genetic tumors were only obtained in Nicotiana glauca, which exhibited tumor features similar to those of N. glauca × N. langsdorffii. Our results suggest that genetic factors may control tumor formation independent of tumor induction dependent on the specific interspecific cross. Genetic tumor formation exhibited high B-type and D-type cyclin expression levels, indicating tumor cells are characterized by an uncontrolled cell cycle.     

The Effects of Density Dependence and Immigration on Local Adaptation and Niche Evolution in a Black-Hole Sink Environment

We examine the effects of density dependence and immigration on local adaptation in a "black-hole sink" habitat, i.e., a habitat in which isolated populations of a species would tend to extinction but where a population is demographically maintained by recurrent one-way migration from a separate source habitat in which the species persists. Using a diploid, one-locus model of a discrete-generation sink population maintained by immigration from a fixed source population, we show that a locally favored allele will spread when rare in the sink if the absolute fitness (or, in some cases, the geometric-mean absolute fitness) of heterozygotes with the favored allele is above one in the sink habitat. With density dependence, the criterion for spread can depend on the rate of immigration, because immigration affects local densities and, hence, absolute fitness. Given the successful establishment of a locally favored allele, it will be maintained by a migration-selection balance and the resulting polymorphic population will be sustained deterministically with either stable or unstable dynamics. The densities of stable polymorphic populations tend to exceed densities that would be maintained in the absence of the favored allele. With strong density regulation, spread of the favored allele may destabilize population dynamics. Our analyses show that polymorphic populations which form subsequent to the establishment of favorable alleles have the capacity to persist deterministically without immigration. Finally, we examined the probabilistic rate at which new favored alleles arise and become established in a sink population. Our results suggest that favored alleles are established most readily at intermediate levels of immigration.     

Models for transmission of disease with immigration of infectives

Simple models for disease transmission that include immigration of infective individuals and variable population size are constructed and analyzed. A model with a general contact rate for a disease that confers no immunity admits a unique endemic equilibrium that is globally stable. A model with mass action incidence for a disease in which infectives either die or recover with permanent immunity has the same qualitative behavior. This latter result is proved by reducing the system to an integro-differential equation. If mass action incidence is replaced by a general contact rate, then the same result is proved locally for a disease that causes fatalities. Threshold-like results are given, but in the presence of immigration of infectives there is no disease-free equilibrium. A considerable reduction of infectives is suggested by the incorporation of screening and quarantining of infectives in a model for HIV transmission in a prison system.     

A semi-symmetric two-locus model

The two-locus symmetric viability model characterized by its invariance with respect to the exchange of alleles at each locus, is a well-studied model of classical two-locus theory. The symmetric model introduced by Lewontin and Kojima is among the few multi-locus models with epistatic interactions between loci for which a polymorphism with linkage equilibrium can be stable and this happens when recombination is sufficiently large. We show that an analogous property holds true for a different model, in which symmetry need exist at only one locus. The properties of this new semi-symmetric model are compared with those of the classical symmetric model. For tight linkage, two classes of polymorphisms are possible, depending on the magnitude of additive epistasis. The recombination rate above which linkage equilibrium becomes stable is derived analytically. As in the symmetric model, intervals of recombination in which no polymorphism is stable are possible, and stable polymorphisms can coexist with stable fixations.     

More on recombination and selection in the modifier theory of sex-ratio distortion

G. Maffi and S.D. Jayakar suggested a model for the two-locus control of sex determination in the mosquito Aedes aegypti (1981, Theor. Pop. Biol. 19, 19-36). This model was extended to multiple alleles and analyzed in mathematical detail by S. Lessard (1987, Theor. Pop. Biol. 31, 339-358). The model supposes that males are "Mm" and females "mm" but the transmission from males is controlled by a second gene with alleles Ai. We show that in addition to the equilibrium in which mAi in females, MAi from males and mAi from males all have the same frequencies, a second class of polymorphic equilibria exists and can be stable. The former class was shown by Lessard to be stable for intermediate and/or loose linkage. The new class of equilibria may be stable for tight linkage under the conditions that preclude stability of the former. We also develop the theory of linkage modification from the neighborhood of the new equilibrium. Successful modifiers of recombination may either reduce or increase the recombination fraction with the outcome depending on the linkage of the modifier to the major genes.     

The deviation from linkage equilibrium with multiple loci varying in a stepping-stone cline

The balance between the creation of associations between alleles at different loci by immigration and the convergence to linkage equilibrium due to the recombination process is studied in a theoretical model. The geographical structure of the model is a stepping-stone chain of populations linking two genetically constant source populations. The model assumes an arbitrary number of autosomal loci and considers genetic variation (two alleles at each locus) that is not subject to natural selection. The gene frequencies at each locus will then show a linear cline through the stepping-stone chain of populations. The deviation from linkage equilibrium through the stepping-stone cline is characterized by an equation for linear measures that provide the linkage disequilibrium measures for a given set of loci in terms of the gene frequencies and the linkage disequilibria in the source populations and in terms of the linkage disequilibrium measures through the cline for lower numbers of loci. Numerical examples of this iterative solution are given, and it is shown that the build-up of the higher order Bennett-disequilibria through the cline is considerably more pronounced than the build-up of two-locus disequilibria.     

Resource-dependent selection.

This paper presents a resource-dependent viability selection differential equation model of continuously reproducing diploid population with two alleles at one locus for a single limiting resource. This model assumes that the genotypic fitness is only a function of the limiting resource. The conditions that the interior equilibrium point of the system exists are that the heterozygote fitness is positive and the homozygote fitness is negative, or the heterozygote fitness is negative and the homozygote fitness is positive at the point. The sufficient and necessary conditions of locally asymptotical stability of the interior equilibrium point are that the heterozygote fitness is positive at the point, or the locally asymptotically stable equilibrium corresponds to the point at which the level of the limiting resource is locally minimized on the zero mean fitness curve, f = 0.     

Some circumstances assuring monomorphism in subdivided populations

It is well known that in a subdivided population subject to soft selection with two alleles at one locus, instability of both fixation states (a “protected polymorphism”) entails at least one stable polymorphic equilibrium. Although stable polymorphic and monomorphic equilibria can coexist in general, a stable fixation state (monomorphic equilibrium) precludes the existence of any polymorphic equilibrium under the circumstances of haploid or submultiplicative diploid viabilities. This provides that a stable monomorphism is robust against random fluctuations in allele frequencies. It also increases the known circumstances where there is a unique globally attracting stable equilibrium, i.e., where allele frequencies are determined by the selection-migration structure independent of the history of the system.     

EPISTASIS AND THE INCREASE IN ADDITIVE GENETIC VARIANCE: IMPLICATIONS FOR PHASE 1 OF WRIGHT'S SHIFTING-BALANCE PROCESS

Central to Wright's shifting-balance theory is the idea that genetic drift and selection in systems with gene interaction can lead to the formation of “adaptive gene complexes.” The theory of genetic drift has been well developed over the last 60 years; however, nearly all of this theory is based on the assumption that only additive gene effects are acting. Wright's theory was developed recognizing that there was a “universality of interaction effects,” which implies that additive theory may not be adequate to describe the process of differentiation that Wright was considering. The concept of an adaptive gene complex implies that an allele that is favored by individual selection in one deme may be removed by selection in another deme. In quantitative genetic terms, the average effects of an allele relative to other alleles changes from deme to deme. The model presented here examines the variance in local breeding values (LBVs) of a single individual and the covariance in the LBVs of a pair of individuals mated in the same deme relative to when they are mated in different demes. Local breeding value is a measure of the average effects of the alleles that make up that individual in a particular deme. I show that when there are only additive effects the covariance between the LBVs of individuals equals the variance in the LBV of an individual. As the amount of epistasis in the ancestral population increases, the variance in the LBV of an individual increases and the covariance between the LBVs of a pair of individuals decreases. The divergence in these two values is a measure of the extent to which the LBV of an individual varies independently of the LBVs of other individuals. When this value is large, it means that the relative ordering of the average effects of alleles will change from deme to deme. These results confirm an important component of Wright's shifting-balance theory: When there is gene interaction, genetic drift can lead to the reordering of the average effects of alleles and when coupled with selection this will lead to the formation of the adaptive gene complexes.     

The equilibrium and stability for n alleles under the density-dependent selection

The problem of the equilibrium under the density-dependent selection for n-alleles belonging to one locus is considered. The new “measure of quality” of the population φ is introduced and it is shown that the equilibrium points are the constrained stationary points of the function φ and all locally stable equilibriums are its local maximum points. The analoque of the Fisher theorem for density-dependent selection is considered.     

HYBRID ZONES WITH DOBZHANSKY-TYPE EPISTATIC SELECTION

Dobzhansky's model of epistatic selection assumes that viable genotypes form “clusters” in genotype space so that populations can evolve from one state to a reproductively isolated state following a “ridge” of well-fit genotypes without crossing any deep adaptive valleys. Recently, the importance of Dobzhansky-type models in evolutionary studies has been reemphasized by Gavrilets (1997a) and Gavrilets and Gravner (1997) who argue that the existence of “ridges” of well-fit genotypes connecting reproductively isolated genotypes is actually a general property of multidimensional adaptive landscapes. Using rigorous techniques and numerical simulations, I analyze clines in the frequencies of selected and neutral alleles maintained by a balance of migration and Dobzhansky-type epistatic selection acting on two diallelic loci. I show that Dobzhansky-type epistatic selection can build up a very strong barrier to neutral gene flow. I describe properties of clines that are indicative of Dobzhansky-type selection.     

Discrete Modeling of Dynamics of Zooplankton Community at the Different Stages of an Antropogeneous Eutrophication

Mathematical modeling is a convenient way for characterization of complex ecosystems. This approach was applied to study the dynamics of zooplankton in Lake Sevan (Armenia) at different stages of anthropogenic eutrophication with the use of a novel method called discrete modeling of dynamical systems with feedback (DMDS). Simulation demonstrated that the application of this method helps in characterization of inter- and intra-component relationships in a natural ecosystem. This method describes all possible pairwise inter-component relationships like “plus–plus,” “minus–minus,” “plus–minus,” “plus–zero,” “minus–zero,” and “zero–zero” that occur in most ecosystems. Based on the results, a working hypothesis was formulated. It was found that the sensitivity to weak external influence in zooplanktons was the greatest during the mid period of eutrophication in Lake Sevan, whereas in the final stages of eutrophication, an outbreak in the biomass production of cyanobacteria was evident. To support this approach, a weak external disturbance in the form of magnetic storm was used to see its effect on species Daphnia longispina sevanica. A statistically significant correlation between the frequency of magnetic storms and the number of this species was revealed and an increase in the number of toxic cyanobacteria species as a consequence of eutrophication. This paper, for the first time, suggests a DMDS method, to diagnose impact of anthropogenic eutrophication on environment.     

Stable linkage disequilibrium without epistasis in subdivided populations

In a large random mating population stable linkage disequilibrium occurs only when there is epistasis. However if a population is divided into a number of subpopulations among which migration occurs, stable linkage disequilibrium in each subpopulation may be produced without epistasis. In the case of two subpopulations a necessary condition for linkage equilibrium in the absence of epistasis is that at least at one of the two loci under consideration the gene frequency must be the same for the two populations. This condition is rather severe and any violation of this will lead to stable linkage disequilibrium. A similar conclusion can be made with more than two populations. In general the presence of linkage disequilibrium does not necessarily imply the existence of epistasis even in equilibrium populations.     

The stationary distribution of an asymmetrical model of selection in a random environment

The stationary distribution for the asymmetrical form of the SAS-CFF model of selection in a random environment is presented. Also presented are the conditions for the stable coexistence of K alleles. These conditions are the same as the conditions obtained from the classical constant-fitness model with the formal substitution of geometric mean fitnesses for the constant fitnesses of the classical model. Two examples are explored. In the “equally spaced” example, increases in the degree of asymmetry raise the homozygosity, which is accompanied by loss of alleles from the population. In the “best allele” example, increases in the degree of asymmetry raise the homozygosity without the loss of alleles. In both cases the frequency spectra are altered by the changes in the degree of asymmetry.     

A SIMULATION STUDY OF MULTILOCUS CLINES

Fisher's method of junctions is used to investigate the degree of association between selected alleles in a cline, in the limit where there is divergence between very many genes. A computer model is used to simulate one of a pair of infinite demes that exchange individuals each generation. Selection is on haploids; it is additive and is equivalent to heterozygote disadvantage. Recombination is uniform over a single chromosome. A “critical value” of selection exists at equilibrium, below which loci act independently and above which they act in association (Barton 1983). Starting with secondary contact, simulation results contrast markedly with the equilibrium solution. The “critical value” is not apparent in the simulated clines, even after many generations. Rather, loci remain associated to some extent under all degrees of selection. The simulation is consistent with the equilibrium analysis in all other respects, and therefore indicates that under weak selection the approach to equilibrium is very slow. This is borne out by further numerical calculations. The slow approach to equilibrium enables us to estimate the time since contact between two demes under idealized conditions. Extending this work toward natural hybrid zones is discussed.     

Allee effects, immigration, and the evolution of species' niches

Theoretical studies of adaptation to sink environments (with conditions outside the niche requirements of a species) have shown that immigration from source habitats can either facilitate or inhibit local adaptation. Here, we examine the influence of immigration on the evolution of local adaptation, given an Allee effect (i.e., at low densities, absolute fitness increases with population density). We consider a deterministic model for evolution at a haploid locus, and a stochastic individual-based model for evolution of a quantitative trait, and several kinds of Allee effects. We demonstrate that increased immigration can greatly facilitate adaptive evolution in the sink; with greater immigration, local population sizes rise, and because of the Allee effect, there is a positive indirect effect of immigration on local fitness. This makes it easier for alleles of modest effect to be captured by natural selection, transforming the sink into a locally adapted population that can persist without immigration.     

Mutation Accumulation May Only Be a Minor Force in Shaping Life-History Traits,Even when Reproduction Is Sexual

In a previous theoretical study we investigated whether adaptive or non-adaptive processes are more important in the evolution of senescence. We built a model that combined both processes and found that mutation accumulation is important only at those ages where mortality has a negligible impact on fitness. This model, however, was limited to haploid organisms. Here we extend our model by introducing diploidy and sexual reproduction. We assume that only recessive (mutated) homozygotes experience detrimental effects. Our results corroborate our previous conclusions, confirming that life histories are largely determined by adaptive processes. We also found that the equilibrium frequencies of mutated alleles are at higher values than in haploid model, because mutations in heterozygotes are hidden for directional selection. Nevertheless, the equilibrium frequencies of recessive homozygotes that make mutations visible to selection are very similar to the equilibrium frequencies of these alleles in our haploid model. Diploidy and sexual reproduction with recombination slows down approaching selection-mutation balance.