Rise of oceanographic barriers in continuous populations of a cetacean: the genetic structure of harbour porpoises in Old World waters
- Michaël C Fontaine1, 2Email author,
- Stuart JE Baird2,
- Sylvain Piry2,
- Nicolas Ray3,
- Krystal A Tolley4, 5,
- Sarah Duke6,
- Alexei BirkunJr7,
- Marisa Ferreira8,
- Thierry Jauniaux9,
- Ángela Llavona10,
- Bayram Öztürk11,
- Ayaka A Öztürk11,
- Vincent Ridoux12,
- Emer Rogan13,
- Marina Sequeira14,
- Ursula Siebert15,
- Gísli A Vikingsson16,
- Jean-Marie Bouquegneau1 and
- Johan R Michaux2, 17
© Fontaine et al; licensee BioMed Central Ltd. 2007
Received: 29 December 2006
Accepted: 25 July 2007
Published: 25 July 2007
Understanding the role of seascape in shaping genetic and demographic population structure is highly challenging for marine pelagic species such as cetaceans for which there is generally little evidence of what could effectively restrict their dispersal. In the present work, we applied a combination of recent individual-based landscape genetic approaches to investigate the population genetic structure of a highly mobile extensive range cetacean, the harbour porpoise in the eastern North Atlantic, with regards to oceanographic characteristics that could constrain its dispersal.
Analyses of 10 microsatellite loci for 752 individuals revealed that most of the sampled range in the eastern North Atlantic behaves as a 'continuous' population that widely extends over thousands of kilometres with significant isolation by distance (IBD). However, strong barriers to gene flow were detected in the south-eastern part of the range. These barriers coincided with profound changes in environmental characteristics and isolated, on a relatively small scale, porpoises from Iberian waters and on a larger scale porpoises from the Black Sea.
The presence of these barriers to gene flow that coincide with profound changes in oceanographic features, together with the spatial variation in IBD strength, provide for the first time strong evidence that physical processes have a major impact on the demographic and genetic structure of a cetacean. This genetic pattern further suggests habitat-related fragmentation of the porpoise range that is likely to intensify with predicted surface ocean warming.
In the marine realm, pelagic species that have large geographic range and high dispersal capabilities represent a serious challenge to the idea of allopatric divergence (i.e., a large continuous population broken up into smaller units by extrinsic barriers) and to speciation processes in a seemingly continuous environment . The high mobility of these species and the dearth of barriers to gene flow in oceans might be expected to limit the division of species' ranges and, as a result, even distant regions might be connected genetically [1, 2]. Although examples of genetic homogeneity over large distances are common in marine systems, there are also many examples of surprising population structure in marine species with high dispersal potential [1, 3–7].
Cetaceans are good examples of this kind of species. Despite their broad range and their high dispersal capabilities, many cetaceans often show substantial genetic structure at regional or even fine scale, although the extent varies among species . It is generally argued that these patterns, not always correlated with geographic features, are related to a combination of complex behaviours, such as philopatry, specialisations for local resources, or social organisation into kinship groups [8, 9]. On the other hand, while the dispersal and segregation of populations of terrestrial mammals are frequently influenced by geographic features or climatic characteristics, few such obvious barriers are expected to restrict cetacean dispersal and gene flow in the world's oceans [10, 11]. Variation in oceanographic properties of the water column, such as depth, temperature, currents and winds, are known as important factors in the life of these animals, most obviously in conditioning the availability of their food (for example, see ), but their effect on cetacean dispersal and on population structure remains enigmatic.
Here, we provide strong and clear evidence that seascape imposes major constraints on the demographic and genetic structure of a cetacean, and thus on its dispersal. This finding is of general interest in the context of climate change and habitat fragmentation for marine species, as ecosystems in the eastern North Atlantic are shifting toward a warmer dynamic equilibrium with significant changes already detected in plankton and fish assemblages.
We applied two complementary Bayesian clustering algorithms, namely Structure v.2.1 [18, 19] and Geneland v.1.0.7 , to infer population structure (i.e., a number of clusters, K) and to assign individuals (probabilistically) to populations (or clusters) based on individual multilocus genotypes and, for the second algorithm, also on individual spatial origins. Both of these approaches assume that populations are panmictic units with distinct allele frequencies. To test whether individual dispersal is restricted in space, we analysed the pattern of isolation by distance (IBD) using the individual-based approach developed by Rousset . This involves regression of an index of genetic differentiation on marine geographic distance among pairs of individuals (see Methods). Finally, recent migration among populations (within the last few generations) was assessed using a Bayesian model implemented in BayesAss v.1.3 . This algorithm requires few assumptions for assigning individual genotypes to population of origin and, in particular, relaxes the key assumption of Hardy-Weinberg (HW) equilibrium within populations.
While Structure uses only the individual multilocus genotype data to infer the population structure, Geneland also exploits the spatial positions of the individual samples as a supplemental parameter in the analysis . Interesting features of the Geneland model that further distinguish it from that of the Structure model are its ability (a) to deal with an unknown number of populations simultaneously with other parameters, (b) to locate population boundaries across space, and (c) to account for uncertainty in the positioning of sampled individuals (see Methods and  for further details). This last feature is particularly useful in the present context as the locations of sampled harbour porpoises, composed of by-caught and stranded animals, might be poorly representative of the normal range of individuals.
Genetic diversity and differentiation among inferred populations
Genetic variation at the 10 microsatellite loci for populations inferred from the cluster analyses
H o /H e
H o /H e
H o /H e
Tests of departure from HW equilibrium (Table 1) show no significant deviation for porpoises from Iberia and the Black Sea, but a significant deficiency in heterozygosity at 9 of the 10 loci analysed in the NAt cluster. This slight heterozygote deficiency recorded at almost all loci in porpoises of northern Atlantic waters and the failure of the Structure model to assign these individuals in comparison to the Geneland model suggest that a subtle spatial structure (i.e., Wahlund effect) with a continuous gradation in allele frequencies across regions and/or isolation by distance could occur [23, 25].
Isolation by distance analyses
Isolation by distance within the North Atlantic cluster (NAt). Three different scales defined from the latitudinal subdivision of the global range in two (Medium scale) and in three (Small scale) parts were analysed. See Figure 1 for the delimitations of the NAt subdivisions.
Mean (max) marine distance (km)
Slope ± SE
Intercept ± SE
4πDσ 2 (1/slope)
0.0037 ± 0.0015
0.0011 ± 0.0148
0.0080 ± 0.0025
-0.0566 ± 0.0293
0.0028 ± 0.0018
0.0153 ± 0.0154
0.0100 ± 0.0031
-0.0867 ± 0.0347
0.0025 ± 0.0025
0.0244 ± 0.0292
0.0030 ± 0.0020
0.0167 ± 0.0174
We found a significant positive relationship between the index of genetic differentiation (a r ) and the marine geographic distance among porpoises in the NAt cluster at all scales considered (Table 2) except one: the region NAt-3B. This latter corresponds to the area where the sample size is the lowest (n = 141), where the sampling is the most spatially heterogeneous (Figure 1), and also where the marine distances among porpoises are the shortest (Table 2). Therefore, the absence of significant evidence in this region likely results from the low power of the analysis to detect IBD (see, for example, ).
Rousset [21, 26] demonstrated that the regression slope is proportional to 1/4πDσ 2, where D is the effective density of individuals and σ2 the second moment of axial dispersal distance, best described as the mean squared parent-offspring axial dispersal distance. σ2 can be understood as a measure of the speed at which two gene lineages issuing from an ancestor move away from each other, as it is the rate at which the mean squared axial distance between these two lineage increases per time unit . The comparison among subset areas at the medium and at small scale showed significant north-south variation in the parameters of the regression for the 10 microsatellite loci (Table 2). The slope (or 1/4πDσ 2) in the south part of the NAt cluster was significantly higher than that in northern parts at medium scale (Wilcoxon paired-sample test, 2A-2B: p = 0.037) and at small scale (Wilcoxon paired-sample test, 3A-3B: p = 0.005; 3A-3C: p = 0.046; 3B-3C: p = 0.399), suggesting that either density (D), dispersal (σ2), or both are reduced in the south part compared to the north.
Recent migration rates among populations
Mean ± SD (95% CI) posterior distributions for migration rates among harbour porpoise populations. Values along the diagonal (bold) are the proportion of individuals derived from the source population (or non-migrant) each generation.
Migration rate from
0.996 ± 0.004
0.002 ± 0.003
0.002 ± 0.003
(n = 78)
0.010 ± 0.011
0.978 ± 0.017
0.011 ± 0.011
(n = 31)
0.003 ± 0.003
0.031 ± 0.012
0.965 ± 0.013
(n = 303)
Virtually all porpoises from the Black Sea were identified as non-migrant (Table 3). Although this result is not surprising, as the Black Sea population is now geographically isolated from the Atlantic populations by the Mediterranean Sea, this result can be useful as reference to assess the status of the Iberian population. Almost all porpoises from Iberian waters were also identified as non-migrant (98% of the individuals and the 95% CI upper limit including 1), while the NAt cluster showed a slightly lower non-migrant proportion (96%; Table 3). The migration rates between Iberia and the NAt cluster were low (m ≤ 0.03) with the lower 95% CI bounds not different from 0, except in one case: the migration rate from Iberia to the NAt cluster appeared slightly higher than the reverse, but the large overlap of 95% CIs did not allow us to conclude there was asymmetry in migration rates.
The individual-based approaches we used here revealed that most of the harbour porpoise range in the central and eastern North Atlantic behaves as a 'continuous' population that widely extends over thousands of kilometres from the French coasts of the Bay of Biscay northwards to the arctic waters of Norway and Iceland, with significant isolation by distance. This striking result is concordant with the low but sometimes significant level of genetic differentiation previously reported at microsatellite loci between arbitrarily defined groups in the North Sea and adjacent waters [31, 32]. However, strong barriers to gene flow in the south-eastern North Atlantic range isolate, on a relatively small scale, porpoises from Iberian waters and on a larger scale porpoises from the Black Sea.
The total isolation of harbour porpoises from the Black Sea has long been suggested on the basis of the lack of field observation of porpoises in the Mediterranean Sea , of private mtDNA alleles reported in that population , and of morphological differences . Our results lend further support to this hypothesis. The pronounced genetic footprint of this isolation left at nuclear and mtDNA loci suggest this is an ancient isolation that might date back to the last Ice Age ( and Fontaine, unpublished results). The genetic differentiation detected at microsatellite loci between the Iberian porpoises and those further north was not apparent at the mtDNA control region previously analysed . The lack of mitochondrial lineage sorting and of private microsatellite alleles suggests that the differentiation we observed with microsatellite analyses is caused by a more recent isolation process than that of the Black Sea.
Still further north, depth increases towards Nordic Seas (Figure 1), but waters remain cold and highly productive . However, the bathymetric change does not seem to restrict gene flow in Nordic Seas, consistent with sightings of some porpoises reported far offshore in deep water . While this suggest there are few, if any, potential barriers to dispersal of porpoises from the northern Bay of Biscay up to Arctic waters, the heterozygosity deficit related to the detected IBD shows nevertheless that porpoises do not mate randomly over that extended area and that gene flow is spatially restricted. We observed a north-south variation in the IBD pattern with higher IBD slope at the southern end of this range compared to northern parts (Table 2). One could argue that this north-south variation in IBD pattern might reflect drift disequilibrium  in northern areas associated with the postglacial porpoise recolonisation of Nordic waters in contrast to the southern habitats, which likely remained more stable in time. However, simulation-based sensitivity analysis of current Dσ 2 estimation to demographic instability in time and space conducted by Leblois et al  showed that spatial expansion with constant density does not significantly affect present-time Dσ 2 estimation, especially when the spatial expansion occurred 20 or more generations ago, as it is the case for postglacial recolonisation. Consequently, the higher IBD slope detected in the southern area (NAt-3A) compared to that in waters further north (NAt-3B and NAt-3C) most likely represents a lower current-time Dσ 2. Although we cannot exclude variation in σ2, a lower porpoise density in southern waters is supported by field estimates based on aerial and ship surveys conducted in the North Sea and adjacent waters [45, 46]. These variations in density (and maybe in dispersal patterns) likely reflect variation in habitat. The southern part of the 'continuous' population (i.e., the northern part of the Bay of Biscay, the English Channel and the southernmost part of the North Sea) borders the barrier to gene flow detected in the southern Bay of Biscay and should thus display sub-optimal conditions for porpoises while the middle (i.e., the central and northern North Sea) and northern areas (i.e., the Nordic Seas) would be more optimal for a cold water species such as the harbour porpoise.
The Mediterranean Sea displays similar characteristics to those encountered in the southern Bay of Biscay but at a larger scale. The Mediterranean is composed mostly of deep basins and narrow continental shelves with warm oligotrophic surface waters prevailing most of the year . These characteristics are likely quite unfavourable for cold water species and might explain why the harbour porpoise is absent from this area. The oceanographic conditions in the Black Sea are, by contrast, more suitable for harbour porpoises with low salinity, colder and more nutrient rich surface waters than in the Mediterranean Sea . There are however reports of porpoise strandings in the northern Aegean Sea . This can be understood with regard to oceanographic features in that area. The subdivision of the Aegean into two basins has long been recognised. The northern basin is under the influence of cold, low salinity waters that pour out of the Black Sea. This water is entrained into a cyclonic circulation affecting the northern and western parts of the Aegean, causing an ecological isolation of the northern basin from the southern basin . In the southern basin the continental shelf is very limited and the waters become quickly characteristic of Mediterranean waters , unfavourable for harbour porpoises.
To summarise, surface water temperature and primary production seem to be the factors that best characterise the nature of barriers to gene flow encountered across the harbour porpoise range, their population structure, and their geographic distribution. It is worth noting however that in oceanography, these two parameters are often linked . Indeed, the sea surface temperature acts as a useful proxy for other physical processes, such as vertical stratification and nutrient contents, regulating the size structure, taxonomic composition, and abundance of the phytoplankton community, and thus the food availability for top predators [52, 53]. These results reinforce previous ecological studies on harbour porpoises that reported significant relationships between abundances and movements with sea surface temperature and food availability [54, 55]. Although bathymetry can be important in harbour porpoise ecology [56, 57], we showed that this factor alone seems not to restrict gene flow in northern waters of the sampling range.
While the proximal causes of porpoise dependence on these habitat characteristics are beyond the scope of this paper, the ultimate underlying mechanism is likely related to the high energetic constraints this small cetacean has to face in order to survive. As one of the smallest endothermic marine predators, and furthermore with limited energy storage capacity, it is currently assumed that harbour porpoises must feed frequently without prolonged periods of fasting [16, 58]. Their distribution, their movements, and in sum their overall biology should therefore be closely related to those of their prey and thus to nutrient rich waters.
In the marine realm, community structure is shaped heavily by physical processes (see, for example, [47, 59]). In this study we provide for the first time strong evidence that physical processes determining food availability have major impacts on the demographic and genetic structure of a cetacean. The small body size of harbour porpoises undoubtedly has profound consequences at all levels of their biology and makes this species particularly sensitive to habitat variation. We can however reasonably expect that this will be also applicable to other cetaceans of similar body size, habitat and thermoregulation constraints. However, these constraints could be reduced for larger cetaceans, leading to more complex patterns of population structure not necessarily correlated to seascape features (see, for example, ). Ecosystems in the eastern North Atlantic are shifting toward a warmer dynamic equilibrium with significant changes detected in plankton and fish assemblages [51, 60–62], but the consequences for marine mammals remain to date unclear . Although further analyses would be require to address the demographic trends of these populations, the genetic pattern highlighted here (i.e., the ancient isolation of harbour porpoises in the Black Sea), the more recent isolation of those in Iberian waters, and the higher IBD in the southern end of the northern Atlantic continuum, suggests that habitat-related fragmentation of harbour porpoise range is under way and that it is likely to continue with the predicted changes in climate.
Tissue samples were taken from by-caught and stranded harbour porpoises. A total of 752 animals distributed along the eastern North Atlantic range of the harbour porpoise and in the Black Sea were analysed (Figure 1). Out of these, 515 samples were analysed in this study and 237 samples from Iceland and Ireland were analysed by Duke .
Most of the individuals were geo-referenced using GPS coordinates recorded at the time or deduced from the reported location where the animal was found. These coordinates are naturally rough approximations to the normal locations of animals, especially for stranded animals, but this error can be considered negligible at the scale of the study range. This source of error can also be taken into account in some of the spatial analyses (see below).
DNA extraction and microsatellite analysis
Total genomic DNA was extracted from tissues using the DNeasy™ Tissue Kit (Qiagen) following the manufacturer's recommendations. Samples were genotyped at 10 microsatellite loci using the multiplex sets defined in . Polymerase chain reaction conditions were as reported in . Amplified DNA was analysed for length variations on an automated 96 capillary MegaBace-1000 DNA Analyser (Amersham Biosciences) using Genetic Profiler v.1.5 (Amersham Biosciences).
Data on habitat characteristics across the study range with respect to salinity and sea surface temperature were taken from the National Oceanographic Data Centre (NODC) . Bathymetric data were extracted from the ETOPO2 dataset available on the US National Geophysical Data Centre (NGDC)  and the data on chlorophyll concentration were taken from the NASA Sea-viewing Wide Field-of-view Sensor database (SeaWIFS) .
We applied two Bayesian model-based clustering algorithms to infer population structure and to assign individuals (probabilistically) to clusters without a priori knowledge of population units and limits.
The first approach, implemented in Structure v.2.1, uses individual multilocus genotype data to cluster individuals into K groups while minimising Hardy-Weinberg disequilibrium and gametic phase disequilibrium between loci within groups [18, 19]. The estimation procedure consists of running trial values of the number of populations K and then comparing the estimated log probability of data under each K, Ln [Pr(X|K)]. We conducted a series of independent runs with different proposals for K, testing all values from 1 to 10. Each runs used 106 iterations after a burn-in of length 4 × 104, testing different models: (a) with or without admixture, and (b) correlated or uncorrelated allele frequencies. To check for convergence of the Markov chain Monte Carlo (MCMC), we performed 10 replicates for each value of K and then checked the consistency of results. The estimated number of clusters (K) was taken to be the value of K with the highest Pr(X|K) .
The second algorithm, implemented in Geneland v.1.0.7, differs from that of Pritchard et al  mainly by taking into account explicitly the spatial dependence of individuals expected for species whose range is much larger than the average intergeneration movement of individuals. This model aims at inferring and locating genetic discontinuities between populations in space from individual geo-referenced multilocus genotypes, while taking into account uncertainty in the location of sampled individuals [20, 69]. All the parameters (including K) are co-estimated simultaneously by the MCMC algorithm. However, for technical reasons discussed in , it is better to proceed in two steps: a first run to infer K, and a second run with K fixed at the modal value to estimate the other parameters (mainly the assignment of individuals to the inferred populations). The first step was replicated 10 times to check for convergence, allowing K to vary from 1 to 10 clusters and using the following run parameters: 106 MCMC iterations, maximum rate of Poisson process fixed at 700, maximum number of nuclei in the Poisson-Voronoi tessellation fixed to 500, and an uncertainty associated with the spatial coordinates of 50 km. We used the Dirichlet model of allelic frequencies as it has been demonstrated to perform better than the alternative model . We inferred the number of clusters (K) from the modal value of K for these 10 runs, and then ran the MCMC again 100 times with K fixed for this value, 5 × 105 MCMC iterations, and the other parameters unchanged. We calculated the mean logarithm of posterior probability of the data (PPD) for each of the 100 runs and selected the 10 with the highest PPD. These 10 runs were then post-processed (with a burn-in of 5 × 104 iterations) in order to obtain posterior probabilities of population membership for each individual and each pixel of the spatial domain (174 pixels along the X axis and 143 along the Y axis corresponding to a pixel size of 25 km side). We finally checked visually for the consistency of results across these 10 runs.
Descriptive statistics among clusters
The allelic richness, corrected for difference in sample size, the observed (H o ) and expected (H e ) heterozygosity (or genetic diversity), and F IS values were calculated within each cluster using Fstat v.2.9.3 . To test whether genetic diversity was significantly different between clusters, we applied a Wilcoxon paired-sample test  on the 10 single locus values of the statistics of interest.
Level of genetic differentiation at microsatellite loci among clusters was estimated as F ST after Weir and Cockerham  using Fstat v.2.9.3 . The 95% confidence interval was calculated using 15000 bootstrap resamplings . We conducted exact tests to assess deviations from Hardy-Weinberg equilibrium and test for population differentiation using Genepop v.3.4 .
Isolation by distance analysis
In continuous populations, an isolation by distance pattern occurs when genetic differentiation among individuals increases with their geographic distance . Here we consider the statistic a r , a multilocus estimator of an F ST/(1-F ST) analogue between pairs of individuals . When a continuous population is represented by a two dimensional lattice (i.e., fixed individual positions and no spatial density heterogeneity), a r is approximately linearly related to the logarithm of the geographic distance between individuals (r), a r ≈ (Ln(r)/4πDσ 2) + C, where D is the effective density of individuals, σ2 is the second moment of the dispersal distance distribution, and C is the value of the linear approximation at r = 1 length unit. Values of a r were regressed against the log of the marine geographical distance (see below) between paired individuals, as described in Rousset et al . Significance of the regression slope was tested by 105 random permutations of individual locations (similar to a Mantel test) using the computer program SPAGeDi v.1.2 . Assuming low mutation rate, the inverse of the regression slope provides an estimate of the product 4πDσ 2 [21, 26]. To test whether the regression slopes significantly differed between the different parts of a same scale, we used a Wilcoxon paired-sample test  applied on the 10 single locus values of the regression slope.
In the marine realm, the Euclidean distance between individuals might not be representative of the effective geographic distance separating them. Therefore, we computed an effective marine geographic distance between individuals using the least-cost path (LCP) algorithm implemented in the Pathmatrix extension  of the geographical information system software ArcView v.3.X (ESRI, Redlands, CA, USA). This algorithm computes a deterministic LCP between a source point and a target point by using a friction (or resistance) layer. The friction layer is a raster map where each cell (landscape unit) expresses the relative difficulty (or cost) of moving through that cell. A LCP minimises the sum of costs of all cells along the path (for detailed description and discussion of the algorithm, see ). In the present study, a uniform cost was attributed to all sea cells, while land cells harboured an "infinite" cost. This allowed us to compute effective distances avoiding landmasses. The sea/land map was obtained by rasterising (at 2 km resolution) a polygon version  of the GSHHS shoreline dataset v.1.3 . The computations were performed using a gnomonic projection around the centroid of the sampled localities, which minimises the map deformation in planar distances induced by the curvature of the earth (Baird personal communication). Finally, the length of pairwise LCP (in meters) was introduced as the geographic distance matrix separating pairs of individuals in the regression analyses described above.
Migration rates among clusters
Evidence of recent migration events across clusters was assessed using the Bayesian multilocus genotyping procedure implemented with MCMC methods in BayesAss v.1.3 . This approach does not require populations to be in either migration-drift or Hardy-Weinberg equilibrium. To examine the strength of the information in the porpoise microsatellite data set, 95% confidence intervals were determined for migration rates and compared to a scenario where all proposed changes throughout the Markov chain are accepted (thereby simulating the situation where any information that could exist in the data is insufficient to affect the posterior distribution of migration rates, as suggested by the authors). The MCMC was run for a total of 3 × 106 iterations, with the first 106 discarded as a burn-in to allow the chain to reach stationarity. Samples were collected every 2000 iterations to infer posterior probability distributions of parameters of interest.
We are grateful to all the fishermen, stranding networks, and volunteers that contributed in the collection of samples used in the present study. Specifically, we thank D. Bloch (Museum of Natural History, Faroe Islands), MJ Addink and C Smeenk (National Museum of Natural History, Leiden, the Netherlands), N Øien (Institute of Marine Research, Bergen, Norway), and W Dabin (Centre de Recherche sur les Mammifères Marins, La Rochelle, France). We thank also M Galan and A Loiseau for their help in laboratory work, and S Gobert and K Das for administrative support. R Streiff, A Estoup and R Leblois provided helpful assistance and critical comments on the data analyses, the interpretation of results, and on the manuscript. We also thank P Beaubrun, F Bonhomme, PJ Palsbøll, and C Moritz for their critical comments on the manuscript.
This work was funded by the Belgian Office for Scientific, Technical and Cultural Affairs (Contract EV/12/46A). MCF and JRM were supported by research fellowships from the Belgian National Fund for Scientific Research (F.R.S.FNRS, mandate 'Aspirant' and 'Chercheur Qualifié'), NR by a Swiss National Science Foundation grant (No 3100A0–112072), and KAT by the Norwegian Research Council. This is a MARE publication, no. 113.
- Palumbi SR: Genetic divergence, reproductive isolation and marine speciation. Annu Rev Ecol Syst. 1994, 25: 547-572. 10.1146/annurev.es.25.110194.002555.View ArticleGoogle Scholar
- Thorrold SR: Ocean ecology: don't fence me in. Curr Biol. 2006, 16: R638-R640. 10.1016/j.cub.2006.07.034.View ArticlePubMedGoogle Scholar
- Barber PH, Palumbi SR, Erdman MV, Moosa MK: A marine Wallace's line?. Nature. 2000, 406: 692-693. 10.1038/35021135.View ArticlePubMedGoogle Scholar
- Taylor MS, Hellberg ME: Genetic evidence for local retention of pelagic larvae in Caribbean reef fish. Science. 2003, 299: 107-109. 10.1126/science.1079365.View ArticlePubMedGoogle Scholar
- Bekkevold D, Andre C, Dahlgren TG, Clausen LA, Torstensen E, Mosegaard H, Carvalho GR, Christensen TB, Norlinder E, Ruzzante DE: Environmental correlates of population differentiation in Atlantic herring. Evolution. 2005, 59: 2656-2668.View ArticlePubMedGoogle Scholar
- Jorgensen HBH, Hansen MM, Bekkevold D, Ruzzante DE, Loeschcke V: Marine landscapes and population genetic structure of herring (Clupea harengus L.) in the Baltic Sea. Mol Ecol. 2005, 14: 3219-3234. 10.1111/j.1365-294X.2005.02658.x.View ArticlePubMedGoogle Scholar
- Kenchington EL, Patwary MU, Zouros E, Bird CJ: Genetic differentiation in relation to marine landscape in a broadcast-spawning bivalve mollusc (Placopecten magellanicus). Mol Ecol. 2006, 15: 1781-1796. 10.1111/j.1365-294X.2006.02915.x.View ArticlePubMedGoogle Scholar
- Hoelzel AR, Goldsworthy SD, Fleischer RC: Population genetics. Marine Mammal Biology: An Evolutionary Approach. Edited by: Hoelzel AR. 2002, Oxford: Blackwell Science, 325-352.Google Scholar
- Hoelzel AR: Genetic structure of cetacean populations in sympatry, parapatry, and mixed assemblages: implication for conservation policy. J Hered. 1998, 89: 451-458. 10.1093/jhered/89.5.451.View ArticleGoogle Scholar
- Palsbøll PJ, Clapham PJ, Mattila DK, Larsen F, Sears R, Siegismund HR, Sigurùnsson J, Vasquez O, Arctander P: Distribution of mtDNA haplotypes in North Atlantic humpback whales: the influence of behaviour on population structure. Mar Ecol Prog Ser. 1995, 116: 1-10.View ArticleGoogle Scholar
- Natoli A, Birkun A, Aguilar A, Lopez A, Hoelzel AR: Habitat structure and dispersal of male and female bottlenose dolphins (Tursiops truncatus). Proc R Soc B. 2005, 272: 1217-1226. 10.1098/rspb.2005.3076.PubMed CentralView ArticlePubMedGoogle Scholar
- Berta A, Sumich JL: Marine Mammals: Evolutionary Biology. 1999, London: Academic PressGoogle Scholar
- Read AJ: Porpoises, overview. Encyclopedia of Marine Mammals. Edited by: Perrin WF, Würsig B, Thewissen JGM. 2002, London: Academic Press, 982-985.Google Scholar
- Koopman HN, Pabst DA, McLellan WA, Dillaman RM, Read AJ: Changes in blubber distribution and morphology associated with starvation in harbour porpoise (Phocoena phocoena): evidence for regional variation in blubber structure and function. Physiol Biochem Zool. 2002, 75: 498-512. 10.1086/342799.View ArticlePubMedGoogle Scholar
- Gaskin DE: The harbour porpoise Phocoena phocoena (L.): regional populations, status, and information on direct and indirect catches. Rep Int Whal Commn. 1984, 34: 569-584.Google Scholar
- Read AJ: Harbour porpoise (Phocoena phocoena). Handbook of Marine Mammals. Edited by: Ridgway SH, Harrison R. 1999, London: Academic Press, 323-350.Google Scholar
- Frantzis A, Gordon J, Hassidis G, Komenou A: The enigma of harbour porpoise presence in the Mediterranean Sea. Mar Mamm Sci. 2001, 17: 937-943. 10.1111/j.1748-7692.2001.tb01307.x.View ArticleGoogle Scholar
- Pritchard JK, Stephens M, Donnelly P: Inference of population structure using multilocus genotype data. Genetics. 2000, 155: 945-959.PubMed CentralPubMedGoogle Scholar
- Falush D, Stephens M, Pritchard JK: Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies. Genetics. 2003, 164: 1567-1587.PubMed CentralPubMedGoogle Scholar
- Guillot G, Estoup A, Mortier F, Cosson J-F: A spatial model for landscape genetics. Genetics. 2005, 170: 1261-1280. 10.1534/genetics.104.033803.PubMed CentralView ArticlePubMedGoogle Scholar
- Rousset F: Genetic differentiation between individuals. J Evol Biol. 2000, 13: 58-62. 10.1046/j.1420-9101.2000.00137.x.View ArticleGoogle Scholar
- Wilson GA, Rannala B: Bayesian inference of recent migration rates using multilocus genotypes. Genetics. 2003, 163: 1177-1191.PubMed CentralPubMedGoogle Scholar
- Rosenberg NA, Pritchard JK, Weber JL, Cann HM, Kidd KK, Zhivotovsky LA, Feldman MW: Genetic structure of human populations. Science. 2002, 298: 2381-2385. 10.1126/science.1078311.View ArticlePubMedGoogle Scholar
- Weir BS, Cockerham CC: Estimating F-statistics for the analysis of population structure. Evolution. 1984, 38: 1358-1370. 10.2307/2408641.View ArticleGoogle Scholar
- Rosenberg NA, Mahajan S, Ramachandran S, Zhao C, Pritchard JK, Feldman MW: Clines, clusters, and the effect of study design on the inference of human population structure. PLoS Genet. 2005, 1: e70-10.1371/journal.pgen.0010070.PubMed CentralView ArticlePubMedGoogle Scholar
- Rousset F: Genetic differentiation and estimation of gene flow from F-statistics under isolation by distance. Genetics. 1997, 145: 1219-1228.PubMed CentralPubMedGoogle Scholar
- Slatkin M: Isolation by distance in equilibrium and non-equilibrium populations. Evolution. 1993, 47: 264-279. 10.2307/2410134.View ArticleGoogle Scholar
- Leblois R, Estoup A, Rousset F: Influence of mutational and sampling factors on the estimation of demographic parameters in a "Continuous" population under isolation by distance. Mol Biol Evol. 2003, 20: 491-502. 10.1093/molbev/msg034.View ArticlePubMedGoogle Scholar
- Palumbi SR: Marine reserves and ocean neighborhoods: the spatial scale of marine populations and their management. Annu Rev Environ Resourc. 2004, 29: 31-68. 10.1146/annurev.energy.29.062403.102254.View ArticleGoogle Scholar
- Rousset F: Genetic Structure and Selection in Subdivided Populations (Monographs in population biology edn). 2004, Princeton: Princeton University PressGoogle Scholar
- Andersen LW, Ruzzante DE, Walton M, Berggren P, Bjørge A, Lockyer C: Conservation genetics of harbour porpoises, Phocoena phocoena, in eastern and central North Atlantic. Conserv Genet. 2001, 2: 309-324. 10.1023/A:1012534212853.View ArticleGoogle Scholar
- Andersen LW: Harbour porpoises (Phocoena phocoena) in the North Atlantic: distribution and genetic population structure. Harbour Porpoises in the North Atlantic. Edited by: Tromsø: NAMMCO Scientific Publications. 2003, Haug T, Desportes G, Vikingsson GA, Witting L, 11-30.Google Scholar
- Rosel P, Dizon AE, Haygood MG: Variability of the mitochondrial control region in populations of the harbour porpoise, Phocoena phocoena, on interoceanic and regional scales. Can J Fish Aquat Sci. 1995, 52: 1210-1219.View ArticleGoogle Scholar
- Gol'din PE: Growth and body size of the harbour porpoise, Phocoena phocoena (Cetacea, Phocoenidae), in the Sea of Azov and the Black Sea. Vestnik Zoologii. 2004, 38: 59-73.Google Scholar
- Tolley KA, Rosel PE: Population structure and historical demography of eastern North Atlantic harbour porpoises inferred through mtDNA sequences. Mar Ecol Prog Ser. 2006, 297-308. 10.3354/meps327297.Google Scholar
- Koutsikopoulos C, Beillois P, Leroy C, Taillefer F: Temporal trends and spatial structure of the sea surface temperature in the Bay of Biscay. Oceanol Acta. 1998, 21: 335-344. 10.1016/S0399-1784(98)80020-0.View ArticleGoogle Scholar
- Koutsikopoulos C, Le Cann B: Physical processes and hydrological structures related to the Bay of Biscay anchovy. Sci Mar. 1996, 60: 9-19.Google Scholar
- Fiùza AFG: Upwelling patterns off Portugal. Coastal Upwelling: Its Sediment Records (part A). Edited by: Suess E, Thiede J. 1983, New York: Plenum Press, 85-98.View ArticleGoogle Scholar
- Tenore KR, Alonso-Noval M, Alvarez-Ossorio M, Atkinson LP, Cabanas JM, Cal RM, Campos HJ, Castillejo F, Chesney EJ, Gonzalez N, et al: Fisheries and oceanography off Galicia, NW Spain: mesoscale spatial and temporal changes in physical processes and resultant patterns of biological productivity. J Geophys Res. 1995, 100: 10943-10966. 10.1029/95JC00529.View ArticleGoogle Scholar
- Southward AJ, Hawkins SJ, Burrows MT: Seventy years' observations of changes in distribution and abundance of zooplankton and intertidal organisms in the western English Channel in relation to rising sea temperature. J Therm Biol. 1995, 20: 127-155. 10.1016/0306-4565(94)00043-I.View ArticleGoogle Scholar
- OSPAR Commission: Quality Status Report 2000: Region I – Arctic Water. 2000, London: OSPAR CommissionGoogle Scholar
- Donovan GP, Bjørge A: Harbour porpoises in the North Atlantic: edited extract from the report of the IWC scientific committee, Dublin 1995. Report of the International Whaling Commission. 1995, 3-25. Special Issue 16Google Scholar
- Hutchison DW, Templeton AR: Correlation of pairwise genetics and geographic distance measures: inferring the relative influences of gene flow and drift on the distribution of genetic variability. Evolution. 1999, 53: 1898-1914. 10.2307/2640449.View ArticleGoogle Scholar
- Leblois R, Rousset F, Estoup A: Influence of spatial and temporal heterogeneities on the estimation of demographic parameters in a continuous population using individual microsatellite data. Genetics. 2004, 166: 1081-1092. 10.1534/genetics.166.2.1081.PubMed CentralView ArticlePubMedGoogle Scholar
- Hammond PS, Berggren P, Benke H, Borchers DL, Collet A, Heide-Jorgensen MP, Heimlich S, Hiby AR, Leopold MF, Øien N: Abundance of harbour porpoises and other cetaceans in the North Sea and adjacent waters. J Applied Ecol. 2002, 39: 361-376. 10.1046/j.1365-2664.2002.00713.x.View ArticleGoogle Scholar
- SCANS-II : Quarterly newsletter for the small cetaceans in the European Atlantic and North Sea project. Special Issue 9: Survey of the SCANS-II project. [http://biology.st-andrews.ac.uk/scans2/documents/issue9_Dec06.pdf]
- Longhurst AR: Ecological Geography of the Sea. 1998, London: Academic PressGoogle Scholar
- Özsoy E, Ünlüata Ü: Oceanography of the Black Sea: a review of some recent results. Earth Sci Rev. 1997, 42: 231-272. 10.1016/S0012-8252(97)81859-4.View ArticleGoogle Scholar
- Theocharis A, Georgopoulos D, Lascaratos A, Nittis K: Water masses circulation in the central region of the Eastern Mediterranean: Eastern Ionian, South Aegean, and Northwest Levantin 1986–1987. Deep-Sea Res II. 1993, 40: 1121-1142. 10.1016/0967-0645(93)90064-T.View ArticleGoogle Scholar
- Moraitou-Apostolopoulou M: The zooplankton communities of the Eastern Mediterranean (Levantine basin, Aegean Sea); influence of man-made factors. Mediterranean Marine Ecosystems. Edited by: Moraitou-Apostolopoulou M, Kiortsis V. 1985, New York: Plenum Press, 303-331.View ArticleGoogle Scholar
- Richardson AJ, Schoeman DS: Climate impact on plankton ecosystems in the Northeast Atlantic. Science. 2004, 305: 1609-1612. 10.1126/science.1100958.View ArticlePubMedGoogle Scholar
- Sathyendranath S, Cota G, Stuart V, Maass H, Platt T: Remote sensing of phytoplankton pigments: a comparison of empirical and theoretical approaches. Int J Remote Sens. 2001, 22: 249-273. 10.1080/014311601449925.View ArticleGoogle Scholar
- Behrenfeld MJ, O'Malley RT, Siegel DA, McClain CR, Sarmiento JL, Feldman GC, Millingan AJ, Falkowski PG, Letelier RM, Boss ES: Climate-driven trends in contemporary ocean productivity. Nature. 2006, 444: 752-755. 10.1038/nature05317.View ArticlePubMedGoogle Scholar
- Forney KA: Trends in harbour porpoise abundance off central California, 1986–95: evidence for interannual changes in distribution?. J Cetacean Res Manage. 1999, 1: 73-80.Google Scholar
- Johnston DW, Westgate AJ, Read AJ: Effects of fine-scale oceanographic features on the distribution and movements of harbour porpoises Phocoena phocoena in Bay of Fundy. Mar Ecol Prog Ser. 2005, 295: 279-293. 10.3354/meps295279.View ArticleGoogle Scholar
- Carretta JV, Taylor BL, Chivers S: Abundance and depth distribution of harbor porpoise (Phocoena phocoena) in northern California determined from a 1995 ship survey. Fish Bull. 2001, 99: 29-39.Google Scholar
- Fontaine MC, Tolley KA, Siebert U, Gobert S, Lepoint G, Bouquegneau JM, Das K: Long-term feeding ecology and habitat use in harbour porpoises Phocoena phocoena from Scandinavian waters inferred from trace elements and stable isotopes. BMC Ecol. 2007, 7: 1-10.1186/1472-6785-7-1.PubMed CentralView ArticlePubMedGoogle Scholar
- MacLeod CD, Santos MB, Reid RJ, Scott BE, Pierce GJ: Linking sandeel consumption and the likelihood of starvation in harbour porpoises in the Scottish North Sea: could climate change mean more starving porpoises?. Biol Lett. 2007, 3: 185-188. 10.1098/rsbl.2006.0588.PubMed CentralView ArticlePubMedGoogle Scholar
- Li WKW: Macroecological patterns of phytoplankton in the northwestern North Atlantic. Nature. 2002, 419: 154-157. 10.1038/nature00994.View ArticlePubMedGoogle Scholar
- Beaugrand G, Reid PC, Ibanez F, Lindley JA, Edwards M: Reorganization of North Atlantic marine copepod biodiversity and climate. Science. 2002, 296: 1692-1694. 10.1126/science.1071329.View ArticlePubMedGoogle Scholar
- Perry AL, Low PL, Ellis JR, Reynolds JD: Climate change and distribution shifts in marine fishes. Science. 2005, 308: 1912-1915. 10.1126/science.1111322.View ArticlePubMedGoogle Scholar
- Attrill MJ, Power M: Climatic influence on marine fish assemblage. Nature. 2002, 417: 275-278. 10.1038/417275a.View ArticlePubMedGoogle Scholar
- Learmonth JA, MacLeod CD, Santos MB, Pierce GJ, Crick HQP, Robinson RA: Potential effects of climate change on marine mammals. Oceanogr Mar Biol. 2006, 44: 431-464.Google Scholar
- Duke S: The Population and the Social Structure of Harbour Porpoise (Phocoena phocoena) from Around the Coasts of Iceland and Ireland. 2003, Dublin: PhD thesis, Department of Zoology, University College DublinGoogle Scholar
- Fontaine MC, Galan M, Bouquegneau JM, Michaux JR: Efficiency of fluorescent multiplex polymerase chain reactions (PCRs) for rapid genotyping of harbour porpoises Phocoena phocoena with 11 microsatellite loci. Aquat Mamm. 2006, 32: 301-304. 10.1578/AM.32.3.2006.301.View ArticleGoogle Scholar
- National Oceanographic Data Centre. [http://www.nodc.noaa.gov/]
- National Geophysical Data Centre. [http://www.ngdc.noaa.gov/products/ngdc_products.html]
- NASASea-viewing Wide Field-of-view Sensor database. [http://oceancolor.gsfc.nasa.gov/]
- Guillot G, Mortier F, Estoup A: GENELAND: a computer package for landscape genetics. Mol Ecol Notes. 2005, 5: 712-715. 10.1111/j.1471-8286.2005.01031.x.View ArticleGoogle Scholar
- Goudet J: FSTAT, a program to estimate and test gene diversities and fixation indices (version 2.9.3). [http://www2.unil.ch/popgen/softwares/fstat.htm]
- Zar JH: Biostatistical analysis. 1999, Upper Saddle River, NJ: Prentice-Hall, 4Google Scholar
- Raymond M, Rousset F: GENEPOP (Version 1.2): population genetics software for exact tests and eucumenism. J Hered. 1995, 86: 248-249.Google Scholar
- Wright S: Isolation by distance. Genetics. 1943, 28: 114-138.PubMed CentralPubMedGoogle Scholar
- Hardy OJ, Vekemans X: SPAGeDi: a versatile computer program to analyse spatial genetic structure at the individual of population levels. Mol Ecol Notes. 2002, 2: 618-620. 10.1046/j.1471-8286.2002.00305.x.View ArticleGoogle Scholar
- Ray N: PATHMATRIX: a geographical information system tool to compute effective distances among samples. Mol Ecol Notes. 2005, 5: 177-180. 10.1111/j.1471-8286.2004.00843.x.View ArticleGoogle Scholar
- Adriaensen F, Chardon JP, De Blust G, Swinnen E, Villalba S, Gulinck H, Matthysen E: The application of 'least-cost' modelling as a functional landscape model. Land Urb Plan. 2003, 64: 233-257. 10.1016/S0169-2046(02)00242-6.View ArticleGoogle Scholar
- Nelson A, Turner A: Processing global self-consistent hierarchical high resolution shoreline data version 1.2 Into ESRI ArcGIS vector and raster data. 2004, Leeds: CCG, School of Geography Working Paper, University of LeedsGoogle Scholar
- Wessel P, Smith WHF: A global self-consistent, hierarchical, high-resolution shoreline database. J Geophys Res. 1996, 101: 8741-8743. 10.1029/96JB00104.View ArticleGoogle Scholar
- Rosenberg NA: DISTRUCT: a program for graphical display for population structure. Mol Ecol Notes. 2004, 4: 137-138. 10.1046/j.1471-8286.2003.00566.x.View ArticleGoogle Scholar
- Feldman GC, McClain CR: Seasonal, monthly, and weekly climatologies. Ocean Color Web. Edited by: Kuring N, Bailey SW, Thomas D, Franz BF, Meister G, Werdell PJ, Eplee RE, MacDonald M, Rubens M. 2006, Greenbelt, MD: NASA Goddard Space Flight Center, [http://oceancolor.gsfc.nasa.gov/cgi/climatologies.pl]Google Scholar
This article is published under license to 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 cited.