Socially flexible female choice differs among populations of the Pacific field cricket: geographical variation in the interaction coefficient psi (Ψ)

Indirect genetic effects (IGEs) occur when genes expressed in one individual affect the phenotype of a conspecific. Theoretical models indicate that the evolutionary consequences of IGEs critically depend on the genetic architecture of interacting traits, and on the strength and direction of phenotypic effects arising from social interactions, which can be quantified by the interaction coefficient Ψ. In the context of sexually selected traits, strong positive Ψ tends to exaggerate evolutionary change, whereas negative Ψ impedes sexual trait elaboration. Despite its theoretical importance, whether and how Ψ varies among geographically distinct populations is unknown. Such information is necessary to evaluate the potential for IGEs to contribute to divergence among isolated or semi-isolated populations. Here, we report substantial variation in Ψ for a behavioural trait involved in sexual selection in the field cricket Teleogryllus oceanicus: female choosiness. Both the strength and direction of Ψ varied among geographically isolated populations. Ψ also changed over time. In a contemporary population of crickets from Kauai, experience of male song increased female choosiness. In contrast, experience of male song decreased choosiness in an ancestral population from the same location. This rapid change corroborates studies examining the evolvability of Ψ and demonstrates how interpopulation variation in the interaction coefficient might influence sexual selection and accelerate divergence of traits influenced by IGEs that contribute to reproductive isolation in nascent species or subspecies.     

Evolution of animal genitalia: morphological correlates of fitness components in a water strider

Rapid divergence of male genitalia is one of the most general evolutionary trends in animals with internal fertilization, but the mechanisms of genital evolution are poorly understood. The current study represents the first comprehensive attempt to test the main hypotheses that have been suggested to account for genital evolution (the lock-and-key, sexual selection and pleiotropy hypotheses) with intraspecific data. We measure multivariate phenotypic selection in a water strider species, by relating five different components of fitness (mating frequency, fecundity, egg hatching rate, offspring survival rate and offspring growth rate) to a suite of genital and non-genital morphological traits (in total 48). Body size had a series of direct effects in both sexes. Large size in females was positively related to both fecundity and egg hatching rate. There was positive sexual selection for large size in males (mating frequency), which to some extent was offset by a reduced number of eggs laid by females mated to large males. Male genitalic morphology influenced male mating frequency, but the detected directional selection on genitalia was due to indirect selection on phenotypically correlated non-intromittent traits. Further, we found no assortative mating between male intromittent genitalia and female morphology. Neither did we find any indications of male genitalia conveying information of male genetic quality. Several new insights can be gained from our study. Most importantly, our results are in stark disagreement with the long standing lock-and-key hypothesis of genital evolution, as well as with certain models of sexual selection. Our results are, however, in agreement with other models of sexual selection as well as with the pleiotropy hypothesis of genital evolution. Fluctuating asymmetry of bilaterally symmetrical traits, genital as well as non-genital, had few effects on fitness. Females with low fluctuating asymmetry in leg length produced offspring with a higher survival rate, a pattern most proba bly caused by direct phenotypic maternal effects. We also discuss the relevance of our results to sexual conflict over mating, and the evolution of sexual traits by coevolutionary arms races between the sexes.     

The quantitative genetics of indirect genetic effects: a selective review of modelling issues

Indirect genetic effects (IGE) occur when the genotype of an individual affects the phenotypic trait value of another conspecific individual. IGEs can have profound effects on both the magnitude and the direction of response to selection. Models of inheritance and response to selection in traits subject to IGEs have been developed within two frameworks; a trait-based framework in which IGEs are specified as a direct consequence of individual trait values, and a variance-component framework in which phenotypic variance is decomposed into a direct and an indirect additive genetic component. This work is a selective review of the quantitative genetics of traits affected by IGEs, with a focus on modelling, estimation and interpretation issues. It includes a discussion on variance-component vs trait-based models of IGEs, a review of issues related to the estimation of IGEs from field data, including the estimation of the interaction coefficient Ψ (psi), and a discussion on the relevance of IGEs for response to selection in cases where the strength of interaction varies among pairs of individuals. An investigation of the trait-based model shows that the interaction coefficient Ψ may deviate considerably from the corresponding regression coefficient when feedback occurs. The increasing research effort devoted to IGEs suggests that they are a widespread phenomenon, probably particularly in natural populations and plants. Further work in this field should considerably broaden our understanding of the quantitative genetics of inheritance and response to selection in relation to the social organisation of populations.     

Phenotypic plasticity in female mate choice behavior is mediated by an interaction of direct and indirect genetic effects in Drosophila melanogaster

Female mate choice is a complex decision‐making process that involves many context‐dependent factors. In Drosophila melanogaster, a model species for the study of sexual selection, indirect genetic effects (IGEs) of general social interactions can influence female mate choice behaviors, but the potential impacts of IGEs associated with mating experiences are poorly understood. Here, we examined whether the IGEs associated with a previous mating experience had an effect on subsequent female mate choice behaviors and quantified the degree of additive genetic variation associated with this effect. Females from 21 different genetic backgrounds were housed with males from one of two distinct genetic backgrounds for either a short (3 hr) or long (48 hr) exposure period and their subsequent mate choice behaviors were scored. We found that the genetic identity of a previous mate significantly influenced a female's subsequent interest in males and preference of males. Additionally, a hemiclonal analysis revealed significant additive genetic variation associated with experience‐dependent mate choice behaviors, indicating a genotype‐by‐environment interaction for both of these parameters. We discuss the significance of these results with regard to the evolution of plasticity in female mate choice behaviors and the maintenance of variation in harmful male traits.     

EXPERIMENTAL EVIDENCE FOR THE EVOLUTION OF INDIRECT GENETIC EFFECTS: CHANGES IN THE INTERACTION EFFECT COEFFICIENT,PSI (Ψ), DUE TO SEXUAL SELECTION

Indirect genetics effects (IGEs)—when the genotype of one individual affects the phenotypic expression of a trait in another—may alter evolutionary trajectories beyond that predicted by standard quantitative genetic theory as a consequence of genotypic evolution of the social environment. For IGEs to occur, the trait of interest must respond to one or more indicator traits in interacting conspecifics. In quantitative genetic models of IGEs, these responses (reaction norms) are termed interaction effect coefficients and are represented by the parameter psi (Ψ). The extent to which Ψ exhibits genetic variation within a population, and may therefore itself evolve, is unknown. Using an experimental evolution approach, we provide evidence for a genetic basis to the phenotypic response caused by IGEs on sexual display traits in Drosophila serrata. We show that evolution of the response is affected by sexual but not natural selection when flies adapt to a novel environment. Our results indicate a further mechanism by which IGEs can alter evolutionary trajectories—the evolution of interaction effects themselves.     

Genetic and environmental sources of egg size, fecundity and body size in the migrant skipper, Parnara guttata guttata (Lepidoptera: Hesperiidae)

Genetic and environmental sources of egg size, fecundity and body size (forewing length) were examined in the butterfly, Parnara guttata guttata. Phenotypic and genetic correlation and heritability were estimated for these traits under different day-length and temperature conditions. Egg size and fecundity had relatively high heritabilities, and body sizes in males and females had moderate and high heritability, respectively. Negative phenotypic and genetic correlations between egg size and fecundity were estimated in treatments corresponding to the natural conditions during larval development of the first and second generations. Phenotypic and genetic correlations between body size and egg size differed considerably between insects reared under long and short day-lengths. Next, genotype–environment interactions were estimated by comparing reaction norms to day-length or temperature of these traits among families. ANOVA analysis revealed significant genotype–environment interactions in egg size and forewing length in both sexes for day-length and temperature. These results suggested that a large additive genetic variance for egg size might have been maintained by a genetic trade-off and/or by genotype–environment interactions in P. g. guttata.     

Maternal and paternal contributions to egg size and egg number variation in the blackfin pearl killifish <Emphasis Type="Italic">Austrolebias nigripinnis</Emphasis>

Reproductive effort, egg number and egg size are traditionally considered to be ‘female’ life history traits. However, females often adjust the amount of resources allocated to reproduction depending on their mate, causing male environmental effects on life history traits. If females respond to male traits which are genetically variable, then male environmental effects contain indirect genetic effects. Estimates of how much of the total variation in life history traits originates from female effects versus male environmental effects, seems mostly lacking. We have investigated variation in rates of egg production and in egg size in the annual Argentinian blackfin pearl killifish Austrolebias nigripinnis, in a crossed design where males were exchanged repeatedly between females. Our analysis of phenotypic variance components of reproductive effort, egg size and egg number indicates that the amount of variation contributed by male environmental effects is equal (egg size, reproductive effort) or larger (egg number) than that between females. For egg size and number, we find that male environmental effects consist of a male random effect representing the average response of females to male phenotype, plus a female-male interaction term. This term can be understood as the deviation from the population mean of an individual female’s response. For reproductive effort, we find that the male environmental effect consists of an interaction term only. Random effects on egg size and number additionally vary in magnitude depending on the weekday where we collected eggs, probably due to cyclic variation in experimental conditions. Since we find that both male phenotype and environmental conditions affect egg size and number as determined by females, our results suggest that selection on these life history traits will be frequency-dependent.     

Indirect genetic effects and sexual conflicts: Partner genotype influences multiple morphological and behavioral reproductive traits in a flatworm

The expression of an individual's phenotypic traits can be influenced by genes expressed in its social partners. Theoretical models predict that such indirect genetic effects (IGEs) on reproductive traits should play an important role in determining the evolutionary outcome of sexual conflict. However, empirical tests of (i) whether reproductive IGEs exist, (ii) how they vary among genotypes, and (iii) whether they are uniform for different types of reproductive traits are largely lacking. We addressed this in a series of experiments in the simultaneously hermaphroditic flatworm Macrostomum lignano. We found strong evidence for IGEs on both morphological and behavioral reproductive traits. Partner genotype had a significant impact on the testis size of focal individuals—varying up to 2.4‐fold—suggesting that IGEs could mediate sexual conflicts that target the male sex function. We also found that time to first copulation was affected by a genotype × genotype interaction between mating partners, and that partner genotype affected the propensity to copulate and perform the postcopulatory suck behavior, which may mediate conflicts over the fate of received ejaculate components. These findings provide clear empirical evidence for IGEs on multiple behavioral and morphological reproductive traits, which suggests that the evolutionary dynamics of these traits could be altered by genes contained in the social environment.     

Good genes, genetic compatibility and the evolution of polyandry: use of the diallel cross to address competing hypotheses

Genetic benefits can enhance the fitness of polyandrous females through the high intrinsic genetic quality of females' mates or through the interaction between female and male genes. I used a full diallel cross, a quantitative genetics design that involves all possible crosses among a set of genetically homogeneous lines, to determine the mechanism through which polyandrous female decorated crickets (Gryllodes sigillatus) obtain genetic benefits. I measured several traits related to fitness and partitioned the phenotypic variance into components representing the contribution of additive genetic variance ('good genes'), nonadditive genetic variance (genetic compatibility), as well as maternal and paternal effects. The results reveal a significant variance attributable to both nonadditive and additive sources in the measured traits, and their influence depended on which trait was considered. The lack of congruence in sources of phenotypic variance among these fitness-related traits suggests that the evolution and maintenance of polyandry are unlikely to have resulted from one selective influence, but rather are the result of the collective effects of a number of factors.     

THE GENETIC BASIS OF LIFE-HISTORY CHARACTERS IN A POLYCHAETE EXHIBITING PLANKTOTROPHY AND LECITHOTROPHY

The polychaete Streblospio benedicti is unusual in that several field populations consist of individuals that exhibit either planktotrophic or lecithotrophic larval development. Planktotrophy in this species involves production of many small ova that develop into feeding larvae with a two- to three-week planktonic period. Lecithotrophy involves production of fewer, larger ova that develop into nonfeeding larvae that are brooded longer and have a brief planktonic stage. Reciprocal matings were performed to investigate genetic variance components and the correlation structure of life-history traits associated with planktotrophy and lecithotrophy. Our objective was to better understand persistence of this developmental dichotomy in Streblospio benedicti, and among marine invertebrates in general. Substantial additive genetic variation (75–98% of total) was detected for the following characters at first reproduction: female length; position of the first gametogenic setiger and first brood pouch; ovum diameter; three traits related to fecundity (numbers of ova per ovary, larvae per brood pouch, and larvae per brood); median larval planktonic period and the presence of larval swimming setae; but not for total number of brood pouches; larval length; larval feeding; and larval survivorship. Based on the unusual geographic distribution of development modes in this species, we hypothesize that the developmental traits have evolved in allopatry and have only recently come into contact in North Carolina. The high additive contribution to variance observed for many traits may be inflated due to (a) nonrandom breeding in nature, and (b) examination of only one component of an age-structured population at one time. Nuclear interaction variance and maternal variance accounted for 84% of the total variation in larval survivorship. This observation supports other empirical studies and theoretical predictions that nonadditive components of variance will increase in importance in individual traits that are most closely tied to fitness. A network of life-history trait correlations was observed that defines distinct planktotrophic and lecithotrophic trait complexes. Negative genetic correlations were present between fecundity and egg size, between fecundity and position of the first gametes, and between larval survivorship and median planktonic period. Positive genetic correlations were detected between fecundity and female size at first reproduction and between planktonic period and the presence of swimming setae. Intergenerational product-moment correlations were negative for larval length and fecundity, planktonic period and egg size, female size and larval survivorship, and fecundity and larval survivorship. If the genetic correlation structure observed in the laboratory persists in the field, it may constrain responses of individual characters to directional selection, and indirectly perpetuate the dichotomies associated with planktotrophy and lecithotrophy.     

HOW TO MEASURE INDIRECT GENETIC EFFECTS: THE CONGRUENCE OF TRAIT-BASED AND VARIANCE-PARTITIONING APPROACHES

Indirect genetic effects (IGEs), which occur when phenotypic expression in one individual is influenced by genes in another conspecific individual, may have a drastic effect on evolutionary response to selection. General evolutionary models of IGEs have been developed using two distinct theoretical frameworks derived from maternal effects theory. The first framework is trait-based and focuses on how phenotypes are influenced by specific traits in a social partner, with the strength of interactions defined by the matrix Ψ. The second framework partitions total genetic variance into components representing direct effects, indirect effects, and the covariance between them, without identifying specific social traits responsible for IGEs. The latter framework has been employed more commonly by empiricists because the methods for estimating variance components are relatively straightforward. Here, we show how these two theoretical frameworks are related to each other and derive equations that can be used to translate between them. This translation leads to a generalized method that can be used to estimate Ψ via standard quantitative genetic breeding designs or pedigrees from natural populations. This method can be used in a very general set of circumstances and is widely applicable to all IGEs, including maternal effects and other interactions among relatives.     

Bottleneck Effects on Genetic Variance for Courtship Repertoire

Bottleneck effects on evolutionary potential in mating behavior were addressed through assays of additive genetic variances and resulting phenotypic responses to drift in the courtship repertoires of six two-pair founder-flush lines and two control populations of the housefly. A simulation addressed the complication that an estimate of the genetic variance for a courtship trait (e.g., male performance vigor or the female requirement for copulation) must involve assays against the background behavior of the mating partners. The additive ``environmental' effect of the mating partner's phenotype simply dilutes the net parent-offspring covariance for a trait. However, if there is an interaction with this ``environmental' component, negative parent-offspring covariances can result under conditions of high incompatibility between the population's distributions for male performance and female choice requirements, despite high levels of genetic variance. All six bottlenecked lines exhibited significant differentiation from the controls in at least one measure of the parent-offspring covariance for male performance or female choice (estimated by 50 parent-son and 50 parent-daughter covariances for 10 courtship traits per line) which translated to significant phenotypic drift. However, the average effect across traits or across lines did not yield a significant net increase in genetic variance due to bottlenecks. Concerted phenotypic differentiation due to the founder-flush event provided indirect evidence of directional dominance in a subset of traits. Furthermore, indirect evidence of genotype-environment interactions (potentially producing genotype-genotype effects) was found in the negative parent-offspring covariances predicted by the male-female interaction simulation and by the association of the magnitude of phenotypic drift with the absolute value of the parent-offspring covariance. Hence, nonadditive genetic effects on mating behavior may be important in structuring genetic variance for courtship, although most of the increases in genetic variance would be expected to reflect inbreeding depression with relatively rare situations representing the facilitation of speciation by bottlenecks.     

Indirect genetic effects and the evolution of aggression in a vertebrate system

Aggressive behaviours are necessarily expressed in a social context, such that individuals may be influenced by the phenotypes, and potentially the genotypes, of their social partners. Consequently, it has been hypothesized that indirect genetic effects (IGEs) arising from the social environment will provide a major source of heritable variation on which selection can act. However, there has been little empirical scrutiny of this to date. Here we test this hypothesis in an experimental population of deer mice (Peromyscus maniculatus). Using quantitative genetic models of five aggression traits, we find repeatable and heritable differences in agonistic behaviours of focal individuals when presented with an opponent mouse. For three of the traits, there is also support for the presence of IGEs, and estimated correlations between direct and indirect genetic (rAO,F) effects were high. As a consequence, any selection for aggression in the focal individuals should cause evolution of the social environment as a correlated response. In two traits, strong positive rAO,F will cause the rapid evolution of aggression, while in a third case changes in the phenotypic mean will be constrained by negative covariance between direct and IGEs. Our results illustrate how classical analyses may miss important components of heritable variation, and show that a full understanding of evolutionary dynamics requires explicit consideration of the genetic component of the social environment.     

Genetic variation in male effects on female reproduction and the genetic covariance between the sexes

Abstract.— Males of many insect species increase the fecundity and/or egg size of their mates through the amount or composition of their nuptial gifts or ejaculate. The genetic bases of such male effects on fecundity or egg size are generally unknown, and thus their ability to evolve remains speculative. Likewise, the genetic relationship between male and female investment into reproduction in dioecious species, which is expected to be positive if effects on fecundity are controlled by at least some of the same genes in males and females, is also unknown. Males of the seed beetle Stator limbatus contribute large ejaculates to females during mating, and the amount of donated ejaculate is positively correlated with male body mass. Females mated to large males lay more eggs in their lifetime than females mated to small males. We describe an experiment in which we quantify genetic variation in the number of eggs sired by males (mated to a single female) and found that a significant proportion of the phenotypic variance in the number of eggs sired by males was explained by their genotype. Additionally, the number of eggs sired by a male was highly positively genetically correlated with his body mass. The between-sex genetic correlation, that is, the genetic correlation between the number of eggs sired by males and the number of eggs laid by females, was highly positive when eggs were laid on Acacia greggii seeds. This indicates that males that sire many eggs have sisters that lay many eggs. Thus, some of the genes that control male ejaculate size (or some other fecundity-enhancing factor) when expressed in males appear to control fecundity when expressed in females. We found no significant interaction between male and female genotype on fecundity.     

Evolutionary ecology of egg size and number in a seed beetle: genetic trade-off differs between environments

Abstract In many organisms, large offspring have improved fitness over small offspring, and thus their size is under strong selection. However, due to a trade-off between offspring size and number, females producing larger offspring necessarily must produce fewer unless the total amount of reproductive effort is unlimited. Because differential gene expression among environments may affect genetic covariances among traits, it is important to consider environmental effects on the genetic relationships among traits. We compared the genetic relationships among egg size, lifetime fecundity, and female adult body mass (a trait linked to reproductive effort) in the seed beetle, Stator limbatus , between two environments (host-plant species Acacia greggii and Cercidium floridum ). Genetic correlations among these traits were estimated through half-sib analysis, followed with artificial selection on egg size to observe the correlated responses of lifetime fecundity and female body mass. We found that the magnitude of the genetic trade-off between egg size and lifetime fecundity differed between environments–a strong trade-off was estimated when females laid eggs on C. floridum seeds, yet this trade-off was weak when females laid eggs on A. greggii seeds. Also differing between environments was the genetic correlation between egg size and female body mass–these traits were positively genetically correlated for egg size on A. greggii seeds, yet uncorrelated on C. floridum seeds. On A. greggii seeds, the evolution of egg size and traits linked to reproductive effort (such as female body mass) are not independent from each other as commonly assumed in life-history theory.     

Genetic basis of variation in morphological and life-history traits of a wild population of pink salmon

Understanding the genetic basis of phenotypic variation is essential for predicting the direction and rate of phenotypic evolution. We estimated heritabilities and genetic correlations of morphological (fork length, pectoral and pelvic fin ray counts, and gill arch raker counts) and life-history (egg number and individual egg weight) traits of pink salmon (Oncorhynchus gorbuscha) from Likes Creek, Alaska, in order to characterize the genetic basis of phenotypic variation in this species. Families were created from wild-caught adults, raised to the fry stage in the lab, released into the wild, and caught as returning adults and assigned to families using microsatellite loci and a growth hormone locus. Morphological traits were all moderately to highly heritable, but egg number and egg weight were not heritable, suggesting that past selection has eliminated additive genetic variation in egg number and egg weight or that there is high environmental variance in these traits. Genetic correlations were similar for nonadjacent morphological traits and adjacent traits. Genetic correlations predicted phenotypic correlations fairly accurately, but some pairs of traits with low genetic correlations had high phenotypic correlations, and vice versa, emphasizing the need to use caution when using phenotypic correlations as indices of genetic correlations. This is one of only a handful of studies to estimate heritabilities and genetic correlations for a wild population.     

Female crickets trade offspring viability for fecundity

A growing number of studies are suggesting that females can improve the viability of their embryos by mating with multiple males. However, the reason why females should have low rates of embryo viability is puzzling. Here we conduct a quantitative genetic study of maternal effects on embryo viability in the field cricket Teleogryllus oceanicus. After controlling for female body size, we find significant additive genetic variance for ovary weight, a measure of fecundity, and egg hatching success, a measure of embryo viability. Moreover, we show a genetic trade-off between these traits that is predicted from life-history theory. High rates of embryo mortality in this highly fecund species might therefore be explained by selection favouring an optimum balance between fecundity and embryo viability that maximizes maternal fitness. Paternal effects on female fecundity and embryo viability are often seen as benefits driving the evolution of polyandrous behaviour. However, we raise the alternative possibility that paternal effects might shift females from their naturally selected optimum, and present some support for the notion that sexual conflict over a female's optimal fecundity and embryo viability might generate antagonistic coevolution between the sexes.     

Effects of pathogen exposure on life‐history variation in the gypsy moth (Lymantria dispar)

Investment in host defences against pathogens may lead to trade‐offs with host fecundity. When such trade‐offs arise from genetic correlations, rates of phenotypic change by natural selection may be affected. However, genetic correlations between host survival and fecundity are rarely quantified. To understand trade‐offs between immune responses to baculovirus exposure and fecundity in the gypsy moth (Lymantria dispar), we estimated genetic correlations between survival probability and traits related to fecundity, such as pupal weight. In addition, we tested whether different virus isolates have different effects on male and female pupal weight. To estimate genetic correlations, we exposed individuals of known relatedness to a single baculovirus isolate. To then evaluate the effect of virus isolate on pupal weight, we exposed a single gypsy moth strain to 16 baculovirus isolates. We found a negative genetic correlation between survival and pupal weight. In addition, virus exposure caused late‐pupating females to be identical in weight to males, whereas unexposed females were 2–3 times as large as unexposed males. Finally, we found that female pupal weight is a quadratic function of host mortality across virus isolates, which is likely due to trade‐offs and compensatory growth processes acting at high and low mortality levels, respectively. Overall, our results suggest that fecundity costs may strongly affect the response to selection for disease resistance. In nature, baculoviruses contribute to the regulation of gypsy moth outbreaks, as pathogens often do in forest‐defoliating insects. We therefore argue that trade‐offs between host life‐history traits may help explain outbreak dynamics.     

The role of indirect genetic effects in the evolution of interacting reproductive behaviors in the burying beetle,Nicrophorus vespilloides

Social interactions can give rise to indirect genetic effects (IGEs), which occur when genes expressed in one individual affect the phenotype of another individual. The evolutionary dynamics of traits can be altered when there are IGEs. Sex often involves indirect effects arising from first‐order (current) or second‐order (prior) social interactions, yet IGEs are infrequently quantified for reproductive behaviors. Here, we use experimental populations of burying beetles that have experienced bidirectional selection on mating rate to test for social plasticity and IGEs associated with focal males mating with a female either without (first‐order effect) or with (second‐order effect) prior exposure to a competitor, and resource defense behavior (first‐order effect). Additive IGEs were detected for mating rate arising from (first‐order) interactions with females. For resource defense behavior, a standard variance partitioning analysis provided no evidence of additive genetic variance—either direct or indirect. However, behavior was predicted by focal size relative to that of the competitor, and size is also heritable. Assuming that behavior is causally dependent on relative size, this implies that both DGEs and IGEs do occur (and may potentially interact). The relative contribution of IGEs may differ among social behaviors related to mating which has consequences for the evolutionary trajectories of multivariate traits.     

On indirect genetic effects in structured populations

Indirect genetic effects (IGEs) occur when the phenotype of an individual, and possibly its fitness, depends, at least in part, on the genes of its social partners. The effective result is that environmental sources of phenotypic variance can themselves evolve. Simple models have shown that IGEs can alter the rate and direction of evolution for traits involved in interactions. Here we expand the applicability of the theory of IGEs to evolution in metapopulations by including nonlinear interactions between individuals and population genetic structure. Although population subdivision alone generates some dramatic and nonintuitive evolutionary dynamics for interacting phenotypes, the combination of nonlinear interactions with subdivision reveals an even greater importance of IGEs. The presence of genetic structure links the evolution of interacting phenotypes and the traits that influence their expression ("effector traits") even in the absence of genetic correlations. When nonlinear social effects occur in subdivided populations, evolutionary response is altered and can even oppose the direction expected due to direct selection. Because population genetic structure allows for multilevel selection, we also investigate the role of IGEs in determining the response to individual and group selection. We find that nonlinear social effects can cause interference between levels of selection even when they act in the same direction. In some cases, interference can be so extreme that the actual evolutionary response to multilevel selection is opposite in direction to that predicted by summing selection at each level. This theoretical result confirms empirical data that show higher levels of selection cannot be ignored even when selection acts in the same direction at all levels.