Genotypes selected for early and late avian lay date differ in their phenotype, but not fitness, in the wild.

Global warming has shifted phenological traits in many species, but whether species are able to track further increasing temperatures depends on the fitness consequences of additional shifts in phenological traits. To test this, we measured phenology and fitness of great tits (Parus major) with genotypes for extremely early and late egg lay dates, obtained from a genomic selection experiment. Females with early genotypes advanced lay dates relative to females with late genotypes, but not relative to nonselected females. Females with early and late genotypes did not differ in the number of fledglings produced, in line with the weak effect of lay date on the number of fledglings produced by nonselected females in the years of the experiment. Our study is the first application of genomic selection in the wild and led to an asymmetric phenotypic response that indicates the presence of constraints toward early, but not late, lay dates.

A partial migrant relies upon a range-wide cue set but uses population-specific weighting for migratory timing.

BACKGROUND: Many birds species range over vast geographic regions and migrate seasonally between their breeding and overwintering sites. Deciding when to depart for migration is one of the most consequential life-history decisions an individual may make. However, it is still not fully understood which environmental cues are used to time the onset of migration and to what extent their relative importance differs across a range of migratory strategies. We focus on departure decisions of a songbird, the Eurasian blackbird Turdus merula, in which selected Russian and Polish populations are full migrants which travel relatively long-distances, whereas Finnish and German populations exhibit partial migration with shorter migration distances. METHODS: We used telemetry data from the four populations (610 individuals) to determine which environmental cues individuals from each population use to initiate their autumn migration. RESULTS: When departing, individuals in all populations selected nights with high atmospheric pressure and minimal cloud cover. Fully migratory populations departed earlier in autumn, at longer day length, at higher ambient temperatures, and during nights with higher relative atmospheric pressure and more supportive winds than partial migrants; however, they did not depart in higher synchrony. Thus, while all studied populations used the same environmental cues, they used population-specific and locally tuned thresholds to determine the day of departure. CONCLUSIONS: Our data support the idea that migratory timing is controlled by general, species-wide mechanisms, but fine-tuned thresholds in response to local conditions.

Recent natural variability in global warming weakened phenological mismatch and selection on seasonal timing in great tits (Parus major).

Climate change has led to phenological shifts in many species, but with large variation in magnitude among species and trophic levels. The poster child example of the resulting phenological mismatches between the phenology of predators and their prey is the great tit (Parus major), where this mismatch led to directional selection for earlier seasonal breeding. Natural climate variability can obscure the impacts of climate change over certain periods, weakening phenological mismatching and selection. Here, we show that selection on seasonal timing indeed weakened significantly over the past two decades as increases in late spring temperatures have slowed down. Consequently, there has been no further advancement in the date of peak caterpillar food abundance, while great tit phenology has continued to advance, thereby weakening the phenological mismatch. We thus show that the relationships between temperature, phenologies of prey and predator, and selection on predator phenology are robust, also in times of a slowdown of warming. Using projected temperatures from a large ensemble of climate simulations that take natural climate variability into account, we show that prey phenology is again projected to advance faster than great tit phenology in the coming decades, and therefore that long-term global warming will intensify phenological mismatches.

Species-specific effects of thermal stress on the expression of genetic variation across a diverse group of plant and animal taxa under experimental conditions.

Assessing the genetic adaptive potential of populations and species is essential for better understanding evolutionary processes. However, the expression of genetic variation may depend on environmental conditions, which may speed up or slow down evolutionary responses. Thus, the same selection pressure may lead to different responses. Against this background, we here investigate the effects of thermal stress on genetic variation, mainly under controlled laboratory conditions. We estimated additive genetic variance (VA), narrow-sense heritability (h2) and the coefficient of genetic variation (CVA) under both benign control and stressful thermal conditions. We included six species spanning a diverse range of plant and animal taxa, and a total of 25 morphological and life-history traits. Our results show that (1) thermal stress reduced fitness components, (2) the majority of traits showed significant genetic variation and that (3) thermal stress affected the expression of genetic variation (VA, h2 or CVA) in only one-third of the cases (25 of 75 analyses, mostly in one clonal species). Moreover, the effects were highly species-specific, with genetic variation increasing in 11 and decreasing in 14 cases under stress. Our results hence indicate that thermal stress does not generally affect the expression of genetic variation under laboratory conditions but, nevertheless, increases or decreases genetic variation in specific cases. Consequently, predicting the rate of genetic adaptation might not be generally complicated by environmental variation, but requires a careful case-by-case consideration.

Temperature has a causal and plastic effect on timing of breeding in a small songbird.

Phenotypic plasticity is an important mechanism by which an individual can adapt its seasonal timing to predictable, short-term environmental changes by using predictive cues. Identification of these cues is crucial to forecast the response of species to long-term environmental change and to study their potential to adapt. Individual great tits (Parus major) start reproduction early under warmer conditions in the wild, but whether this effect is causal is not well known. We housed 36 pairs of great tits in climate-controlled aviaries and 40 pairs in outdoor aviaries, where they bred under artificial contrasting temperature treatments or in semi-natural conditions, respectively, for two consecutive years, using birds from lines selected for early and late egg laying. We thus obtained laying dates in two different thermal environments for each female. Females bred earlier under warmer conditions in climate-controlled aviaries, but not in outdoor aviaries. The latter was inconsistent with laying dates from our wild population. Further, early selection line females initiated egg laying consistently ∼9 days earlier than late selection line females in outdoor aviaries, but we found no difference in the degree of plasticity (i.e. the sensitivity to temperature) in laying date between selection lines. Because we found that temperature causally affects laying date, climate change will lead to earlier laying. This advancement is, however, unlikely to be sufficient, thereby leading to selection for earlier laying. Our results suggest that natural selection may lead to a change in mean phenotype, but not to a change in the sensitivity of laying dates to temperature.

Opinion: Is gene mapping in wild populations useful for understanding and predicting adaptation to global change?

Changing environmental conditions will inevitably alter selection pressures. Over the long term, populations have to adapt to these altered conditions by evolutionary change to avoid extinction. Quantifying the 'evolutionary potential' of populations to predict whether they will be able to adapt fast enough to forecasted changes is crucial to fully assess the threat for biodiversity posed by climate change. Technological advances in sequencing and high-throughput genotyping have now made genomic studies possible in a wide range of species. Such studies, in theory, allow an unprecedented understanding of the genomics of ecologically relevant traits and thereby a detailed assessment of the population's evolutionary potential. Aimed at a wider audience than only evolutionary geneticists, this paper gives an overview of how gene-mapping studies have contributed to our understanding and prediction of evolutionary adaptations to climate change, identifies potential reasons why their contribution to understanding adaptation to climate change may remain limited, and highlights approaches to study and predict climate change adaptation that may be more promising, at least in the medium term.

Quantifying individual variation in reaction norms: Mind the residual.

Phenotypic plasticity is a central topic in ecology and evolution. Individuals may differ in the degree of plasticity (individual-by-environment interaction (I × E)), which has implications for the capacity of populations to respond to selection. Random regression models (RRMs) are a popular tool to study I × E in behavioural or life-history traits, yet evidence for I × E is mixed, differing between species, populations, and even between studies on the same population. One important source of discrepancies between studies is the treatment of heterogeneity in residual variance (heteroscedasticity). To date, there seems to be no collective awareness among ecologists of its influence on the estimation of I × E or a consensus on how to best model it. We performed RRMs with differing residual variance structures on simulated data with varying degrees of heteroscedasticity and plasticity, sample size and environmental variability to test how RRMs would perform under each scenario. The residual structure in the RRMs affected the precision of estimates of simulated I × E as well as statistical power, with substantial lack of precision and high false-positive rates when sample size, environmental variability and plasticity were small. We show that model comparison using information criteria can be used to choose among residual structures and reinforce this point by analysis of real data of two study populations of great tits (Parus major). We provide guidelines that can be used by biologists studying I × E that, ultimately, should lead to a reduction in bias in the literature concerning the statistical evidence and the reported magnitude of variation in plasticity.

Comparing two measures of phenological synchrony in a predator-prey interaction: Simpler works better.

Global climate change has sparked a vast research effort into the demographic and evolutionary consequences of mismatches between consumer and resource phenology. Many studies have used the difference in peak dates to quantify phenological synchrony (match in dates, MD), but this approach has been suggested to be inconclusive, since it does not incorporate the temporal overlap between the phenological distributions (match in overlap, MO). We used 24 years of detailed data on the phenology of a predator-prey system, the great tit (Parus major) and the main food for its nestlings, caterpillars, to estimate MD and MO at the population and brood levels. We compared the performance of both metrics on two key demographic parameters: offspring recruitment probability and selection on the timing of reproduction. Although MD and MO correlated quadratically as expected, MD was a better predictor for both offspring recruitment and selection on timing than MO. We argue-and verify through simulations-that this is because quantifying MO has to be based on nontrivial, difficult-to-verify assumptions that likely render MO too inaccurate as a proxy for food availability in practice. Our results have important implications for the allocation of research efforts in long-term population studies in highly seasonal environments.

Exploration of tissue-specific gene expression patterns underlying timing of breeding in contrasting temperature environments in a song bird.

BACKGROUND: Seasonal timing of breeding is a life history trait with major fitness consequences but the genetic basis of the physiological mechanism underlying it, and how gene expression is affected by date and temperature, is not well known. In order to study this, we measured patterns of gene expression over different time points in three different tissues of the hypothalamic-pituitary-gonadal-liver axis, and investigated specifically how temperature affects this axis during breeding. We studied female great tits (Parus major) from lines artificially selected for early and late timing of breeding that were housed in two contrasting temperature environments in climate-controlled aviaries. We collected hypothalamus, liver and ovary samples at three different time points (before and after onset of egg-laying). For each tissue, we sequenced whole transcriptomes of 12 pools (n = 3 females) to analyse gene expression. RESULTS: Birds from the selection lines differed in expression especially for one gene with clear reproductive functions, zona pellucida glycoprotein 4 (ZP4), which has also been shown to be under selection in these lines. Genes were differentially expressed at different time points in all tissues and most of the differentially expressed genes between the two temperature treatments were found in the liver. We identified a set of hub genes from all the tissues which showed high association to hormonal functions, suggesting that they have a core function in timing of breeding. We also found ample differentially expressed genes with largely unknown functions in birds. CONCLUSIONS: We found differentially expressed genes associated with selection line and temperature treatment. Interestingly, the latter mainly in the liver suggesting that temperature effects on egg-laying date may happen down-stream in the physiological pathway. These findings, as well as our datasets, will further the knowledge of the mechanisms of tissue-specific avian seasonality in the future.

Between- and Within-Individual Variation of Maternal Thyroid Hormone Deposition in Wild Great Tits (Parus major).

Maternal hormones are often considered a mediator of anticipatory maternal effects; namely, mothers adjust maternal hormone transfer to prepare the offspring for the anticipated environment. The flexibility for mothers to adjust hormone transfer is therefore a prerequisite for such anticipatory maternal effects. Nevertheless, previous studies have focused only on the average differences of maternal hormone transfer between groups and neglected the substantial individual variation, despite the fact that individual plasticity in maternal hormone transfer is actually the central assumption. In this study, we studied the between- and within-individual variation of maternal thyroid hormones (THs) in egg yolk of wild great tits (Parus major) and estimated the individual plasticity of maternal yolk THs across environmental temperature, clutch initiation dates, and egg laying order using linear mixed effects models. Interestingly, our models provide statistical evidence that the two main THs-the main biologically active hormone T3 and T4, which is mostly considered a prohormone-exhibited different variation patterns. Yolk T3 showed significant between-individual variation on the average levels, in line with its previously reported moderate heritability. Yolk T4, however, showed significant between-clutch variation in the pattern over the laying sequence, suggesting a great within-individual plasticity. Our findings suggest that the role and function of the hormone within the endocrine axis likely influences its flexibility to respond to environmental change. Whether the flexibility of T4 deposition brings a fitness advantage should be examined along with its potential effects on offspring, which remain to be further investigated.

Fine-tuning of seasonal timing of breeding is regulated downstream in the underlying neuro-endocrine system in a small songbird.

The timing of breeding is under selection in wild populations as a result of climate change, and understanding the underlying physiological processes mediating this timing provides insight into the potential rate of adaptation. Current knowledge on this variation in physiology is, however, mostly limited to males. We assessed whether individual differences in the timing of breeding in females are reflected in differences in candidate gene expression and, if so, whether these differences occur in the upstream (hypothalamus) or downstream (ovary and liver) parts of the neuroendocrine system. We used 72 female great tits from two generations of lines artificially selected for early and late egg laying, which were housed in climate-controlled aviaries and went through two breeding cycles within 1 year. In the first breeding season we obtained individual egg-laying dates, while in the second breeding season, using the same individuals, we sampled several tissues at three time points based on the timing of the first breeding attempt. For each tissue, mRNA expression levels were measured using qPCR for a set of candidate genes associated with the timing of reproduction and subsequently analysed for differences between generations, time points and individual timing of breeding. We found differences in gene expression between generations in all tissues, with the most pronounced differences in the hypothalamus. Differences between time points, and early- and late-laying females, were found exclusively in the ovary and liver. Altogether, we show that fine-tuning of the seasonal timing of breeding, and thereby the opportunity for adaptation in the neuroendocrine system, is regulated mostly downstream in the neuro-endocrine system.

Genomic selection on breeding time in a wild bird population.

Artificial selection experiments are a powerful tool in evolutionary biology. Selecting individuals based on multimarker genotypes (genomic selection) has several advantages over phenotype-based selection but has, so far, seen very limited use outside animal and plant breeding. Genomic selection depends on the markers tagging the causal loci that underlie the selected trait. Because the number of necessary markers depends, among other factors, on effective population size, genomic selection may be in practice not feasible in wild populations as most wild populations have much higher effective population sizes than domesticated populations. However, the current possibilities of cost-effective high-throughput genotyping could overcome this limitation and thereby make it possible to apply genomic selection also in wild populations. Using a unique dataset of about 2000 wild great tits (Parus major), a small passerine bird, genotyped on a 650 k SNP chip we calculated genomic breeding values for egg-laying date using the so-called GBLUP approach. In this approach, the pedigree-based relatedness matrix of an "animal model," a special form of the mixed model, is replaced by a marker-based relatedness matrix. Using the marker-based relatedness matrix, the model seemed better able to disentangle genetic and permanent environmental effects. We calculated the accuracy of genomic breeding values by correlating them to the phenotypes of individuals whose phenotypes were excluded from the analysis when estimating the genomic breeding values. The obtained accuracy was about 0.20, with very little effect of the used genomic relatedness estimator but a strong effect of the number of SNPs. The obtained accuracy is lower than typically seen in domesticated species but considerable for a trait with low heritability (∼0.2) as avian breeding time. Our results show that genomic selection is possible also in wild populations with potentially many applications, which we discuss here.

Short-term, but not long-term, increased daytime workload leads to decreased night-time energetics in a free-living song bird.

Reproduction is energetically expensive and to obtain sufficient energy, animals can either alter their metabolic system over time to increase energy intake (increased-intake hypothesis) or reallocate energy from maintenance processes (compensation hypothesis). The first hypothesis predicts a positive relationship between basal metabolic rate (BMR) and energy expenditure (DEE) because of the higher energy demands of the metabolic system at rest. The second hypothesis predicts a trade-off between different body functions, with a reduction of the BMR as a way to compensate for increased daytime energetic expenditure. We experimentally manipulated the workload of wild pied flycatchers by adding or removing chicks when chicks were 2 and 11 days old. We then measured the feeding frequency (FF), DEE and BMR at day 11, allowing us to assess both short- and long-term effects of increased workload. The manipulation at day 2 caused an increase in FF when broods were enlarged, but no response in DEE or BMR, while the manipulation at day 11 caused an increase in FF, no change in DEE and a decrease in BMR in birds with more chicks. Our results suggest that pied flycatchers adjust their workload but that this does not lead to a higher BMR at night (no support for the increased-intake hypothesis). In the short term, we found that birds reallocate energy with a consequent reduction of BMR (evidence for the compensation hypothesis). Birds thus resort to short-term strategies to increase energy expenditure, which could explain why energy expenditure and hard work are not always correlated in birds.

Evolutionary and demographic consequences of phenological mismatches.

Climate change has often led to unequal shifts in the seasonal timing (phenology) of interacting species, such as consumers and their resource, leading to phenological 'mismatches'. Mismatches occur when the time at which a consumer species's demands for a resource are high does not match with the period when this resource is abundant. Here, we review the evolutionary and population-level consequences of such mismatches and how these depend on other ecological factors, such as additional drivers of selection and density-dependent recruitment. This review puts the research on phenological mismatches into a conceptual framework, applies this framework beyond consumer-resource interactions and illustrates this framework using examples drawn from the vast body of literature on mismatches. Finally, we point out priority questions for research on this key impact of climate change.

Phenological mismatch drives selection on elevation, but not on slope, of breeding time plasticity in a wild songbird.

Phenotypic plasticity is an important mechanism for populations to respond to fluctuating environments, yet may be insufficient to adapt to a directionally changing environment. To study whether plasticity can evolve under current climate change, we quantified selection and genetic variation in both the elevation (RNE ) and slope (RNS ) of the breeding time reaction norm in a long-term (1973-2016) study population of great tits (Parus major). The optimal RNE (the caterpillar biomass peak date regressed against the temperature used as cue by great tits) changed over time, whereas the optimal RNS did not. Concordantly, we found strong directional selection on RNE , but not RNS , of egg-laying date in the second third of the study period; this selection subsequently waned, potentially due to increased between-year variability in optimal laying dates. We found individual and additive genetic variation in RNE but, contrary to previous studies on our population, not in RNS . The predicted and observed evolutionary change in RNE was, however, marginal, due to low heritability and the sex limitation of laying date. We conclude that adaptation to climate change can only occur via micro-evolution of RNE, but this will necessarily be slow and potentially hampered by increased variability in phenotypic optima.

Environmental coupling of heritability and selection is rare and of minor evolutionary significance in wild populations.

Predicting the rate of adaptation to environmental change in wild populations is important for understanding evolutionary change. However, predictions may be unreliable if the two key variables affecting the rate of evolutionary change-heritability and selection-are both affected by the same environmental variable. To determine how general such an environmentally induced coupling of heritability and selection is, and how this may influence the rate of adaptation, we made use of freely accessible, open data on pedigreed wild populations to answer this question at the broadest possible scale. Using 16 populations from 10 vertebrate species, which provided data on 50 traits (relating to body mass, morphology, physiology, behaviour and life history), we found evidence for an environmentally induced relationship between heritability and selection in only 6 cases, with weak evidence that this resulted in an increase or decrease in the expected selection response. We conclude that such a coupling of heritability and selection is unlikely to strongly affect evolutionary change, even though both heritability and selection are commonly postulated to be dependent on the environment.