- Research article
- Open Access
Genomic evidence of demographic fluctuations and lack of genetic structure across flyways in a long distance migrant, the European turtle dove
BMC Evolutionary Biology volume 16, Article number: 237 (2016)
Understanding how past climatic oscillations have affected organismic evolution will help predict the impact that current climate change has on living organisms. The European turtle dove, Streptopelia turtur, is a warm-temperature adapted species and a long distance migrant that uses multiple flyways to move between Europe and Africa. Despite being abundant, it is categorized as vulnerable because of a long-term demographic decline. We studied the demographic history and population genetic structure of the European turtle dove using genomic data and mitochondrial DNA sequences from individuals sampled across Europe, and performing paleoclimatic niche modelling simulations.
Overall our data suggest that this species is panmictic across Europe, and is not genetically structured across flyways. We found the genetic signatures of demographic fluctuations, inferring an effective population size (Ne) expansion that occurred between the late Pleistocene and early Holocene, followed by a decrease in the Ne that started between the mid Holocene and the present. Our niche modelling analyses suggest that the variations in the Ne are coincident with recent changes in the availability of suitable habitat.
We argue that the European turtle dove is prone to undergo demographic fluctuations, a trait that makes it sensitive to anthropogenic impacts, especially when its numbers are decreasing. Also, considering the lack of genetic structure, we suggest all populations across Europe are equally relevant for conservation.
Climate oscillations during the Quaternary strongly affected organismic evolution worldwide (e.g. [1–3]). The responses of different species to these oscillations have been diverse, and dependent on their environmental requirements and adaptations [4, 5]. In Europe, for example, most species that are adapted to warm-temperate climates underwent southward range shifts and reductions in their abundance during glacial periods, recovering their northern distributions and numbers during interglacial periods [6–8]. Similarly, current global warming and human activities are also affecting the evolution and distribution of living organisms worldwide (e.g. [9, 10]). Long distance migratory bird species are particularly sensitive to global climate change; their migration routes and timing have shifted drastically during the past decades [11, 12], mainly as a consequence of phenological responses to the rising temperatures . Moreover, long distance migrants have suffered stronger population declines than other bird species, because their survival depends on the quality of habitats that are often geographically distant (i.e. breeding, stopover and wintering territories) [11, 14, 15].
Around 2 billion birds migrate between Europe and Africa each year , crossing the Mediterranean Sea and the Sahara Desert. Even though migration of some species occurs broadly throughout the Mediterranean Sea, higher numbers of birds migrate over its eastern and western edges , thus reducing over-sea flying distances. As a consequence, alternative flyways have arisen in some species (e.g., White storks, Ciconia ciconia, ; Montagu’s harrier, Circus pygargus, . Similar strategies have been observed in other long distance migratory birds across the world (e.g. North American’s Swainson’s thrush, Catharus ustulatus, ; the Asian greenish warbler, Phylloscopus trochiloides, ). In some cases, such alternative flyways might be mirrored in the breeding grounds in the form of migratory divides, defined as narrow regions of contact between populations with different migratory pathways . The presence of alternative flyways and migratory divides has been postulated as a mechanism promoting intraspecific genetic diversification .
The European turtle dove, Streptopelia turtur (hereafter turtle dove), is the only columbid that takes part in the African-European migration system . The species is warm-climate adapted and requires patchy areas composed by shrubby woodlands and open grasslands to breed . The turtle dove breeds across the western Palearctic and spends winters in the sub-Saharan African region, in latitudes between 10°N and 20°N . Mark-recovery data suggests that turtle doves use three main flyways during their spring and autumn journeys . Birds that breed in western Europe fly over the Iberian Peninsula, those breeding in central Europe fly through Italy and Malta, and eastern European breeders fly over Greece (Fig. 1a). Despite being an abundant species, the turtle dove has been undergoing a long-term demographic decline, and as a result is considered a vulnerable species throughout its European distribution . Both the role of the alternative migratory flyways on genetically structuring turtle dove populations, and how the long-term demographic decline has impacted the species’ genetic diversity remain unknown.
Here we studied the demographic history and population genetic structure of the turtle dove across Europe. We performed a genome wide scan using double digest restriction-site associated DNA sequencing (ddRADseq) and obtained mitochondrial DNA sequences (cytochrome b). Demographic history analyses were complemented with niche modelling analysis to test if variations in the effective population size are correlated with past changes in the availability of suitable habitat. Based on the ecological requirements of our focal species, we expected to find evidence of large effective population sizes during warm interglacial periods, followed by small effective sizes during cold glacial periods. We would expect to observe a similar trend in the availability of suitable habitat. We also hypothesized that the current demographic decline in the turtle dove would be strong enough to show the genetic signature of a decrease in its effective population size. Finally, we tested whether the alternative flyways described for the turtle dove actually correspond to migratory divides promoting intraspecific genetic diversification. In that case we expect to find three genetically differentiated groups associated to each flyway (i.e. East, Central and West).
We analysed samples from 104 adult individuals, representing the three main migratory flyways of the turtle dove (Fig. 1a). The western flyway (WF) was represented with samples from the United Kingdom (n = 18, breeding season), Spain (n = 10, breeding/migrating), Germany (n = 1, breeding season) and France (n = 19, breeding); the Central flyway (CF) with samples from Italy (n = 15, migrating) and Malta (n = 11, migrating); and the Eastern flyway (EF) with samples from Bulgaria (n = 14, breeding/migrating) and Greece (n = 16, breeding/migrating). Sampling took place between 2011 and 2014. Tissue samples were taken from wings of hunted individuals, and blood samples were taken from the brachial vein of individuals either caught at their breeding sites or during migration while on a flyway. Tissue samples were preserved in ethanol, and blood samples were preserved either in ethanol or on FTA cards (Whatman International Ltd., Maidstone, UK). A table showing more details of the analysed samples can be found in the additional files section (Additional file 1).
Genomic DNA was isolated using the DNeasy blood and tissue kit (Qiagen, CA) and quantified with a Qubit fluorometer (Invitrogen, CA). Double-digest restriction-site associated DNA markers (ddRADtags) were generated following Peterson et al., . We also obtained 95 sequences of the cytochrome b mitochondrial gene, 892 base pairs long (GenBank accession numbers: KU588290-KU588384). Primers used were specifically designed for this study: F_cytb_St: 5′-TGATAACTCAAATCCTAACTGGTC-3′ and R_cytb_St: 5′-TTGTTTTCTAGGGCTCCGAT-3′. Further details on the ddRADtags and cytochrome b amplification laboratory methods are provided in the additional files section (Additional file 2).
RAD-sequence pre-processing and single nucleotide polymorphism (SNP) calling
Filtering and quality control
The samples were sequenced on two Illumina HiSeq 2500 lanes that produced a total of ~304 million, 151 bp single-end reads. After assessing quality using FASTQC (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/), we trimmed the final 29 bp to discard lower-quality base calls using FASTX-Toolkit (http://hannonlab.cshl.edu/fastx_toolkit/). We filtered reads with at least a single base call with a Phred quality score below 10 (90 % call accuracy), and/or more than 5 % below a Phred score of 20 (99 % call accuracy). We demultiplexed reads with the process_radtags program from the STACKS v1.20  bioinformatics pipeline. We discarded reads that did not contain an SbfI cut site or a unique P1 barcode in the 5′ end, as well as those with adapter contamination in the 3′ end. Finally, because the barcodes used for demultiplexing ranged between 5 and 7 bp, we trimmed all reads to 115 bp, the length of the shortest sequences once barcodes were removed.
Reference-based RAD loci assembly and variant calling
We obtained 1.03 ± 0.46 million, 115 bp quality-filtered sequences per individual, which were aligned to the Columba livia genome (assembly 1.0) with BOWTIE2 v2.2.2 . The average alignment rate was 85.7 %. We conducted a referenced-based assembly using the pstacks/cstacks/sstacks pipeline controlled by the ref_map script in STACKS, followed by the rxstacks error correction module and a final iteration of cstacks/sstacks. This pipeline will identify a locus when a minimum coverage is reached: set to 20 with the m parameter. A catalogue of all loci is generated, and we allowed up to five differences among aligned loci of different individuals (parameter n). The error correction module examines SNP and haplotype calls across all individuals and removes loci that match more than one locus in the catalogue (presumably from repetitive areas of the genome), eliminates excess haplotypes that are uncommon and thus probably erroneous, and eliminates loci below a user established likelihood threshold (set to −200 based on the distribution across all loci). Once errors are corrected a new catalogue is generated and SNPs are identified. The pipeline produced a catalogue consisting of 18,026 loci. We exported SNP and sequence data for downstream analysis using the program populations, retaining loci that were present in every sampling locality, had less than 20 % missing data and a coverage of at least 20x. These settings produced 3,012 RAD loci, 2,934 of which had at least one SNP. When a RAD locus had more than one SNP, we selected the first one for downstream analyses to avoid using tightly linked SNPs.
Analysis of population genetic structure
Summary statistics of levels of polymorphism (i.e. haplotype and nucleotide diversity) were obtained with DNASP v5.10.1 . We obtained a median joining haplotype network, using the software POPART . Hierarchical analyses of molecular variance (AMOVAs) were performed in ARLEQUIN v220.127.116.11 , with significance values obtained after 10,000 permutations. We tested if grouping individuals by migration flyway explained significant levels of molecular variance, a signal of genetic differentiation. We tested for the possibility of three genetically differentiated groups (i.e. flyways), a Western group including samples from United Kingdom, France, Germany and Spain; a Central group containing the samples from Italy and Malta, and an Eastern group with the samples from Bulgaria and Greece.
Genetic structure was investigated using two kinds of approaches; model and non-model based methods. For a model-based approach we employed the variational Bayesian inference method implemented in fastSTRUCTURE , to investigate the most likely number of genetic clusters (K) that better explained our data. Variational Bayesian inference aims to solve the problem of inference as an optimization rather than a sampling problem, a key advantage of this method is that it runs faster while achieving similar levels of accuracy than other available algorithms . We performed a first exploratory run using the simple prior testing from K = 1 to K = 10. Afterwards, three independent runs with different starting random seeds were performed with the logistic prior, testing from K = 1 to K = 3; the usage of the latter prior is advised when genetic structure is difficult to determine or weak . During both runs all other settings were left as default. Mean admixture proportions (Q-values) for the value of K that better explained our data were graphed with DISTRUCT v1.1 . For non-model based methods we performed a Principal Component Analysis using the function dudi.pca, from the ADE4 R package . An AMOVA was conducted with ARLEQUIN v18.104.22.168  following the grouping scheme described above (see mtDNA analyses). We also searched for loci under selection among the groups described above using BAYESCANR  included in the R package PopGenome .
We searched for signatures of past population expansion or contraction in the turtle dove effective population size (Ne). The Ne is defined as the population size that is evolutionary relevant (i.e., the number of breeding individuals) . We used DNASP v5  to investigate our mtDNA data set, calculating Fu’s F s and Tajima’s D; the significance for each test was estimated after 1,000 permutations.
We also explored the demographic history of the turtle dove by performing an Extended Bayesian Skyline plot (EBSP), as implemented in BEAST version 1.8.2 . We did not find evidence of population genetic structure (Fig. 1b and c), and thus were able to pool all our samples into one population. Because of the computationally intensive nature of this analysis, we conducted runs including only 30 loci per individual. We used our cytochrome b data and sequences from 29 randomly selected nuclear loci. All partitions were unlinked and we implemented a strict clock model, with a substitutions rate of 2.1 % per million years as estimated by Weir and Schluter  for the cytochrome b gene in birds. We ran the analysis using the “Coalescent: Extended Bayesian Skyline” tree prior for 100 million generations, sampling every 4,000. We inspected the results in Tracer version 1.6 , checking for convergence, sufficiently large effective sample sizes for parameter estimates and discarding the initial 10 % of samples as burn-in. This analysis was run three times using independent sets of 29 randomly selected nuclear loci in combination with the cytochrome b sequences; all three runs produced similar results.
An Approximate Bayesian Computation (ABC) analysis was performed using DIYABC v2 . ABC methods bypass exact likelihood calculation by using summary statistics and massive simulations, allowing them to handle large data sets . Thus, we included the complete SNP dataset (2954 loci) for all genotyped individuals (104) and pooled our samples into a single population (as previously explained for the EBSP). We compared five alternative demographic scenarios, with different combinations of Ne variation in time (t), see Fig. 2 for a schematic representation of the scenarios. All priors were set as normally distributed, time (t) expressed in million years before present and effective population size (N) in million individuals (Fig. 2). Priors N1, N2, t1 and t2 were informed by the EBSP (see results). Prior N4 at t3 was set equal to N1 based on a similar availability of suitable habitat estimated for the last interglacial (LIG) and the present (see results). Prior N3 was based on the demographic decline reported for the turtle dove . Simulated scenarios were generated using all the available summary statistics in DIYABC for a single population analysis based on SNPs: proportion of monomorphic loci, mean gene diversity across polymorphic loci, variance of gene diversity across polymorphic loci and mean gene diversity across all loci. Runs consisted of 5 ×106 simulations (1 million simulations per tested scenario). The PCA function embedded in DIYABC was used to test the goodness-of-fit of the observed data against that generated by the simulated scenarios. The relative posterior probability was then estimated for each scenario with the direct and logistic regression approaches, on 0.1 and 1 % respectively, of the simulations closest to the observed data set. Finally, we estimated the posterior distribution of parameters for the most likely scenario; we assumed a generation time of one year for the turtle dove  to scale the estimated parameters.
The realized niche of the present distribution of the turtle dove was modelled and extrapolated using MaxEnt 3.3.3 k (maximum entropy; ). We combined data obtained from GBIF (Global Biodiversity Information Facility, http://www.gbif.org/), with the Bioclim data sets (http://worldclim.org/) for the present, the Mid-Holocene (MH, 6000 ybp), the Last Glacial Maximum (LGM, 22,000 ybp), and the last interglacial (LIG; ~120,000–140,000 ybp). MaxEnt t is a program for modelling ecological niches from presence-only species records (further details in Additional file 2).
Summary statistics and population genetic structure
Cytochrome b sequences were highly polymorphic, with 73 segregating sites and 40 haplotypes found among 95 turtle dove individuals (haplotype diversity (Hd): 0.92; nucleotide diversity (π): 0.0075). Our median joining network suggested two main haplogroups, differentiated by six mutational steps (Fig. 1b). However, the observed haplogroups did not cluster samples by migration flyway or their geographic origin. Moreover, haplotypes from individuals sampled in the same locality appeared scattered across the network (Fig. 1b). Similarly, the AMOVA based on mtDNA sequences did not show evidence that migration flyways explain population genetic structure in the turtle doves. The highest percentage of the variance was explained by the components “among individuals, within population” and “within individuals” (Additional file 3).
SNP-based analyses with fastSTRUCTURE -using the simple prior- suggested that the number of components used to explain the genetic structure in our data (K * Øc) and the model complexity that maximizes the marginal likelihood (K * ε) were both 1. When the logistic prior was used, K * Øc was still 1, however K * ε = 3. Even though K = 3 was suggested with the latter prior, the greatest proportion of the ancestry of all individuals was explained by only one cluster (mean Q = 0.999) (shown in Additional file 4). Hence, K = 1 seems to be the most likely partition to explain genetic structure in our data (Additional file 4). The AMOVA result based on SNPs resembled the results obtained for the mtDNA sequences, with the “among groups” component explaining the lowest percentage of the variance when testing for the differences among flyways (Additional file 3). The PCA based on SNP data did not show any clustering of individuals (Fig. 1c). Finally, we did not find loci under selection in BayescanR (p > 0.05).
Demographic history and niche modelling
The mtDNA dataset showed signatures of a recent demographic expansion in the turtle dove population: F s = −20.95 (p = 0.001) and Tajima’s D = −1.70 (p = 0.02). The EBSP analysis that combined mtDNA sequences with 29 randomly selected nuDNA loci also supported a population expansion. Based on our assumption of molecular clock calibration for cytochrome b (2.1 per million years), this expansion started ~14,000 years before present (Fig. 3).
Our DIYABC analysis supported a compatible, yet more complex demographic scenario than that suggested by the EBSP. From the five simulated scenarios, scenario four (Fig. 2) showed the best goodness-of-fit to the observed data (Additional file 5) and had the highest posterior probability, 0.99 (CI 95 % = 0.97–1) when using the direct estimation method and 1 (CI 95 %–1–1) with the logistic regression method (Additional file 6). Regardless of the method used, the confidence intervals estimated from the competing scenarios did not overlap. Scenario four assumes that the turtle dove population suffered two variations in its effective population size. The Ne during the mid-late Pleistocene (~78,300 ybp) was estimated to be ~84,400 individuals. Following a demographic expansion during the early-mid Holocene (~7,600 ybp), the population reached an estimated Ne of ~3,160,000 doves. Finally, we found evidence of a recent contraction that started between the mid-Holocene and the present, resulting in a current Ne estimate of ~172,000 individuals (see Table 1 for more details on the estimated parameters). It is possible that the recent contraction in the Ne was not detected by the EBSP analysis because of lack of statistical power in the smaller dataset (30 loci). None of our analyses recovered evidence of a large Ne during the last interglacial (~120–140,000 ybp).
MaxEnt models achieved a good fit to the data of the present distribution (AUC value of 0.71 ± 0.01 SD). Bio1 (Annual mean temperature) was identified as the most important parameter (Percent contribution, PC: 47.9), and also achieved the highest importance in the permutation test (Permutation importance, PI: 45.5). Other key parameters were Bio18 (Precipitation of Warmest Quarter, PC: 24.2, PI: 22.3) and Bio12 (Annual Precipitation, PC: 23.4, PI: 19.5). Bio8 (Mean Temperature of Wettest Quarter) was of minor importance (PC: 4.5, PI: 12.7). The 10th percentile training presence logistic threshold was 0.36 and this value was applied to map and calculate the suitable habitat. The resulting fits were extrapolated to the past climate. According to the global climate model MIROC-ESM, the potentially suitable habitat during the breeding season was found to be largest during the mid Holocene (41 % of the total land area covered in the analysis) and smallest during the Last Glacial Maximum (19 %) (Fig. 4). Large surfaces of suitable habitat were also observed during the last interglacial (37 %) and present (39 %) (Fig. 4). Results obtained with alternative climate models showed similar outputs and are provided in the additional files section (Additional file 7).
We found genetic evidence of historical fluctuations in the effective population size of the turtle dove. As predicted, we inferred a drastic demographic expansion, estimated to have taken place during the early-mid Holocene. Our analyses also supported a relatively recent demographic contraction, which took place between the mid Holocene and the present. The demographic scenario supported by our genetic data is correlated with changes in the availability of suitable habitat for the turtle dove, as inferred by our ecological niche modelling results. We did not, however, detect evidence of population genetic structure associated to the alternative migratory flyways. Finally, we argue that the turtle dove’s fluctuating demography along with its decreasing effective population size make it particularly sensitive to current anthropogenic threats.
A fluctuating demographic history
The turtle dove is a warm-climate adapted species, with specific breeding habitat requirements that include a mixture of grasslands and woodlands . The climate oscillations that occurred during the Quaternary period strongly affected the distribution and abundance of this kind of environment . Therefore, we expected to find changes in the turtle dove effective population size in response to those shifts.
The last interglacial (~140–120,000 ybp) was a warm period, with climate conditions that were similar to those of the Holocene ; therefore, we expected to observe a large Ne during that time. Although our niche modelling showed a large surface of suitable turtle dove habitat during this period (Fig. 4), the genetic data did not support our prediction or was not informative about this period in time. We also expected to observe a drastic demographic contraction during the late Pleistocene, when a cooling process started ~75,000 ybp and lasted until the last glacial maximum . This cooling was coupled with shifts in precipitation patterns, creating an east–west gradient  and modifying habitat composition [5, 48]. The western Palearctic in general and the northwest Palearctic in particular, became dominated by dry and open-tundra like environments , unsuitable for most woodland dependent bird species. In agreement with the described scenario, we found a drastic reduction in the surface of available suitable habitat for the turtle dove (Fig. 4) and a very small estimated Ne (Table 1, Fig. 2). It is possible that this period represented a bottleneck for the species. If very few gene trees coalesce prior to the bottleneck, the demographic signal before the last glacial maximum would be erased. Such a bottleneck is a plausible explanation of why our data do not support models that predict a large Ne during the last interglacial.
The beginning of the Holocene is characterized by a rise in the average temperatures , which generated an increment in the suitable available habitat for most warm-temperature adapted species such as the turtle dove (Fig. 4, ). In agreement with this, our Tajima’s D and Fs analysis suggest a scenario of a recent demographic expansion, which, according to the EBSP analysis started during the early Holocene and reached the highest Ne in the present (Fig. 3). Our DIYABC analysis, however, supports a scenario in which the turtle dove went through a strong demographic expansion during the early-mid Holocene and reached its highest Ne during the mid Holocene (Fig. 2 and Table 1). This analysis also inferred the occurrence of a recent contraction in the Ne, indicating that the current Ne for the turtle dove is relatively low and closer to that observed during the late Pleistocene (Table 1). This finding fits our predictions; because it is coherent with the long-term abundance decline reported for the turtle dove . However, our estimate of when this Ne contraction started is imprecise and depends on several assumptions, including mutation rates, generation time and the parameters in our model. The observed contraction in the turtle dove effective size could have started at any point between the mid Holocene and the present (Table 1). Therefore, we are unable to assert that it is the strict consequence of modern human activities; moreover it also coincides with a slight reduction in the availability of suitable habitat for the species (projected data, Fig. 4).
Our analyses complement each other, and the observed differences are probably a consequence of the amount of genetic data each method is able to incorporate, its assumptions, and differences in the time scale at which mitochondrial and nuclear markers are able to provide resolution. Nonetheless, all our analyses coincide in showing a history of demographic fluctuations in the turtle dove. Overall, it is likely that this last contraction in the Ne is the outcome of the turtle dove’s fluctuating demographic history responding to natural climatic changes, currently enhanced by anthropogenic actions.
Population genetic structure: different flyways, the same gene pool
We hypothesized that the alternative flyways used by migrating turtle doves between Europe and Africa, would generate genetic differentiation among populations, a phenomenon previously described for other long distance migratory bird species worldwide (e.g. [20, 21, 52]). However, we did not detect any evidence of population genetic structure (Fig. 1b, c and Additional file 4). The turtle doves’ breeding range expands east of Bulgaria, through Russia, reaching China . Therefore we cannot comment on the possibility of genetic structure beyond the three European flyways, that were the focus of our study.
We propose two possible explanations for not finding genetic structure associated to the migratory flyways. 1) The described flyways for the turtle dove are coupled with a weak migratory connectivity (as defined in ), meaning that individuals from the same breeding population do not share the same wintering ground and vice versa. Therefore, these flyways would not constitute a barrier to gene flow. 2) Migratory connectivity could be strong, as suggested by the very low probability of recapturing western breeders in the eastern flyway and vice versa . If the latter held true, then flyways could generate genetic differentiation that we might not be able to detect, because divergence is too recent and individuals still share most of their alleles (i.e., high levels of incomplete lineage sorting). Our genetic data will not allow us to distinguish between these alternatives. However, the observation of one individual that was ringed as a western breeder but recaptured on eastern breeding grounds, together with the moderate probability of western and eastern breeders to be recaptured in the central flyway  would support the first explanation.
Implications for conservation
Our results suggest that the turtle dove’s demography make it particularly vulnerable to anthropogenic threats. Firstly, we found that it is undergoing a phase of effective population size contraction (Fig. 2, Table 1). This is of particular interest because the turtle dove has been recently categorized as vulnerable by the IUCN Red List (http://www.iucnredlist.org/), as a consequence of a global abundance decline (30–49 % in 15 years, ). This trend has been particularly dramatic in the United Kingdom (97 % decline in 40 years, ). The reported decline is attributed to anthropogenic degradation of its breeding  and wintering territories , also game hunting and poaching [57, 58]. Even though we cannot attribute our finding of a recent Ne contraction to modern human actions alone (discussed above), this is still relevant to turtle dove’s conservation. Many bird species that also suffered recent declines in abundance and are categorized as endangered (sensu IUCN Red List) show genomic signatures of long-term demographic declines that predate the appearance of modern anthropogenic threats . We also found signatures of a demographic history of fluctuations for the turtle dove. It has been observed that in demographically stable wild species the ratio of the effective and census population size (Ne/Nc) is around 0.1 , however several factors can produce deviations from this ratio . The estimated Nc for the turtle dove in Europe is between 6 and 12 million individuals  and our mean Ne estimate for the present is ~172,000 individuals (Table 1). Using the latter Nc and Ne parameters we observe a lower ratio than expected, Ne/Nc = 0.028–0.014, which is usually the result of drastic demographic fluctuations [62, 63]. Because factors such as sampling bias might affect the estimation of demographic parameters obtained using different analytical methods (e.g. ), we suggest to interpret our demographic parameters with caution and instead to focus on the general pattern. In this regard, species that exhibit the described drastic natural fluctuations could be particularly sensitive to anthropogenic threats, especially when they are going through a contraction phase (e.g. North American passenger pigeon, Ectopistes migratorius, , and Rocky Mountain grasshopper, Melanoplus spretus, ).
The described demographic scenario (i.e. long-term decline and drastic fluctuations) could increase species’ vulnerability to the impacts of macroecological changes resulting from natural earth cycles, which are currently enhanced by anthropogenic activities [67, 68]. To further illustrate this, a paradigmatic case is the North American passenger pigeon; which went extinct almost a century ago, although it was considered one of the most abundant bird species in the world . A genomic study of this species suggested that its naturally fluctuating demographic history, coupled with a downward trend in its Ne that coincided with the appearance of human threats were the main causes of its sudden extinction .
Our finding that the turtle dove shows the genomic signatures of a fluctuating demographic history coupled with a recent effective population size contraction, provides support for its vulnerable conservation status. Finally, because we found no evidence of population genetic structure, we prompt that all turtle dove populations across Europe (the portion of the range that was sampled) are equally relevant for conservation.
Hewitt G. The genetic legacy of the Quaternary ice ages. Nature. 2000;405:907–13. doi:10.1038/35016000.
Nyári ÁS, Joseph L. Comparative phylogeography of Australo-Papuan mangrove-restricted and mangrove-associated avifaunas. Biol J Lin Soc. 2013;109:574–98. doi:10.1111/bij.12082.
Carnaval AC, Waltari E, Rodrigues MT, Rosauer D, et al. Prediction of phylogeographic endemism in an environmentally complex biome. Proc R Soc B. 2014;281:20141461. doi:10.1098/rspb.2014.1461.
Stewart JR, Lister AM, Barnes I, Dalén L. Refugia revisited: individualistic responses of species in space and time. Proc R Soc B. 2010;277:661–71. doi:10.1098/rspb.2009.1272.
Holm SR, Svenning JC. 180,000 years of climate change in Europe: Avifaunal responses and vegetation implications. PLoS One. 2014;9:e94021. doi:10.1371/journal.pone.0094021.
Hewitt G. Genetic consequences of climatic oscillations in the Quaternary. Phil Trans R Soc Lond B. 2004;359:183–95. doi:10.1098/rstb.2003.1388.
Schmitt T. Molecular biogeography of Europe: Pleistocene cycles and postglacial trends. Front Zool. 2007;4:1–13. doi:10.1186/1742-9994-4-11.
Bennet KD, Provan J. What do we mean by ‘refugia’? Quaternary Sci Rev. 2008;27:2449–55. doi:10.1016/j.quascirev.2008.08.019.
Hitch AT, Leberg PL. Breeding distributions of North American bird species moving north as a result of climate change. Conserv Biol. 2007;21:534–9. doi:10.1111/j.1523-1739.2006.00609.x.
Peterson AT, Navarro-Sigüenza AG, Martínez-Meyer E, Cuervo-Robayo AP, Berlanga H, Soberón J. Twentieth century turnover of Mexican endemic avifaunas: Landscape change versus climate drivers. Sci Adv. 2015;1:e1400071. doi:10.1126/sciadv.1400071.
Cox GW. Bird migration and global change. Washington, DC: Island Press; 2010.
Newton I. The migration ecology of birds. London: Academic; 2010.
Fontaine JJ, Stutzman RJ, Gannes LZ, Fontaine JJ, Stutzman RJ, Gannes LZ. Leaps, chains and climate change for western migratory songbirds. Phenological synchrony and bird migration: changing climate and seasonal resources in north America. Boca Raton: CRC Press; 2015. p. 8–11.
Newton I. Population limitation in migrants. Ibis. 2004;146:197–226. doi:10.1111/j.1474-919X.2004.00293.x.
Sanderson FJ, Donald PF, Pain DJ, Burfield IJ, van Bommel FPJ. Long-term population declines in Afro-Palearctic migrant birds. Biol Conserv. 2006;131:93–105. doi:10.1016/j.biocon.2006.02.008.
Hahn S, Bauer S, Liechti F. The natural link between Europe and Africa – 2.1 billion birds on migration. Oikos. 2009;118:624–6. doi:10.1111/j.1600-0706.2008.17309.x.
Bruderer B, Liechti F. Bird migration across the Mediterranean. Proc Int Ornithol Congr. 1999;22:1983–99.
Bairlein F. Results of bird ringing in the study of migration routes and behaviour. Ardea. 2001;89:7–19.
Trierweiler C, Klaassen RH, Drent RH, Exo KM, Komdeur J, Bairlein F, Koks BJ. Migratory connectivity and population-specific migration routes in a long-distance migratory bird. Proc R Soc B. 2014;281:20132897. doi:10.1098/rspb.2013.2897.
Ruegg KC, Smith TB. Not as the crow flies: a historical explanation for circuitous migration in Swainson’s thrush (Catharus ustulatus). Proc R Soc Lond B. 2002;269:1375–81. doi:10.1098/rspb.2002.2032.
Alcaide M, Scordato ES, Price TD, Irwin DE. Genomic divergence in a ring species complex. Nature. 2014;511:83–5. doi:10.1038/nature13285.
Ruegg K. Genetic, morphological, and ecological characterization of a hybrid zone that spans a migratory divide. Evolution. 2008;62:452–66. doi:10.1111/j.1558-5646.2007.00263.x.
Irwin DE, Irwin JH. Siberian migratory divides: the role of seasonal migration in speciation. In: Greenberg R, Marra PP, editors. Birds of two worlds: the ecology and evolution of bird migration. Baltimore: Johns Hopkins University Press; 2005. p. 27–40.
Browne SJ, Aebischer NJ. Habitat use, foraging ecology and diet of Turtle Doves Streptopelia turtur in Britain. Ibis. 2003;145:572–82. doi:10.1046/j.1474-919X.2003.00185.x.
Jarry G. Tourterelle des Bois Streptopelia turtur. In: Yeatman-Berthelot D, Jarry G, editors. Nouvel atlas des oiseaux nicheurs de France. Paris: Société Ornithologique de France; 1995. p. 380–3.
Marx M, Korner-Nievergelt F, Quillfeldt P. Analysis of Ring Recoveries of European Turtle Doves Streptopelia turtur-Flyways, Migration Timing and Origin Areas of Hunted Birds. Acta Ornithologica. 2016;51:55–70. http://dx.doi.org/10.3161/00016454AO2016.51.1.005.
BirdLife International. Streptopelia turtur. The IUCN Red List of Threatened Species 2015: e.T22690419A82781167. 2015 http://dx.doi.org/10.2305/IUCN.UK.2015-4.RLTS.T22690419A82781167.en.
Peterson BK, Weber JN, Kay EH, et al. Double digest RADseq: an inexpensive method for de novo SNP discovery and genotyping in model and non-model species. PLoS One. 2012;7:e37135. doi:10.1371/journal.pone.0037135.
Catchen J, Amores A, Hohenlohe P, et al. Stacks: building and genotyping loci de novo from short-read sequences. G3 (Bethesda). 2011;1:171–82. doi:10.1534/g3.111.000240.
Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10:R25. doi:10.1186/gb-2009-10-3-r25.
Librado P, Rozas J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009;25:1451–2. doi:10.1093/bioinformatics/btp187.
Leigh J, Bryant D. Popart: full-feature software for haplotype network construction. Methods Ecol Evol. 2015;6:1110–6.
Excoffier L, Lischer HE. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Res. 2010;10:564–7. doi:10.1111/j.1755-0998.2010.02847.x.
Raj A, Stephens M, Pritchard JK. fastSTRUCTURE: variational inference of population structure in large SNP data sets. Genetics. 2014;197:573–89. doi:10.1534/genetics.114.164350.
Rosenberg NA. DISTRUCT: a program for the graphical display of population structure. Mol Ecol Notes. 2004;4:137–8. doi:10.1046/j.1471-8286.2003.00566.x.
Dray S, Dufour AB. The ade4 package: implementing the duality diagram for ecologists. J Stat Softw. 2007;22:1–20. doi: 10.18637/jss.v022.i04.
Foll M, Gaggiotti O. A genome-scan method to identify selected loci appropriate for both dominant and codominant markers: a Bayesian perspective. Genetics. 2008;180:977–93. doi:10.1534/genetics.108.092221.
Pfeifer B, Wittelsbürger U, Ramos-Onsins SE, Lercher MJ. PopGenome: An Efficient Swiss Army Knife for Population Genomic Analyses in R. Mol Biol Evol. 2014;31:1929–36. doi:10.1093/molbev/msu136.
Wright S. Evolution in mendelian populations. Genetics. 1930;16:97–158.
Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012;29:1969–73. doi:10.1093/molbev/mss075.
Weir JT, Schluter D. Calibrating the avian molecular clock. Mol Ecol. 2008;17:2321–8. doi:10.1111/j.1365-294X.2008.03742.x.
Rambaut A, Suchard MA, Xie D, Drummond AJ. Tracer v1.6. 2014. Available from: http://beast.bio.ed.ac.uk/Tracer.
Cornuet JM, Pudlo P, Veyssier J, Dehne-Garcia A, et al. DIYABC v2. 0: A software to make approximate Bayesian computation inferences about population history using single nucleotide polymorphism, DNA sequence and microsatellite data. Bioinformatics. 2014;30:1187–9. doi:10.1093/bioinformatics/btt763.
Cornuet JM, Ravigné V, Estoup A. Inference on population history and model checking using DNA sequence and microsatellite data with the software DIYABC (v1. 0). BMC Bioinformatics. 2010;11:401. doi:10.1186/1471-2105-11-401.
Glutz Von Blotzheim UN. Handbuch der Vögel Mitteleuropas. Wiesbaden: Akademische Verlagsgesellschaft; 1980.
Phillips SJ, Anderson RP, Schapire RE. Maximum entropy modeling of species geographic distributions. Ecol Model. 2006;190:231–59. doi:10.1016/j.ecolmodel.2005.03.026.
Prentice IC, Jolly D, et al. Mid-Holocene and glacial-maximum vegetation geography of the northern continents and Africa. J Biogeogr. 2000;27:507–19. doi:10.1046/j.1365-2699.2000.00425.x.
Van Andel TH, Tzedakis PC. Palaeolithic landscapes of Europe and environs, 150,000–25,000 years ago: an overview. Quart Sci Rev. 1996;15:481–500. doi:10.1016/0277-3791(96)00028-5.
Ehlers J, Gibbard PL, Hughes PD. Quaternary glaciations-extent and chronology: a closer look. Amsterdam: Elsevier; 2011. p. 1–13.
Dansgaard W, Johnsen SJ, Clausen HB, et al. Evidence of general instability of past climate from a 250-kyr ice-core record. Nature. 1993;364:218–20. doi:10.1038/364218a0.
Mackay A, Battarbee R, Birks J, Oldfield F. Global change in the Holocene. London: Arnold; 2003. p. 1–3.
Bearhop S, Fiedler W, Furness RW, Votier SC, Waldron S, Newton J, Bowen GJ, Berthold P, Farnsworth K. Assortative mating as a mechanism for rapid evolution of a migratory divide. Science. 2005;310:502–4. doi:10.1126/science.1115661.
Webster MS, Marra PP, Haig SM, Bensch S, Holmes RT. Links between worlds: unraveling migratory connectivity. Trends Ecol Evol. 2002;17:76–83. doi:10.1016/S0169-5347(01)02380-1.
Hayhow DB, Bond AL, Eaton MA, Grice PV, Hall C, Hall J, Harris SJ, Hearn RD, Holt CA, Noble DG, Stroud DA, Wotton S. The state of the UK’s birds 2015. Sandy: RSPB, BTO, WWT, JNCC, NE, NIEA, NRW and SNH; 2015.
Browne SJ, Aebischer NJ. Temporal changes in the breeding ecology of European Turtle Doves Streptopelia turtur in Britain, and implications for conservation. Ibis. 2004;146:125–37. doi:10.1111/j.1474-919X.2004.00235.x.
Eraud C, Boutin J, Riviere M, Brun J, Barbraud C, Lormee H. Survival of Turtle Doves Streptopelia turtur in relation to western Africa environmental conditions. Ibis. 2009;151:186–90.
Browne SJ, Aebischer NJ. Studies of West Palearctic birds: turtle dove. British Birds. 2005;98:58–72.
Lutz M. Management plan for Turtle Dove (Streptopelia turtur) 2007–2009. Luxemburg: Office for Official Publications of the European Communities; 2007.
Nadachowska-Brzyska K, Li C, Smeds L, Zhang G, Ellegren H. Temporal dynamics of avian populations during Pleistocene revealed by whole-genome sequences. Curr Biol. 2015;25:1375–80. doi:10.1016/j.cub.2015.03.047.
Frankham R. Effective population size/adult population size ratios in wildlife: a review. Genet Res. 1995;66:95–107. doi:10.1017/S0016672300034455.
Turner TF, Osborne MJ, Moyer GR, Benavides MA, Alo D. Life history and environmental variation interact to determine effective population to census size ratio. Proc R Soc B. 2006;273:3065–73. doi:10.1098/rspb.2006.367.
Nei MMT, Chakraborty R. The bottleneck effect and genetic variability in populations. Evolution. 1975;29:1–10. doi:10.2307/2407137.
Motro U, Thomson G. On heterozygosity and the effective size of populations subject to size changes. Evolution. 1982;36:1059–66. doi:10.2307/2408083.
Heller R, Chikhi L, Siegismund HR. The confounding effect of population structure on Bayesian skyline plot inferences of demographic history. PLoS One. 2013;8:e62992. http://dx.doi.org/10.1371/journal.pone.0062992.
Hung CM, Shaner PJL, Zink RM, Liu WC, Chu TC, Huang WS, Li SH. Drastic population fluctuations explain the rapid extinction of the passenger pigeon. Proc Natl Acad Sci U S A. 2014;111:10636–41. doi:10.1073/pnas.1401526111.
Lockwood JA, DeBrey LD. A solution for the sudden and unexplained extinction of the Rocky Mountain grasshopper (Orthoptera: Acrididae). Environ Entomol. 1990;19:1194–205.
Parmesan C. Ecological and evolutionary responses to recent climate change. Ann Rev Ecol Evol Syst. 2006;37:637–69.
Chen IC, Hill JK, Ohlemüller R, Roy DB, Thomas CD. Rapid range shifts of species associated with high levels of climate warming. Science. 2011;333:1024–6. doi:10.1126/science.1206432.
Schorger AW. The passenger pigeon: its natural history and extinction. Madison: Univ. of Wisconsin Press; 1955.
Would like to thank Tanja Romeike, Steffen Koschkar and Fernando Spina for help in the field. We are indebted to the members of the Fuller Evolutionary Biology program, especially Dr. Irby Lovette and Dr. Bronwyn Butcher for access to the research facilities at the Cornell Lab of Ornithology and for providing valuable feedback for this project.
This study was supported by the Deutsche Forschungsgemeinschaft (DFG, Germany, Heisenberg grant to PQ, QU148/5), and also by a grant from the Research fund (Forschungsförderfonds) of the Justus Liebig University (Germany).
LuC carried out laboratory work to generate the cytochrome b and SNP datasets, also performed bioinformatics analyses on both datasets, and wrote the manuscript; LeC provided guidance and participated in laboratory work to generate the SNP dataset and performed bioinformatics analyses on that data, also collaborated in writing the manuscript; TW helped with project design, provided the facilities to generate the cytochrome b dataset and drafted the manuscript. LH, EC, DJ, RG, ZP, BDE, MB, CJG and MM conducted field work and provided blood and tissue samples for this work; they also drafted the manuscript. PQ designed and coordinated the project, participated in the field work, performed the niche modeling simulation analysis and helped write the manuscript. All authors gave their final approval for publication.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval approval and consent to participate
We used mist nets to catch individuals for blood sampling, under the licences obtained from the national authorities of each country. The birds were bled, ringed and released as fast as possible in order to reduce stress. Birds’ wings hunted during migration were provided by hunting estates. No birds were hunted for the purpose of the present study.
Data used in MaxEnt was downloaded from open access sources, namely: http://www.worldclim.org (climate data) and http://www.gbif.org (occurrence data). mtDNA sequences are available in GenBank (accession numbers for mtDNA: KU588290-KU588384) and Rad-seq data is available in DRYAD (doi:10.5061/dryad.p5cf4).
Additional file 1:
Details of the samples used. A table with additional information on the blood and tissue samples obtained for our work. (DOC 143 kb)
Additional file 2:
Further details on materials and methods. Additional information on the laboratory methods related to the RAD-seq procedures and mtDNA sequencing. Also, further details on the demographic history modeling performed with DIYABC, and the niche modeling analysis conducted in MAXENT. (DOC 30 kb)
Additional file 3:
Analysis of molecular variance (AMOVAs). Results of the AMOVAs performed, testing for significant genetic differences among the three main flyways (using mitochondrial and nuclear DNA data). (DOC 41 kb)
Additional file 4:
fastSTRUCTURE analysis indicating the number of genetic clusters suggested for the turtle dove, K = 1. Barplot graphics representing K = 1 and 3 are shown. (DOC 393 kb)
Additional file 5:
Principal components analysis to test the goodness-of-fit of each evaluated scenario against the simulated data in DIYABC. PCA plots show that Scenario 4 had the best goodness-of-fit among all tested scenarios. (DOC 630 kb)
Additional file 6:
Posterior probability of scenario choice in DIYABC. Direct estimation approach that uses 0.1 % of the closest simulations to the observed data. (b) The logistic regression uses 1 % of the closest simulations. (DOC 211 kb)
Additional file 7:
MaxEnt habitat suitability models for European Turtle doves. Niche modelling analysis based on additional climate models: CC (CCSM4), MR (MIROC-ESM) and MI (MPI-ESM-P). Models available on http://www.worldclim.org/. (DOC 729 kb)
Rights and permissions
Open Access This 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.
About this article
Cite this article
Calderón, L., Campagna, L., Wilke, T. et al. Genomic evidence of demographic fluctuations and lack of genetic structure across flyways in a long distance migrant, the European turtle dove. BMC Evol Biol 16, 237 (2016). https://doi.org/10.1186/s12862-016-0817-7
- Population genetic structure
- Migratory birds
- Paleoclimatic niche modelling
- Climate change