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.     

Characterizing selection in black‐throated blue warblers using a sexual network approach

Our understanding of trait evolution is built upon studies that examine the correlation between traits and fitness, most of which implicitly assume all individuals experience similar selective environments. However, accounting for differences in selective pressures, such as variation in the social environment, can advance our understanding of how selection shapes individual traits and subsequent fitness. In this study, we test whether variation in the social environment affects selection on individual phenotype. We apply a new sexual network framework to quantify each male's social environment as the mean body size of his primary competitors. We test for direct and social selection on male body size using a 10‐year data set on black‐throated blue warblers (Setophaga caerulescens), a territorial species for which body size is hypothesized to mediate competition for mates. We found that direct selection on body size was weak and nonsignificant, as was social selection via the body size of the males' competitors. Analysing both types of selection simultaneously allows us to firmly reject a role for body size in competitive interactions between males and subsequent male fitness in this population. We evaluate the application of the sexual network approach to empirical data and suggest that other phenotypic traits such as song characteristics and plumage may be more relevant than body size for male–male competition in this small passerine bird.     

On the standardization of fitness and traits in comparative studies of phenotypic selection

Comparisons of the strength and form of phenotypic selection among groups provide a powerful approach for testing adaptive hypotheses. A central and largely unaddressed issue is how fitness and phenotypes are standardized in such studies; standardization across or within groups can qualitatively change conclusions whenever mean fitness differs between groups. We briefly reviewed recent relevant literature, and found that selection studies vary widely in their scale of standardization, but few investigators motivated their rationale for chosen standardization approaches. Here, we propose that the scale at which fitness should be relativized should reflect whether selection is likely to be hard or soft; that is, the scale at which populations (or hypothetical populations in the case of a contrived experiment) are regulated. We argue that many comparative studies of selection are implicitly or explicitly focused on soft selection (i.e., frequency and density‐dependent selection). In such studies, relative fitness should preferably be calculated using within‐group means, although this approach is taken only occasionally. Related difficulties arise for the standardization of phenotypes. The appropriate scale at which standardization should take place depends on whether groups are considered to be fixed or random. We emphasize that the scale of standardization is a critical decision in empirical studies of selection that should always warrant explicit justification.     

An integrated analysis of phenotypic selection on insect body size and development time

Most studies of phenotypic selection do not estimate selection or fitness surfaces for multiple components of fitness within a unified statistical framework. This makes it difficult or impossible to assess how selection operates on traits through variation in multiple components of fitness. We describe a new generation of aster models that can evaluate phenotypic selection by accounting for timing of life‐history transitions and their effect on population growth rate, in addition to survival and reproductive output. We use this approach to estimate selection on body size and development time for a field population of the herbivorous insect, Manduca sexta (Lepidoptera: Sphingidae). Estimated fitness surfaces revealed strong and significant directional selection favoring both larger adult size (via effects on egg counts) and more rapid rates of early larval development (via effects on larval survival). Incorporating the timing of reproduction and its influence on population growth rate into the analysis resulted in larger values for size in early larval development at which fitness is maximized, and weaker selection on size in early larval development. These results illustrate how the interplay of different components of fitness can influence selection on size and development time. This integrated modeling framework can be readily applied to studies of phenotypic selection via multiple fitness components in other systems.     

The causes of selection on flowering time through male fitness in a hermaphroditic annual plant

Flowering is a key life‐history event whose timing almost certainly affects both male and female fitness, but tests of selection on flowering time through male fitness are few. Such selection may arise from direct effects of flowering time, and indirect effects through covariance between flowering time and the environment experienced during reproduction. To isolate these intrinsically correlated associations, we staggered planting dates of Brassica rapa families with known flowering times, creating populations in which age at flowering (i.e., flowering time genotype) and Julian date of flowering (i.e., flowering time environment) were positively, negatively, or uncorrelated. Genetic paternity analysis revealed that male fitness was not strongly influenced by seasonal environmental changes. Instead, when age and date were uncorrelated, selection through male fitness strongly favored young age at flowering. Strategic sampling offspring for paternity analysis rejected covariance between sire age at flowering and dam quality as the cause of this selection. Results instead suggest a negative association between age at flowering and pollen competitive ability. The manipulation also revealed that, at least in B. rapa, the often‐observed correlation between flowering time and flowering duration is environmental, not genetic, in origin.     

Fine‐scale temporal adaptation within a salmonid population: mechanism and consequences

We demonstrate a clear example of local adaptation of seasonal timing of spawning and embryo development. The consequence is a population of pink salmon that is segmented into spawning groups that use the same limited habitat. We synthesize published observations with results of new analyses to demonstrate that genetic variation of these traits results in survival differentials related to that variation, and that density‐dependent embryo mortality and seasonally variable juvenile mortality are a mechanism of selection. Most examples of local adaptation in natural systems depend on observed correlations between environments and fitness traits, but do not fully demonstrate local adaptation: that the trait is genetically determined, exhibits different fitness in common environments or across different environments, and its variation is mechanistically connected to fitness differences. The geographic or temporal scales of local adaptation often remain obscure. Here, we show that heritable, fine‐scale differences of timing of reproductive migration in a pink salmon (Oncorhynchus gorbuscha) resulted in temporal structure that persisted several generations; the differences enable a density‐dependent population to pack more spawners into limited spawning habitat, that is, enhance its fitness. A balanced trade‐off of survivals results because embryos from early‐migrating fish have a lower freshwater survival (harsh early physical conditions and disturbance by late spawners), but emigrant fry from late‐migrating fish have lower marine survivals (timing of their vernal emergence into the estuarine environment). Such fine‐scale local adaptations increase the genetic portfolio of the populations and may provide a buffer against the impacts of climate change.     

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.     

Pervasive,yet idiosyncratic,epistatic pleiotropy during adaptation in a behaviourally complex microbe

Understanding how multiple mutations interact to jointly impact multiple ecologically important traits is critical for creating a robust picture of organismal fitness and the process of adaptation. However, this is complicated by both environmental heterogeneity and the complexity of genotype‐to‐phenotype relationships generated by pleiotropy and epistasis. Moreover, little is known about how pleiotropic and epistatic relationships themselves change over evolutionary time. The soil bacterium Myxococcus xanthus employs several distinct social traits across a range of environments. Here, we use an experimental lineage of M. xanthus that evolved a novel form of social motility to address how interactions between epistasis and pleiotropy evolve. Specifically, we test how mutations accumulated during selection on soft agar pleiotropically affect several other social traits (hard agar motility, predation and spore production). Relationships between changes in swarming rate in the selective environment and the four other traits varied greatly over time in both direction and magnitude, both across timescales of the entire evolutionary lineage and individual evolutionary time steps. We also tested how a previously defined epistatic interaction is pleiotropically expressed across these traits. We found that phenotypic effects of this epistatic interaction were highly correlated between soft and hard agar motility, but were uncorrelated between soft agar motility and predation, and inversely correlated between soft agar motility and spore production. Our results show that ‘epistatic pleiotropy’ varied greatly in magnitude, and often even in sign, across traits and over time, highlighting the necessity of simultaneously considering the interacting complexities of pleiotropy and epistasis when studying the process of adaptation.     

Paths to selection on life history loci in different natural environments across the native range of Arabidopsis thaliana

Selection on quantitative trait loci (QTL) may vary among natural environments due to differences in the genetic architecture of traits, environment‐specific allelic effects or changes in the direction and magnitude of selection on specific traits. To dissect the environmental differences in selection on life history QTL across climatic regions, we grew a panel of interconnected recombinant inbred lines (RILs) of Arabidopsis thaliana in four field sites across its native European range. For each environment, we mapped QTL for growth, reproductive timing and development. Several QTL were pleiotropic across environments, three colocalizing with known functional polymorphisms in flowering time genes (CRY2, FRI and MAF2‐5), but major QTL differed across field sites, showing conditional neutrality. We used structural equation models to trace selection paths from QTL to lifetime fitness in each environment. Only three QTL directly affected fruit number, measuring fitness. Most QTL had an indirect effect on fitness through their effect on bolting time or leaf length. Influence of life history traits on fitness differed dramatically across sites, resulting in different patterns of selection on reproductive timing and underlying QTL. In two oceanic field sites with high prereproductive mortality, QTL alleles contributing to early reproduction resulted in greater fruit production, conferring selective advantage, whereas alleles contributing to later reproduction resulted in larger size and higher fitness in a continental site. This demonstrates how environmental variation leads to change in both QTL effect sizes and direction of selection on traits, justifying the persistence of allelic polymorphism at life history QTL across the species range.     

Experimental demonstration of the impact of hard and soft selection regimes on polymorphism maintenance in spatially heterogeneous environments†

Predicting and managing contemporary adaption requires a proper understanding of the determinants of genetic variation. Spatial heterogeneity of the environment may stably maintain polymorphism when habitat contribution to the next generation can be considered independent of the degree of adaptation of local populations within habitats (i.e., under soft selection). In contrast, when habitats contribute proportionally to the mean fitness of the populations they host (hard selection), polymorphism is not expected to be maintained by selection. Although mathematically established decades ago, this prediction had never been demonstrated experimentally. Here, we provide an experimental test in which polymorphic populations of Escherichia coli growing in heterogeneous habitats were exposed to hard and soft selection regimes. As predicted by theory, polymorphism was preserved longer under soft selection. Complementary tests established that soft selection slowed fixation processes and could even protect polymorphism in the long term by providing a systematic advantage to rare genotypes.     

The effects of competition on the strength and softness of selection

Deleterious alleles are constantly introduced into populations due to mutation. In subdivided populations, the impact of these mutations depends on the strength of selection as well as the softness of selection, that is, the extent to which fitness is governed by local rather than global competition. It is widely appreciated that the intensity and type of competition will affect selection on deleterious mutations but most empirical work has focused solely on the effects of competition on selection strength. However, competition has rarely been studied in the context of selection ‘softness’ even though competition is at the conceptual root of soft selection. All other things being equal, theory predicts that inter‐ and intraspecific competitions have opposing effects on the softness of selection. Using Drosophila melanogaster, we estimated the strength and softness of selection in a ‘baseline’ competitive environment as well as two additional competitive environments characterized by either additional intra‐ or interspecific competitors. We found that competitive environment had little effect on the average strength of selection. While the softness of selection was affected by the type of competition, the direction of change varied across tests of different genes, contrary to expectation. Although the ‘hard/soft’ selection paradigm implicitly assumes that all individuals are equally sensitive to the local competitive environment, we found this not to be the case. Wild‐type individuals were more sensitive to changes in the genetic quality of their local competitors than mutant individuals.     

LIFE-CYCLE COMPONENTS OF SELECTION IN ERIGERON ANNUUS: I. PHENOTYPIC SELECTION

The magnitude and direction of phenotypic selection on emergence date and seedling size in Erigeron annuus was measured to determine the heterogeneity of selection among sites and the proportion of fitness variance explained by seedling size and emergence date. Three disturbance treatments (open, annual vegetation, perennial vegetation) were imposed to test the hypothesis of stronger selection on seedlings in competitive environments. Selection was most heterogeneous early in the life cycle, with significant spatial heterogeneity in the magnitude of selection on a local scale. The disturbance treatments affected only fecundity selection gradients and selection was strongest in open plots. Significant variation in the sign of selection differentials on emergence date was observed for establishment and fall viability selection episodes; at later stages selection varied in magnitude but not direction. Seedlings in the earliest cohort experienced high mortality during establishment, but increased size and fecundity later in the life cycle. Both stabilizing and disruptive selection on emergence date were observed during establishment, but in general selection was purely directional. At Stony Brook most selection on emergence date operated indirectly through seedling size, whereas at the Weld Preserve direct selection was stronger. There were persistent effects of both seedling emergence date and rosette diameter on adult fitness components, and October rosette diameter explained 18% of the total phenotypic variance in fecundity. Overall, viability fitness components were much more important than fecundity selection. Winter survivorship was the single most important episode of selection.     

Phenotypic selection and covariation in the life‐history traits of elephant seals: heavier offspring gain a double selective advantage

Early developmental conditions contribute to individual heterogeneity of both phenotypic traits and fitness components, ultimately affecting population dynamics. Although the demographic consequences of ontogenic growth are best quantified using an integrated measure of fitness, most analyses to date have instead studied individual fitness components in isolation. Here, we estimated phenotypic selection on weaning mass in female southern elephant seals Mirounga leonina by analyzing individual‐based data collected between 1986 and 2016 with capture–recapture and matrix projection models. In support of a hypothesis predicting a gradual decrease of weaning mass effects with time since weaning (the replacement hypothesis), we found that the estimated effects of weaning mass on future survival and recruitment probability was of intermediate duration (rather than transient or permanent). Heavier female offspring had improved odds of survival in early life and a higher probability to recruit at an early age. The positive link between weaning mass and recruitment age is noteworthy, considering that pre‐recruitment mortality already imposed a strong selective filter on the population, leaving only the most ‘robust’ individuals to reproduce. The selection gradient on asymptotic population growth rate, a measure of mean absolute fitness, was weaker than selection on first‐year survival and recruitment probabilities. Weaker selection on mean fitness occurs because weaning mass has little impact on adult survival, the fitness component to which the population growth of long‐lived species is most sensitive. These results highlight the need to interpret individual variation in phenotypic traits in a context that considers the demographic pathways between the trait and an inclusive proxy of individual fitness. Although variation in weaning mass do not translate to permanent survival differences among individuals in adulthood, it explains heterogeneity and positive covariation between survival and breeding in early life, which contribute to between‐individual variation in fitness.     

Very low levels of direct additive genetic variance in fitness and fitness components in a red squirrel population

A trait must genetically correlate with fitness in order to evolve in response to natural selection, but theory suggests that strong directional selection should erode additive genetic variance in fitness and limit future evolutionary potential. Balancing selection has been proposed as a mechanism that could maintain genetic variance if fitness components trade off with one another and has been invoked to account for empirical observations of higher levels of additive genetic variance in fitness components than would be expected from mutation–selection balance. Here, we used a long‐term study of an individually marked population of North American red squirrels (Tamiasciurus hudsonicus) to look for evidence of (1) additive genetic variance in lifetime reproductive success and (2) fitness trade‐offs between fitness components, such as male and female fitness or fitness in high‐ and low‐resource environments. “Animal model” analyses of a multigenerational pedigree revealed modest maternal effects on fitness, but very low levels of additive genetic variance in lifetime reproductive success overall as well as fitness measures within each sex and environment. It therefore appears that there are very low levels of direct genetic variance in fitness and fitness components in red squirrels to facilitate contemporary adaptation in this population.     

Population density mediates the interaction between pre‐ and postmating sexual selection

When females mate with more than one male, sexual selection acts both before and after mating. The interaction between pre‐ and postmating episodes of selection is expected to be context dependent, but few studies have investigated how total sexual selection changes under different ecological conditions. We examined how population density mediates the interaction between pre‐ and postmating sexual selection by establishing replicate populations of the horned dung beetle Onthophagus taurus at low, medium, and high densities, and then using microsatellite‐based parentage analyses to measure male fitness. We found that mating success and fertilization success were positively correlated at all three densities, but the strength of the correlation decreased with increasing density. We also found a shift from negative to positive linear selection on testes mass as density increased, and opposing selection on weapons and testes at high densities. These patterns suggest that the importance of postmating processes increases with increasing population density, which reduces the selective advantage of weapons for premating contest competition, and increases the selective advantage of large ejaculates for postmating sperm competition. We expect that density‐dependent selection on testes mass has contributed to the phenotypic variation observed between natural populations of O. taurus that differ in density.     

Among‐year variation in selection during early life stages and the genetic basis of fitness in Arabidopsis thaliana

Incomplete information regarding both selection regimes and the genetic basis of fitness limits our understanding of adaptive evolution. Among‐year variation in the genetic basis of fitness is rarely quantified, and estimates of selection are typically based on single components of fitness, thus potentially missing conflicting selection acting during other life‐history stages. Here, we examined among‐year variation in selection on a key life‐history trait and the genetic basis of fitness covering the whole life cycle in the annual plant Arabidopsis thaliana. We planted freshly matured seeds of >200 recombinant inbred lines (RILs) derived from a cross between two locally adapted populations (Italy and Sweden), and both parental genotypes at the native site of the Swedish population in three consecutive years. We quantified selection against the nonlocal Italian genotype, mapped quantitative trait loci (QTL) for fitness and its components, and quantified selection on timing of germination during different life stages. In all 3 years, the local Swedish genotype outperformed the nonlocal Italian genotype. However, both the contribution of early life stages to relative fitness, and the effects of fitness QTL varied among years. Timing of germination was under conflicting selection through seedling establishment vs. adult survival and fecundity, and both the direction and magnitude of net selection varied among years. Our results demonstrate that selection during early life stages and the genetic basis of fitness can vary markedly among years, emphasizing the need for multiyear studies considering the whole life cycle for a full understanding of natural selection and mechanisms maintaining local adaptation.     

Woodland caribou habitat selection patterns in relation to predation risk and forage abundance depend on reproductive state

The ideal free distribution assumes that animals select habitats that are beneficial to their fitness. When the needs of dependent offspring differ from those of the parent, ideal habitat selection patterns could vary with the presence or absence of offspring. We test whether habitat selection depends on reproductive state due to top‐down or bottom‐up influences on the fitness of woodland caribou (Rangifer tarandus caribou), a threatened, wide‐ranging herbivore. We combined established methods of fitting resource and step selection functions derived from locations of collared animals in Ontario with newer techniques, including identifying calf status from video collar footage and seasonal habitat selection analysis through latent selection difference functions. We found that females with calves avoided predation risk and proximity to roads more strongly than females without calves within their seasonal ranges. At the local scale, females with calves avoided predation more strongly than females without calves. Females with calves increased predation avoidance but not selection for food availability upon calving, whereas females without calves increased selection for food availability across the same season. These behavioral responses suggest that habitat selection by woodland caribou is influenced by reproductive state, such that females with calves at heel use habitat selection to offset the increased vulnerability of their offspring to predation risk.     

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.     

Dynamic changes in host–virus interactions associated with colony founding and social environment in fire ant queens (Solenopsis invicta)

The dynamics of host – parasite interactions can change dramatically over the course of a chronic infection as the internal (physiological) and external (environmental) conditions of the host change. When queens of social insects found a colony, they experience changes in both their physiological state (they develop their ovaries and begin laying eggs) and the social environment (they suddenly stop interacting with the other members of the mother colony), making this an excellent model system for examining how these factors interact with chronic infections. We investigated the dynamics of host – viral interactions in queens of Solenopsis invicta (fire ant) as they transition from mating to colony founding/brood rearing to the emergence of the first workers. We examined these dynamics in naturally infected queens in two different social environments, where queens either founded colonies as individuals or as pairs. We hypothesized that stress associated with colony founding plays an important role in the dynamics of host – parasite interactions. We also hypothesized that different viruses have different modalities of interaction with the host that can be quantified by physiological measures and genomic analysis of gene expression in the host. We found that the two most prevalent viruses, SINV‐1 and SINV‐2, are associated with different fitness costs that are mirrored by different patterns of gene expression in the host. In fact SINV‐2, the virus that imposes the significant reduction of a queen's reproductive output is also associated with larger changes of global gene expression in the host. These results show the complexity of interactions between S. invicta and two viral parasites. Our findings also show that chronic infections by viral parasites in insects are dynamic processes that may pose different challenges in the host, laying the groundwork for interesting ecological and evolutionary considerations.