Suzanne Bonamour
email: suzanne.bonamour-du-tartre@cefe.cnrs.fr

Luis-Miguel Chevin

Anne Charmantier

Céline Teplitsky

Phenotypic plasticity in

Keywords: phenotypic plasticity, climate change, cue reliability, blue tits, laying date

Introduction

Global change encompasses modifications of the environment both at a global scale (e.g. climate) and at a local scale but in so large proportion that the whole planet is affected (e.g. urbanization, invasive species). It thus represents complex modifications of the environments in which populations face radically new conditions (e.g. pesticides, roads), but also new modalities of historically known environments (e.g. extreme climatic events, invasive predator). Studying whether and how wild organisms adapt to these rapid environmental changes is both an opportunity for in-depth evolutionary ecology scrutiny, and a societal challenge. Phenotypic plasticity is a major mechanism of response to environmental variability, which may allow organisms to cope with rapid environmental changes, including global change. It has indeed been identified as the main mechanism of phenotypic change in response to climate change [START_REF] Gienapp | Climate change and evolution: disentangling environmental and genetic responses[END_REF][START_REF] Merilä | Climate change, adaptation, and phenotypic plasticity: The problem and the evidence[END_REF][START_REF] Matesanz | Global change and the evolution of phenotypic plasticity in plants[END_REF] and other human-induced rapid changes such as urbanization [START_REF] Alberti | Urban driven phenotypic changes: Empirical observations and theoretical implications for eco-evolutionary feedback[END_REF]. Phenotypic plasticity is also suspected to have an important role in colonization of new environments, geographic range shifts, and the success of invasive species [START_REF] Davidson | Do invasive species show higher phenotypic plasticity than native species and, if so, is it adaptive? A meta-analysis[END_REF][START_REF] Lande | Evolution of phenotypic plasticity in colonizing species[END_REF]. However, whether or not plasticity will aid adaptation and population persistence in a new environment depends on whether it is adaptive there. When it is adaptive, phenotypic plasticity can increase the probability of population persistence [START_REF] Chevin | Adaptation, plasticity, and extinction in a changing environment: towards a predictive theory[END_REF], as exemplified in two great tit (Parus major) populations where phenotypic plasticity is predicted to increase the likelihood of population persistence under various scenarios of climate change [START_REF] Vedder | Quantitative assessment of the importance of phenotypic plasticity in adaptation to climate change in wild bird populations[END_REF][START_REF] Gienapp | Predicting demographically sustainable rates of adaptation: can great tit breeding time keep pace with climate change?[END_REF]. However, these predictions are based on models that generally consider a "static" estimate of plasticity, without taking into account potential changes in reaction norms due to plasticity evolution, or the possibility that the adaptive nature of plasticity may be altered under new environmental conditions [START_REF] Ghalambor | Adaptive versus non-adaptive phenotypic plasticity and the potential for contemporary adaptation in new environments[END_REF], as is the case under global change.

A core prerequisite for adaptive predictive plasticity is the existence of reliable cues, i.e. accurate environmental information about future selection on the expressed plastic phenotype [START_REF] Moran | The evolutionary maintenance of alternative phenotypes[END_REF][12][START_REF] Botero | Evolutionary tipping points in the capacity to adapt to environmental change[END_REF][START_REF] Scheiner | The genetics of phenotypic plasticity: XII: Temporal and spatial heterogeneity[END_REF]. While the definition of cues is fairly straightforward in theory, their empirical characterization may be a very arduous task (electronic supplementary material, Box S1). Yet it is a crucial step in order to understand how global change can affect adaptive plasticity and its consequences on population biology. For instance, theory predicts that if cue reliability decreases, plasticity may be selected to decrease [START_REF] Gavrilets | The genetics of phenotypic plasticity. V. Evolution of reaction norm shape[END_REF], or else may drive populations to extinction by increasing the expected load caused by phenotypic mismatch with the optimum in a fluctuating environment [START_REF] Reed | Phenotypic plasticity and population viability: the importance of environmental predictability[END_REF][START_REF] Chevin | Phenotypic plasticity in evolutionary rescue experiments[END_REF][START_REF] Ashander | Predicting evolutionary rescue via evolving plasticity in stochastic environments[END_REF]. This is just one possible scenario, but it illustrates the point that one of the most central tasks for understanding and predicting the ecological and evolutionary roles of phenotypic plasticity under global change is to decipher which environmental cues are used by organisms, why (i.e. which selective pressures), and the extent to which global change is likely to affect cue reliability.

To address these questions, we will start with a brief overview of the nature and evolution of cues in wild populations. We will then outline the potentially complex effects of global change on plasticity, through its effect on the availability and reliability of cues. Finally, we will illustrate some of these aspects with a case study comparing plasticity of phenology in four blue tits populations from contrasted habitats.

Nature and evolution of cues for phenotypic plasticity

For plasticity to be adaptive, the environment influencing the development and/or expression of a particular phenotypic trait (or set thereof) needs to be a reliable cue for the selective pressure on this trait. In a variable environment, this reliability implies that the cue(s) should predict the selective environment, i.e. the environmental conditions driving natural selection on the expressed plastic trait (figure 1 and [START_REF] Moran | The evolutionary maintenance of alternative phenotypes[END_REF]12,[START_REF] Scheiner | The genetics of phenotypic plasticity: XII: Temporal and spatial heterogeneity[END_REF]), so that the plastic response leads to increased fitness [START_REF] Ghalambor | Adaptive versus non-adaptive phenotypic plasticity and the potential for contemporary adaptation in new environments[END_REF]. Reliable cues do not require that the same environmental variable affects both the expression of a trait and selection on this trait: any environmental variable that is correlated with the environment of selection through space or time can act as a reliable cue for phenotypic plasticity (e.g. red/infra-red ratio signaling competitor presence [START_REF] Griffith | Shade tolerance plasticity in response to neutral vs green shade cues in Polygonum species of contrasting ecological breadth[END_REF]). In fact, environmental variables that are most informative about the selective environment at the time of phenotypic determination may differ from the selective environment itself at this time [START_REF] Chevin | Evolution of environmental cues for phenotypic plasticity[END_REF]. For instance, rain may be a better predictor of the peak of food abundance at time t+ (where  is the lag between phenotypic determination and selection on the trait) than is food abundance itself at time t. More generally, organisms are expected to respond to the linear combination of environmental variables that best predicts variation in the phenotypic optimum [START_REF] Gavrilets | The genetics of phenotypic plasticity. V. Evolution of reaction norm shape[END_REF][START_REF] Chevin | Evolution of environmental cues for phenotypic plasticity[END_REF]. Very few case studies have investigated how an optimum phenotype is predicted by one environmental variable [START_REF] Chevin | Estimating the variation, autocorrelation, and environmental sensitivity of phenotypic selection[END_REF], and probably even fewer address the joint effect of several environmental variables. In contrast, many empirical studies have focused on the multidimensional aspect of plasticity, and highlighted that organisms respond to different environmental factors in an integrated way [START_REF] Westneat | Complex interactions among temporal variables affect the plasticity of clutch size in a multi-brooded bird[END_REF][START_REF] Sgrò | What can plasticity contribute to insect responses to climate change?[END_REF][START_REF] Valladares | Ecological limits to plant phenotypic plasticity[END_REF][START_REF] Polgar | Leaf-out phenology of temperate woody plants: From trees to ecosystems[END_REF]. Responding to multiple cues allows organisms to fine tune their phenotype in complex environments: for example, preys can be simultaneously confronted to different types of predators [START_REF] Black | Predator-induced phenotypic plasticity in Daphnia pulex: Life history and morphological responses to Notenecta and Chaoborus[END_REF][START_REF] Teplitsky | Hierarchical responses of tadpoles to multiple predators[END_REF], or plants face the need to respond to herbivory as well as to competition for light and water availability [START_REF] Valladares | Ecological limits to plant phenotypic plasticity[END_REF].

The aspect of cue reliability that influences the evolution of plasticity (in models of linear reaction norms) is the regression slope of the environment of selection on the environment of development [START_REF] Gavrilets | The genetics of phenotypic plasticity. V. Evolution of reaction norm shape[END_REF][START_REF] De Jong | Unpredictable selection in a structured population leads to genetic differenciation in evolved reaction norms[END_REF]. When the time lag between these two environments increases, cues are expected to become less informative [START_REF] Bradshaw | Evolutionary significance of phenotypic plasticity in plants[END_REF][START_REF] Padilla | Plastic inducible morphologies are not always adaptive: The importance of time delays in a stochastic environment[END_REF], to an extent that depends on the time scale of environmental fluctuations: strong temporal autocorrelation of the environment may allow the evolution of plastic responses to early cues. Similarly, dispersal between cue perception and selection causes a spatial lag analog to the time lag [START_REF] Moran | The evolutionary maintenance of alternative phenotypes[END_REF], with an effect on plasticity evolution that depends on the degree of spatial autocorrelation over the typical dispersal distance. However, even a cue that reliably predicts the environment of selection can lead to maladaptive plastic responses, if there is a development lag (time between cue perception and phenotype expression [START_REF] Scheiner | The genetics of phenotypic plasticity: XII: Temporal and spatial heterogeneity[END_REF]), such that the "set value" corresponding to that cue is not reached when selection operates on the trait. When the environment fluctuates on time scales longer than the generation time and dispersal is low, information about parental environment can be a reliable cue for offspring, and transgenerational plasticity is expected to evolve (figure 1, [START_REF] Mousseau | The adaptive significance of maternal effects[END_REF][START_REF] Leimar | The evolution of transgenerational integration of information in heterogeneous environments[END_REF][START_REF] Mcnamara | Detection vs. selection: integration of genetic, epigenetic and environmental cues in fluctuating environments[END_REF]). Transgenerational effects may be especially important for traits that are expressed and fixed early in the ontogeny (including many morphological traits), because offspring may lack sensory abilities during development [START_REF] Uller | Developmental plasticity and the evolution of parental effects[END_REF], and have fewer time to integrate information about the environment. For example, in the nematode Caenorhabditis elegans, maternal provisioning in glycogen increases offspring fitness under anoxia. Using experimental evolution, Dey et al [START_REF] Dey | Adaptation to temporally fluctuating environments by the evolution of maternal effects[END_REF] showed that when normoxia (i.e. normal levels of oxygen) and anoxia (<1% oxygen) were predictably alternating, mothers experiencing normoxia increased glycogen provisioning to their embryos (and vice versa for anoxia), demonstrating the evolution of maternal effects under predictable variation of the environment. In marine sticklebacks (Gasterosteus aculeatus) where body size is a key component of fitness, maternal effects could help adaptation to increased temperature under climate change, as offspring from mothers acclimated to 21°C were larger in this warm environment than conspecifics from mothers acclimated to 17°C, thanks to different mitochondrial performances [START_REF] Shama | Transgenerational plasticity in marine sticklebacks: Maternal effects mediate impacts of a warming ocean[END_REF].

Recently, there has been a growing interest in understanding theoretically how cue variability and reliability at different time scales interact to shape plastic responses. Organisms are predicted to integrate various sources of information, including cues from the current environment (for labile traits that change continuously in life), earlier cues perceived during development (early-life or carry-over effects), and transgenerational cues from the environment to which previous generations were exposed (figure 1, [START_REF] Mcnamara | Detection vs. selection: integration of genetic, epigenetic and environmental cues in fluctuating environments[END_REF]), including the parental phenotype itself [START_REF] Kuijper | When to rely on maternal effects and when on phenotypic plasticity?[END_REF]. The relative weight of these different types of cues is expected to vary depending on their predictive power over different time scales (related to the temporal autocorrelation of the environment), and the strength of selection [START_REF] Leimar | The evolution of transgenerational integration of information in heterogeneous environments[END_REF][START_REF] Mcnamara | Detection vs. selection: integration of genetic, epigenetic and environmental cues in fluctuating environments[END_REF][START_REF] Kuijper | When to rely on maternal effects and when on phenotypic plasticity?[END_REF].

Cue variation through space and/or time can lead to variation in plasticity, with an alternation of different ecotypes. A telling example of within-population variation in plasticity comes from the wild mosquitofish (Gambusia holbrooki). Fish born in the early spring (development at cool temperature) are expected to experience more variable environments than fish born in early summer (development at warm temperature). Accordingly, metabolic rate is plastic in fish born in early spring, but not in early summer fish [START_REF] Seebacher | Regulation of thermal acclimation varies between generations of the short-lived mosquitofish that developed in different environmental conditions[END_REF]. Following on from such encouraging case studies, more work is needed to assess, both theoretically and empirically, whether and how much variation in cue variability and predictability drives within-species variation in plasticity.

How may global change influence cues and plasticity?

Under global change, several aspects of the environment are expected to be modified, notably through the emergence of new environments (due to local habitat change or dispersal in new habitats), alteration of covariances among environmental variables, and increased climatic variability. These environmental changes can be gradual, but abrupt changes are also increasingly likely because of anthropogenic pressures [START_REF] Ratajczak | Abrupt change in ecological systems: inference and diagnosis[END_REF], e.g. in the form of extreme climate events such as floods, droughts, or storms. Altogether, environmental changes may affect the expression and evolution of plasticity through changes in cue reliability, cue perception and interpretation, and development of the phenotype.

(a) Changes in the availability and reliability of cues

The simplest way environmental change may affect cue use is by disrupting their detection, either by degrading the signal or by disturbing organisms' sensory abilities. For example, juvenile damselfish (Pomacentrus wardi) no longer respond to predator cues in degraded environments (dead coral reefs), and the likely mechanism is a degradation of the conspecific alarm cue [START_REF] Lönnstedt | Degraded environments alter prey risk assessment[END_REF]. Similarly, eutrophication can impair visual signals, leading to maladaptive plasticity in male sticklebacks (Gasterosteus aculeatus), as their increased investment in courtship behavior is not correlated with female interest [START_REF] Candolin | Population responses to anthropogenic disturbance: Lessons from threespined sticklebacks Gasterosteus aculeatus in eutrophic habitats[END_REF]. Pervasive electrosmog, i.e. human-made electromagnetic noise, completely disrupts the magnetic compass orientation in the migratory European robin Erithacus rubecula [START_REF] Engels | Anthropogenic electromagnetic noise disrupts magnetic compass orientation in a migratory bird[END_REF].

When cues are still available, their reliability can be affected if the new environment is similar in some ways to a known environment. In the most extreme cases, a new environment is generating a signal similar to a previously known cue, but completely uncorrelated with both the original cue and the environment of selection. For example asphalt and ponds polarize light in the same way, leading mayflies to lay their eggs on the road rather than ponds [START_REF] Kriska | Why do mayfly lay their eggs en masse on dry asphalt roads?[END_REF], and insecticides can trigger costly morphological defenses in Daphnia ambigua in the absence of predators [START_REF] Hanazato | Pesticides as chemical agents inducing helmet formation in Daphnia ambigua[END_REF][START_REF] Riessen | Costs of predator-induced morphological defences in Daphnia[END_REF]. Such maladaptive responses based on cues that were previously reliable are termed evolutionary traps [START_REF] Robertson | Ecological novelty and the emergence of evolutionary traps[END_REF][START_REF] Schlaepfer | Ecological and evolutionary traps[END_REF]. These traps will be all the more effective as the previous cue was highly reliable [START_REF] Robertson | Ecological novelty and the emergence of evolutionary traps[END_REF].

On the other end, the cue can be sampled from the same historical environmental variable, but the correlation between this variable and the selective environment may have decreased or entirely vanished. This is expected under climate change, for instance because of increased climate stochasticity within years, or because the average temperature is not changing at the same rate for all seasons [START_REF] Marrot | Current spring warming as a driver of selection on reproductive timing in a wild passerine[END_REF][START_REF] Williams | Cold truths: How winter drives responses of terrestrial organisms to climate change[END_REF]. For example, the yellow bellied marmot (Marmota flaviventris) has gradually emerged earlier from hibernation (study led between 1975 and 1999) because of warmer air temperature earlier in the spring. However, the date of snowmelt has not changed in the Rocky Mountains, leading to lower foraging opportunities at emergence. While warmer air was historically a good predictor of snow melt, climate change has led to a mismatch between air temperature and snow cover [START_REF] Inouye | Climate change is affecting altitudinal migrants and hibernating species[END_REF].

Because species are embedded in networks of ecological interactions, loss of cue reliability may arise from species responding differently to environmental changes, or interacting species responding to different cues. Changes in phenology are the most ubiquitous responses to climate change [START_REF] Parmesan | Ecological and evolutionary responses to recent climate change[END_REF], and the best documented cases of disruption of ecological networks are based on a loss of synchrony between ecological interactors, such as predators and their preys, or plants and their pollinators [START_REF] Visser | Shifts in phenology due to global climate change: the need for a yardstick[END_REF][START_REF] Byers | Studying plant-pollinator interactions in a changing climate: A review of approaches[END_REF]. If interacting species are responding to the same cue, their phenological mismatch may increase because their reaction norms are different. This is generally the case in trophic interactions, where consumers tend to display shallower (flatter) reaction norms to climate than their resources [START_REF] Thackeray | Phenological sensitivity to climate across taxa and trophic levels[END_REF]. This occurs because consumers generally have longer generation times than resources (due to larger body size), and thus longer developmental lags for the expression of the plastic trait, which reduces their ability to predict the environment of selection, as temporal autocorrelation affecting cue reliability decreases on average with longer time lags [START_REF] Gienapp | Why climate change will invariably alter selection pressures on phenology[END_REF]. Another possibility is that interacting species are actually using different cues that undergo different time trends. For example, four tree species of Prunus in Japan are flowering earlier due to increased temperatures during the time window most related to their phenology, yet temperature remains unchanged during the time window related to the phenology of the butterfly Pieris rapae that uses them as host plants. This leads to strong plasticity of phenology in trees but not in butterflies [START_REF] Doi | Heterogeneous intra-annual climatic changes drive different phenological responses at two trophic levels[END_REF], resulting in an increasing plantinsect phenology mismatch due to warming (see also our case study below). Some mechanisms may mitigate the loss of reliability in some cues. First, the potentially negative effects of loss of reliability of a particular environmental variable as a cue could be greatly alleviated if organisms are using multiple cues, because partial redundancy in information from different environmental sources may increase the robustness of the signal [START_REF] Dore | The role of complex cues in social and reproductive plasticity[END_REF]. For example in the above-mentioned case of the coral reef damselfish, the loss of olfactory cue in dead coral water can be compensated by the use of visual cues [START_REF] Lönnstedt | Degraded environments alter prey risk assessment[END_REF]. Similarly, we could expect that the mismatch in trophic webs could be alleviated if species are using each other's presence, abundance, or phenotype as cues (as investigated empirically by Phillimore et al. [START_REF] Phillimore | Dissecting the contributions of plasticity and local adaptation to the phenology of a butterfly and its host plants[END_REF] with butterflies and their host plants), rather than using a unique climatic variable such as temperature alone. Second, the expression of maladaptive responses to an unreliable cue could be reduced through learning, for example if females are able to recalibrate their cue use based on past experience [START_REF] Grieco | Evidence for the effect of learning on timing of reproduction in blue tits[END_REF].

(b) Consequences on the expression and evolution of plasticity

It is difficult to predict in general whether and how much plasticity should be adaptive in novel environments. It has been suggested that reaction norms should have random and erratic shapes in environments that were never or rarely encountered previously (as is expected to occur under global change), due to the absence of past selection in these environments [START_REF] Ghalambor | Adaptive versus non-adaptive phenotypic plasticity and the potential for contemporary adaptation in new environments[END_REF]. However, more quantitative arguments can be made by taking into account two realistic features of plasticity in new environments.

First, environments are seldom entirely new; instead, major environmental shifts mostly mean that previously (perhaps very) rare environments have become common. How rare a new environment has been prior to the shift determines how relaxed selection has been in this environment, and the opportunity there has been for genetic drift to produce erratic reaction norms there [START_REF] Chevin | Evolution of phenotypic plasticity in extreme environments[END_REF]. Second, reaction norms do not evolve completely freely, even in environments where they have not been under selection, because their shapes are generally constrained to some extent by genetic correlations of trait values across environments. Considering these two points makes it likely that reaction norms remain at least partly adaptive over the new environmental range in the absence of major changes in the adaptive landscape, and simple quantitative statements can be made about this question based on the frequency of extremes and the correlations of trait values across environments [START_REF] Chevin | Evolution of phenotypic plasticity in extreme environments[END_REF].

Retaining partially adaptive plastic responses in novel environments also requires that cues remain reliable. In the specific case of extreme climatic events, whether such cues exist is still unknown. One recent study tackling this issue found no evidence that the Eurasian oystercatchers (Haematopus ostralegus) uses cues (lunar cycle, past and current water height) to avoid flooding of their nests during extreme precipitation events [START_REF] Bailey | No phenotypic plasticity in nest-site selection in response to extreme flooding events[END_REF]. Furthermore, global change generally involves variation of more than one environmental variable, and an alteration of the correlation structure of these variables within and between times. This may trigger conflicting ecological demands on organisms [START_REF] Valladares | Ecological limits to plant phenotypic plasticity[END_REF], such that new combinations of environmental variables become the best predictors of the optimum phenotype. The ability of organisms to track these altered patterns of change in the optimum phenotype is predicted to depend on their ability to evolve a new combination of environmental variables used as a cue [START_REF] Chevin | Evolution of environmental cues for phenotypic plasticity[END_REF]. This theoretical prediction has not been tested yet, but a couple of studies suggest that cue use can evolve relatively fast. A study in the pitcher plant mosquito (Wyeomia smithii) suggests that the critical photoperiod for diapause induction has shifted in 5 years towards shorter day length (i.e. later in the season) due to longer growing seasons [START_REF] Bradshaw | Genetic shift in photoperiodic response correlated with global warming[END_REF]. Antipredator responses of the Iberian waterfrog (Pelophylax perezi) tadpoles to an invasive species (crayfish, Procambarus clarkia) is also a likely case of cue evolution. Tadpoles from this species have the ability to build morphological defenses against native predators [START_REF] Nunes | Morphological and life-history responses of anurans to predation by an invasive crayfish: An integrative approach[END_REF] but only tadpoles from populations coexisting with P. clarkia since 10-15 generations display anti-predator response to this new predator [START_REF] Nunes | Rapid evolution of constitutive and inducible defenses against an invasive predator[END_REF].

Changes in cue reliability are expected to have a major impact on the evolution of plasticity, including a transition between predictive plasticity and bet hedging as adaptive strategies [START_REF] Botero | Evolutionary tipping points in the capacity to adapt to environmental change[END_REF][START_REF] Tufto | Genetic evolution, plasticity, and bet-hedging as adaptive responses to temporally autocorrelated fluctuating selection: A quantitative genetic model[END_REF]. In the context of climate change and the associated increased variability of temperature, decreased predictability of the environment could for example impact the adaptiveness of transgenerational plasticity, and select for a strategy of bet hedging, where parental effects increase variance in offspring phenotype rather than alter their mean phenotypes [START_REF] Donelson | Transgenerational plasticity and climate change experiments: Where do we go from here?[END_REF].

In terms of predictive plasticity, the reliability of the cue determines the steepness of the optimal reaction norm: the less reliable the cue, the shallower the reaction norm that evolves at equilibrium, as compared to "perfect" plasticity, defined by how the optimum phenotype is affected by the environment [12,[START_REF] Lande | Adaptation to an extraordinary environment by evolution of phenotypic plasticity and genetic assimilation[END_REF]. In the case of an abrupt change in cue reliability, causing previously adaptive plasticity to become maladaptive (for instance because it causes overshoots of the optimum as in [START_REF] Reed | Phenotypic plasticity and population viability: the importance of environmental predictability[END_REF]), plasticity is expected to evolve towards a shallower reaction norm that matches the current level of cue reliability. Such plasticity shift may cause evolutionary rescue by reducing the load caused by stochastic fluctuations in the environment [START_REF] Chevin | Phenotypic plasticity in evolutionary rescue experiments[END_REF][START_REF] Ashander | Predicting evolutionary rescue via evolving plasticity in stochastic environments[END_REF].

Case study: local variation of cues in Mediterranean blue tits and their preys

Investigating the mechanisms of plastic response to climate change and their limits in the wild is notoriously difficult (electronic supplementary material, Box S1). However, an exceptionally wellcharacterized research model of adaptive plasticity in wild populations is the earlier egg-laying date of temperate forest insectivorous passerines during warmer springs. This earlier breeding phenology allows birds to track the phenology of the caterpillars they use as main food source for their offspring [START_REF] Chevin | Estimating the variation, autocorrelation, and environmental sensitivity of phenotypic selection[END_REF][START_REF] Visser | Shifts in caterpillar biomass phenology due to climate change and its impact on the breeding biology of an insectivorous bird[END_REF]. In the following section, we investigate in some detail how plastic phenological change relates to environmental cues and their reliability in four Mediterranean populations of blue tits (Cyanistes caeruleus).

The four wild populations of blue tits have been monitored for 20 to 42 years (electronic supplementary material, table S1). Two populations breed in forests dominated by the Deciduous downy oak Quercus pubescens (the D-Rouvière and the D-Muro populations, in the French mainland close to Montpellier and Corsica respectively, where 'D' stands for deciduous) and two in forests dominated by the evergreen holm oak Quercus ilex (the E-Muro and the E-Pirio populations, both in Corsica; 'E' stands for evergreen forest [START_REF] Blondel | A thirtyyear study of phenotypic and genetic variation of Blue tits in Mediterranean habitat mosaics[END_REF]). These ecological characteristics have induced local adaptation in birds [START_REF] Charmantier | Mediterranean blue tits as a case study of local adaptation[END_REF], particularly in their life-history traits. Average phenology of two of these populations has been changing over time: mean laying date is increasingly early (more than three days per decade) for the populations breeding in evergreen, but not in deciduous forests (figure 2a and electronic supplementary material, table S2a), despite similar trends of climate warming in the four sites (electronic supplementary material, table S2c and figure S1; note that table S1 reports trends in temperature over the spring -from April to end-June -while figure S1 presents trends in temperature over the entire year). This pattern is all the more striking as the E-Muro and D-Muro populations are only 5.6 km apart [START_REF] Charmantier | Mediterranean blue tits as a case study of local adaptation[END_REF], and connected by gene flow [START_REF] Szulkin | Population genomic footprints of fine scale differentiation between habitats in Mediterranean blue tits[END_REF].

One plausible explanation behind the different time trends in laying date could be local variation of environmental cues for phenotypic plasticity. Females from different populations could be sensitive to different climate variables (e.g. amount of rain, mean temperature or presence of extreme climatic events), or to the same variables but at different time periods during the year (e.g. early vs. late spring). In order to assess the population-specific cues for laying date plasticity in blue tits, we implemented sliding windows analyses using the package climwin [START_REF] Bailey | climwin : An R toolbox for climate window analysis[END_REF][START_REF] Van De Pol | Identifying the best climatic predictors in ecology and evolution[END_REF]. For each population, we tested several weather variables: mean, maximum, minimum for both daily temperature and daily rain, and positive Extreme Climatic Events of temperature (ECE hereafter, electronic supplementary material Box S2 for details on the methods). Climate data were obtained from the national meteorological stations of Saint Martin de Londres for the mainland (about 24 Km from D-Rouvière) and Calvi for Corsica (9-19 km from the three Corsican study sites). Temperatures from the meteorological stations are highly correlated to local temperature in study sites (electronic supplementary material, table S3).

In three of the populations (D-Rouvière, D-Muro and E-Pirio), the climatic factor that most influences variation of mean laying date is the mean temperature (electronic supplementary materials, table S4) explaining between 54 and 75% of variability in laying date (table 1a). "Evergreen" populations are also sensitive to the number of ECE. In E-Pirio, this effect is additive with the mean temperature effect, but not highly robust. In turn, the number of ECE within a window is the best predictor of average laying date in the E-Muro population (electronic supplementary materials, table S4), but the simultaneous use of information from mean temperature cannot be excluded (mean temperature: -3.47 ± 0.88, t =-3.94, p = 0.001, ECE: -1.59 ± 0.42, t = -3.79, p = 0.002, r² = 0.75). We found no effect of other climate variables such as rainfall (electronic supplementary materials, table S4). Hence, although some combination of other factors (e.g. bud development [START_REF] Bourgault | Spring vegetation phenology is a robust predictor of breeding date across broad landscapes: A multi-site approach using the Corsican blue tit (Cyanistes caeruleus)[END_REF] or population density [START_REF] Bourret | Multidimensional environmental influences on timing of breeding in a tree swallow population facing climate change[END_REF]) likely play a role, temperature is a key driver behind the plasticity of laying date. This is in line with other correlative studies in great tits, showing that the average spring temperature explains more than 50% of the variability of laying date [START_REF] Visser | Shifts in caterpillar biomass phenology due to climate change and its impact on the breeding biology of an insectivorous bird[END_REF][START_REF] Charmantier | Adaptive phenotypic plasticity in response to climate change in a wild bird population[END_REF], and with experiments showing that laying date responds to temperature treatments [START_REF] Visser | Temperature has a causal effect on avian timing of reproduction[END_REF]. The underlying physiological mechanisms remain elusive, as no effect of temperature on the neuroendocrine system linked to reproduction has been detected yet [START_REF] Schaper | Individual variation in avian reproductive physiology does not reliably predict variation in laying date[END_REF][START_REF] Caro | Is microevolution the only emergency exit in a warming world? Temperature influences egg laying but not its underlying mechanisms in great tits[END_REF]. Despite many studies on bird breeding phenology, knowledge about mechanisms driving phenological plasticity are still scarce, and the relative roles of direct and indirect effects of temperature (e.g. for the latter, a constraint on the timing of egg laying mediated by food abundance or quality) are unknown [START_REF] Caro | The Case of the Missing Mechanism: How Does Temperature Influence Seasonal Timing in Endotherms[END_REF].

The populations differ in the specific time window of mean temperature influencing laying date. Females from populations in deciduous habitats use a cue based on a long time period encompassing winter and early spring, while females from the E-Pirio population use a shorter and later time period corresponding to a month in spring (figure 3, electronic supplementary materials, table S5). The pattern is qualitatively similar but not robust in E-Muro (figure 3, electronic supplementary material, table S4 andS5). Larger windows in deciduous populations could be related to the correlation between the best window and other possible windows remaining high across a larger range of time periods than in evergreen populations (electronic supplementary materials Figure S2). A formal statistical comparison between the windows used in D-Rouvière and E-Pirio shows significant differences between the two populations in terms of ordinal calendar (electronic supplementary materials Box S2).

The use of different time windows in the four bird populations leads to contrasted rates of change of the cue over the study period. In the "Evergreen" populations, where the window includes only spring temperature, the mean temperature inside the cue window increases over the 1991-2017 period (E-Muro: 0.4 ± 0.07 °C/decade, t-value = 6.18, p-value <0.0001; E-Pirio: 0.6 ± 0.08, t-value = 7.81, p-value <0.0001, figure 2b), but there is no such increase in the two other sites where the windows include information from spring and winter (D-Rouvière: 0.03 ± 0.3 °C/decade, t-value = 0.11, p-value = 0.91; D-Muro 0.22 ± 0.12, t-value = 1.76, p-value = 0.09). This is in line with the very slow or absent rate of climate change in winter as compared to spring in these Mediterranean areas (electronic supplementary material, figure S1). Climate change is also associated with increasing frequency of ECE, and their number increases in the window used by the E-Muro population (0.49 ± 0.16 positive ECE/decade, t-value = 2.98, p-value = 0.005, see electronic supplementary material table S6 for other populations). Hence overall, even though spring warming has been pervasive in all populations (electronic supplementary materials, table S2c, figure S1), we only detect climate warming during the cue windows influencing bird phenology in E-Pirio and E-Muro. This difference of windows across habitats explains laying dates are advancing in the evergreen avian populations but not in the deciduous populations.

It is also possible to investigate the reliability of cues in this system. The fitness consequences of phenological responses of blue tits depend on the phenology of caterpillars of Tortrix viridana, the main food resource for nestlings [START_REF] Bańbura | Nestling diet variation in an insular Mediterranean population of blue tits Parus caeruleus: effects of years, territories and individuals[END_REF], thought to influence reproductive success through offspring survival. The synchronization of the nesting period (especially when nestlings reach the peak of their energetic needs, around 9 days after hatching [START_REF] Dias | Local specialization and maladaptation in the Mediterranean blue tit (Parus caeruleus)[END_REF]) with date of peak abundance in caterpillars depends on how reliably the environmental cue used by birds to modify their phenology predicts caterpillar peak date. We can thus measure the linear regression slope of the environment affecting selection on timing of breeding on the environment affecting timing of breeding itself, in order to quantify the reliability of the cue used by birds in a way that matches theoretical predictions [e.g. [START_REF] Gavrilets | The genetics of phenotypic plasticity. V. Evolution of reaction norm shape[END_REF][START_REF] De Jong | Unpredictable selection in a structured population leads to genetic differenciation in evolved reaction norms[END_REF]. Caterpillar phenology is studied in the four sites thanks to coprometers (i.e. 50*50 cm cloth square collecting the frass of caterpillars under oak trees). The mass of frass collected during a given time period allows us to estimate the abundance of caterpillars throughout the season. The peak date of caterpillar abundance is the day with the highest collected quantity of caterpillar frass. To evaluate the cue used by caterpillars, we used the same sliding window analysis as for birds. Dataset sizes allow caterpillar sliding windows analyses only in E-Pirio and D-Muro (at least 20 years of data are required [START_REF] Van De Pol | Identifying the best climatic predictors in ecology and evolution[END_REF]). Temperature explained ca. 60% of the variation in the date of caterpillar peak in both sites (table 1b, figure 3, electronic supplementary materials table S5). The reliability of the cue was high in both populations: 1.02 ± 0.04 (t-value = 25.17, p-value < 0.0001) in D-Muro and 0.70 ± 0.08 in E-Pirio (tvalue = 9.071, p-value < 0.0001; with lower reliability in E-Pirio than D-Muro (negative blue tit cue*E-Pirio interaction: -0.32 ± 0.09, t-value = -3.62, p-value = 0.0005).

Another quantity of interest is how well the cue used by birds predict the food peak itself. The slope of the regression of the caterpillar date on the cue used by blue tits is the reaction norm that is favoured by natural selection, if caterpillar date is taken as the optimum phenotype each year. This relationship is significantly negative in all sites (figure 2e, electronic supplementary material, table S7): early prey abundance is associated with high temperature in the cue window. However, the proportion of variation in the food peak that is captured by its relationship with the cue for plasticity is lower in E-Pirio (r² = 0.21) than in D-Muro (r² = 0.59, see discussion), and also moderate in the two other populations (D-Rouvière: r² = 0.33, E-Muro r² = 0.29). The strength of this relationship in D-Muro is similar as that found for a Dutch great tit population in Hoge Veluwe [START_REF] Visser | Shifts in caterpillar biomass phenology due to climate change and its impact on the breeding biology of an insectivorous bird[END_REF]. If the food peak correctly predicts the optimum laying date, this suggests that plasticity in response to temperature in the time windows we identified only allows tracking a moderate proportion of temporal fluctuations in this optimum. However, note that caterpillar abundance data in the D-Rouvière site should be used with caution: data collection stopped in 2002 because of frequent rainfalls, and 10 years of data may not be sufficient for analyses to be reliable.

A plasticity-mediated phenological mismatch between predators and their food source can arise from use of different cues, or from different responses to the same cues. Regardless of the mechanism, theory predicts that the slope of the consumer should be shallower than the slope of the resource, because cue reliability is expected to be lower in consumers [START_REF] Gienapp | Why climate change will invariably alter selection pressures on phenology[END_REF]. Accordingly, a wide scale analysis over 812 taxa revealed that secondary consumers have lower climate sensitivity than other groups (e.g. primary producers or consumers [START_REF] Thackeray | Phenological sensitivity to climate across taxa and trophic levels[END_REF]). In line with expectations, in E-Pirio birds are less sensitive to temperature than caterpillars (figure 2c,d, table 1, species*temperature interaction: -4.06 (± 1.28), t-value = -3.17, p-value = 0.002). However, in the D-Muro site, the slopes of the reaction norms are similar in birds and caterpillars (figure 2c,d,table 1, species*temperature interaction: -1.32 (± 1.02), t-value = -1.29, p-value = 0.20). The abundance of caterpillars in D-Muro is (much) higher than in the three other populations (electronic supplementary material, table S1). This make its likely that in other sites such as E-Pirio, birds are exploiting other resources, especially when caterpillars are rare [START_REF] Bańbura | Nestling diet variation in an insular Mediterranean population of blue tits Parus caeruleus: effects of years, territories and individuals[END_REF][START_REF] Tremblay | The effects of habitat quality on foraging patterns, provisionning rate and nestling growth in corsican blue tits Parus caeruleus[END_REF]. When food sources are diversified, flatter reaction norms in birds than in caterpillars may not imply that the reaction norms are not adaptive, but that birds are responding to a more multidimensional environment [START_REF] Chevin | Evolution of environmental cues for phenotypic plasticity[END_REF]. An important next step would thus be to estimate how the cue for avian phenology and the date of caterpillars jointly predict the annual optimum laying date [START_REF] Chevin | Estimating the variation, autocorrelation, and environmental sensitivity of phenotypic selection[END_REF].

However, this may not be sufficient to get a full understanding of the role of plasticity in adaptation in this system. For example, even in D-Muro where average laying date is tracking the caterpillar date accurately, birds are always late compared to the food peak (figure 2f, electronic supplementary material table S8), and there is strong directional selection for earlier laying date [START_REF] Porlier | Variation in phenotypic plasticity and selection patterns in blue tit breeding time: between-and withinpopulation comparisons[END_REF]. An integrative estimate of optimal laying date, using multiple components of fitness such as survival or recruitment success (as in e.g. [START_REF] Price | Directional selection and the evolution of breeding date in birds[END_REF][START_REF] Reed | Population Growth in a Wild Bird Is Buffered Against Phenological Mismatch[END_REF]), could help understand this lag, and whether or not it is adaptive.

Discussion and perspectives

Our literature review and case study make it clear that, if we wish to understand how global change will affect phenotypic plasticity and its consequences for population biology, we need a deeper understanding of the environmental mechanisms shaping variation in plasticity within species.

Perhaps the most difficult and needed measurements concern spatio-temporal changes in phenotypic selection , which underlie the evolution of phenotypic cues for plasticity. There have been repeated calls for measuring the environmental sensitivity of selection, notably as a way to identify the causes of natural selection [START_REF] Chevin | Adaptation, plasticity, and extinction in a changing environment: towards a predictive theory[END_REF][START_REF] Wade | The causes of natural selection[END_REF][START_REF] Maccoll | The ecological causes of evolution[END_REF]. This becomes a necessity when investigating the adaptiveness and evolution of plasticity, but the available methods and their applications are still limited [START_REF] Chevin | Estimating the variation, autocorrelation, and environmental sensitivity of phenotypic selection[END_REF][START_REF] Hunter | Quantification and decomposition of environment-selection relationships[END_REF]. In particular, we are not aware of an attempt to use multiple environmental variables as predictors of changes in an optimum phenotype, as a way to investigate selection on environmental cues for plasticity, consistent with predictions from theory [START_REF] Gavrilets | The genetics of phenotypic plasticity. V. Evolution of reaction norm shape[END_REF][START_REF] Chevin | Evolution of environmental cues for phenotypic plasticity[END_REF]. Manipulative experimental approaches are also a powerful yet underused tool for measuring selection on phenotypic plasticity and understanding the adaptive role of environmental cues. For instance, Schmitt et al. [START_REF] Schmitt | A test of the adaptive plasticity hypothesis using transgenic and mutant plants disabled in phytochrome-mediated elongation responses to neighbors[END_REF] experimentally shut down plastic stem elongation in response to crowding in plants, by constructing genetically modified lines of Brassica rapa that lack the photoreceptor phytochrome A involved in the detection of the relevant cue (ratio of wavelength characteristic of shading by other plants). They showed that the resulting lack of plastically induced elongation is detrimental in a crowded environment, while constitutive expression of stem elongation is detrimental in uncrowded environments. Other experiments of this kind, where cue perception/use is disrupted, or decoupled from the selective pressure, would yield extremely useful information, but are still too rare. Furthermore, evaluating whether phenotypic plasticity will remain adaptive in the face of global change requires an integrative approach at different levels. The first level of integration is across life stages. Environmental changes can affect phenotypic expression at different stages of the life cycle, so that the effects of global change can be mediated not only by current cues, but also by transgenerational, early-life, and carry-over effects (figure 1). The investigation of how these mechanisms of response to the environment interact is just starting [START_REF] Donelson | Transgenerational plasticity and climate change experiments: Where do we go from here?[END_REF][START_REF] Beaman | Evolution of Plasticity: Mechanistic Link between Development and Reversible Acclimation[END_REF], but it will provide important results on the dynamics of expression of plasticity. Second, there is also a strong need for an integrative approach at the community and/or ecosystem level. Interactions among species are a major source of selection for adaptive plasticity (e.g. synchrony of phenology, antipredator defences, competitive interactions). For example, because the migrant pied and collared flycatchers display shallower reaction norms to temperature than blue tits, they breed later in warmer environments, and suffer from competition with resident tit species [START_REF] Samplonius | Phenological sensitivity to climate change is higher in resident than in migrant bird populations among European cavity breeders[END_REF]. An important step would be to increase the number of studies comparing plasticity of interacting species in wild populations, notably examining how climate and environmental change are affecting cues and the expression of plasticity across trophic levels.

Comparing populations can provide insights as to how phenotypic plasticity and cue use are shaped by evolution. In our case study, geographically very close, and connected, populations occupying different habitats differed extensively in the cues affecting laying date plasticity, as well as the reliability of these cues as predictors of the putative selective pressure, namely the peak of abundance of the main food source (caterpillars). The populations from evergreen and deciduous habitats differed in the time windows -and to a lesser extent, in the environmental variables -that they use. The difference in time windows among populations is consistent with results from another study showing shorter windows for populations breeding in "late" habitats, such as high latitudes [START_REF] Phillimore | Passerines may be sufficiently plastic to track temperature-mediated shifts in optimum lay date[END_REF]. This suggests that these differences in time windows can be adaptive. However, detection of cues in wild populations relies on a correlative approach, and we cannot completely exclude that another unmeasured environmental factor correlated to temperature drives the plasticity of laying date, as well as local variation of cue windows. The ecological reality of these populations is certainly highly complex, and more information on environmental heterogeneity would be needed to assess which aspect of environmental variability can be linked to the difference in cue use. For instance, the focal populations vary in the abundance of caterpillars, suggesting that the required use of additional preys may partly explain the apparent difference among populations in the reliability of cues for predicting the selective environment. Altogether, understanding the phenology of blue tits (and other insectivorous passerines in temperate forests) requires integrating knowledge on -at least -three trophic levels: an insectivorous bird, a herbivorous caterpillar and an autotroph tree. In temperate regions, spring temperature seems to be the main environmental factor influencing hatching and development period in insects [START_REF] Van Asch | Phenology of Forest Caterpillars and Their Host Trees: The Importance of Synchrony[END_REF], following a period of low metabolic rates caused by ambient temperature during winter (or even a diapause period mainly driven by photoperiod [START_REF] Tauber | Insect seasonal cycles: genetics and evolution[END_REF]). Similarly, tree phenology often includes a chilling period (i.e. a period of cold temperature required to lift dormancy after winter [START_REF] Polgar | Leaf-out phenology of temperate woody plants: From trees to ecosystems[END_REF]). Chilling requirement and cold tolerance vary among species, more specifically deciduous oaks seem more sensitive to cold temperatures and freezing stress than evergreen oaks [START_REF] Cavender-Bares | Hydraulic properties and freezing-induced cavitation in sympatric evergreen and deciduous oaks with contrasting habitats[END_REF][START_REF] Cavender-Bares | Summer and winter sensitivity of leaves and xylem to minimum freezing temperatures: A comparison of cooccurring Mediterranean oaks that differ in leaf lifespan[END_REF]. Cold sensitivity of deciduous oaks suggest that winter temperature could partly predict tree and thus insect phenology, offering a potential explanation as to why winter temperatures are included in the cue window of birds breeding in deciduous oak forests. Finally, the use of extreme climatic events as an additional source of information in E-Muro (ECE) also requires further investigation before clear interpretations can be made. Our results overall highlight that ecological (biotic and abiotic components) variation can have strong effects on plasticity, but further empirical work is needed to be able to discuss our results from a more mechanistic point of view.

We find little evidence for multiple redundant climatic cues, but we cannot exclude that some other, unmeasured cues, such as tree phenology, could also be important in determining laying date [START_REF] Bourgault | Spring vegetation phenology is a robust predictor of breeding date across broad landscapes: A multi-site approach using the Corsican blue tit (Cyanistes caeruleus)[END_REF]. Although partial redundancy should increase the reliability of the information [START_REF] Dore | The role of complex cues in social and reproductive plasticity[END_REF], information acquisition is also expected to be costly [START_REF] Dewitt | Costs and limits of phenotypic plasticity[END_REF]. An open remaining question is thus the circumstances under which we can expect organisms to use several, partly redundant sources of information, and whether this could mitigate the effects of global change. For example, while temperature in late winter/early spring may not remain a reliable cue, tree phenology could provide more robust information on the long term. Detecting redundant environmental cues is challenging because of collinearity among factors leading to statistical issues (as is probably the case in our analysis), especially in short time series (but see [START_REF] Morrissey | Multiple Regression Is Not Multiple Regressions: The Meaning of Multiple Regression and the Non-Problem of Collinearity[END_REF]).

If a cue is no longer informative, there could be evolution of cue use. In general, we know little about the evolutionary potential of cue use but a couple of studies suggest some genetic variation exists for this. For example, in Daphnia magna, the integration of different cues for plasticity of life

Tables & Figures

Table 1 Average population reaction norms in response to the cues estimated with the sliding window analyses (see figure 3, table S5) for a) laying date of blue tits, and b) caterpillar date of peak abundance.

Depending on the best predictors of phenology detected by sliding window analysis, the reaction norms were estimated using either the mean temperature (MeanT) or the number of ECE during the best window. Sample sizes are given in the electronic supplementary material Table S1 S7); f) Synchrony between blue tits and caterpillars illustrated as the correlation between mean laying date of blue tits (± se) and caterpillar peak date (see table S8).

In all panels, the "Evergreen" sites are in green (dark green for E-Pirio, light green for E-Muro), and "Deciduous" sites are in blue (dark blue for D-Rouvière, light blue for D-Muro). All dates are ordinal dates with 1= 1 st January. In panel f, the red dashed line represents the expected optimal relationship between laying date and caterpillar date. In some cases, error bars are not visible because they are very small compared to the figure scale. [START_REF] Bailey | climwin : An R toolbox for climate window analysis[END_REF][START_REF] Van De Pol | Identifying the best climatic predictors in ecology and evolution[END_REF]. A good agreement between best and median window is an indicator that the window is precisely estimated.

  history traits varies among clones [100]. In A. thaliana the knowledge of genetic pathways involved in determining flowering time allows to understand variability in plastic responses and how multiple cues are integrated, and thus to predict the effects of environmental changes on plasticity [101,102]. Similar approaches in the wild seem very challenging for now, but laboratory experiments could provide very useful results in this yet little explored area.

Figure 1

 1 Figure 1Summarized diagram of the relationships between the environment, cues, and phenotypes in the framework of phenotypic plasticity. The multivariate environment can be used as a cue to initiate a plastic response in an individual (thin black arrows). The focal phenotype (on which selection is acting in our example, black thunderbolt) can be affected by current environmental conditions, earlier life events and the preceding parental generation. Information about past environments (either parental or from early life) includes different sources, such as the phenotype itself (at the same or another trait), or epigenetics. For phenotypic plasticity to be adaptive, cues in the successive periods and generations have to be correlated with the future selective environment (white arrows).

Figure 2 a

 2 Figure 2 a) Temporal trends in annual mean laying dates in the four Blue tit populations (mean ± se); b) Temporal trends in mean temperature (in °C ± se) for the cue windows used by each population; c) Blue tit reaction norms: annual mean laying date (± se) as a function of the estimated cue (in °C ± se); d) Caterpillar reaction norm (only for E-Pirio and D-Muro): annual caterpillar peak date as a function of the estimated cue (°C ± se); e) Cue reliability for blue tits measured as the correlation between the caterpillar peak date and the cue (tableS7); f) Synchrony between blue tits and caterpillars illustrated

Figure 3

 3 Figure 3Estimated mean temperature cue windows for blue tits (a, c, d, f) and caterpillars (b, e) in each population. The open circles represent the mean phenology (mean laying date or mean caterpillar peak date). Open and close day of the window are given in ordinal day, with 1= 1 st January. Coloured lines represent the best windows (from the best model, green for evergreen and blue for deciduous); black lines represent the median windows from the 95% confidence set (see electronic supplementary material, Box S2 and[START_REF] Bailey | climwin : An R toolbox for climate window analysis[END_REF][START_REF] Van De Pol | Identifying the best climatic predictors in ecology and evolution[END_REF]. A good agreement between best and median window is an indicator that

  .

	a)	Cue		Intercept			Slope		r²
			Estimates	t-value p-value	Estimates	t-value p-value
			(± se)			(± se)
	D-Rouvière MeanT		131.78	34.51	<0.0001	-4.42	-8.91	<0.0001	0.75
	(17 th Jan. -1 st April)	(± 3.82)			(± 0.50)
	D-Muro MeanT		137.61	20.31	<0.0001	-3.67	-5.92	<0.0001	0.59
	(4 th Feb. -7 th April)	(± 6.78)			(± 0.62)
	E-Muro MeanT		175.33	1.14	<0.0001	-5.13	-5.11	<0.0001	0.57
	(15 th March -27 th April)	(± 13.34)			(± 1.00)
	ECE		112.21	92.90	<0.0001	-2.41	-4.95	0.0001	0.55
	(23 th March -15 th April)	(± 1.21)			(± 0.49)
	E-Pirio MeanT		200.53	16.82	<0.0001	-5.10	-6.07	<0.0001	0.54
	(31th March -7 th May)	(± 11.92)			(± 0.84)
	b)
	D-Muro MeanT		171.22	19.01	<0.0001	-5.00	-6.15	<0.0001	0.61
	(11 th Feb. -7 th April)	(± 9.01)			(± 0.81)
	E-Pirio MeanT		283.38	14.14	<0.0001	-8.47	-6.46	<0.0001	0.58
	(13 th March -7 th June)	(± 20.04)			(± 1.31)

Acknowledgments

We thank Ally Phillimore and an anonymous reviewer for their constructive comments on the manuscript. We also thank Anssi Laurila, Samuel Caro and Pierre-Yves Henry for discussions, and Martijn van de Pol and Liam Bailey for useful advice about windows comparison. We thank all the people who helped maintain the study site and conduct the blue tit monitoring, in particular Jacques Blondel, Philippe Perret, Christophe de Franceschi, Annick Lucas, Pablo Giovannini, Samuel Caro, Claire Doutrelant and Marcel Lambrechts.

Funding

This project was funded by the CNRS and the European Research Council (Starting grant ERC-2013-StG-337365-SHE to A.C., and starting grant StG-678140-FluctEvol to L.M.C.).

Ethics

The blue tit monitoring protocol was approved by the Animal Care and Use committee Languedoc-Rousillon (CEEA-LR-12066) as well as by Regional Institutions (bylaw issued by the Prefecture on 15/06/2012 n° 2012167-0003). Captures of breeding birds were performed under ringing permits delivered by the CRBPO (Centre de Recherches sur la Biologie des Populations d'Oiseaux, Paris; ringing permit number 1907 to AC, program permit number 369).

Data accessibility

Data will be available from the Dryad Digital Repository after acceptance of the article. Meteorological data are accessible from https://donneespubliques.meteofrance.fr/ Author's contributions

Competing interests

We have no competing interests.