Differential migration from high fitness demes in the shining fungus beetle, Phalacrus substriatus

Abstract.— Using data from three years (1994–1996), I tested whether differential migration occurs from demes of high mean fitness in the shining fungus beetle, Phalacrus substriatus . The results show evidence for differential migration, thus providing evidence from a natural population for a critical demographic assumption of many interdemic selection models. To predict the evolutionary response to interdemic selection through differential migration, the genetic basis of the variation among demes in mean fitness must be known because the observed patterns could also be explained by some demes having an intrinsically favorable habitat. Thus, the importance of differential migration through interdemic selection in natural populations cannot be unequivocally answered without experiments specifically addressing the question of what causes differences in mean fitness among demes.     

THE ROLE OF DEME SIZE,REPRODUCTIVE PATTERNS,AND DISPERSAL IN THE DYNAMICS OF t-LETHAL HAPLOTYPES

The t-lethal haplotypes (t) found in house mouse (Mus musculus) populations are recessive lethals favored by gametic selection whereby male heterozygotes exhibit a non-Mendelian transmission ratio of about 95% t. The expected equilibrium frequency is 0.385; however, empirical values are lower, averaging close to 0.13. We examined the hypothesis that interdemic selection is the cause of the low empirical values by using a deme-structured simulation model that included overlapping generations, a realistic breeding system, differential deme productivity, and a large total population. We found that under some conditions interdemic selection could lower t frequency below 0.13 in the face of immigration rates up to 5%. Low frequencies were correlated with effective deme size (n_e), regardless of whether n_e was changed through changing deme size (n) or through changing the proportion of breeding adults. Earlier workers showed how the first two phases of interdemic selection (random genetic differentiation and mass selection) interacted to reduce the haplotype frequency, but here we show the importance of the third phase (differential productivity of demes) once demes are linked by dispersal. The effect of this phase is not due to the (negative) covariation between deme productivity and haplotype frequency, but occurs when differential deme productivity generates a difference in t frequency between the population of juveniles recruited into their natal deme and the population of juvenile dispersers. This difference was maximized when the average productivity of demes was low, either because few adult females bred at any one time and/or because fecundity was low. Contrary to an earlier prediction, male-biased dispersal also reduced haplotype frequency, and this probably stems from the relative excess of wild-type genotypes among dispersers compared to the deme residents. Another unexpected finding was that the randomly generated excess of heterozygotes (F_IS < 0) found in small demes favored t haplotypes; however, the effect was only seen when the more powerful influence of the third phase of interdemic selection was removed. Simulations of neutral polymorphisms showed that a deme structure giving F_ST ≤ 0.6 is inconsistent with a haplotype frequency below 0.13. Based on current empirical estimates of F_ST (about 0.2), we concluded that immigration rates in the field are too high for interdemic selection alone to cause the observed deficit of lethal haplotypes. One factor that could combine with population structure effects is the observation that the transmission ratio is lowered to around 0.6 in litters produced from postpartum estrus (PPE). Incorporating this factor, we showed that interdemic selection could be effective in lowering the frequency of t below 0.13 when F_ST was above 0.43 even when migration rates were up to 10%. These results suggest that if empirical haplotype and F_ST estimates are accurate, then additional factors such as a lowered fitness of heterozygotes may be involved.     

THE EFFECTIVE SIZE OF A HIERARCHICALLY STRUCTURED POPULATION

A knowledge of the effective size of a population (N_e) is important in understanding its current and future evolutionary potential. Unfortunately, the effective size of a hierarchically structured population is not, in general, equal to the sum of its parts. In particular, the inbreeding structure has a major influence on N_e. Here I link N_e to Wright's hierarchical measures of inbreeding, F_IS and F_ST, for an island-structured population (or metapopulation) of size N_T. The influence of F_ST depends strongly on the degree to which island productivity is regulated. In the absence of local regulation (the interdemic model), interdemic genetic drift reduces N_e. When such drift is combined with local inbreeding under otherwise ideal conditions, the effects of F_IS and F_ST are identical: increasing inbreeding either within or between islands reduces N_e, with N_e = N_T/[(1 + F_IS)(1 + F_ST) ? 2F_ISF_ST]. However, if islands are all equally productive because of local density regulation (the traditional island model), then N_e = N_T/[(1 + F_IS)(1 –F_ST)] and the effect of F_ST is reversed. Under the interdemic model, random variation in the habitat quality (and hence productivity) of islands act to markedly decrease N_e. This variation has no effect under the island model because, by definition, all islands are equally productive. Even when no permanent island structure exists, spatial differences in habitat quality can significantly increase the overall variance in reproductive success of both males and females and hence lower N_e. Each of these basic results holds when other nonideal factors are added to the model. These factors, deviations from a 1:1 sex ratio, greater than Poisson variance in female reproductive success, and variation in male mating success due to polygynous mating systems, all act to lower N_e. The effects of male and female variance on N_e have important differences because only females affect island productivity. Finally, it is noted that to use these relationships, F_IS and F_ST must be estimated according to Wright's definition (and corrected to have a zero expectation under the null model). A commonly used partitioning (θ, θ_g) can be biased if either island size or the number of islands is small.     

THE INFLUENCE OF DISPERSAL PATTERNS AND MATING SYSTEM ON GENETIC DIFFERENTIATION WITHIN AND BETWEEN POPULATIONS OF THE RED HOWLER MONKEY (ALOUATTA SENICULUS)

The relationship between social structure and partitioning of genetic variance was examined in two red howler monkey populations (W and G) in Venezuela, one of which (G) was undergoing rapid growth through colonization by new troops. Rates and patterns of gene flow had been determined through radiotelemetry and direct observation data on solitary migrants, and 10 years of troop censusing. Standard electrophoresis techniques were used to examine 29 loci in blood samples taken from 137 of the study animals. Analysis of genetic variance demonstrated: (1) a significantly high level of genetic variation among troops within populations (F_ST = 0.225 for W and 0.142 for G), and (2) a significant excess of heterozygosity within troops relative to expected (F_IS = -0.136 for W and -0.064 for G), despite relatively high levels of observed and inferred inbreeding in W. Differences between the populations in F_ST values conformed to those predicted based on differences in colonization rate. Comparison of partitioning of genetic variance among different genealogical subsets of troops demonstrated that the pattern of genetic differentiation observed among troops within populations was promoted by an essentially single-male harem breeding structure, a very low rate of random exchange of breeding males among troops, and a high degree of relatedness among troop females. Between-troop genetic differentiation (F_ST) was thereby increased relative to that expected from other types of social organization, while the correlation between uniting gametes within troops (F_IS) was decreased. Genetic differentiation between populations (2%) corresponded to that predicted from migration rates. Such a mosaic of genetic variation, combined with differences in reproductive success observed among troops and a high troop failure rate, create conditions in which interdemic selection could result in more rapid spread of advantageous gene combinations than would be expected in a panmictic population, particularly in a colonizing situation in which the founder population is small.     

Is local selection so widespread in river organisms? Fractal geometry of river networks leads to high bias in outlier detection

Identifying local adaptation is crucial in conservation biology to define ecotypes and establish management guidelines. Local adaptation is often inferred from the detection of loci showing a high differentiation between populations, the so‐called F_ST outliers. Methods of detection of loci under selection are reputed to be robust in most spatial population models. However, using simulations we showed that F_ST outlier tests provided a high rate of false‐positives (up to 60%) in fractal environments such as river networks. Surprisingly, the number of sampled demes was correlated with parameters of population genetic structure, such as the variance of F_STs, and hence strongly influenced the rate of outliers. This unappreciated property of river networks therefore needs to be accounted for in genetic studies on adaptation and conservation of river organisms.     

HIERARCHICAL ORGANIZATION OF GENETIC VARIATION IN THE SAILFIN MOLLY,POECILIA LATIPINNA (PISCES: POECILIIDAE)

Twenty-two percent of all allozyme variation documented in the sailfin molly (Poecilia latipinna) was attributable to regional differences, while only 3% was attributable to differences among demes within regions. Of the variation documented in a given region, 6–12% was attributable to variation among demes. Cluster analyses supported these conclusions quantitatively. Spatial-autocorrelation analyses offered more explicit support: demes separated by increasingly greater distances were increasingly dissimilar. Analyses using F statistics and rare alleles suggest “effective gene flow rates” (the product of effective population size and gene flow rate) of approximately 4, a level more than sufficient to prevent local independence of gene-frequency dynamics. These results, taken together, suggest that mollies do not have a population structure conducive to the operation of Wright's shifting-balance process and make the striking patterns of interdemic variation in body size and sexual behavior observed in this species all the more interesting.     

TEMPORAL FLUCTUATIONS IN DEMOGRAPHIC PARAMETERS AND THE GENETIC VARIANCE AMONG POPULATIONS

Contrary to assumptions commonly made in the study of population genetics, the demographic properties of many populations are not always constant. Important characteristics of populations such as migration rate and population size may vary in time and space. Moreover, local populations often come and go; the rate of extinction and the properties of colonization may also vary. In this paper, the approach to equilibrium following a disturbance in the genetic variance among populations is described. The rate of migration is shown to be critical in determining the extent to which extinction and recolonization affects genetic differentiation. Perturbations and variations through time and space in demographic parameters such as population size and migration rate are shown to be important in determining the partitioning of genetic variance. Equations are given to predict the average through time of genetic differentiation among populations in the event of a single disturbance or in constant fluctuations in the pertinent demographic parameters. In general, these fluctuations increase the F_ST of a species. Spatial demographic variation affects F_STmuch more than temporal variation. These demographic properties make some species unsuitable for the empirical analysis of migration with indirect genetic measures. Demographic instability may play a large role in the evolution of genetic variation.     

The Effects of Background and Interference Selection on Patterns of Genetic Variation in Subdivided Populations

It is well known that most new mutations that affect fitness exert deleterious effects and that natural populations are often composed of subpopulations (demes) connected by gene flow. To gain a better understanding of the joint effects of purifying selection and population structure, we focus on a scenario where an ancestral population splits into multiple demes and study neutral diversity patterns in regions linked to selected sites. In the background selection regime of strong selection, we first derive analytic equations for pairwise coalescent times and F_ST as a function of time after the ancestral population splits into two demes and then construct a flexible coalescent simulator that can generate samples under complex models such as those involving multiple demes or nonconservative migration. We have carried out extensive forward simulations to show that the new methods can accurately predict diversity patterns both in the nonequilibrium phase following the split of the ancestral population and in the equilibrium between mutation, migration, drift, and selection. In the interference selection regime of many tightly linked selected sites, forward simulations provide evidence that neutral diversity patterns obtained from both the nonequilibrium and equilibrium phases may be virtually indistinguishable for models that have identical variance in fitness, but are nonetheless different with respect to the number of selected sites and the strength of purifying selection. This equivalence in neutral diversity patterns suggests that data collected from subdivided populations may have limited power for differentiating among the selective pressures to which closely linked selected sites are subject.     

NEUTRAL THEORY OF QUANTITATIVE GENETIC VARIANCE IN AN ISLAND MODEL WITH LOCAL EXTINCTION AND COLONIZATION

Additive genetic variance maintained by mutation in a selectively neutral quantitative character is analyzed for an ideal population distributed on n islands, each with local effective size N, that exchange migrants at a small rate, m. In a stable population structure, the expected genetic variance maintained within islands is identical to that in a panmictic population of the same total size, regardless of the migration rate (m > 0). This result contrasts with Wright's classical conclusion, based on inbreeding coefficients, that at least one immigrant per island every other generation (Nm > ½) is necessary for the genetic variance within local populations to approach that under panmixia. The expected genetic variance maintained among islands is inversely proportional to m and increases with the number of islands, but is independent of N. Local extinction and colonization diminish the genetic variance maintained within islands by reducing the effective size of island populations through the founder effect, although the expected genetic variance within islands is nearly as large as that in a panmictic population of the same total effective size. If the founders of new colonies originate from more than one island, rates of local extinction and colonization larger than about twice the migration rate will substantially reduce the genetic variance maintained among islands. These results indicate the importance of mutation and migration in maintaining quantitative genetic variance within small local populations.     

Determinants of Genetic Structure in a Nonequilibrium Metapopulation of the Plant Silene latifolia

Population genetic differentiation will be influenced by the demographic history of populations, opportunities for migration among neighboring demes and founder effects associated with repeated extinction and recolonization. In natural populations, these factors are expected to interact with each other and their magnitudes will vary depending on the spatial distribution and age structure of local demes. Although each of these effects has been individually identified as important in structuring genetic variance, their relative magnitude is seldom estimated in nature. We conducted a population genetic analysis in a metapopulation of the angiosperm, Silene latifolia, from which we had more than 20 years of data on the spatial distribution, demographic history, and extinction and colonization of demes. We used hierarchical Bayesian methods to disentangle which features of the populations contributed to among population variation in allele frequencies, including the magnitude and direction of their effects. We show that population age, long-term size and degree of connectivity all combine to affect the distribution of genetic variance; small, recently-founded, isolated populations contributed most to increase F _ST in the metapopulation. However, the effects of population size and population age are best understood as being modulated through the effects of connectivity to other extant populations, i.e. F _ST diminishes as populations age, but at a rate that depends how isolated the population is. These spatial and temporal correlates of population structure give insight into how migration, founder effect and within-deme genetic drift have combined to enhance and restrict genetic divergence in a natural metapopulation.     

Adaptive divergence despite strong genetic drift: genomic analysis of the evolutionary mechanisms causing genetic differentiation in the island fox (Urocyon littoralis)

The evolutionary mechanisms generating the tremendous biodiversity of islands have long fascinated evolutionary biologists. Genetic drift and divergent selection are predicted to be strong on islands and both could drive population divergence and speciation. Alternatively, strong genetic drift may preclude adaptation. We conducted a genomic analysis to test the roles of genetic drift and divergent selection in causing genetic differentiation among populations of the island fox (Urocyon littoralis). This species consists of six subspecies, each of which occupies a different California Channel Island. Analysis of 5293 SNP loci generated using Restriction‐site Associated DNA (RAD) sequencing found support for genetic drift as the dominant evolutionary mechanism driving population divergence among island fox populations. In particular, populations had exceptionally low genetic variation, small N_e (range = 2.1–89.7; median = 19.4), and significant genetic signatures of bottlenecks. Moreover, islands with the lowest genetic variation (and, by inference, the strongest historical genetic drift) were most genetically differentiated from mainland grey foxes, and vice versa, indicating genetic drift drives genome‐wide divergence. Nonetheless, outlier tests identified 3.6–6.6% of loci as high F_ST outliers, suggesting that despite strong genetic drift, divergent selection contributes to population divergence. Patterns of similarity among populations based on high F_ST outliers mirrored patterns based on morphology, providing additional evidence that outliers reflect adaptive divergence. Extremely low genetic variation and small N_e in some island fox populations, particularly on San Nicolas Island, suggest that they may be vulnerable to fixation of deleterious alleles, decreased fitness and reduced adaptive potential.     

A SIMULATION OF WRIGHT'S SHIFTING-BALANCE PROCESS: MIGRATION AND THE THREE PHASES

Wright partitioned the shifting-balance process into three phases. Phase one is the shift of a deme within a population to the domain of a higher adaptive peak from that of the historical peak. Phase two is mass selection within a deme towards that higher peak. Phase three is the conversion of additional demes to the higher peak. The migration rate between demes is critical for the existence of phases one and three. Phase one requires small effective population sizes, hence low migration rates. Phase three is optimal under high migration rates that spread the most-fit genotype from deme to deme. Thus, a population-wide peak shift requires intermediate levels of migration. By altering the rates of phases one and three, migration affects the predominant direction of mass selection within a population. This study examines the degree to which migration, through its effects on phases one and three, determines the probability of a simulated population arriving at its genotypic optimum after 12,000 generations. These simulations reveal that there is a range of migration rates for which an entire population might be expected to shift to a higher peak. Below m = 0.001 peak shifts occur frequently (phases I and II) but are not successfully exported out of subpopulations (phase III), and above 0.01 peak shifts within demes (phase I and II), required to initiate phase III, become increasingly uncommon. Because it is unlikely that real populations will have uniform migration rates from generation to generation, the probable effects of varying migration rates on broadening the range of conditions producing peak shifts are discussed.     

Generation of a neutral FST baseline for testing local adaptation on gill raker number within and between European whitefish ecotypes in the Baltic Sea basin

Divergent selection at ecologically important traits is thought to be a major factor driving phenotypic differentiation between populations. To elucidate the role of different evolutionary processes shaping the variation in gill raker number of European whitefish (Coregonus lavaretus sensu lato) in the Baltic Sea basin, we assessed the relationships between genetic and phenotypic variation among and within three whitefish ecotypes (sea spawners, river spawners and lake spawners). To generate expected neutral distribution of F_ST and to evaluate whether highly variable microsatellite loci resulted in deflated F_ST estimates compared to less variable markers, we performed population genetic simulations under finite island and hierarchical island models. The genetic divergence observed among (F_CT = 0.010) and within (F_ST = 0.014–0.041) ecotypes was rather low. The divergence in gill raker number, however, was substantially higher between sea and river spawners compared to observed microsatellite data and simulated neutral baseline (P_CT > F_CT). This suggests that the differences in gill raker number between sea and river spawners are likely driven by divergent natural selection. We also found strong support for divergent selection on gill raker number among different populations of sea spawners (P_ST > F_ST), most likely caused by highly variable habitat use and diverse diet. The putative role of divergent selection within lake spawners initially inferred from empirical microsatellite data was not supported by simulated F_ST distributions. This work provides a first formal test of divergent selection on gill raker number in Baltic whitefish, and demonstrates the usefulness of population genetic simulations to generate informative neutral baselines for P_ST–F_ST analyses helping to disentangle the effects of stochastic evolutionary processes from natural selection.     

GENE FLOW IN GROUND BEETLES (COLEOPTERA: CARABIDAE) OF DIFFERING HABITAT PREFERENCE AND FLIGHT-WING DEVELOPMENT

Estimates of gene flow vary 100-fold among five carabid species, ranging from the winged lowland subtropical Agonum elongatulum to the flightless montane temperate Platynus angustatus. Results based on Wright's (1943) F_ST method, and Slatkin's (1981) graphical and (1985a) private-allele methods are concordant. Genetic heterogeneity, measured by Wright's F_ST, is not correlated with degree of flight-wing development; one fully winged species exhibits heterogeneity of the same order as a vestigially winged species. Genetic heterogeneity is positively correlated with the average elevation of collection sites for these species. Lower levels of gene flow associated with greater genetic subdivision may occur in upland areas because of habitat fragmentation (due to topographic diversity) and habitat persistence (leading to a lower extinction rate for populations). In at least one species, the distribution of stable infraspecific polymorphisms indicates that the high estimate of present-day gene flow is likely to be due to historical gene flow and not to present-day conditions.     

Population structure and interdemic selection in the cooperative spider Anelosimus eximius

The degree of relatedness among interacting individuals helps determine the fitness consequences of particular behaviors, whereas the partitioning (and amount) of genetic variation among and within groups controls the level at which selection will act most effectively. Three criteria are considered necessary for selection to act at the group or interdemic level: high rate of group initiation/extinction; differential survival and reproduction among groups; and highly subdivided population structure. The first two criteria have been demonstrated by earlier studies of Anelosimus eximius colonies. This study employs hierarchical analysis of allozyme polymorphisms to demonstrate the third criterion, subdivided population structure. Anelosimus eximius were collected from Suriname, Panama, Ecuador, Peru and Trinidad. Seven of 40 scorable enzyme loci revealed variation; 4 of these were polymorphic within colonies or regions. Expected heterozygosities were low, ranging from 0 (Ecuador, Peru) to ～0.03 (Suriname). For each polymorphic locus, hierarchical F -statistics were used to partition overall genetic variation into among-region (or among-population; F _rt ), among-colony ( F _sr ), and within-colony ( F _is ) components. Samples from Suriname (43 colonies, 4 local populations) were the most informative; lack of scorable variation limited the inferences that could be drawn from other regions. A. eximius colonies are highly inbred: negative estimates of F _is imply very small effective colony sizes (～6.5 for Suriname samples). By contrast, estimates of F _sr were very high: the mean for Suriname samples was 0.890, indicating neglibible gene flow among established colonies. Inbreeding within colonies, and genetic differentiation among colonies are consistent with demographic and behavioral observations of A. eximius . We suggest that interdemic selection is probable in this species and other cooperative spiders with this type of social system, and that mutual tolerance and absence of nest-mate recognition, as well as female-biased sex ratios, may have arisen by interdemic selection.     

MAXIMUM-LIKELIHOOD ESTIMATION OF POPULATION DIVERGENCE TIMES AND POPULATION PHYLOGENY IN MODELS WITHOUT MUTATION

In this paper we present a method for estimating population divergence times by maximum likelihood in models without mutation. The maximum-likelihood estimator is compared to a commonly applied estimator based on Wright's F_ST statistic. Simulations suggest that the maximum-likelihood estimator is less biased and has a lower variance than the F_ST-based estimator. The maximum-likelihood estimator provides a statistical framework for the analysis of population history given genetic data. We demonstrate how maximum-likelihood estimates of the branching pattern of divergence of multiple populations may be obtained. We also describe how the method may be applied to test hypotheses such as whether populations have maintained equal population sizes. We illustrate the method by applying it to two previously published sets of human restriction fragment length polymorphism (RFLP) data.     

CORRELATION OF PAIRWISE GENETIC AND GEOGRAPHIC DISTANCE MEASURES: INFERRING THE RELATIVE INFLUENCES OF GENE FLOW AND DRIFT ON THE DISTRIBUTION OF GENETIC VARIABILITY

Attempts to relate estimates of regional F_ST to gene flow and drift via Wright's (1931) equation F_ST ≈ 1/ (4Nm + 1) are often inappropriate because most natural sets of populations probably are not at equilibrium (McCauley 1993), as assumed by the island model upon which the equation is based, or ineffective because the influences of gene flow and drift are confounded in the product Nm. Evaluations of the association between genetic (F_ST) and geographic distances separating all pairwise populations combinations in a region allows one to test for regional equilibrium, to evaluate the relative influences of gene flow and drift on population structure both within and between regions, and to visualize the behavior of the association across all degrees of geographic separation. Tests of the model using microsatellite data from 51 populations of eastern collared lizards (Crotaphytus collaris collaris) collected from four distinct geographical regions gave results highly consistent with predicted patterns of association based on regional differences in various historical and ecological factors that affect the amount of drift and gene flow. The model provides a prerequisite for and an alternative to regional F_ST analyses, which often simply assume regional equilibrium, thus potentially leading to erroneous and misleading inferences regarding regional population structure.     

F(ST) and Q(ST) under neutrality

A commonly used test for natural selection has been to compare population differentiation for neutral molecular loci estimated by F_ST and for the additive genetic component of quantitative traits estimated by Q_ST. Past analytical and empirical studies have led to the conclusion that when averaged over replicate evolutionary histories, Q_ST = F_ST under neutrality. We used analytical and simulation techniques to study the impact of stochastic fluctuation among replicate outcomes of an evolutionary process, or the evolutionary variance, of Q_ST and F_ST for a neutral quantitative trait determined by n unlinked diallelic loci with additive gene action. We studied analytical models of two scenarios. In one, a pair of demes has recently been formed through subdivision of a panmictic population; in the other, a pair of demes has been evolving in allopatry for a long time. A rigorous analysis of these two models showed that in general, it is not necessarily true that mean Q_ST = F_ST (across evolutionary replicates) for a neutral, additive quantitative trait. In addition, we used finite-island model simulations to show there is a strong positive correlation between Q_ST and the difference Q_ST − F_ST because the evolutionary variance of Q_ST is much larger than that of F_ST. If traits with relatively large Q_ST values are preferentially sampled for study, the difference between Q_ST and F_ST will also be large and positive because of this correlation. Many recent studies have used tests of the null hypothesis Q_ST = F_ST to identify diversifying or uniform selection among subpopulations for quantitative traits. Our findings suggest that the distributions of Q_ST and F_ST under the null hypothesis of neutrality will depend on species-specific biology such as the number of subpopulations and the history of subpopulation divergence. In addition, the manner in which researchers select quantitative traits for study may introduce bias into the tests. As a result, researchers must be cautious before concluding that selection is occurring when Q_ST ≠ F_ST.     

A HIERARCHICAL ANALYSIS OF GENETIC DIFFERENTIATION IN A MONTANE LEAF BEETLE CHRYSOMELA AENEICOLLIS (COLEOPTERA: CHRYSOMELIDAE)

Herbivorous insects that use the same host plants as larvae and adults can have a subdivided population structure that corresponds to the distribution of their hosts. Having a subdivided population structure favors local adaptation of subpopulations to small-scale environmental differences and it may promote their genetic divergence. In this paper, I present the results of a hierarchical study of population structure in a montane willow leaf beetle, Chrysomela aeneicollis (Coleoptera: Chrysomelidae). This species spends its entire life associated with the larval host (Salix spp.), which occurs in patches along high-elevation streams and in montane bogs. I analyzed the genetic differentiation of C. aeneicollis populations along three drainages in the Sierra Nevada mountains of California at five enzyme loci: ak-1, idh-2, mpi-1, pgi-1, and pgm-1, using recent modifications of Wright's F-statistics. My results demonstrated significant differentiation (F_ST = 0.043) among drainages that are less than 40 kilometers apart. One locus, pgi-1, showed much greater differentiation than the other four (F_ST = 0.412), suggesting that it is under natural selection. C. aeneicollis populations were also subdivided within drainages, with significant differentiation 1) among patches of willows (spanning less than three kilometers) and 2) in some cases, among trees within a willow patch. My results demonstrate that this species has the capacity to adapt to local environmental variation at small spatial scales.