Genomic evidence of demographic fluctuations and lack of genetic structure across flyways in a long distance migrant, the European turtle dove
© The Author(s). 2016
Received: 1 August 2016
Accepted: 28 October 2016
Published: 7 November 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.
KeywordsPopulation genetic structure Demography Migratory birds Genomics Conservation Paleoclimatic niche modelling Climate change
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 .
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.
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
Demographic parameters inferred with DIYABC
Late Pleistocene (t2)
Mid Holocene (t1)
Late Pleistocene Ne (N2)
Mid Holocene Ne (N1)
Present Ne (N3)
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.
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:https://doi.org/10.5061/dryad.p5cf4).
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. doi:https://doi.org/10.1038/35016000.View ArticlePubMedGoogle Scholar
- 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:https://doi.org/10.1111/bij.12082.View ArticleGoogle Scholar
- 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:https://doi.org/10.1098/rspb.2014.1461.View ArticlePubMedPubMed CentralGoogle Scholar
- 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:https://doi.org/10.1098/rspb.2009.1272.View ArticlePubMedGoogle Scholar
- Holm SR, Svenning JC. 180,000 years of climate change in Europe: Avifaunal responses and vegetation implications. PLoS One. 2014;9:e94021. doi:https://doi.org/10.1371/journal.pone.0094021.View ArticlePubMedPubMed CentralGoogle Scholar
- Hewitt G. Genetic consequences of climatic oscillations in the Quaternary. Phil Trans R Soc Lond B. 2004;359:183–95. doi:https://doi.org/10.1098/rstb.2003.1388.View ArticleGoogle Scholar
- Schmitt T. Molecular biogeography of Europe: Pleistocene cycles and postglacial trends. Front Zool. 2007;4:1–13. doi:https://doi.org/10.1186/1742-9994-4-11.View ArticleGoogle Scholar
- Bennet KD, Provan J. What do we mean by ‘refugia’? Quaternary Sci Rev. 2008;27:2449–55. doi:https://doi.org/10.1016/j.quascirev.2008.08.019.View ArticleGoogle Scholar
- 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:https://doi.org/10.1111/j.1523-1739.2006.00609.x.View ArticlePubMedGoogle Scholar
- 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:https://doi.org/10.1126/sciadv.1400071.View ArticlePubMedPubMed CentralGoogle Scholar
- Cox GW. Bird migration and global change. Washington, DC: Island Press; 2010.Google Scholar
- Newton I. The migration ecology of birds. London: Academic; 2010.Google Scholar
- 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.Google Scholar
- Newton I. Population limitation in migrants. Ibis. 2004;146:197–226. doi:https://doi.org/10.1111/j.1474-919X.2004.00293.x.View ArticleGoogle Scholar
- 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:https://doi.org/10.1016/j.biocon.2006.02.008.View ArticleGoogle Scholar
- 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:https://doi.org/10.1111/j.1600-0706.2008.17309.x.View ArticleGoogle Scholar
- Bruderer B, Liechti F. Bird migration across the Mediterranean. Proc Int Ornithol Congr. 1999;22:1983–99.Google Scholar
- Bairlein F. Results of bird ringing in the study of migration routes and behaviour. Ardea. 2001;89:7–19.Google Scholar
- 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:https://doi.org/10.1098/rspb.2013.2897.View ArticlePubMedPubMed CentralGoogle Scholar
- 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:https://doi.org/10.1098/rspb.2002.2032.View ArticleGoogle Scholar
- Alcaide M, Scordato ES, Price TD, Irwin DE. Genomic divergence in a ring species complex. Nature. 2014;511:83–5. doi:https://doi.org/10.1038/nature13285.View ArticlePubMedGoogle Scholar
- Ruegg K. Genetic, morphological, and ecological characterization of a hybrid zone that spans a migratory divide. Evolution. 2008;62:452–66. doi:https://doi.org/10.1111/j.1558-5646.2007.00263.x.View ArticlePubMedGoogle Scholar
- 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.Google Scholar
- Browne SJ, Aebischer NJ. Habitat use, foraging ecology and diet of Turtle Doves Streptopelia turtur in Britain. Ibis. 2003;145:572–82. doi:https://doi.org/10.1046/j.1474-919X.2003.00185.x.View ArticleGoogle Scholar
- 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.Google Scholar
- 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:https://doi.org/10.1371/journal.pone.0037135.View ArticlePubMedPubMed CentralGoogle Scholar
- 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:https://doi.org/10.1534/g3.111.000240.View ArticleGoogle Scholar
- 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:https://doi.org/10.1186/gb-2009-10-3-r25.View ArticlePubMedPubMed CentralGoogle Scholar
- Librado P, Rozas J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009;25:1451–2. doi:https://doi.org/10.1093/bioinformatics/btp187.View ArticlePubMedGoogle Scholar
- Leigh J, Bryant D. Popart: full-feature software for haplotype network construction. Methods Ecol Evol. 2015;6:1110–6.View ArticleGoogle Scholar
- 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:https://doi.org/10.1111/j.1755-0998.2010.02847.x.View ArticleGoogle Scholar
- Raj A, Stephens M, Pritchard JK. fastSTRUCTURE: variational inference of population structure in large SNP data sets. Genetics. 2014;197:573–89. doi:https://doi.org/10.1534/genetics.114.164350.View ArticlePubMedPubMed CentralGoogle Scholar
- Rosenberg NA. DISTRUCT: a program for the graphical display of population structure. Mol Ecol Notes. 2004;4:137–8. doi:https://doi.org/10.1046/j.1471-8286.2003.00566.x.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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:https://doi.org/10.1534/genetics.108.092221.View ArticlePubMedPubMed CentralGoogle Scholar
- 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:https://doi.org/10.1093/molbev/msu136.View ArticlePubMedPubMed CentralGoogle Scholar
- Wright S. Evolution in mendelian populations. Genetics. 1930;16:97–158.Google Scholar
- 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:https://doi.org/10.1093/molbev/mss075.View ArticlePubMedPubMed CentralGoogle Scholar
- Weir JT, Schluter D. Calibrating the avian molecular clock. Mol Ecol. 2008;17:2321–8. doi:https://doi.org/10.1111/j.1365-294X.2008.03742.x.View ArticlePubMedGoogle Scholar
- 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:https://doi.org/10.1093/bioinformatics/btt763.View ArticlePubMedGoogle Scholar
- 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:https://doi.org/10.1186/1471-2105-11-401.View ArticlePubMedPubMed CentralGoogle Scholar
- Glutz Von Blotzheim UN. Handbuch der Vögel Mitteleuropas. Wiesbaden: Akademische Verlagsgesellschaft; 1980.Google Scholar
- Phillips SJ, Anderson RP, Schapire RE. Maximum entropy modeling of species geographic distributions. Ecol Model. 2006;190:231–59. doi:https://doi.org/10.1016/j.ecolmodel.2005.03.026.View ArticleGoogle Scholar
- 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:https://doi.org/10.1046/j.1365-2699.2000.00425.x.View ArticleGoogle Scholar
- 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:https://doi.org/10.1016/0277-3791(96)00028-5.View ArticleGoogle Scholar
- Ehlers J, Gibbard PL, Hughes PD. Quaternary glaciations-extent and chronology: a closer look. Amsterdam: Elsevier; 2011. p. 1–13.View ArticleGoogle Scholar
- 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:https://doi.org/10.1038/364218a0.View ArticleGoogle Scholar
- Mackay A, Battarbee R, Birks J, Oldfield F. Global change in the Holocene. London: Arnold; 2003. p. 1–3.Google Scholar
- 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:https://doi.org/10.1126/science.1115661.View ArticlePubMedGoogle Scholar
- Webster MS, Marra PP, Haig SM, Bensch S, Holmes RT. Links between worlds: unraveling migratory connectivity. Trends Ecol Evol. 2002;17:76–83. doi:https://doi.org/10.1016/S0169-5347(01)02380-1.View ArticleGoogle Scholar
- 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.Google Scholar
- 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:https://doi.org/10.1111/j.1474-919X.2004.00235.x.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Browne SJ, Aebischer NJ. Studies of West Palearctic birds: turtle dove. British Birds. 2005;98:58–72.Google Scholar
- Lutz M. Management plan for Turtle Dove (Streptopelia turtur) 2007–2009. Luxemburg: Office for Official Publications of the European Communities; 2007.Google Scholar
- 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:https://doi.org/10.1016/j.cub.2015.03.047.View ArticlePubMedPubMed CentralGoogle Scholar
- Frankham R. Effective population size/adult population size ratios in wildlife: a review. Genet Res. 1995;66:95–107. doi:https://doi.org/10.1017/S0016672300034455.View ArticleGoogle Scholar
- 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:https://doi.org/10.1098/rspb.2006.367.View ArticlePubMedPubMed CentralGoogle Scholar
- Nei MMT, Chakraborty R. The bottleneck effect and genetic variability in populations. Evolution. 1975;29:1–10. doi:https://doi.org/10.2307/2407137.View ArticleGoogle Scholar
- Motro U, Thomson G. On heterozygosity and the effective size of populations subject to size changes. Evolution. 1982;36:1059–66. doi:https://doi.org/10.2307/2408083.View ArticleGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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:https://doi.org/10.1073/pnas.1401526111.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticleGoogle Scholar
- Parmesan C. Ecological and evolutionary responses to recent climate change. Ann Rev Ecol Evol Syst. 2006;37:637–69.View ArticleGoogle Scholar
- 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:https://doi.org/10.1126/science.1206432.View ArticlePubMedGoogle Scholar
- Schorger AW. The passenger pigeon: its natural history and extinction. Madison: Univ. of Wisconsin Press; 1955.Google Scholar