Being cosmopolitan: evolutionary history and phylogeography of a specialized raptor, the Osprey Pandion haliaetus
© Monti et al. 2015
Received: 27 July 2015
Accepted: 9 November 2015
Published: 17 November 2015
The Osprey (Pandion haliaetus) is one of only six bird species with an almost world-wide distribution. We aimed at clarifying its phylogeographic structure and elucidating its taxonomic status (as it is currently separated into four subspecies). We tested six biogeographical scenarios to explain how the species’ distribution and differentiation took place in the past and how such a specialized raptor was able to colonize most of the globe.
Using two mitochondrial genes (cyt b and ND2), the Osprey appeared structured into four genetic groups representing quasi non-overlapping geographical regions. The group Indo-Australasia corresponds to the cristatus ssp, as well as the group Europe-Africa to the haliaetus ssp. In the Americas, we found a single lineage for both carolinensis and ridgwayi ssp, whereas in north-east Asia (Siberia and Japan), we discovered a fourth new lineage. The four lineages are well differentiated, contrasting with the low genetic variability observed within each clade. Historical demographic reconstructions suggested that three of the four lineages experienced stable trends or slight demographic increases. Molecular dating estimates the initial split between lineages at about 1.16 Ma ago, in the Early Pleistocene.
Our biogeographical inference suggests a pattern of colonization from the American continent towards the Old World. Populations of the Palearctic would represent the last outcomes of this colonization. At a global scale the Osprey complex may be composed of four different Evolutionary Significant Units, which should be treated as specific management units. Our study brought essential genetic clarifications, which have implications for conservation strategies in identifying distinct lineages across which birds should not be artificially moved through exchange/reintroduction schemes.
The modern distribution of living organisms has been shaped by multiple processes that had profound effects on the dispersal, genetic structure and evolutionary histories of plant and animal populations. Movements of land-masses and successive multiple glacial events that occurred during the Pleistocene caused severe habitat changes which confined many species to warmer refugia and led other taxa to experience demographic reductions or complete extinction . Favourable periods during climatic fluctuations allowed successive population expansions, together with the recolonization of portions of the ancient ranges .
Despite the high potential dispersive power of flying birds, it is striking that only a few taxa did colonize most of the world. Excluding seabirds, for which the distribution pattern depends more on ocean basins than on the major land-masses , only six landbird species (out of ca. 10,000 species) are known to be cosmopolitan, breeding in each biogeographical region of the world, except Antarctica. This group includes the Great Egret Ardea alba, the Cattle Egret Bubulcus ibis, the Glossy Ibis Plegadis falcinellus, the Barn Owl Tyto alba, the Peregrine Falcon Falco peregrinus and the Osprey Pandion haliaetus.
The Osprey is a medium-sized raptor with flexible breeding habitat requirements across its range. Despite its specialization as a piscivore, it is an opportunistic forager that can feed in both freshwater and marine environments . Northern populations are known to migrate long distances [4, 5], whereas individuals from lower latitudes (e.g. Caribbean, Atlantic islands and Mediterranean basin) are mostly sedentary, or perform small-scale interbreeding movements [4, 6]. One could therefore predict that such broad habitat tolerances and high mobility should homogenize genomes, limiting genetic differentiation across populations at a continental level, as described in other widespread raptors (e.g. Haliaeetus albicilla: ; Falco peregrinus: ). However, adult ospreys tend to return to their natal area to breed  and such strong philopatry may have played in favour of genetic structuring among populations across the extensive range. Similarly, genetic differences may be expected between long-distance migratory and partially migratory/resident populations.
In this context, using molecular markers is a powerful method for inferring the evolutionary history of the Osprey. Only two genetic studies have been carried out on this species [14, 15], but they did not investigate all subspecies, in the entire distributional range. Here, we carried out an extensive phylogeographic study based on mitochondrial DNA sequences (cytochrome b and ND2 genes) from samples covering the worldwide distribution of the Osprey. This exhaustive dataset allowed addressing specific questions. a) Does the Osprey show any phylogeographic structure in mitochondrial DNA at the continental level? b) How was such a specialized raptor able to colonize the entire world? We propose a hypothetical evolutionary scenario explaining how the species’ distribution and differentiation took place in the past. In the light of these new genetic clarifications, we discuss the potential implications for revisions of the taxonomy and for conservation (e.g. defining Evolutionary Significant Units (ESUs; ) to design relevant conservation strategies at the appropriate spatial scale.
Sampling, DNA sequencing and alignment
Fresh samples (n = 114; both blood and feathers) were obtained from wild ospreys at the nest during ringing activities in different breeding sites. For each individual, about 0.5 ml of blood was taken by venepuncture from the wing and stored either on filter blotting papers or in 70 % ethanol in Eppendorf tubes. In some cases, growing body feathers containing traces of blood within the calamus were collected and stored in envelopes.
For areas where it was not possible to collect fresh samples, we completed our sampling using 90 museum specimens. We collected small pieces of skin (about 2 mm from the toepad) from museum study Osprey-skins following the procedure described by Mundy et al. . From each museum, we mainly choose museum study skins of various origins and collected during the breeding season, so excluding possible vagrants, wintering birds or dispersing animals.
From both fresh and museum specimens DNA was extracted and amplified by PCR (see technical procedures below) for the mitochondrial Cytochrome b (cyt b) and NADH dehydrogenase subunit 2 (ND2). For the cyt b, five sequences from the GenBank database were also included, leading to a total of 209 sequences (see Additional file 1). For ND2 a subset of only 37 individuals returned high-quality PCR products to which we added one sequence from GenBank, thus leading to a dataset of 38 individuals (Additional file 1). Finally, four other raptor species were used as outgroups: White-tailed eagle (Haliaeetus albicilla), Tawny eagle (Aquila rapax), Black-breasted buzzard (Hamirostra melanosternon) and Gray-headed kite (Leptodon cayanensis). The choice of these outgroups was motivated because 1) both cyt b and ND2 sequences were available in GenBank (see  for outgroup accession numbers), and 2) they represent four lineages of the Accipitridae family, the sister taxa of Pandionidae (including only Pandion haliaetus; ). The number of individuals and the length of sequences used for each analysis are reported in supplementary information (Additional file 2).
Total genomic DNA was extracted from fresh and historical samples using a Qiagen DNeasy Tissue kit, following the manufacturer’s instructions (Qiagen, Hilden, Germany). To avoid contamination with contemporary DNA , all extractions from museum specimens were performed using the facilities of the platform “degraded DNA” (Labex CeMEB, Montpellier, France) dedicated to degraded DNA experiments, where we adopted the following specific safety measures. Equipment, consumables and Qiagen DNeasy Tissue kits used in the platform were purchased new, while the room was regularly cleaned and exposed to UV overnight after each DNA extraction cycle, in order to destroy possible traces of DNA between successive extractions. Experimentators wore protective clothing and footwear. As a further precaution and following Bantock et al. , we worked with a maximum of 12 samples during each series of extraction to reduce the risks of cross-contamination and possible handling errors with tubes. We only used half of each foot-skin sample (about 10 mg of tissue) from which total DNA was extracted. Samples were incubated at least for one night at 56 °C to be digested during the lysis. The tissue was digested in 180 μL buffer ATL/20 μL proteinase K solution for 20-h at 55 °C; other reagents and the spin column were used according to the manufacturer's instructions ("Tissue protocol"), and final DNA elution was performed with 2 × 100 μL of 10 mM Tris, 0.5 mM HC1 pH 9.0 preheated to 70 °C. Multiple negative extraction and amplification controls were carried out simultaneously, using the same instruments and reagents, to detect possible contamination.
Portions of the mitochondrial cyt b and ND2 were amplified by PCR. Specific external and internal mitochondrial cyt b and ND2 primers were designed in this study for Pandion haliaetus (Additional file 3). PCR was performed using 1-μl (fresh samples) to 3-μl (museum specimens) of total DNA extracted, 5-μl of the Qiagen Multiplex PCR kit (containing HotStarTaq DNA polymerase, DNTPs and Multiplex PCR buffer), 1-μl for each primer at 2 pm and 2-μl of purified water. PCR reaction was performed using a MasterCycler Eppendorf thermocycler and began with an initial denaturation of 15 min at 95 °C, followed by 30 cycles of 30 s denaturation at 94 °C, 90 s annealing at 54 °C, 1 min extension at 72 °C and a 30 min final extension at 60 °C. A mitochondrial cytochrome b 1040 nucleotides fragment was amplified with PANHF1 and PANHR5 primers; F13 and PHND2-R1 primers were used to amplify a ND2 nucleotide fragment of 1100 bp. In case of degraded DNA, we used internal primers to amplify cyt b and ND2 in 300 to 500 nucleotides overlapping fragments (see Additional file 3). Screening of the PCR products was performed by running on a 1 % agarose gel using GelRed TM nucleic Acid gel stain (Biotium). Size products have been compared to long fragments ladder from Eurogentec, Smart LadderTM. After band sizes were determined, PCR products were sequenced at the Genoscope (Evry, France).
Electrophoregrams were read using CodonCode Aligner 4.0.4 software and sequences were aligned by eye using Seaview 4 software . Sequences were also translated into amino acids to check for any stop codons and possible amplification of pseudo-genes. Consensus sequences obtained for cyt b and ND2 from both fresh and museum samples were deposited in GenBank under accession numbers given in Additional file 1.
Partitioning and phylogenetic analyses
Phylogenetic relationships were inferred from the cyt b alone or from the concatenated cyt b + ND2 datasets. We determined both the best-fit partitioning scheme and the best models of sequence evolution using PartitionFinder 1.0.1 .
Phylogenetic trees were reconstructed using two probabilistic methods: Bayesian inference (BI) and maximum likelihood (ML). Bayesian analyses were performed with MrBayes 3.2  using the partitioning strategy obtained with PartitionFinder (see also Results). Two separate runs of five million generations (sampled every 250 generations) were conducted simultaneously. Tracer 1.5  was used to check the convergence between the two runs and to determine the burn-in period. On this basis, the first 2000 phylogenetic trees were discarded (10 %), and the remaining 18000 trees were used to estimate posterior parameters and probability distributions. ML tree was constructed with RaxMl 8.0.17 . As GTR is the only nucleotide substitution model available in RaxMl, GTR + G was applied to all partitions previously determined by PartitionFinder. The robustness of nodes was evaluated with 1000 bootstrap replicates with the option –b. The consensus tree was obtained using the program Consense of the Phylip 3.69 package .
Relationships between haplotypes were also visualized as a minimum spanning network, using the Median-Joining (MJ) network algorithm implemented in the program network 18.104.22.168 . In order to avoid artefactual groupings due to missing data, the MJ network was built considering the most complete dataset in terms of nucleotides and individuals, which is a fragment of 661 bp of the cyt b on 146 samples.
Genetic diversity, demographic history and molecular dating
We used Dnasp 5.10  to compute the number of haplotype (nH), haplotype diversity (H), nucleotide diversity (π) as well as the average number of nucleotide differences (k). Mean genetic distances within and between groups were computed using the p-distance and a pairwise deletion for the gaps/missing data treatment, as implemented in the Mega 5.10 software .
Demographic history of the haplogroups and the whole dataset was determined with different methods. Firstly, R2 , Fu’s Fs  statistics and their significance were calculated with DnaSP. Ramos-Onsins & Rozas  recommended using R2 when population sizes are small (~10) and Fs when sample sizes are large (~50). The historical demography of main haplogroups was also estimated based on the cyt b dataset using the skyline plot method (BSP; [30, 31] implemented in Beast 1.8.0. BSP analyses were performed on each group (including all cyt b sequences) separately with the cyt b partitioned according to codon position and using the HKY + G model as sequence evolution. The likelihood-ratio test performed with tree-puzzle 5.2  rejected the strict molecular clock hypothesis (p < 0.05). BSP analyses were thus conducted using a lognormal-relaxed molecular clock with a substitution rate of 0.01973 per lineage per million years as estimated by Nabholz et al.  for the Osprey cytochrome b. Although this rate has been obtained for the third codon position we apply it because most substitutions (83.3 %) are located at this position in our dataset. Analyses were run for 50 million generations, sampled every 1000 generations, after discarding the first 10 % as burn-in. We used Tracer 1.5 to analyse the results and draw the BSPs.
Time of the most recent common ancestor (TMRCA) was estimated with Beast 1.8.0 based on 23 sequences (19 Osprey haplotypes and 4 outgroups) of the cytochrome b (661 positions). The whole alignment was partitioned according to the three codon positions using a HKY + G model of sequence evolution. As previously, analyses were performed using a lognormal-relaxed molecular clock using a substitution rate (ucld. mean) following a normal prior distribution of mean 0.01973  and a standard deviation of 0.005 to take into account some rate uncertainty. Two runs were performed, each of 50 million generations, sampled every 1000 generations, and a 10 % burn-in was applied. The resulting tree files were combined with LogCombiner 1.8.0 and the maximum clade credibility tree (mean height) was obtained with Tree-Annotator 1.8.0.
Historical biogeography reconstruction
Probabilistic inference of ancestral range was performed using BioGeoBEARS  as implemented in R. We performed six different models [34, 35], to reconstruct the biogeographic history of Osprey across continents. Each analysis allows for a different subset of biogeographical possibilities, such as dispersal, vicariance and extinction (see [34, 35]). Accordingly, we tested the following models: Dispersal-Extinction Cladogenesis Model (DEC), Dispersal-Vicariance Analysis (DIVA), Bayesian inference of historical biogeography for discrete areas (BayArea) and the same models including the founder-event speciation parameter (‘j’), DECj, DIVAj and BayAreaj. Finally, statistical fit of the six different models were compared using a model choice procedure by means of the Akaike Information Criterion implemented in the R package BioGeoBEARS [34, 35].
For the mitochondrial cyt b dataset (213 individuals = 209 ospreys + 4 outgroups), the best-fit scheme was a partitioning according to codon position with the models HKY + I for position 1 (367 nucleotides), TrN for position 2 (368) and K80 + G for position 3 (368). As the TrN substitution model was not available in MrBayes, the parameter Nst was set to 6. The partitioned ML analysis was performed with 1000 bootstrap replicates using a GTR + G substitution model for each codon position in RaxMl software. The average Bayesian posterior probabilities (pp) and bootstrap values (BP) for supported clades are shown on the tree in Fig. 2. The phylogenetic tree including the totality of individuals (n = 209) is reported in supplementary information (Additional file 4).
The second cluster (IND-AUS; pp = 1.00, BP = 100) is composed of 3 haplotypes including 16 individuals retrieved from the Indo-Australasian area: 6 from various islands of Indonesia and 10 from Australia.
A third group was composed of 2 haplotypes and 11 ospreys from Asia (ASIA; pp = 1.00, BP = 71): 6 from the Pacific coast of Siberia (e.g. regions of Magadan, Khabarovsky, Primorskii and the Kurile islands) and five from Japan (Fig. 3). In addition, two samples that belonged to this haplogroup were collected in other geographic areas: one from continental Asia (Mongolia) and another one from the Indo-Australasian region (New Guinea).
Finally, the largest clade (EUR-AFR; pp = 1, BP = 96) was formed by 102 individuals representing 9 haplotypes, mainly belonging to the Western Palearctic area, with a few exceptions. 93 of these samples were collected in Europe along a latitudinal gradient scattered from northern Europe (Fenno-Scandia and western Russia), central Europe (Germany, France), to southern localities in the Mediterranean area (Corsica, Balearics, Italy, Portugal). Samples from North African coasts (e.g. Morocco) and Atlantic islands (e.g. Canary and Cape Verde) were included in this haplogroup, together with ospreys from the Red Sea and Persian Gulf areas. Three single haplotypes identified particular geographical regions (Fig. 3): Canary Islands (4 individuals), Persian Gulf (2 ind.) and Red Sea (3 ind.). The remaining samples of the EUR-AFR group (93 ind.) mostly belonged to two frequent haplotypes. Interestingly, three geographical exceptions were recorded in this group: two individuals from central Siberia (Tuva and Baikal regions) and one from India.
In total, three potential mixing areas were detected between phylogenetic lineages in the Old world (see Fig. 7): a) one in central Siberia between EUR-AFR and ASIA; b) one in Indonesia between ASIA and IND-AUS and c) a third one between west Indonesia and India between EUR-AFR and IND-AUS.
Uncorrected pairwise (p)-distance expressed as percentage (%) for cyt b within (in bold) and between clades in Pandion haliaetus. The number of individuals analzed for each group is given into parentheses
Relationships between the four haplogroups did not appear well resolved. The structure of the cyt b phylogenetic tree (Fig. 2a) revealed that the three Old World lineages (EUR-AFR, ASIA and IND-AUS) formed a rather well supported grouping (pp = 0.99, BP = 66) leaving the AMER group as the basal clade. On the other hand, the sister group relationships between EUR-AFR and ASIA is poorly supported (pp = 0.33, BP = 35). In order to improve resolution between groups we sequenced the ND2 gene for a subsample of 39 individuals (Additional file 4b). The concatenation of the two mitochondrial fragments represented 2037-bp and provided fourteen haplotypes from thirty-eight sequences (Additional file 1). Four partitions were obtained for the best-fit scheme: one for the position 1 of cyt b (with the model K80 + G), one for the cyt b-position 2 and ND2-position 3 (model HKY + I), one for the cyt b-position 3 and the ND2-position1 (model TrN) and one for the ND2-position 2 (model K81uf + G). As previously stated, the TrN model was approximated with Nst = 6 in MrBayes. The analysis carried out on the two genes (Fig. 2b for the 14 haplotypes and Additional file 4c for the 38 individuals) reinforced the support for the four main lineages, although it should be noted that the ASIA group was now represented only by five samples from Japan. Even though the number of nucleotides has been doubled, the only noticeable gain is a slight increased support for the node EUR-AFR/ASIA (pp = 0.66, BP = 47) whereas other relationships, including the position of the root, remained unresolved.
Network, genetic variability and demography
The cyt b network (Fig. 3) confirmed that four major groups which were included in 19 unique mtDNA haplotypes. The EUR-AFR clade (n = 102) resulted in nine haplotypes differing by only one nucleotide change. Two out of the nine haplotypes were most frequent, and shared by the majority of the individuals (58 and 25 individuals, respectively). Despite remarkable differences in breeding and movement ecology, Osprey populations of lower latitudes within the EUR-AFR did not show notable haplotypic differences when compared to the northern and central European birds. The four samples from the Canary Islands shared a single haplotype (Fig. 3). Within the IND-AUS group (n = 16) only three haplotypes were found, differing by only one nucleotide position. Five haplotypes were observed within the AMER group and 13 samples out of 17 showed the same haplotype, which was shared by ospreys from both western and eastern USA and from the Caribbean. Finally, within the ASIA group (n = 11) two haplotypes were recorded. The five samples from Japan were characterized by a single haplotype, which was shared with 2 samples from East Siberia.
Despite slight variation within each group (haplotypes were mainly distant by only one or two positions), a larger number of nucleotide differences were recorded between clades. The AMER group recorded the greatest genetic distance with EUR-AFR (a minimum of 15 nucleotides changed), whereas the minimum number of changes is 11 positions between IND-AUS and ASIA (Fig. 3).
Estimates of across and within-population variability of cyt b sequences of Osprey mtDNA
Molecular dating and biogeographical inference
Comparison of the six biogeographical reconstruction models for three different scenarios (“S0”, “S1”, “S2”; see text for definition)
DEC + J
DIVALIKE + J
BAYAREALIKE + J
Details of parameter estimations (dispersal (d), extinction (e), and founder-event speciation (j)) from the analyses of biogeographic inference of scenario S1 in BioGeoBEARS.The model with lowest AIC value is marked in bold font
Number of Parameters
DEC + J
DIVA + J
BayArea + J
Species diversity and demography
Our study revealed that the Osprey is structured into four main genetic groups, representing quasi non-overlapping geographical regions. Each lineage, though including birds from distant geographic areas, showed low internal genetic variability as revealed by the very low values (range: 0.1-0.2 %; Table 1) observed for the mean intra-group genetic distances. Haplotypic and nucleotide diversities were also weak, with only 9 haplotypes observed for the largest group including 102 ospreys (EUR-AFR), and values were even smaller for the other three lineages (Table 2). The two most common haplotypes in the clade EUR-AFR showed high overlap in ranges (e.g. individuals from very distant populations, like Finnish and Moroccan birds, were included in the same haplotype). A marked geographic substructuring was recorded only for individuals of the Canary Islands, Red Sea and Persian Gulf areas, each one representing a single haplotype (even if sample size is low; Fig. 3).
Despite the low variation within each group, the network (Fig. 3) revealed a high number of nucleotide differences between the four distinct clades. Overall, the mean sequence divergence across all populations (the entire Pandion mix) is 1.2 %, which is higher than the values recorded for the Red kite Milvus milvus (0.75 % observed for the mitochondrial DNA control region;  or the white-tailed eagle Haliaeetus albicilla (0.7 % for mitochondrial DNA control region; ). The mean genetic distance between groups (range: 1.5-2.6 %) is comparable to, or even greater than, those observed for the mitochondrial cyt b in several members of closely related sister eagle species from the genera Aquila, Hieraaetus (range: 1.7-2.1 %; ) and Haliaeetus (range: 0.3-9.8 %; ).
Populations within each group were poorly differentiated, suggesting that they might have experienced a reduction of genetic variation. Such low values are usually related to populations that experienced demographic crashes or remained isolated in fragmented habitats [38, 39]. However, the historical demographic reconstruction for each lineage suggested that they did not experience any strong bottleneck phases (Fig. 4), but rather underwent stable trends or slight increases.
For the EUR-AFR clade, Bayesian skyline curves (Fig. 4) suggest that populations encountered a recent expansion, which probably started about 10,000 years ago. This fits well with a recent review of the Holocene fossils of Osprey in central Europe from 10,000 years BP . In Europe, the Osprey expanded its distribution area in the following centuries until the beginning of the 20th century. Then, despite declines during the 1960-70s, populations were able to recover, showing positive demographic trends in most areas . On the other hand, the ASIA clade seems to be the only one that potentially suffered some demographic decline (Fig. 4). The same trend is also suggested by very low nucleotide diversity and a positive Fs value. However, such values might also be related to the small sample size of this group (15 samples). This result needs hence to be confirmed by further analyses with more extensive sampling from East Asia.
Why then does each clade show such low genetic variability? Low levels of genetic variation can be the consequence of recent population declines, or represent an ancestral condition, inherited from an ancient evolutionary history . To bring support to these two hypotheses, we compared haplotype diversity between museum and fresh specimens among the EUR-AFR group (n = 102). We found that for 10 oldest museum specimens (dated between 1872 and 1959; see Additional file 1), most haplotypes (7 out of 9) are also represented in the fresh samples, thus suggesting that there has been no recent loss of haplotypic diversity. Several studies have also reported stable genetic diversity despite declines in population size (e.g. ). For example, no obvious loss of genetic diversity was detected among Canadian peregrine falcons (Falco peregrinus) despite a population bottleneck in the 20th century . In our case, one possible explanation is the colonization of new areas by few individuals carrying only a few mitochondrial lineages from the genetic pool (founder effect). Furthermore, source populations could have experienced a reduction in genetic variability due to climate fluctuations during the Quaternary; remnant populations in refugia represented the genetic source for the following recolonization (see below).
The cyt b phylogenetic tree (Fig. 2a and Additional file 4a) suggests that the AMER haplogroup constitutes the first genetic group that emerged among the four groups evidenced for Pandion haliaetus. This relationship is more supported by the cyt b alone (1103 bp) than by the ND2 alone (1078 bp; Additional file 4b). Moreover, the combination of cyt b and ND2 (2037 bp, Fig. 2b) did not increase markedly the resolution, thus suggesting that the divergence of the four major clades occurred over a relatively short time period. Another explanation could be that irresolution arose as a misplacement of the root of the tree due to the use of too divergent outgroups. Indeed, there are no optimal outgroups available for the phylogeny rooting of Pandion, since it constitutes a long branch in the phylogeny of raptors, distant from its sister Accipitridae family [11, 45]. In preliminary analyses we also included other outgroups such as Sagittarius serpentarius and/or different species from the Cathartidae family, but no evident improvements in branching were recorded. Similar problems of rooting have been reported in other phylogeography studies [46, 47].
The major result yielded by phylogeographic inferences is a split of Pandion in four clades that appear to have diverged from each other within relatively short time. Because of the uncertainty with regard to the relative branching order of these clades (see Results), patterns of colonization phases should be considered as tentative.
The best biogeographical scenario supports America as the most likely area of origin of Pandion (step 1 Fig. 7). The selected model (DIVAj) favours an expansion of Osprey populations that underwent diversification by subsequent founder events followed by vicariance processes. Synthetically, from America, Osprey should have colonized the Indo-Australasian region passing through the Pacific coast of Asia (step 2 Fig. 7). From refugia located in Indonesia-Oceania, rapid range expansion allowed the settlement in eastern Asia and the Western Palearctic (Fig. 5, Fig. 6, step 3 Fig. 7). As a supporting case, cyt b data suggest that Japan was likely colonized after Russia and not directly from Australia. Such a scenario would explain the current distribution of Pandion across all continents and the genetic divergence found among the four lineages. This hypothesis is corroborated by the p-distances which displayed lower values between AMER and IND-AUS (p = 0.020) than those recorded between AMER and EUR-AFR (0.026). This result indicates that the EUR-AFR group did not originate directly from the AMER group, as also suggested by the lowest likelihood values yielded by the alternative biogeographical reconstruction (scenario “S2” in which EUR-AFR would have been colonized from AMER; Table 3).
Implications for taxonomy and conservation
The four genetic groups that we found do not entirely correspond to the four subspecies described, based on morphological characters [4, 13]. The IND-AUS lineage fully matched geographically with the subspecies cristatus . Contrary to current taxonomy, North American carolinensis ospreys did not differ from Caribbean ridgwayi birds as there is no evident structure in the phylogeographic tree and in the network. On the other hand, in Eurasia, we discovered that the subspecies haliaetus was actually composed of 2 lineages (EUR-AFR and ASIA) that are undistinguishable morphologically. This reflects the poor knowledge of the species in Asia where detailed information about biology and distribution are needed . More samples should be collected to clarify the geographical limits of these lineages, especially in the regions where we found a zone of overlap (e.g. central Siberia with EUR-AFR clade and in Indonesia with IND-AUS clade).
Overall, genetic distances between Osprey clades (Table 1) are in a range which has already been used by taxonomists for designating distinct raptor species (e.g. ). However, we detected a relatively low overall nucleotide diversity (1.0 %) compared to other large raptor species with a similar wide distribution (e.g. Gypaetus barbatus, 2.9 % for the control region; ; but see  for Haliaetus albicilla, 0.7 %). It might be objected that more variability would potentially have been found by including more geographic samples for the different clades. However, if we consider the EUR-AFR clade (131 individuals covering a large geographic area), the variability was as weak as for other groups that are represented by a much smaller number of birds (15 to 37). The decision for splitting ospreys into different species (e.g. the split of P. h. cristatus as a full species called Eastern Osprey; ) should integrate also other factors besides morphology and mtDNA differences; e.g. nuclear genes and behavioural aspects, including migrations that could play an important role as reproductive barriers between distant populations .
A first step towards a sound global management and conservation plan is to define Management Units (MUs) and Evolutionary Significant Units (ESUs; ). As a matter of fact, subspecies have often been used as proxies for units of conservation in absence of a genetic data indicating distinct evolutionary units . Our results evidenced four different lineages that may deserve the status of ESUs for specific management actions. A better knowledge of the distribution range of each lineage is strongly needed in the near future. In particular, the ASIA lineage should be a priority target for multiple reasons: 1) this lineage has never been described before; 2) it relied on a limited number of samples (n = 15) from only a few areas; 3) the majority of the samples analysed were museums specimens, so the current presence of this lineage in East Russia and Indonesia needs to be confirmed; 4) it is the only one clade showing signals of demographic decline; and 5) no clear morphological characteristics can presently help identification.
Our study revealed that, at the global scale, the particular evolutionary history of the Osprey has partitioned the species into 4 distinct clades with clear genetic differentiation. In addition, our findings indicate that broad habitat requirements and high mobility of Osprey were more important factors than philopatry in shaping genetic diversity at the intra-clade level. Further genetic study using microsatellite markers are ongoing (Monti et al., in prep.) in order to reveal more recent differences and to quantify gene flow among populations that show differences in migratory or reproductive behaviours (e.g. North and Central European populations vs Mediterranean, Canary Islands, Cape Verde and Red Sea).
Even though the Osprey is currently globally listed as “Least Concern” according to IUCN criterions , it is considered a priority species for conservation across its distributional range. Indeed, the Osprey is taxonomically unique and conserving its phylogenetic diversity should be a priority [53–55]. The species experienced a severe decline during the 19th and 20th century that led to significant demographic declines and local extinctions [56–58]. Consequently, the Osprey has emerged as an important flagship species and during recent decades has been involved in 25 reintroduction projects in the USA [9, 58] and in Europe [57, 59]. Our results suggest that future reintroduction projects should be conducted using source populations from the same lineage (e.g. within Europe birds chosen for translocation should originate from the Western Palearctic, avoiding individuals belonging to the other lineages). However, before concluding that no restriction needs to be adopted for translocations between populations within the Western Palearctic, other ecological, demographical or behavioural variables should be considered.
A specific ethical committee was not required for this study. Bird handling was performed under animal experimentation permits 34–369 (David Grémillet) delivered by the French ‘Direction Départementale de la Protection des Populations’ and under the licence of Olivier Duriez from the Centre for Bird Population Studies (CRBPO) of the Natural History Museum (MNHN, Paris): according to the French law of 22 September 2008, the CRBPO has the delegation by the Ministry of Ecology, Energy, Sustainable Development and Land Settlement for allowing the owners of a general bird ringing licence to capture and handle birds from protected species, and collect samples or mark them (with rings or any other device like GPS units).
Availability of supporting data
The data set supporting the results of this article is included within the article and its additional files “Additional files 1, 2, 3 and 4.docx”. Nucleotide sequences have been submitted to GenBank and accession numbers are provided in Additional file 1. Alignments used for each analysis (please refer to the table in Additional file 2 for matching) have been included as Additional file 5.
We thank all the colleagues who kindly provided Osprey samples: Pertti Saurola, Harri Koskinen and Juhani Koivu (Finland), Urmas Sellis (Estonia), Aigars Kalvāns (Latvia), Miroslav V. Babushkin and Vladimir V. Ivanovsky (Russia), Martin Haase and Daniel Schmidt (Germany), Rolf Wahl, Franҫois Gossmann and Eric Pasquet (France), Luis Palma and Ursula Hofle Hansen (Portugal), Manuel Siverio (Canary islands), Pedro López-Suárez (Cape Verde), Houssine Nibani (Morocco), Rafel Triay Bagur (Balearic islands). Special thanks to Yula Kapetanakos and Alan Poole (Cornell lab of Ornithology, Ithaca, USA) for collecting and performing free analyses of 10 samples from the USA. Many thanks are also addressed to the following museums and birds’ collection curators who allowed us to gather samples from their collection: Mark Adams (Natural History Museum of London; UK), Steven van der Mije and Bechy Desjardins (Netherlands Centre for Biodiversity Naturalis of Leiden), Isao Nishiumi (National Museum of Nature and Science of Amakubo; Japan), Mikhail Kalyakin, Pavel Tomkovitch, Yaroslav Red’kin (Zoological Museum of Moscow; Russia), Guillermo Delgado (Museo de la Naturaleza y el Hombre of Tenerife Island; Spain) and Enrico Borgo (Museo Civico di Storia Naturale Giacomo Doria di Genova; Italy). We also thank the Portuguese Navy and in particular the contra-almirante José Luís Branco Seabra de Melo (Head of Cabinet of the Navy Chief of Staff) and the capitão-de-mar-e-guerra José Jaime Gonçalves Ribeiro (Director of the Vasco da Gama Aquarium) for providing historical samples of Portuguese ospreys. All extractions from museum specimens were performed using the facilities of the platform “ADN dégradé” (Labex CeMEB, Montpellier, France) coordinated by Christelle Tougard. Sequencing was supported by the project SEQMED attributed to CM by the “Bibliothèque du Vivant” sequencing project (CNRS-INRA-MNHN). Alan Poole provided valuable comments and improved the English of earlier drafts. This study was funded thanks to a grant from the Foundation Prince Albert II de Monaco and the Associazone italiana delle Fondation Prince Albert II de Monaco ONLUS. It also benefits from the support of the Parc Naturel Régional de Corse, Maremma Regional Park and Conservatoire du Littoral (through the Mediterranean Small Islands Initiative PIM). F. Monti is supported by a grant of the Maremma Regional Park (Italy) and by a mobility grant of the Università Italo Francese/Université Franco Italienne (UIF/UFI).
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.
- Hewitt G. The genetic legacy of the Quaternary ice ages. Nature. 2000;405:907–13.View ArticlePubMedGoogle Scholar
- Hewitt GM. Genetic consequences of climatic oscillations in the Quaternary. Philos Trans R Soc B Biol Sci. 2004;359:183–95.View ArticleGoogle Scholar
- Newton I. Speciation and Biogeography of Birds. London & San Diego: Academic; 2003.Google Scholar
- Poole AF. Ospreys: a Natural and Unnatural History. New York: Cambridge University Press; 1989.Google Scholar
- Alerstam T, Hake M, Kjellén N. Temporal and spatial patterns of repeated migratory journeys by ospreys. Anim Behav. 2006;71:555–66.View ArticleGoogle Scholar
- Thibault JC, Triay R, Beabrun P, Boukhalfa D, Dominici JM, Torre A. Osprey (Pandion haliaetus) in the Mediterranean: characteristics of a resident population with a patchy distribution. In: Muntaner J, Mayol J, editors. Biologia y Conservacion de las Rapaces Mediterraneas. Madrid: Monografias, n°4. SEO; 1996. p. 135–44.Google Scholar
- Hailer F, Helander B, Folkestad AO, Ganusevich SA, Garstad S, Hauff P, et al. Phylogeography of the white-tailed eagle, a generalist with large dispersal capacity. J Biogeogr. 2007;34:1193–206.View ArticleGoogle Scholar
- Bell DA, Griffiths CS, Caballero IC, Hartley RR, Lawson RH. Genetic Evidence for Global Dispersal in the Peregrine Falcon (Falco peregrinus) and Affinity with the Taita Falcon (Falco fasciinucha). Journal of Raptor Research. 2014;48:44–53.View ArticleGoogle Scholar
- Martell MS, Englund JV, Tordoff HB. An urban Osprey population established by translocation. Journal of Raptor Research. 2002;36(2):91–6.Google Scholar
- Wink M, Sauer-Gürth H. Phylogenetic relationships in diurnal raptors based on nucleotide sequences of mitochondrial and nuclear marker genes. In: Raptors Worldwide. Edited by Chancelor RD, Meyburg BU. WWGBP, Berlin; 2004.Google Scholar
- Lerner HRL, Mindell DP. Phylogeny of eagles, Old World vultures, and other Accipitridae based on nuclear and mitochondrial DNA. Mol Phylogenet Evol. 2005;37:327–46.View ArticlePubMedGoogle Scholar
- Griffiths CS, Barrowclough GF, Groth JG, Mertz LA. Phylogeny, diversity, and classification of the Accipitridae based on DNA sequences of the RAG-1 exon. J Avian Biol. 2007;38:587–602.View ArticleGoogle Scholar
- Strandberg R. Ageing, sexing and subspecific identification of Osprey, and two WP records of American Osprey. Dutch Birding. 2013;35:69–87.Google Scholar
- Wink M, Sauer-Gürth H, Witt HH. Phylogenetic differentiation of the Osprey (Pandion haliaetus) inferred from nucleotide sequences of the mitochondrial cytochrome b gene. In: Raptors Worldwide. Edited by Chancelor RD, Meyburg BU. WWGBP, Berlin; 2004.Google Scholar
- Helbig AJ, Schmidt D, Seibold I. Mitochondrial DNA sequences reveal differentiation between Nearctic and Palearctic Osprey (Pandion haliaetus) populations [abstract]. Biol Cons Fauna. 1998;102:224.Google Scholar
- Moritz C. Defining “Evolutionarily Significant Units” for conservation. Trends Ecol Evol. 1994;9:373–5.View ArticlePubMedGoogle Scholar
- Mundy NI, Unitt P, Woodruff DS. Skin from feet of museum specimens as a Non-destructive source of DNA for avian genotyping. The Auk. 1997;114:126–9.View ArticleGoogle Scholar
- Bantock TM, Prys-Jones RP, Lee PLM. New and improved molecular sexing methods for museum bird specimens. Mol Ecol Resour. 2008;8:519–28.View ArticlePubMedGoogle Scholar
- Galtier N, Gouy M, Gautier C. SEAVIEW and PHYLO_WIN: two graphic tools for sequence alignment and molecular phylogeny. Comput Appl Biosci. 1996;12:543–8.PubMedGoogle Scholar
- Lanfear R, Calcott B, Ho SYW, Guindon S. PartitionFinder: combined selection of partitioning schemes and substitution models for phylogenetic analyses. Mol Biol Evol. 2012;29:1695–701.View ArticlePubMedGoogle Scholar
- Ronquist F, Huelsenbeck J, MrBayes 3. Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003;19:1572–4.View ArticlePubMedGoogle Scholar
- Rambaut A, Drummond AJ. Tracer: MCMC Trace Analysis Tool. Edinburgh: Institute of Evolutionary Biology, University of Edinburgh; 2007.Google Scholar
- Stamatakis A. RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics. 2006;22:2688–90.View ArticlePubMedGoogle Scholar
- Felsenstein J. PHYLIP (Phylogeny Inference Package) version 3.6. Distributed by the author. Seattle: Department of Genome Sciences, University of Washington; 2005.Google 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
- Rozas J, Sánchez-DelBarrio JC, Messeguer X, Rozas R. DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics. 2003;19:2496–7.View ArticlePubMedGoogle Scholar
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28:2731–9.PubMed CentralView ArticlePubMedGoogle Scholar
- Ramos-Onsins SE, Rozas J. Statistical properties of new neutrality tests against population growth. Mol Biol Evol. 2002;19:2092–100.View ArticlePubMedGoogle Scholar
- Fu Y-X. Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics. 1997;147:915–25.PubMed CentralPubMedGoogle Scholar
- Drummond AJ, Rambaut A, Shapiro B, Pybus OG. Bayesian coalescent inference of past population dynamics from molecular sequences. Mol Biol Evol. 2005;22:1185–92.View ArticlePubMedGoogle Scholar
- Drummond AJ, Rambaut A. BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007;7:214.PubMed CentralView ArticlePubMedGoogle Scholar
- Schmidt HA, Strimmer K, Vingron M, von Haeseler A. TREE-PUZZLE: maximum likelihood phylogenetic analysis using quartets and parallel computing. Bioinformatics. 2002;18:502–4.View ArticlePubMedGoogle Scholar
- Nabholz B, Glémin S, Galtier N. The erratic mitochondrial clock: variations of mutation rate, not population size, affect mtDNA diversity across birds and mammals. BMC Evol Biol. 2009;9:54.PubMed CentralView ArticlePubMedGoogle Scholar
- Matzke NJ. BioGeoBEARS: BioGeography with Bayesian (and Likelihood) Evolutionary Analysis in R Scripts. R package, version 0.2.1, published July 27, 2013 at https://cran.r-project.org/web/packages/BioGeoBEARS/index.html.
- Matzke NJ. Thesis abstract: probabilistic historical biogeography: new models for founder-event speciation, imperfect detection, and fossils allow improved accuracy and model-testing. Front Biogeogr. 2013;5(4):242–8.Google Scholar
- Roques S, Negro JJ. MtDNA genetic diversity and population history of a dwindling raptorial bird, the Red Kite (Milvus milvus). Biol Conserv. 2005;126:41–50.View ArticleGoogle Scholar
- Wink M, Heidrich P, Fentzlo C. A mtDNA phylogeny of sea eagles (genus Haliaeetus) based on nucleotide sequences of the cytochrome b gene. Biochem Syst Ecol. 1996;24:783–91.View ArticleGoogle Scholar
- Lifjeld J, Bjornstad G, Steen O, Nesje M. Reduced genetic variation in Norwegian Peregrine Falcons Falco peregrinus indicated by minisatellite DNA fingerprinting. Heredity. 2002;144:19–26.Google Scholar
- Martínez-Cruz B, Godoy JA, Negro JJ. Population genetics after fragmentation: the case of the endangered Spanish imperial eagle (Aquila adalberti). Mol Ecol. 2004;13:2243–55.View ArticlePubMedGoogle Scholar
- Zachos FE, Schmölcke U. Archeozoological records and distribution history of the osprey (Pandion haliaetus) in central Europe. Journal of Ornithology. 2006;147:565–8.View ArticleGoogle Scholar
- Schmidt-Rothmund D, Dennis R, Saurola P. The osprey in the western palearctic: breeding population size and trends in the early 21 st century. Journal of Raptor Research. 2014;48:375–86.View ArticleGoogle Scholar
- Wandeler P, Hoeck PEA, Keller LF. Back to the future: museum specimens in population genetics. Trends Ecol Evol. 2007;22:634–42.View ArticlePubMedGoogle Scholar
- Pertoldi C, Hansen MM, Loeschcke V, Madsen AB, Jacobsen L, Baagoe H. Genetic consequences of population decline in European Otter Lutra lutra: An assessment of microsatellite DNA variation in Danish otters from 1883 to 1993. Proc R Soc Lond B. 2001;268:1775–81.View ArticleGoogle Scholar
- Brown JW, de Groot PJV, Birt TP, Seutin G, Boag PT, Friesen VL. Appraisal of the consequences of the DDT-induced bottleneck on the level of geographic distribution of neutral genetic variation in Canadian peregrine falcons, Falco peregrinus. Mol Ecol. 2007;16:327–43.View ArticlePubMedGoogle Scholar
- Helbig AJ, Kocum A, Seibold I, Braun M. A multi-gene phylogeny of aquiline eagles (Aves: Accipitriformes) reveals extensive paraphyly at the genus level. Mol Phylogenet Evol. 2005;35:147–64.View ArticlePubMedGoogle Scholar
- Godoy JA, Negro JJ, Hiraldo F, Donázar JA. Phylogeography, genetic structure and diversity in the endangered bearded vulture (Gypaetus barbatus, L.) as revealed by mitochondrial DNA. Mol Ecol. 2004;13:371–90.View ArticlePubMedGoogle Scholar
- Igea J, Aymerich P, Fernández-González A, González-Esteban J, Gómez A, Alonso R, et al. Phylogeography and postglacial expansion of the endangered semi-aquatic mammal Galemys pyrenaicus. BMC Evol Biol. 2013;13:115.PubMed CentralView ArticlePubMedGoogle Scholar
- Christidis L, Boles WE. Systematics and taxonomy of Australian birds. Collingwood, Australia: CSIRO Publishing; 2008.Google Scholar
- Shoji A, Sugiyma A, Brazil MA. The status and breeding biology of Ospreys in Hokkaido, Japan. Condor. 2011;113:762–7.View ArticleGoogle Scholar
- Helbig AJ, Knox AG, Parkin DT, Sangester G, Collinson M. Guidelines for assigning species rank. Ibis. 2002;144:518–25.View ArticleGoogle Scholar
- Zink RM. The role of subspecies in obscuring avian biological diversity and misleading conservation policy. Proc R Soc Lond B. 2004;271:561–4.View ArticleGoogle Scholar
- BirdLife International (2015). Pandion haliaetus. The IUCN Red List of Threatened Species. Version 2015-3. <www.iucnredlist.org>. Downloaded on 12 November 2015.
- Winter M, Devictor V, Schweiger O. Phylogenetic diversity and nature conservation: where are we? Trends Ecol Evol. 2013;28:199–204.View ArticlePubMedGoogle Scholar
- Jetz W, Thomas GH, Joy JB, Redding DW, Hartmann K, Mooers AO. Global distribution and conservation of evolutionary distinctness in birds. Curr Biol. 2014;24:919–30.View ArticlePubMedGoogle Scholar
- Nunes LA, Turvey ST, Rosindell J. The price of conserving avian phylogenetic diversity: a global prioritization approach. Philos Trans R Soc Lond B Biol Sci. 2015;370:20140004.PubMed CentralView ArticlePubMedGoogle Scholar
- Cramp S, Simmons KEL. Handbook of the birds of Europe, vol. 2. Oxford: Oxford University Press; 1980. p. 265–77.Google Scholar
- Dennis R, Dixon H. The experimental reintroduction of Ospreys Pandion haliaetus from Scotland to England. Vogelwelt. 2001;122:147–54.Google Scholar
- Bierregaard RO, David AB, Gibson L, Kennedy RS, Poole AF, Scheibel MS, et al. Post-DDT recovery of osprey (Pandion haliaetus) populations in Southern New England and Long Island, New York, 1970–2013. Journal of Raptor Research. 2014;48:361–74.View ArticleGoogle Scholar
- Monti F, Sforzi A, Dominici JM. Post-fledging dependence period of ospreys Pandion haliaetus released in central Italy: home ranges, space use and aggregation. Ardeola. 2012;59:17–30.View ArticleGoogle Scholar