On the neutrality of molecular genetic markers: pedigree analysis of genetic variation in fragmented populations

Many studies employ molecular markers to infer ecological and evolutionary processes, assuming that variation found at genetic loci offers a reliable representation of stochastic events in natural populations. Increasingly, evidence emerges that molecular markers might not always be selectively neutral. However, only a few studies have analysed how deviations from neutrality could affect estimates of genetic variation, using populations with known genealogy. We monitored changes in allozyme variation over eight generations in captive metapopulations of the butterfly Bicyclus anynana. Population demography was recorded by individually marking 35 000 butterflies and constructing pedigrees. We designed a computer program that simulated the inheritance of founder allozyme alleles in butterfly pedigrees. We thus tested whether the observed transmission of allozyme alleles could be explained by random genetic drift alone, or whether there was evidence for positive or negative selection. This analysis showed that in the smallest metapopulations the loss of allozyme variation exceeded the neutral rate. Possibly, linkage disequilibria between deleterious mutations and marker alleles resulted in background selection and a faster erosion of allozyme variation. In larger metapopulations, one locus (MDH) showed a significant heterozygote excess and smaller than expected loss in heterozygosity, observations consistent with (associative) overdominance. This study demonstrates that the neutrality of molecular markers cannot always be assumed, particularly in small populations with a high mutation load.     

The mode of evolution of molecular markers in populations of house mice under artificial selection for locomotor behavior

A complete understanding of the mode of evolution of molecular markers is important for making inferences about different population genetic parameters, especially because a number of studies have reported patterns of allelic variation at molecular markers that are not in agreement with neutral evolutionary expectations. In the present study, house mice (Mus domesticus) from the fourteenth generation of a selection experiment for increased voluntary wheel-running activity were used to test how selection on a complex behavior affects the distribution of allelic variation by examining patterns of variation at six microsatellite and four allozyme loci. This population had a hierarchical structure that allowed for simultaneous testing of the effects of selection and genetic drift on the distribution of allelic variation by comparing observed patterns of allele frequencies and estimates of genetic divergence at multiple hierarchical levels to expectations under models of neutral evolution. The levels of genetic divergence among replicate lines and between selection groups, estimated from microsatellite data or pooled microsatellite and allozyme data, were not significantly different from expectations under neutral evolution. Furthermore, the pattern of change of allele frequencies between the base population and generation 14 was largely in agreement with expectations under neutral evolution (although the PGM locus exhibited a pattern of change within populations that was difficult to explain under neutral evolution). Overall the results generally provide support for the neutral evolution of molecular markers.     

Levels of genetic polymorphism: marker loci versus quantitative traits.

Species are the units used to measure ecological diversity and alleles are the units of genetic diversity. Genetic variation within and among species has been documented most extensively using allozyme electrophoresis. This reveals wide differences in genetic variability within, and genetic distances among, species, demonstrating that species are not equivalent units of diversity. The extent to which the pattern observed for allozymes can be used to infer patterns of genetic variation in quantitative traits depends on the forces generating and maintaining variability. Allozyme variation is probably not strictly neutral but, nevertheless, heterozygosity is expected to be influenced by population size and genetic distance will be affected by time since divergence. The same is true for quantitative traits influenced by many genes and under weak stabilizing selection. However, the limited data available suggest that allozyme variability is a poor predictor of genetic variation in quantitative traits within populations. It is a better predictor of general phenotypic divergence and of postzygotic isolation between populations or species, but is only weakly correlated with prezygotic isolation. Studies of grasshopper and planthopper mating signal variation and assortative mating illustrate how these characters evolve independently of general genetic and morphological variation. The role of such traits in prezygotic isolation, and hence speciation, means that they will contribute significantly to the diversity of levels of genetic variation within and among species.     

Population genetic consequences of feeding habits in some forest lepidoptera

By surveying variation at allozyme loci in several phytophagous lepidopteran species (Geometridae), we have tested two hypotheses about the relationship of genetic variation to environmental heterogeneity: (1) that allozyme polymorphisms may exist because of associations between genotypes and "niches" (different host plants, in this instance), and (2) that the overall genetic variation of a species is correlated with environmental heterogeneity (or breadth of the species' overall ecological niche).—Genetic differentiation among samples of oligophagous or polyphagous species taken from different host species was observed in one of three species, at only one of seven polymorphic loci. The data thus provide no evidence for pronounced genetic substructuring, or "host race" formation in these sexually reproducing species, although host plant-genotype associations in a parthenogenetic moth give evidence of the potential for diversifying selection.—In a comparison of allozyme variation in polyphagous ("generalized") and oligophagous ("specialized") species, heterozygosity appeared to be higher in specialized species, at all polymorphic loci but one. It is possible that this unexpected result arises from a functional relation between breadth of diet and genetic variation.     

Natural selection at linked sites in humans

Theoretical and empirical work indicates that patterns of neutral polymorphism can be affected by linked, selected mutations. Under background selection, deleterious mutations removed from a population by purifying selection cause a reduction in linked neutral diversity. Under genetic hitchhiking, the rise in frequency and fixation of beneficial mutations also reduces the level of linked neutral polymorphism. Here we review the evidence that levels of neutral polymorphism in humans are affected by selection at linked sites. We then discuss four approaches for distinguishing between background selection and genetic hitchhiking based on (i) the relationship between polymorphism level and recombination rate for neutral loci with high mutation rates, (ii) relative levels of variation on the X chromosome and the autosomes, (iii) the frequency distribution of neutral polymorphisms, and (iv) population-specific patterns of genetic variation. Although the evidence for selection at linked sites in humans is clear, current methods and data do not allow us to clearly assess the relative importance of background selection and genetic hitchhiking in humans. These results contrast with those obtained for Drosophila, where the signals of positive selection are stronger.     

Inter- and intraspecific genetic and morphological variation in a sibling pair of carabid species

Background

Pogonus littoralis and Pogonus chalceus are very close related species with quite different ecological preferences within salt marshes. We study the evolutionary processes in and between these presumably young species. Therefore, we compare the variation in ecologically relevant characters and the genetic variation within one of the species (intraspecific differentiation) with the variation of the two types of characters between the two species (interspecific variation). Data are compared between two independent sets of populations, one set at a small geographical scale (the ecologically diverse Guérande area in France) and the other set at a Atlantic-Mediterranean scale.

Results

Body and relative wing size and IDH1 allozyme data show that the intraspecific variation in P. chalceus is high and in the same range as the interspecific variation (P. chalceus versus P. littoralis). Based on neutral markers (other allozymes and mitochondrial DNA) on the other hand, the intraspecific variation in P. chalceus is much lower in comparison to the interspecific variation.

Conclusion

The different ecotypes in the highly polytypic species P. chalceus are as highly differentiated in ecological characters as true species, but are not recognised as such by screening neutral DNA polymorphisms. This can be interpreted as a case of ongoing speciation driven by natural selection adapting each ecotype to its respective ecological niche. The same ecological process can be recognised in the differentiation between the two sister species, where en plus reproductive isolation between the two gene pools occurred, allowing independent drift and mutation accumulation in neutral genetic characters.     

Background Selection as Baseline for Nucleotide Variation across the Drosophila Genome

The constant removal of deleterious mutations by natural selection causes a reduction in neutral diversity and efficacy of selection at genetically linked sites (a process called Background Selection, BGS). Population genetic studies, however, often ignore BGS effects when investigating demographic events or the presence of other types of selection. To obtain a more realistic evolutionary expectation that incorporates the unavoidable consequences of deleterious mutations, we generated high-resolution landscapes of variation across the Drosophila melanogaster genome under a BGS scenario independent of polymorphism data. We find that BGS plays a significant role in shaping levels of variation across the entire genome, including long introns and intergenic regions distant from annotated genes. We also find that a very large percentage of the observed variation in diversity across autosomes can be explained by BGS alone, up to 70% across individual chromosome arms at 100-kb scale, thus indicating that BGS predictions can be used as baseline to infer additional types of selection and demographic events. This approach allows detecting several outlier regions with signal of recent adaptive events and selective sweeps. The use of a BGS baseline, however, is particularly appropriate to investigate the presence of balancing selection and our study exposes numerous genomic regions with the predicted signature of higher polymorphism than expected when a BGS context is taken into account. Importantly, we show that these conclusions are robust to the mutation and selection parameters of the BGS model. Finally, analyses of protein evolution together with previous comparisons of genetic maps between Drosophila species, suggest temporally variable recombination landscapes and, thus, local BGS effects that may differ between extant and past phases. Because genome-wide BGS and temporal changes in linkage effects can skew approaches to estimate demographic and selective events, future analyses should incorporate BGS predictions and capture local recombination variation across genomes and along lineages.     

Hitchhiking effect of a beneficial mutation spreading in a subdivided population

A central problem in population genetics is to detect and analyze positive natural selection by which beneficial mutations are driven to fixation. The hitchhiking effect of a rapidly spreading beneficial mutation, which results in local removal of standing genetic variation, allows such an analysis using DNA sequence polymorphism. However, the current mathematical theory that predicts the pattern of genetic hitchhiking relies on the assumption that a beneficial mutation increases to a high frequency in a single random-mating population, which is certainly violated in reality. Individuals in natural populations are distributed over a geographic space. The spread of a beneficial allele can be delayed by limited migration of individuals over the space and its hitchhiking effect can also be affected. To study this effect of geographic structure on genetic hitchhiking, we analyze a simple model of directional selection in a subdivided population. In contrast to previous studies on hitchhiking in subdivided populations, we mainly investigate the range of sufficiently high migration rates that would homogenize genetic variation at neutral loci. We provide a heuristic mathematical analysis that describes how the genealogical structure at a neutral locus linked to the locus under selection is expected to change in a population divided into two demes. Our results indicate that the overall strength of genetic hitchhiking--the degree to which expected heterozygosity decreases--is diminished by population subdivision, mainly because opportunity for the breakdown of hitchhiking by recombination increases as the spread of the beneficial mutation across demes is delayed when migration rate is much smaller than the strength of selection. Furthermore, the amount of genetic variation after a selective sweep is expected to be unequal over demes: a greater reduction in expected heterozygosity occurs in the subpopulation from which the beneficial mutation originates than in its neighboring subpopulations. This raises a possibility of detecting a "hidden" geographic structure of population by carefully analyzing the pattern of a selective sweep.     

Amplified fragment length polymorphism (AFLP) analysis of the genetic structure of the zebra mussel, Dreissena polymorpha, in the Mississippi River

1. We predicted that zebra mussel, Dreissena polymorpha (Pallas), genetic structure in the Mississippi River would follow a model of invasive species genetics, which predicts low genetic structure among populations of recently established species. This prediction was upheld in our previous genetic study using allozymes, however, one locus yielded anomalous results. 2. We employed amplified fragment length polymorphism (AFLP) analysis as a neutral marker to assess the amount of genetic structure within and among populations, and as a test of expected population structure from both invasion genetic theory, and the results from our previous study. 3. There was greater spatial differentiation, as measured by F_st, observed using AFLP's than for allozymes (P < 0.001). There was no evidence that AFLP variation conformed to an isolation by distance model, and genetic relationships of populations, as measured by AFLP markers, were not similar to those detected in our allozyme survey. 4. The lack of concordance between these two genetic marker systems probably reflects their differential responses to drift, migration, and selection occurring during this rapid invasion. Strong population structure is counter to predictions that populations of invasive species will not be differentiated, as with observations based on allozyme markers. Therefore, newly established species may require genetic surveys using multiple marker systems to evaluate population structure.     

Variation after a selective sweep in a subdivided population

The effect of genetic hitchhiking on neutral variation is analyzed in subdivided populations with differentiated demes. After fixation of a favorable mutation, the consequences on particular subpopulations can be radically different. In the subpopulation where the mutation first appeared by mutation, variation at linked neutral loci is expected to be reduced, as predicted by the classical theory. However, the effect in the other subpopulations, where the mutation is introduced by migration, can be the opposite. This effect depends on the level of genetic differentiation of the subpopulations, the selective advantage of the mutation, the recombination frequency, and the population size, as stated by analytical derivations and computer simulations. The characteristic outcomes of the effect are three. First, the genomic region of reduced variation around the selected locus is smaller than that predicted in a panmictic population. Second, for more distant neutral loci, the amount of variation increases over the level they had before the hitchhiking event. Third, for these loci, the spectrum of gene frequencies is dominated by an excess of alleles at intermediate frequencies when compared with the neutral theory. At these loci, hitchhiking works like a system that takes variation from the between-subpopulation component and introduces it into the subpopulations. The mechanism can also operate in other systems in which the genetic variation is distributed in clusters with limited exchange of variation, such as chromosome arrangements or genomic regions closely linked to targets of balancing selection.     

植物种群分子进化中对生境的适应——对荒漠植物遗传多样性研究结果的初步分析

长期以来,自然选择理论与中性理论对生物分子进化中的环境适应机理存在着激烈争论。目前,在植物种群分子进化中对生境适应的研究中正面临着一些难题：中性突变是分子水平进化的唯一原因,自然选择发挥主要作用的适应性进化是否存在于分子水平,选择与中性两种学说两种机制完全不同,如何才能将两者联系和统一起来,部分学者利用建立各种模型来描述自然选择对分子标记位点以及连锁序列的直接作用,如生态位宽度变异假设等。本研究小组以新疆阜康荒漠植物为研究对象,通过对两种重要荒漠植物遗传多样性的研究,分析两种植物各亚种群不同生境的生态因子与其遗传变异的关系,讨论生态位宽度变异假设,揭示遗传变异的产生与维持。中性论者与选择论者都试图解释生物环境适应与分子变异之间的关系。中性论和选择论是反映进化的两个侧面,它们不是绝对的,可以相互转化。     

Soft sweeps: molecular population genetics of adaptation from standing genetic variation

A population can adapt to a rapid environmental change or habitat expansion in two ways. It may adapt either through new beneficial mutations that subsequently sweep through the population or by using alleles from the standing genetic variation. We use diffusion theory to calculate the probabilities for selective adaptations and find a large increase in the fixation probability for weak substitutions, if alleles originate from the standing genetic variation. We then determine the parameter regions where each scenario-standing variation vs. new mutations-is more likely. Adaptations from the standing genetic variation are favored if either the selective advantage is weak or the selection coefficient and the mutation rate are both high. Finally, we analyze the probability of "soft sweeps," where multiple copies of the selected allele contribute to a substitution, and discuss the consequences for the footprint of selection on linked neutral variation. We find that soft sweeps with weaker selective footprints are likely under both scenarios if the mutation rate and/or the selection coefficient is high.     

Analysis of F ST outliers at allozyme loci in Pacific salmon: implications for natural selection

Natural selection has been invoked to explain the observed geographic distribution of allozyme allele frequencies for a number of teleost species. The effects of selection on allozyme loci in three species of Pacific salmon were tested. A simulation-based approach to estimate the null distribution of population differentiation (F _ST) and test for F _ST outliers was used. This approach showed that a majority of allozyme loci conform to neutral expectations predicted by the simulation model, with relatively few F _ST outliers found. No consistent F _ST outlier loci were found across species. Analysis of population sub-groups based on geography and genetic identity reduced the number of outlier loci for some species, indicating that large geographic groups may include genetically divergent populations and/or that there is geographic heterogeneity in selection pressure upon allozyme loci. Two outlier allozyme loci found in this analysis, lactate dehydrogenase-B and malic enzyme, have been shown to be influenced by selection in other teleost species. This approach is also useful in identifying allozyme loci (or other genetic markers) that meet assumptions for population genetic study.     

Regional subdivision in wild barley allozyme variation: adaptive or neutral?

We examined the adaptive importance of allozyme variation in wild barley (Hordeum spontaneum). The test involved a nested sampling design with four population groups, each representing a different environment, and a comparison of observed allozyme variation with that expected under the assumption that allozymes are not neutral. Measurements of plant fitness in indigenous and alien environments in reciprocal introductions of seeds and seedlings in the four environments provided a guideline for the expected pattern of allozyme variation. The results showed considerable variation in both the degree of regional and population subdivision and the pattern of the subdivision among loci. The observed pattern of variation was ambiguous. Although two alleles exhibited a pattern of distribution that cannot be explained by genetic drift as a function of geographic distance, we failed to detect either a significant relationship between genetic distance and environmental similarity or any favored epistatic allele combinations across the four environments. Our results suggest that interpretation of allozyme variation in wild barley as adaptive and directly related to local environment still needs justification. Although we could not reject the null hypothesis, a proposed methodology seeking a concordance between observed and "adaptive" (i.e., expected under hypothesis that allozymes are not neutral) allozyme variation may prove to be effective in resolving the neutralist-selectionist debate when applied to other species.     

Comparative losses of quantitative and molecular genetic variation in finite populations of Drosophila melanogaster

Quantitative genetic variation, the main determinant of the ability to evolve, is expected to be lost in small populations, but there are limited data on the effect, and controversy as to whether it is similar to that for near neutral molecular variation. Genetic variation for abdominal and sternopleural bristle numbers and allozyme heterozygosity were estimated in 23 populations of Drosophila melanogaster maintained at effective population sizes of 25, 50, 100, 250 or 500 for 50 generations, as well as in 19 highly inbred populations and the wild outbred base population. Highly significant negative regressions of proportion of initial genetic variation retained on inbreeding due to finite population size were observed for both quantitative characters (b = -0.67 +/- 0.14 and -0.58 +/- 0.11) and allozyme heterozygosity (b = -0.79 +/- 0.10), and the regression coefficients did not differ significantly. Thus, quantitative genetic variation is being lost at a similar rate to molecular genetic variation. However, genetic variation for all traits was lost at rates significantly slower than predicted by neutral theory, most likely due to associative overdominance. Positive, but relatively low correlations were found among the different measures of genetic variation, but their low magnitudes were attributed to large sampling errors, rather than differences in the underlying processes of loss.     

One species or several? Discordant patterns of geographic variation between allozymes and mtDNA sequences among spiders in the genus Metepeira (Araneae: Araneidae)

Paradoxically, an allozyme study of Metepeira "spinipes" (sensu lato) demonstrated extensive gene flow among four populations whose members are nevertheless morphologically and behaviorally distinct. Initially, the authors tentatively concluded that the populations exhibited panmixis and suggested that local environmental effects accounted for the apparent morphological and behavioral differences. However, they later concluded that such differences were too great to be accounted for by the environment alone and that the four populations actually represented three different species. To confirm that the allozyme results were, in fact, artifactual, we reexamined the relationships among these populations by sequencing a portion of the 12S mtDNA ribosomal subunit. In contrast to the allozyme result, our results demonstrate good agreement between patterns of genetic and morphological/behavioral variation. We suggest (1) that the allozyme allele frequencies are homogenized by balancing selection, not gene flow as was previously concluded, and therefore (2) that this study provides another instance in which inferences about population structures from allozyme data are misleading.     

MICROSATELLITE EVOLUTION IN VERTEBRATES: INFERENCE FROM AC DINUCLEOTIDE REPEATS

Abstract We analyze published data from 592 AC microsatellite loci from 98 species in five vertebrate classes including fish, reptiles, amphibians, birds, and mammals. We use these data to address nine major questions about microsatellite evolution. First, we find that larger genomes do not have more microsatellite loci and therefore reject the hypothesis that microsatellites function primarily to package DNA into chromosomes. Second, we confirm that microsatellite loci are relatively rare in avian genomes, but reject the hypothesis that this is due to physical constraints imposed by flight. Third, we find that microsatellite variation differs among species within classes, possibly relating to population dynamics. Fourth, we reject the hypothesis that microsatellite structure (length, number of alleles, allele dispersion, range in allele sizes) differs between poikilotherms and homeotherms. The difference is found only in fish, which have longer microsatellites and more alleles than the other classes. Fifth, we find that the range in microsatellite allele size at a locus is largely due to the number of alleles and secondarily to allele dispersion. Sixth, length is a major factor influencing mutation rate. Seventh, there is a directional mutation toward an increase in microsatellite length. Eighth, at the species level, microsatellite and allozyme heterozygosity covary and therefore inferences based on large-scale studies of allozyme variation may also reflect microsatellite genetic diversity. Finally, published microsatellite loci (isolated using conventional hybridization methods) provide a biased estimate of the actual mean repeat length of microsatellites in the genome.     

Evidence of a reproductive barrier between two forms of the marine periwinkle Littorina fabalis (Gastropoda)

Studies of allozyme variation may reveal unexpected patterns of genetic variation which challenge earlier conclusions of species delimitations based on morphological data. However, allozyme variation alone may not be sufficient to resolve this kind of problem. For example, populations of the marine intertidal snail Littorina fabalis (=Littorina mariae) from wave exposed parts end from protected parts of the same shores are distinguished by different alleles of arginine kinase (Ark) while indifferent, or very nearly so, in another 29 loci. Intermediate populations have large deficiencies of exposed/sheltered heterozygote classes of Ark and we have earlier suggested habitat-related selection in this locus as the explanation. In this study we estimated growth rate of individual snails of different Ark-genotypes in three different habitats (exposed, sheltered and intermediate). In all habitats the snails homozygous for alleles of ‘exposed’ type grew faster and matured at a larger size than did snails homozygous for alleles of ‘sheltered’ types. This relationship was indirectly confirmed in three additional sites of intermediate exposure where exposed AA-genotypes dominated among large (>8 mm) snails while the sheltered genotypes dominated among small (<5 mm) snails of truly sympatric samples. We furthermore found small differences in allele frequencies of two other loci (Pgi and Pgm-2) and in shell colour frequencies, comparing sympatric snails of exposed and sheltered Ark-homozygotes. Although we found no signs of habitat-related selection among snails of different Ark-genotype, or selection against heterozygotes, we cannot reject selection in Ark, as our experiments only covered one island, one season and grown-up snails. The coupling between allozyme and phenotypic characters in strictly sympatric samples of snails suggests the presence of two gene pools. Perhaps the large and small forms of L. fabalis represent very closely related cryptic taxa. However, introgression between them seems a possible explanation for the striking similarities in the vast majority of morphological and allozyme characters.     

THE DIVERGENCE OF NEUTRAL QUANTITATIVE CHARACTERS AMONG PARTIALLY ISOLATED POPULATIONS

The neutral model of phenotypic evolution has yielded several simple predictions about the long-term rates of between-population divergence of polygenic traits and about the equilibrium level of within-population variance when mutation and random genetic drift are the sole evolutionary forces. These conclusions must be modified if populations are only partially isolated. A quantitative model is presented for the development of within- and between-population variance for neutral quantitative characters in pairs of populations with arbitrary effective sizes and migration rates. Both the variance in the base population and subsequent variance generated by mutation are considered, and several dynamical and equilibrium properties are shown to be adequately described by simple approximations. The resultant formulations provide some insight into the sensitivity of measures of morphological distance to gene flow, the necessity of isolation for the accumulation of variation between incipient species, and the consequences of gene flow into captive populations of endangered species.     

Functional Effects of Pgi Allozymes in ESCHERICHIA COLI

Five alleles representing three electromorphs of phosphoglucose isomerase (PGI) have been transferred from natural isolates of E. coli into the genetic background of E. coli K12 and examined for their effect on growth rate in chemostats limited for glucose or fructose. With glucose limitation, all alleles are selectively neutral or nearly neutral within the limit of resolution of the technique, whether the genetic background is nonmutant or whether it contains a deletion of the locus of glucose-6-phosphate dehydrogenase, the enzyme that provides an alternative metabolic pathway for the substrate of PGI. With fructose limitation, one of the naturally occurring alleles has a small but reproducible detrimental effect on growth rate. A kinetic difference in this detrimental allozyme, apparently relating to an inhibition constant, has been observed in some, but not all, lots of substrate, and a similar difference has also been noted in one of the rare electromorphs that could not be transferred into E. coli K12. These results support a model of genetic variation in which the alleles are neutral or nearly neutral in the prevailing environment but have a potential for selection that can be expressed under the appropriate conditions of environment or genetic background. This hypothesis is discussed in the context of allozyme polymorphisms observed in other organisms.