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.     

Constrained evolution of the sex comb in Drosophila simulans

Male fitness is dependent on sexual traits that influence mate acquisition (precopulatory sexual selection) and paternity (post‐copulatory sexual selection), and although many studies have documented the form of selection in one or the other of these arenas, fewer have done it for both. Nonetheless, it appears that the dominant form of sexual selection is directional, although theoretically, populations should converge on peaks in the fitness surface, where selection is stabilizing. Many factors, however, can prevent populations from reaching adaptive peaks. Genetic constraints can be important if they prevent the development of highest fitness phenotypes, as can the direction of selection if it reverses across episodes of selection. In this study, we examine the evidence that these processes influence the evolution of the multivariate sex comb morphology of male Drosophila simulans. To do this, we conduct a quantitative genetic study together with a multivariate selection analysis to infer how the genetic architecture and selection interact. We find abundant genetic variance and covariance in elements of the sex comb. However, there was little evidence for directional selection in either arena. Significant nonlinear selection was detected prior to copulation when males were mated to nonvirgin females, and post‐copulation during sperm offence (again with males mated to nonvirgins). Thus, contrary to our predictions, the evolution of the D. simulans sex comb is limited neither by genetic constraints nor by antagonistic selection between pre‐ and post‐copulatory arenas, but nonlinear selection on the multivariate phenotype may prevent sex combs from evolving to reach some fitness maximizing optima.     

Spatiotemporal environmental heterogeneity and the maintenance of the tailspot polymorphism in the variable platyfish (Xiphophorus variatus)

Genetic variation is critical for adaptive evolution. Despite its importance, there is still limited evidence in support of some prominent theoretical models explaining the maintenance of genetic polymorphism within populations. We examined 84 populations of Xiphophorus variatus, a livebearing fish with a genetic polymorphism associated with physiological performance, to test: (1) whether niche differentiation explains broad‐scale maintenance of polymorphism, (2) whether polymorphism is maintained among populations by local adaptation and migration, or (3) whether heterogeneity in explicit environmental variables could be linked to levels of polymorphism within populations. We found no evidence of climatic niche differentiation that could generate or maintain broad geographic variation in polymorphism. Subsequently, hierarchical partitioning of genetic richness and partial mantel tests revealed that 76% of the observed genetic richness was partitioned within populations with no effect of geographic distance on polymorphism. These results strongly suggest a lack of migration‐selection balance in the maintenance of polymorphism, and model selection confirmed a significant relationship between environmental heterogeneity and genetic richness within populations. Few studies have demonstrated such effects at this scale, and additional studies in other taxa should examine the generality of gene‐by‐environment interactions across populations to better understand the dynamics and scale of balancing selection.     

Experimental evidence for multivariate stabilizing sexual selection

Stabilizing selection is a fundamental concept in evolutionary biology. In the presence of a single intermediate optimum phenotype (fitness peak) on the fitness surface, stabilizing selection should cause the population to evolve toward such a peak. This prediction has seldom been tested, particularly for suites of correlated traits. The lack of tests for an evolutionary match between population means and adaptive peaks may be due, at least in part, to problems associated with empirically detecting multivariate stabilizing selection and with testing whether population means are at the peak of multivariate fitness surfaces. Here we show how canonical analysis of the fitness surface, combined with the estimation of confidence regions for stationary points on quadratic response surfaces, may be used to define multivariate stabilizing selection on a suite of traits and to establish whether natural populations reside on the multivariate peak. We manufactured artificial advertisement calls of the male cricket Teleogryllus commodus and played them back to females in laboratory phonotaxis trials to estimate the linear and nonlinear sexual selection that female phonotactic choice imposes on male call structure. Significant nonlinear selection on the major axes of the fitness surface was convex in nature and displayed an intermediate optimum, indicating multivariate stabilizing selection. The mean phenotypes of four independent samples of males, from the same population as the females used in phonotaxis trials, were within the 95% confidence region for the fitness peak. These experiments indicate that stabilizing sexual selection may play an important role in the evolution of male call properties in natural populations of T. commodus.     

Genomics,environment and balancing selection in behaviourally bimodal populations: The caribou case

Selection forces that favour different phenotypes in different environments can change frequencies of genes between populations along environmental clines. Clines are also compatible with balancing forces, such as negative frequency‐dependent selection (NFDS), which maintains phenotypic polymorphisms within populations. For example, NFDS is hypothesized to maintain partial migration, a dimorphic behavioural trait prominent in species where only a fraction of the population seasonally migrates. Overall, NFDS is believed to be a common phenomenon in nature, yet a scarcity of studies were published linking naturally occurring allelic variation with bimodal or multimodal phenotypes and balancing selection. We applied a Pool‐seq approach and detected selection on alleles associated with environmental variables along a North–South gradient in western North American caribou, a species displaying partially migratory behaviour. On 51 loci, we found a signature of balancing selection, which could be related to NFDS and ultimately the maintenance of the phenotypic polymorphisms known within these populations. Yet, remarkably, we detected directional selection on a locus when our sample was divided into two behaviourally distinctive groups regardless of geographic provenance (a subset of GPS‐collared migratory or sedentary individuals), indicating that, within populations, phenotypically homogeneous groups were genetically distinctive. Loci under selection were linked to functional genes involved in oxidative stress response, body development and taste perception. Overall, results indicated genetic differentiation along an environmental gradient of caribou populations, which we found characterized by genes potentially undergoing balancing selection. We suggest that the underlining balancing force, NFDS, plays a strong role within populations harbouring multiple haplotypes and phenotypes, as it is the norm in animals, plants and humans too.     

ESTIMATING AND VISUALIZING FITNESS SURFACES USING MARK–RECAPTURE DATA

Understanding how selection operates on a set of phenotypic traits is central to evolutionary biology. Often, it requires estimating survival (or other fitness‐related life‐history traits) which can be difficult to obtain for natural populations because individuals cannot be exhaustively followed. To cope with this issue of imperfect detection, we advocate the use of mark‐recapture data and we provide a general framework for both the estimation of linear and nonlinear selection gradients and the visualization of fitness surfaces. To quantify the strength of selection, the standard second‐order polynomial regression method is integrated in mark‐recapture models. To visualize the form of selection, we use splines to display selection acting on multivariate phenotypes in the most flexible way. We employ Markov chain Monte Carlo sampling in a Bayesian framework to estimate model parameters, assessing traits relevance and calculating the optimal amount of smoothing. We illustrate our approach using data from a wild population of Common blackbirds (Turdus merula) to investigate survival in relation to morphological traits, and provide evidence for correlational selection using the new methodology. Overall, the framework we propose will help in exploring the full potential of mark‐recapture data to study natural selection.     

Evidence of adaptive divergence in plasticity: density- and site-dependent selection on shade-avoidance responses in Impatiens capensis

We investigated the conditions under which plastic responses to density are adaptive in natural populations of Impatiens capensis and determined whether plasticity has evolved differently in different selective environments. Previous studies showed that a population that evolved in a sunny site exhibited greater plasticity in response to density than did a population that evolved in a woodland site. Using replicate inbred lines in a reciprocal transplant that included a density manipulation, we asked whether such population differentiation was consistent with the hypothesis of adaptive divergence. We hypothesized that plasticity would be more strongly favored in the sunny site than in the woodland site; consequently, we predicted that selection would be more strongly density dependent in the sunny site, favoring the phenotype that was expressed at each density. Selection on internode length and flowering date was consistent with the hypothesis of adaptive divergence in plasticity. Few costs or benefits of plasticity were detected independently from the expressed phenotype, so plasticity was selected primarily through selection on the phenotype. Correlations between phenotypes and their plasticity varied with the environment and would cause indirect selection on plasticity to be environment dependent. We showed that an appropriate plastic response even to a rare environment can greatly increase genotypic fitness when that environment is favorable. Selection on the measured characters contributed to local adaptation and fully accounted for fitness differences between populations in all treatments except the woodland site at natural density.     

Quantification and decomposition of environment‐selection relationships

In nature, selection varies across time in most environments, but we lack an understanding of how specific ecological changes drive this variation. Ecological factors can alter phenotypic selection coefficients through changes in trait distributions or individual mean fitness, even when the trait‐absolute fitness relationship remains constant. We apply and extend a regression‐based approach in a population of Soay sheep (Ovis aries) and suggest metrics of environment‐selection relationships that can be compared across studies. We then introduce a novel method that constructs an environmentally structured fitness function. This allows calculation of full (as in existing approaches) and partial (acting separately through the absolute fitness function slope, mean fitness, and phenotype distribution) sensitivities of selection to an ecological variable. Both approaches show positive overall effects of density on viability selection of lamb mass. However, the second approach demonstrates that this relationship is largely driven by effects of density on mean fitness, rather than on the trait‐fitness relationship slope. If such mechanisms of environmental dependence of selection are common, this could have important implications regarding the frequency of fluctuating selection, and how previous selection inferences relate to longer term evolutionary dynamics.     

Social traits,social networks and evolutionary biology

The social environment is both an important agent of selection for most organisms, and an emergent property of their interactions. As an aggregation of interactions among members of a population, the social environment is a product of many sets of relationships and so can be represented as a network or matrix. Social network analysis in animals has focused on why these networks possess the structure they do, and whether individuals’ network traits, representing some aspect of their social phenotype, relate to their fitness. Meanwhile, quantitative geneticists have demonstrated that traits expressed in a social context can depend on the phenotypes and genotypes of interacting partners, leading to influences of the social environment on the traits and fitness of individuals and the evolutionary trajectories of populations. Therefore, both fields are investigating similar topics, yet have arrived at these points relatively independently. We review how these approaches are diverged, and yet how they retain clear parallelism and so strong potential for complementarity. This demonstrates that, despite separate bodies of theory, advances in one might inform the other. Techniques in network analysis for quantifying social phenotypes, and for identifying community structure, should be useful for those studying the relationship between individual behaviour and group‐level phenotypes. Entering social association matrices into quantitative genetic models may also reduce bias in heritability estimates, and allow the estimation of the influence of social connectedness on trait expression. Current methods for measuring natural selection in a social context explicitly account for the fact that a trait is not necessarily the property of a single individual, something the network approaches have not yet considered when relating network metrics to individual fitness. Harnessing evolutionary models that consider traits affected by genes in other individuals (i.e. indirect genetic effects) provides the potential to understand how entire networks of social interactions in populations influence phenotypes and predict how these traits may evolve. By theoretical integration of social network analysis and quantitative genetics, we hope to identify areas of compatibility and incompatibility and to direct research efforts towards the most promising areas. Continuing this synthesis could provide important insights into the evolution of traits expressed in a social context and the evolutionary consequences of complex and nuanced social phenotypes.     

Social phenotype extended to communities: Expanded multilevel social selection analysis reveals fitness consequences of interspecific interactions

In social species, fitness consequences are associated with both individual and social phenotypes. Social selection analysis has quantified the contribution of conspecific social traits to individual fitness. There has been no attempt, however, to apply a social selection approach to quantify the fitness implications of heterospecific social phenotypes. Here, we propose a novel social selection based approach integrating the role of all social interactions at the community level. We extended multilevel selection analysis by including a term accounting for the group phenotype of heterospecifics. We analyzed nest activity as a model social trait common to two species, the lesser kestrel (Falco naumanni) and jackdaw (Corvus monedula), nesting in either single‐ or mixed‐species colonies. By recording reproductive outcome as a measure of relative fitness, our results reveal an asymmetric system wherein only jackdaw breeding performance was affected by the activity phenotypes of both conspecific and heterospecific neighbors. Our model incorporating heterospecific social phenotypes is applicable to animal communities where interacting species share a common social trait, thus allowing an assessment of the selection pressure imposed by interspecific interactions in nature. Finally, we discuss the potential role of ecological limitations accounting for random or preferential assortments among interspecific social phenotypes, and the implications of such processes to community evolution.     

SOLVING THE PARADOX OF STASIS: SQUASHED STABILIZING SELECTION AND THE LIMITS OF DETECTION

Despite the potential for rapid evolution, stasis is commonly observed over geological timescales—the so‐called “paradox of stasis.” This paradox would be resolved if stabilizing selection were common, but stabilizing selection is infrequently detected in natural populations. We hypothesize a simple solution to this apparent disconnect: stabilizing selection is hard to detect empirically once populations have adapted to a fitness peak. To test this hypothesis, we developed an individual‐based model of a population evolving under an invariant stabilizing fitness function. Stabilizing selection on the population was infrequently detected in an “empirical” sampling protocol, because (1) trait variation was low relative to the fitness peak breadth; (2) nonselective deaths masked selection; (3) populations wandered around the fitness peak; and (4) sample sizes were typically too small. Moreover, the addition of negative frequency‐dependent selection further hindered detection by flattening or even dimpling the fitness peak, a phenomenon we term “squashed stabilizing selection.” Our model demonstrates that stabilizing selection provides a plausible resolution to the paradox of stasis despite its infrequent detection in nature. The key reason is that selection “erases its traces”: once populations have adapted to a fitness peak, they are no longer expected to exhibit detectable stabilizing selection.     

Evolutionary genetics of maternal effects

Maternal genetic effects (MGEs), where genes expressed by mothers affect the phenotype of their offspring, are important sources of phenotypic diversity in a myriad of organisms. We use a single‐locus model to examine how MGEs contribute patterns of heritable and nonheritable variation and influence evolutionary dynamics in randomly mating and inbreeding populations. We elucidate the influence of MGEs by examining the offspring genotype‐phenotype relationship, which determines how MGEs affect evolutionary dynamics in response to selection on offspring phenotypes. This approach reveals important results that are not apparent from classic quantitative genetic treatments of MGEs. We show that additive and dominance MGEs make different contributions to evolutionary dynamics and patterns of variation, which are differentially affected by inbreeding. Dominance MGEs make the offspring genotype‐phenotype relationship frequency dependent, resulting in the appearance of negative frequency‐dependent selection, while additive MGEs contribute a component of parent‐of‐origin dependent variation. Inbreeding amplifies the contribution of MGEs to the additive genetic variance and, therefore enhances their evolutionary response. Considering evolutionary dynamics of allele frequency change on an adaptive landscape, we show that this landscape differs from the mean fitness surface, and therefore, under some condition, fitness peaks can exist but not be “available” to the evolving population.     

Heterosis in age‐specific selected populations of a seed beetle: Sex differences in longevity and reproductive behavior

We tested mutation accumulation hypothesis for the evolution of senescence using short‐lived and long‐lived populations of the seed‐feeding beetle, Acanthoscelides obtectus (Say), obtained by selection on early‐ and late‐life for many generations. The expected consequence of the mutation accumulation hypothesis is that in short‐lived populations, where the force of natural selection is the strongest early in life, the late‐life fitness traits should decline due to genetic drift which increases the frequency of mutations with deleterious effects in later adult stages. Since it is unlikely that identical deleterious mutations will increase in several independent populations, hybrid vigor for late‐life fitness is expected in offspring obtained in crosses among populations selected for early‐life fitness traits. We tested longevity of both sexes, female fecundity and male reproductive behavior for hybrid vigor by comparing hybrid and nonhybrid short‐lived populations. Hybrid vigor was confirmed for male virility, mating speed and copulation duration, and longevity of both sexes at late ages. In contrast to males, the results on female fecundity in short‐lived populations did not support mutation accumulation as a genetic mechanism for the evolution of this trait. Contrary to the prediction of this hypothesis, male mating ability indices and female fecundity in long‐lived populations exhibited hybrid vigor at all assayed age classes. We demonstrate that nonhybrid long‐lived populations diverged randomly regarding female and male reproductive fitness, indicating that sexually antagonistic selection, when accompanied with genetic drift for female fecundity and male virility, might be responsible for overriding natural selection in the independently evolving long‐lived populations.     

WHY WE ARE NOT DEAD ONE HUNDRED TIMES OVER

The possibility of pervasive weak selection at tens or hundreds of millions of sites across the genome, suggested by recent studies of silent site DNA sequence variation and divergence, raises the problem of the survival of the population in the face of the large genetic load that may result. Two alternative resolutions of this problem are presented for populations where recombination is sufficiently frequent that different sites under selection evolve independently. One invokes weak stabilizing selection, of the magnitude compatible with abundant silent site variability. This can be shown to produce only a modest genetic load, due to the effectiveness of even weak stabilizing selection in keeping the trait mean close to the optimum. The other invokes soft selection, whereby individuals compete for a limiting resource whose abundance determines the absolute fitness of the population. Weak purifying selection at a large number of sites produces only a small variance in fitness among individuals within the population, due to the fact that most sites are fixed rather than polymorphic. Even when it produces a large genetic load, it is compatible with the observations on fitness variance when selection is soft. It may be very difficult to distinguish between these two possibilities.     

SOCIAL CAUSES OF CORRELATIONAL SELECTION AND THE RESOLUTION OF A HERITABLE THROAT COLOR POLYMORPHISM IN A LIZARD

Abstract When selection acts on social or behavioral traits, the fitness of an individual depends on the phenotypes of its competitors. Here, we describe methods and statistical inference for measuring natural selection in small social groups. We measured selection on throat color alleles that arises from microgeographic variation in allele frequency at natal sites of side‐blotched lizards (Uta stansburiana). Previous game‐theoretic analysis indicates that two color morphs of female side‐blotched lizards are engaged in an offspring quantity‐quality game that promotes a density‐and frequency‐dependent cycle. Orange‐throated females are r‐strategists. They lay large clutches of small progeny, which have poor survival at high density, but good survival at low density. In contrast, yellow‐throated females are K‐strategists. They lay small clutches of large progeny, which have good survival at high density. We tested three predictions of the female game: (1) orange progeny should have a fitness advantage at low density; (2) correlational selection acts to couple color alleles and progeny size; and (3) this correlational selection arises from frequency‐dependent selection in which large hatchling size confers an advantage, but only when yellow alleles are rare. We also confirmed the heritability of color, and therefore its genetic basis, by producing progeny from controlled matings. A parsimonious cause of the high heritability is that three alleles (o, b, y) segregate as one genetic factor. We review the physiology of color formation to explain the possible genetic architecture of the throat color trait. Heritability of color was nearly additive in our breeding study, allowing us to compute a genotypic value for each individual and thus predict the frequency of progeny alleles released on 116 plots. Rather than study the fitness of individual progeny, we studied how the fitness of their color alleles varied with allele frequency on plots. We confirmed prediction 1: When orange alleles are present in female progeny, they have higher fitness at low density when compared to other alleles. Even though the difference in egg size of the female morphs was small (0.02 g), it led to knife‐edged survival effects for their progeny depending on local social context. Selection on hatchling survival was not only dependent on color alleles, but on a fitness interaction between color alleles and hatchling size, which confirmed prediction 2. Sire effects, which are not confounded by maternal phenotype, allowed us to resolve the frequency dependence of correlational selection on egg size and color alleles and thereby confirmed prediction 3. Selection favored large size when yellow sire alleles were rare, but small size when they were common. Correlational selection promotes the formation of a self‐reinforcing genetic correlation between the morphs and life‐history variation, which causes selection in the next density and frequency cycle to be exacerbated. We discuss general conditions for the evolution of self‐reinforcing genetic correlations that arise from social selection associated with frequency‐dependent sexual and natural selection.     

Integrating large‐scale geographic patterns in flight morphology,flight characteristics and sexual selection in a range‐expanding damselfly

While geographic trait variation along environmental clines is widespread, associated patterns in sexual selection remain largely unexplored. Geographic patterns in sexual selection may be expected if 1) phenotypes vary geographically and sexual selection is dependent on the local phenotypes in the population, and if 2) sexual selection is influenced by geographically structured environmental conditions. We quantified geographic variation in flight‐related traits and flight performance in mated and unmated males and tested for geographic variation in sexual selection on these traits in the poleward range‐expanding damselfly Coenagrion scitulum across a set of eleven core and edge populations ordered along thermal gradients in the larval and in the adult stage. We found little support for trait differentiation between core and edge populations, instead we found considerable geographic trait variation along the larval and adult thermal gradients. As expected under time constraints, body mass decreased with shorter larval growth seasons. Lower temperatures during the adult flight period were associated with a higher body mass, a higher flight speed and a higher fat content; these traits likely evolved to buffer flight ability at suboptimal temperatures and to optimize starvation resistance. Across the large geographic scale, we found a consistent higher flight duration in mated males. Instead, sexual selection for higher fat content was stronger in populations with lower adult flight temperatures and sexual selection for lower body mass acted only in edge populations. Our results indicate sexual selection on flight performance to be consistent over a large geographic scale and this despite the clear geographic patterns in sexual selection on the underlying morphological traits. Our results highlight that to fully understand the fitness implications of geographically changing trait patterns, researchers should consider the entire phenotype–performance–fitness axis and incorporate effects of geographically structured life‐stage specific environmental conditions on this axis.     

When maladaptive gene flow does not increase selection

Populations receiving high maladaptive gene flow are expected to experience strong directional selection—because gene flow pulls mean phenotypes away from local fitness peaks. We tested this prediction by means of a large and replicated mark‐recapture study of threespine stickleback (Gasterosteus aculeatus) in two stream populations. One of the populations (outlet) experiences high gene flow from the lake population and its morphology is correspondingly poorly adapted. The other population (inlet) experiences very low gene flow from the lake population and its morphology is correspondingly well adapted. Contrary to the above prediction, selection was not stronger in the outlet than in the inlet, a result that forced us to consider potential reasons for why maladaptive gene flow might not increase selection. Of particular interest, we show by means of a simple population genetic model that maladaptive gene flow can—under reasonable conditions—decrease the strength of directional selection. This outcome occurs when immigrants decrease mean fitness in the resident population, which decreases the strength of selection against maladapted phenotypes. We argue that this previously unrecognized effect of gene flow deserves further attention in theoretical and empirical studies.     

Contribution of maternal effects to dietary selection in Mediterranean fruit flies

Individual responses to dietary variation represent a fundamental component of fitness, and nutritional adaptation can occur over just a few generations. Maternal effects can show marked proximate responses to nutrition, but whether they contribute to longer term dietary adaptation is unclear. Here, we tested the hypotheses that maternal effects: (i) contribute to dietary adaptation, (ii) diminish when dietary conditions are constant between generations, (iii) are trait‐specific and (iv) interact with high‐ and low‐quality food. We used experimental evolution regimes in the medfly (Ceratitis capitata) to test these predictions by subjecting an outbred laboratory‐adapted population to replicated experimental evolution on either constant high calorie sugar (‘A’) or low‐calorie starch (‘S’) larval diets, with a standard adult diet across both regimes. We measured the contribution of maternal effects by comparing developmental and adult phenotypes of individuals reared on their own diet with those swapped onto the opposite diet for either one or two generations (high and low maternal effect conditions, respectively), both at the start and after 30 generations of selection. Initially, there were strong maternal effects on female body mass and male mating success but not larval survival. Interestingly, the initial maternal effects observed in female body mass and male mating success showed sex‐specific interactions when individuals from high calorie regimes were tested on low calorie diets. However, as populations responded to selection, the effects of maternal provisioning on all traits diminished. The results broadly supported the predictions. They show how the contribution of maternal effects to dietary responses evolves in a context‐dependent manner, with significant variation across different fitness‐related traits. We conclude that maternal effects can evolve during nutritional adaptation and hence may be an important life history trait to measure, rather than to routinely minimize.     

Scale dependence of sex ratio in wild plant populations: implications for social selection

Social context refers to the composition of an individual''s social interactants, including potential mates. In spatially structured populations, social context can vary among individuals within populations, generating the opportunity for social selection to drive differences in fitness functions among individuals at a fine spatial scale. In sexually polymorphic plants, the local sex ratio varies at a fine scale and thus has the potential to generate this opportunity. We measured the spatial distribution of two wild populations of the gynodioecious plant Silene vulgaris and show that there is fine‐scale heterogeneity in the local distribution of the sexes within these populations. We demonstrate that the largest variance in sex ratio is among nearest neighbors. This variance is greatly reduced as the spatial scale of social interactions increases. These patterns suggest the sex of neighbors has the potential to generate fine‐scale differences in selection differentials among individuals. One of the most important determinants of social interactions in plants is the behavior of pollinators. These results suggest that the potential for selection arising from sex ratio will be greatest when pollen is shared among nearest neighbors. Future studies incorporating the movement of pollinators may reveal whether and how this fine‐scale variance in sex ratio affects the fitness of individuals in these populations.