Genetic structure reflects natal dispersal movements at different spatial scales in the blue tit, Cyanistes caeruleus

Fine-scale genetic structure and phenotypic divergence of a passerine bird population inhabiting a continuous Mediterranean woodland

Genetic differentiation between populations inhabiting ecologically different habitats might appear because of limited dispersal and gene flow, which may lead to patterns of phenotypic divergence and local adaptation. In this study, we use dispersal, genotypic (24 microsatellite loci) and phenotypic (body size and clutch size) data to analyse patterns of genetic structuring and phenotypic divergence in a blue tit (Cyanistes caeruleus) population inhabiting a continuous and heterogeneous woodland along a valley. The two slopes of the valley differ in their forest formations and environmental conditions. Findings showed that most blue tits reproduced within their natal slope. Accordingly, microsatellite analyses revealed that populations of blue tits established in the two slopes show subtle genetic differentiation. The two genetic populations diverged in clutch size, exceeding the level of differentiation expected based on genetic drift, hence suggesting divergent selection (or other processes promoting divergence) on this life-history trait. Our findings reveal that restricted dispersal and spatial heterogeneity may lead to genetic differentiation among bird populations at a surprisingly small scale. In this respect, it is worth highlighting that such differentiation occurs for an organism with high dispersal capacity and within a continuous woodland. Moreover, we show that small-scale ecological differences, together with limited gene flow, can result in selection favouring different phenotypes even within the same continuum population.

Replicated, urban-driven exposure to metallic trace elements in two passerines

While there are increasing examples of phenotypic and genotypic differences between urban and non-urban populations of plants and animals, few studies identified the mechanisms explaining those dissimilarities. The characterization of the urban landscape, which can only be achieved by measuring variability in relevant environmental factors within and between cities, is a keystone prerequisite to understand the effects of urbanization on wildlife. Here, we measured variation in bird exposure to metal pollution within 8 replicated urbanization gradients and within 2 flagship bird species in urban evolutionary ecology: the blue tit (Cyanistes caeruleus) and the great tit (Parus major). We report on a highly significant, positive linear relationship between the magnitude of urbanization-inferred as either tree cover, impervious surface cover, or an urbanization score computed from several environmental variables, and copper, zinc and lead concentrations in bird feathers. The reverse relationship was measured in the case of mercury, while cadmium and arsenic did not vary in response to the urbanization level. This result, replicated across multiple cities and two passerine species, strongly suggests that copper, zinc, lead and mercury pollution is likely to trigger the emergence of parallel responses at the phenotypic and/or genotypic level between urban environments worldwide.

Genetic structure of urban and non-urban populations differs between two common parid species

Landscape conversions induced by human activities can affect dispersal patterns of various bird species and, as a result, affect genetic structure of their populations. Genetic differentiation of bird populations may be enhanced by habitat variation, especially in urban-non-urban systems. The majority of population genetic studies focus on single species, which inflicts limitations for direct comparisons of genetic responses of avian populations to urbanization. Here, we used a set of microsatellite markers to examine genetic diversity, gene flow and population structure in two common parid species, great tits Parus major and blue tits Cyanistes caeruleus occupying three sites in habitats with contrasting urbanization level in central Poland. We found low but significant divergence of urban park population with both suburban and non-urban forest great tit populations, while no differentiation was found between suburban forest and non-urban forest populations. In contrast, no evidence for genetic differentiation was found between blue tit populations from the urban park, suburban forest and non-urban forest sites. We conclude that great tits and blue tits respond to urbanization-related changes in a different way, which may be a result of different rates of migration and/or dispersal, likely higher in blue tits. Some impact may be also induced by interspecific competition. We suggest that changing the focus of urban genetic research from single to multiple species may provide novel insights into how natural populations respond to the processes of urbanization.

Demographic history and genomics of local adaptation in blue tit populations

Understanding the genomic processes underlying local adaptation is a central aim of modern evolutionary biology. This task requires identifying footprints of local selection but also estimating spatio‐temporal variations in population demography and variations in recombination rate and in diversity along the genome. Here, we investigated these parameters in blue tit populations inhabiting deciduous versus evergreen forests, and insular versus mainland areas, in the context of a previously described strong phenotypic differentiation. Neighboring population pairs of deciduous and evergreen habitats were weakly genetically differentiated (F_ST = 0.003 on average), nevertheless with a statistically significant effect of habitat type on the overall genetic structure. This low differentiation was consistent with the strong and long‐lasting gene flow between populations inferred by demographic modeling. In turn, insular and mainland populations were moderately differentiated (F_ST = 0.08 on average), in line with the inference of moderate ancestral migration, followed by isolation since the end of the last glaciation. Effective population sizes were large, yet smaller on the island than on the mainland. Weak and nonparallel footprints of divergent selection between deciduous and evergreen populations were consistent with their high connectivity and the probable polygenic nature of local adaptation in these habitats. In turn, stronger footprints of divergent selection were identified between long isolated insular versus mainland birds and were more often found in regions of low recombination, as expected from theory. Lastly, we identified a genomic inversion on the mainland, spanning 2.8 Mb. These results provide insights into the demographic history and genetic architecture of local adaptation in blue tit populations at multiple geographic scales.

Network Analysis Shows Asymmetrical Flows within a Bird Metapopulation

How the spatial expansion of a species changes at a human time scale is a process difficult to determine. We studied the dispersal pattern of the French white stork population, using a 21-year ringing/resighting dataset. We used the graph-theory to investigate the strength of links between 5 populations (North-East, North-West, Centre, West, and South) and to determine factors important for the birds’ movements. Two clusters of populations were identified within the metapopulation, with most frequent movements of individuals between North-Eastern and Centre populations, and between North-Western and Western populations. Exchanges of individuals between populations were asymmetrical, where North-Eastern and North-Western populations provided more emigrants than they received immigrants. Neither the geographical distance between populations, nor the difference in densities influenced the number of individuals exchanging between populations. The graph-theory approach provides a dynamic view of individual movements within a metapopulation and might be useful for future population studies in the context of conservation.

The strength of the association between heterozygosity and probability of interannual local recruitment increases with environmental harshness in blue tits

The extent of inbreeding depression and the magnitude of heterozygosity–fitness correlations (HFC) have been suggested to depend on the environmental context in which they are assayed, but little evidence is available for wild populations. We combine extensive molecular and capture–mark–recapture data from a blue tit (Cyanistes caeruleus) population to (1) analyze the relationship between heterozygosity and probability of interannual adult local recruitment and (2) test whether environmental stress imposed by physiologically suboptimal temperatures and rainfall influence the magnitude of HFC. To address these questions, we used two different arrays of microsatellite markers: 14 loci classified as neutral and 12 loci classified as putatively functional. We found significant relationships between heterozygosity and probability of interannual local recruitment that were most likely explained by variation in genomewide heterozygosity. The strength of the association between heterozygosity and probability of interannual local recruitment was positively associated with annual accumulated precipitation. Annual mean heterozygosity increased over time, which may have resulted from an overall positive selection on heterozygosity over the course of the study period. Finally, neutral and putatively functional loci showed similar trends, but the former had stronger effect sizes and seemed to better reflect genomewide heterozygosity. Overall, our results show that HFC can be context dependent, emphasizing the need to consider the role of environmental heterogeneity as a key factor when exploring the consequences of individual genetic diversity on fitness in natural populations.

Heterogeneous road networks have no apparent effect on the genetic structure of small mammal populations

Roads are widely recognized to represent a barrier to individual movements and, conversely, verges can act as potential corridors for the dispersal of many small mammals. Both barrier and corridor effects should generate a clear spatial pattern in genetic structure. Nevertheless, the effect of roads on the genetic structure of small mammal populations still remains unclear. In this study, we examine the barrier effect that different road types (4-lane highway, 2-lane roads and single-lane unpaved roads) may have on the population genetic structure of three species differing in relevant life history traits: southern water vole Arvicola sapidus, the Mediterranean pine vole Microtus duodecimcostatus and the Algerian mouse Mus spretus. We also examine the corridor effect of highway verges on the Mediterranean pine vole and the Algerian mouse. We analysed the population structure through pairwise estimates of FST among subpopulations bisected by roads, identified genetic clusters through Bayesian assignment approaches, and used simple and partial Mantel tests to evaluate the relative barrier or corridor effect of roads. No strong evidences were found for an effect of roads on population structure of these three species. The barrier effect of roads seems to be site-specific and no corridor effect of verges was found for the pine vole and Algerian mouse populations. The lack of consistent results among species and for each road type lead us to believe that the ability of individual dispersers to use those crossing structures or the habitat quality in the highway verges may have a relatively higher influence on gene flow among populations than the presence of crossing structures per se. Further research should include microhabitat analysis and the estimates of species abundance to understand the mechanisms that underlie the genetic structure observed at some sites.

Population genomic footprints of fine-scale differentiation between habitats in Mediterranean blue tits

Linking population genetic variation to the spatial heterogeneity of the environment is of fundamental interest to evolutionary biology and ecology, in particular when phenotypic differences between populations are observed at biologically small spatial scales. Here, we applied restriction-site associated DNA sequencing (RAD-Seq) to test whether phenotypically differentiated populations of wild blue tits (Cyanistes caeruleus) breeding in a highly heterogeneous environment exhibit genetic structure related to habitat type. Using 12 106 SNPs in 197 individuals from deciduous and evergreen oak woodlands, we applied complementary population genomic analyses, which revealed that genetic variation is influenced by both geographical distance and habitat type. A fine-scale genetic differentiation supported by genome- and transcriptome-wide analyses was found within Corsica, between two adjacent habitats where blue tits exhibit marked differences in breeding time while nesting < 6 km apart. Using redundancy analysis (RDA), we show that genomic variation remains associated with habitat type when controlling for spatial and temporal effects. Finally, our results suggest that the observed patterns of genomic differentiation were not driven by a small proportion of highly differentiated loci, but rather emerged through a process such as habitat choice, which reduces gene flow between habitats across the entire genome. The pattern of genomic isolation-by-environment closely matches differentiation observed at the phenotypic level, thereby offering significant potential for future inference of phenotype-genotype associations in a heterogeneous environment.

Mediterranean blue tits as a case study of local adaptation

While the study of the origins of biological diversity across species has provided numerous examples of adaptive divergence, the realization that it can occur at microgeographic scales despite gene flow is recent, and scarcely illustrated. We review here evidence suggesting that the striking phenotypic differentiation in ecologically relevant traits exhibited by blue tits Cyanistes caeruleus in their southern range‐edge putatively reflects adaptation to the heterogeneity of the Mediterranean habitats. We first summarize the phenotypic divergence for a series of life history, morphological, behavioural, acoustic and colour ornament traits in blue tit populations of evergreen and deciduous forests. For each divergent trait, we review the evidence obtained from common garden experiments regarding a possible genetic origin of the observed phenotypic differentiation as well as evidence for heterogeneous selection. Second, we argue that most phenotypically differentiated traits display heritable variation, a fundamental requirement for evolution to occur. Third, we discuss nonrandom dispersal, selective barriers and assortative mating as processes that could reinforce local adaptation. Finally, we show how population genomics supports isolation – by – environment across landscapes. Overall, the combination of approaches converges to the conclusion that the strong phenotypic differentiation observed in Mediterranean blue tits is a fascinating case of local adaptation.

Individual genetic diversity and probability of infection by avian malaria parasites in blue tits (Cyanistes caeruleus)

Understanding the importance of host genetic diversity for coping with parasites and infectious diseases is a long-standing goal in evolutionary biology. Here, we study the association between probability of infection by avian malaria (Plasmodium relictum) and individual genetic diversity in three blue tit (Cyanistes caeruleus) populations that strongly differ in prevalence of this parasite. For this purpose, we screened avian malaria infections and genotyped 789 blue tits across 26 microsatellite markers. We used two different arrays of markers: 14 loci classified as neutral and 12 loci classified as putatively functional. We found a significant relationship between probability of infection and host genetic diversity estimated at the subset of neutral markers that was not explained by strong local effects and did not differ among the studied populations. This relationship was not linear, and probability of infection increased up to values of homozygosity by locus (HL) around 0.15, reached a plateau at values of HL from 0.15 to 0.40 and finally declined among a small proportion of highly homozygous individuals (HL > 0.4). We did not find evidence for significant identity disequilibrium, which may have resulted from a low variance of inbreeding in the study populations and/or the small power of our set of markers to detect it. A combination of subtle positive and negative local effects and/or a saturation threshold in the association between probability of infection and host genetic diversity in combination with increased resistance to parasites in highly homozygous individuals may explain the observed negative quadratic relationship. Overall, our study highlights that parasites play an important role in shaping host genetic variation and suggests that the use of large sets of neutral markers may be more appropriate for the study of heterozygosity-fitness correlations.

Large-scale patterns of genetic variation in a female-biased dispersing passerine: the importance of sex-based analyses

Dispersal affects the distribution, dynamics and genetic structure of natural populations, and can be significantly different between sexes. However, literature records dealing with the dispersal of migratory birds are scarce, as migratory behaviour can notably complicate the study of dispersal. We used the barn swallow Hirundo rustica as model taxon to investigate patterns of genetic variability in males and in females of a migratory species showing sex-biased dispersal. We collected blood samples (n = 186) over the period 2006 to 2011 from adults (H. r. rustica subspecies) nesting in the same breeding site at either high (Ireland, Germany and Russia) or low (Spain, Italy and Cyprus) latitude across Europe. We amplified the Chromo Helicase DNA gene in all birds in order to warrant a sex-balanced sample size (92 males, 94 females). We investigated both uniparental (mitochondrial ND2 gene) and biparental (microsatellite DNA: 10 loci) genetic systems. The mtDNA provided evidence for demographic expansion yet no significant partition of the genetic variability was disclosed. Nevertheless, a comparatively distant Russian population investigated in another study, whose sequences were included in the present dataset, significantly diverged from all other ones. Different to previous studies, microsatellites highlighted remarkable genetic structure among the studied populations, and pointed to the occurrence of differences between male and female barn swallows. We produced evidence for non-random patterns of gene flow among barn swallow populations probably mediated by female natal dispersal, and we found significant variability in the philopatry of males of different populations. Our data emphasize the importance of taking into account the sex of sampled individuals in order to obtain reliable inferences on species characterized by different patterns of dispersal between males and females.

Breed locally, disperse globally: fine-scale genetic structure despite landscape-scale panmixia in a fire-specialist

An exciting advance in the understanding of metapopulation dynamics has been the investigation of how populations respond to ephemeral patches that go 'extinct' during the lifetime of an individual. Previous research has shown that this scenario leads to genetic homogenization across large spatial scales. However, little is known about fine-scale genetic structuring or how this changes over time in ephemeral patches. We predicted that species that specialize on ephemeral habitats will delay dispersal to exploit natal habitat patches while resources are plentiful and thus display fine-scale structure. To investigate this idea, we evaluated the effect of frequent colonization of ephemeral habitats on the fine-scale genetic structure of a fire specialist, the black-backed woodpecker (Picoides arcticus) and found a pattern of fine-scale genetic structure. We then tested for differences in spatial structure between sexes and detected a pattern consistent with male-biased dispersal. We also detected a temporal increase in relatedness among individuals within newly burned forest patches. Our results indicate that specialist species that outlive their ephemeral patches can accrue significant fine-scale spatial structure that does not necessarily affect spatial structure at larger scales. This highlights the importance of both spatial and temporal scale considerations in both sampling and data interpretation of molecular genetic results.

Habitat-linked population genetic differentiation in the blue tit Cyanistes caeruleus

Although the recent emergence of the field of landscape genetics has led to several studies investigating the effects of habitat composition between populations on genetic differentiation, much less is known on the impact of within-habitat ecological characteristics on levels of gene flow and genetic differentiation among populations. Using data on 840 individuals sampled in 8 sites in Corsica and 1 in southern France and analyzed at 10 microsatellite loci, we assessed the spatial and temporal population genetic structure of blue tits (Cyanistes caeruleus) in Corsica and identified the ecological factors, both between- and within-sampling sites, responsible for the observed genetic structure. We found temporally stable fine-scale genetic structure within Corsica, with genetic differentiation values among populations corresponding to ~25% of levels observed between Corsica and the mainland. This structure was not explained by the geographic distance among populations or by the presence of physical barriers but was instead related to local habitat types (deciduous or evergreen oaks). Our results are thus consistent with previously documented phenotypic differences among habitats in morphological and reproductive traits. These findings suggest that although individuals have high dispersal ability, local adaptation might reduce gene flow among populations located in different habitats.