The consequences of fluctuating selection for isozyme polymorphisms in Daphnia

Temporal sequences of allele frequencies in natural populations of Daphnia are analyzed to obtain the mean and variance of the selection coefficient for both asexual and sexual phases. In general, the alleles at enzyme loci appear to be quasi-neutral. Although significant variation exists for the estimated selection coefficients, the means are in all cases close to zero. Estimates of the variance of selection intensity are applied to existing models to demonstrate the implications of fluctuating selection for the spatial and temporal distribution of gene frequencies in Daphnia. The empirical and analytical results are shown to provide a possible solution to some previously puzzling aspects of Daphnia population genetic surveys. Neither genetic drift nor diversifying selection are necessary conditions for the local diversification of gene frequencies.     

Allozyme variation of Cyclobalanopsis championii (Fagaceae), a narrowly distributed species in southern Taiwan

Allozyme genetic variability in five natural populations of Cyclobalanopsis championii (Fagaceae) in Taiwan was investigated using 12 loci from 9 enzyme systems. The average values of parameters describing within-population variation, expected heterozygosity (He = 0.151), the percentage of polymorphic loci per individual (P = 50%), the average number of alleles per locus (A = 1.7), effective number of alleles per locus (Ae = 1.25), and the average number of alleles per polymorphic loci (AP = 2.2) are comparable to those of other long-lived woody plants. The overall fixation index (Fis = 0.208) indicates a significant deficiency of heterozygotes at the population level. Allelic frequency deviation from Hardy-Weinberg equilibrium was found for different loci in different populations. An exact test for population differentiation using the Tools for Population Genetic Analyses program also indicates that allelic frequencies among populations are significantly different (P < .001). Among-population variation, Gst, accounted for 9.2% of the total heterozygosity. The population at Shouchia and the southernmost population Nanjenshan had higher inbreeding coefficients (0.177 and 0.153, respectively) than did the northern populations. Genetic drift is supported by the observations of the variance components of linkage disequilibrium and a large proportion of loci in Nanjenshan and Shouchia that show pairwise locus disequilibrium. We believe continuous genetic drift in the southern populations will increase genetic divergence among populations of C. championii in Taiwan. Significant correlation was found between elevation and expected heterozygosity. We therefore inferred that temperature is the most important ecological factor to influence the genetic diversity of C. championii.     

Gene action,genetic variation,and GWAS: A user-friendly web tool

Fisher’s partitioning of genotypic values and genetic variance is highly relevant in the current era of genome-wide association studies (GWASs). However, despite being more than a century old, a number of persistent misconceptions related to nonadditive genetic effects remain. We developed a user-friendly web tool, the Falconer ShinyApp, to show how the combination of gene action and allele frequencies at causal loci translate to genetic variance and genetic variance components for a complex trait. The app can be used to demonstrate the relationship between a SNP effect size estimated from GWAS and the variation the SNP generates in the population, i.e., how locus-specific effects lead to individual differences in traits. In addition, it can also be used to demonstrate how within and between locus interactions (dominance and epistasis, respectively) usually do not lead to a large amount of nonadditive variance relative to additive variance, and therefore, that these interactions usually do not explain individual differences in a population.     

Geographic components of linkage disequilibrium in natural populations of Escherichia coli

Geographic variation in the genetic structure of natural enteric populations of Escherichia coli was assessed at both the single-locus and dilocus levels from allozyme genotypes at 12 enzyme loci in 178 cell lines isolated from human hosts in Sweden, Iowa, and Tonga. Although there was significant heterogeneity in allele frequencies at six of the 12 loci, geographic variation accounted for only 2.0% of the total genetic diversity (HT = 0.518). Ohta's D-statistics were used to partition the total variance of dilocus linkage disequilibrium into within-population and between-population components. The observed total variance in disequilibrium (0.0339), averaged over 66 locus-pairs, was significantly greater than would be expected (0.0103) if alleles were randomly associated in an unstructured total population; and both within-locality and between-locality components made substantial contributions to the total variance. Half the locus-pairs exhibited the specific dual relationship among components expected when random factors are generating disequilibrium, but 20% of the locus-pairs showed the opposite relationship, reflecting systematic allele associations. The magnitude of dilocus disequilibrium apparently is unrelated to the chromosomal distance between loci. This and other evidence indicates that substitutive recombination rates in natural populations are sufficiently low to permit indirect periodic selection to play a prominent role in generating multilocus genetic structure.      

Fluxes and metabolic pools as model traits for quantitative genetics. I. The L-shaped distribution of gene effects

The fluxes through metabolic pathways can be considered as model quantitative traits, whose QTL are the polymorphic loci controlling the activity or quantity of the enzymes. Relying on metabolic control theory, we investigated the relationships between the variations of enzyme activity along metabolic pathways and the variations of the flux in a population with biallelic QTL. Two kinds of variations were taken into account, the variation of the average enzyme activity across the loci, and the variation of the activity of each enzyme of the pathway among the individuals of the population. We proposed analytical approximations for the flux mean and variance in the population as well as for the additive and dominance variances of the individual QTL. Monte Carlo simulations based on these approximations showed that an L-shaped distribution of the contributions of individual QTL to the flux variance (R(2)) is consistently expected in an F(2) progeny. This result could partly account for the classically observed L-shaped distribution of QTL effects for quantitative traits. The high correlation we found between R(2) value and flux control coefficients variance suggests that such a distribution is an intrinsic property of metabolic pathways due to the summation property of control coefficients.     

The effect of neutral nonadditive gene action on the quantitative index of population divergence

For neutral additive genes, the quantitative index of population divergence (Q(ST)) is equivalent to Wright's fixation index (F(ST)). Thus, divergent or convergent selection is usually invoked, respectively, as a cause of the observed increase (Q(ST) > F(ST)) or decrease (Q(ST) < F(ST)) of Q(ST) from its neutral expectation (Q(ST) = F(ST)). However, neutral nonadditive gene action can mimic the additive expectations under selection. We have studied theoretically the effect of consecutive population bottlenecks on the difference F(ST) - Q(ST) for two neutral biallelic epistatic loci, covering all types of marginal gene action. With simple dominance, Q(ST) < F(ST) for only low to moderate frequencies of the recessive alleles; otherwise, Q(ST) > F(ST). Additional epistasis extends the condition Q(ST) < F(ST) to a broader range of frequencies. Irrespective of the type of nonadditive action, Q(ST) < F(ST) generally implies an increase of both the within-line additive variance after bottlenecks over its ancestral value (V(A)) and the between-line variance over its additive expectation (2F(ST)V(A)). Thus, both the redistribution of the genetic variance after bottlenecks and the F(ST) - Q(ST) value are governed largely by the marginal properties of single loci. The results indicate that the use of the F(ST) - Q(ST) criterion to investigate the relative importance of drift and selection in population differentiation should be restricted to pure additive traits.     

The effect of epistasis on the excess of the additive and nonadditive variances after population bottlenecks

The effect of population bottlenecks on the components of the genetic variance generated by two neutral independent epistatic loci has been studied theoretically (VA, additive; VD, dominant; VAA, additive x additive; VAD, additive x dominant; VDD; dominant x dominant components of variance). Nonoverdominance and overdominance models were considered, covering all possible types of marginal gene action at the single locus level. The variance components in an infinitely large panmictic population (ancestral components) were compared with their expected values at equilibrium, after t consecutive bottlenecks of equal size N (derived components). Formulae were obtained in terms of allele frequencies and effects at each locus and the corresponding epistatic value. An excess of VA after bottlenecks can be assigned to two sources: (1) the spatiotemporal changes in the marginal average effects of gene substitution alpha(i), which are equal to zero only for additive gene action within and between loci; and (2) the covariance between alpha2(i) and the heterozygosity at the loci involved, which is generated by dominance, with or without epistasis. Numerical examples were analyzed, indicating that an increase in VA after bottlenecks will only occur if its ancestral value is minimal or very small. For the nonoverdominance model with weak reinforcing epistasis, that increase has been detected only for extreme frequencies of the negative allele at one or both loci. With strong epistasis, however, this result can be extended to a broad range of intermediate frequencies. With no epistasis, the same qualitative results were found, indicating that dominance can be considered as the primary cause of an increase in VA following bottlenecks. In parallel, the derived total nonadditive variance exceeded its ancestral value (V(NA) = V(D) + V(AA) + V(AD) + V(DD)) for a range of combinations of allele frequencies covering those for an excess of VA and for very large frequencies of the negative allele at both loci. For the overdominance model, an increase in V(A) and V(NA) was respectively observed for equilibrium (intermediate) frequencies at one or both loci or for extreme frequencies at both loci. For all models, the magnitude of the change of V(A) and V(NA) was inversely related to N and t. At low levels of inbreeding, the between-line variance was not affected by the type of gene action. For the models considered, the results indicate that it is unlikely that the rate of evolution may be accelerated after population bottlenecks, in spite of occasional increments of the derived V(A) over its ancestral value.     

Multiple Allele Approach to the Study of Genes in DROSOPHILA MELANOGASTER That Are Involved in Imaginal Disc Development

The phenotypes of five different lethal mutants of Drosophila melanogaster that have small imaginal discs were analyzed in detail. From these results, we inferred whether or not the observed imaginal disc phenotype resulted exclusively from a primary imaginal disc defect in each mutant. To examine the validity of these inferences, we employed a multiple-allele method. Lethal alleles of the five third-chromosome mutations were identified by screening EMS-treated chromosomes for those which fail to complement with a chromosome containing all five reference mutations. Twenty-four mutants were isolated from 13,197 treated chromosomes. Each of the 24 was then tested for complementation with each of the five reference mutants. There was no significant difference in the mutation frequencies at these five loci. The stage of lethality and the imaginal disc morphology of each mutant allele were compared to those of its reference allele in order to examine the range of defects to be found among lethal alleles of each locus. In addition, hybrids of the alleles were examined for intracistronic complementation. For two of the five loci, we detected no significant phenotypic variation among lethal alleles. We infer that each of the mutant alleles at these two loci cause expression of the null activity phenotype. However, for the three other loci, we did detect significant phenotypic variation among lethal alleles. In fact, one of the mutant alleles at each of these three loci causes no detectable imaginal disc defect. This demonstrates that attempting to assess the developmental role of a gene by studying a single mutant allele may lead to erroneous conclusions. As a byproduct of the mutagenesis procedure, we have isolated two dominant, cold-sensitive mutants.     

Parallel effects of genetic variation in ACE activity in baboons and humans

Like humans, savannah baboons (Papio sp.) show heritable interindividual variation in complex physiological phenotypes. One prominent example of such variation involves production of the homeostatic regulator protein angiotensin converting enzyme (ACE), which shows heritable variation in both baboons and humans. In humans, this phenotypic variation is associated with an Alu insertion-deletion polymorphism in the ACE gene, which explains approximately half of the variation in serum ACE activity. We identified a similar Alu insertion-deletion polymorphism in the baboon ACE homologue and measured its frequency in a wild population and a captive population of baboons. We also analyzed the contribution of ACE genotype at this indel to variation in serum ACE activity in the captive population. When conditioned on weight, a known factor affecting ACE activity in humans, age and ACE genotype both accounted for variance in ACE activity; in particular, we identified a significant nonadditive interaction between age and genotype. A model incorporating this interaction effect explained 21.6% of the variation in residual serum ACE activity. Individuals homozygous for the deletion mutation exhibited significantly higher levels of ACE activity than insertion-deletion heterozygotes at younger ages (10-14 years), but showed a trend towards lower levels of ACE activity compared with heterozygotes at older ages (> or =15 years). These results demonstrate an interesting parallel between the genetic architecture underlying ACE variation in humans and baboons, suggesting that further attention should be paid in humans to the relationship between ACE genetic variation and aging.     

The effect of genetic drift on the variance/covariance components generated by multilocus additive x additive epistatic systems

The effect of population bottlenecks on the components of the genetic variance/covariance generated by n neutral independent additive x additive loci has been studied theoretically. In its simplest version, this situation can be modelled by specifying the allele frequencies and homozygous effects at each locus, and an additional factor measuring the strength of the n-th order epistatic interaction. The variance/covariance components in an infinitely large panmictic population (ancestral components) were compared with their expected values at equilibrium over replicates randomly derived from the base population, after t bottlenecks of size N (derived components). Formulae were obtained giving the derived components (and the between-line variance) as functions of the ancestral ones (alternatively, in terms of allele frequencies and effects) and the corresponding inbreeding coefficient F(t). The n-th order derived component of the genetic variance/covariance is continuously eroded by inbreeding, but the remaining components may increase initially until a critical F(t) value is attained, which is inversely related to the order of the pertinent component, and subsequently decline to zero. These changes can be assigned to the between-line variances/covariances of gene substitution and epistatic effects induced by drift. Numerical examples indicate that: (1) the derived additive variance/covariance component will generally exceed its ancestral value unless epistasis is weak; (2) the derived epistatic variance/covariance components will generally exceed their ancestral values unless allele frequencies are extreme; (3) for systems showing equal ancestral additive and total non-additive variance/covariance components, those including a smaller number of epistatic loci may generate a larger excess in additive variance/covariance after bottlenecks than others involving a larger number of loci, provided that F(t) is low. Our results indicate that it is unlikely that the rate of evolution may be significantly accelerated after population bottlenecks, in spite of occasional increments of the derived additive variance over its ancestral value.     

Influence of Gene Interaction on Complex Trait Variation with Multilocus Models

Although research effort is being expended into determining the importance of epistasis and epistatic variance for complex traits, there is considerable controversy about their importance. Here we undertake an analysis for quantitative traits utilizing a range of multilocus quantitative genetic models and gene frequency distributions, focusing on the potential magnitude of the epistatic variance. All the epistatic terms involving a particular locus appear in its average effect, with the number of two-locus interaction terms increasing in proportion to the square of the number of loci and that of third order as the cube and so on. Hence multilocus epistasis makes substantial contributions to the additive variance and does not, per se, lead to large increases in the nonadditive part of the genotypic variance. Even though this proportion can be high where epistasis is antagonistic to direct effects, it reduces with multiple loci. As the magnitude of the epistatic variance depends critically on the heterozygosity, for models where frequencies are widely dispersed, such as for selectively neutral mutations, contributions of epistatic variance are always small. Epistasis may be important in understanding the genetic architecture, for example, of function or human disease, but that does not imply that loci exhibiting it will contribute much genetic variance. Overall we conclude that theoretical predictions and experimental observations of low amounts of epistatic variance in outbred populations are concordant. It is not a likely source of missing heritability, for example, or major influence on predictions of rates of evolution.     

Estimation of relatedness by DNA fingerprinting

The recent discovery of hypervariable VNTR (variable number of tandem repeat) loci has led to much excitement among population biologists regarding the feasibility of deriving individual estimates of relatedness in field populations by DNA fingerprinting. It is shown that unbiased estimates of relatedness cannot be obtained at the individual level without knowledge of the allelic distributions in both the individuals of interest and the base population unless the proportion of shared marker alleles between unrelated individuals is essentially zero. Since the latter is usually on the order of 0.1-0.5 and since there are enormous practical difficulties in obtaining the former, only an approximate estimator for the relatedness can be given. The bias of this estimator is individual specific and inversely related to the number of marker loci and frequencies of marker alleles. Substantial sampling variance in estimates of relatedness arises from variation in identity by descent within and between loci and, with finite numbers of alleles, from variation in identity in state between genes that are not identical by descent. In the extreme case of 25 assayed loci, each with an effectively infinite number of alleles, the standard error of a relatedness estimate is no less than 14%, 20%, 35%, and 53% of the expectation for full sibs and second-, third-, and fourth-order relationships, respectively. Attempts to ascertain relatedness by means of DNA fingerprinting should proceed with caution.      

The Aymara of Western Bolivia. IV. Gene frequencies for eight blood groups and 19 protein and erythrocyte enzyme systems.

A total of 315 individuals, mainly of Aymara origin, from western Bolivia were examined for genetic variation at eight red cell antigen and 19 serum protein and red cell enzyme loci. The gene frequencies for polymorphic loci and the discovery of several rare variants are discussed in terms of previous work among the Aymara and the closely related Quechua. The effect of inclusion of related individuals in the sample on gene frequency, variance of gene frequency and genetic distance, is discussed.     

GENETIC DIVERSITY AND POPULATION STRUCTURE IN CAMELLIA JAPONICA L. (THEACEAE)

Camellia japonica is a widespread and morphologically diverse tree native to parts of Japan and adjacent islands. Starch gel electrophoresis was used to score allelic variation at 20 loci in seeds collected from 60 populations distributed throughout the species range. In comparison with other plant species, the level of genetic diversity within C. japonica populations is very high: 66.2% of loci were polymorphic on average per population, with a mean number of 2.16 alleles per locus; the mean observed and panmictic heterozygosities were 0.230 and 0.265, respectively. Genotypic proportions at most loci in most populations fit Hardy-Weinberg expectations. However, small heterozygote deficiencies were commonly observed (mean population fixation index = 0.129). It is suggested that the most likely cause of the observed deficiencies is population subdivision into genetically divergent subpopulations. The overall level of population differentiation is greater than is typically observed in out-breeders: The mean genetic distance and identity (Nei's D and I) between pairs of populations were 0.073 and 0.930, respectively, and Wright's Fst was 0.144. Differences among populations appeared to be manifested as variation in gene frequencies at many loci rather than variation in allelic composition per se. However, the patterns of variation were not random. Reciprocal clinal variation of gene frequencies was observed for allele pairs at six loci. In addition, principal components analysis revealed that populations tended to genetically cluster into four regions representing the geographic areas Kyushu, Shikoku, western Honshu, and eastern Honshu. There was a significant relationship between genetic and geographic distance (r = 0.61; P < 0.01). Analysis of variance on allozyme frequencies showed that there was approximately four times as much differentiation among populations within regions, as among regions. It is likely that the observed patterns of population relationships result from the balance between genetic drift in small subpopulations and gene flow between them.     

Heritability of heterozygosity offers a new way of understanding why dominant gene action contributes to additive genetic variance

Whenever allele frequencies are unequal, nonadditive gene action contributes to additive genetic variance and therefore the resemblance between parents and offspring. The reason for this has not been easy to understand. Here, we present a new single‐locus decomposition of additive genetic variance that may give greater intuition about this important result. We show that the contribution of dominant gene action to parent–offspring resemblance only depends on the degree to which the heterozygosity of parents and offspring covary. Thus, dominant gene action only contributes to additive genetic variance when heterozygosity is heritable. Under most circumstances this is the case because individuals with rare alleles are more likely to be heterozygous, and because they pass rare alleles to their offspring they also tend to have heterozygous offspring. When segregating alleles are at equal frequency there are no rare alleles, the heterozygosities of parents and offspring are uncorrelated and dominant gene action does not contribute to additive genetic variance.     

Population structure of the African savannah elephant inferred from mitochondrial control region sequences and nuclear microsatellite loci

Two hundred and thirty-six mitochondrial DNA nucleotide sequences were used in combination with polymorphism at four nuclear microsatellite loci to assess the amount and distribution of genetic variation within and between African savannah elephants. They were sampled from 11 localities in eastern, western and southern Africa. In the total sample, 43 haplotypes were identified and an overall nucleotide diversity of 2.0% was observed. High levels of polymorphism were also observed at the microsatellite loci both at the level of number of alleles and gene diversity. Nine to 14 alleles per locus across populations and 44 alleles in the total sample were found. The gene diversity ranged from 0.51 to 0.72 in the localities studied. An analysis of molecular variance showed significant genetic differentiation between populations within regions and also between regions. The extent of subdivision between populations at the mtDNA control region was approximately twice as high as shown by the microsatellite loci (mtDNA F(ST) = 0.59; microsatellite R(ST) = 0.31). We discuss our results in the light of Pleistocene refugia and attribute the observed pattern to population divergence in allopatry accompanied by a recent population admixture following a recent population expansion.     

Frequencies of null alleles at enzyme Loci in natural populations of ponderosa and red pine

Pinus ponderosa and P. resinosa population samples have mean frequencies of enzymatically inactive alleles of 0.0031 and 0.0028 at 29 and 27 enzyme loci, respectively. Such alleles are rare and are apparently maintained by selection-mutation balance. Ponderosa pine have much higher amounts of allozymic and polygenic phenotypic variation than red pine, yet both species have similar frequencies of null alleles. Thus, null alleles apparently do not contribute to polygenic variation, as has been suggested. The concordance between allozymic and polygenic variation adds support to the view that allozyme studies may be valuable in predicting the relative amount of polygenic variation in populations.     

Data and theory point to mainly additive genetic variance for complex traits

The relative proportion of additive and non-additive variation for complex traits is important in evolutionary biology, medicine, and agriculture. We address a long-standing controversy and paradox about the contribution of non-additive genetic variation, namely that knowledge about biological pathways and gene networks imply that epistasis is important. Yet empirical data across a range of traits and species imply that most genetic variance is additive. We evaluate the evidence from empirical studies of genetic variance components and find that additive variance typically accounts for over half, and often close to 100%, of the total genetic variance. We present new theoretical results, based upon the distribution of allele frequencies under neutral and other population genetic models, that show why this is the case even if there are non-additive effects at the level of gene action. We conclude that interactions at the level of genes are not likely to generate much interaction at the level of variance.     

The Inheritance of Metabolic Flux: Expressions for the within-Sibship Mean and Variance Given the Parental Genotypes

Recent developments have related quantitative trait expression to metabolic flux. The present paper investigates some implications of this for statistical aspects of polygenic inheritance. Expressions are derived for the within-sibship genetic mean and genetic variance of metabolic flux given a pair of parental, diploid, n-locus genotypes. These are exact and hold for arbitrary numbers of gene loci, arbitrary allelic values at each locus, and for arbitrary recombination fractions between adjacent gene loci. The within-sibship, genetic variance is seen to be simply a measure of parental heterozygosity plus a measure of the degree of linkage coupling within the parental genotypes. Approximations are given for the within-sibship phenotypic mean and variance of metabolic flux. These results are applied to the problem of attaining adequate statistical power in a test of association between allozymic variation and inter-individual variation in metabolic flux. Simulations indicate that statistical power can be greatly increased by augmenting the data with predictions and observations on progeny statistics in relation to parental allozyme genotypes. Adequate power may thus be attainable at small sample sizes, and when allozymic variation is scored at a only small fraction of the total set of loci whose catalytic products determine the flux.     

Patterns of Gene Variation in Central and Marginal Populations of DROSOPHILA ROBUSTA

The central and marginal populations of D. robusta differ greatly in the level of inversion polymorphism; the marginal populations are monomorphic or nearly so and the central populations are highly polymorphic. This paper presents the frequencies of alleles at forty gene loci in various populations of D. robusta, studied by electrophoresis of proteins and enzymes. Population samples were obtained from eight widely separated populations of D. robusta which included the central, the extreme marginal and the intervening populations between the center and the margins. We find that the proportion of polymorphic loci and average heterozygosity per individual is slightly higher in the marginal populations than the central populations. In D. robusta on an average, 39% of the loci are polymorphic and the average proportion of loci heterozygous per individual is 11%. A breakdown of loci in three categories, viz, hydrolytic enzymes and some other enzymes, larval proteins and glycolytic and Kreb's cycle enzymes, shows that in all populations the level of polymorphism is highest in the hydrolytic enzymes, intermediate in larval proteins and least in the glycolytic and Kreb's cycle enzymes. On the average, the proportion of loci heterozygous per individual for three groups of loci is: hydrolytic enzymes and others (.164), larval proteins (.115) and glycolytic and Kreb's cycle enzymes (.037). We also observe that in all populations the level of polymorphism on the X chromosome is far less than the expected 38%; in salivary gland cells the euchromatic length of the X chromosome is 38% of the entire genome. Lower levels of polymorphism for the X chromosome loci are explained due to low probability of balanced polymorphisms for the X-linked loci since the conditions for establishment of balanced polymorphism for X-linked loci are more restrictive than for the autosomal loci.-The polymorphic loci can be grouped according to pattern of allele frequencies in different populations as follows: (1) The allele frequencies are similar in all populations at the XDH, Pep-1 and Hex-1 loci. (2) The alleles at the Est-1, Est-2, Amy loci and the AP-4(1.0) and the LAP-1(.90) alleles show north south clinal change in frequency. (3) There is north south and east west differentiation at the Pt-5, Pt-8 and Pt-9 loci and the allele AP-4(.81). (4) Polymorphism at loci such as Fum, B.Ox, Hex-8, Pep-2 and Pep-3 are restricted to only one or two of the populations. (5) Allele frequencies at the MDH and ODH loci fluctuate between populations. (6) Allele frequencies at many polymorphic loci such as Est-1, Est-2, LAP-1, AP-4, Pt-5, Pt-8, Pt-9, Pt-16, MDH, Fum change clinally within a gene arrangement. The pattern of gene variation in D. robusta is very complex and cannot be easily explained due to migration of neutral alleles between once-isolated populations or to semi-isolation of neutral alleles. The observations of the pattern of allele variation in different populations, high levels of polymorphism in the marginal populations which have small population size and low levels of polymorphism of the X chromosome loci all support the argument in favor of balancing selection as the main mechanism for the maintenance of these polymorphisms. Environmental factors must play a role in the maintenance of a great deal of these polymorphisms, since we observe clinal allele frequency changes even within a given inversion type.