Genomic organization of molecular differentiation in Norway spruce (Picea abies)

Diversity and differentiation among three populations representing the geographical domains commonly recognized within the natural distribution area of Picea abies were analysed by using a set of 292 AFLP (amplified fragment length polymorphism), SSR (single sequence repeat) and ESTP (expressed sequence tags polymorphism) markers. As usually observed in forest trees, results showed high within-population diversity (H(S) reaching 0.79) and low among-population differentiation (G(ST) approximately 2%). The genomic organization of differentiation was then investigated on the basis of a subsample of 150 AFLP, SSR and ESTP mapped markers. The number of the loci differentiating the Baltico-Nordic from the central European populations (25 loci) and, within the central European populations, the Alpine from the Hercyno-Carpathian populations (12 loci), were different. These 37 differentiated loci, with individual G(ST) values ranging from 0.008 to 0.20, were evenly distributed on all linkage groups and mostly followed the neutral expectations, suggesting genome-wide effects on differentiation. Nine of them however behave as 'outlier' loci indicating possible locus-specific selective effects. Contribution of ongoing evolutionary forces and historical effects to the geographical differentiation of the species are discussed.     

Evidence for contrasting modes of selection at interacting globin genes in the European rabbit (Oryctolagus cuniculus)

In hybrid zones between genetically differentiated populations, variation in locus-specific rates of introgression may reflect adaptation to different environments or adaptation to different genetic backgrounds. The European rabbit, Oryctolagus cuniculus, is well-suited to studies of such hybrid zone dynamics because it is composed of two genetically divergent subspecies that hybridize in a zone of secondary contact in central Iberia. A species-wide survey of allozyme variation revealed a broad range of locus-specific divergence levels (F(ST) ranged from 0 to 0.54, mean F(ST)=0.16). Interestingly, the two loci that fell at opposite ends of the distribution of F(ST) values, haemoglobin alpha-chain (HBA) and haemoglobin beta-chain (HBB), encode interacting subunits of the haemoglobin protein. The contrasting patterns of spatial variation at these two loci could not be reconciled under a neutral model of population structure. The HBA gene exhibited higher-than-expected levels of population differentiation, consistent with a history of spatially varying selection. The HBB gene exhibited lower-than-expected levels of population differentiation, consistent with some form of spatially uniform selection. Patterns of linkage disequilibrium and allele frequency variation do not appear to fit any simple model of two-locus epistatic selection.     

On the distribution of temporal variations in allele frequency: consequences for the estimation of effective population size and the detection of loci undergoing selection

The effective population size (Ne) is frequently estimated using temporal changes in allele frequencies at neutral markers. Such temporal changes in allele frequencies are usually estimated from the standardized variance in allele frequencies (Fc). We simulate Wright-Fisher populations to generate expected distributions of Fc and of Fc (Fc averaged over several loci). We explore the adjustment of these simulated Fc distributions to a chi-square distribution and evaluate the resulting precision on the estimation of Ne for various scenarios. Next, we outline a procedure to test for the homogeneity of the individual Fc across loci and identify markers exhibiting extreme Fc-values compared to the rest of the genome. Such loci are likely to be in genomic areas undergoing selection, driving Fc to values greater (or smaller) than expected under drift alone. Our procedure assigns a P-value to each locus under the null hypothesis (drift is homogeneous throughout the genome) and simultaneously controls the rate of false positive among loci declared as departing significantly from the null. The procedure is illustrated using two published data sets: (i) an experimental wheat population subject to natural selection and (ii) a maize population undergoing recurrent selection.     

Adaptation to a steep environmental gradient and an associated barrier to gene exchange in Littorina saxatilis

Steep environmental gradients offer important opportunities to study the interaction between natural selection and gene flow. Allele frequency clines are expected to form at loci under selection, but unlinked neutral alleles may pass easily across these clines unless a generalized barrier evolves. Here we consider the distribution of forms of the intertidal gastropod Littorina saxatilis, analyzing shell shape and amplified fragment length polymorphism (AFLP) loci on two rocky shores in Britain. On the basis of previous work, the AFLP loci were divided into differentiated and undifferentiated groups. On both shores, we have shown a sharp cline in allele frequencies between the two morphs for differentiated AFLP loci. This is coincident with a habitat transition on the shore where the two habitats (cliff and boulder field) are immediately contiguous. The allele frequency clines coincide with a cline in shell morphology. In the middle of the cline, linkage disequilibrium for the differentiated loci rises in accordance with expectation. The clines are extremely narrow relative to dispersal, probably as a result of both strong selection and habitat choice. An increase in F(ST) for undifferentiated AFLPs between morphs, relative to within-morph comparisons, is consistent with there being a general barrier to gene flow across the contact zone. These features are consistent either with an episode of allopatric divergence followed by secondary contact or with primary, nonallopatric divergence. Further data will be needed to distinguish between these alternatives.     

Evolutionary inference from QST

Q(ST) is a commonly used metric of the degree of genetic differentiation among populations displayed by quantitative traits. Typically, Q(ST) is compared to F(ST) measured on putatively neutral loci; if Q(ST)=F(ST), this is taken as evidence of spatially heterogeneous and diversifying selection. This paper reviews the uses, assumptions and statistics of Q(ST) and F(ST) comparisons. Unfortunately, Q(ST)/F(ST) comparisons are statistically challenging. For a single trait, Q(ST) must be compared not to the mean F(ST) but to the distribution of F(ST) values. The sources of biases and sampling error for Q(ST) are reviewed, and a new method for comparing Q(ST) and F(ST) is suggested. Simulation results suggest that the distribution of neutral F(ST) and Q(ST) values are little affected by various deviations from the island model. Consequently, the distributions of Q(ST) and F(ST) are well approximated by the Lewontin-Krakauer prediction, even with realistic deviations from the island-model assumptions.     

Natural selection on protein polymorphism in the rodent genus Peromyscus: evidence from interlocus contrasts

The effects of natural selection are generally locus-specific, whereas migration, drift, and inbreeding are expected to have relatively uniform effects across the entire genome. This suggests that multilocus surveys of multiple populations can be used to distinguish selection from demographic effects. The purpose of this study was to test for evidence of selection on protein polymorphism in natural populations of mice in the genus Peromyscus. We analyzed published data from geographic surveys of allozyme variation and used a coalescent-based simulation model to identify specific loci that deviated from neutral expectations. Observed F(ST) values generally exhibited a remarkably close fit to the expected neutral distributions, indicating that the majority of loci are simply tracking stochastic demographic processes. A smaller number of loci exhibited highly significant departures from the expectations of the neutral model and thus appear to be tracking the direct or indirect effects of selection. Most departures from neutrality were characterized by F(ST) values that far exceeded neutral expectations and were therefore attributable to spatially varying selection. Interestingly, the albumin locus was implicated as a candidate gene for local adaptation in four different species of Peromyscus. The results also demonstrate that selection can severely bias marker-based estimates of neutral parameters.     

Empirical Bayes inference of pairwise F(ST) and its distribution in the genome

Populations often have very complex hierarchical structure. Therefore, it is crucial in genetic monitoring and conservation biology to have a reliable estimate of the pattern of population subdivision. F(ST)'s for pairs of sampled localities or subpopulations are crucial statistics for the exploratory analysis of population structures, such as cluster analysis and multidimensional scaling. However, the estimation of F(ST) is not precise enough to reliably estimate the population structure and the extent of heterogeneity. This article proposes an empirical Bayes procedure to estimate locus-specific pairwise F(ST)'s. The posterior mean of the pairwise F(ST) can be interpreted as a shrinkage estimator, which reduces the variance of conventional estimators largely at the expense of a small bias. The global F(ST) of a population generally varies among loci in the genome. Our maximum-likelihood estimates of global F(ST)'s can be used as sufficient statistics to estimate the distribution of F(ST) in the genome. We demonstrate the efficacy and robustness of our model by simulation and by an analysis of the microsatellite allele frequencies of the Pacific herring. The heterogeneity of the global F(ST) in the genome is discussed on the basis of the estimated distribution of the global F(ST) for the herring and examples of human single nucleotide polymorphisms (SNPs).     

Comparative analysis of population genetic structure in Athyrium distentifolium (Pteridophyta) using AFLPs and SSRs from anonymous and transcribed gene regions

To examine the performance and information content of different marker systems, comparative assessment of population genetic diversity was undertaken in nine populations of Athyrium distentifolium using nine genomic and 10 expressed sequence tag (EST) microsatellite (SSR) loci, and 265 amplified fragment length polymorphism (AFLP) loci from two primer combinations. In range-wide comparisons (European vs. North American populations), the EST-SSR loci showed more reliable amplification and produced more easily scorable bands than genomic simple sequence repeats (SSRs). Genomic SSRs showed significantly higher levels of allelic diversity than EST-SSRs, but there was a significant correlation in the rank order of population diversities revealed by both marker types. When AFLPs, genomic SSRs, and EST-SSRs are considered, comparisons of different population diversity metrics/markers revealed a mixture of significant and nonsignificant rank-order correlations. However, no hard incongruence was detected (in no pairwise comparison of populations did different marker systems or metrics detect opposingly significant different amounts of variation). Comparable population pairwise estimates of F(ST) were obtained for all marker types, but whilst absolute values for genomic and EST-SSRs were very similar (F(ST) = 0.355 and 0.342, respectively), differentiation was consistently higher for AFLPs in pairwise and global comparisons (global AFLP F(ST) = 0.496). The two AFLP primer combinations outperformed 18 SSR loci in assignment tests and discriminatory power in phenetic cluster analyses. The results from marker comparisons on A. distentifolium are discussed in the context of the few other studies on natural plant populations comparing microsatellite and AFLP variability.     

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.     

A population genetics study of Anopheles darlingi (Diptera: Culicidae) from Colombia based on random amplified polymorphic DNA-polymerase chain reaction and amplified fragment lenght polymorphism markers

The genetic variation and population structure of three populations of Anopheles darlingi from Colombia were studied using random amplified polymorphic markers (RAPDs) and amplified fragment length polymorphism markers (AFLPs). Six RAPD primers produced 46 polymorphic fragments, while two AFLP primer combinations produced 197 polymorphic fragments from 71 DNA samples. Both of the evaluated genetic markers showed the presence of gene flow, suggesting that Colombian An. darlingi populations are in panmixia. Average genetic diversity, estimated from observed heterozygosity, was 0.374 (RAPD) and 0.309 (AFLP). RAPD and AFLP markers showed little evidence of geographic separation between eastern and western populations; however, the F ST values showed high gene flow between the two western populations (RAPD: F ST = 0.029; Nm: 8.5; AFLP: F ST = 0.051; Nm: 4.7). According to molecular variance analysis (AMOVA), the genetic distance between populations was significant (RAPD:phiST = 0.084; AFLP:phiST = 0.229, P < 0.001). The F ST distances and AMOVAs using AFLP loci support the differentiation of the Guyana biogeographic province population from those of the Chocó-Magdalena. In this last region, Chocó and Córdoba populations showed the highest genetic flow.     

Genetic diversity analysis using lowly polymorphic dominant markers: the example of AFLP in pigs

DNA markers are commonly used for large-scale evaluation of genetic diversity in farm animals, as a component of the management of animal genetic resources. AFLP markers are useful for such studies as they can be generated relatively simply; however, challenges in analysis arise from their dominant scoring and the low level of polymorphism of some markers. This paper describes the results obtained with a set of AFLP markers in a study of 59 pig breeds. AFLP fingerprints were generated using four primer combinations (PC), yielding a total of 148 marker loci, and average harmonic mean of breed sample size was 37.3. The average proportion of monomorphic populations was 63% (range across loci: 3%-98%). The moment-based method of Hill and Weir (2004, Mol Ecol 13:895-908) was applied to estimate gene frequencies, gene diversity (F(ST)), and Reynolds genetic distances. A highly significant average F(ST) of 0.11 was estimated, together with highly significant PC effects on gene diversity. The variance of F(ST) across loci also significantly exceeded the variance expected under the hypothesis of AFLP neutrality, strongly suggesting the sensitivity of AFLP to selection or other forces. Moment estimates were compared to estimates derived from the square root estimation of gene frequency, as currently applied for dominant markers, and the biases incurred in the latter method were evaluated. The paper discusses the hypotheses underlying the moment estimations and various issues relating to the biallelic, dominant, and lowly polymorphic nature of this set of AFLP markers and to their use as compared to microsatellites for measuring genetic diversity.     

Adaptive divergence for a fitness-related trait among invasive Ambrosia artemisiifolia populations in France

The impact of natural selection on the adaptive divergence of invasive populations can be assessed by testing the null hypothesis that the extent of quantitative genetic differentiation (Q(ST) ) would be similar to that of neutral molecular differentiation (F(ST) ). Using eight microsatellite loci and a common garden approach, we compared Q(ST) and F(ST) among ten populations of an invasive species Ambrosia artemisiifolia (common ragweed) in France. In a common garden study with varying water and nutrient levels, we measured Q(ST) for five traits (height, total biomass, reproductive allocation, above- to belowground biomass ratio, and days to flowering). Although low F(ST) indicated weak genetic structure and strong gene flow among populations, we found significant diversifying selection (Q(ST) > F(ST) ) for reproductive allocation that may be closely related to fitness. It suggests that abiotic conditions may have exerted selection pressure on A. artemisiifolia populations to differentiate adaptively, such that populations at higher altitude or latitude evolved greater reproductive allocation. As previous studies indicate multiple introductions from various source populations of A. artemisiifolia in North America, our results suggest that the admixture of introduced populations may have increased genetic diversity and additive genetic variance, and in turn, promoted the rapid evolution and adaptation of this invasive species.     

Among-locus variation in Fst: fish, allozymes and the Lewontin-Krakauer test revisited.

Variation among loci in the distribution of allele frequencies among subpopulations is well known; how to tell when the variation exceeds that expected when all loci are subject to uniform evolutionary processes is not well known. If locus-specific effects are important, the ability to detect those effects should vary with the level of gene flow. Populations with low gene flow should exhibit greater variation among loci in Fst than populations with high gene flow, because gene flow acts to homogenize allele frequencies among subpopulations. Here I use Lewontin and Krakauer's k statistic to describe the variance among allozyme loci in 102 published data sets from fishes. As originally proposed, k > 2 was considered evidence that the variation in Fst among loci is greater than expected from neutral evolution. Although that interpretation is invalid, large differences in k in different populations suggest that locus-specific forces may be important in shaping genetic diversity. In these data, k is not greater for populations with expected low levels of gene flow than for populations with expected high levels of gene flow. There is thus no evidence that locus-specific forces are of general importance in shaping the distribution of allele frequencies at enzyme loci among populations of fishes.     

Identifying and reducing AFLP genotyping error: an example of tradeoffs when comparing population structure in broadcast spawning versus brooding oysters

Phylogeographic inferences about gene flow are strengthened through comparison of co-distributed taxa, but also depend on adequate genomic sampling. Amplified fragment length polymorphisms (AFLPs) provide a rapid and inexpensive source of multilocus allele frequency data for making genomically robust inferences. Every AFLP study initially generates markers with a range of locus-specific genotyping error rates and applies criteria to select a subset for analysis. However, there has been very little empirical evaluation of the best tradeoff between culling all but the lowest-error loci to minimize overall genotyping error versus the potential for increasing population genetic signal by retaining more loci. Here, we used AFLPs to compare population structure in co-distributed broadcast spawning (Crassostrea virginica) and brooding (Ostrea equestris) oyster species. Using existing methods for almost entirely automated marker selection and scoring, genotyping error tradeoffs were evaluated by comparing results across a nested series of data sets with mean mismatch errors of 0, 1, 2, 3, 4 and >4%. Artifactual population structure was diagnosed in high-error data sets and we assessed the low-error point at which expected population substructure signal was lost. In both species, we identified substructure patterns deemed to be inaccurate at average mismatch error rates 2 and >4%. In the species comparison, the optimum data sets showed higher gene flow for the brooding oyster with more oceanic salinity tolerances. AFLP tradeoffs may differ among studies, but our results suggest that important signal may be lost in the pursuit of 'acceptable' error levels and our procedures provide a general method for empirically exploring these tradeoffs.     

Controlling for the effects of history and nonequilibrium conditions in gene flow estimates in northern bullfrog (Rana catesbeiana) populations

Nonequilibrium conditions due to either allopatry followed by secondary contact or recent range expansion can confound measurements of gene flow among populations in previously glaciated regions. We determined the scale at which gene flow can be estimated among breeding aggregations of bullfrogs (Rana catesbeiana) at the northern limit of their range in Ontario, Canada, using seven highly polymorphic DNA microsatellite loci. We first identified breeding aggregations that likely share a common history, determined from the pattern of allelic richness, factorial correspondence analysis, and a previously published mtDNA phylogeography, and then tested for regional equilibrium by evaluating the association between pairwise F(ST) and geographic distance. Regional breeding aggregations in eastern Ontario separated by <100 km were determined to be at or near equilibrium. High levels of gene flow were measured using traditional F-statistics and likelihood estimates of Nm. Similarly high levels of recent migration (past one to three generations) were estimated among the breeding aggregations using nonequilibrium methods. We also show that, in many cases, breeding aggregations separated by up to tens of kilometers are not genetically distinct enough to be considered separate genetic populations. These results have important implications both for the identification of independent "populations" and in assessing the effect of scale in detecting patterns of genetic equilibrium and gene flow.     

Isolation by distance in the eastern oyster, Crassostrea virginica, in Chesapeake Bay

Intensive efforts are underway to restore depleted stocks of Crassostrea virginica in Chesapeake Bay. However, the extent of gene flow among local populations, an important force mediating the success of these endeavors, is poorly understood. Spatial and temporal population structures were examined in C. virginica from Chesapeake Bay using eight microsatellite loci. Deficits in heterozygosity relative to Hardy-Weinberg expectations were seen at all loci and were best explained by null alleles. Permutation tests indicated that heterozygote deficiency reduced power in tests of differentiation. Nonetheless, genotypic exact tests demonstrated significant levels of geographic differentiation overall, and a subtle pattern of isolation by distance (IBD) was observed. Comparisons between age classes failed to show differences in genotype frequencies, allelic richness, gene diversity, or differentiation as measured by F(ST), contrary to predictions made by the sweepstakes hypothesis. The IBD pattern could reflect an evolutionary equilibrium established because local gene flow predominates, or be influenced in either direction by recent anthropogenic activities. An evolutionary interpretation appears justified as more parsimonious, implying that local efforts to restore oyster populations will have local demographic payoffs, perhaps at the scale of tributaries or regional subestuaries within Chesapeake Bay.     

Patterns of differentiation in a colour polymorphism and in neutral markers reveal rapid genetic changes in natural damselfly populations

The existence and mode of selection operating on heritable adaptive traits can be inferred by comparing population differentiation in neutral genetic variation between populations (often using F(ST) values) with the corresponding estimates for adaptive traits. Such comparisons indicate if selection acts in a diversifying way between populations, in which case differentiation in selected traits is expected to exceed differentiation in neutral markers [F(ST )(selected) > F(ST )(neutral)], or if negative frequency-dependent selection maintains genetic polymorphisms and pulls populations towards a common stable equilibrium [F(ST) (selected) < F(ST) (neutral)]. Here, we compared F(ST) values for putatively neutral data (obtained using amplified fragment length polymorphism) with estimates of differentiation in morph frequencies in the colour-polymorphic damselfly Ischnura elegans. We found that in the first year (2000), population differentiation in morph frequencies was significantly greater than differentiation in neutral loci, while in 2002 (only 2 years and 2 generations later), population differentiation in morph frequencies had decreased to a level significantly lower than differentiation in neutral loci. Genetic drift as an explanation for population differentiation in morph frequencies could thus be rejected in both years. These results indicate that the type and/or strength of selection on morph frequencies in this system can change substantially between years. We suggest that an approach to a common equilibrium morph frequency across all populations, driven by negative frequency-dependent selection, is the cause of these temporal changes. We conclude that inferences about selection obtained by comparing F(ST) values from neutral and adaptive genetic variation are most useful when spatial and temporal data are available from several populations and time points and when such information is combined with other ecological sources of data.     

Historical introgression and the role of selective vs. neutral processes in structuring nuclear genetic variation (AFLP) in a circumpolar marine fish, the capelin (Mallotus villosus)

The capelin (Mallotus villosus) is a widespread marine fish species for which previous work has identified geographically distinct mtDNA clades, the frontiers of which are well within adult and larval dispersal capabilities. Here, we use AFLPs to test for the presence of nuclear gene flow among clades. In addition, we evaluate genetic structuring within one clade, the Northwest Atlantic (NWA). We found that each of the mtDNA clades corresponds with a unique nuclear DNA genetic cluster. Within the NWA clade, we detected individuals with small but significant amounts of genetic ancestry from other clades, likely due to historical introgression. Further support for historical introgression comes from analyses of variance in locus-specific differentiation, which support introgression between some clades and divergence without gene flow between others. Within the NWA, we identified two genetic clusters that correspond to sites in geographically adjacent areas. However, these clusters differ primarily at 'outlier' loci, and a genetic subdivision (K=2) was not supported by genetic clustering programs using neutral loci. Significant neutral F(ST) differentiation was found only between sites that otherwise differed at outlier loci. Thus, these populations may be in the initial stages of 'isolation by adaptation'. These results suggest strong between-clade reproductive isolation despite opportunities for gene flow and support the hypothesis that selection can contribute to divergence in otherwise 'open' systems.     

Genomic variation in cline shape across a hybrid zone

Hybrid zones are unique biological interfaces that reveal both population level and species level evolutionary processes. A genome‐scale approach to assess gene flow across hybrid zones is vital, and now possible. In Mexican towhees (genus Pipilo), several morphological hybrid gradients exist. We completed a genome survey across one such gradient (9 populations, 140 birds) using mitochondrial DNA, 28 isozyme, and 377 AFLP markers. To assess variation in introgression among loci, cline parameters (i.e., width, center) for the 61 clinally varying loci were estimated and compiled into genomic distributions for tests against three empirical models spanning the range of observed cline shape. No single model accounts for observed variation in cline shape among loci. Numerous backcross individuals near the gradient center confirm a hybrid origin for these populations, contrary to a previous hypothesis based on social mimicry and character displacement. In addition, the observed variation does not bin into well‐defined categories of locus types (e.g., neutral vs. highly selected). Our multi‐locus analysis reveals cross‐genomic variation in selective constraints on gene flow and locus‐specific flexibility in the permeability of the interspecies membrane.     

A multilocus sequence survey in Arabidopsis thaliana reveals a genome-wide departure from a neutral model of DNA sequence polymorphism

The simultaneous analysis of multiple genomic loci is a powerful approach to studying the effects of population history and natural selection on patterns of genetic variation of a species. By surveying nucleotide sequence polymorphism at 334 randomly distributed genomic regions in 12 accessions of Arabidopsis thaliana, we examined whether a standard neutral model of nucleotide sequence polymorphism is consistent with observed data. The average nucleotide diversity was 0.0071 for total sites and 0.0083 for silent sites. Although levels of diversity are variable among loci, no correlation with local recombination rate was observed, but polymorphism levels were correlated for physically linked loci (<250 kb). We found that observed distributions of Tajima's D- and D/D(min)- and of Fu and Li's D-, D*- and F-, F*-statistics differed significantly from the expected distributions under a standard neutral model due to an excess of rare polymorphisms and high variances. Observed and expected distributions of Fay and Wu's H were not different, suggesting that demographic processes and not selection at multiple loci are responsible for the deviation from a neutral model. Maximum-likelihood comparisons of alternative demographic models like logistic population growth, glacial refugia, or past bottlenecks did not produce parameter estimates that were more consistent with observed patterns. However, exclusion of highly polymorphic "outlier loci" resulted in a fit to the logistic growth model. Various tests of neutrality revealed a set of candidate loci that may evolve under selection.