Carry-over effects of seasonal migration on reproductive success through breeding site retention in a partially migratory bird.

Understanding the maintenance and dynamics of phenotypic polymorphisms requires unpicking key ecological mechanisms shaping the fitness costs and benefits of expressing alternative phenotypes, generating selection. Seasonal migration versus year-round residence expressed in partially migratory populations represents one common polymorphism that can experience strong selection through differential reproductive success. Yet, key hypothesised pathways that could generate such selection remain to be empirically tested. One hypothesis is that migratory tactics affect subsequent reproductive success through carry-over effects on breeding site retention and resulting breeding dispersal. By remaining in breeding areas all year round, residents could retain their preferred breeding site between years, and consequently have higher reproductive success. Conversely, migrants that escape harsh non-breeding season conditions could return in better condition, with high resource holding potential, and outcompete residents to retain their site. Such effects could further depend on migration timing and vary between years. Yet, such pathways have not been quantified, precluding empirical parameterisation of partial migration theory. We used 4 years of breeding and non-breeding season data from partially migratory European shags (Gulosus aristotelis) to test whether the three most frequent migratory tactics in this population (full resident, early migrant departing soon after breeding, and late migrant departing in late autumn) differed in their breeding site retention; whether site retention predicted reproductive success; and hence whether effects of migratory tactic on reproductive success were explicable through site retention. Overall, residents were much more likely to retain their breeding site between years than both early and late migrants, and site retention was associated with increased reproductive success. Yet, these effects varied somewhat among years: late migrants were always least likely to retain their site but had variable relative reproductive success. Path analyses revealed that effects of migratory tactic on reproductive success were only partly attributable to breeding site retention. These results indicate that multiple mechanisms underlie reproductive selection on migratory tactics, potentially contributing to maintaining behavioural polymorphisms. Yet, the clear associations between migratory tactics and local breeding dispersal reveal that these movements can be strongly interlinked across seasons, shaping overall spatioseasonal dynamics in partially migratory systems.

Multigenerational Fitness Effects of Natural Immigration Indicate Strong Heterosis and Epistatic Breakdown in a Wild Bird Population.

AbstractThe fitness of immigrants and their descendants produced within recipient populations fundamentally underpins the genetic and population dynamic consequences of immigration. Immigrants can in principle induce contrasting genetic effects on fitness across generations, reflecting multifaceted additive, dominance, and epistatic effects. Yet full multigenerational and sex-specific fitness effects of regular immigration have not been quantified within naturally structured systems, precluding inference on underlying genetic architectures and population outcomes. We used four decades of song sparrow (Melospiza melodia) life history and pedigree data to quantify fitness of natural immigrants, natives, and their F1, F2, and backcross descendants and test for evidence of nonadditive genetic effects. Values of key fitness components (including adult lifetime reproductive success and zygote survival) of F1 offspring of immigrant-native matings substantially exceeded their parent mean, indicating strong heterosis. Meanwhile, F2 offspring of F1-F1 matings had notably low values, indicating surprisingly strong epistatic breakdown. Furthermore, magnitudes of effects varied among fitness components and differed between female and male descendants. These results demonstrate that strong nonadditive genetic effects on fitness can arise within weakly structured and fragmented populations experiencing frequent natural immigration. Such effects will substantially affect the net degree of effective gene flow and resulting local genetic introgression and adaptation.

Multi-generation genetic contributions of immigrants reveal cryptic elevated and sex-biased effective gene flow within a natural meta-population.

Impacts of immigration on micro-evolution and population dynamics fundamentally depend on net rates and forms of resulting gene flow into recipient populations. Yet, the degrees to which observed rates and sex ratios of physical immigration translate into multi-generational genetic legacies have not been explicitly quantified in natural meta-populations, precluding inference on how movements translate into effective gene flow and eco-evolutionary outcomes. Our analyses of three decades of complete song sparrow (Melospiza melodia) pedigree data show that multi-generational genetic contributions from regular natural immigrants substantially exceeded those from contemporary natives, consistent with heterosis-enhanced introgression. However, while contributions from female immigrants exceeded those from female natives by up to three-fold, male immigrants' lineages typically went locally extinct soon after arriving. Both the overall magnitude, and the degree of female bias, of effective gene flow therefore greatly exceeded those which would be inferred from observed physical arrivals, altering multiple eco-evolutionary implications of immigration.

Sexual selection and mate limitation shape evolution of species' range limits.

Understanding what processes shape the formation of species' geographic range limits is one central objective linking ecology and evolutionary biology. One potentially key process is sexual selection; yet, theory examining how sexual selection could shape eco-evolutionary dynamics in marginal populations is still lacking. In species with separate sexes, range limits could be shaped by limitations in encountering mates at low densities. Sexual selection could therefore modulate mate limitation and resulting extinction-colonization dynamics at range margins, through evolution of mate encounter ability and/or mate competition traits, and their demographic consequences. We use a spatially explicit eco-genetic model to reveal how different forms of sexual selection can variably affect emerging range limits. Larger ranges emerged when sexual selection acted exclusively on traits increasing mate encounter probability, thus reducing female's mate limitation toward the range margins. In contrast, sexual selection via mate competition narrowed range limits due to increased trait-dependent mortality in males and elevated mate limitation for females. When mate encounter coevolved with mate competition, their combined effects on range limits depended on the mating system (polygyny vs. monogamy). Our results demonstrate that evolution of species' ranges may be importantly shaped by feedbacks between sexual selection and spatial population demography and dynamics.

Environmental variance in male mating success modulates the positive versus negative impacts of sexual selection on genetic load.

Sexual selection on males is predicted to increase population fitness, and delay population extinction, when mating success negatively covaries with genetic load across individuals. However, such benefits of sexual selection could be counteracted by simultaneous increases in genome-wide drift resulting from reduced effective population size caused by increased variance in fitness. Resulting fixation of deleterious mutations could be greatest in small populations, and when environmental variation in mating traits partially decouples sexual selection from underlying genetic variation. The net consequences of sexual selection for genetic load and population persistence are therefore likely to be context dependent, but such variation has not been examined. We use a genetically explicit individual-based model to show that weak sexual selection can increase population persistence time compared to random mating. However, for stronger sexual selection such positive effects can be overturned by the detrimental effects of increased genome-wide drift. Furthermore, the relative strengths of mutation-purging and drift critically depend on the environmental variance in the male mating trait. Specifically, increasing environmental variance caused stronger sexual selection to elevate deleterious mutation fixation rate and mean selection coefficient, driving rapid accumulation of drift load and decreasing population persistence times. These results highlight an intricate balance between conflicting positive and negative consequences of sexual selection on genetic load, even in the absence of sexually antagonistic selection. They imply that environmental variances in key mating traits, and intrinsic genetic drift, should be properly factored into future theoretical and empirical studies of the evolution of population fitness under sexual selection.

Additive genetic and environmental variation interact to shape the dynamics of seasonal migration in a wild bird population.

Dissecting joint micro-evolutionary and plastic responses to environmental perturbations requires quantifying interacting components of genetic and environmental variation underlying expression of key traits. This ambition is particularly challenging for phenotypically discrete traits where multiscale decompositions are required to reveal nonlinear transformations of underlying genetic and environmental variation into phenotypic variation, and when effects must be estimated from incomplete field observations. We devised a joint multistate capture-recapture and quantitative genetic animal model, and fitted this model to full-annual-cycle resighting data from partially-migratory European shags (${Gulosus~{}aristotelis}$) to estimate key components of genetic, environmental and phenotypic variance in the ecologically critical discrete trait of seasonal migration versus residence. We demonstrate non-negligible additive genetic variance in latent liability for migration, resulting in detectable micro-evolutionary responses following two episodes of strong survival selection. Further, liability-scale additive genetic effects interacted with substantial permanent individual and temporary environmental effects to generate complex nonadditive effects on expressed phenotypes, causing substantial intrinsic gene-by-environment interaction variance on the phenotypic scale. Our analyses therefore reveal how temporal dynamics of partial seasonal migration arise from combinations of instantaneous micro-evolution and within-individual phenotypic consistency, and highlight how intrinsic phenotypic plasticity could expose genetic variation underlying discrete traits to complex forms of selection.

Immune genotypes, immune responses, and survival in a wild bird population.

Individuals vary in their immune genotype, inbreeding coefficient f, immune responses, survival to adulthood, and adult longevity. However, whether immune genes predict survival or longevity, whether such relationships are mediated through immune responses, and how f affects immune genotype remain unclear. We use a wild song sparrow (Melospiza melodia) population in which survival to adulthood, adult longevity, and f were measured precisely, and in which immune responses have previously been assessed. We investigate four toll-like receptor (TLR) and the major histocompatibility complex (MHC) class IIB exon 2 genes. We test whether immune genes predict fitness (survival to adulthood or adult longevity); whether immune genes predict immune response; whether immune response predicts fitness and whether fitness, immune responses, or immune genotypes are correlated with f. We find that survival to adulthood is not associated with immune gene variation, but adult longevity is decreased by high MHC allele diversity (especially in birds that were relatively outbred), and by the presence of a specific MHC supertype. Immune responses were affected by specific immune genotypes. Survival to adulthood and adult longevity were not predicted by immune response, implying caution in the use of immune response as a predictor for fitness. We also found no relationship between f and immune genotype. This finding indicates that immune gene associations with longevity and immune response are not artefacts of f, and suggests that pathogen-mediated selection at functional loci can slow the loss of genetic variation arising from genetic drift and small population size.

Hierarchical Variation in Phenotypic Flexibility across Timescales and Associated Survival Selection Shape the Dynamics of Partial Seasonal Migration.

AbstractPopulation responses to environmental variation ultimately depend on within-individual and among-individual variation in labile phenotypic traits that affect fitness and resulting episodes of selection. Yet complex patterns of individual phenotypic variation arising within and between time periods, as well as associated variation in selection, have not been fully conceptualized or quantified. We highlight how structured patterns of phenotypic variation in dichotomous threshold traits can theoretically arise and experience varying forms of selection, shaping overall phenotypic dynamics. We then fit novel multistate models to 10 years of band-resighting data from European shags to quantify phenotypic variation and selection in a key threshold trait underlying spatioseasonal population dynamics: seasonal migration versus residence. First, we demonstrate substantial among-individual variation alongside substantial between-year individual repeatability in within-year phenotypic variation ("flexibility"), with weak sexual dimorphism. Second, we demonstrate that between-year individual variation in within-year phenotypes ("supraflexibility") is structured and directional, consistent with the threshold trait model. Third, we demonstrate strong survival selection on within-year phenotypes-and hence on flexibility-that varies across years and sexes, including episodes of disruptive selection representing costs of flexibility. By quantitatively combining these results, we show how supraflexibility and survival selection on migratory flexibility jointly shape population-wide phenotypic dynamics of seasonal movement.

Eco-evolutionary extinction and recolonization dynamics reduce genetic load and increase time to extinction in highly inbred populations.

Understanding how genetic and ecological effects can interact to shape genetic loads within and across local populations is key to understanding ongoing persistence of systems that should otherwise be susceptible to extinction through mutational meltdown. Classic theory predicts short persistence times for metapopulations comprising small local populations with low connectivity, due to accumulation of deleterious mutations. Yet, some such systems have persisted over evolutionary time, implying the existence of mechanisms that allow metapopulations to avoid mutational meltdown. We first hypothesize a mechanism by which the combination of stochasticity in the numbers and types of mutations arising locally (genetic stochasticity), resulting local extinction, and recolonization through evolving dispersal facilitates metapopulation persistence. We then test this mechanism using a spatially and genetically explicit individual-based model. We show that genetic stochasticity in highly structured metapopulations can result in local extinctions, which can favor increased dispersal, thus allowing recolonization of empty habitat patches. This causes fluctuations in metapopulation size and transient gene flow, which reduces genetic load and increases metapopulation persistence over evolutionary time. Our suggested mechanism and simulation results provide an explanation for the conundrum presented by the continued persistence of highly structured populations with inbreeding mating systems that occur in diverse taxa.

Conceptualizing the evolutionary quantitative genetics of phenological life-history events: Breeding time as a plastic threshold trait.

Successfully predicting adaptive phenotypic responses to environmental changes, and predicting resulting population outcomes, requires that additive genetic (co)variances underlying microevolutionary and plastic responses of key traits are adequately estimated on appropriate quantitative scales. Such estimation in turn requires that focal traits, and their underlying quantitative genetic architectures, are appropriately conceptualized. Here, we highlight that directly analyzing observed phenotypes as continuously distributed quantitative traits can potentially generate biased and misleading estimates of additive genetic variances and individual-by-environment and gene-by-environment interactions, and hence of forms of plasticity and genetic constraints, if in fact the underlying biology is best conceptualized as an environmentally sensitive threshold trait. We illustrate this scenario with particular reference to the key phenological trait of seasonal breeding date, which has become a focus for quantifying joint microevolutionary, plastic, and population responses to environmental change, but has also become a focus for highlighting that predicted adaptive outcomes are not always observed. Specifically, we use simple simulations to illustrate how potentially misleading inferences on magnitudes of additive genetic variance, and forms of environmental interactions, can arise by directly analyzing observed breeding dates if the transition to breeding in fact represents a threshold trait with latent-scale plasticity. We summarize how existing and new datasets could be (re)analyzed, potentially providing new insights into how critical microevolutionary and plastic phenological responses to environmental variation and change can arise and be constrained.

Intrinsic emergence and modulation of sex-specific dominance reversals in threshold traits.

Sex-specific dominance reversals (SSDRs) in fitness-related traits, where heterozygotes' phenotypes resemble those of alternative homozygotes in females versus males, can simultaneously maintain genetic variation in fitness and resolve sexual conflict and thereby shape key evolutionary outcomes. However, the full implications of SSDRs will depend on how they arise and the resulting potential for evolutionary, ecological and environmental modulation. Recent field and laboratory studies have demonstrated SSDRs in threshold(-like) traits with dichotomous or competitive phenotypic outcomes, implying that such traits could promote the emergence of SSDRs. However, such possibilities have not been explicitly examined. I show how phenotypic SSDRs can readily emerge in threshold traits given genetic architectures involving large-effect loci alongside sexual dimorphism in the mean and variance in polygenic liability. I also show how multilocus SSDRs can arise in line-cross experiments, especially given competitive reproductive systems that generate nonlinear fitness outcomes. SSDRs can consequently emerge in threshold(-like) traits as functions of sexual antagonism, sexual dimorphism and reproductive systems, even with purely additive underlying genetic effects. Accordingly, I identify theoretical and empirical advances that are now required to discern the basis and occurrence of SSDRs in nature, probe forms of (co-)evolutionary, ecological and environmental modulation, and evaluate net impacts on sexual conflict.

Modelling the responses of partially migratory metapopulations to changing seasonal migration rates: From theory to data.

Among-individual and within-individual variation in expression of seasonal migration versus residence is widespread in nature and could substantially affect the dynamics of partially migratory metapopulations inhabiting seasonally and spatially structured environments. However, such variation has rarely been explicitly incorporated into metapopulation dynamic models for partially migratory systems. We, therefore, lack general frameworks that can identify how variable seasonal movements, and associated season- and location-specific vital rates, can control system persistence. We constructed a novel conceptual framework that captures full-annual-cycle dynamics and key dimensions of metapopulation structure for partially migratory species inhabiting seasonal environments. We conceptualize among-individual variation in seasonal migration as two variable vital rates: seasonal movement probability and associated movement survival probability. We conceptualize three levels of within-individual variation (i.e. plasticity), representing seasonal or annual variation in seasonal migration or lifelong fixed strategies. We formulate these concepts as a general matrix model, which is customizable for diverse life-histories and seasonal landscapes. To illustrate how variable seasonal migration can affect metapopulation growth rate, demographic structure and vital rate elasticities, we parameterize our general models for hypothetical short- and longer-lived species. Analyses illustrate that elasticities of seasonal movement probability and associated survival probability can sometimes equal or exceed those of vital rates typically understood to substantially influence metapopulation dynamics (i.e. seasonal survival probability or fecundity), that elasticities can vary non-linearly, and that metapopulation outcomes depend on the level of within-individual plasticity. We illustrate how our general framework can be applied to evaluate the consequences of variable and changing seasonal movement probability by parameterizing our models for a real partially migratory metapopulation of European shags Gulosus aristotelis assuming lifelong fixed strategies. Given observed conditions, metapopulation growth rate was most elastic to breeding season adult survival of the resident fraction in the dominant population. However, given doubled seasonal movement probability, variation in survival during movement would become the primary driver of metapopulation dynamics. Our general conceptual and matrix model frameworks, and illustrative analyses, thereby highlight complex ways in which structured variation in seasonal migration can influence dynamics of partially migratory metapopulations, and pave the way for diverse future theoretical and empirical advances.

Genetic variance in fitness indicates rapid contemporary adaptive evolution in wild animals.

The rate of adaptive evolution, the contribution of selection to genetic changes that increase mean fitness, is determined by the additive genetic variance in individual relative fitness. To date, there are few robust estimates of this parameter for natural populations, and it is therefore unclear whether adaptive evolution can play a meaningful role in short-term population dynamics. We developed and applied quantitative genetic methods to long-term datasets from 19 wild bird and mammal populations and found that, while estimates vary between populations, additive genetic variance in relative fitness is often substantial and, on average, twice that of previous estimates. We show that these rates of contemporary adaptive evolution can affect population dynamics and hence that natural selection has the potential to partly mitigate effects of current environmental change.

Adaptation to climate change through seasonal migration revealed by climatic versus demographic niche models.

Predicting the geographic range of species and their response to climatic variation and change are entwined goals in conservation and evolutionary ecology. Species distribution models (SDMs) are foundational in this effort and used to visualize the geographic range of species as the spatial representation of its realized niche. SDMs are also used to forecast range shifts under climate change, but often in the absence of empirical evidence that climate limits population growth. We explored the influence of climate on demography, seasonal migration, and the extent of the geographic range in song sparrows (Melospiza melodia), a species thought to display marked local adaptation to regional climate. To do so, we developed SDMs to predict the demographic and climate niches of migratory and resident song sparrows across our study area in western North America from California to Alaska, using 48 years of demographic data from a focal population in British Columbia and 1.2 million continental-scale citizen science observations. Spatial agreement of our demographic and climate niche models in the region of our focal population was strong (76%), supporting the hypothesis that demographic performance and the occurrence of seasonal migration varied predictably with climatic conditions. In contrast, agreement at the northern (58%) and southern (40%) extents of our study area was lower, as expected if the factors limiting population growth vary regionally. Our results support the hypothesis that local climate drives spatial variation in the occurrence of seasonal migration in song sparrows by limiting the fitness of year-round residents, and suggest that climate warming has favored range expansions and facilitated an upward shift in elevational range song sparrows that forgo seasonal migration. Our work highlights the potential role of seasonal migration in climate adaptation and limits on the reliability of climate niche models not validated with demographic data.

Strong spatial population structure shapes the temporal coevolutionary dynamics of costly female preference and male display.

Female mating preferences for exaggerated male display traits are commonplace. Yet, comprehensive understanding of the evolution and persistence of costly female preference through indirect (Fisherian) selection in finite populations requires some explanation for the persistence of additive genetic variance (Va ) underlying sexual traits, given that directional preference is expected to deplete Va in display and hence halt preference evolution. However, the degree to which Va , and hence preference-display coevolution, may be prolonged by spatially variable sexual selection arising solely from limited gene flow and genetic drift within spatially structured populations has not been examined. Our genetically and spatially explicit model shows that spatial population structure arising in an ecologically homogeneous environment can facilitate evolution and long-term persistence of costly preference given small subpopulations and low dispersal probabilities. Here, genetic drift initially creates spatial variation in female preference, leading to persistence of Va in display through "migration-bias" of genotypes maladapted to emerging local sexual selection, thus fueling coevolution of costly preference and display. However, costs of sexual selection increased the probability of subpopulation extinction, limiting persistence of high preference-display genotypes. Understanding long-term dynamics of sexual selection systems therefore requires joint consideration of coevolution of sexual traits and metapopulation dynamics.

Properties of phenotypic plasticity in discrete threshold traits.

Forms of phenotypic plasticity in key traits, and forms of selection on and genetic variation in such plasticity, fundamentally underpin phenotypic, population dynamic, and evolutionary responses to environmental variation and directional change. Accordingly, numerous theoretical and empirical studies have examined properties and consequences of plasticity, primarily considering traits that are continuously distributed on observed phenotypic scales with linear reaction norms. However, many environmentally sensitive traits are expressed as discrete alternative phenotypes and are appropriately characterized as quantitative genetic threshold traits. Here, we highlight that forms of phenotypic plasticity, genetic variation, and inheritance in plasticity, and outcomes of selection on plasticity, could differ substantially between threshold traits and continuously distributed traits (as are typically considered). We thereby highlight theoretical developments that are required to rationalize and predict phenotypic and microevolutionary dynamics involving plastic threshold traits, and outline how intrinsic properties of such traits could provide relatively straightforward explanations for apparently idiosyncratic observed patterns of phenotypic variation. We summarize how key quantitative genetic parameters underlying threshold traits can be estimated, and thereby set the scene for embedding dynamic discrete traits into theoretical and empirical understanding of the role of plasticity in driving phenotypic, population, and evolutionary responses to environmental variation and change.

Integrating advances in population and evolutionary ecology with conservation strategy through long-term studies of red-billed choughs.

Conceptual and methodological advances in population and evolutionary ecology are often pursued with the ambition that they will help identify demographic, ecological and genetic constraints on population growth rate (λ), and ultimately facilitate evidence-based conservation. However, such advances are often decoupled from conservation practice, impeding translation of scientific understanding into effective conservation and of conservation-motivated research into wider conceptual understanding. We summarise key outcomes from long-term studies of a red-billed chough Pyrrhocorax pyrrhocorax population of conservation concern, where we proactively aimed to achieve the dual and interacting objectives of advancing population and evolutionary ecology and advancing effective conservation. Estimation of means, variances and covariances in key vital rates from individual-based demographic data identified temporal and spatial variation in subadult survival as key constraints on λ, and simultaneously provided new insights into how vital rates can vary as functions of demographic structure, natal conditions and parental life history. Targeted analyses showed that first-year survival increased with prey abundance, implying that food limitation may constrain λ. First-year survival then decreased dramatically, threatening population viability and prompting emergency supplementary feeding interventions. Detailed evaluations suggested that the interventions successfully increased first-year survival in some years and additionally increased adult survival and successful reproduction, thereby feeding back to inform intervention refinements and understanding of complex ecological constraints on λ. Genetic analyses revealed novel evidence of expression of a lethal recessive allele, and demonstrated how critically small effective population size can arise, thereby increasing inbreeding and loss of genetic variation. Population viability analyses parameterised with all available demographic and genetic data showed how ecological and genetic constraints can interact to limit population viability, and identified ecological management as of primacy over genetic management to ensure short-term persistence of the focal population. This case study demonstrates a full iteration through the sequence of primary science, evidence-based intervention, quantitative evaluation and feedback that is advocated in conservation science but still infrequently achieved. It thereby illustrates how pure science advances informed conservation actions to ensure the (short-term) stability of the target population, and how conservation-motivated analyses fed back to advance fundamental understanding of population processes.

Are immigrants outbred and unrelated? Testing standard assumptions in a wild metapopulation.

Immigration into small recipient populations is expected to alleviate inbreeding and increase genetic variation, and hence facilitate population persistence through genetic and/or evolutionary rescue. Such expectations depend on three standard assumptions: that immigrants are outbred, unrelated to existing natives at arrival, and unrelated to each other. These assumptions are rarely explicitly verified, including in key field systems in evolutionary ecology. Yet, they could be violated due to non-random or repeated immigration from adjacent small populations. We combined molecular genetic marker data for 150-160 microsatellite loci with comprehensive pedigree data to test the three assumptions for a song sparrow (Melospiza melodia) population that is a model system for quantifying effects of inbreeding and immigration in the wild. Immigrants were less homozygous than existing natives on average, with mean homozygosity that closely resembled outbred natives. Immigrants can therefore be considered outbred on the focal population scale. Comparisons of homozygosity of real or hypothetical offspring of immigrant-native, native-native and immigrant-immigrant pairings implied that immigrants were typically unrelated to existing natives and to each other. Indeed, immigrants' offspring would be even less homozygous than outbred individuals on the focal population scale. The three standard assumptions of population genetic and evolutionary theory were consequently largely validated. Yet, our analyses revealed some deviations that should be accounted for in future analyses of heterosis and inbreeding depression, implying that the three assumptions should be verified in other systems to probe patterns of non-random or repeated dispersal and facilitate precise and unbiased estimation of key evolutionary parameters.

Episodes of opposing survival and reproductive selection cause strong fluctuating selection on seasonal migration versus residence.

Quantifying temporal variation in sex-specific selection on key ecologically relevant traits, and quantifying how such variation arises through synergistic or opposing components of survival and reproductive selection, is central to understanding eco-evolutionary dynamics, but rarely achieved. Seasonal migration versus residence is one key trait that directly shapes spatio-seasonal population dynamics in spatially and temporally varying environments, but temporal dynamics of sex-specific selection have not been fully quantified. We fitted multi-event capture-recapture models to year-round ring resightings and breeding success data from partially migratory European shags (Phalacrocorax aristotelis) to quantify temporal variation in annual sex-specific selection on seasonal migration versus residence arising through adult survival, reproduction and the combination of both (i.e. annual fitness). We demonstrate episodes of strong and strongly fluctuating selection through annual fitness that were broadly synchronized across females and males. These overall fluctuations arose because strong reproductive selection against migration in several years contrasted with strong survival selection against residence in years with extreme climatic events. These results indicate how substantial phenotypic and genetic variation in migration versus residence could be maintained, and highlight that biologically important fluctuations in selection may not be detected unless both survival selection and reproductive selection are appropriately quantified and combined.

Immigration counter-acts local micro-evolution of a major fitness component: Migration-selection balance in free-living song sparrows.

Ongoing adaptive evolution, and resulting "evolutionary rescue" of declining populations, requires additive genetic variation in fitness. Such variation can be increased by gene flow resulting from immigration, potentially facilitating evolution. But, gene flow could in fact constrain rather than facilitate local adaptive evolution if immigrants have low additive genetic values for local fitness. Local migration-selection balance and micro-evolutionary stasis could then result. However, key quantitative genetic effects of natural immigration, comprising the degrees to which gene flow increases the total local additive genetic variance yet counteracts local adaptive evolutionary change, have not been explicitly quantified in wild populations. Key implications of gene flow for population and evolutionary dynamics consequently remain unclear. Our quantitative genetic analyses of long-term data from free-living song sparrows (Melospiza melodia) show that mean breeding value for local juvenile survival to adulthood, a major component of fitness, increased across cohorts more than expected solely due to drift. Such micro-evolutionary change should be expected given nonzero additive genetic variance and consistent directional selection. However, this evolutionary increase was counteracted by negative additive genetic effects of recent immigrants, which increased total additive genetic variance but prevented a net directional evolutionary increase in total additive genetic value. These analyses imply an approximate quantitative genetic migration-selection balance in a major fitness component, and hence demonstrate a key mechanism by which substantial additive genetic variation can be maintained yet decoupled from local adaptive evolutionary change.