Evidence of spatial genetic structure in a California bunchgrass population

We investigated the scale of genetic variation of purple needlegrass (Nassella pulchra), a species commonly used in California for grassland restoration. Common garden and field data revealed evidence of genetic differentiation between two intermixed microhabitats characterized by differences in soil depth and community composition. We assessed the genetic variation within a single population using randomly amplified polymorphic DNA (RAPD) data collected from clusters of five individuals in 40 locations. We found no evidence for genetic structure at the whole population level. At smaller spatial scales, however, we found strong evidence that genetic subdivision of the population occurs at the level of the maternal neighborhood. We suggest that the interaction between widespread pollen dispersal and restricted seed dispersal may be the primary factor generating these results; panmictic pollen dispersal will make detection of genetic patterning difficult at larger spatial scales while limited seed dispersal will generate local genetic structure. As a result, the detection of population genetic structure will depend on the spatial scale of analysis. Local selection gradients related to topography and soil depth are also likely to play a role in structuring local genetic variation. Since N. pulchra is widely used in California in grassland and woodland habitat restoration, we suggest that, as a general rule, care should be exercised in transferring germplasm for the purposes of conservation when little is known about the within-population genetic subdivision of a plant species. Received: 23 December 1996 / Accepted: 20 May 1997     

Fine-scale genetically based differentiation of life-history traits in the perennial shrub Lupinus arboreus

Across large spatial scales, plants often exhibit genetically based differentiation in traits that allow adaptation to local sites. At smaller spatial scales, sharp boundaries between edaphic conditions also can create strong gradients in selection that counteract gene flow and result in local adaptation. Few studies, however, have examined the degree to which continuous populations of perennial plants exhibit genetically based differentiation in life-history traits over small spatial scales. We quantified the degree of genetically based differentiation in adaptive traits among bush lupine (Lupinus arboreus) from nearby dune and grassland sites (sites separated by < 0.75 km) that formed part of a larger continuous population of L. arboreus. We also investigated the spatial genetic structure of bush lupine by examining how genetic structure differed between seeds and juvenile plants that were less than two years old. We calculated F-statistics from gel electrophoresis of 10 polymorphic loci. We then used these values to infer levels of gene flow. To examine differentiation in adaptive traits, we created full-sibling/half-sibling families of lupine within each area and established reciprocal common gardens at each site. Across two years, we measured canopy volume, flowering time, seed set, and mortality of progeny planted in each garden. Spatial genetic structure among seeds was virtually nonexistent (F(ST) = 0.002), suggesting that gene flow between the three areas could be quite high. However, genetic structure increased 20-fold among juvenile plants (F(ST) = 0.041). We found strong evidence for fine-scale genetically based differentiation and local adaptation in adaptive traits such as plant size, flowering phenology, fecundity, and mortality. Thus, it is likely that strong but differing selection regimes within each area drive spatial differentiation in lupine life-history traits.     

Cuckoo parasitism and productivity in different magpie subpopulations predict frequencies of the 457bp allele: a mosaic of coevolution at a small geographic scale

The level of defense against great spotted cuckoo (Clamator glandarius) parasitism in different European populations of magpie (Pica pica) depends on selection pressures due to parasitism and gene flow between populations, which suggests the existence of coevolutionary hot spots within a European metapopulation. A mosaic of coevolution is theoretically possible at small geographical scales and with strong gene flow, because, among other reasons, plots may differ in productivity (i.e., reproductive success of hosts in the absence of parasitism) and defensive genotypes theoretically should be more common in plots of high productivity. Here, we tested this prediction by exploring the relationship between parasitism rate, level of defense against parasitism (estimated as both rejection rate and the frequency of the 457bp microsatellite allele associated with foreign egg rejection in magpies), and some variables related to the productivity (average laying date, clutch size, and number of hatchlings per nest) of magpies breeding in different subpopulations. We found that both estimates of defensive ability (egg rejection rate and frequency of the 457bp allele) covaried significantly with between-plot differences in probability of parasitism, laying date, and number of hatchlings per nest. Moreover, the parasitism rate was larger in more productive plots. These results confirm the existence of a mosaic of coevolution at a very local geographical scale, and the association between laying date and number of hatchlings with variables related to defensive ability and the selection pressure arising from parasitism supports the prediction of coevolutionary gradients in relation to host productivity.     

Isolation by Distance in a Quantitative Trait

Random genetic drift in a quantitative character is modeled for a population with a continuous spatial distribution in an infinite habitat of one or two dimensions. The analysis extends Wright's concept of neighborhood size to spatially autocorrelated sampling variation in the expected phenotype at different locations. Weak stabilizing selection is assumed to operate toward the same optimum phenotype in every locality, and the distribution of dispersal distances from parent to offspring is a (radially) symmetric function. The equilibrium pattern of geographic variation in the expected local phenotype depends on the neighborhood size, the genetic variance within neighborhoods, and the strength of selection, but is nearly independent of the form of the dispersal function. With all else equal, geographic variance is smaller in a two-dimensional habitat than in one dimension, and the covariance between expected local phenotypes decreases more rapidly with the distance separating them in two dimensions than in one. The equilibrium geographic variance is less than the phenotypic variance within localities, unless the neighborhood size is small and selection is extremely weak, especially in two dimensions. Nevertheless, dispersal of geographic variance created by random genetic drift is an important mechanism maintaining genetic variance within local populations. For a Gaussian dispersal function it is shown that, even with a small neighborhood size, a population in a two-dimensional habitat can maintain within neighborhoods most of the genetic variance that would occur in an infinite panmictic population.     

Fine‐grained adaptive divergence in an amphibian: genetic basis of phenotypic divergence and the role of nonrandom gene flow in restricting effective migration among wetlands

Adaptive ecological differentiation among sympatric populations is promoted by environmental heterogeneity, strong local selection and restricted gene flow. High gene flow, on the other hand, is expected to homogenize genetic variation among populations and therefore prevent local adaptation. Understanding how local adaptation can persist at the spatial scale at which gene flow occurs has remained an elusive goal, especially for wild vertebrate populations. Here, we explore the roles of natural selection and nonrandom gene flow (isolation by breeding time and habitat choice) in restricting effective migration among local populations and promoting generalized genetic barriers to neutral gene flow. We examined these processes in a network of 17 breeding ponds of the moor frog Rana arvalis, by combining environmental field data, a common garden experiment and data on variation in neutral microsatellite loci and in a thyroid hormone receptor (TRβ) gene putatively under selection. We illustrate the connection between genotype, phenotype and habitat variation and demonstrate that the strong differences in larval life history traits observed in the common garden experiment can result from adaptation to local pond characteristics. Remarkably, we found that haplotype variation in the TRβ gene contributes to variation in larval development time and growth rate, indicating that polymorphism in the TRβ gene is linked with the phenotypic variation among the environments. Genetic distance in neutral markers was correlated with differences in breeding time and environmental differences among the ponds, but not with geographical distance. These results demonstrate that while our study area did not exceed the scale of gene flow, ecological barriers constrained gene flow among contrasting habitats. Our results highlight the roles of strong selection and nonrandom gene flow created by phenological variation and, possibly, habitat preferences, which together maintain genetic and phenotypic divergence at a fine‐grained spatial scale.     

Detecting clinal and balanced selection using spatial autocorrelation analysis under kin-structured migration

Recently spatial autocorrelation has been employed to infer microevolutionary processes from patterns of genetic variation. In theory, different processes should show characteristic signature correlograms; e. g., clinal selection should produce correlograms decreasing from positive to negative autocorrelation, whereas uniform balanced selection should lead to no spatial autocorrelation. The ability of a statistical method such as spatial autocorrelation analysis to distinguish between these selective regimes or even to detect departures from neutrality is dependent on the strength of the evolutionary force and the population structure. Weak selection or migration will not be apparent against the expected background of stochastic noise. Moreover, the population structure may generate sufficient stochastic variation such that even strong evolutionary forces may fail to be detected. This study uses computer simulation to examine the effects of kin-structured migration and three different selective regimes on the shape of spatial correlograms to assess the ability of this technique to detect different microevolutionary processes. Genetic variation among 8 loci is simulated in a linear set of 25 artificial populations. Kin-structured stepping-stone migration among adjacent populations is modeled; directional, balanced, and clinal selection, as well as neutral loci are considered. These experiments show that strong selection produces correlograms of the predicted shape. However, with an anthropologically reasonable population structure, considerable stochastic variation among correlograms for different alleles may still exist. This suggests the need for caution in inferring genetic process from spatial patterns. © 1994 Wiley-Liss, Inc.     

Local and regional scale habitat heterogeneity contribute to genetic adaptation in a commercially important marine mollusc (Haliotis rubra) from southeastern Australia

Characterising adaptive genetic divergence among conspecific populations is often achieved by studying genetic variation across defined environmental gradients. In marine systems this is challenging due to a paucity of information on habitat heterogeneity at local and regional scales and a dependency on sampling regimes that are typically limited to broad longitudinal and latitudinal environmental gradients. As a result, the spatial scales at which selection processes operate and the environmental factors that contribute to genetic adaptation in marine systems are likely to be unclear. In this study we explore patterns of adaptive genetic structuring in a commercially‐ harvested abalone species (Haliotis rubra) from southeastern Australia, using a panel of genome‐wide SNP markers (5,239 SNPs), and a sampling regime informed by marine LiDAR bathymetric imagery and 20‐year hindcasted oceanographic models. Despite a lack of overall genetic structure across the sampling distribution, significant genotype associations with heterogeneous habitat features were observed at local and regional spatial scales, including associations with wave energy, ocean current, sea surface temperature, and geology. These findings provide insights into the potential resilience of the species to changing marine climates and the role of migration and selection on recruitment processes, with implications for conservation and fisheries management. This study points to the spatial scales at which selection processes operate in marine systems and highlights the benefits of geospatially‐informed sampling regimes for overcoming limitations associated with marine population genomic research.     

Genome-wide association to fine-scale ecological heterogeneity within a continuous population of Biscutella laevigata (Brassicaceae)

Gene flow, drift and selection can be detected through different signatures across the genome and the landscape. Genetic discontinuities along with their correlation to environmental features can be used to tease out isolation-by-distance and isolation-by-time from processes related to selection. Using spatial statistics (spatial autocorrelation methods, canonical correspondence analysis and partial Mantel tests) dealing with genome-wide amplified fragment length polymorphism (AFLP) under unlikely Hardy-Weinberg assumptions, this study investigates 124 individuals within a continuous population of the autopolyploid Biscutella laevigata (Brassicaceae). Fine-scale spatial genetic structure was strong and the mosaic-like distribution of AFLP genotypes was consistently associated with habitat factors, even when controlled for geographical distances. The use of multivariate analyses enabled separation of the factors responsible for the repartition of the genetic variance and revealed a composite effect of isolation by distance, phenological divergence and local adaptation to habitats characterised by different solar radiation regimes. These results suggest that the immigrant inviability barrier facilitated the maintenance of adapted subpopulations to distinct environmental conditions at the local scale.     

Signatures of environmental genetic adaptation pinpoint pathogens as the main selective pressure through human evolution

Previous genome-wide scans of positive natural selection in humans have identified a number of non-neutrally evolving genes that play important roles in skin pigmentation, metabolism, or immune function. Recent studies have also shown that a genome-wide pattern of local adaptation can be detected by identifying correlations between patterns of allele frequencies and environmental variables. Despite these observations, the degree to which natural selection is primarily driven by adaptation to local environments, and the role of pathogens or other ecological factors as selective agents, is still under debate. To address this issue, we correlated the spatial allele frequency distribution of a large sample of SNPs from 55 distinct human populations to a set of environmental factors that describe local geographical features such as climate, diet regimes, and pathogen loads. In concordance with previous studies, we detected a significant enrichment of genic SNPs, and particularly non-synonymous SNPs associated with local adaptation. Furthermore, we show that the diversity of the local pathogenic environment is the predominant driver of local adaptation, and that climate, at least as measured here, only plays a relatively minor role. While background demography by far makes the strongest contribution in explaining the genetic variance among populations, we detected about 100 genes which show an unexpectedly strong correlation between allele frequencies and pathogenic environment, after correcting for demography. Conversely, for diet regimes and climatic conditions, no genes show a similar correlation between the environmental factor and allele frequencies. This result is validated using low-coverage sequencing data for multiple populations. Among the loci targeted by pathogen-driven selection, we found an enrichment of genes associated to autoimmune diseases, such as celiac disease, type 1 diabetes, and multiples sclerosis, which lends credence to the hypothesis that some susceptibility alleles for autoimmune diseases may be maintained in human population due to past selective processes.     

COMPARING STRENGTHS OF DIRECTIONAL SELECTION: HOW STRONG IS STRONG?

The fundamental equation in evolutionary quantitative genetics, the Lande equation, describes the response to directional selection as a product of the additive genetic variance and the selection gradient of trait value on relative fitness. Comparisons of both genetic variances and selection gradients across traits or populations require standardization, as both are scale dependent. The Lande equation can be standardized in two ways. Standardizing by the variance of the selected trait yields the response in units of standard deviation as the product of the heritability and the variance-standardized selection gradient. This standardization conflates selection and variation because the phenotypic variance is a function of the genetic variance. Alternatively, one can standardize the Lande equation using the trait mean, yielding the proportional response to selection as the product of the squared coefficient of additive genetic variance and the mean-standardized selection gradient. Mean-standardized selection gradients are particularly useful for summarizing the strength of selection because the mean-standardized gradient for fitness itself is one, a convenient benchmark for strong selection. We review published estimates of directional selection in natural populations using mean-standardized selection gradients. Only 38 published studies provided all the necessary information for calculation of mean-standardized gradients. The median absolute value of multivariate mean-standardized gradients shows that selection is on average 54% as strong as selection on fitness. Correcting for the upward bias introduced by taking absolute values lowers the median to 31%, still very strong selection. Such large estimates clearly cannot be representative of selection on all traits. Some possible sources of overestimation of the strength of selection include confounding environmental and genotypic effects on fitness, the use of fitness components as proxies for fitness, and biases in publication or choice of traits to study.     

Genomic signature of natural and anthropogenic stress in wild populations of the waterflea Daphnia magna: validation in space, time and experimental evolution

Natural populations are confronted with multiple selection pressures resulting in a mosaic of environmental stressors at the landscape level. Identifying the genetic underpinning of adaptation to these complex selection environments and assigning causes of natural selection within multidimensional selection regimes in the wild is challenging. The water flea Daphnia is a renowned ecological model system with its well-documented ecology, the possibility to analyse subfossil dormant egg banks and the short generation time allowing an experimental evolution approach. Capitalizing on the strengths of this model system, we here link candidate genome regions to three selection pressures, known to induce micro-evolutionary responses in Daphnia magna: fish predation, parasitism and land use. Using a genome scan approach in space, time and experimental evolution trials, we provide solid evidence of selection at the genome level under well-characterized environmental gradients in the wild and identify candidate genes linked to the three environmental stressors. Our study reveals differential selection at the genome level in Daphnia populations and provides evidence for repeatable patterns of local adaptation in a geographic mosaic of environmental stressors fuelled by standing genetic variation. Our results imply high evolutionary potential of local populations, which is relevant to understand the dynamics of trait changes in natural populations and their impact on community and ecosystem responses through eco-evolutionary feedbacks.     

On humans and wildlife in Mediterranean islands

Aim To investigate the effects of human‐induced landscape changes in Mediterranean islands on the ecological and evolutionary responses of bird communities and populations. The combination of mass extinction of large mammals and massive deforestation by humans was hypothesized to produce new selection regimes to which organisms were likely to respond. Habitat selection and niche breadth have been investigated at the scale of species, and phenotypic variation at the scale of local populations. Location The study was carried out along habitat gradients and in habitat mosaics at different spatial scales on the island of Corsica and in areas of similar size and structure in continental France. Methods Two sets of gradients have been used for investigating habitat selection and niche breadth: gradients of altitude, and gradients of vegetation structure. Population studies focused on the blue tit, Cyanistes caeruleus. Large samples of breeding attempts by this species in 10 habitats provided detailed data on phenotypic variation of fitness‐related traits both on Corsica and on the mainland. Results The extent of niche space used by birds differed substantially depending on which habitat gradient was considered. Many species have been found to contract their habitat niche along the elevation gradient on Corsica compared with the mainland, whereas all species in the vegetation gradient broadened their niche on the island. Breeding patterns of the blue tit differed considerably depending on whether they settle in deciduous oaks (Quercus humilis) or in evergreen sclerophyllous oaks (Quercus ilex). Phenotypic variation of breeding traits was much higher on the island, where more populations were correctly timed for the best breeding period than on the mainland, a pattern that is likely to result from lower dispersal of organisms on the island. Main conclusions The differences in observed niche breadth between the two series of habitat gradients is explained both by the species‐specific ecology of the species and the human‐induced environmental history of Corsica. Large‐scale landscape changes provided new opportunities for island colonization by non‐forest species, which are isolated as small, ‘fugitive’ local populations. In both gradients, forest species that are typical components of the Corsican bird fauna definitely expanded their niche and occupied a wider range of habitats on Corsica than on the mainland. At the population scale, landscapes included habitat patches with contrasted selection regimes, which resulted in high phenotypic variation for many fitness‐related traits. Reduced dispersal of birds on the island resulted in a much higher degree of local differentiation on Corsica than on the mainland.     

Genetic analysis of differentiation among breeding ponds reveals a candidate gene for local adaptation in Rana arvalis

One of the main questions in evolutionary and conservation biology is how geographical and environmental features of the landscape shape neutral and adaptive genetic variation in natural populations. The identification of genomic polymorphisms that account for adaptive variation can aid in finding candidate loci for local adaptation. Consequently, a comparison of spatial patterns in neutral markers and loci under selection may help disentangle the effects of gene flow, genetic drift and selection at the landscape scale. Many amphibians breed in wetlands, which differ in environmental conditions and in the degree of isolation, enhancing the potential for local adaptation. We used microsatellite markers to measure genetic differentiation among 17 local populations of Rana arvalis breeding in a network of wetlands. We found that locus RC08604 deviated from neutral expectations, suggesting that it is a good candidate for directional selection. We used a genetic network analysis to show that the allele distribution in this locus is correlated with habitat characteristics, whereas this was not the case at neutral markers that displayed a different allele distribution and population network in the study area. The graph approach illustrated the genomic heterogeneity (neutral loci vs. the candidate locus for directional selection) of gene exchange and genetic divergence among populations under directional selection. Limited gene flow between wetlands was only observed at the candidate genomic region under directional selection. RC08604 is partially located inside an up‐regulated thyroid‐hormone receptor (TRβ) gene coordinating the expression of other genes during metamorphosis and appears to be linked with variation in larval life‐history traits found among R. arvalis populations. We suggest that directional selection on genes coding larval life‐history traits is strong enough to maintain the divergence in these genomic regions, reducing the effective recombination of locally adapted alleles but not in other regions of the genome. Integrating this knowledge into conservation plans at the landscape scale will improve the design of management strategies to preserve adaptive genetic diversity in wetland networks.     

Geographical variation in egg mass and egg content in a passerine bird

Reproductive, phenotypic and life-history traits in many animal and plant taxa show geographic variation, indicating spatial variation in selection regimes. Maternal deposition to avian eggs, such as hormones, antibodies and antioxidants, critically affect development of the offspring, with long-lasting effects on the phenotype and fitness. Little is however known about large-scale geographical patterns of variation in maternal deposition to eggs. We studied geographical variation in egg components of a passerine bird, the pied flycatcher (Ficedula hypoleuca), by collecting samples from 16 populations and measuring egg and yolk mass, albumen lysozyme activity, yolk immunoglobulins, yolk androgens and yolk total carotenoids. We found significant variation among populations in most egg components, but ca. 90% of the variation was among individuals within populations. Population however explained 40% of the variation in carotenoid levels. In contrast to our hypothesis, we found geographical trends only in carotenoids, but not in any of the other egg components. Our results thus suggest high within-population variation and leave little scope for local adaptation and genetic differentiation in deposition of different egg components. The role of these maternally-derived resources in evolutionary change should be further investigated.     

Geographic variation of wing morphology of great fruit‐eating bats (Artibeus lituratus): environmental,genetic and spatial correlates of phenotypic differences

Responses of species to environmental gradients are important and frequent determinants of geographic phenotypic variation that can drive adaptive processes. Nonetheless, random genetic processes such as drift can also result in geographic variation in phenotypes, and should be evaluated before implicating selection as the process driving phenotypic change. We examined geographic variation in wing morphology of Artibeus lituratus among 18 different sites distributed across interior Atlantic Forest of Paraguay and Argentina. Moreover, we contrasted geographic variation with environmental, spatial, and genetic variation to test hypotheses related to selection and drift and their impacts on wing morphology. For A. lituratus distributed across interior Atlantic Forest, significant differences among sites characterized variation in wing morphology. Geographic variation was significantly related to climatic variables but not spatial or genetic distances. Such a pattern suggests that phenotypic variation is related to selection for particular environmental regimes, and not genetic drift. Four significant dimensions of phenotypic variation were determined. Three dimensions were related to variation among individuals in terms of wing tips, whereas one was related to overall body size. Wing tips are important for manoeuverability during flight and differences among sites likely reflect differences in forest and vegetation structure that must be managed during foraging. Although climate provides good surrogates for environmental variation, it is probably only an indirect cue of selection regimes that determine variation in wing morphology. Future studies should evaluate more direct environmental measures such as vegetation structure when attempting to interpret geographical variation in wing morphology.     

When can refuges mediate the genetic effects of fire regimes? A simulation study of the effects of topography and weather on neutral and adaptive genetic diversity in fire‐prone landscapes

Understanding how landscape heterogeneity mediates the effects of fire on biodiversity is increasingly important under global changes in fire regimes. We used a simulation experiment to investigate how fire regimes interact with topography and weather to shape neutral and selection‐driven genetic diversity under alternative dispersal scenarios, and to explore the conditions under which microrefuges can maintain genetic diversity of populations exposed to recurrent fire. Spatial heterogeneity in simulated fire frequency occurred in topographically complex landscapes, with fire refuges and fire‐prone “hotspots” apparent. Interannual weather variability reduced the effect of topography on fire patterns, with refuges less apparent under high weather variability. Neutral genetic diversity was correlated with long‐term fire frequency under spatially heterogeneous fire regimes, being higher in fire refuges than fire‐prone areas, except under high dispersal or low fire severity (low mortality). This generated different spatial genetic structures in fire‐prone and fire‐refuge components of the landscape, despite similar dispersal. In contrast, genetic diversity was only associated with time since the most recent fire in flat landscapes without predictable refuges and hotspots. Genetic effects of selection driven by fire‐related conditions depended on selection pressure, migration distance and spatial heterogeneity in fire regimes. Allele frequencies at a locus conferring higher fitness under successional environmental conditions followed a pattern of “temporal adaptation” to contemporary conditions under strong selection pressure and high migration. However, selected allele frequencies were correlated with spatial variation in long‐term mean fire frequency (relating to environmental predictability) under weak dispersal, low selection pressure and strong spatial heterogeneity in fire regimes.     

Molecular divergence in tropical tree populations occupying environmental mosaics

Unveiling the genetic basis of local adaptation to environmental variation is a major goal in molecular ecology. In rugged landscapes characterized by environmental mosaics, living populations and communities can experience steep ecological gradients over very short geographical distances. In lowland tropical forests, interspecific divergence in edaphic specialization (for seasonally flooded bottomlands and seasonally dry terra firme soils) has been proven by ecological studies on adaptive traits. Some species are nevertheless capable of covering the entire span of the gradient; intraspecific variation for adaptation to contrasting conditions may explain the distribution of such ecological generalists. We investigated whether local divergence happens at small spatial scales in two stands of Eperua falcata (Fabaceae), a widespread tree species of the Guiana Shield. We investigated Single Nucleotide Polymorphisms (SNP) and sequence divergence as well as spatial genetic structure (SGS) at four genes putatively involved in stress response and three genes with unknown function. Significant genetic differentiation was observed among sub‐populations within stands, and eight SNP loci showed patterns compatible with disruptive selection. SGS analysis showed genetic turnover along the gradients at three loci, and at least one haplotype was found to be in repulsion with one habitat. Taken together, these results suggest genetic differentiation at small spatial scale in spite of gene flow. We hypothesize that heterogeneous environments may cause molecular divergence, possibly associated to local adaptation in E. falcata.     

Spatial and temporal dynamics in a sexual selection mosaic

Selective regimes and phenotypic optima could either change smoothly and in a clinal fashion or be spatially organized in a more unpredictable mosaic pattern over the geographic landscape. When natural or sexual selection is driven by intra- or interspecific biotic interactions, fine-grained spatial variation in selective regimes could result in selection mosaics rather than clinal variation in selection. We investigated temporal variation and spatial organization in sexual selection on male body size along an ecological coastal-inland gradient of a polymorphic damselfly Ischnura elegans. Body size increased in a clinal fashion along this gradient: animals were smaller in size at the coast, but became larger in the inland areas. In contrast, the sexual selection regimes on male body size showed evidence of more fine-grained spatial organization with no evidence for a clinal pattern and low spatial autocorrelations between populations. These spatially fine-grained sexual selection regimes varied in sign and magnitude and were driven by a combination of the densities of heritable female color morphs and local female body sizes. We suggest that the spatial organization of the selective regimes can be interpreted as a sexual selection mosaic that is influenced by highly localized density- and frequency-dependent social interactions.     

An analysis of selection on a colour polymorphism in the northern leopard frog

In this study, we investigated the role of selection in the maintenance of a dorsal colour polymorphism in natural populations of the northern leopard frog, Rana pipiens. We determined genetic structure both spatially and temporally from a suite of putatively neutral molecular markers and tested whether or not the colour locus exhibited patterns of genetic variation that differed from those of the neutral loci. Spatial genetic structure at the colour locus was indistinguishable from structure at neutral loci [95% confidence intervals of F(ST) (neutral) = (0.07, 0.35), F(ST) (colour locus) = 0.114]. In the temporal analysis, we found that the variance among populations in the change in allele frequency at the colour locus (equal to 0.004) lies within the 95% confidence intervals for the variance among populations in changes in allele frequencies at neutral loci. In light of our inability to show evidence for the selective maintenance of the colour polymorphism, we used computer simulations to infer the power of our spatial and temporal techniques to detect selection. The computer simulations showed that although the strength of selection (s) would need to be relatively strong to have been detected by the temporal approach (s = 0.1-0.4, depending on the model tested), the spatial analysis would have detected all but weak selection (s = 0.01-0.04, depending on the model tested). This study illustrates the importance of using a locus comparison approach to detect evidence for selective maintenance before conducting studies to measure the selective mechanisms maintaining a polymorphism.