Email updates

Keep up to date with the latest news and content from BMC Ecology and BioMed Central.

Open Access Highly Accessed Research article

Population genetic structure of gray wolves (Canis lupus) in a marine archipelago suggests island-mainland differentiation consistent with dietary niche

Astrid V Stronen12*, Erin L Navid3, Michael S Quinn4, Paul C Paquet56, Heather M Bryan567 and Christopher T Darimont567*

Author Affiliations

1 Mammal Research Institute, Polish Academy of Sciences, ul. Waszkiewicza 1, Białowieża 17-230, Poland

2 Department of Biotechnology, Chemistry and Environmental Engineering, Aalborg University, Sohngaardsholmsvej 57, Aalborg 9000, Denmark

3 Faculty of Environmental Design, University of Calgary, 2500 University Dr. NW, Calgary, Alberta T2N 1 N4, Canada

4 Institute for Environmental Sustainability, Mount Royal University, 4825 Mount Royal Gate SW, Calgary, Alberta T3E 6 K6, Canada

5 Department of Geography, University of Victoria, PO Box 3060 STN CSC, Victoria, British Columbia V8W 3R4, Canada

6 Raincoast Conservation Foundation, PO Box 86 Denny Island, British Columbia V0T 1B0, Canada

7 Hakai Beach Institute, Box 309, Heriot Bay, British Columbia V0P 1H0, Canada

For all author emails, please log on.

BMC Ecology 2014, 14:11  doi:10.1186/1472-6785-14-11


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1472-6785/14/11


Received:25 September 2013
Accepted:15 April 2014
Published:10 June 2014

© 2014 Stronen et al.; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. 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.

Abstract

Background

Emerging evidence suggests that ecological heterogeneity across space can influence the genetic structure of populations, including that of long-distance dispersers such as large carnivores. On the central coast of British Columbia, Canada, wolf (Canis lupus L., 1758) dietary niche and parasite prevalence data indicate strong ecological divergence between marine-oriented wolves inhabiting islands and individuals on the coastal mainland that interact primarily with terrestrial prey. Local holders of traditional ecological knowledge, who distinguish between mainland and island wolf forms, also informed our hypothesis that genetic differentiation might occur between wolves from these adjacent environments.

Results

We used microsatellite genetic markers to examine data obtained from wolf faecal samples. Our results from 116 individuals suggest the presence of a genetic cline between mainland and island wolves. This pattern occurs despite field observations that individuals easily traverse the 30 km wide study area and swim up to 13 km among landmasses in the region.

Conclusions

Natal habitat-biased dispersal (i.e., the preference for dispersal into familiar ecological environments) might contribute to genetic differentiation. Accordingly, this working hypothesis presents an exciting avenue for future research where marine resources or other components of ecological heterogeneity are present.

Keywords:
Canis lupus; Ecological divergence; Marine resources; Niche; Population genetic structure; Traditional ecological knowledge; Wolf

Background

Recent evidence indicates that ecological and environmental variation can result in genetic differentiation within many taxa, including highly mobile species. Examples include sea turtles (reviewed in Bowen and Karl [1]), fish species such as herring (Clupea harengus L., 1758; André et al. [2]) and hake (Merluccius merluccius L., 1758; Milano et al. [3]), and mammal species including the orca (Orcinus orca L., 1758; Hoelzel et al. [4]), cougar (Puma concolor L., 1771; McRae et al., [5]), lynx (Lynx canadensis Kerr, 1792; Rueness et al., [6]), coyote (Canis latrans Say, 1823; Sacks et al. [7]), and wolves (C. lupus L., 1758; Musiani et al. [8]; Pilot et al. [9]; Weckworth et al. [10]–[12]). For example, Muñoz-Fuentes et al. [13] showed strong genetic divergence over distances less than 500 km between wolves of coastal and interior regions of British Columbia (BC), Canada. Ecological and environmental dimensions such as climate and prey availability between areas, not distance, best explained population structure. These patterns arise because individuals may be more likely to survive and reproduce within their natal habitats (Davis and Stamps [14], Nosil et al. [15], Edelaar et al. [16]), which, in turn, can influence population genetic structure. A prediction from this body of work is that genetic divergence might be detected even over short geographical distances, and for highly mobile animals, should there be a sharp gradient in environmental conditions.

Such sharp ecological transitions occur between mainland and adjacent island environments within coastal BC. Although distances between mainland and neighbouring islands are small (<1500 m), the environments have striking geological and ecological differences. The mainland is topographically rugged, contains less shoreline for a given area and is relatively species-rich. In contrast, the neighbouring islands are less mountainous, have more complex shorelines, and host fewer species; notably absent are grizzly bears, (Ursus acrtos horribilis Ord, 1815), which compete with wolves for marine resources (Darimont and Paquet [17]; Paquet et al. [18]). Owing to these different environments, analyses of faeces and stable isotope data have identified distinctly different realized niches. Wolves from island populations rely on marine resources for up to 85% of their diet, whereas mainland conspecifics rarely include more than 30% (Darimont et al. [19,20]). Additionally, the coastal mainland supports moose (Alces alces L., 1758) and mountain goats (Oreamnus americanus Blainville, 1816) that are absent or rare on coastal islands. Consequently, these major prey items are commonly detected in wolf diet in mainland areas and only very rarely on islands (Darimont et al. [21]). Moreover, likely reflecting these distinct habitat and dietary niches, parasite prevalence also differs between areas; there is higher faecal prevalence of Giardia sp. infections on islands and a lower prevalence of Diphyllobothrium sp. relative to mainland sites (Bryan et al. [22]).

Our objective was to examine genetic data from wolves of coastal BC over a limited geographic area (~2000 km2, with a generally east–west mainland-island axis of <30 km) to test the hypothesis that ecological heterogeneity can drive population genetic structure of a highly mobile animal within a small area. We note that this prediction was also informed by holders of traditional ecological knowledge (TEK) in the Heiltsuk First Nation area, who distinguish between mainland “timber wolf” and island “coastal wolf” forms. Given these scholarly- and TEK-informed hypotheses and the sharp environmental gradients on the BC coast, we expected mainland-island genetic differentiation that mirrors ecological differences among neighbouring social groups.

Methods

Study area

The central coast of BC is a remote network of islands and naturally fragmented mainland landmasses with limited (but increasing) industrial anthropogenic disturbance. The area is characterized by a wet and temperate climate, and annual precipitation typically exceeds 350 cm (Darimont and Paquet [17]). A core area (~2000 km2) centered on Bella Bella (52°10’ N, 128° 09’ W) served as the location for this study (Additional file 1). This landscape is surrounded by ocean, which separates a mainland landmass (823 km2) and five main islands ranging in size from 150–250 km2. Distances from island to mainland range from 250 m to 1450 m. Observational and genetic data (Darimont et al. [19]; Navid [23]) suggest that wolf packs, defined by the multi-year association of genetically and morphologically distinct individuals, have either island or mainland home ranges. However, one group (Yeo-Coldwell [YC]) primarily uses island habitat but also a portion of the adjacent mainland. Other units are either mainland groups (Upper Roscoe [UR], Lower Roscoe [LR] or island groups (Cunningham-Chatfield [CC], Denny-Campbell [DC]). Moreover, wolves are commonly observed swimming among landmasses, and home ranges of social groups often include multiple islands or mainland landmasses (e.g. peninsulas; Paquet et al. [18]; Darimont [24]; McAllister and Darimont [25]).

Additional file 1. Map of the study area on the central coast of British Columbia, Canada. Shown are estimated home ranges of five wolf (Canis lupus) social groups.

Format: DOC Size: 3.6MB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

Sampling

One thousand and seventy-four (1074) wolf faecal samples were collected between winter 2003 and winter 2004. We collected the following number of samples per season: spring: n = 416 summer: n = 297 fall: n = 292 winter: n = 69. Sampling areas included wildlife trails, logging roads, and electrical power rights-of-way. We preserved each sample in a 50-ml Falcon tube with 95% ethanol. We selected samples for genetic analysis based on characteristics of the samples and collection sites that best predicted amplification success (minimal physical decay, high moisture content, canopy cover; Navid [23]). We extracted DNA from faecal samples with Qiagen QIAamp® DNA Stool Mini Kits and the ‘Protocol for isolation of DNA from larger amounts of stool’ (QIAamp® DNA Stool Mini Kit handbook, http://www.qiagen.com/literature/ webcite). We performed DNA extractions in a room physically separated from amplified PCR products and used exclusively for this study to reduce the risk of contamination. Final purified extracts were refrigerated at +4°C until use.

Microsatellite amplification

We amplified a panel of 14 microsatellite markers (13 autosomal and one Y chromosome marker). These were FH2001, FH2010, FH2017, FH2054, FH2088, FH2096, FH2422 (Breen et al. [26]), FH3313, FH3725 (Guyon et al. [27]), PEZ06, PEZ08, PEZ15, PEZ19 (Halverson J. in Neff et al. [28]), and the Y-chromosome marker MS41B (Sundquist et al. [29]). We genotyped n = 477 faecal samples. Polymerase chain reaction (PCR) conditions optimized for the markers, based on the Qiagen multiplexing kit, were: initial denaturation at 95°C for 15 minutes, then 35 cycles of denaturation at 94°C for 30 sec, annealing at 58°C for 90 sec, extension at 72°C for 60 sec, with final extension at 60°C for 30 min. Organisation of markers into multiplexes is shown in Additional file 2. Amplified PCR product was loaded into a 6.5% denaturing polyacrylamide gel, and run on a LICOR4300s DNA analyzer. Genotyping was done with a LICOR’s SAGA GT version 3.3 microsatellite analysis software.

Additional file 2. Multiplex combinations of 14 microsatellite markers for genetic analyses of wolves from the central coast of British Columbia, Canada.

Format: DOC Size: 29KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

We accepted for further analyses samples that amplified at least 9/14 loci, and used the Excel Microsatellite Toolkit (Park [30]) to test for the presence of matching profiles. We consolidated matches (i.e., profiles with ≥ 75% matching alleles to account for uncertainty in genotyping) into one profile and retained the profile with the highest amplification rate. Matches were tested across all samples, but only observed within wolf groups and 33 profiles were removed.

N = 116 individual profiles were identified. Only 28 wolves (24%) were identified as males, and we observed six alleles for marker MS41B (209, 211, 213, 217, 219, 221). We used MICRO-CHECKER 2.2.3 (van Oosterhout et al. [31]) to assess possibilities of null alleles, large allele dropout, and scoring errors due to stutter peaks. We repeated genotyping of 50 samples collected during the fall season (deemed to represent 50 different individuals from all four groups [CC, DC, LR, YC] based on the abovementioned criteria) to evaluate data quality and estimate genotyping error. Here we estimated per-locus error rates based on the percentage of loci that did not show the same result twice (Additional file 3). Loci for which we obtained the same results twice were accepted as duplicated loci. From these results we obtained duplicate genotypes comprising five or more loci for 18 individuals (i.e. every locus in each of these 18 genotypes provided consistent results when re-tested). Based on amplification and error rates (Additional file 3) we removed MS41B, PEZ08, FH2017, and FH3313 from further analyses.

Additional file 3. Calculation of amplification and error rates and assessment of null alleles, large allele dropout, and stutter peaks for wolf samples (n = 116) from the central coast of British Columbia, Canada.

Format: DOC Size: 44KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

Statistical analyses

We calculated allelic diversity and observed and expected heterozygosity (with correction for sample size bias; Nei [32]) per locus in GENEPOP 3.4 (Raymond and Rousset [33]) and Genetix 4.05.2 (Belkhir et al. [34]), and FIS according to Weir and Cockerham [35], for mainland and island wolves. We tested for departures from Hardy-Weinberg equilibrium per locus in GENEPOP 3.4 with the Markov chain method (Guo and Thompson [36]). The results were adjusted to account for multiple comparisons with the false discovery rate (FDR, Verhoeven et al. [37]). Subsequently, we performed centered and scaled principal component analyses (PCA) with the adegenet-package (Jombart [38]) in R 2.14.2 (R development Core Team [39]). The PCA approach does not assume genetic equilibrium conditions and is well-suited for identifying spatial patterns such as genetic clines (i.e., gradients rather than separate clusters or complete admixture) that can be difficult to detect (Jombart et al. [40]). We repeated the PCA with the 18 individuals for which we had duplicated genotypes (and thus higher confidence), and all wolves identified as males to confirm the presence of the observed cline. On average, male wolves may disperse longer distances than females and are more likely to join new packs. As our data set comprised relatively few confirmed males, we tested these results separately to check if the island-mainland gradient remained consistent. A high proportion of females in our sample might otherwise have contributed to the observed gradient if females disperse less frequently and/or shorter distances than males. We then performed a spatial PCA (henceforth sPCA; Jombart et al. [40]), which also takes spatial sampling information into account. As multiple samples were at times collected from the same location, we added 100 m of jitter (small amount of noise) to the UTM coordinates. We performed a spatial autocorrelation in GenAlEx (Peakall and Smouse [41] and references therein) to examine the possible existence of isolation-by-distance in our data set. We used distance classes of 5 km to obtain fine-scale results for our study area. Finally, we performed a partial Mantel-test in R with the Vegan package (Oksanen et al. [42]) to examine the relationship between genetic distance and island-mainland habitat type while controlling for geographic distance. This allowed us to test whether there was an effect of habitat type on fine-scale genetic structure after accounting for the effect of geographic distance. Geographic distance and habitat may be co-linear and their effects could be difficult to separate. Consequently, we also examined the relationship between geographic distance and habitat type. For these tests we incorporated co-dominant genotypic and Euclidean geographic distance matrices exported from GenAlEx and a third matrix with island-mainland habitat designations. We used Pearson’s correlation coefficient with n = 999 permutations.

Results

The average number of alleles per locus was 5.8 for mainland wolves and 6.8 for island wolves (Table 1). For mainland wolves, expected heterozygosity was 0.632 and five loci showed departures from Hardy-Weinberg equilibrium with observed levels of heterozygosity lower than expected. FIS results were positive for all loci with a mean value of 0.264. For island wolves, expected heterozygosity was 0.690 and seven loci showed departures from Hardy-Weinberg equilibrium (four of these were consistent between mainland and island wolves). FIS results for island wolves were positive for all except two loci, with a mean value of 0.211. We identified possible null alleles and stutter peaks for the overall sample, but dropout of large alleles was not detected (Additional file 3).

Table 1. Genetic diversity measures for wolves ( Canis lupus ) from the central coast of British Columbia, Canada

PCA results indicated the presence of a genetic cline between island and mainland wolves (Figure 1a, c). Although overlap was extensive, the results suggested an east–west gradient in profiles across < 30 km. Examination of genetic profiles based on the known wolf groups in the area (UR was not represented in the second analysis) suggested limited overlap between LR (Mainland) and YC (Island) wolves (Figure 1b, d). The CC and DC island groups occupy an intermediate position, along with the UR group from the mainland. Colour plots (Additional file 4) show the individual genetic profiles throughout the study area, and display a similar east–west gradient from the mainland to the islands. The PCA results for individuals identified as males (n = 28) were consistent with island-mainland differentiation (Additional file 5). For the sPCA, one global structure (and no local structure) was apparent (Additional file 6). When mapped across the geographic space, the global structure revealed an east–west gradient where YC and LR were the most differentiated groups (Figure 2). The partial Mantel test gave a correlation coefficient of 0.011 (p-value 0.351) between genetic distance and habitat matrices. The test between geographic distance and habitat matrices produced a correlation coefficient of 0.568 (p-value 0.001). Spatial autocorrelation results were positive for the first 17 km, negative from approximately 17–45 km, and subsequently positive (though this may be considered as zero autocorrelation at the larger distance classes with wide confidence intervals; Additional file 7).

thumbnailFigure 1. Principal component analyses of wolves (Canis lupus) from the central coast of British Columbia, Canada showing geographic distribution of individuals. a) Individual (n = 116) profiles based on ≥ 10 microsatellite loci labelled according to mainland (MA) and island (IS) sample locations. b) Individual profiles (n = 116) based on ≥ 10 microsatellite loci labeled according to membership in five wolf family groups: Upper Roscoe (UR) and Lower Roscoe (LR) on the mainland, and Yeo-Coldwell (YC), Cunningham-Chatfield (CC), and Denny-Campbell (DC) islands. Note that the label for DC (green colour) is overlapped by UR (red colour). c) A subsample of individual profiles (n = 18) with duplicated genotypes based on ≥ 5 loci labelled according to mainland and island sample locations. d) Individual profiles (n = 18) with duplicated genotypes based on ≥ 5 loci labelled according to membership in four wolf family groups LR, YC, CC, and DC (none from UR).

Additional file 4. Colour plot of wolf profiles from the central coast of British Columbia, Canada. a) Individual profiles (n = 116) based on ≥ 10 microsatellite loci. The first axis represents 6.1% of the variation, the second axis 5.1%. b) A subsample of individual profiles (n = 18) based on ≥ 5 duplicated loci. Genetic diversity is represented by distance and colour; individuals further apart and/or labelled with more dissimilar colours have more divergent genotypes. The first axis represents 18.8% of the variation, the second axis 13.9%.

Format: DOC Size: 50KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

Additional file 5. Principal component analysis (PCA) of male wolves from the central coast of British Columbia, Canada, showing island (IS, n = 19) and mainland (MA, n = 9) individuals. The first axis represents 11.9% of the variation, the second axis 9.7%. PCA is based on the 10 loci retained for final analyses (Table 1 and Additional file 2).

Format: DOC Size: 54KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

Additional file 6. Eigenvalues from a spatial principal component analysis (sPCA) on 10 microsatellite loci from 116 wolves from the central coast of British Columbia, Canada. Positive values (left side) represent global structures and negative values (right side) show local patterns. Tests for local and global structure revealed the presence of one global structure, which was subsequently interpreted.

Format: DOC Size: 39KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

thumbnailFigure 2. Spatial principal component analysis of wolves (Canis lupus) from the central coast of British Columbia, Canada, showing the first global structure mapped across the study area. Individual profiles (n = 116) are based on ≥ 10 microsatellite loci and originate from five wolf family groups: Upper Roscoe (UR) and Lower Roscoe (LR) on the mainland, and Yeo-Coldwell (YC), Cunningham-Chatfield (CC), and Denny-Campbell (DC) islands.

Additional file 7. Spatial autocorrelation analysis of wolf samples (n = 116) from the central coast of British Columbia, Canada, using 5 km distance classes. The Y axis shows the kinship coefficient (r), and U and L are the upper and lower limits for the 95% confidence interval of no spatial structure occurring in the data set after permutation (n = 999). Error bars show the 95% confidence interval around r as determined by bootstrap resampling (n = 999).

Format: DOC Size: 50KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

Discussion

Genetic variation

Allelic diversity and expected heterozygosity for island wolves (6.8, 0.690) and mainland wolves (5.8, 0.63) were relatively high and comparable to values reported for island populations of wolves on the Pacific Coast in southeast Alaska (5, 0.52; Weckworth et al. [10]) and coastal island populations in Arctic Canada (4.2, 0.61; Carmichael et al. [43]). Allelic diversity and expected heterozygosity were somewhat lower for wolves on the mainland portion of our study area, although this might, at least in part, reflect chance effects of our relatively small sample sizes. Comparison with FIS values from southeast Alaska islands (0.05) and coastal islands in Arctic Canada (0.181) suggest a higher degree of mating among relatives in mainland (0.264) and island (0.211) wolves from coastal BC. However, Carmichael et al. [43] also observed high FIS values on Victoria Island (0.427, n = 52) and on islands in the High Arctic (0.629, n = 11). Based on the findings from Alaska wolves, continental wolves appear to have higher genetic diversity. We would also expect a similar situation for our study area, as mainland wolves have a wider surrounding area from which to receive immigrants. However, there are known wolf groups on neighbouring islands not included in this study and we cannot exclude the possibility that immigration from these areas may have augmented the diversity in our sample of island wolves.

Non-invasive sampling and genotyping

Allelic dropout in non-invasive sampling (Santini et al. [44]) could, at least in part, explain the lower values for observed heterozygosity, the high number of loci not in Hardy-Weinberg equilibrium, and the positive FIS values. Our results could also have been influenced by the presence of null alleles. When most of the loci indicate null alleles, however, the MICRO-CHECKER program warns there may not be random mating in the population (panmixia). The PCA and sPCA findings of island-mainland differentiation suggest absence of panmixia in our study area. We therefore believe that island-mainland structure contributed to the frequent reports of null alleles. The study area is difficult to access, and many samples may have been several weeks old and thus affected by exposure to the humid climate (Santini et al. [44]; Navid [23]). Our results are based on analyses of faecal material, where duplicated genotypes were obtained for 15% (18 of 116) individuals. Error rates were high, but we do not expect any consistent bias between areas. Results from the duplicated genotypes accord with the larger dataset, although further sampling and multiple-tube analyses (e.g. Santini et al. [44]) would be necessary for accurate identification of individual wolves and to confirm dispersal events in our study area.

Evolutionary ecology and genetic differentiation between mainland and island wolves

The partial Mantel test showed no significant relationship between genetic distance and island-mainland habitat type when accounting for geographic distance. However, there was a significant correlation between geographic distance and island-mainland habitat type, suggesting that the two matrices are collinear and their effects cannot be differentiated. The spatial autocorrelation indicated negative autocorrelation from approximately 17–45 km. These results appear to contrast with those of Muñoz-Fuentes et al. [13] who reported that geographic distance was unlikely to explain the spatial structure of wolf mtDNA haplotypes in a broader study of coastal and central BC. Wolves are highly capable dispersers able to travel > 70 km/day (Mech and Boitani [45]), and it seems unlikely that geographic distance alone can explain the island-mainland structure suggested by the sPCA. In such a situation, we would expect the spatial autocorrelation results to show consistent (and increasing) negative kinship-values with geographic distance. In contrast, the 45–50 distance class that represents wolves in the northern- and southernmost parts of our study, which are farthest apart in geographic distance, showed positive values (or, more likely, no autocorrelation). Multiple interacting factors, including distance, water, terrain ruggedness etc., may affect genetic structure in our study area. Although it is essential to evaluate the possible influence of physiography on the differences observed between island and mainland wolves, the observed correlation between geographic distance and habitat type combined with the physical complexity of the landscape make it problematic and potentially misleading to use linear distances for estimating wolf movement.

Water barriers between the mainland and islands might restrict dispersal and gene flow. For example, captive wolves released on Coronation Island in Alaska did not swim 900 m to nearby habitat with abundant food (Klein [46]). In our study, we reject this hypothesis because we commonly observe wolves swimming among landmasses and distances among islands (including the multiple landmasses used by some groups) are often larger than the distances between islands and the mainland (Darimont et al., unpublished data). Immigrants from outside the study area could also influence the observed east–west gradient in genetic profiles. The differentiation seen in YC profiles, for example, may be explained by gene flow from unsampled wolves on the outer islands father west. Similarly, profiles from the UR group, which showed considerable overlap with island wolves, might result from immigration by one or more island wolves with high reproductive success. Furthermore, the presence of intermediate profiles in the north and south of our study area implies an island-mainland gradient. A strict island-mainland dichotomy may thus be simplistic and should be evaluated on a broader geographic scale. Without genetic data from a larger spatial extent, however, we cannot evaluate these hypotheses.

Family group structure might also have influenced our results, especially for long-lived animals for which the genetic influence of one successful breeder can be detected for many generations. Difficulties with amplification of MS41B likely reduced our ability to identify male wolves. A possible higher prevalence of females in the sample might nevertheless exacerbate genetic structuring in species where males are more likely to disperse. However, male wolf profiles we assessed showed a similar island-mainland gradient. Observational and tracking data suggest that wolf group size in the study area was ≤ 10 individuals (Darimont [24]), and it appears unlikely that the observed gradient in genetic profiles could be explained by social structure (i.e., wolf pack membership) alone. The identification of 116 individuals in the study area appears reasonable for a sampling period that included 2 litters, winter pup mortality that may exceed 50%, and the likelihood that 20% of individuals would be solitary or extra-territorial dispersers (Mech and Boitani [47]).

Despite the above-mentioned uncertainties, we offer the working hypothesis that the sharp ecological gradient between island and mainland locations, as revealed by the landscape characteristics and the dietary and parasitic data from wolves in our study area, can influence population genetic structure. Although our study must be interpreted with caution, and should be repeated with genetic profiles of higher quality, the results appear consistent with an increasing body of literature reporting genetic differentiation in wolves and other highly mobile species (see Introduction) influenced by ecological and environmental factors. Dispersal rates and gene flow might differ substantially between island and mainland sub-populations, and the extent to which populations are demographically independent could help define management units (Palsbøll et al. [48]) along the Pacific coast.

Associated morphological or other characteristics observed over time might have allowed TEK knowledge holders to recognize these dissimilar wolf forms. Such intra-specific nomenclature is common among indigenous knowledge holders (Turner et al. [49]). Indeed, in adjacent southeast Alaska, the frequency of the black colour phase among wolves killed by trappers is ~50% on the mainland and only ~20% on the islands (Person et al. [50]). Additional morphological differences among wolves of coastal BC might have led to mainland-island classification by local people.

The evolutionary influence of marine resources, which are pronounced on islands in our study area, can be dramatic for terrestrial wildlife. For example, polar bears (U. maritimus Phipps, 1774) are thought to have evolved from grizzly bears in peripheral areas where marine resources were abundant (Shields et al. [51]). Moreover, wolves of coastal BC (mainland and island populations) were thought to have diverged from interior populations in part because of marine resource availability in coastal zones (Muñoz-Fuentes et al. [13]). Individuals born in this distinct environment are likely better able to survive and reproduce within, compared to beyond, these conditions.

Conclusions

Our results indicate the presence of a genetic cline between island and mainland wolves. Although overlap was extensive, the results suggest an east–west gradient in profiles across < 30 km. We hypothesize that adaptive responses to heterogeneity in food resources can influence genetic differentiation. Accordingly, this line of inquiry presents an exciting avenue for future research where marine resources or other components of ecological heterogeneity are present.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

ELN carried out laboratory work, and contributed to fieldwork, statistical analyses, manuscript writing and revisions. AVS performed statistical analyses and manuscript writing. MSQ supervised the study and helped revise the manuscript. CTD participated in study design, fieldwork, and manuscript writing. PCP contributed to study design, field work and manuscript revisions. HMB created the maps and helped revise the manuscript. All authors read and approved the final manuscript.

Acknowledgements

We thank the Heiltsuk First Nation for their permission to study an important cultural animal in their traditional territory. We also thank J. Carpenter, J. Gordon-Walker, B. Mann, M. Musiani, M. Niedziałkowska, G. Pflueger, Y. Plante, and C. Service for field, laboratory and/or intellectual contributions. We are especially grateful to Chester Starr (Lone Wolf) of the Heiltsuk Nation for the knowledge he shared with us, which inspired this work. Raincoast Conservation Foundation, National Geographic Society, Wilburforce Foundation, and the Tula Foundation supported this work. ELN was supported by a Natural Science and Engineering Research Council Graduate Scholarship. CTD was supported by Natural Science and Engineering Research Council Discovery Grant no. 435683. J. Olano-Marin and B. Martinez-Cruz offered constructive comments on an earlier version of the text. We thank the editor and two anonymous reviewers for their assistance in improving our manuscript.

References

  1. Bowen BW, Karl SA: Population genetics and phylogeography of sea turtles.

    Mol Ecol 2007, 16:4886-4907. PubMed Abstract | Publisher Full Text OpenURL

  2. André C, Larsson LC, Laikre L, Bekkevold D, Brigham D, Carvalho GR, Dahlgren TG, Hutchinson WF, Mariani S, Mudde K, Ruzzante DE, Ryman N: Detecting population structure in a high gene-flow species, Atlantic herring (Clupea harengus): direct, simultaneous evaluation of neutral vs putatively selected loci.

    Heredity 2011, 106:270-280. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  3. Milano I, Babbucci M, Cariani A, Atanassova M, Bekkevold D, Carvalho GR, Espiñeira M, Fiorentino F, Garofalo G, Geffen AJ, Hansen JH, Helyar SJ, Nielsen EE, Ogden R, Patarnello T, Stagioni M, Tinti F, Bargelloni L, FishPopTrace Consortium: Outlier markers reveal fine-scale genetic structuring across European hake populations (Merluccius merluccius).

    Mol Ecol 2014, 23:118-135. PubMed Abstract | Publisher Full Text OpenURL

  4. Hoelzel AR, Hey J, Dahlheim M, Nicholson C, Burkanov V, Black N: Evolution of population structure in a highly social top predator, the killer whale.

    Mol Biol Evol 2007, 24:1407-1415. PubMed Abstract | Publisher Full Text OpenURL

  5. McRae BH, Beier P, Dewald LE, Huynh LY, Keim P: Habitat barriers limit gene flow and illuminate historical events in a wide-ranging carnivore, the American puma.

    Mol Ecol 2005, 14:1965-1977. PubMed Abstract | Publisher Full Text OpenURL

  6. Rueness EK, Stenseth NC, O’Donoghue M, Boutin S, Ellegren H, Jacobsen KS: Ecological and genetic spatial structuring in the Canadian lynx.

    Nature 2003, 425:69-72. PubMed Abstract | Publisher Full Text OpenURL

  7. Sacks BN, Brown SK, Ernest HB: Population structure of California coyotes corresponds to habitat-specific breaks and illuminates species history.

    Mol Ecol 2004, 13:1265-1275. PubMed Abstract | Publisher Full Text OpenURL

  8. Musiani M, Leonard JA, Cluff HD, Gates CC, Mariani S, Paquet PC, Vilà C, Wayne RK: Differentiation of tundra and boreal coniferous forest wolves: genetics, coat color and association with migratory caribou.

    Mol Ecol 2007, 16:4149-4170. PubMed Abstract | Publisher Full Text OpenURL

  9. Pilot M, Jędrzejewski W, Sidorovich VE, Meier-Augustein W, Hoelzel AR: Dietary differentiation and the evolution of population genetic structure in a highly mobile carnivore.

    PLoS ONE 2012, 7:e39341. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  10. Weckworth BV, Talbot S, Sage GK, Person DK, Cook J: A signal for independent coastal and continental histories among North American wolves.

    Mol Ecol 2005, 14:917-931. PubMed Abstract | Publisher Full Text OpenURL

  11. Weckworth BV, Talbot SL, Cook JA: Phylogeography of wolves (Canis lupus) in the Pacific Northwest.

    J Mammal 2010, 91:363-375. Publisher Full Text OpenURL

  12. Weckworth BV, Dawson NG, Talbot SL, Flamme MJ, Cook JA: Going Coastal: Shared evolutionary history between Coastal British Columbia and Southeast Alaska wolves (Canis lupus).

    PLoS ONE 2011, 6:e19582. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  13. Muñoz-Fuentes V, Darimont CT, Wayne RK, Paquet PC, Leonard JA: Ecological factors drive differentiation in wolves from British Columbia.

    J Biogeogr 2009, 36:1516-1531. Publisher Full Text OpenURL

  14. Davis JM, Stamps JA: The effect of natal experience on habitat preferences.

    Trends Ecol Evol 2004, 19:411-416. PubMed Abstract | Publisher Full Text OpenURL

  15. Nosil P, Vines TH, Funk DJ: Reproductive isolation caused by natural selection against immigrants from divergent habitats.

    Evolution 2005, 59:705-719. PubMed Abstract OpenURL

  16. Edelaar P, Siepielski AM, Clobert J: Matching habitat choice causes direct gene flow: a neglected dimension in evolution and ecology.

    Evolution 2008, 62:2462-2472. PubMed Abstract | Publisher Full Text OpenURL

  17. Darimont CT, Paquet PC: The gray wolves, Canis lupus, of British Columbia’s central and north coast: distribution and conservation assessment.

    Can Field-Nat 2002, 116:416-422. OpenURL

  18. Paquet PC, Alexander SM, Swan PL, Darimont CT: The influence of natural landscape fragmentation and resource availability on connectivity and distribution of marine gray wolf populations on the Central Coast, British Columbia, Canada. In Connectivity Conservation. Edited by Crooks K, Sanjayan MA. Cambridge: Cambridge University Press; 2006:130-156. OpenURL

  19. Darimont CT, Paquet PC, Reimchen TE: Spawning salmon disrupt trophic coupling between wolves and ungulate prey in coastal British Columbia.

    BMC Ecol 2008, 8:14. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  20. Darimont CT, Paquet PC, Reimchen TE: Landscape heterogeneity and marine subsidy generate extensive intrapopulation niche diversity in a large terrestrial vertebrate.

    J Anim Ecol 2009, 78:126-133. PubMed Abstract | Publisher Full Text OpenURL

  21. Darimont CT, Price MHH, Winchester NN, Gordon-Walker J, Paquet PC: Predators in natural fragments: foraging ecology of wolves in British Columbia’s central and north coast archipelago.

    J Biogeogr 2004, 31:1867-1877. Publisher Full Text OpenURL

  22. Bryan HM, Darimont CT, Hill JE, Paquet PC, Thompson RCA, Wagner B, Smits JE: Seasonal and biogeographical patterns of gastrointestinal parasites in large carnivores: wolves in a coastal archipelago.

    Parasitology 2012, 139:781-790. PubMed Abstract | Publisher Full Text OpenURL

  23. Navid EL: Examination of genotyping success and population genetic structure of wolves in coastal British Columbia using non-invasive molecular markers. University of Calgary, Calgary: Faculty of Environmental Design; 2009. [MSc thesis] OpenURL

  24. Darimont CT: Ecological and evolutionary causes and consequences of intra-population variability in foraging niche; predator and prey in a marine archipelago. University of Victoria, Victoria: Department of Biology; 2007. [PhD dissertation] OpenURL

  25. McAllister I, Darimont CT: The Last Wild Wolves. Berkeley: University of California Berkeley Press; 2007. [Ghosts of the Rainforest] OpenURL

  26. Breen M, Jouquand S, Renier C, Mellersh CS, Hitte C, Holmes NG, Chéron A, Suter N, Vignaux F, Bristow AE, Priat C, McCann E, André C, Boundy S, Gitsham P, Thomas R, Bridge WL, Spriggs HF, Ryder EJ, Curson A, Sampson J, Ostrander EA, Binns MM, Galibert F: Chromosome-specific single-locus FISH probes allow anchorage of an 1800-marker integrated radiation-hybrid/linkage map of the domestic dog genome to all chromosomes.

    Genome Res 2001, 11:1784-1795. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  27. Guyon R, Lorentzen TD, Hitte C, Kim L, Cadieu E, Parker HG, Quignon P, Lowe JK, Renier C, Gelfenbeyn B, Vignaux F, DeFrance HB, Gloux S, Mahairas GG, André C, Galibert F, Ostrander EA: A 1-Mb resolution radiation hybrid map of the canine genome.

    Proc Natl Acad Sci USA 2003, 100:5296-5301. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  28. Neff MW, Broman KW, Mellersh CS, Ray K, Acland GM, Aguirre GD, Ziegle JS, Ostrander EA, Rine J: A second-generation genetic linkage map of the domestic dog, Canis familiaris.

    Genetics 1999, 151:803-820. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  29. Sundquist AK, Ellegren H, Olivier M, Vilà C: Y-chromosome haplotyping in Scandinavian wolves (Canis lupus) based on microsatellite markers.

    Mol Ecol 2001, 10:1959-1966. PubMed Abstract | Publisher Full Text OpenURL

  30. Park SDE: Trypanotolerance in West African cattle and the population genetic effects of selection. University of Dublin, Dublin: Smurfit Institute of Genetics; 2001. [PhD dissertation] OpenURL

  31. van Oosterhout C, Hutchinson WF, Wills DPM, Shipley P: MICRO-CHECKER: software for identifying and correcting genotyping errors in microsatellite data.

    Mol Ecol Notes 2004, 4:535-538. Publisher Full Text OpenURL

  32. Nei M: Estimation of average heterozygosity and genetic distance from a small number of individuals.

    Genetics 1978, 89:583-590. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  33. Raymond M, Rousset F: Genepop (version 1.2): Population genetics software for exact tests and ecuminiscism.

    J Hered 1995, 86:248-249. OpenURL

  34. Belkhir K, Borsa P, Chikhi L, Raufaste N, Bonhomme F: GENETIX 4.05, logiciel sous Windows TM pour la génétique des populations. Montpellier (France): Laboratoire Génome, Populations, Interactions, CNRS UMR 5000, Université de Montpellier II; 2004. OpenURL

  35. Weir BS, Cockerham CC: Estimating F-statistics for the analyses of population structure.

    Evolution 1984, 38:1358-1370. Publisher Full Text OpenURL

  36. Guo SW, Thompson EA: Performing the exact test of Hardy-Weinberg proportion for multiple alleles.

    Biometrics 1992, 48:361-372. PubMed Abstract | Publisher Full Text OpenURL

  37. Verhoeven KJF, Simonsen KL, McIntyre LM: Implementing false discovery rate control: Increasing your power.

    Oikos 2005, 108:643-647. Publisher Full Text OpenURL

  38. Jombart T: adegenet: a R package for the multivariate analysis of genetic markers.

    Biochem Genet 2008, 10:149-163. OpenURL

  39. R Development Core Team: R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2012.

    ISBN 3-900051-07-0. [http://www.R-project.org/ webcite]

    OpenURL

  40. Jombart T, Devillard S, Durfour A-B: Revealing cryptic spatial patterns in genetic variability by a new multivariate method.

    Heredity 2008, 101:92-103. PubMed Abstract | Publisher Full Text OpenURL

  41. Peakall R, Smouse PE: GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research – an update.

    Bioinformatics 2012, 28:2537-2539. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  42. Oksanen J, Blanchet FG, Kindt R, Legendre P, Minchin PR, O’Hara RB, Simpson GL, Solymos P, Henry M, Stevens H, Wagner H: Vegan: Community Ecology Package.

    2013.

    R-package version 2.0-10. http://CRAN.R-project.org/package=vegan webcite

  43. Carmichael LE, Krizan J, Nagy JA, Dumond M, Johnson D, Veitch A, Strobeck C: Northwest passages: conservation genetics of Arctic Island wolves.

    Conserv Genet 2008, 9:879-892. Publisher Full Text OpenURL

  44. Santini A, Lucchini V, Fabbri E, Randi E: Ageing and environmental factors affect PCR success in wolf (Canis lupus) extremental DNA samples.

    Mol Ecol Notes 2007, 7:955-961. Publisher Full Text OpenURL

  45. Mech LD, Boitani L: Introduction. In Wolves: Behaviour, Ecology and Conservation . Edited by Mech LD, Boitani L. Chicago: University of Chicago Press; 2003.

    xv-xvii

    OpenURL

  46. Klein DR: The introduction, increase, and demise of wolves on Coronation Island, Alaska. In Ecology and conservation of wolves in a changing world. Edited by Carbyn LN, Fritts SH, Seip DR. Edmonton: Canadian Circumpolar Institute, Occasional Publication No. 35; 1995:275-280. OpenURL

  47. Mech LD, Boitani L: Wolf social ecology. In Wolves: Behaviour, Ecology and Conservation . Edited by Mech LD, Boitani L. Chicago: University of Chicago Press; 2003:1-34. OpenURL

  48. Palsbøll PJ, Bérubé M, Allendorf FW: Identification of management units using population genetic data.

    Trends Ecol Evol 2006, 22:11-16. PubMed Abstract | Publisher Full Text OpenURL

  49. Turner NJ, Ignace MB, Ignace R: Traditional ecological knowledge and wisdom of aboriginal peoples in British Columbia.

    Ecol Appl 2000, 10:1275-1287. Publisher Full Text OpenURL

  50. Person DK, Kirchhoff M, Van Ballenberghe V, Iverson GC, Grossman E: The Alexander Archipelago wolf: a conservation assessment. Portland, OR: U.S: Department of Agriculture, Forest Service, Pacific Research Station; 1996. [Gen. Tech. Rep. PNW-GTR-384] OpenURL

  51. Shields GF, Adams D, Garner G, Labelle M, Pietsch J, Ramsay M, Schwartz C, Titus K, Williamson S: Phylogeography of Mitochondrial DNA Variation in Brown Bears and Polar Bears.

    Mol Phylogenet Evol 2000, 15:319-326. PubMed Abstract | Publisher Full Text OpenURL