- Research article
- Open Access
Does genetic structure reflect differences in non-breeding movements? A case study in small, highly mobile seabirds
© The Author(s). 2017
- Received: 3 April 2017
- Accepted: 27 June 2017
- Published: 5 July 2017
In seabirds, the extent of population genetic and phylogeographic structure varies extensively among species. Genetic structure is lacking in some species, but present in others despite the absence of obvious physical barriers (landmarks), suggesting that other mechanisms restrict gene flow. It has been proposed that the extent of genetic structure in seabirds is best explained by relative overlap in non-breeding distributions of birds from different populations. We used results from the analysis of microsatellite DNA variation and geolocation (tracking) data to test this hypothesis. We studied three small (130–200 g), very abundant, zooplanktivorous petrels (Procellariiformes, Aves), each sampled at two breeding populations that were widely separated (Atlantic and Indian Ocean sectors of the Southern Ocean) but differed in the degree of overlap in non-breeding distributions; the wintering areas of the two Antarctic prion (Pachyptila desolata) populations are separated by over 5000 km, whereas those of the blue petrels (Halobaena caerulea) and thin-billed prions (P. belcheri) show considerable overlap. Therefore, we expected the breeding populations of blue petrels and thin-billed prions to show high connectivity despite their geographical distance, and those of Antarctic prions to be genetically differentiated.
Microsatellite (at 18 loci) and cytochrome b sequence data suggested a lack of genetic structure in all three species. We thus found no relationship between genetic and spatial structure (relative overlap in non-breeding distributions) in these pelagic seabirds.
In line with other Southern Ocean taxa, geographic distance did not lead to genetic differences between widely spaced populations of Southern Ocean petrel species.
- Falkland / Malvinas Islands
- Genetic structure
- Kerguelen Islands
- Non-breeding distribution
- South Georgia
- Spatial distribution
The Southern Ocean (water masses south of the Subtropical Front) is characterized by an environmental gradient from subantarctic to polar waters, and the fauna reflect adaptations to the thermal and latitudinal structure in this region (e.g. [1, 2]). Seabirds in the Southern Ocean are faced with a limited availability of habitable islands, resulting in disjunctive breeding distributions, with populations often located at island groups that are several thousand kilometres apart . Because geographical barriers are absent across the Southern Ocean, physical isolation by geographic landmarks may not restrict gene flow. However, distance as well as intrinsic barriers, such as site fidelity, can reduce dispersal and lead to genetic structuring of seabird populations [4, 5]. Geographical distance can result in local differentiation (isolation by distance), the extent of which increases as connectivity decreases [4, 6]. However, the movements of pelagic seabirds, which are extensive during the non-breeding season, can counteract the geographical isolation of their breeding colonies. Hence, integrative studies of both population genetic structure and at-sea distributions are required to understand evolutionary processes in these taxa. An overlap of nonbreeding areas may facilitate contact and pair-bonds between individuals of different populations, even when their breeding areas are far apart. A direct study of site fidelity and movements between colonies in these species, e.g. using mark-recapture techniques, is difficult due to the remoteness and size of the breeding colonies, which can be as large as two million pairs (e.g. Thin-billed prions Pachyptila belcheri on New Island, Falkland Islands). Thus, genetic methods are required.
Several studies that have analysed patterns of genetic diversity in seabird species with circumpolar distributions in the Southern Ocean, have found little genetic and phylogeographic structure despite the large distances between colonies (e.g. Adélie penguins Pygoscelis adeliae, ; grey-headed albatrosses Thalassarche chrysostoma, ; wandering albatrosses Diomedea exulans, ). In contrast, the black-browed albatross T. melanophris displays genetic structure despite the absence of physical barriers, potentially because of its stronger association with continental shelf habitats; this conclusion was based on the comparison with the closely-related grey-headed albatross, which forages mainly at frontal systems which are circumpolar . Similarly, genetic structure has been found in other species associated with foraging on land or over continental shelves (giant petrels ; skuas ). Alternatively, in white-chinned petrels the Atlantic and Indian Ocean populations are not separated genetically, but New Zealand forms a separate cluster . To explain this variability, Friesen et al.  suggested that the extent of genetic structure in seabirds could be explained by the overlap of non-breeding areas, as most species with two or more population-specific non-breeding areas showed phylogeographic structure.
Due to their high mutation rates, microsatellite loci (i.e. regions across the genome where short nucleotide sequences are repeated in tandem) provide powerful tools in population genetics. They allow the testing of hypotheses related to within- and between-population genetic variability and the estimation of other demographic parameters such as gene flow and effective population size [13, 14]. Recently, Moodley et al.  used next-generation sequencing to develop microsatellite markers for thin-billed prions (Pachyptila belcheri), and other Pachyptila species (Aves, Procellariiformes; albatrosses and petrels). Here we use this newly developed marker set, together with tracking (geolocation) data to test the hypothesis that genetic structure is dependent on the non-breeding distribution pattern in thin-billed prions, Antarctic prions Pachyptila desolata and the closely-related blue petrel Halobaena caerulea. These small petrels have wide breeding distributions, including islands in different ocean basins, and are wide-ranging from subtropical to Antarctic waters [2, 16–19], making them ideal biological models to test for a correspondence between genetic structure and non-breeding distribution. The aim of the present study was thus to test the hypothesis that the extent of genetic structure is explained by the degree of overlap in non-breeding distributions of birds from different populations, by comparing three closely-related species with different migration strategies. We expected higher gene flow in the species in which populations showed greatest overlap during the non-breeding season.
Study species and sample collection
Blue petrels, thin-billed prions and Antarctic prions breed on islands in the Southern Ocean, and in each species, we sampled two large breeding populations, one in the south-west Atlantic Ocean (Falkland Islands and South Georgia) and the other in the Indian Ocean (Kerguelen Islands) . The three petrel species migrate away from their breeding grounds during the non-breeding season, where they segregate by choosing temperate (Antarctic prions), subantarctic (thin-billed prions) or polar (blue petrels) waters [2, 17–19]. Blue petrels, thin-billed prions and Antarctic prions show the typical procellariiform pattern of high longevity (several decades), and an extended period (several years) before they return to the colony to breed . Both sexes take part in incubation, and in feeding of their single, slow-growing chick [22, 23]. These three species are zooplanktivorous, with a preference for crustaceans [24, 25].
During 2010–2012, blood samples were obtained from adult birds breeding in colonies of thin-billed prions located at Île Mayès, Kerguelen (49°28′S, 69°57′E; N = 34) and New Island, Falkland/Malvinas Islands (51°43′S, 61°18′W; N = 42); Antarctic prions at Île Verte, Kerguelen (49°30′S, 70°02′E; N = 38) and Bird Island, South Georgia (54°00′S, 38°02′W; N = 35); and blue petrels at Île Mayès (N = 30) and Bird Island (N = 19). Blood (0.2–0.4 ml) was sampled by puncture of the wing vein and collection using heparinized capillaries, or using 0.3 ml syringes, and immediately transferred to FTA classic cards (Whatman International Ltd). These breeding colonies represent major breeding sites of all three species, with estimates at each site of several million birds for prions, and >100,000 individuals for the blue petrel [20, 26].
Analyses of distribution
Data on non-breeding distributions were obtained using 1 g leg-mounted geolocators, or Global Location Sensors (GLS loggers; model MK10, supplied by British Antarctic Survey and Biotrack UK), as described previously [2, 17, 18, 27]. Briefly, we attached geolocators to 12 to 25 individuals from each population during the breeding seasons 2009/10 (New Island, Bird Island) and 2011/12 (Kerguelen). Repeatability of migration routes across years is high in these species . The geolocators weighed 1 g, equivalent to <1% of the mean body mass of blue petrels (ca. 190 g), thin-billed prions (ca. 130 g) and Antarctic prions (ca. 160 g), and were fixed to plastic leg bands. Tagged individuals were fitted with standard metal rings on the other leg. Burrows were revisited and devices retrieved in the following two seasons (New I., N = 20 thin-billed prions; Bird I., N = 9 Antarctic prions, N = 11 blue petrels; Kerguelen, N = 10 Antarctic prions, N = 16 blue petrels and N = 19 thin-billed prions). We analysed distribution data from one non-breeding season per bird. Processing of geolocator data was carried out as described previously . Changes in distribution by month were examined using kernel analysis of filtered locations . The non-parametric fixed kernel density estimator was used to determine the 95% density contour. Kernel densities do not require serial independence of observations when estimating foraging ranges . Kernel analyses were performed in a Lambert azimuthal equal-area projection centred on the South Pole using ARCGIS 9.3 (ESRI, Redlands, CA, USA) using HAWTH’S ANALYSIS TOOLS . The overlap in 95% density polygons between populations were calculated for each species using the INTERSECT and UNION tools of ARCGIS 9.3.
Genetic variation and tests for selection and/or population expansion. Parameters are based on 18 polymorphic microsatellites recorded in Atlantic and Indian Ocean populations of the blue petrel Halobaena caerulea, thin-billed prion Pachyptila belcheri, and Antarctic prion P. desolata
Mitochondrial DNA cytochrome b 889 bp
7.2 ± 3.7
0.53 ± 0.24
0.66 ± 0.28
7.7 ± 4.9
0.52 ± 0.28
0.62 ± 0.30
8.4 ± 2.6
0.71 ± 0.10
0.74 ± 0.11
8.5 ± 2.8
0.74 ± 0.12
0.76 ± 0.11
8.5 ± 2.6
0.69 ± 0.12
0.75 ± 0.13
9.4 ± 3.5
0.71 ± 0.16
0.74 ± 0.16
Genetic structure was analysed using population-model based, and non-model based approaches. We used ARLEQUIN 3.5  to perform model-free analyses of molecular variance (AMOVAs) and tested for significance using 1000 permutations. As a model based analysis, we used the Bayesian clustering algorithm implemented in STRUCTURE 2.3.4 . We ran STRUCTURE under the admixture ancestry model and correlated allele frequencies without locality as priors. We explored values of K = 1–5 (10 iterations per K value), each consisting of 400,000 MCMC repetitions after a burn-in of 100,000. The most likely K value was determined using STRUCTURE HARVESTER  in accordance with the methods described in Evanno et al. . Different iterations of the optimal K value were combined in CLUMPP 1.1.2  and displayed using DISTRUCT 1.1 . Genetic structure was further analysed using factorial correspondence analysis (FCA) of population multilocus scores carried out in GENETIX 4.05 . Using the allele frequencies for all loci, FCA allows the visualisation in space of the genetic differences between individuals. We inferred gene flow in BAYESASS 3.0 using bi-directional migration rates (m) . This method was successful in the study of immigration and emigration in populations that are not in equilibrium e.g. [4, 43, 44]. BAYESASS estimates the posterior probability of an individual’s history and allows an estimation of the rate and direction of recent dispersal . The acceptance rates for the main parameters (i.e. ‘migration’ rate, inbreeding coefficient and allele frequencies) were adjusted during several preliminary runs. Convergence was assessed by checking the trace files in TRACER 1.6 . Final parameter estimates were obtained after performing three independent runs using different starting random seed numbers. The MCMC was run for 50,000,000 iterations with a burn-in period of 10,000,000 and a sampling frequency of 5000 iterations.
Genetic analyses using cytochrome b sequences
An 889-bp fragment of the cytochrome b (cyt b) gene was amplified using specific primers (CytB_Pri_F: 5′ -CTAGCTATACACTACACCGC-3′ and CytB_Pri_R: 5′ -CTAGTTGGCCGATGATGATG-3′) . PCRs were conducted in 20 μl reaction volumes containing 100 ng DNA template, 10 mM of each primer, 10 mM dNTPs (Roth, Karlsruhe), 2 mM MgCl and 0.5 U Taq DNA polymerase (BioLabs Taq DNA polymerase) in a 1× PCR reaction buffer. Thermocycling included initial denaturation at 94 °C for 2 min, 30 cycles of denaturation at 94 °C for 30 s, annealing at 60 °C for 45 s and extension at 72 °C for 1 min, followed by a final extension step of 5 min at 72 °C. PCR products were purified by digestion with exonuclease-shrimp alkaline phosphatase (from Fermentas Life Sciences), following the manufacturer’s specifications. PCR products were then sequenced in both directions using Big Dye chemistry (Applied Biosystems) and run on an AB 3130xl genetic analyser (for cyt b). Sequences were assembled and aligned in CLC Main Workbench 6.9.2 (CLC bio, Aarhus, Denmark).
A Fairy prion Pachyptila turtur sample was included as an outgroup. The best-fit model of nucleotide substitution for each of our sequence alignment was determined using Bayesian Information Criterion scores and corrected Akaike Information Criterion values in MEGA 6.06 . The cyt b trees were inferred by using the Maximum Likelihood method based on the Hasegawa-Kishino-Yano model  in MEGA, applying the Neighbor-Joining method for the heuristic search. Branch support was assessed with 10,000 pseudoreplicates (bootstrapping). Median-joining haplotype networks were estimated using NETWORK 5 . Genetic structure was also investigated by AMOVA, calculating the variance component distributed between populations (F ST) following Weir & Cockerham  in ARLEQUIN 3.5 . The significance of F ST values were tested with 1000 permutations.
Wintering range overlap
Population genetic structure – microsatellite data
The inbreeding coefficient (F IS ) based on 25 polymorphic microsatellites suggested a deviation from Hardy-Weinberg equilibrium due to an excess of homozygotes in all six populations. To test if this was due to non-random mating (either inbreeding or the existence of population genetic structure, i.e. the Wahlund effect ), or the existence of null alleles, we examined whether the effect was locus-specific. We found a significant deviation in seven of the 25 loci. Independently of the population analysed, always the same seven loci showed Hardy-Weinberg disequilibrium, which suggest the existence of null alleles. We therefore removed these seven loci from further analyses.
The genetic diversity found at the two major colonies of blue petrels, thin-billed prions, and Antarctic prions is summarised in Table 1. For most parameters considered, no differences were observed between populations from the Atlantic or Indian Oceans (Table 1). However, the inbreeding coefficient (F IS ) was higher in the Atlantic Ocean populations for all three species (Table 1).
In addition, classical estimators of population differentiation were very low between Atlantic and Indian Ocean populations for all species. AMOVA results indicated a higher population differentiation between the populations of Antarctic prions (F ST = 0.008, P < 0.001), than those of blue petrels (F ST = 0.001, P < 0.001), and thin-billed prions (F ST = −0.0002, P = 0.002).
Posterior mean migration rates and standard deviation of the marginal posterior distribution for each estimate. Mean migration rates (m) as a proportion from 0 to 1 and standard deviation (s.d.) were calculated between Atlantic and Indian Ocean populations of the blue petrel Halobaena caerulea, thin-billed prion Pachyptila belcheri, and Antarctic prion P. desolata
m[Atlantic] ← [Indian]
m[Indian] ← [Atlantic]
Population genetic structure – cytochrome b sequences
High between-population haplotype sharing was also inferred from maximum likelihood reconstructions of species-level phylogenies in blue petrels and Antarctic prions, where all splits included individuals from both colonies and with the highest bootstrap support being 69% (see Additional file 1: Figures S2 to S4). Again the cyt b phylogeny for thin billed prions appeared more structured, with a lower degree of haplotype sharing between populations.
Results of the AMOVA testing for genetic structure in the cytochrome b sequences. Test were carried out between Atlantic and Indian Ocean populations of the blue petrel Halobaena caerulea, thin-billed prion Pachyptila belcheri, and Antarctic prion P. desolata
Source of variation
Sum of squares
Percentage of variation
In the present study, we compare the population genetic structure of three small pelagic petrels in relation to their distribution during the non-breeding season. The three petrels differed in the degree of spatio-temporal segregation of birds from large populations in the Atlantic and Indian Oceans, ranging from no overlap in Antarctic prions, limited spatial overlap for a prolonged period (6 months, max. 17%) in thin-billed prions, to more pronounced overlap but for a shorter time (3 months, max. 56%) in blue petrels. Based on this and on the hypothesis of Friesen et al. , we would expect the least gene flow in the Antarctic prions.
In contrast, high levels of gene flow in all three species meant that these differences in relative segregation of non-breeding ranges were barely detectable through standard Bayesian and multivariate analyses of population genetic structure. Only model-free analyses of hierarchical genetic structure through AMOVA indicated low, but significant levels of population differentiation in both nuclear and mitochondrial marker sets between Atlantic and Indian Ocean populations of Antarctic prions, in line with expectations. In the blue petrel, which showed the highest spatial overlap between populations in the non-breeding season, there was very little between-population differentiation at both molecular levels, which was also in line with our hypotheses.
The lack of population genetic structure of the species studied here is similar to results observed in other Southern Ocean seabirds that feed in pelagic waters e.g. [8, 9]. The Procellariiformes have adapted to pelagic foraging, and only come ashore to breed, usually on remote islands . In the Southern Ocean, there are relatively few island groups suitable for breeding, and a pronounced genetic differentiation between populations might be expected based on the large distances and strong natal philopatry e.g. [52–54]. Alternatively, as these seabirds are highly mobile, they could maintain high levels of gene flow, as suggested by our results. The high migration rates also suggest that natal philopatry is not as strong as previously thought. In the breeding season, petrels are central-place foragers, constrained in their foraging distribution by the necessity to return to their nest site at regular intervals to incubate the egg or feed the chick. Even so, they can exploit vast ocean areas, foraging >1000 km from the nest [13, 28, 55]. In the non-breeding season, when not restricted by breeding duties, many seabird species, including the Procellariiformes, move over even larger distances, sometimes making circumpolar or trans-equatorial migrations [56–58]. The distribution of prions and blue petrels at sea has been the focus of recent stable isotope and tracking studies during the breeding and non-breeding season, which indicate considerable ecological segregation [2, 16, 17, 19, 59]. This segregation may have led to the speciation of prions and blue petrels via environmental specialisation (e.g. ) as reflected by genetic differences among species (Fig. 4). Other seabird species complexes breeding in different ecological (latitudinal) ranges in the Southern Ocean show a similar pattern (e.g. great albatrosses Diomedea spp. ; northern and southern rockhopper penguins Eudyptes moseleyi and E. chrysocome ).
Within our three study species, we found considerable admixture and evidence for high migration rates (Table 2). The lack of population structure even in Antarctic prions with non-overlapping wintering areas implies a higher gene flow than expected from the observed non-breeding distributions of adult birds. This may reflect the greater dispersal of young birds, which may spend several years roaming the oceans before settling to breed and may be much more wide-ranging even than adults during the non-breeding season . Ecological studies on the role of spatial structure in determining the dynamics of populations suggest that migration rates below 0.10 (10% of the individuals between the populations each generation) are indicative of populations behaving independently . In the three studied species, in addition to the high migration rates, BAYESASS analyses indicated a source-sink relationship between populations, with migration in one direction being much higher than in the opposite (Table 2). The reason for these asymmetries would require further detailed study; possible mechanisms include subtle differences in mating behaviour, low natal philopatry, or an active selection of the most productive colony by pre-breeders prior to recruitment. However, chick-provisioning rates observed for thin-billed prions do not support the latter hypothesis: although chicks on Kerguelen (Indian Ocean) are fed less frequently than those in the Falklands/Malvinas , the former received over 30% of immigrants per generation. To understand the reason for the observed asymmetries in the patterns of gene flow across ocean basins, a better understanding of the movements of juveniles and immatures prior to pair formation would be needed. However, prions only start breeding at 4–7 years old , and multi-year remote tracking is not yet feasible in these very small seabirds.
The mitochondrial data in the present study also allow an analysis of population demographic history. Negative Fu’s F S for thin-billed and Antarctic prions, as well as Tajima’s D for thin-billed prions at the Falklands and Antarctic prions at Kerguelen (Table 1), and the star-shaped networks may be considered as signatures of recent population expansions (Fig. 5). Natural climate cycles, such as past glacial/interglacial phases, shape species distributions. Many Southern Ocean species have undergone range restrictions and expansions associated with glacial cycles . The suitable glacial refugia for seabirds in the subantarctic during glacial maxima were restricted to a few lower-latitude islands (Gough Island, the Falkland Islands, or around New Zealand) [67, 68]. The patterns observed in the cytochrome b networks in the prions and blue petrels are compatible with a restriction to glacial refugia and subsequent population expansions when the species recolonized other subantarctic islands. A similar scenario has been postulated for other Southern Ocean seabirds, including king penguins Aptenodytes patagonicus , Emperor penguins Aptenodytes forsteri  and pygoscelid penguins .
In summary, our study highlights the low genetic structure of three abundant pelagic seabirds which show varying degrees of spatio-temporal segregation at sea between different breeding populations. In line with several other Southern Ocean seabirds, geographic distance between breeding colonies in these petrels does not generate genetic differences, regardless of the degree of overlap in non-breeding distributions. Thus, our study did not support the hypothesis of Friesen et al.  who suggested that the extent of genetic structure in seabirds could be explained by the overlap of non-breeding areas. Contrary to other seabird species with two or more population-specific non-breeding areas, the three small petrels in our study did not show population genetic structure. Some of the species investigated in  differ from this study either in being resident at, or near, breeding colonies or having passed through population bottlenecks. The reasons for other species differences (e.g. differing degrees of philopatry) need further investigation.
The authors thank A. Corbeau, J. Ferrer-Obiol, M. Passerault and T. Lacombe for assistance in the field at Kerguelen. We are also grateful to Ruth Brown for assistance with blood-sampling, device deployment and retrieval, to Gopi K. Munimanda for assistance with lab work, and to Theresa L. Cole for lab work and sample genotyping.
Funding was obtained through a grant of the Deutsche Forschungsgemeinschaft (Germany), Heisenberg program (DFG, Qu 148-5). Fieldwork was supported through a CGS award from the Natural Environment Research Council (NERC) and at Kerguelen by the Institut Polaire Français Paul Emile Victor (Programme N°109, HW) and by Terres Australes et Antarctiques Françaises.
Availability of data and materials
All relevant data are within the paper and its Additional file. DNA sequences can be accessed via GenBank.
PQ, JFM and YM conceived and designed the study. HW, YC, KD, RP, JN and JFM carried out the extensive fieldwork. JFM and LC carried out the lab work. PQ, RP and JN analysed the tracking data. JFM carried out the bioinformatics analyses. PQ and JFM drafted the manuscript. All authors reviewed the final draft of the manuscript. All authors read and approved the final manuscript.
The study involved free-ranging animals and was carried out under licence obtained by the Falkland Islands Government and Terres Australes et Antarctiques Françaises.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Knox GA. Biology of the southern ocean. Boca Raton: CRC Press; 2006.Google Scholar
- Quillfeldt P, Cherel Y, Delord K, Weimerkirch H. Cool, cold or colder? Spatial segregation of prions and blue petrels is explained by differences in preferred sea surface temperatures. Biol Lett. 2015;11(4):20141090.View ArticlePubMedPubMed CentralGoogle Scholar
- Ropert-Coudert Y, Hindell MA, Phillips R, Charrassin J-B, Trudelle L, Raymond B. Biogeographic patterns of birds and mammals. The Biogeographic Atlas of the Southern Ocean Scientific Committee on Antarctic Research. 2014:364–87. http://atlas.biodiversity.aq. Accessed 30 June 2017.
- Friesen V, Burg T, McCoy K. Mechanisms of population differentiation in seabirds. Mol Ecol. 2007;16(9):1765–85.View ArticlePubMedGoogle Scholar
- Friesen VL. Speciation in seabirds: why are there so many species… and why aren’t there more? J Ornithol. 2015;156(Suppl 1):S27–39.View ArticleGoogle Scholar
- Levy H, Clucas GV, Rogers AD, Leaché AD, Ciborowski KL, Polito MJ, Lynch HJ, Dunn MJ, Hart T. Population structure and phylogeography of the Gentoo Penguin (Pygoscelis papua) across the Scotia Arc. Ecol Evol. 2016;6:1834–53.Google Scholar
- Roeder AD, Marshall RK, Mitchelson AJ, Visagathilagar T, Ritchie PA, Love DR, Pakai TJ, McPartlan HC, Murray ND, Robinson NA. Gene flow on the ice: genetic differentiation among Adélie penguin colonies around Antarctica. Mol Ecol. 2001;10(7):1645–56.View ArticlePubMedGoogle Scholar
- Burg T, Croxall J. Global relationships amongst black-browed and grey-headed albatrosses: analysis of population structure using mitochondrial DNA and microsatellites. Mol Ecol. 2001;10:2647–60.View ArticlePubMedGoogle Scholar
- Milot E, Weimerskirch H, Bernatchez L. The seabird paradox: dispersal, genetic structure and population dynamics in a highly mobile, but philopatric albatross species. Mol Ecol. 2008;17:1658–73.View ArticlePubMedGoogle Scholar
- Techow NMSM, O’Ryan C, Phillips RA, Gales R, Marin M, Patterson-Fraser D, Quintana F, Ritz MS, Thompson DR, Wanless RM, et al. Speciation and phylogeography of giant petrels Macronectes. Mol Phylogenet Evol. 2010;54(2):472–87.Google Scholar
- Ritz MS, Millar C, Miller GD, Phillips RA, Ryan P, Sternkopf V, Liebers-Helbig D, Peter H-U. Phylogeography of the southern skua complex—rapid colonization of the southern hemisphere during a glacial period and reticulate evolution. Mol Phylogenet Evol. 2008;49:292–303.View ArticlePubMedGoogle Scholar
- Techow NMSM, Ryan PG, O'Ryan C. Phylogeography and taxonomy of White-chinned and Spectacled Petrels. Mol Phylogenet Evol. 2009;52:25–33.View ArticleGoogle Scholar
- Bruford MW, Wayne RK. Microsatellites and their application to population genetic studies. Curr Opin Genet Dev. 1993;3(6):939–43.View ArticlePubMedGoogle Scholar
- Sunnucks P. Efficient genetic markers for population biology. Trends Ecol Evol. 2000;15(5):199–203.View ArticlePubMedGoogle Scholar
- Moodley Y, Masello JF, Cole TL, Calderon L, Munimanda GK, Thali MR, Alderman R, Cuthbert RJ, Marin M, Massaro M, et al. Evolutionary factors affecting the cross-species utility of newly developed microsatellite markers in seabirds. Mol Ecol Resour. 2015;15(5):1046–58.View ArticlePubMedGoogle Scholar
- Navarro J, Votier SC, Aguzzi J, Chiesa JJ, Forero MG, Phillips RA. Ecological segregation in space, time and trophic niche of sympatric planktivorous petrels. PLoS One. 2013;8(4):e62897.View ArticlePubMedPubMed CentralGoogle Scholar
- Navarro J, Cardador L, Brown R, Phillips RA. Spatial distribution and ecological niches of non-breeding planktivorous petrels. Sci Rep. 2015;5:12164.View ArticlePubMedPubMed CentralGoogle Scholar
- Quillfeldt P, Masello JF, Navarro J, Phillips RA. Year-round distribution suggests spatial segregation of two small petrel species in the South Atlantic. J Biogeogr. 2013;40:430–41.View ArticleGoogle Scholar
- Cherel Y, Quillfeldt P, Delord K, Weimerskirch H. Combination of at-sea activity, geolocation and feather stable isotopes documents where and when seabirds moult. Front Ecol Evol. 2016;4Google Scholar
- Weimerskirch H, Zotier R, Jouventin P. The avifauna of the Kerguelen Islands. Emu. 1989;89:15–29.View ArticleGoogle Scholar
- Barbraud C, Weimerskirch H. Climate and density shape population dynamics of a marine top predator. Proc R Soc Lond B Biol Sci. 2003;270(1529):2111–6.View ArticleGoogle Scholar
- Chaurand T, Weimerskirch H. Incubation routine, body mass regulation and egg neglect in the blue petrel Halobaena caerulea. Ibis. 1994;136(3):285–90.Google Scholar
- Quillfeldt P, Strange I, Segelbacher G, Masello JF. Male and female contributions to provisioning rates of Thin-billed prions Pachyptila belcheri in the South Atlantic. J Ornithol. 2007;38:298–308.Google Scholar
- Cherel Y, Bocher P, De Broyer C, Hobson KA. Food and feeding ecology of the sympatric thin-billed Pachyptila belcheri and Antarctic P. desolata prions at Iles Kerguelen, Southern Indian Ocean. Mar Ecol Prog Ser. 2002;228:263–81.Google Scholar
- Cherel Y, Bocher P, Trouvé C, Weimerskirch H. Diet and feeding ecology of blue petrels Halobaena caerulea at Iles Kerguelen, Southern Indian Ocean. Mar Ecol Prog Ser. 2002;228:283–99.Google Scholar
- Barbraud C, Delord K. Population census of blue petrels Halobaena caerulea at Mayes Island, Iles Kerguelen. Antarct Sci. 2006;18(02):199–204.Google Scholar
- Quillfeldt P, Cherel Y, Masello JF, Delord K, McGill RAR, Furness RW, Moodley Y, Weimerskirch H. Half a world apart? Overlap in nonbreeding distributions of Atlantic and Indian Ocean thin-billed prions. PLoS One. 2015;10(5):e0125007.View ArticlePubMedPubMed CentralGoogle Scholar
- Phillips RA, Silk JRD, Phalan B, Catry P, Croxall JP. Seasonal sexual segregation in two Thalassarche albatross species: competitive exclusion, reproductive role specialization or foraging niche divergence? Proc R Soc B. 2004;271:1283–91.Google Scholar
- de Solla SR, Bonduriansky R, Brooks RJ. Eliminating autocorrelation reduces biological relevance of home range estimates. J Anim Ecol. 1999;68:221–34.View ArticleGoogle Scholar
- Beyer HL. Hawth's Analysis Tools for ArcGIS. 2004.Google Scholar
- Merino S, Hennicke J, Martínez J, Ludynia K, Torres R, Work TM, Stroud S, Masello JF, Quillfeldt P. Infection by Haemoproteus parasites in four species of frigatebirds and the description of a new species of Haemoproteus (Haemosporida: Haemoproteidae). J Parasitol. 2012;98:388–97.View ArticlePubMedGoogle Scholar
- van Oosterhout C, Hutchinson WF, Wills DP, Shipley P. Micro-checker: software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol Notes. 2004;4:535–8.View ArticleGoogle Scholar
- Dieringer D, Schlötterer C. Microsatellite analyser (MSA): a platform independent analysis tool for large microsatellite data sets. Mol Ecol Notes. 2003;3:167–9.View ArticleGoogle Scholar
- Excoffier L, Lischer HEL. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour. 2010;10:564–7.View ArticlePubMedGoogle Scholar
- Raymond M, Rousset F. GENEPOP Version 3.1d: population genetics software for exact tests and ecumenism. J Hered. 1995;86:248–9.View ArticleGoogle Scholar
- Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155(2):945–59.PubMedPubMed CentralGoogle Scholar
- Earl DA, von Holdt BM. STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv Genet Resour. 2012;4:359–61.View ArticleGoogle Scholar
- Evanno G, Regnaut S, Goudet J. Detecting the number of clusters of individuals using the software Structure: a simulation study. Mol Ecol. 2005;14(8):2611–20.View ArticlePubMedGoogle Scholar
- Jakobsson M, Rosenberg NA. CLUMPP: a cluster matching and permutation program for dealing with label switching and multimodality in analysis of population structure. Bioinformatics. 2007;23(14):1801–6.View ArticlePubMedGoogle Scholar
- Rosenberg NA. DISTRUCT: a program for the graphical display of population structure. Mol Ecol Notes. 2004;4(1):137–8.View ArticleGoogle Scholar
- Belkhir K, Borsa P, Chikhi L, Raufaste N, Bonhomme F. GENETIX 4.05, logiciel sous Windows TM pour la génétique des populations. 4.05 ed. Montpellier: Laboratoire Génome, Populations, Interactions, CNRS UMR 5171, Université de Montpellier II; 2004.Google Scholar
- Wilson GA, Rannala B. Bayesian inference of recent migration rates using multilocus genotypes. Genetics. 2003;163(3):1177–91.PubMedPubMed CentralGoogle Scholar
- Genovart M, Thibault J-C, Igual JM, MdM B-R, Rabouam C, Bretagnolle V. Population structure and dispersal patterns within and between Atlantic and Mediterranean populations of a large-range pelagic seabird. PLoS One. 2013;8:e70711.View ArticlePubMedPubMed CentralGoogle Scholar
- Calderón L, Quintana F, Cabanne GS, Lougheed SC, Tubaro PL. Phylogeography and genetic structure of two Patagonian shag species (Aves: Phalacrocoracidae). Mol Phylogenet Evol. 2014;72:42–53.View ArticlePubMedGoogle Scholar
- Rambaut A, Drummond AJ. Tracer v1.4. 2007.Google Scholar
- Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: Molecular Evolutionary Genetics Analysis Version 6.0. Mol Biol Evol. 2013;30(12):2725–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Hasegawa M, Kishino H, Yano T. Dating of the human-ape splitting by a molecular clock of mitochondrial DNA. J Mol Evol. 1985;22:160–74.View ArticlePubMedGoogle Scholar
- Bandelt HJ, Forster P, Röhl A. Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999;16:37–48.View ArticlePubMedGoogle Scholar
- Weir BS, Cockerham CC. Estimating F-statistics for the analysis of population structure. Evolution. 1984;38:1358–70.PubMedGoogle Scholar
- Koffi M, De Meeûs T, Séré M, Bucheton B, Simo G, Njiokou F, Salim B, Kaboré J, MacLeod A, Camara M, et al. Population genetics and reproductive strategies of African trypanosomes: revisiting available published data. PLoS Negl Trop Dis. 2015;9(10):e0003985.View ArticlePubMedPubMed CentralGoogle Scholar
- Brooke M. Albatrosses and petrels across the world. Oxford: Oxford University Press; 2004.Google Scholar
- Ovenden J, Wust-Saucy A, Bywater R, Brothers N, White R. Genetic evidence for philopatry in a colonially nesting seabird, the Fairy Prion (Pachyptila turtur). Auk. 1991;108:688–94.Google Scholar
- Steeves TE, Anderson DJ, Friesen VL. A role for nonphysical barriers to gene flow in the diversification of a highly vagile seabird, the masked booby (Sula dactylatra). Mol Ecol. 2005;14:3877–87.View ArticlePubMedGoogle Scholar
- Bicknell A, Knight M, Bilton D, Reid J, Burke T, Votier S. Population genetic structure and long-distance dispersal among seabird populations: Implications for colony persistence. Mol Ecol. 2012;21(12):2863–76.View ArticlePubMedGoogle Scholar
- Quillfeldt P, Phillips RA, Marx M, Masello JF. Colony attendance and at-sea distribution of thin-billed prions during the early breeding season. J Avian Biol. 2014;45(4):315–24.View ArticleGoogle Scholar
- Rayner MJ, Taylor GA, Gummer HD, Phillips RA, Sagar PM, Shaffer SA, Thompson DR. The breeding cycle, year-round distribution and activity patterns of the endangered Chatham Petrel (Pterodroma axillaris). Emu. 2012;112(2):107–16.Google Scholar
- Shaffer SA, Tremblay Y, Weimerskirch H, Scott D, Thompson DR, Sagar PM, Moller H, Taylor GA, Foley DG, Block BA, et al. Migratory shearwaters integrate oceanic resources across the Pacific Ocean in an endless summer. Proc Natl Acad Sci. 2006;103(34):12799–802.View ArticlePubMedPubMed CentralGoogle Scholar
- Weimerskirch H, Delord K, Guitteaud A, Phillips RA, Pinet P. Extreme variation in migration strategies between and within wandering albatross populations during their sabbatical year, and their fitness consequences. 2015;5:8853.Google Scholar
- Quillfeldt P, Masello JF, McGill RAR, Adams M, Furness RW. Moving polewards in winter: a recent change in migratory strategy. Front Zool. 2010;7:15.11.View ArticleGoogle Scholar
- Sexton JP, Hangartner SB, Hoffmann AA. Genetic isolation by environment or distance: which pattern of gene flow is most common? Evolution. 2014;68(1):1–15.View ArticlePubMedGoogle Scholar
- Burg T, Croxall J. Global population structure and taxonomy of the wandering albatross species complex. Mol Ecol. 2004;13:2345–55.View ArticlePubMedGoogle Scholar
- de Dinechin M, Ottvall R, Quillfeldt P, Jouventin P. Speciation chronology of Rockhopper penguins inferred from molecular, geological and palaeoceanographic data. J Biogeogr. 2009;36:693–702.View ArticleGoogle Scholar
- Thiers L, Delord K, Barbraud C, Phillips RA, Pinaud D, Weimerskirch H. Foraging zones of the two sibling species of giant petrels in the Indian Ocean throughout the annual cycle: implication for their conservation. Mar Ecol Prog Ser. 2014;499:233–48.View ArticleGoogle Scholar
- Hastings A. Complex interactions between dispersal and dynamics: lessons from coupled logistic equations. Ecology. 1993;74(5):1362–72.View ArticleGoogle Scholar
- Quillfeldt P, Strange I, Masello JF. Sea surface temperatures, variable food supply and behavioural buffering capacity in Thin-billed prions Pachyptila belcheri: breeding success, provisioning and chick begging. J Avian Biol. 2007;38:298–308.View ArticleGoogle Scholar
- Schreiber EA, Burger J. Biology of marine birds. Boca Raton: CRC Press; 2001.Google Scholar
- Fraser CI, Nikula R, Ruzzante DE, Waters JM. Poleward bound: biological impacts of Southern Hemisphere glaciation. Trends Ecol Evol. 2012;27(8):462–71.View ArticlePubMedGoogle Scholar
- Trucchi E, Gratton P, Whittington JD, Cristofari R, Le Maho Y, Stenseth NC, Le Bohec C. King penguin demography since the last glaciation inferred from genome-wide data. Proc R Soc B, The Royal Society. 2014:20140528.Google Scholar
- Younger JL, Clucas GV, Kooyman G, Wienecke B, Rogers AD, Trathan PN, Hart T, Miller KJ. Too much of a good thing: Sea ice extent may have forced emperor penguins into refugia during the last glacial maximum. Glob Chang Biol. 2015;21(6):2215–26.View ArticlePubMedGoogle Scholar
- Clucas GV, Dunn MJ, Dyke G, Emslie SD, Naveen R, Polito MJ, Pybus OG, Rogers AD, Hart T. A reversal of fortunes: climate change 'winners' and 'losers' in Antarctic Peninsula penguins. Sci Rep. 2014;4:5024.View ArticlePubMedPubMed CentralGoogle Scholar