Extra-pair paternity in two passerine birds breeding in a gradient of urbanisation.

Urbanisation has been increasing worldwide in recent decades, driving environmental change and exerting novel selective pressures on wildlife. Phenotypic differences between urban and rural individuals have been widely documented in several taxa. However, the extent to which urbanisation impacts mating strategies is less known. Here, we investigated extra-pair paternity variation in great tits (Parus major) and blue tits (Cyanistes caeruleus) breeding in nestboxes set in a gradient of urbanisation in Warsaw, Poland, over three breeding seasons. Urbanisation was quantified as the amount of light pollution, noise pollution, impervious surface area (ISA) and tree cover within a 100-m radius around each nestbox. We obtained genotypes for 1213 great tits at 7344 SNP markers and for 1299 blue tits at 9366 SNP markers with a genotyping-by-sequencing method, and inferred extra-pair paternity by computing a genomewide relatedness matrix. We report higher extra-pair paternity in blue tits breeding in more urbanised areas, for example, with higher light pollution and ISA, and lower tree cover. However, no such trend was found in great tits. Late-stage survival of individual nestlings in both species was not associated with paternity or urbanisation proxies, thus we were not able to detect fitness benefits or drawbacks of being an extra-pair offspring in relation to urbanisation. Our results contribute to the growing body of knowledge reporting on the effects of urbanisation on avian ecology and behaviour, and confirm species-specific and population-specific patterns of extra-pair paternity variation.

Projection of current and future distribution of adaptive genetic units in an alpine ungulate.

Climate projections predict major changes in alpine environments by the end of the 21st century. To avoid climate-induced maladaptation and extinction, many animal populations will either need to move to more suitable habitats or adapt in situ to novel conditions. Since populations of a species exhibit genetic variation related to local adaptation, it is important to incorporate this variation into predictive models to help assess the ability of the species to survive climate change. Here, we evaluate how the adaptive genetic variation of a mountain ungulate-the Northern chamois (Rupicapra rupicapra)-could be impacted by future global warming. Based on genotype-environment association analyses of 429 chamois using a ddRAD sequencing approach, we identified genetic variation associated with climatic gradients across the European Alps. We then delineated adaptive genetic units and projected the optimal distribution of these adaptive groups in the future. Our results suggest the presence of local adaptation to climate in Northern chamois with similar genetic adaptive responses in geographically distant but climatically similar populations. Furthermore, our results predict that future climatic changes will modify the Northern chamois adaptive landscape considerably, with various degrees of maladaptation risk.

Linking genetic, morphological, and behavioural divergence between inland island and mainland deer mice.

The island syndrome hypothesis (ISH) stipulates that, as a result of local selection pressures and restricted gene flow, individuals from island populations should differ from individuals within mainland populations. Specifically, island populations are predicted to contain individuals that are larger, less aggressive, more sociable, and that invest more in their offspring. To date, tests of the ISH have mainly compared oceanic islands to continental sites, and rarely smaller spatial scales such as inland watersheds. Here, using a novel set of genome-wide SNP markers in wild deer mice (Peromyscus maniculatus) we conducted a genomic assessment of predictions underlying the ISH in an inland riverine island system: analysing island-mainland population structure, and quantifying heritability of phenotypes thought to underlie the ISH. We found clear genomic differentiation between the island and mainland populations and moderate to high marker-based heritability estimates for overall variation in traits previously found to differ in line with the ISH between mainland and island locations. FST outlier analyses highlighted 12 loci associated with differentiation between mainland and island populations. Together these results suggest that the island populations examined are on independent evolutionary trajectories, the traits considered have a genetic basis (rather than phenotypic variation being solely due to phenotypic plasticity). Coupled with the previous results showing significant phenotypic differentiation between the island and mainland groups in this system, this study suggests that the ISH can hold even on a small spatial scale.

Population genetics reveals divergent lineages and ongoing hybridization in a declining migratory fish species complex.

Deciphering the effects of historical and recent demographic processes responsible for the spatial patterns of genetic diversity and structure is a key objective in evolutionary and conservation biology. Using population genetic analyses, we investigated the demographic history, the contemporary genetic diversity and structure, and the occurrence of hybridization and introgression of two species of anadromous fish with contrasting life history strategies and which have undergone recent demographic declines, the allis shad (Alosa alosa) and the twaite shad (Alosa fallax). We genotyped 706 individuals from 20 rivers and 5 sites at sea in Southern Europe at thirteen microsatellite markers. Genetic structure between populations was lower for the nearly semelparous species A. alosa, which disperses greater distances compared to the iteroparous species, A. fallax. Individuals caught at sea were assigned at the river level for A. fallax and at the region level for A. alosa. Using an approximate Bayesian computation framework, we inferred that the most likely long term historical divergence scenario between both species and lineages involved historical separation followed by secondary contact accompanied by strong population size decline. Accordingly, we found evidence for contemporary hybridization and bidirectional introgression due to gene flow between both species and lineages. Moreover, our results support the existence of at least one distinct species in the Mediterrannean sea: A. agone in Golfe du Lion area, and another divergent lineage in Corsica. Overall, our results shed light on the interplay between historical and recent demographic processes and life history strategies in shaping population genetic diversity and structure of closely related species. The recent demographic decline of these species' populations and their hybridization should be carefully considered while implementing conservation programs.

Riverscape genetics in brook lamprey: genetic diversity is less influenced by river fragmentation than by gene flow with the anadromous ecotype.

Understanding the effect of human-induced landscape fragmentation on gene flow and evolutionary potential of wild populations has become a major concern. Here, we investigated the effect of riverscape fragmentation on patterns of genetic diversity in the freshwater resident European brook lamprey (Lampetra planeri) that has a low ability to pass obstacles to migration. We tested the hypotheses of (i) asymmetric gene flow following water current and (ii) an effect of gene flow with the closely related anadromous river lamprey (L. fluviatilis) ecotype on L. planeri genetic diversity. We genotyped 2472 individuals, including 225 L. fluviatilis, sampled from 81 sites upstream and downstream barriers to migration, in 29 western European rivers. Linear modelling revealed a strong positive relationship between genetic diversity and the distance from the river source, consistent with expected patterns of decreased gene flow into upstream populations. However, the presence of anthropogenic barriers had a moderate effect on spatial genetic structure. Accordingly, we found evidence for downstream-directed gene flow, supporting the hypothesis that barriers do not limit dispersal mediated by water flow. Downstream L. planeri populations in sympatry with L. fluviatilis displayed consistently higher genetic diversity. We conclude that genetic drift and slight downstream gene flow drive the genetic make-up of upstream L. planeri populations whereas gene flow between ecotypes maintains higher levels of genetic diversity in L. planeri populations sympatric with L. fluviatilis. We discuss the implications of these results for the design of conservation strategies of lamprey, and other freshwater organisms with several ecotypes, in fragmented dendritic river networks.

Surface temperatures of non-incubated eggs in great tits (Parus major) are strongly associated with ambient temperature.

Temperature is one of the best investigated environmental factors in ecological life-history studies and is increasingly considered in the contexts of climate change and urbanization. In avian ecology, few studies have examined the associations between thermal dynamics in the nest environment and its neighbouring air. Here, we placed avian nests and non-incubated eggs inside nest boxes at various air temperatures that ranged from 0.3 to 33.1 °C, both in the field and in laboratory conditions. We measured how the design of the boxes, their compass orientation and their location in more or less urbanized environments affected the surface temperature of nests and eggs. We also assessed whether covering the eggs with lining material influenced their surface temperature. Overall, across all performed tests, we found that the surface temperature of nests and eggs strongly reflected the air temperature measured outside of the nest boxes. While the design of the nest boxes had little influence on the temperature of nests and eggs, orienting the nest boxes to the north or to the west significantly decreased their surface temperature. The presence of lining material also kept eggs slightly warmer when air temperatures were low. Altogether these results suggest that non-incubated eggs are not well protected against extreme air temperatures prior to the onset of incubation. From an evolutionary point of view, producers of ectotherm eggs need therefore to time egg-laying appropriately in order to avoid unfavourable thermal nest environments.

Assortative Mating in Animals and Its Role for Speciation.

Evolutionary theory predicts that positive assortative mating-the tendency of similar individuals to mate with each other-plays a key role for speciation by generating reproductive isolation between diverging populations. However, comprehensive tests for an effect of assortative mating on species richness at the macroevolutionary scale are lacking. We used a meta-analytic approach to test the hypothesis that the strength of assortative mating within populations is positively related to species richness across a broad range of animal taxa. Specifically, we ran a phylogenetically independent meta-analysis using an extensive database of 1,447 effect sizes for the strength of assortative mating, encompassing 307 species from 130 families and 14 classes. Our results suggest that there is no relationship between the strength of assortative mating and species richness across and within major taxonomic groups and trait categories. Moreover, our analysis confirms an earlier finding that animals typically mate assortatively (global Pearson correlation coefficient: r=0.36; 95% confidence interval: 0.19-0.52) when accounting for phylogenetic nonindependence. We argue that future advances will rely on a better understanding of the evolutionary causes and consequences of the observed intra- and interspecific variation in the strength of assortative mating.

No evidence of inbreeding depression in fast declining herds of migratory caribou.

Identifying inbreeding depression early in small and declining populations is essential for management and conservation decisions. Correlations between heterozygosity and fitness (HFCs) provide a way to identify inbreeding depression without prior knowledge of kinship among individuals. In Northern Quebec and Labrador, the size of two herds of migratory caribou (Rivière-George, RG and Rivière-aux-Feuilles, RAF) has declined by one to two orders of magnitude in the last three decades. This raises the question of a possible increase in inbreeding depression originating from, and possibly contributing to, the demographic decline in those populations. Here, we tested for the association of genomic inbreeding indices (estimated with 22,073 SNPs) with body mass and survival in 400 caribou sampled in RG and RAF herds between 1996 and 2016. We found no association of individual heterozygosity or inbreeding coefficient with body mass or annual survival. Furthermore, those genomic inbreeding indices remained stable over the period monitored. These results suggest that the rapid and intense demographic decline of the herds did not cause inbreeding depression in those populations. Although we found no evidence for HFCs, if demographic decline continues, it is possible that such inbreeding depression would be triggered.

Gene flow does not prevent personality and morphological differentiation between two blue tit populations.

Understanding the causes and consequences of population phenotypic divergence is a central goal in ecology and evolution. Phenotypic divergence among populations can result from genetic divergence, phenotypic plasticity or a combination of the two. However, few studies have deciphered these mechanisms for populations geographically close and connected by gene flow, especially in the case of personality traits. In this study, we used a common garden experiment to explore the genetic basis of the phenotypic divergence observed between two blue tit (Cyanistes caeruleus) populations inhabiting contrasting habitats separated by 25 km, for two personality traits (exploration speed and handling aggression), one physiological trait (heart rate during restraint) and two morphological traits (tarsus length and body mass). Blue tit nestlings were removed from their population and raised in a common garden for up to 5 years. We then compared adult phenotypes between the two populations, as well as trait-specific Qst and Fst . Our results revealed differences between populations similar to those found in the wild, suggesting a genetic divergence for all traits. Qst -Fst comparisons revealed that the trait divergences likely result from dissimilar selection patterns rather than from genetic drift. Our study is one of the first to report a Qst -Fst comparison for personality traits and adds to the growing body of evidence that population genetic divergence is possible at a small scale for a variety of traits including behavioural traits.

Do genetic drift and accumulation of deleterious mutations preclude adaptation? Empirical investigation using RADseq in a northern lacustrine fish.

Understanding genomic signatures of divergent selection underlying long-term adaptation in populations located in heterogeneous environments is a key goal in evolutionary biology. In this study, we investigated neutral, adaptive and deleterious genetic variation using 7,192 SNPs in 31 Lake Trout (Salvelinus namaycush) populations (n = 673) from Québec, Canada. Average genetic diversity was low, weakly shared among lakes, and positively correlated with lake size, indicating a major role for genetic drift subsequent to lake isolation. Putatively deleterious mutations were on average at lower frequencies than the other SNPs, and their abundance relative to the entire polymorphism in each population was positively correlated with inbreeding, suggesting that the effectiveness of purifying selection was negatively correlated with inbreeding, as predicted from theory. Despite evidence for pronounced genetic drift and inbreeding, several outlier loci were associated with temperature and found in or close to genes with biologically relevant functions notably related to heat stress and immune responses. Outcomes of gene-temperature associations were influenced by the inclusion of the most inbred populations, in which allele frequencies deviated the most from model predictions. This result illustrates challenge in identifying gene-environment associations in cases of high genetic drift and restricted gene flow and suggests limited adaptation in populations experiencing higher inbreeding. We discuss the relevance of these findings for the conservation and management, notably regarding stocking and genetic rescue, of Lake Trout populations and other species inhabiting highly fragmented habitats.

Inferring the demographic history underlying parallel genomic divergence among pairs of parasitic and nonparasitic lamprey ecotypes.

Understanding the evolutionary mechanisms generating parallel genomic divergence patterns among replicate ecotype pairs remains an important challenge in speciation research. We investigated the genomic divergence between the anadromous parasitic river lamprey (Lampetra fluviatilis) and the freshwater-resident nonparasitic brook lamprey (Lampetra planeri) in nine population pairs displaying variable levels of geographic connectivity. We genotyped 338 individuals with RAD sequencing and inferred the demographic divergence history of each population pair using a diffusion approximation method. Divergence patterns in geographically connected population pairs were better explained by introgression after secondary contact, whereas disconnected population pairs have retained a signal of ancient migration. In all ecotype pairs, models accounting for differential introgression among loci outperformed homogeneous migration models. Generating neutral predictions from the inferred divergence scenarios to detect highly differentiated markers identified greater proportions of outliers in disconnected population pairs than in connected pairs. However, increased similarity in the most divergent genomic regions was found among connected ecotype pairs, indicating that gene flow was instrumental in generating parallelism at the molecular level. These results suggest that heterogeneous genomic differentiation and parallelism among replicate ecotype pairs have partly emerged through restricted introgression in genomic islands.

Investigating the extent of parallelism in morphological and genomic divergence among lake trout ecotypes in Lake Superior.

Understanding the emergence of species through the process of ecological speciation is a central question in evolutionary biology which also has implications for conservation and management. Lake trout (Salvelinus namaycush) is renowned for the occurrence of different ecotypes linked to resource and habitat use throughout North America. We aimed to unravel the fine genetic structure of the four lake trout ecotypes in Lake Superior. A total of 486 individuals from four sites were genotyped at 6822 filtered SNPs using RADseq technology. Our results revealed different extent of morphological and genetic differentiation within the different sites. Overall, genetic differentiation was weak but significant and was on average three times higher between sites (mean FST  = 0.016) than between ecotypes within sites (mean FST  = 0.005) indicating higher level of gene flow or a more recent shared ancestor between ecotypes within each site than between populations of the same ecotype. Evidence of divergent selection was also found between ecotypes and/or in association with morphological variation. Outlier loci found in genes related to lipid metabolism and visual acuity were of particular interest in this context of ecotypic divergence. However, we did not find clear indication of parallelism at the genomic level, despite the presence of phenotypic parallelism among some ecotypes from different sampling sites. Overall, the occurrence of different levels of both genomic and phenotypic differentiation between ecotypes within each site with several differentiated loci linked to relevant biological functions supports the presence of a continuum of divergence in lake trout.

Genetic and phenotypic changes in an Atlantic salmon population supplemented with non-local individuals: a longitudinal study over 21 years.

While introductions and supplementations using non-native and potentially domesticated individuals may have dramatic evolutionary effects on wild populations, few studies documented the evolution of genetic diversity and life-history traits in supplemented populations. Here, we investigated year-to-year changes from 1989 to 2009 in genetic admixture at 15 microsatellite loci and in phenotypic traits in an Atlantic salmon (Salmo salar) population stocked during the first decade of this period with two genetically and phenotypically distinct source populations. We detected a pattern of temporally increasing introgressive hybridization between the stocked population and both source populations. The proportion of fish returning to the river after a single winter at sea (versus several ones) was higher in fish assigned to the main source population than in local individuals. Moreover, during the first decade of the study, both single-sea-winter and multi-sea-winter (MSW) fish assigned to the main source population were smaller than local fish. During the second decade of the study, MSW fish defined as hybrids were lighter and smaller than fish from parental populations, suggesting outbreeding depression. Overall, this study suggests that supplementation with non-local individuals may alter not only the genetic diversity of wild populations but also life-history traits of adaptive significance.

RAD genotyping reveals fine-scale genetic structuring and provides powerful population assignment in a widely distributed marine species, the American lobster (Homarus americanus).

Deciphering genetic structure and inferring connectivity in marine species have been challenging due to weak genetic differentiation and limited resolution offered by traditional genotypic methods. The main goal of this study was to assess how a population genomics framework could help delineate the genetic structure of the American lobster (Homarus americanus) throughout much of the species' range and increase the assignment success of individuals to their location of origin. We genotyped 10 156 filtered SNPs using RAD sequencing to delineate genetic structure and perform population assignment for 586 American lobsters collected in 17 locations distributed across a large portion of the species' natural distribution range. Our results revealed the existence of a hierarchical genetic structure, first separating lobsters from the northern and southern part of the range (FCT  = 0.0011; P-value = 0.0002) and then revealing a total of 11 genetically distinguishable populations (mean FST  = 0.00185; CI: 0.0007-0.0021, P-value < 0.0002), providing strong evidence for weak, albeit fine-scale population structuring within each region. A resampling procedure showed that assignment success was highest with a subset of 3000 SNPs having the highest FST . Applying Anderson's (Molecular Ecology Resources, 2010, 10, 701) method to avoid 'high-grading bias', 94.2% and 80.8% of individuals were correctly assigned to their region and location of origin, respectively. Lastly, we showed that assignment success was positively associated with sample size. These results demonstrate that using a large number of SNPs improves fine-scale population structure delineation and population assignment success in a context of weak genetic structure. We discuss the implications of these findings for the conservation and management of highly connected marine species, particularly regarding the geographic scale of demographic independence.

Stable isotope analyses on archived fish scales reveal the long-term effect of nitrogen loads on carbon cycling in rivers.

Stable isotope analysis of organic matter in sediment records has long been used to track historical changes in productivity and carbon cycling in marine and lacustrine ecosystems. While flow dynamics preclude stratigraphic measurements of riverine sediments, such retrospective analysis is important for understanding biogeochemical cycling in running waters. Unique collections of riverine fish scales were used to analyse δ(15) N and δ(13) C variations in the food web of two European rivers that experience different degrees of anthropogenic pressure. Over the past four decades, dissolved inorganic N loading remained low and constant in the Teno River (70°N, Finland); in contrast, N loading increased fourfold in the Scorff River (47°N, France) over the same period. Archived scales of Atlantic salmon parr, a riverine life-stage that feeds on aquatic invertebrates, revealed high δ(15) N values in the Scorff River reflecting anthropogenic N inputs to that riverine environment. A strong correlation between dissolved inorganic N loads and δ(13) C values in fish scales was observed in the Scorff River, whereas no trend was found in the Teno River. This result suggests that anthropogenic N-nutrients enhanced atmospheric C uptake by primary producers and its transfer to fish. Our results illustrate for the first time that, as for lakes and marine ecosystems, historical changes in anthropogenic N loading can affect C cycling in riverine food webs, and confirm the long-term interactions between N and C biogeochemical cycles in running waters.

Parallel and nonparallel genome-wide divergence among replicate population pairs of freshwater and anadromous Atlantic salmon.

Little is known about the genetic basis differentiating resident and anadromous forms found in many salmonid species. Using a medium-density SNP array, we documented genomic diversity and divergence at 2336 genetically mapped loci among three pairs of North American anadromous and freshwater Atlantic salmon populations. Our results show that across the genome, freshwater populations have lower diversity and a smaller proportion of private polymorphism relative to anadromous populations. Moreover, differentiation was more pronounced among freshwater than among anadromous populations at multiple spatial scales, suggesting a large effect of genetic drift in these isolated freshwater populations. Using nonhierarchical and hierarchical genome scans, we identified hundreds of markers spread across the genome that are potentially under divergent selection between anadromous and freshwater populations, but few outlier loci were repeatedly found in all three freshwater­anadromous comparisons. Similarly, a sliding window analysis revealed numerous regions of high divergence that were nonparallel among the three comparisons. These last results show little evidence for the parallel evolution of alleles selected for in freshwater populations, but suggest nonparallel adaptive divergence at many loci of small effects distributed through the genome. Overall, this study emphasizes the important role of genetic drift in driving genome-wide reduction in diversity and divergence in freshwater Atlantic salmon populations and suggests a complex multigenic basis of adaptation to resident and anadromous strategies with little parallelism.

Extra-pair paternity, breeding density, and synchrony in natural cavities versus nestboxes in two passerine birds.

Most of what is known about extra-pair paternity in hole-nesting birds derives from studies using artificial nesting sites, such as nestboxes. However, it has rarely been investigated whether inference drawn from breeding events taking place in nestboxes matches what would be observed under natural conditions, that is, in natural cavities. We here report on a variation in promiscuity in blue tits and great tits nesting in natural cavities and nestboxes in an urban forest in Warsaw, Poland. Specifically, we tested whether local breeding density, local breeding synchrony, and extra-pair paternity (inferred from SNP data generated with a high-throughput genotyping by sequencing method) differed between birds nesting in natural cavities and nestboxes. In both blue tits and great tits, the frequency of extra-pair paternity was similar between the two cavity types. In blue tits, we observed shorter nearest neighbor distance, higher neighbor density, and higher synchronous neighbor density (i.e., density of fertile females) in nestboxes relative to natural cavities. No such pattern was found in great tits. Moreover, we detected a positive relationship between the proportion of extra-pair offspring in the nest and neighbor density around the nest in blue tits. Our results revealed that the provisioning of nestboxes did not change rates of extra-pair paternity, suggesting that conclusions drawn from nestbox studies might adequately represent the natural variation in extra-pair matings in some species or sites. However, the observed differences in spatiotemporal components of breeding dynamics highlight the fact that these parameters should be carefully considered when comparing mating behavior across studies and/or sites.

Epigenetics and the city: Non-parallel DNA methylation modifications across pairs of urban-forest Great tit populations.

Identifying the molecular mechanisms involved in rapid adaptation to novel environments and determining their predictability are central questions in evolutionary biology and pressing issues due to rapid global changes. Complementary to genetic responses to selection, faster epigenetic variations such as modifications of DNA methylation may play a substantial role in rapid adaptation. In the context of rampant urbanization, joint examinations of genomic and epigenomic mechanisms are still lacking. Here, we investigated genomic (SNP) and epigenomic (CpG methylation) responses to urban life in a passerine bird, the Great tit (Parus major). To test whether urban evolution is predictable (i.e. parallel) or involves mostly nonparallel molecular processes among cities, we analysed both SNP and CpG methylation variations across three distinct pairs of city and forest Great tit populations in Europe. Our analyses reveal a polygenic response to urban life, with both many genes putatively under weak divergent selection and multiple differentially methylated regions (DMRs) between forest and city great tits. DMRs mainly overlapped transcription start sites and promotor regions, suggesting their importance in modulating gene expression. Both genomic and epigenomic outliers were found in genomic regions enriched for genes with biological functions related to the nervous system, immunity, or behavioural, hormonal and stress responses. Interestingly, comparisons across the three pairs of city-forest populations suggested little parallelism in both genetic and epigenetic responses. Our results confirm, at both the genetic and epigenetic levels, hypotheses of polygenic and largely nonparallel mechanisms of rapid adaptation in novel environments such as urbanized areas.

Evolving spatial conservation prioritization with intraspecific genetic data.

Spatial conservation prioritization (SCP) is a planning framework used to identify new conservation areas on the basis of the spatial distribution of species, ecosystems, and their services to human societies. The ongoing accumulation of intraspecific genetic data on a variety of species offers a way to gain knowledge of intraspecific genetic diversity and to estimate several population characteristics useful in conservation, such as dispersal and population size. Here, we review how intraspecific genetic data have been integrated into SCP and highlight their potential for identifying conservation area networks that represent intraspecific genetic diversity comprehensively and that ensure the long-term persistence of biodiversity in the face of global change.

Determinants of hierarchical genetic structure in Atlantic salmon populations: environmental factors vs. anthropogenic influences.

Disentangling the effects of natural environmental features and anthropogenic factors on the genetic structure of endangered populations is an important challenge for conservation biology. Here, we investigated the combined influences of major environmental features and stocking with non-native fish on the genetic structure and local adaptation of Atlantic salmon (Salmo salar) populations. We used 17 microsatellite loci to genotype 975 individuals originating from 34 French rivers. Bayesian analyses revealed a hierarchical genetic structure into five geographically distinct clusters. Coastal distance, geological substrate and river length were strong predictors of population structure. Gene flow was higher among rivers with similar geologies, suggesting local adaptation to geological substrate. The effect of river length was mainly owing to one highly differentiated population that has the farthest spawning grounds off the river mouth (up to 900km) and the largest fish, suggesting local adaptation to river length. We detected high levels of admixture in stocked populations but also in neighbouring ones, implying large-scale impacts of stocking through dispersal of non-native individuals. However, we found relatively few admixed individuals suggesting a lower fitness of stocked fish and/or some reproductive isolation between wild and stocked individuals. When excluding stocked populations, genetic structure increased as did its correlation with environmental factors. This study overall indicates that geological substrate and river length are major environmental factors influencing gene flow and potential local adaptation among Atlantic salmon populations but that stocking with non-native individuals may ultimately disrupt these natural patterns of gene flow among locally adapted populations.