Laura Gervais
email: lauragervais51@gmail.com

Pedigree-free quantitative genetic approach provides evidence for heritability of movement tactics in wild roe deer

Keywords: (3-6 keys word): evolutionary potential, repeatability, personality variation, animal model, evolvability, genomic relatedness matrix

Assessing the evolutionary potential of animal populations in the wild is crucial to understanding how they may respond to selection mediated by rapid environmental change (e.g. habitat loss and fragmentation). A growing number of studies have investigated the adaptive role of behaviour, but assessments of its genetic basis in a natural setting remain scarce. We combined intensive biologging technology with genome-wide data and a pedigree-free quantitative genetic approach to quantify repeatability, heritability and evolvability for a suite of behaviours related to the risk avoidanceresource acquisition trade-off in a wild roe deer (Capreolus capreolus) population inhabiting a heterogeneous, human-dominated landscape. These traits, linked to the stress response, movement and space-use behaviour, were all moderately to highly repeatable. Furthermore, the repeatable among-individual component of variation in these traits was partly due to additive genetic variance, with heritability estimates ranging from 0.21±0.08 to 0.70±0.11 and evolvability ranging from 1.1 to 4.3%. Changes in the trait mean can therefore occur under hypothetical directional selection over just a few generations. To the best of our knowledge, this is the first empirical demonstration of additive genetic variation in space-use behaviour in a free-ranging population based on genomic relatedness data. We conclude that wild animal populations may have the potential to adjust their spatial behaviour to human-driven environmental modifications through micro-evolutionary change.

Introduction

Over evolutionary time, all organisms have evolved in relation to the environmental heterogeneity that is a ubiquitous feature in the wild. However, because recent human-induced environmental change is increasingly rapid, organisms are currently experiencing novel and rapidly fluctuating conditions to which they must adjust over very short timescales [START_REF] Parmesan | Ecological and Evolutionary Responses to Recent Climate Change[END_REF]. This adjustment occurs mainly through two non-exclusive mechanisms, by modification of the phenotype (i.e. microevolutionary change, phenotypic plasticity) or by modification of the environment (i.e. selection or adjustment of the environment) [START_REF] Edelaar | Appreciating the Multiple Processes Increasing Individual or Population Fitness[END_REF]. Behaviour plays an essential role in this adjustment since it is a key mediator of the relationship between an organism and its environment. Indeed, it is generally considered to be one of the most plastic components of the phenotype [START_REF] Tuomainen | Behavioural responses to human-induced environmental change[END_REF]. However, behavioural plasticity may not suffice to sustain a directional response, and micro-evolutionary change may be required for a population to adapt to a highly modified environment [START_REF] Both | Adjustment to climate change is constrained by arrival date in a long-distance migrant bird[END_REF][START_REF] Wong | Behavioral responses to changing environments[END_REF].

A genetic basis for behavioural variation among individuals is the raw material for adaptive evolution, assuming that the ability of certain behavioural phenotypes to cope with novel conditions provides them with a fitness benefit. Long considered as statistical noise [START_REF] Magurran | Adaptive individual differences within single populations[END_REF], behavioural variation among individuals is currently a central focus for behavioural ecologists. Many recent studies have highlighted the remarkable consistency of behavioural differences among individuals across time and ecological contexts in a wide range of taxa (also known as "personality traits" [START_REF] Réale | Integrating animal temperament within ecology and evolution[END_REF]. While this consistency (aka repeatability) does not preclude a degree of behavioural plasticity, it implies that an individual does not express the entire range of behavioural responses observed within the population. To predict how animals will behaviourally adjust to a novel environmental context, we need to evaluate the respective importance of personality and plasticity in observed behavioural variation and to quantify the extent to which personality variation is heritable. Adaptive explanations for the existence of consistent among-individual differences assume that phenotypic variation reliably reflects additive genetic variation (i.e. "the phenotypic gambit [START_REF] Grafen | Natural Selection, Kin Selection and Group Selection[END_REF]. Studies under controlled conditions have reported heritability values encompassing the widest range possible (from 0 to 100%), with a 30% average [START_REF] Mousseau | Natural selection and the heritability of fitness components[END_REF][START_REF] Stirling | Selection, structure and the heritability of behaviour[END_REF][START_REF] Van Oers | Quantitative and molecular genetics[END_REF]. [START_REF] Postma | Four decades of estimating heritabilities in wild vertebrate populations: Improved methods, more data, better estimates[END_REF] recently estimated an average heritability of 50% across wild populations. However, this value is probably inflated because additive genetic effects may be confounded with non-genetic permanent environmental effects in field studies [START_REF] Dingemanse | Individual behaviour: Behavioural ecology meets quantitative genetics[END_REF]. As a consequence, whether the phenotypic gambit generally holds for personality variations in the wild

Accepted Article

Version postprint

Comment citer ce document : Gervais, L. (Auteur de correspondance), Hewison, A. M., Morellet, N., Bernard, M., Merlet, J., Cargnelutti, B., Chaval, Y., Pujol, B., Quéméré, E. (Auteur de correspondance) (2020). Pedigreefree quantitative genetic approach provides evidence for heritability of movement tactics in wild roe deer. Journal of Evolutionary Biology. , DOI : 10.1111/jeb.13594

This article is protected by copyright. All rights reserved remains unclear. Hence, it is crucial to evaluate the genetic and environmental basis of amongindividual variation in behaviour which can only be reliably assessed when both heritability and repeatability are estimated simultaneously [START_REF] Dochtermann | The contribution of additive genetic variation to personality variation: Heritability of personality[END_REF].

Measuring behavioural repeatability in free ranging populations is challenging. As a consequence, most studies of wild animals index personality variation using common garden experiments [START_REF] Archard | The importance of wild populations in studies of animal temperament: Animal temperament in wild populations[END_REF], likely constraining an individual's behavioural repertoire. Individual specialisation in spatial behaviour underlies many personality traits, for example, activity or exploratory behaviour, both of which are linked to movement and space-use. Recent advances in biologging technology offer the opportunity to study the interaction between individual spatial behaviours and the environmental context in which these behaviours are expressed. However, the use of tracking technologies to assess facets of personality linked to movement has only recently been exploited [START_REF] Nilson | Individuality in movement: The role of animal personality[END_REF][START_REF] Spiegel | What's your move? Movement as a link between personality and spatial dynamics in animal populations[END_REF][START_REF] Spiegel | When the going gets tough: Behavioural type-dependent space use in the sleepy lizard changes as the season dries[END_REF][START_REF] Villegas-Ríos | Individual level consistency and correlations of fish spatial behaviour assessed from aquatic animal telemetry[END_REF][START_REF] Ward-Fear | The ecological and life history correlates of boldness in free-ranging lizards[END_REF]. Furthermore, while it is well established that some exceptional movements such as natal dispersal or migration are partly genetically based [START_REF] Charmantier | Heritability of short-scale natal dispersal in a large-scale foraging bird, the wandering albatross[END_REF][START_REF] Greenwood | The Role of Dispersal in the Great Tit (Parus major): The Causes, Consequences and Heritability of Natal Dispersal[END_REF][START_REF] Saastamoinen | Genetics of dispersal[END_REF][START_REF] Waser | Heritability of dispersal in banner-tailed kangaroo rats, Dipodomys spectabilis[END_REF], little is known about the genetic variation underpinning routine movement behaviours that drive most space-use behaviour. In order to evaluate their evolutionary potential, it is therefore necessary to use an integrative approach, combining the movement ecology and personality frameworks with quantitative genetics.

Measuring genetic variance and the heritability of a trait in wild populations requires estimates of pairwise relatedness amongst individuals [START_REF] Lynch | Genetics and analysis of quantitative traits[END_REF]. Expected relatedness estimates have traditionally been calculated on the basis of population-level multigenerational pedigrees, which are themselves based on field observations or on parentage assignment inferred from multilocus genotypes [START_REF] Wilson | An ecologist's guide to the animal model[END_REF]. Advances in high-throughput sequencing technologies mean that it is now possible to directly estimate realised pairwise genetic relatedness from genomic data (thousands of SNPs) in free-ranging organisms at a reasonable cost [START_REF] Gervais | RAD-sequencing for estimating GRM-based heritability in the wild: A case study in roe deer[END_REF][START_REF] Gienapp | Genomic Quantitative Genetics to Study Evolution in the Wild[END_REF]. In this study, we combined a pedigree-free quantitative genetic approach with high spatio-temporal resolution monitoring of individual behaviour in a large sample of free ranging roe deer (Capreolus capreolus) inhabiting a highly heterogeneous agricultural landscape. This landscape has been extensively modified by human activities (e.g. agricultural practice, infrastructures) such that high quality foraging resources occur in locations where human disturbance is also high. Hence, animals must trade-off acquisition of high-quality resources against risk avoidance (Bonnot et al.,

Materials & methods

Study area and data collection

This study is based on the long-term monitoring by capture-mark-recapture of a free-ranging Three hundred and thirty four roe deer were caught between 2002 and 2018 during annual winter capture operations using drive netting. When a roe deer was caught in the net, it was tranquilized with an intra-muscular injection of acepromazine (Calmivet-3cc) and placed in a wooden retention box designed to minimize stress and risk of injury. During marking procedures (ca. 10 mins/animal), each animal was sexed, aged as adult (>1 year old) or juvenile (<1 year old) in relation to tooth eruption [START_REF] Ratcliffe | Age determination of roe deer[END_REF]) and ear tagged. Blood and skin were sampled for subsequent analyses of the physiological stress response to capture and genetic information, respectively. The total time elapsed from capture in the net to release was recorded for each individual. Mean ambient temperature (°C) for each capture event was provided by the nearest Meteo-France weather station. All capture and marking procedures were in accordance with French and European law for animal welfare (prefectural order from the Toulouse Administrative Authority to capture and monitor wild roe deer and agreement no. A31113001 approved by the Departmental Authority of Population Protection). For more details on the capture procedure and ethical considerations, see [START_REF] Morellet | The Effect of Capture on Ranging Behaviour and Activity of the European Roe Deer Capreolus capreolus[END_REF] and Bonnot et al. (2018).

Physiological and behavioural responses to capture

Capture can be considered as a form of acute stress for wild animals which should provoke a suite of physiological (e.g. hyperthermia, change in blood parameters) and behavioural [START_REF] Moberg | The biology of animal stress: Basic principles and implications for animal welfare[END_REF][START_REF] Moe | Effects of Handling and Physical Restraint on Rectal Temperature, Cortisol, Glucose and Leucocyte[END_REF] This article is protected by copyright. All rights reserved the nearest 0.1°C during handling and recorded the average of the temperature at the beginning and at the end of the marking procedure as a measure of stress-induced hyperthermia. Individuals with abnormally low body temperature (T°C < 36.9) were removed from the dataset (n=5).

Haematological parameters were assessed on blood samples taken from the jugular vein (10 mL).

Differential white cell counts were performed by visual determination on the first 100 Wright-Giemsa-stained blood smears to quantify the proportion of lymphocytes and neutrophils in relation to the total white cell count. We calculated the ratio of neutrophils to lymphocytes (hereafter NL ratio) as an indicator of the chronic stress response [START_REF] Davis | The use of leukocyte profiles to measure stress in vertebrates: A review for ecologists[END_REF]. In a recent work on three roe deer populations experiencing fluctuating environmental conditions, [START_REF] Carbillet | The neutrophil to lymphocyte ratio indexes individual variation in the behavioural stress response of wild roe deer across fluctuating environmental conditions[END_REF] showed that this measure was consistently repeatable, irrespective of the capture protocol and provided an informative indicator of chronic stress that is partly related to environmental quality (e.g. food availability and quality, human disturbance). Finally, we recorded a behavioural score (hereafter docility) to index the individual's overall response to handling using a simple binary metric: we attributed a score of one if the individual struggled during manipulation (non-docile individuals), or a score of zero if the individual was passive (docile).

Space-use and mobility

Most of the studied deer (286 out of 334) were equipped with GPS collars to monitor individual space-use and movement. Collars were programmed to obtain a GPS fix every four hours (in 2003 and 2004) or every six hours (from 2005 to 2018) for approximately one year (baseline monitoring).

In addition, the collars were programmed to obtain a GPS fix at high frequency (1 location/10 min) over a 24h period approximately once per month (intensive monitoring). We performed differential correction to improve fix accuracy [START_REF] Adrados | GPS approach to study finescale site use by wild red deer during active and inactive behaviors[END_REF] and removed aberrant GPS fixes that would imply unfeasible movements (0.05% of all locations). We also removed all GPS fixes taken during the first 10 days after capture because of the potential alteration of spatial behaviour [START_REF] Morellet | The Effect of Capture on Ranging Behaviour and Activity of the European Roe Deer Capreolus capreolus[END_REF]. Furthermore, as one third of juveniles disperse [START_REF] Debeffe | Condition-dependent natal dispersal in a large herbivore: Heavier animals show a greater propensity to disperse and travel further[END_REF], we only retained locations recorded after the dispersal period (from the end of March to the end of May, [START_REF] Debeffe | Condition-dependent natal dispersal in a large herbivore: Heavier animals show a greater propensity to disperse and travel further[END_REF] for deer classified as dispersers (n=45) to ensure that we analysed routine home range behaviour only. Finally, we also removed those months that had less than 15 days (for baseline monitoring) or less than 20 hours (for intensive monitoring) of data for a given individual.

Based on aerial photographs of the study site (from the IGN's BD Ortho, http://professionnels.ign.fr/bdortho-50cm), we manually digitized homogeneous habitat polygons

Accepted Article

Version postprint

Comment citer ce document : Gervais, L. (Auteur de correspondance), Hewison, A. M., Morellet, N., Bernard, M., Merlet, J., Cargnelutti, B., Chaval, Y., Pujol, B., Quéméré, E. (Auteur de correspondance) (2020). Pedigreefree quantitative genetic approach provides evidence for heritability of movement tactics in wild roe deer. Journal of Evolutionary Biology. , DOI : 10.1111/jeb.13594

This article is protected by copyright. All rights reserved (in ArcView GIS 3.3, Esri, Redlands, CA, U.S.A.). Each polygon was then assigned to a habitat type (e.g. woodland, hedgerow, meadow or crop) that was determined annually by field observations during summer. Meadows and agricultural crops are numerous in the study system and provide abundant high quality resources [START_REF] Abbas | Landscape fragmentation generates spatial variation of diet composition and quality in a generalist herbivore[END_REF][START_REF] Hewison | Landscape fragmentation influences winter body mass of roe deer[END_REF], but are also more exposed to human disturbance and associated risks (e.g. road traffic accidents, agricultural practices, hunting) [START_REF] Coulon | Inferring the effects of landscape structure on roe deer (Capreolus capreolus) movements using a step selection function[END_REF][START_REF] Martin | Temporal shifts in landscape connectivity for an ecosystem engineer, the roe deer, across a multiple-use landscape[END_REF][START_REF] Padié | Roe deer at risk: Teasing apart habitat selection and landscape constraints in risk exposure at multiple scales[END_REF]. Roe deer are known to actively

avoid roads and open habitats during daytime, when human disturbance is high [START_REF] Bonnot | Habitat use under predation risk: Hunting, roads and human dwellings influence the spatial behaviour of roe deer[END_REF][START_REF] Martin | Temporal shifts in landscape connectivity for an ecosystem engineer, the roe deer, across a multiple-use landscape[END_REF][START_REF] Padié | Roe deer at risk: Teasing apart habitat selection and landscape constraints in risk exposure at multiple scales[END_REF].Thus, we assessed roe deer space-use behaviour in terms of risk taking as the probability of being in open habitat (i.e. meadows or crop fields) during daytime (hereafter P OH , [START_REF] Bonnot | Interindividual variability in habitat use: Evidence for a risk management syndrome in roe deer?[END_REF], and the proximity to roads.

P OH was calculated as the average proportion of fixes at 12:00 that occurred in open habitat, ranging between 0 for an individual that was always located in closed habitat during daytime, and 1 for an individual always located in open habitat during daytime. For each individual and each GPS fix at 12:00, we also measured the distance (km) to the nearest road (hereafter distance to roads). P OH and distance to roads were averaged per month in order to describe seasonal variation. Finally, we measured movement behaviour by calculating daily average speed (km.h -1 ) during the intensive monitoring periods as the total distance travelled during 24h (sum of step lengths) divided by the total interval duration (approx. 24h) (package adehabitatLT, [START_REF] Calenge | Analysis of animal movements in R: The adehabitatLT Package[END_REF], and averaged these values per month for each individual.

Genomic relatedness

The DNA of the 334 individuals was extracted from skin samples using DNeasy TM Tissue Kit (Qiagen)

and genotyped using a double-digest RAD sequencing approach [START_REF] Peterson | Double Digest RADseq: An Inexpensive Method for De Novo SNP Discovery and Genotyping in Model and Non-Model Species[END_REF] following the procedure detailed in [START_REF] Gervais | RAD-sequencing for estimating GRM-based heritability in the wild: A case study in roe deer[END_REF]. Briefly, individuals DNA was digested using two restriction enzymes, EcoR1 and MsPI, and size-selected libraries of 270-330 bp insert-size were sequenced on Illumina HiSeq 2500 (V4) (paired-end reads of 125 bp) so as to obtain an average depth of coverage of 76 reads/loci. De novo reconstruction of markers (i.e without a reference genome) and SNP calling was carried out with Stacks 1.35 software [START_REF] Catchen | Stacks: An analysis tool set for population genomics[END_REF]. After removing loci present in less than 80% of individuals and with Minor-Allele-Frequency (MAF) lower than 1%, we retained 14,887 polymorphic loci with a very low rate of genotyping error (<1%) and missing data (< 6 %) (see [START_REF] Gervais | RAD-sequencing for estimating GRM-based heritability in the wild: A case study in roe deer[END_REF] for more details). This final dataset was used to compute the Genomic Relatedness Matrix (GRM) of

Accepted Article

Comment citer ce document : Gervais, L. (Auteur de correspondance), Hewison, A. M., Morellet, N., Bernard, M., Merlet, J., Cargnelutti, B., Chaval, Y., Pujol, B., Quéméré, E. (Auteur de correspondance) (2020). Pedigreefree quantitative genetic approach provides evidence for heritability of movement tactics in wild roe deer. Journal of Evolutionary Biology. , DOI : 10.1111/jeb.13594

This article is protected by copyright. All rights reserved the population using identity by state (IBS) SNP relationships as implemented in the GCTA software tool [START_REF] Yang | GCTA: A Tool for Genome-wide Complex Trait Analysis[END_REF].

Quantitative genetic analyses

Our aim was to partition the phenotypic variance for fitness-related behavioural and physiological traits (Bonnot et al., 2018;[START_REF] Monestier | Is a proactive mum a good mum? A mother's coping style influences early fawn survival in roe deer[END_REF] into among-and within-individual variation, and then partition among-individual variation into its genetic and non-genetic (permanent environmental effects) components. We used Linear-Mixed-effect Models (LMMs) where parameters were calculated using the REstricted Maximum Likelihood (REML) approach [START_REF] Lynch | Genetics and analysis of quantitative traits[END_REF]. This approach, also called the "animal model" in quantitative genetics, partitions the phenotypic variance (V P ) into additive genetic variance (V A ), permanent environmental variance (V PE )

and residual variance (V R ), so that V P= V A +V PE +V R conditional on any fixed effects. We included ageclass and sex as fixed effects in all models because behavioural traits may differ between life history stages and sexes [START_REF] Debeffe | Short-and long-term repeatability of docility in the roe deer: Sex and age matter[END_REF]. We also included a year effect to account for inter-annual differences in individual behaviour and environmental conditions in all models. In addition, in models investigating variation in response to the capture event (body temperature, NL ratio and docility), we included as fixed effects (i) the cumulative number of capture events for each individual, because repeated capture might lead to habituation and an attenuation of the behavioural response [START_REF] Bonnot | Who's afraid of the big bad wolf? Variation in the stress response among personalities and populations in a large wild herbivore[END_REF]) and (ii) the time between the actual capture event and release (i.e. time delay) to control for potential fluctuations in the stress response, assuming a linear relationship with time elapsed [START_REF] Montané | Effect of acepromazine on the signs of capture stress in captive and free-ranging roe deer (Capreolus capreolus)[END_REF]. In the few cases (n=14) where the time delay was not available, we used the population average time delay (248 min). For the body temperature model only, we additionally included the mean ambient temperature (°C) during capture. In models investigating variation in space-use and movement behaviour (P OH , distance to roads and daily average speed), we additionally included month (1-12)

and sex and their two-way interaction as fixed effects to account for potential divergence in temporal variation between the sexes. Finally, as recommended by [START_REF] Wilson | Why h2 does not always equal VA/VP?: Why h2 does not always equal VA/VP?[END_REF], all models were also run with "methodological bias" fixed-effects only (thus, excluding all fixed-effects related to biological or temporal variation i.e. age-class, sex, year of capture and month of observation) (Table S2).

We This article is protected by copyright. All rights reserved of traits was estimated as ind²=V IND /(V IND +V R ) [START_REF] Nakagawa | Repeatability for Gaussian and non-Gaussian data: A practical guide for biologists[END_REF]. As suggested by [START_REF] Dochtermann | The mean matters: Going beyond repeatability to interpret behavioural variation[END_REF], we also reported a mean-squared standardised metric as

𝐼 𝐼 = 𝑉 𝐼𝑁𝐷 𝑋 2 ⁄
, where 𝑋 corresponds to the trait mean, to allow for direct comparison of the magnitude of the among-individual variation across traits (on a ratio or log scales, see below for further details). Second, we estimated additive genetic and non-genetic variation by decomposing among-individual variation into additive genetic variance (V A ) and permanent environmental variance (V PE ), including relatedness between pairs of individuals (GRM) as a random effect on the intercept so that

V IND =V A +V PE with V P= V IND +V R .
Note that the GRM was non-positive definite but could nevertheless be implemented in an REML animal model (see [START_REF] Gilmour | Handling Non Positive Definite relationship matrices in Mixed Models[END_REF][START_REF] Gilmour | Average information residual maximum likelihood in practice[END_REF][START_REF] Meyer | Bending" and beyond: Better estimates of quantitative genetic parameters[END_REF] for implementation of none-positive definite matrix).

Narrow sense heritability of the traits was estimated as h²=V A /(V A +V PE +V R ) and the permanent environmental effect as pe²=V PE /(V A +V PE +V R ), while the within-individual effect was approximated by the residual effect calculated as r²=V R /(V A +V PE +V R ) [START_REF] Wilson | An ecologist's guide to the animal model[END_REF]. Evolvability of the traits was estimated using the additive genetic variance standardized by the square mean of the trait as

𝐼 𝐴 = 𝑉 𝐴 𝑋 2 ⁄
where 𝑋 is the trait mean before selection [START_REF] Houle | Comparing evolvability and variability of quantitative traits[END_REF]. This can be interpreted as the expected rate of trait change if it were submitted to a selection pressure which was as strong as on fitness itself ("a unit strength of selection", [START_REF] Houle | Comparing evolvability and variability of quantitative traits[END_REF]. Evolvability is interpretable for traits on a log-interval or ratio measurement scale (Garcia-Gonzalez, Simmons, Tomkins, Kotiaho, & Evans, 2012), but is not appropriate for traits without a natural zero point, proportions, or traits with a mean of zero [START_REF] Hansen | Heritability is not Evolvability[END_REF]. As a consequence, we do not report evolvability estimates for body temperature, docility or P OH . For log-transformed data (NL ratio, distance to roads, daily average speed, see below), we directly used the additive genetic variance as an approximation of I A on a raw scale (Garcia-Gonzalez, Simmons, Tomkins, Kotiaho, & Evans, 2012; [START_REF] Hansen | Heritability is not Evolvability[END_REF] .

All analyses were run on ASREML-R v4.00 [START_REF] Butler | ASReml-R reference manual (version 3). The State of Queensland[END_REF]. After inspecting the distributions of residuals for the models, we decided to log-transform traits with continuous distributions to meet assumptions of homoscedasticity and normality of residuals (with the exception of body temperature and P OH that already had normally distributed residuals and distance to roads that was log (x+1) transformed). We analysed the binary trait (docility) as if it was Gaussian, assuming that a discontinuous phenotype is determined by an underlying continuous distribution with a threshold above which the behaviour is considered to be expressed and below which it is not [START_REF] Elston | Query: Estimating" Heritability" of a dichotomous trait[END_REF][START_REF] Roff | The evolution of threshold traits: A quantitative genetic analysis of the physiological and life-history correlates of wing dimorphism in the sand cricket[END_REF][START_REF] Van Vleck | Estimation of heritability of threshold characters[END_REF]. We report the heritability estimate for docility on the true data-scale (see [START_REF] Dempster | Heritability of threshold characters[END_REF]. Note that

Accepted Article

Version This article is protected by copyright. All rights reserved REML animal models with a binary data distribution family are known to underestimate all the variance components [START_REF] Butler | ASReml-R reference manual (version 3). The State of Queensland[END_REF] and that estimates from threshold models must also be interpreted with caution [START_REF] De Villemereuil | Quantitative genetic methods depending on the nature of the phenotypic trait[END_REF][START_REF] De Villemereuil | General Methods for Evolutionary Quantitative Genetic Inference from Generalized Mixed Models[END_REF]. Finally, the statistical significance of variance components was assessed using Likelihood Ratio Tests (LRT) [START_REF] Pinheiro | Mixed-effects models in S and S-PLUS[END_REF], assuming that twice the difference in log-likelihood between nested models follows a chi-squared distribution with one degree of freedom.

Results

Data summary

The number of times that each of the 334 individuals was captured ranged from one to six, with a mean of 1.93 captures per individual. Of the 286 individuals equipped with GPS collars, 10 individuals were removed due to an insufficient number of days monitoring. The mean number of monitoring days per individual was 330 (ranging from 57 to 1,268). See Table S1 for more details on sample size, means and phenotypic variance (V P ) of traits. The Genomic Relatedness Matrix included 55,945 pairwise relatedness coefficients, of which 825 pairs had a relatedness coefficient higher than 0.1. The majority of individuals were unrelated (see Figure S1). The total variance of the relatedness coefficient off-diagonal in the population was 0.014.

Accepted Article

Version This article is protected by copyright. All rights reserved

Partitioning trait variation into between-and within-individual components

We found moderate to high levels of consistent among-individual variation (repeatability) for all the behavioural and physiological traits studied (Figure 1, Table 1). First, traits relating to stress management yielded similar levels of repeatability (r²±SE), ranging from 0.26±0.10 for docility, to 0.31±0.12 for the NL ratio, and 0.38±0.11 for body temperature (Figure 1). Second, movement behaviour (daily average speed) was similarly repeatable (r²=0.29±0.03). Finally, individuals showed a high level of consistency in space-use behaviour, with a repeatability of 0.62±0.02 for P OH and 0.76±0.02 for distance to roads (Figure 1). However, although highly repeatable, distance to roads had the lowest value for mean-squared standardised among-individual variation (I I =1.2±0.1%). In contrast, the NL ratio was moderately repeatable, but had the highest level of among-individual variation (I I =8.3±3.3%). Therefore, among individual variation was more marked for the NL ratio than for distance to roads.

Estimating additive genetic and permanent environmental variance

We found statistically significant additive genetic variance for traits related to space-use and movement, but not for traits related to stress management (body temperature, NL ratio, docility)

(Table 1). Narrow sense heritability (h²±SE) differed markedly among traits (Figure 1), ranging from 0.16±0.15 for the NL ratio to 0.70±0.11 for distance to roads, while the permanent environmental effects (pe²) ranged from 0.05±0.10 for distance to roads to 0.20±0.09 for P OH . Evolvability (I A ±SE%) also differed markedly among traits, ranging from 1.1±0.2% to 4.3±4.2% for distance to roads and NL ratio, respectively. The level of among-individual variation in the three stress management behaviours was similarly moderate. However, the models were unable to tease apart the genetic and permanent environmental components for these three traits, despite non-zero additive genetic variance (for NL ratio h²=0.16±0.15, χ²(1,N=301)=1.018, P=0.313, for docility h²=0.17±0.12, χ²(1,N=334) =3.442, P=0.063 and for body temperature h²=0.32±0.16, χ²(1,N=329)=3.525, P=0.060).

In contrast, we were able to estimate both additive genetic and permanent environmental variances for daily average speed, despite a moderate level of repeatability (r²=29%). Lastly, among-individual variation for traits related to movement and space-use behaviours (daily average speed, P OH and distance to roads) was mostly attributed to additive genetic variance (h²=0.21±0.08, 0.41±0.10 and 0.70±0.11, respectively) (Figure 1, Table 1).

Accepted Article

Version postprint

Comment citer ce document : Gervais, L. (Auteur de correspondance), Hewison, A. M., Morellet, N., Bernard, M., Merlet, J., Cargnelutti, B., Chaval, Y., Pujol, B., Quéméré, E. (Auteur de correspondance) (2020). Pedigreefree quantitative genetic approach provides evidence for heritability of movement tactics in wild roe deer. Journal of Evolutionary Biology. , DOI : 10.1111/jeb.13594

This article is protected by copyright. All rights reserved Table 1. Variance component estimates and their associated ratios for stress management, space-use and movement behaviours. Components of phenotypic variance (V PE : permanent environmental variance, V A : additive genetic variance, V R : residual variance) and their associated ratios (ind²: repeatability, pe²: permanent environmental effects, h²: heritability, r²: residual effect, I I : mean-squared standardised among-individual variation and I A : evolvability) are reported with their standard errors in brackets.

V IND : among-individual variance is the sum of V PE and V A . Note that repeatability estimates and associated statistics were derived before the inclusion of the genetic effect.

Significance of ind², h², pe² was assessed with a LRT test (Chi-square statistic (χ²) and p-value (P) with one degree of freedom). I A values are only reported for variables on a log-interval or ratio scale.

Variables ind² V PE V A V R pe²

Discussion

There is a growing interest in how environmental heterogeneity and individual consistency can generate specialisation in movement tactics. However, an understanding of the causes of movement specialisation and its eco-evolutionary consequences is still lacking [START_REF] Spiegel | What's your move? Movement as a link between personality and spatial dynamics in animal populations[END_REF]. Our study provides, as far as we know, the first evidence of substantial heritability and evolvability in space-use behaviour in the wild. It also illustrates the potential of genomic-based relatedness approaches to provide evidence for the evolutionary potential of formerly unexplored fitness-related traits in nonmodel species. In a wild large herbivore population inhabiting a spatially heterogeneous landscape, we found significant evidence of additive genetic variation for a suite of behavioural traits which are directly implicated in the management of the risk avoidance-resource acquisition trade-off and are related to fitness components (Bonnot et al., 2018;[START_REF] Monestier | Is a proactive mum a good mum? A mother's coping style influences early fawn survival in roe deer[END_REF]. Below we discuss the evolutionary implications of this in terms of how roe deer populations may adapt to global changes, such as landscape fragmentation, and highlight potential caveats in assessing the micro-evolutionary potential of space-use and movement behaviour in a heterogeneous landscape.

Presence of personality variation in stress-management, space-use and movement behaviour

Our results revealed the presence of substantial among-individual variation in all the studied behavioural and physiological traits. Average repeatability across traits (ca. 30%) was consistent with mean values observed by [START_REF] Bell | The repeatability of behaviour: A meta-analysis[END_REF][START_REF] Holtmann | Metabolic rates, and not hormone levels, are a likely mediator of between-individual differences in behaviour: A meta-analysis[END_REF]. However, repeatability differed markedly among behavioural traits. In particular, space-use behaviours were twice as consistent at the individual level as the other traits studied. This is not surprising, as roe deer are known to be sedentary and exhibit high spatial fidelity [START_REF] Hewison | Social organisation of European roe deer. The European Roe Deer: The Biology of Success[END_REF]. Studies reporting the presence of consistent among-individual variation in space-use and movement behaviour are on the increase [START_REF] Campos-Candela | A mechanistic theory of personalitydependent movement behaviour based on dynamic energy budgets[END_REF][START_REF] Hertel | Don't poke the bear: Using tracking data to quantify behavioural syndromes in elusive wildlife[END_REF][START_REF] Schirmer | Individuals in space: Personality-dependent space use, movement and microhabitat use facilitate individual spatial niche specialization[END_REF][START_REF] Spiegel | What's your move? Movement as a link between personality and spatial dynamics in animal populations[END_REF]. While these findings have implications for our understanding of how individuals may respond to habitat fragmentation [START_REF] Spiegel | What's your move? Movement as a link between personality and spatial dynamics in animal populations[END_REF], little is currently known about the mechanisms underpinning individual-level consistency in space-use behaviour.

To facilitate comparisons of traits both within and across studies, it is important to estimate meanrather than variance-standardised among-individual variation for each trait [START_REF] Dochtermann | The mean matters: Going beyond repeatability to interpret behavioural variation[END_REF]. Indeed, repeatability estimates are, by definition, strongly dependent on the total phenotypic variance for a given trait which, in the specific case of space-use behaviours, is partly related to inherent differences in habitat structure. Consequently, although space-use behaviours were twice as repeatable compared to other non-spatial traits, mean-squared standardised among-

Accepted Article

Version postprint

Comment citer ce document : Gervais, L. (Auteur de correspondance), Hewison, A. M., Morellet, N., Bernard, M., Merlet, J., Cargnelutti, B., Chaval, Y., Pujol, B., Quéméré, E. (Auteur de correspondance) (2020). Pedigreefree quantitative genetic approach provides evidence for heritability of movement tactics in wild roe deer. Journal of Evolutionary Biology. , DOI : 10.1111/jeb.13594

This article is protected by copyright. All rights reserved individual variation was much lower. This indicates that the higher repeatability of space-use behaviours is not due to more pronounced differences among individuals, but rather to a low level of within-individual variation. There is a clear lack of comparable estimates for mean-squared standardised among-individual variation in behavioural traits from other wild systems. More studies are necessary if we want to understand how and why the magnitude of among-individual differences varies across traits and populations [START_REF] Dochtermann | The mean matters: Going beyond repeatability to interpret behavioural variation[END_REF][START_REF] Holtmann | Metabolic rates, and not hormone levels, are a likely mediator of between-individual differences in behaviour: A meta-analysis[END_REF].

Trait heritabilities, evolvabilities, and evolutionary implications

Understanding the relationship between heritability and repeatability is of critical interest in behavioural ecology [START_REF] Boake | Repeatability: Its role in evolutionary studies of mating behavior[END_REF][START_REF] Dochtermann | The contribution of additive genetic variation to personality variation: Heritability of personality[END_REF]. Indeed, adaptive explanations regarding the origin and maintenance of personality variation are widely assumed when using a phenotypic approach, but most often without investigating the mechanisms underpinning observed repeatability (genetic or environmental) [START_REF] Dingemanse | Recent models for adaptive personality differences: A review[END_REF]. Although it is now widely acknowledged that personality variation has, at least, some genetic basis [START_REF] Dochtermann | The contribution of additive genetic variation to personality variation: Heritability of personality[END_REF], there is still a lack of understanding of how the relative importance of genetic and non-genetic sources of variance might differ among personality traits [START_REF] Réale | Integrating animal temperament within ecology and evolution[END_REF]. In a recent metaanalysis, [START_REF] Dochtermann | The mean matters: Going beyond repeatability to interpret behavioural variation[END_REF] showed that heritability estimates varied from 14 to 45% among behavioural categories (e.g. activity, foraging, mating etc. behaviours). Similarly, we found that heritability estimates varied strongly among behavioural traits ranging from 21±8% for daily average speed to 70±11% for distance to roads. Comparing heritability among traits may lead to erroneous conclusions as it does not solely depend on the absolute magnitude of additive genetic variation.

Houle (1992) therefore suggested using mean-standardize additive genetic variation that is a better measure of trait evolvability. Studies reporting heritability or evolvability of behaviour in wild populations are still scarce in comparison to other types of traits. Our results showed that behavioural traits exhibiting significant heritability also had non-zero evolvability. Therefore, consistent directional selection on those traits should lead to marked changes in the trait mean over a few hundred generations (1.1%-4.6% per generation per unit of selection). This is within the range of expected change reported for behavioural traits by Postma (2014) (3.99±3.76%). However, the realization of this potential for evolutionary change is likely to be limited by constraining biological mechanisms [START_REF] Pujol | The Missing Response to Selection in the Wild[END_REF]. For example, genetic correlations might act as constraints to evolution and it is now acknowledged that personality traits may form a syndrome based on genetic correlations that tie any change in one trait to the changes occurring in other traits [START_REF] Blumstein | Heritability of anti-predatory traits: Vigilance and locomotor performance in marmots[END_REF][START_REF] Careau | Genetic correlation between resting metabolic rate and exploratory behaviour in deer mice (Peromyscus maniculatus)[END_REF][START_REF] Petelle | Heritability and genetic correlations of personality traits in a wild population of yellow-bellied marmots (Marmota flaviventris)[END_REF][START_REF] Taylor | Low heritabilities, but genetic and maternal correlations between red squirrel behaviours[END_REF].

Thus, further studies will be required to understand whether this mechanism may affect the

Accepted Article

Version This article is protected by copyright. All rights reserved evaluation of evolutionary potential in our study system. Moreover, selection pressures are expected to vary in space in this human-dominated landscape. For example, hunting intensity is variable across the landscape and may impose a selection pressure on behaviour [START_REF] Leclerc | Harvesting as a potential selective pressure on behavioural traits[END_REF], with some personality types more prone to be killed than others [START_REF] Ciuti | Human selection of elk behavioural traits in a landscape of fear[END_REF][START_REF] Lone | An adaptive behavioural response to hunting: Surviving male red deer shift habitat at the onset of the hunting season[END_REF]. This may constrain the change in the trait mean distribution, as fluctuating selection might also influence the evolutionary response [START_REF] Kruuk | New answers for old questions: the evolutionary quantitative genetics of wild animal populations[END_REF][START_REF] Pujol | The Missing Response to Selection in the Wild[END_REF]. To test this hypothesis, future work will investigate the predominant selection pressures and the responses of the roe deer in our study system.

The limits of assessing microevolutionary potential in heterogeneous landscapes

One of the novelties of our study is that we investigated the evolutionary potential of space-use and movement behaviour in a heterogeneous landscape. This is challenging because the observed heritability of a given trait may partly be due to environmental similarity when related individuals share close and/or similar habitats [START_REF] Kruuk | How to separate genetic and environmental causes of similarity between relatives: Separating genetic and environmental effects[END_REF][START_REF] Thomson | A guide to using a multiple-matrix animal model to disentangle genetic and nongenetic causes of phenotypic variance[END_REF]. Such environmental similarity effects were observed in a wild red deer population (Cervus elaphus, [START_REF] Stopher | Shared Spatial Effects on Quantitative Genetic Parameters: Accounting for Spatial Autocorrelation and Home Range Overlap Reduces Estimates of Heritability in Wild Red Deer[END_REF], but were absent in Soay sheep (Ovis aries, [START_REF] Regan | Accounting for female space sharing in St. Kilda Soay sheep (Ovis aries) results in little change in heritability estimates[END_REF] and song sparrows (Melospiza melodia, [START_REF] Germain | Direct and indirect genetic and finescale location effects on breeding date in song sparrows[END_REF]. Whether this might affect the heritability estimates that we report for our study population remains to be explored, but cannot be ignored. Relatives are likely to share similar environments because the majority of juveniles are philopatric in our study population (~60%, [START_REF] Debeffe | Condition-dependent natal dispersal in a large herbivore: Heavier animals show a greater propensity to disperse and travel further[END_REF].

Furthermore, although dispersers may travel relatively long distances (mean of 12.3±10.5 km [START_REF] Debeffe | Condition-dependent natal dispersal in a large herbivore: Heavier animals show a greater propensity to disperse and travel further[END_REF] in our study site) compared to the size of a home range (mean of 0.76 km² [START_REF] Morellet | Seasonality, weather and climate affect home range size in roe deer across a wide latitudinal gradient within Europe[END_REF], they may nevertheless select a settlement range which is as similar as possible to their natal range (Natal Habitat Preference Induction, [START_REF] Davis | The effect of natal experience on habitat preferences[END_REF]. In the absence of additional information to account for maternal or space-sharing effects, caution must be taken when interpreting additive genetic variance for habitat-dependent traits because estimates might be thus inflated (see [START_REF] Andrews | A bioinformatic pipeline for identifying informative SNP panels for parentage assignment from RADseq data[END_REF][START_REF] Huisman | Pedigree reconstruction from SNP data: Parentage assignment, sibship clustering and beyond[END_REF][START_REF] Thomson | A guide to using a multiple-matrix animal model to disentangle genetic and nongenetic causes of phenotypic variance[END_REF].

Accounting for environmental similarity between individuals implies that habitat use has no genetic basis and influences behaviour in an unidirectional way. However, causal relationships may be hard to disentangle under natural conditions if habitat use and behaviours share a confounded genetic and environmental basis. Home-range establishment during early life potentially covaries with personality-dependent spatial behaviour [START_REF] Spiegel | What's your move? Movement as a link between personality and spatial dynamics in animal populations[END_REF][START_REF] Villegas-Ríos | Personalities influence spatial responses to environmental fluctuations in wild fish[END_REF]. Thus, different personality types might consistently differ in the environmental stimuli that they experience [START_REF] Spiegel | What's your move? Movement as a link between personality and spatial dynamics in animal populations[END_REF], which may generate among-individual variation in

Accepted Article

Version postprint

Comment citer ce document : Gervais, L. (Auteur de correspondance), Hewison, A. M., Morellet, N., Bernard, M., Merlet, J., Cargnelutti, B., Chaval, Y., Pujol, B., Quéméré, E. (Auteur de correspondance) (2020). Pedigreefree quantitative genetic approach provides evidence for heritability of movement tactics in wild roe deer. Journal of Evolutionary Biology. , DOI : 10.1111/jeb.13594

This article is protected by copyright. All rights reserved behaviour [START_REF] Bonnot | Habitat use under predation risk: Hunting, roads and human dwellings influence the spatial behaviour of roe deer[END_REF][START_REF] Bonnot | Interindividual variability in habitat use: Evidence for a risk management syndrome in roe deer?[END_REF][START_REF] Coulon | Inferring the effects of landscape structure on roe deer (Capreolus capreolus) movements using a step selection function[END_REF][START_REF] Padié | Roe deer at risk: Teasing apart habitat selection and landscape constraints in risk exposure at multiple scales[END_REF]. This might be the case in roe deer because they use sub-seasonal home-ranges to track spatio-temporal variation in risk and resources [START_REF] Couriot | Truly sedentary? The multi-range tactic as a response to resource heterogeneity and unpredictability in a large herbivore[END_REF] and actively avoid exposure to human disturbance (e.g. roads, dwellings, agricultural crops) through adjustment of their circadian rhythm [START_REF] Bonnot | Habitat use under predation risk: Hunting, roads and human dwellings influence the spatial behaviour of roe deer[END_REF][START_REF] Coulon | Inferring the effects of landscape structure on roe deer (Capreolus capreolus) movements using a step selection function[END_REF][START_REF] Martin | Temporal shifts in landscape connectivity for an ecosystem engineer, the roe deer, across a multiple-use landscape[END_REF][START_REF] Padié | Roe deer at risk: Teasing apart habitat selection and landscape constraints in risk exposure at multiple scales[END_REF]. In addition, because roe deer are highly sedentary [START_REF] Hewison | Social organisation of European roe deer. The European Roe Deer: The Biology of Success[END_REF], habitat use during adulthood is strongly constrained by habitat selection during post-dispersal settlement. Hence, the habitat selection decision an individual takes at around one year of age [START_REF] Debeffe | Condition-dependent natal dispersal in a large herbivore: Heavier animals show a greater propensity to disperse and travel further[END_REF] will have long-term consequences on the behavioural tactic it employs to trade-off risk avoidance and resource acquisition.

Behavioural tactics employed during the juvenile stage may also influence habitat selection decisions during post-dispersal settlement and, more generally, behaviours during adulthood. In this situation, experimental approaches in semi-natural systems might be necessary to assess how personality variation is expressed when individuals are raised under identical or contrasting environmental conditions. This kind of approach might help clarify whether environmental drivers shape personality variation or personalities select particular environments.

Genomic relatedness and spatial behaviour in wild non-model organisms: promises and caveats

Genomic tools offer the opportunity for behavioural ecologists to study non-model organisms in natural populations [START_REF] Bengston | Genomic tools for behavioural ecologists to understand repeatable individual differences in behaviour[END_REF][START_REF] Gienapp | Genomic Quantitative Genetics to Study Evolution in the Wild[END_REF]. GRM-based approaches have already been shown to have suitable power for estimating trait heritabilities in free ranging populations [START_REF] Bérénos | Estimating quantitative genetic parameters in wild populations: A comparison of pedigree and genomic approaches[END_REF][START_REF] Perrier | Heritability estimates from genome wide relatedness matrices in wild populations: Application to a passerine, using a small sample size[END_REF].

Our study demonstrated that a GRM-based approach can be used with a relatively small sample size (e.g. a few hundred individuals). However, the critical points to consider appear to be variance in relatedness and the number of individual observations required to ensure enough statistical power to tease apart the genetic and non-genetic components of trait variation. Variance in relatedness in our data (1.4 × 10 -3 ) was within the range reported in Soay sheep (1.3 × 10 -3 ) [START_REF] Bérénos | Estimating quantitative genetic parameters in wild populations: A comparison of pedigree and genomic approaches[END_REF] and blue tits (4 × 10 -3 ) (C. Perrier, pers. comm.), which is 100 times higher than the predicted values for unrelated humans (Vinkhuyzen et al., 2013, 2 x 10 -5 ).

The use of repeated measures is known to provide better accuracy when decomposing phenotypic variation [START_REF] Martin | Measuring individual differences in reaction norms in field and experimental studies: a power analysis of random regression models: Power analysis of random regression models[END_REF]Pol, 2012) and also helps to distinguish permanent environmental effects from residual effects [START_REF] Wilson | An ecologist's guide to the animal model[END_REF]. Our results indirectly support the idea that the use of several observations per individual greatly improves statistical power to detect heritability [START_REF] Åkesson | Estimating Heritabilities and Genetic Correlations: Comparing the 'Animal Model' with Parent-Offspring Regression Using Data from a Natural Population[END_REF]. Indeed, we were able to detect a This article is protected by copyright. All rights reserved heritability of 0.21 (±0.08) using 205 individuals and 1856 observations for daily average speed, while we were unable to detect significant heritability for docility (0.17±0.12) using 249 individuals and 334 observations, despite comparable levels of repeatability. Although this conclusion is not restricted to the GRM-approach, it is likely an important consideration for GRM-based heritability assessments based on small sample sizes that often lack adequate statistical power. Yet, to the best of our knowledge, no theoretical studies have investigated the statistical power of GRM-based approaches as a function of the number of individual observations, although similar issues have emerged in Genome Wide Association Studies [START_REF] Rönnegård | Increasing the power of genome wide association studies in natural populations using repeated measuresevaluation and implementation[END_REF][START_REF] Silva | Insights into the genetic architecture of morphological traits in two passerine bird species[END_REF]. Biologging techniques that provide multiple observations of a number of behaviours across time therefore appear as a powerful tool to overcome these limitations.

Accepted Article

Conclusion

By combining intensive biologging technology and genome-wide data, we provided an empirical demonstration of substantial heritability and evolvability in movement and space-use tactics related to the management of the risk/resource trade-off. In a recent review, [START_REF] Tucker | Moving in the Anthropocene: Global reductions in terrestrial mammalian movements[END_REF] showed that human activities have profound consequences on wild animal movements, but were unable to address the evolutionary consequences for their long-term sustainability. Our findings imply that wild populations have the potential to respond to selection pressure generated by human activity through micro-evolutionary change in spatial behaviours. The lack of comparable data from other wild systems highlights the need for more studies on the heritability and evolvability of stress management, space-use and movement behaviour under natural conditions. We conclude that the synergy between data-rich technologies that monitor the spatial locations of individuals and genome wide analysis paves the way for such advances.

Data accessibility:

Genomic Relatedness Matrix, spatial behaviour data and associated R-script are available from the Dryad Digital Repository: https:// doi:10.5061/dryad.dbrv15dxg

ACKNOWLEDGMENTS:

We are grateful to the many co-workers, students and volunteers who helped during field data collection, in particular Jean-Marc Angibault, Jean-Luc Rames, Bruno Lourtet, Denis Picot, Nicolas

Cebe and Hélène Verheyden. We also thank Edinburgh genomics for RADseq data and the Genotoul Bioinfromatics Plateform Toulouse Midi-Pyrénées for providing help, computing and storage of RADseq data. We acknowledge John Hunt, Pim Edelaar and one anonymous reviewer for their comments that helped to improve the quality of this manuscript. We would also like to warmly

Accepted Article

Version Accepted Article

  decomposed phenotypic variance conditioned on fixed effects, into within-(V R ) and amongindividual variance (V IND ), including individual identity as a random effect on the intercept. V R represents within-individual variance due to plasticity and/or measurement error. The repeatability Accepted Article Comment citer ce document : Gervais, L. (Auteur de correspondance), Hewison, A. M., Morellet, N., Bernard, M., Merlet, J., Cargnelutti, B., Chaval, Y., Pujol, B., Quéméré, E. (Auteur de correspondance) (2020). Pedigreefree quantitative genetic approach provides evidence for heritability of movement tactics in wild roe deer. Journal of Evolutionary Biology. , DOI : 10.1111/jeb.13594

  postprint Comment citer ce document : Gervais, L. (Auteur de correspondance), Hewison, A. M., Morellet, N., Bernard, M., Merlet, J., Cargnelutti, B., Chaval, Y., Pujol, B., Quéméré, E. (Auteur de correspondance) (2020). Pedigreefree quantitative genetic approach provides evidence for heritability of movement tactics in wild roe deer. Journal of Evolutionary Biology. , DOI : 10.1111/jeb.13594

  postprint Comment citer ce document : Gervais, L. (Auteur de correspondance), Hewison, A. M., Morellet, N., Bernard, M., Merlet, J., Cargnelutti, B., Chaval, Y., Pujol, B., Quéméré, E. (Auteur de correspondance) (2020). Pedigreefree quantitative genetic approach provides evidence for heritability of movement tactics in wild roe deer. Journal of Evolutionary Biology. , DOI : 10.1111/jeb.13594

  postprint Comment citer ce document : Gervais, L. (Auteur de correspondance), Hewison, A. M., Morellet, N., Bernard, M., Merlet, J., Cargnelutti, B., Chaval, Y., Pujol, B., Quéméré, E. (Auteur de correspondance) (2020). Pedigreefree quantitative genetic approach provides evidence for heritability of movement tactics in wild roe deer. Journal of Evolutionary Biology. , DOI : 10.1111/jeb.13594

  postprint Comment citer ce document : Gervais, L. (Auteur de correspondance), Hewison, A. M., Morellet, N., Bernard, M., Merlet, J., Cargnelutti, B., Chaval, Y., Pujol, B., Quéméré, E. (Auteur de correspondance) (2020). Pedigreefree quantitative genetic approach provides evidence for heritability of movement tactics in wild roe deer. Journal of Evolutionary Biology. , DOI : 10.1111/jeb.13594 This article is protected by copyright. All rights reserved thank Ingrid David, Caroline Thomson and Isabel Winney for their advice on statistical analysis and Pierre de Villemereuil and Denis Réale for many helpful discussions on the project. This project was financially supported by the French National Institute for Agricultural Research (INRA), the UK Natural Environment Research Council (NERC), the Biomolecular Analysis Facility at the University of Edinburgh (NBAF901), the "PATCH" RPDOC ANR project (ANR-12-PDOC-0017-01) to Cécile Vanpé, and the "Mov-It" ANR project (ANR-16-CE02-0010-02) to Mark Hewison and Nicolas Morellet, the Centre National de la Recherche Scientifique (CNRS) core funding to Benoit Pujol, the EU within the framework of the Marie-Curie FP7 COFUND People Programme through the award of an Agreenskills fellowship to Erwan Quéméré, and the Federal University of Toulouse through the PhD scholarship of Laura Gervais (IDEX grant "GENEMOV" to Benoit Pujol and Erwan Quéméré).

Figures:Figure 1 :

 1 Figures:

  responses. Since 2009, we continuously measured rectal temperature to

	Accepted Article Version postprint
	Comment citer ce document :
	Gervais, L. (Auteur de correspondance), Hewison, A. M., Morellet, N., Bernard, M., Merlet, J.,
	Cargnelutti, B., Chaval, Y., Pujol, B., Quéméré, E. (Auteur de correspondance) (2020).
	Pedigreefree quantitative genetic approach provides evidence for heritability of movement tactics in
	wild roe deer. Journal of Evolutionary Biology. , DOI : 10.1111/jeb.13594