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.

Increased parental effort fails to buffer the cascading effects of warmer seas on common guillemot demographic rates.

Climate warming can reduce food resources for animal populations. In species exhibiting parental care, parental effort is a 'barometer' of changes in environmental conditions. A key issue is the extent to which variation in parental effort can buffer demographic rates against environmental change. Seabirds breed in large, dense colonies and globally are major predators of small fish that are often sensitive to ocean warming. We explored the causes and consequences of annual variation in parental effort as indicated by standardised checks of the proportions of chicks attended by both, one or neither parent, in a population of common guillemots Uria aalge over four decades during which there was marked variation in marine climate and chick diet. We predicted that, for parental effort to be an effective buffer, there would be a link between environmental conditions and parental effort, but not between parental effort and demographic rates. Environmental conditions influenced multiple aspects of the prey delivered by parents to their chicks with prey species, length and energy density all influenced by spring sea surface temperature (sSST) in the current and/or previous year. Overall, the mean annual daily energy intake of chicks declined significantly when sSST in the current year was higher. In accordance with our first prediction, we found that parental effort increased with sSST in the current and previous year. However, the increase was insufficient to maintain chick daily energy intake. In contrast to our second prediction, we found that increased parental effort had major demographic consequences such that growth rate and fledging success of chicks, and body mass and overwinter survival of breeding adults all decreased significantly. Common guillemot parents were unable to compensate effectively for temperature-mediated variation in feeding conditions through behavioural flexibility, resulting in immediate consequences for breeding population size because of lower adult survival and potentially longer-term impacts on recruitment because of lower productivity. These findings highlight that a critical issue for species' responses to future climate change will be the extent to which behavioural buffering can offer resilience to deteriorating environmental conditions.

Site-dependent regulation of breeding success: Evidence for the buffer effect in the common guillemot, a colonially breeding seabird.

Density-dependent regulation can offer resilience to wild populations experiencing fluctuations in environmental conditions because, at lower population sizes, the average quality of habitats or resources is predicted to increase. Site-dependent regulation is a mechanism whereby individuals breed at the highest quality, most successful, sites, leaving poorer quality, less successful sites vacant. As population size increases, higher quality sites become limiting but when populations decline, lower quality sites are vacated first, offering resilience. This process is known as the 'buffer effect'. However, few studies have tested whether such regulation operates in populations experiencing changes in size and trend. We used data from a population of common guillemots Uria aalge, a colonially breeding seabird, to investigate the relationship between site occupancy probability, site quality and population size and trend. These data were collected at five sub-colonies spanning a 38-year period (1981-2018) comprising phases of population increase, decrease and recovery. We first tested whether site quality and population size in sub-colonies explained which sites were occupied for breeding, and if this was robust to changes in sub-colony trend. We then investigated whether disproportionate use of higher quality sites drove average site quality and breeding success across sub-colony sizes and trends. Finally, we tested whether individuals consistently occupied higher quality sites during periods of decline and recovery. Higher quality sites were disproportionality used when sub-colony size was smaller, resulting in higher average site quality and breeding success at lower population sizes. This relationship was unaffected by changes in sub-colony trend. However, contrary to the predictions of the buffer effect, new sites were established at a similar rate to historically occupied sites during sub-colony decline and recovery despite being of lower quality. Our results provide support for the buffer effect conferring resilience to populations, such that average breeding success was consistently higher at lower population size during all phases of population change. However, this process was tempered by the continued establishment of new, lower quality, sites which could act to slow population recovery after periods when colony size was low.

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.

Variation and correlation in the timing of breeding of North Atlantic seabirds across multiple scales.

Timing of breeding, an important driver of fitness in many populations, is widely studied in the context of global change, yet despite considerable efforts to identify environmental drivers of seabird nesting phenology, for most populations we lack evidence of strong drivers. Here we adopt an alternative approach, examining the degree to which different populations positively covary in their annual phenology to infer whether phenological responses to environmental drivers are likely to be (a) shared across species at a range of spatial scales, (b) shared across populations of a species or (c) idiosyncratic to populations. We combined 51 long-term datasets on breeding phenology spanning 50 years from nine seabird species across 29 North Atlantic sites and examined the extent to which different populations share early versus late breeding seasons depending on a hierarchy of spatial scales comprising breeding site, small-scale region, large-scale region and the whole North Atlantic. In about a third of cases, we found laying dates of populations of different species sharing the same breeding site or small-scale breeding region were positively correlated, which is consistent with the hypothesis that they share phenological responses to the same environmental conditions. In comparison, we found no evidence for positive phenological covariation among populations across species aggregated at larger spatial scales. In general, we found little evidence for positive phenological covariation between populations of a single species, and in many instances the inter-year variation specific to a population was substantial, consistent with each population responding idiosyncratically to local environmental conditions. Black-legged kittiwake Rissa tridactyla was the exception, with populations exhibiting positive covariation in laying dates that decayed with the distance between breeding sites, suggesting that populations may be responding to a similar driver. Our approach sheds light on the potential factors that may drive phenology in our study species, thus furthering our understanding of the scales at which different seabirds interact with interannual variation in their environment. We also identify additional systems and phenological questions to which our inferential approach could be applied.

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.

Meeting Paris agreement objectives will temper seabird winter distribution shifts in the North Atlantic Ocean.

We explored the implications of reaching the Paris Agreement Objective of limiting global warming to <2°C for the future winter distribution of the North Atlantic seabird community. We predicted and quantified current and future winter habitats of five North Atlantic Ocean seabird species (Alle alle, Fratercula arctica, Uria aalge, Uria lomvia and Rissa tridactyla) using tracking data for ~1500 individuals through resource selection functions based on mechanistic modeling of seabird energy requirements, and a dynamic bioclimate envelope model of seabird prey. Future winter distributions were predicted to shift with climate change, especially when global warming exceed 2°C under a "no mitigation" scenario, modifying seabird wintering hotspots in the North Atlantic Ocean. Our findings suggest that meeting Paris agreement objectives will limit changes in seabird selected habitat location and size in the North Atlantic Ocean during the 21st century. We thereby provide key information for the design of adaptive marine-protected areas in a changing ocean.

Strong survival selection on seasonal migration versus residence induced by extreme climatic events.

Elucidating the full eco-evolutionary consequences of climate change requires quantifying the impact of extreme climatic events (ECEs) on selective landscapes of key phenotypic traits that mediate responses to changing environments. Episodes of strong ECE-induced selection could directly alter population composition, and potentially drive micro-evolution. However, to date, few studies have quantified ECE-induced selection on key traits, meaning that immediate and longer-term eco-evolutionary implications cannot yet be considered. One widely expressed trait that allows individuals to respond to changing seasonal environments, and directly shapes spatio-seasonal population dynamics, is seasonal migration versus residence. Many populations show considerable among-individual phenotypic variation, resulting in 'partial migration'. However, variation in the magnitude of direct survival selection on migration versus residence has not been rigorously quantified, and empirical evidence of whether seasonal ECEs induce, intensify, weaken or reverse such selection is lacking. We designed full annual cycle multi-state capture-recapture models that allow estimation of seasonal survival probabilities of migrants and residents from spatio-temporally heterogeneous individual resightings. We fitted these models to 9 years of geographically extensive year-round resighting data from partially migratory European shags Phalacrocorax aristotelis. We thereby quantified seasonal and annual survival selection on migration versus residence across benign and historically extreme non-breeding season (winter) conditions, and tested whether selection differed between females and males. We show that two of four observed ECEs, defined as severe winter storms causing overall low survival, were associated with very strong seasonal survival selection against residence. These episodes dwarfed the weak selection or neutrality evident otherwise, and hence caused selection through overall annual survival. The ECE that caused highest overall mortality and strongest selection also caused sex-biased mortality, but there was little overall evidence of sex-biased selection on migration versus residence. Our results imply that seasonal ECEs and associated mortality can substantially shape the landscape of survival selection on migration versus residence. Such ECE-induced phenotypic selection will directly alter migrant and resident frequencies, and thereby alter immediate spatio-seasonal population dynamics. Given underlying additive genetic variation, such ECEs could potentially cause micro-evolutionary changes in seasonal migration, and thereby cause complex eco-evolutionary population responses to changing seasonal environments.

No evidence for fitness signatures consistent with increasing trophic mismatch over 30 years in a population of European shag Phalacrocorax aristotelis.

As temperatures rise, timing of reproduction is changing at different rates across trophic levels, potentially resulting in asynchrony between consumers and their resources. The match-mismatch hypothesis (MMH) suggests that trophic asynchrony will have negative impacts on average productivity of consumers. It is also thought to lead to selection on timing of breeding, as the most asynchronous individuals will show the greatest reductions in fitness. Using a 30-year individual-level dataset of breeding phenology and success from a population of European shags on the Isle of May, Scotland, we tested a series of predictions consistent with the hypothesis that fitness impacts of trophic asynchrony are increasing. These predictions quantified changes in average annual breeding success and strength of selection on timing of breeding, over time and in relation to rising sea surface temperature (SST) and diet composition. Annual average (population) breeding success was negatively correlated with average lay date yet showed no trend over time, or in relation to increasing SST or the proportion of principal prey in the diet, as would be expected if trophic mismatch was increasing. At the individual level, we found evidence for stabilising selection and directional selection for earlier breeding, although the earliest birds were not the most productive. However, selection for earlier laying did not strengthen over time, or in relation to SST or slope of the seasonal shift in diet from principal to secondary prey. We found that the optimum lay date advanced by almost 4 weeks during the study, and that the population mean lay date tracked this shift. Our results indicate that average performance correlates with absolute timing of breeding of the population, and there is selection for earlier laying at the individual level. However, we found no fitness signatures of a change in the impact of climate-induced trophic mismatch, and evidence that shags are tracking long-term shifts in optimum timing. This suggests that if asynchrony is present in this system, breeding success is not impacted. Our approach highlights the advantages of examining variation at both population and individual levels when assessing evidence for fitness impacts of trophic asynchrony.

Interactions between Environmental Contaminants and Gastrointestinal Parasites: Novel Insights from an Integrative Approach in a Marine Predator.

Environmental contaminants and parasites are ubiquitous stressors that can affect animal physiology and derive from similar dietary sources (co-exposure). To unravel their interactions in wildlife, it is thus essential to quantify their concurring drivers. Here, the relationship between blood contaminant residues (11 trace elements and 17 perfluoroalkyl substances) and nonlethally quantified gastrointestinal parasite loads was tested while accounting for intrinsic (sex, age, and mass) and extrinsic factors (trophic ecology inferred from stable isotope analyses and biologging) in European shags Phalacrocorax aristotelis. Shags had high mercury (range 0.65-3.21 µg g-1 wet weight, ww) and extremely high perfluorooctanoic acid (PFOA) and perfluorononanoic acid (PFNA) residues (3.46-53 and 4.48-44 ng g-1 ww, respectively). Males had higher concentrations of arsenic, mercury, PFOA, and PFNA than females, while the opposite was true for selenium, perfluorododecanoic acid (PFDoA), and perfluooctane sulfonic acid (PFOS). Individual parasite loads (Contracaecum rudolphii) were higher in males than in females. Females targeted pelagic-feeding prey, while males relied on both pelagic- and benthic-feeding organisms. Parasite loads were not related to trophic ecology in either sex, suggesting no substantial dietary co-exposure with contaminants. In females, parasite loads increased strongly with decreasing selenium:mercury molar ratios. Females may be more susceptible to the interactive effects of contaminants and parasites on physiology, with potential fitness consequences.

The energetic cost of parasitism in a wild population.

Parasites have profound fitness effects on their hosts, yet these are often sub-lethal, making them difficult to understand and quantify. A principal sub-lethal mechanism that reduces fitness is parasite-induced increase in energetic costs of specific behaviours, potentially resulting in changes to time and energy budgets. However, quantifying the influence of parasites on these costs has not been undertaken in free-living animals. We used accelerometers to estimate energy expenditure on flying, diving and resting, in relation to a natural gradient of endo-parasite loads in a wild population of European shags Phalacrocorax aristotelis We found that flight costs were 10% higher in adult females with higher parasite loads and these individuals spent 44% less time flying than females with lower parasite loads. There was no evidence for an effect of parasite load on daily energy expenditure, suggesting the existence of an energy ceiling, with the increase in cost of flight compensated for by a reduction in flight duration. These behaviour specific costs of parasitism will have knock-on effects on reproductive success, if constraints on foraging behaviour detrimentally affect provisioning of young. The findings emphasize the importance of natural parasite loads in shaping the ecology and life-history of their hosts, which can have significant population level consequences.

The role of parasitism in the energy management of a free-ranging bird.

Parasites often prompt sub-lethal costs to their hosts by eliciting immune responses. These costs can be hard to quantify but are crucial to our understanding of the host's ecology. Energy is a fundamental currency to quantify these costs, as energetic trade-offs often exist between key fitness-related processes. Daily energy expenditure (DEE) comprises of resting metabolic rate (RMR) and energy available for activity, which are linked via the energy management strategy of an organism. Parasitism may play a role in the balance between self-maintenance and activity, as immune costs can be expressed in elevated RMR. Therefore, understanding energy use in the presence of parasitism enables mechanistic elucidation of potential parasite costs. Using a gradient of natural parasite load and proxies for RMR and DEE in a wild population of breeding European shags (Phalacrocorax aristotelis), we tested the effect of parasitism on maintenance costs as well as the relationship between proxies for RMR and DEE. We found a positive relationship between parasite load and our RMR proxy in females but not males, and no relationship between proxies for RMR and DEE. This provides evidence for increased maintenance costs in individuals with higher parasite loads and suggests the use of an allocation energy management strategy, whereby an increase to RMR creates restrictions on energy allocation to other activities. This is likely to have fitness consequences as energy allocated to immunity is traded off against reproduction. Our findings demonstrate that understanding energy management strategies alongside fitness drivers is central to understanding the mechanisms by which these drivers influence individual fitness.

Breeding density, fine-scale tracking, and large-scale modeling reveal the regional distribution of four seabird species.

Population-level estimates of species' distributions can reveal fundamental ecological processes and facilitate conservation. However, these may be difficult to obtain for mobile species, especially colonial central-place foragers (CCPFs; e.g., bats, corvids, social insects), because it is often impractical to determine the provenance of individuals observed beyond breeding sites. Moreover, some CCPFs, especially in the marine realm (e.g., pinnipeds, turtles, and seabirds) are difficult to observe because they range tens to ten thousands of kilometers from their colonies. It is hypothesized that the distribution of CCPFs depends largely on habitat availability and intraspecific competition. Modeling these effects may therefore allow distributions to be estimated from samples of individual spatial usage. Such data can be obtained for an increasing number of species using tracking technology. However, techniques for estimating population-level distributions using the telemetry data are poorly developed. This is of concern because many marine CCPFs, such as seabirds, are threatened by anthropogenic activities. Here, we aim to estimate the distribution at sea of four seabird species, foraging from approximately 5,500 breeding sites in Britain and Ireland. To do so, we GPS-tracked a sample of 230 European Shags Phalacrocorax aristotelis, 464 Black-legged Kittiwakes Rissa tridactyla, 178 Common Murres Uria aalge, and 281 Razorbills Alca torda from 13, 20, 12, and 14 colonies, respectively. Using Poisson point process habitat use models, we show that distribution at sea is dependent on (1) density-dependent competition among sympatric conspecifics (all species) and parapatric conspecifics (Kittiwakes and Murres); (2) habitat accessibility and coastal geometry, such that birds travel further from colonies with limited access to the sea; and (3) regional habitat availability. Using these models, we predict space use by birds from unobserved colonies and thereby map the distribution at sea of each species at both the colony and regional level. Space use by all four species' British breeding populations is concentrated in the coastal waters of Scotland, highlighting the need for robust conservation measures in this area. The techniques we present are applicable to any CCPF.

Reproductive performance of resident and migrant males, females and pairs in a partially migratory bird.

Quantifying among-individual variation in life-history strategies, and associated variation in reproductive performance and resulting demographic structure, is key to understanding and predicting population dynamics and life-history evolution. Partial migration, where populations comprise a mixture of resident and seasonally migrant individuals, constitutes a dimension of life-history variation that could be associated with substantial variation in reproductive performance. However, such variation has rarely been quantified due to the challenge of measuring reproduction and migration across a sufficient number of seasonally mobile males and females. We used intensive winter (non-breeding season) resightings of colour-ringed adult European shags (Phalacrocorax aristotelis) from a known breeding colony to identify resident and migrant individuals. We tested whether two aspects of annual reproductive performance, brood hatch date and breeding success, differed between resident and migrant males, females and breeding pairs observed across three consecutive winters and breeding seasons. The sex ratios of observed resident and migrant shags did not significantly differ from each other or from 1:1, suggesting that both sexes are partially migratory and that migration was not sex-biased across surveyed areas. Individual resident males and females hatched their broods 6 days earlier and fledged 0.2 more chicks per year than migrant males and females on average. Resident individuals of both sexes therefore had higher breeding success than migrants. Hatch date and breeding success also varied with a pair's joint migratory strategy such that resident-resident pairs hatched their broods 12 days earlier than migrant-migrant pairs, and fledged 0.7 more chicks per year on average. However, there was no evidence of assortative pairing with respect to migratory strategy: observed frequencies of migrant-migrant and resident-resident pairs did not differ from those expected given random pairing. These data demonstrate substantial variation in two key aspects of reproductive performance associated with the migratory strategies of males, females and breeding pairs within a partially migratory population. These patterns could reflect direct and/or indirect mechanisms, but imply that individual variation in migratory strategy and variation in pairing among residents and migrants could influence selection on migration and drive complex population and evolutionary dynamics.

Protracted treatment with corticosterone reduces breeding success in a long-lived bird.

Determining the physiological mechanisms underpinning life-history decisions is essential for understanding the constraints under which life-history strategies can evolve. In long-lived species, where the residual reproductive value of breeders is high, adult survival is a key contributor to lifetime reproductive success. We therefore expect that when adult survival is compromised during reproduction, mechanisms will evolve to redirect resources away from reproduction, with implications for reproductive hormones, adult body mass, nest attendance behaviour and breeding success. We investigated whether manipulating corticosterone, to simulate exposure to an environmental stressor, affected the secretion of prolactin and breeding success in the black-legged kittiwake Rissa tridactyla. We used implanted Alzet® osmotic pumps to administer corticosterone to incubating kittiwakes at a constant rate over a period of approximately 8days. Manipulated birds were compared with sham implanted birds and control birds, which had no implants. There was no significant difference in the body mass of captured individuals at the time of implantation and implant removal. Corticosterone-implanted males showed lower nest attendance during the chick rearing period compared to sham-implanted males; the opposite pattern was found in females. Corticosterone treated birds showed a marginally significant reduction in breeding success compared to sham-implanted individuals, with all failures occurring at least 1week after implant removal. However, prolactin concentrations at implant removal were not significantly different from initial values. We were unable to measure the profile of change in corticosterone during the experiment. However, our results suggest a delayed effect of elevated corticosterone on breeding success rather than an immediate suppression of prolactin concentrations causing premature failure.

Genetic and biochemical differences in populations bred for extremes in maize grain methionine concentration.

BACKGROUND: Methionine is an important nutrient in animal feed and several approaches have been developed to increase methionine concentration in maize (Zea mays L.) grain. One approach is through traditional breeding using recurrent selection. Using divergent selection, genetically related populations with extreme differences in grain methionine content were produced. In order to better understand the molecular mechanisms controlling grain methionine content, we examined seed proteins, transcript levels of candidate genes, and genotypes of these populations. RESULTS: Two populations were selected for high or low methionine concentration for eight generations and 40 and 56% differences between the high and low populations in grain methionine concentration were observed. Mean values between the high and low methionine populations differed by greater than 1.5 standard deviations in some cycles of selection. Other amino acids and total protein concentration exhibited much smaller changes. In an effort to understand the molecular mechanisms that contribute to these differences, we compared transcript levels of candidate genes encoding high methionine seed storage proteins involved in sulfur assimilation or methionine biosynthesis. In combination, we also explored the genetic mechanisms at the SNP level through implementation of an association analysis. Significant differences in methionine-rich seed storage protein genes were observed in comparisons of high and low methionine populations, while transcripts of seed storage proteins lacking high levels of methionine were unchanged. Seed storage protein levels were consistent with transcript levels. Two genes involved in sulfur assimilation, Cys2 and CgS1 showed substantial differences in allele frequencies when two selected populations were compared to the starting populations. Major genes identified across cycles of selection by a high-stringency association analysis included dzs18, wx, dzs10, and zp27. CONCLUSIONS: We hypothesize that transcriptional changes alter sink strength by altering the levels of methionine-rich seed storage proteins. To meet the altered need for sulfur, a cysteine-rich seed storage protein is altered while sulfur assimilation and methionine biosynthesis throughput is changed by selection for certain alleles of Cys2 and CgS1.

Seabirds show foraging site and route fidelity but demonstrate flexibility in response to local information.

BACKGROUND: Fidelity to a given foraging location or route may be beneficial when environmental conditions are predictable but costly if conditions deteriorate or become unpredictable. Understanding the magnitude of fidelity displayed by different species and the processes that drive or erode it is therefore vital for understanding how fidelity may shape the demographic consequences of anthropogenic change. In particular, understanding the information that individuals may use to adjust their fidelity will facilitate improved predictions of how fidelity may change as environments change and the extent to which it will buffer individuals against such changes. METHODS: We used movement data collected during the breeding season across eight years for common guillemots, Atlantic puffins, razorbills, and black-legged kittiwakes breeding on the Isle of May, Scotland to understand: (1) whether foraging site/route fidelity occurred within and between years, (2) whether the degree of fidelity between trips was predicted by personal foraging effort, and (3) whether different individuals made more similar trips when they overlapped in time at the colony prior to departure and/or when out at sea suggesting the use of the same local environmental cues or information on the decisions made by con- and heterospecifics. RESULTS: All species exhibited site and route fidelity both within- and between-years, and fidelity between trips in guillemots and razorbills was related to metrics of foraging effort, suggesting they adjust fidelity to their personal foraging experience. We also found evidence that individuals used local environmental cues of prey location or availability and/or information gained by observing conspecifics when choosing foraging routes, particularly in puffins, where trips of individuals that overlapped temporally at the colony or out at sea were more similar. CONCLUSIONS: The fidelity shown by these seabird species has the potential to put them at greater risk in the face of environmental change by driving individuals to continue using areas being degraded by anthropogenic pressures. However, our results suggest that individuals show some flexibility in their fidelity, which may promote resilience under environmental change. The benefits of this flexibility are likely to depend on numerous factors, including the rapidity and spatial scale of environmental change and the reliability of the information individuals use to choose foraging sites or routes, thus highlighting the need to better understand how organisms combine cues, prior experience, and other sources of information to make movement decisions.

Quantifying the impacts of predation by Great Black-backed Gulls Larus marinus on an Atlantic Puffin Fratercula arctica population: Implications for conservation management and impact assessments.

The management of predator-prey conflicts can be a key aspect of species conservation. For management approaches to be effective, a robust understanding of the predator-prey relationship is needed, particularly when both predator and prey are species of conservation concern. On the Isle of May, Firth of Forth, Scotland, numbers of breeding Great Black-backed Gulls Larus marinus, a generalist predator, have been increasing since the 1980s, which has led to increasing numbers of sympatrically breeding Atlantic Puffins Fratercula arctica being predated during the breeding season. This may have consequences for species management on the Isle of May and impact assessments of offshore windfarms in the wider Firth of Forth area. We used population viability analysis to quantify under what predation pressure the Atlantic Puffin population may decline and become locally extinct over a three-generation period. The predation level empirically estimated in 2017 (1120 Puffins per year) was not sufficient to drive a decline in the Puffin population. Rather, an increase to approximately 3000 Puffins per year would be required to cause a population decline, and >4000 to drive the population to quasi-extinction within 66 years. We discuss the likelihood of such a scenario being reached on the Isle of May, and we recommend that where predator-prey conflicts occur, predation-driven mortality should be regularly quantified to inform conservation management and population viability analyses associated with impact assessments.

Consequences of cross-season demographic correlations for population viability.

Demographic correlations are pervasive in wildlife populations and can represent important secondary drivers of population growth. Empirical evidence suggests that correlations are in general positive for long-lived species, but little is known about the degree of variation among spatially segregated populations of the same species in relation to environmental conditions. We assessed the relative importance of two cross-season correlations in survival and productivity, for three Atlantic puffin (Fratercula arctica) populations with contrasting population trajectories and non-overlapping year-round distributions. The two correlations reflected either a relationship between adult survival prior to breeding on productivity, or a relationship between productivity and adult survival the subsequent year. Demographic rates and their correlations were estimated with an integrated population model, and their respective contributions to variation in population growth were calculated using a transient-life table response experiment. For all three populations, demographic correlations were positive at both time lags, although their strength differed. Given the different year-round distributions of these populations, this variation in the strength population-level demographic correlations points to environmental conditions as an important driver of demographic variation through life-history constraints. Consequently, the contributions of variances and correlations in demographic rates to population growth rates differed among puffin populations, which has implications for-particularly small-populations' viability under environmental change as positive correlations tend to reduce the stochastic population growth rate.

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.