Natal habitat imprinting counteracts the diversifying effects of phenotype-dependent dispersal in a spatially structured population
© The Author(s). 2016
Received: 24 May 2016
Accepted: 18 July 2016
Published: 8 August 2016
Habitat selection may have profound evolutionary consequences, but they strongly depend on the underlying preference mechanism, including genetically-determined, natal habitat and phenotype-dependent preferences. It is known that different mechanisms may operate at the same time, yet their relative contribution to population differentiation remains largely unexplored empirically mainly because of the difficulty of finding suitable study systems. Here, we investigate the role of early experience and genetic background in determining the outcome of settlement by pied flycatchers (Ficedula hypoleuca) breeding in two habitat patches between which dispersal and subsequent reproductive performance is influenced by phenotype (body size). For this, we conducted a cross-fostering experiment in a two-patch system: an oakwood and a conifer plantation separated by only 1 km.
Experimental birds mostly returned to breed in the forest patch where they were raised, whether it was that of their genetic or their foster parents, indicating that decisions on where to settle are determined by individuals’ experience in their natal site, rather than by their genetic background. Nevertheless, nearly a third (27.6 %) moved away from the rearing habitat and, as previously observed in unmanipulated individuals, dispersal between habitats was phenotype-dependent. Pied flycatchers breeding in the oak and the pine forests are differentiated by body size, and analyses of genetic variation at microsatellite loci now provide evidence of subtle genetic differentiation between the two populations. This suggests that phenotype-dependent dispersal may contribute to population structure despite the short distance and widespread exchange of birds between the study plots.
Taken together, the current and previous findings that pied flycatchers do not always settle in the habitat to which they are best suited suggest that their strong tendency to return to the natal patch regardless of their body size might lead to maladaptive settlement decisions and thus constrain the potential of phenotype-dependent dispersal to promote microgeographic adaptation.
Selection of breeding environments is an important determinant of individual fitness and has therefore been a topic of considerable study in the realm of the ecology and evolution of dispersal [1–5]. Recent literature indicates that dispersal is typically nonrandom, and there is accumulating evidence that individuals often display a preference for a specific habitat type [6, 7]. Habitat preferences can have different underlying causes that include both genetic and environmental factors (e.g. [8–10]). First, experience in the natal patch can shape habitat preferences in adulthood, so that individuals tend to return to their birthplace or to other breeding habitats that resemble those they encountered at an early age, a phenomenon called natal habitat preference induction [11, 12]. Second, habitat preferences can be genetically determined [13, 14]. Finally, individuals may modify habitat selection according to their particular phenotype to settle in the habitats they are best suited to .
Numerous theoretical studies suggest that individual variation in habitat preferences can play an important role in population differentiation and ultimately also sympatric speciation (reviewed by ). Nevertheless, recent simulation-based studies of the evolutionary consequences of habitat preferences suggest that these may be largely contingent on the underlying preference mechanism [16, 17]. For example, natal habitat preference induction could promote population differentiation because the offspring of dispersers settling in a new, previously unused habitat may become reproductively isolated very quickly from the original source population. However, when the imprinting mechanism entails a substantial cost (e.g. energetic or nutritional costs of information processing and storing; ), population differentiation is more likely to occur through genetically-determined preferences . It should be noted that different preference mechanisms may act redundantly in nature, producing similar dispersal patterns, and even operate synergistically to promote population differentiation . For example, phenotype-dependent and natal habitat preferences may reinforce each other, since the former can facilitate local adaptation while the latter typically contributes to reproductive isolation [15, 19].
Different mechanisms generating habitat preference and nonrandom dispersal have been well characterized from a theoretical standpoint [6, 7, 11]. In addition, some efforts have been made to quantify the heritability of the propensity to disperse [14, 20] and of dispersal distance [14, 21] (but see ), and to assess the influence of early life experience [23–26] and phenotypic traits [27, 28] on dispersal and settlement decisions. However, almost no studies have been conducted to determine the relative contribution of different preference mechanisms to the expression of breeding site selection within a biologically realistic framework (but see ).
In this study, we tested the role of early life experience and genetic background in the settlement patterns of pied flycatchers (Ficedula hypoleuca) breeding in two contrasting habitats: a coniferous forest and a nearby deciduous forest separated by only 1 km. Males, and to a lesser degree, females breeding in the deciduous forest are morphologically different from those breeding in the nearby coniferous forest, the latter being smaller in size [28, 29]. Over 25 % of pied flycatchers returning to breed for the first time in the study area change habitats through dispersal and, from that moment, they rarely change the patch where they first bred . Male pied flycatchers disperse between the two forests according to body size, so that individuals moving from the coniferous to the deciduous patch are larger than those moving the other way round, and also than those that remain in the coniferous patch. We have also shown that, contrary to the deciduous forest, where male size does not determine fitness, the latter increases non-linearly with body size in the coniferous forest. Nevertheless, the observed patterns of phenotype-dependent dispersal and settlement do not translate into fitness benefits . Neither the proportion of individuals that return to their natal patch is what one would expect based on nest-site availability in each patch , suggesting that additional preference mechanisms (e.g. natal habitat preference induction or genetically-determined preference) must be operating.
Cross-fostering experiments between alternative habitats provide an essential tool to assess whether individuals settle to breed in one forest type or another owing to their genetic background or either mostly on the basis of the environment they experienced during sensitive periods along their ontogeny [30, 31]. Here, we capitalize on this approach to separate these two types of effects. Note that, since only one patch per habitat type is considered in this study, the terms “natal patch” and “natal habitat” are used interchangeably in the context of this work. No specific assessment is therefore made to separate the potential natal habitat effects on settlement decisions from those attributable to the specific features of each patch. As non-random, phenotype-dependent movements sustained across time may contribute to population structuring even at small spatial scales [27, 32, 33], we additionally explore whether there is detectable genetic divergence between pied flycatchers breeding in the two adjacent patches. Finally, we integrate our findings with earlier work on the same system concerning size-dependent dispersal and fitness differences among differently-sized individuals and discuss the relative importance of concurrent preference mechanisms in the spatial assortment of genotypes and phenotypes.
Fieldwork was conducted between 2006 and 2015 in a two-patch system: a mature oak (Quercus pyrenaica) forest of 9.3 ha, and a nearby (1.1 km) mixed pine (mostly Pinus sylvestris) plantation of 4.8 ha separated by unsuitable breeding habitat for pied flycatchers (i.e. a mixture of rock outcrops and riverside vegetation holding few suitable holes). Many aspects of the ecology of pied flycatchers have been studied in the deciduous and the coniferous forest since wooden nest-boxes (156 and 81) were provided in 1984 and 1988, respectively (e.g. ). Pied flycatchers are small (ca. 12 g) long-distance migrants overwintering south of the Sahara and breeding across Europe . Males often arrive from spring migration before females, search for a suitable cavity for nesting and announce themselves to females by singing actively [34, 35]. Pied flycatchers exhibit strong natal and breeding site fidelity , and local recruitment rates are the highest reported in the literature on the species [37, 38]. Males, and to a lesser degree females, may postpone breeding until their second or, more rarely, third year of life .
After arrival of the earliest males from spring migration, nest-boxes were checked at 1 − 3-day intervals to determine exact laying dates, clutch sizes, hatching dates and fledging success. Breeding birds were caught on day 8 (hatch day = 1) using a nest box trap, and measured for tarsus length (to the nearest 0.05 mm) and body mass (to the nearest 0.1 g). Fledgling mass and tarsus length were measured on day 13, when nestlings have already attained the definitive adult size [39, 40]. A small (<50 μL) blood sample was taken from the brachial vein of nestlings for molecular sexing and other genetic analyses (see e.g. ). Both adults and fledglings were individually marked with numbered metal rings. Nestling mortality was controlled from hatching to fledging by considering both the number of young surviving at day 13 and those found dead in the nest 18–20 days after hatching.
Cross-fostering was conducted in 2006–2009 (complete clutches and half broods; 2–3 nestlings) and 2012–2013 (half broods; 2–3 nestlings). Note that recruits stemming from complete and half clutches or broods swapped between nests are equally suitable for the purpose of this study. We cross-fostered birds both within and between habitats in order to test, on one hand, how likely it is for a bird raised in its habitat of origin to return there to breed (i.e. control group) and, on the other hand, how likely it is for a bird experimentally raised in a habitat other than its habitat of origin to move back to its habitat of origin (i.e. experimental group). Cross-fostering was carried out 2–3 days after the onset of incubation (clutches) or on day 2–3 after hatching (broods). All transfers were performed between nests of similar phenology (i.e. matching hatching dates, and a difference of ± 1 day for eggs). Nestlings were individually marked for individual identification with a unique combination of colour markings on underparts (legs and wings) made with non-toxic felt pens and repainted every 2–3 days. To avoid skewing the masses of experimental nestlings in any direction during partial cross-fostering and thus minimize non-deliberate effects of competition between native and foster birds, all nestlings were weighed, ranked for weight, and then sequentially assigned to be either cross-fostered or remaining in the source nest.
Experimental manipulations had no apparent effect on recruitment rates (see ‘results’). For all recruits stemming from cross-fostering experiments that were subsequently captured as first-time breeders, we recorded whether they returned to breed in their rearing patch or, conversely, moved into the adjacent one. It should be noted that, because the two habitat patches sampled in this study are isolated from other patches by unsuitable habitat, pied flycatchers returning to breed in the study area must choose between settling in their natal or the adjacent plot, or risk not breeding at all .
To examine genetic differences between pied flycatchers breeding in the two plots we followed the rationale described by . Our dataset included 351 males and females breeding in 2005 and 2006 in the coniferous and the deciduous forest (108 and 243 birds, respectively), genotyped at 15 microsatellite loci (see ). Test for linkage disequilibrium were performed in Genepop 4.0  and subsequently adjusted with Bonferroni corrections for multiples tests. We also used Genepop to test for deviation from Hardy-Weinberg equilibrium across loci and thus for the occurrence of inbreeding (Fis, estimated following ; Markov chain parameters were: 10,000 dememorisations, 1000 batches and 10,000 iterations per batch). Within each population, genetic diversity was calculated across loci using FSTAT  and compared using a Wilcoxon signed-rank test.
We used GENETIX to estimate the extent of genetic differentiation (F ST) between study plots (5000 permutations; ). Further, we tested for genotypic differentiation between populations in Genepop (“exact G test” option; Markov chain parameters were: 10,000 dememorisations, 1000 batches and 10,000 iterations per batch). Finally, we used a Bayesian approach as implemented in the program STRUCTURE  to infer the number of genetic clusters (K) in the population. Simulations were run assuming the admixture model with correlated allele frequencies. Four independent runs (K = 4), with twenty replicates for each K, were performed with a 105 burning period followed by 106 MCMC repeats after burning. To obtain the true value of K based on the method described by  we used the STRUCTURE HARVESTER web .
To investigate the roles of early experience and genetic variation in determining selection of either habitat by adult pied flycatchers in the two-patch system, we fitted a generalized linear mixed model (GLMM; binomial error distribution and logit link function) including birds’ propensity to return to the rearing patch (0 = return to breed in the forest patch of their genetic parents; 1 = return to breed in the forest patch of their foster parents) as the response variable. As fixed effects, we included sex, the type of experimental treatment (cross-fostering within or between patches) as a control term, and the life stage of cross-fostering (egg or nestling) to account for possible ontogenetic effects on imprinting. Laying date at first reproduction, which is an accurate proxy for the arrival date in our study population , was also included as a covariate to account for the potential effect of arrival time (e.g. decreasing nest-site availability as the season progresses) on bird distribution across the two forests . First-time breeders aged 2–3 years could have visited the study area in one or two seasons before being detected . However, when restricting the analysis to birds breeding for the first time at age 1 (n = 43), results were exactly the same (details not shown), suggesting that the age at first breeding does not affect the settlement pattern of pied flycatchers.
It is possible that some fledglings had visited (and possibly imprinted on) the adjacent forest prior to departure to the wintering areas, since post-fledging explorations may exceed the minimum distance separating the two plots (see [52, 53]. As the direction of exploratory trips is assumed to be random , the likelihood of visiting the adjacent plot should increase with proximity of the latter to the rearing nest, thus obscuring the interpretation of our results. Accordingly, we used GPS coordinates to calculate the minimum linear distance between the rearing nest and the adjacent forest and included it in the model as a covariate.
Migratory birds may rely on a host of environmental cues and constraints to guide their settlement decisions after arrival, such as nest-site  and food availability , predation risk  or conspecific attraction . Nest-site availability, measured as the annual number of nest-boxes not occupied by other hole-nesting species, differs between the two study patches and, furthermore, the magnitude of such differences may change over time . Likewise, it is likely that other settlement cues, although not measured in this population, may vary over time and space, and this could severely bias the settlement decisions of pied flycatchers. To account for all these factors (i.e. between-year variability in habitat heterogeneity), we included in the model return year as a random effect. Nest identity was also fitted as a random effect to account for repeated measurements of experimental nests that produced ≥2 recruits (17 %, n = 99).
Finally, we tested the interaction between tarsus length and rearing patch according to the notion by  that, when phenotype-dependent dispersal occurs, an individual’s propensity to disperse into a non-natal habitat may depend on the interaction between its morphology and its rearing habitat.
GLMMs were fitted in R 2.14.0 (http://www.R-project.org) using the function lmer in the package ‘lme4’ . Selection of the final model (containing only statistically significant terms) was carried out by dropping non-significant terms from a fully saturated model (containing all main effects and interactions) in a hierarchical way, starting with the least significant order terms.
Between 2006 and 2013, a total of 445 individuals from 116 nests (58 dyads) and 496 individuals from 120 nests (60 dyads) were cross-fostered between and within patches, respectively. Overall recruitment rates of cross-fostered young from experimental and control nests were, respectively, 14.7 % (57 recruits from 389 confirmed fledged offspring) and 11.8 % (48 recruits from 406 confirmed fledged offspring), the difference being not statistically significant (Pearson’s Chi-squared test with Yates’ continuity correction, χ 2 = 1.152, d.f. = 1, P = 0.28). Recruitment rates did not differ between cross-fostered and non-manipulated fledglings either (respectively, 13.2 % (n = 795) and 11.7 % (n = 1658); χ 2 = 1.003, d.f. = 1, P = 0.32), indicating that cross-fostering had no effect on the natural recruitment rates of pied flycatchers. Nests subjected to egg or nestling translocations showed similar recruitments rates (respectively, 13 % (n = 502 eggs) and 13.7 % (n = 293 nestlings); χ 2 = 0.030, d.f. = 1, P = 0.83), indicating that the life stage at which the experiment was carried out did not bias our estimates. Overall recruitment rates did not differ either between cross-fostered males and females (respectively, 38 out of 335 vs. 44 out of 286 genetically-sexed nestlings; χ 2 = 1.860, d.f. = 1, P = 0.17).
Results of the GLMM (binomial error distribution and logit link function) analyzing the effects of the experimental treatment (cross-fostering within and between habitats), life stage of cross-fostering (egg and nestling), sex, proximity to the adjacent habitat patch, breeding date, rearing patch, body size (tarsus length) and the interaction between rearing patch and body size on birds’ propensity to return to the habitat where they had been raised (0 = return to breed in the non-rearing patch; 1 = return to breed in the rearing patch)
Life stage of cross-fostering
Proximity to adjacent patch
Tarsus length x Rearing patch
Results from the GLMM showed no effect of sex, breeding date, age at first reproduction, or proximity to the adjacent forest patch on birds’ propensity to return to the same patch where they were fostered (Table 1). The interaction between body size (tarsus length) and rearing patch had no statistically significant effect. Body size or patch effects alone were not influential either, indicating that there was no overall tendency for a particular phenotype to exchange between patches more than others, nor were individuals reared in one particular patch more likely to move into the adjacent patch (Table 1).
F ST values revealed a low but statistically significant genetic differentiation between pied flycatchers breeding in the coniferous and the deciduous forest (F ST = 0.002, P < 0.001), which was corroborated by the exact test of genotypic differentiation (χ 2 = 53.99, d.f. = 30, P = 0.004). However, unlike in the coniferous habitat (χ 2 = 22.82, d.f. = 30, P = 0.82), the deciduous forest population showed a deviation from H-W equilibrium (χ 2 = 43.9, d.f. = 30, P = 0.05) due to an excess of homozygotes in three markers (Fhy 301, Fhy 339 and Fhy 401), which might affect the F ST estimates. Further analyses showed that the excess of homozygotes in those loci was due to rare, private alleles appearing at low frequencies. After removing from the analyses the individuals carrying rare alleles in the deciduous patch, the two populations showed no evidence for H-W disequilibrium (deciduous patch: χ 2 = 37.79, d.f. = 30, P = 0.15; coniferous patch: χ 2 = 25.6, d.f. = 30, P = 0.69), had similar genetic diversity (Z = 0.78, P = 0.43) and, although marginally non-significant, F ST values were of similar magnitude as that obtained from the entire data set (F ST = 0.001, P = 0.06). These results suggest that rare alleles in some individuals from the deciduous forest likely contributed to, but were not fully responsible for, the genetic differentiation of pied flycatchers between the study plots. Bayesian clustering approach did not provide support for two genetic clusters. This is not surprising since, as shown by simulations , STRUCTURE may not appropriately summarize population structure in scenarios with a low number of genetic clusters (n = 2) and low genetic differentiation, as occurs in our system.
By following individuals that had been cross-fostered to the same or a different habitat patch type, we provided evidence that pied flycatchers returned mostly to breed as adults in the area where they had fledged, regardless of their sex and origin. This is consistent with the hypothesis that individuals’ decisions on where to settle are based on the environment they experienced at an early age, rather than on their genetic background. Birds breeding in the coniferous and the deciduous forest differ in body size [28, 29] and, as shown here, there is also evidence of genetic differentiation between the two populations. Combined with our previous work demonstrating that bird movement between the study plots is affected by body size , our results suggest that phenotype-dependent dispersal may act concurrently with natal site fidelity and thus contribute to the observed population differentiation of pied flycatchers over a very small geographic scale.
Spatial and temporal scale of imprinting
Our study supports earlier evidence from experimental [61, 62] and observational studies  showing that pied flycatchers tend to return as adults to their natal patch, and that this is not genetically pre-programmed. However, it is not known whether the returning birds seek for a natal-like habitat or simply for their birthplace, as they only had one patch of each habitat type to choose between but, according to other studies on this species [61, 62], the former seems more likely. On the other hand, the observed propensity of pied flycatchers to settle in their natal patch could alternatively be interpreted in terms of spatial constraints − rather than choice − if, for example, nest-site availability away from the natal site is limited  or if it is easier to return to the same site after migration . Nevertheless, it is unlikely that pied flycatchers are constrained to choose between the two study plots since, on one hand, the distance between them is extremely short, thus facilitating bird exchange between both sites and, on the other hand, nest-site availability does not yet seem to be a limiting factor for pied flycatchers in the study area [28, 29].
It is known that experience in the natal site can strongly influence later habitat choice in the pied flycatcher [54, 62], but evidence from translocation experiments suggests that imprinting may also occur well beyond fledging.  transferred caged 5–6 weeks old fledglings 250 km away from their birthplace in northern Germany and found that all recruits returned to the area where they had been released, not to their natal patch, indicating that imprinting may extend over the fledgling stage. However, such tendency to return to the site of release after translocation might alternatively result from the typically high costs of finding the way home .
Besides the duration of the imprinting process, the spatial scale of post-fledging exploration may affect settlement patterns in the subsequent years. Male pied flycatchers commonly settle closer to their natal sites than females, although both sexes seem to imprint on a similarly small area (i.e. several kilometres in diameter; ). Mean distance of post-fledging dispersal has been estimated between 0.6 − 1.4 km, depending on the study population [52, 53]. Exploratory trips by fledglings that extend beyond the distance separating the two forests (1.1 km) might facilitate imprinting on the patch adjacent to their natal patch, and thus determine future habitat shifts. However, even though detailed data on prospecting behaviour by fledglings are not available for the study area, this seems unlikely to bias our results, as the likelihood of changing habitats was apparently not determined by proximity of the rearing nest to the neighbour forest.
Some individuals could have visited the study area in one or more seasons before reproducing for the first time and explored future territories away from the area they explored as a fledgling . But if this occurred, it certainly did not confound our estimates, as there were no differences between first-year and older first breeders in the propensity to change habitats. In addition, the timing of arrival from spring migration, which is known to have a profound effect on the breeding phenology and success of pied flycatchers , might influence territory acquisition. Late migrants could be time-constrained to choose their breeding site and thus occupy the remaining free territories, whatever the habitat type . This may be particularly true for first-year breeders, which usually show delayed arrival dates [35, 37]. However, we found no phenological or age effects on birds’ propensity to return to their rearing patch. Finally, an additional potentially confounding factor is that some individuals could perceive certain forest patches as offering better foraging and nest-site opportunities than others based on e.g. asymmetry in breeding density and presence of competitors, and settle accordingly . However, we think this is unlikely to bias our estimates for two main reasons. First, as shown here, males and females are equally likely to return to breed in the area where they were raised, despite the fact that only the former compete for nest sites . Second, earlier results on the same system indicate that nest-site availability and population density do not play an important role in the local distribution of pied flycatchers [28, 29].
Natal habitat imprinting vs. phenotype-dependent dispersal
So far, few studies have attempted to tease apart the relative effects of natal patch preferences and phenotype-dependent dispersal on population differentiation.  also reported a strong tendency in lake and stream sticklebacks to return to their natal area and, as in pied flycatchers, their movement between habitats was phenotype-dependent. In our study population, around a quarter of the (unmanipulated) natal dispersers do not return to the forest patch they presumably imprinted on. Conversely, they move naturally into the adjacent habitat patch according to their body size, so that individuals dispersing from the coniferous to the deciduous forest are larger than those moving the other way round (respectively, 19.47 ± 0.07 mm and 19.39 ± 0.08 mm; mean ± SE; ). Likewise, we observed qualitatively similar differences between cross-fostered individuals moving from the coniferous to the deciduous forest and those moving the other way round (respectively, 19.52 ± 0.08 mm and 19.28 ± 0.15 mm; mean ± SE). Even though this difference (1.1 %) might seem modest, it is well within the range reported in other studies made at small spatial scales (5–10 km; 0.9–1.2 % difference; [33, 64, 65]. Nevertheless, in contrast to the large sample of unmanipulated individuals, the effect of the interaction between body size and rearing patch on the propensity to change habitats did not reach significance in our experimental birds (see ‘results’), possibly because the sample size of ‘non-philopatric’ recruits stemming from the experiment is very limited compared with the observational study (29 vs. 115). Indeed, a power analysis based on the observed differences between the two groups of dispersers (effect size d = 0.52) showed that, for α = 0.05, an overall sample size of n = 120 would be required to have 80 % power, which means that, according to the actual recruitment rate (13.2 %) and proportion of recruits moving between habitats (27.6 %), over 3200 individuals would need to be cross-fostered. Still, taken together, the results of this experiment and our previous studies [28, 29] give support to the notion by  that phenotype-dependent dispersal and natal site fidelity may act concurrently.
More pied flycatchers disperse between the two habitats compared with the sticklebacks study (30 % vs. 10 %) but, unlike in the stickleback populations, phenotype-dependent dispersal and natal habitat preferences of pied flycatchers do not appear to act synergistically but partially cancel out each other.  showed that morphologically different sticklebacks preferentially settle in the habitat conferring a fitness advantage (ultimately promoting microgeographic adaptation). By contrast, decisions on where to settle and fitness appear to be decoupled in our study population , possibly because of the strong tendency to return to the natal habitat patch, not to that in which their particular phenotype performs best (see [26, 66]; but see . Thus, the strong (maladaptive) reluctance of the majority of birds to move away from either of the two habitat patches might constrain, rather than promote, microgeographic adaptation and population divergence. Even with this constraint, the level of genetic differentiation between both pied flycatcher populations is of the same order of magnitude as that reported between the less mobile lake and stream sticklebacks (F ST = 0.008; ), and also similar to those found over much more extensive regions in other highly mobile organisms, such as the house sparrow (Passer domesticus) (F ST = 0.004 across Finland; ) or the Eurasian reed warblers (Acrocephalus scirpaceus) (F ST = 0.013 across Europe; ). Accordingly, it seems that the diversifying effects of phenotype-dependent dispersal alone might be strong enough to generate a detectable signal of population differentiation at exceedingly small geographic scales [27, 33].
Early experience in the natal patch may play a crucial role in determining subsequent dispersal and settlement in the pied flycatcher ([54, 61, 69], this study) and possibly also in a broad range of animals . Taken together, the results presented herein support previous evidence indicating that phenotype-dependent dispersal between the two plots might contribute to the observed morphological  and genetic (this study) differentiation of pied flycatchers. However, the potential of phenotype-dependent dispersal to increase the magnitude of divergence between these two populations might be constrained by the strong natal site fidelity, as individual performance in each forest patch is strongly determined by morphology . We suggest that the heretofore largely neglected − but likely widespread − interplay between early experience in the natal site and individual phenotype should be fully taken into account in future studies investigating the mechanisms underlying non-random dispersal and habitat selection.
We are grateful to Inés Valencia, María Cuenca, Alba Ruiz, Marta Guntiñas, Tomás Redondo, Gregorio Moreno-Rueda, Gerardo Jiménez, Carlos Marfil, Guadalupe Macías, Fran Romero and David Ochoa for their assistance with fieldwork and to Alex Centeno for help with genetic analyses.
CC received financial support from the Spanish Ministry of Economy and Competitiveness, through the Severo Ochoa Programme for Centres of Excellence in R&D&I (SVP-2013-067686). DC was supported by project CGL 2009-10652. Since 1987, JP’s work has been funded by the Spanish Governments, most recently by projects CGL2015-70639-P, CGL2014-55969-P and CGL2012-35232.
Availability of data and materials
The authors confirm that all data underlying the findings will be made publicly available at https://digital.csic.es/handle/10261/135062 in case the paper is accepted.
CC, DC and JP conceived and performed the cross-fostering experiment, collected and analyzed the data and wrote the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Research was conducted under permits from the authorities in Madrid and Castilla-La Mancha autonomies. The Ethics Committee on Animal Experimentation from Doñana Biological Station (CSIC) approved the experimental procedure, which complies with the current laws of Spain.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Johnson ML, Gaines MS. Evolution of dispersal: theoretical models and empirical tests using birds and mammals. Ann Rev Ecol Syst. 1990;21:449–80.View ArticleGoogle Scholar
- Szulkin M, Sheldon BC. Dispersal as a means of inbreeding avoidance in a wild bird population. Proc R Soc Lond B. 2008;275:703–11.View ArticleGoogle Scholar
- Duckworth RA, Badyaev AV. Coupling of dispersal and aggression facilitates the rapid range expansion of a passerine bird. Proc Natl Acad Sci U S A. 2007;104:15017–22.View ArticlePubMedPubMed CentralGoogle Scholar
- Duckworth RA, Kruuk LEB. Evolution of genetic integration between dispersal and colonization ability in a bird. Evolution. 2009;63:968–77.View ArticlePubMedGoogle Scholar
- Clobert J, Baguette M, Benton TG, Bullock JM. Dispersal ecology and evolution Oxford University press Oxford. 2012.View ArticleGoogle Scholar
- Clobert J, Le Galliard L, Cote J, Meylan S, Massot M. Informed dispersal heterogeneity in animal dispersal syndromes and the dynamics of spatially structured populations. Ecol Lett. 2009;12:197–209.View ArticlePubMedGoogle Scholar
- Edelaar P, Bolnick DI. Non-random gene flow: an underappreciated force in evolution and ecology. Trends Ecol Evol. 2012;27:659–65.View ArticlePubMedGoogle Scholar
- Massot M, Clobert J. Processes at the origin of similarities in dispersal behavior among siblings. J Evol Biol. 2000;13:707–19.View ArticleGoogle Scholar
- Cote J, Clobert J. Social personalities influence natal dispersal in a lizard. Proc R Soc Lond B. 2007;274:383–90.View ArticleGoogle Scholar
- Duckworth RA. Maternal effects and range expansion: a key factor in a dynamic process? Phil Trans R Soc B. 2009;364:1075–86.View ArticlePubMedPubMed CentralGoogle Scholar
- Davis JM, Stamps JA. The effect of natal experience on habitat preferences. Trends Ecol Evol. 2004;19:411–6.View ArticlePubMedGoogle Scholar
- Stamps JA, Krishnan VV, Willits NH. How different types of natal experience affect habitat preference. Am Nat. 2009;174:623–30.View ArticlePubMedGoogle Scholar
- Jaenike J, Holt RD. Genetic variation for habitat preference: evidence and explanations. Am Nat. 1991;137:S67–90.View ArticleGoogle Scholar
- Doligez B, Gustafsson L, Pärt T. ‘Heritability’ of dispersal propensity in a patchy population. Proc R Soc Lond B. 2009;276:2829–36.View ArticleGoogle Scholar
- Edelaar P, Siepielski AM, Clobert J. Matching habitat choice causes directed gene flow: a neglected dimension in evolution and ecology. Evolution. 2008;62:2462–72.View ArticlePubMedGoogle Scholar
- Beltman JB, Metz JAJ. Speciation: more likely through a genetic or through a learned habitat preference? Proc R Soc Lond B. 2005;272:1455–63.View ArticleGoogle Scholar
- Berner D, Thibert‐Plante X. How mechanisms of habitat preference evolve and promote divergence with gene flow. J Evol Biol. 2015;28:1641–55.View ArticlePubMedGoogle Scholar
- Dukas R. Costs of memory: ideas and predictions. J Theor Biol. 1999;197:41–50.View ArticlePubMedGoogle Scholar
- Bolnick DI, Snowberg LK, Patenia C, Stutz WE, Ingram T, Lau OL. Phenotype-dependent native habitat preference facilitates divergence between parapatric lake and stream stickleback. Evolution. 2009;63:2004–16.View ArticlePubMedGoogle Scholar
- Hansson B, Bensch S, Hasselquist D. Heritability of dispersal in the great reed warbler. Ecol Lett. 2003;6:290–4.View ArticleGoogle Scholar
- Greenwood PJ, Harvey PH, Perrins CM. The role of dispersal in the great tit (Parus major): the causes consequences and heritability of natal dispersal. J Anim Ecol. 1979;48:123–42.View ArticleGoogle Scholar
- van Noordwijk AJ. Problems in the analysis of dispersal and a critique on its ‘heritability’ in the great tit. J Anim Ecol. 1984;53:533–44.View ArticleGoogle Scholar
- Alonso JC, Martín E, Alonso JA, Morales MB. Proximate and ultimate causes of natal dispersal in the great bustard Otis tarda. Behav Ecol. 1998;9:243–52.View ArticleGoogle Scholar
- Ekman J, Eggers S, Griesser M. Fighting to stay: the role of sibling rivalry for delayed dispersal. Anim Behav. 2002;64:453–9.View ArticleGoogle Scholar
- Cote J, Clobert J, Brodin T, Fogarty S, Sih A. Personality-dependent dispersal: characterization ontogeny and consequences for spatially structured populations. Phil Trans R Soc B. 2010;65:4065–76.View ArticleGoogle Scholar
- Piper WH, Palmer MW, Banfield N, Meyer MW. Can settlement in natal-like habitat explain maladaptive habitat selection? Proc R Soc Lond B. 2013;280(1765):20130979.View ArticleGoogle Scholar
- Garant D, Kruuk LEB, Wilkin TA, McCleery RH, Sheldon BC. Evolution driven by differential dispersal within a wild bird population. Nature. 2005;433:60–5.View ArticlePubMedGoogle Scholar
- Camacho C, Canal D, Potti J. Non-random dispersal drives phenotypic divergence within a bird population. Ecol Evol. 2013;3:4841–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Camacho C, Canal D, Potti J. Testing the matching habitat choice hypothesis: phenotype-environment correlation and fitness in a wild bird population. Evol Ecol. 2015;29:873–86.View ArticleGoogle Scholar
- Merilä J, Sheldon BC. Avian quantitative genetics. Curr Ornithol. 2001;16:179–255.Google Scholar
- Stinchcombe JR. Cross-pollination of plants and animals: wild quantitative genetics and plant evolutionary genetics. In: Charmantier A, Garant D, Kruuk LEB, editors. Quantitative genetics in the wild. Oxford: Oxford University Press; 2014. p. 128–46.View ArticleGoogle Scholar
- Postma E, van Noordwijk AJ. Gene flow maintains a large genetic difference in clutch size at a small spatial scale. Nature. 2005;433:65–8.View ArticlePubMedGoogle Scholar
- Senar JC, Borrás A, Cabrera J, Cabrera T, Björklund M. Local differentiation in the presence of gene flow in the citril finch Serinus citrinella. Biol Lett. 2006;2:85–7.View ArticlePubMedGoogle Scholar
- Lundberg A, Alatalo RV. The pied flycatcher T. London: AD Poyser; 1992.Google Scholar
- Potti J. Arrival time from spring migration in male pied flycatchers: individual consistency and familial resemblance. Condor. 1998;100:702–8.View ArticleGoogle Scholar
- Montalvo S, Potti J. Breeding dispersal in Spanish pied flycatchers Ficedula hypoleuca. Ornis Scand. 1992;23:491–8.View ArticleGoogle Scholar
- Potti J, Montalvo S. Return rate age at first breeding and natal dispersal of pied flycatchers Ficedula hypoleuca in central Spain. Ardea. 1991;79:419–28.Google Scholar
- Canal D, Serrano D, Potti J. Exploring heterozygosity-survival correlations in a wild songbird population: contrasting effects between juvenile and adult stages. PLoS One. 2014;9(8):e105020.View ArticlePubMedPubMed CentralGoogle Scholar
- Alatalo RV, Lundberg A. Heritability and selection on tarsus length in the pied flycatcher (Ficedula hypoleuca). Evolution. 1986;40:574–83.View ArticleGoogle Scholar
- Potti J, Merino S. Heritability estimates and maternal effects on tarsus length in pied flycatchers Ficedula hypoleuca. Oecologia. 1994;100:331–8.View ArticleGoogle Scholar
- Postma E, Den Tex RJ, van Noordwijk AJ, Mateman AC. Neutral markers mirror small‐scale quantitative genetic differentiation in an avian island population. Biol J Linn Soc. 2009;97:867–75.View ArticleGoogle Scholar
- Raymond M, Rousset F. Genepop (version 1 2) population genetics software for exact tests and ecumenicism. J Hered. 1995;86:248–9.Google Scholar
- Weir BS, Cockerham CC. Estimating F-statistics for the analysis of population structure. Evolution. 1984;38:1358–70.View ArticleGoogle Scholar
- Goudet J. FSTAT (version 1 2): a computer program to calculate F-statistics. J Hered. 1995;86:485–6.Google Scholar
- Belkhir K, Borsa P, Chikhi L, Raufaste N, Bonhomme F. GENETIX 4 05: logiciel sous Windows TM pour la génétique des populations Laboratoire Génome Populations Interactions CNRS UMR 5171: Université de Montpellier II Montpellier France; 2004Google Scholar
- Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155:945–59.PubMedPubMed CentralGoogle Scholar
- Evanno G, Regnaut S, Goudet J. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol. 2005;14:2611–20.View ArticlePubMedGoogle Scholar
- Earl DA, VonHoldt BM. STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv Genet Resour. 2012;4:359–61.View ArticleGoogle Scholar
- Potti J, Montalvo S. Male arrival and female mate choice in Pied Flycatchers Ficedula hypoleuca in central Spain. Ornis Scand. 1991;22:45–54.View ArticleGoogle Scholar
- Smith RJ, Moore FR. Arrival timing and seasonal fitness in a long-distance migratory landbird. Behav Ecol Sociobiol. 2005;57:231–9.View ArticleGoogle Scholar
- Becker PH, Dittmann T, Ludwigs JD, Limmer B, Ludwig SC, Bauch C, Braasch A, Wendeln H. Timing of initial arrival at the breeding site predicts age at first reproduction in a long-lived migratory bird. Proc Natl Acad Sci U S A. 2008;105:12349–52.View ArticlePubMedPubMed CentralGoogle Scholar
- Van Balen JH. Observations on the post-fledging dispersal of the pied flycatcher Ficedula hypoleuca. Ardea. 1979;67:134–7.Google Scholar
- Vysotsky VG, Bardin AV, Sokolov LV. Distance of post-fledging dispersal of the Pied Flycatcher Ficedula hypoleuca. Proc Zool Inst Leningrad. 1990;210:35–45.Google Scholar
- Chernetsov N, Sokolov LV, Kosarev V, Leoke D, Markovets M, Tsvey A, Shapoval AP. Sex-related natal dispersal of pied flycatchers: how far away from home? Condor. 2006;108:711–7.View ArticleGoogle Scholar
- Mänd R, Tilgar V, Lõhmus A. Providing nest boxes for hole-nesting birds–does habitat matter? Biodivers Conserv. 2005;14:1823–40.View ArticleGoogle Scholar
- Marshall MR, Cooper RJ. Territory size of a migratory songbird in response to caterpillar density and foliage structure. Ecology. 2004;85:432–45.View ArticleGoogle Scholar
- Morosinotto C, Thomson RL, Korpimäki E. Habitat selection as an antipredator behaviour in a multi‐predator landscape: all enemies are not equal. J Anim Ecol. 2010;79:327–33.View ArticlePubMedGoogle Scholar
- Doligez B, Danchin E, Clobert J, Gustafsson L. The use of conspecific reproductive success for breeding habitat selection in a non-colonial hole-nesting species the collared flycatcher. J Anim Ecol. 1999;68:1193–206.View ArticleGoogle Scholar
- Bates D, Maechler M, Bolker B. lme4: linear mixed-effects models using S4 classes R package version 0 999375–42. 2011.Google Scholar
- Kalinowski ST. The computer program STRUCTURE does not reliably identify the main genetic clusters within species: simulations and implications for human population structure. Heredity. 2011;106:625–32.View ArticlePubMedGoogle Scholar
- Berndt R, Winkel W. Verfrachtungs-Experimente zur Frage der Geburtsortsprägung beim Trauerschnäpper (Ficedula hypoleuca). J Ornithol. 1979;120:41–5.View ArticleGoogle Scholar
- Vallin N, Qvarnström A. Learning the hard way: imprinting can enhance enforced shifts in habitat choice. Int J Ecol. 2011;2011:287532.View ArticleGoogle Scholar
- Pärt T. The importance of local familiarity and search costs for age-and sex-biased philopatry in the collared flycatcher. Anim Behav. 1995;49:1029–38.View ArticleGoogle Scholar
- García-Navas V, Ferrer ES, Sanz JJ, Ortego J. The role of immigration and local adaptation on fine-scale genotypic and phenotypic population divergence in a less mobile passerine. J Evol Biol. 2014;27:1590–603.View ArticlePubMedGoogle Scholar
- Camacho C, Sáez P, Sánchez S, Palacios S, Molina C, Potti J. The road to opportunities: landscape change promotes body size divergence in a highly mobile species. Curr Zool. 2016;62:7–14.View ArticleGoogle Scholar
- Fletcher RJ, Robertson EP, Wilcox RC, Reichert BE, Austin JD, Kitchens WM. Affinity for natal environments by dispersers impacts reproduction and explains geographical structure of a highly mobile bird. Proc R Soc Lond B. 2015;282:20151545.View ArticleGoogle Scholar
- Kekkonen J, Seppä P, Hanski IK, Jensen H, Väisänen RA, Brommer JE. Low genetic differentiation in a sedentary bird: house sparrow population genetics in a contiguous landscape. Heredity. 2011;106:183–90.View ArticlePubMedGoogle Scholar
- Procházka P, Stokke BG, Jensen H, Fainová D, Bellinvia E, Fossøy F, Vikan JR, Bryja J, Soler M. Low genetic differentiation among reed warbler Acrocephalus scirpaceus populations across Europe. J Avian Biol. 2011;42:103–13.View ArticleGoogle Scholar
- Burger C, Both C. Translocation as a novel approach to study effects of a new breeding habitat on reproductive output in wild birds. PLoS One. 2011;6(3):e18143.View ArticlePubMedPubMed CentralGoogle Scholar