Next Article in Journal
Influence of Landscape Factors on Amphibian Roadkills at the National Level
Next Article in Special Issue
High Genetic Diversity among Breeding Red-Backed Shrikes Lanius collurio in the Western Palearctic
Previous Article in Journal
After the Fall: Legacy Effects of Biogenic Structure on Wind-Generated Ecosystem Processes Following Mussel Bed Collapse
Previous Article in Special Issue
Winter Roost Tree Selection and Phenology of the Long-Eared Owl (Asio otus) in Crimea
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Pliocene Origin, Ice Ages and Postglacial Population Expansion Have Influenced a Panmictic Phylogeography of the European Bee-Eater Merops apiaster

1
Department of Biology, Institute of Pharmacy and Molecular Biotechnology, Heidelberg University, 69120 Heidelberg, Germany
2
Bee-Eater Study Group of the DO-G, Geschwister-Scholl-Str. 15, 67304 Kerzenheim, Germany
*
Authors to whom correspondence should be addressed.
Diversity 2019, 11(1), 12; https://doi.org/10.3390/d11010012
Submission received: 20 August 2018 / Revised: 29 November 2018 / Accepted: 26 December 2018 / Published: 15 January 2019

Abstract

:
Oscillations of periods with low and high temperatures during the Quaternary in the northern hemisphere have influenced the genetic composition of birds of the Palearctic. During the last glaciation, ending about 12,000 years ago, a wide area of the northern Palearctic was under lasting ice and, consequently, breeding sites for most bird species were not available. At the same time, a high diversity of habitats was accessible in the subtropical and tropical zones providing breeding grounds and refugia for birds. As a result of long-term climatic oscillations, the migration systems of birds developed. When populations of birds concentrated in refugia during ice ages, genetic differentiation and gene flow between populations from distinct areas was favored. In the present study, we explored the current genetic status of populations of the migratory European bee-eater. We included samples from the entire Palearctic-African distribution range and analyzed them via mitochondrial and nuclear DNA markers. DNA data indicated high genetic connectivity and panmixia between populations from Europe, Asia and Africa. Negative outcomes of Fu’s Fs and Tajima’s D tests point to recent expansion events of the European bee-eater. Speciation of Merops apiaster started during the Pliocene around three million years ago (Mya), with the establishment of haplotype lineages dated to the Middle Pleistocene period circa 0.7 Mya. M. apiaster, which breed in Southern Africa are not distinguished from their European counterparts, indicating a recent separation event. The diversification process of the European bee-eater was influenced by climatic variation during the late Tertiary and Quaternary. Bee-eaters must have repeatedly retracted to refugia in the Mediterranean and subtropical Africa and Asia during ice ages and expanded northwards during warm periods. These processes favored genetic differentiation and repeated lineage mixings, leading to a genetic panmixia, which we still observe today.

1. Introduction

Abiotic and biotic conditions have an important influence on the genetic structure of resident and long-distance migratory birds. When physical and/or ecological barriers limit gene flow between populations, phylogeographic differentiation can occur. Without barriers, gene flow will take place between populations. When differentiated lineages meet in refugia, lineages can mix, leading to genetic panmixia. These scenarios have been observed in phylogeographic studies of several migratory birds [1,2,3]. Philopatry, habitat fragmentation, specific migratory flyways, and climatic conditions can influence the speciation of migratory taxa by limiting gene flow between populations and promoting genetic divergence [4]. Species with differing migratory routes and destinations are usually genetically distinct, while species with similar breeding and wintering areas show weak genetic structure [2,5].
In the last two million years, a cyclical change of climate was observed in the northern hemisphere, with long warm and long cold periods. The Pleistocene climatic oscillations have led to compressions and expansions of the distribution range of many Palearctic bird taxa. During ice ages, populations retracted to distinct refugia in the Mediterranean, African and Asian subtropics. During warm periods, they left the refugia and distributed themselves over wide areas of the Palearctic again. These frequent cycles promoted genetic differentiation, gene flow and even panmixia [6]. These environmental conditions have had important impacts not only on the bird richness and distribution but also on their migration systems [7]. Thus, the migration systems of birds developed in response to long-term climatic cycles. They determined the phylogeographic pattern in many species [8].
The bee-eaters (family Meropidae) consist of widely distributed species in Africa, Asia, Europe, Australia, and Indonesia, inhabiting a variety of habitats from arid areas to wet forests [9]. Phylogenetic relationships between the bee-eaters have been explored on the basis of DNA data [10] and morphological diagnosis [11]. The phylogenetic analyses of all bee-eaters, based on mitochondrial and nuclear data, revealed two main clades, which differ in their migratory behavior: one lineage comprising basically sedentary species and one including mostly migratory species [10].
The long-distance migrant European bee-eater Merops apiaster exhibits a wide distribution in the Palearctic. Its breeding range covers a large area in Europe, extending to Western Asia, as well as the northwest and the extreme south of Africa [12]. Although this species usually inhabits mainly dry areas, populations of bee-eaters have expanded towards Central/Southeast Europe, north of the Alps and Great Britain during the last decades, probably due to the occurrence of warmer summers [13].
The phylogenetic position and the genetic variability of the European bee-eater have been studied before [9,10,11,14]. A low level of genetic inter-population divergence in M. apiaster had been observed [10]. However, the detailed evolutionary history of M. apiaster and the clarification of the phylogenetic relationships between the bee-eaters is still a demand. Therefore, in this study we investigated sequences of mitochondrial marker genes and nuclear microsatellite markers: (i) to explore the evolutionary history of the migratory European bee-eater, (ii) to trace the genetic connectivity between breeding populations in Europe, Asia and Africa, and (iii) to estimate the speciation time of M. apiaster. Special emphasis was placed on populations in Germany (Rhineland-Palatinate, Baden-Württemberg, Saxony and Saxony-Anhalt), which have been established during the last two decades, and on populations which breed in Southern Africa.

2. Materials and Methods

2.1. Taxon Sampling, PCR Amplification, and Sequencing

We sampled a total of 279 individuals of Merops apiaster from 26 localities (Table S1, Genbank accession numbers MH217684—MH218120). All the samples were deposited in the bird tissue collection of the Institute of Pharmacy and Molecular Biotechnology (IPMB) or loaned by the following collections and researchers: Claire Spottiswood; Field Museum of Natural History, USA (FMNH); Louisiana State University Museum of Natural History, USA (LSUMNS); Natural History Museum of Denmark, Copenhagen, Denmark (ZMUC); Zoological Museum of Bashkir University, Russia (ZMBU); Zoological Museum of Belarusian University, Azerbaijan. Additional DNA sequences from GenBank were included in the analysis (Accession numbers: KU984741-KU984911, KU957811-KU957981, EU021509, EU021514-EU021517, EU021519-EU021521, EU021524, EU021527, EU021528, EU021530, EU021536, EU021537, KJ455498, KJ455499).
The DNA was extracted from tissue, feathers, blood and swab samples by the standard phenol-chloroform protocol [15]. The mitochondrial cytochrome b (Cyt b) and cytochrome oxidase I (COI) genes were amplified using the primers published from previous studies [3,16,17]. PCR conditions were set as follows: 94 °C denaturation for 5 min, 52–62 °C for annealing temperature for 45 s, and 72 °C for extension for 2 min, running in 38 cycles. A final extension step was conducted at 72 °C for 10 min. The reaction was executed in a 40 µL final volume, each reaction containing 60 µg of DNA, 38.7 µL of H2O, 100 µM of dNTP, PCR Solution Buffer Complete (10x), 10 pmol of forward and reverse primers, and 2 units of Taq Polymerase.
The PCR product was purified by precipitation in 4 M NH4Ac and absolute ethanol (1:1:10), centrifuged at 13,000 rpm for a period of 30 min, followed by centrifugation with 70% ethanol using the same settings and later dissolved in 25 µL of sterile H2O. For sequencing, 1.0 µL of sequencing primer (10 pmol/µL) was combined with 7.0 µL of PCR products. The Sanger sequencing was executed on an ABI 3730 automated capillary sequencer (Applied Biosystems, Carlsbad, CA, USA) with the ABI Prism Big Dye Terminator Cycle Sequencing Ready Reaction Kit 3.1 (carried out by STARSEQ GmbH, Mainz, Germany). The nucleotide sequences were aligned using the program CodonCode aligner version 5 (http://www.codoncode.com/aligner). Additionally, we refined the alignment by eye and checked for stop codons in the sequences to exclude the possibility that nuclear copies had been amplified.

2.2. Demographic Analysis

DNAsp 5.0 [18] was used to access the demographic parameters for each population: number of haplotypes (H), number of polymorphic (segregating) sites (s), haplotype diversity (Hd), nucleotide diversity (Pi), Fu and Li’s F, Tajima’s D and Fu’s Fs statistics. To test for signals of sudden demographic expansion [19], the mismatch distribution was estimated using the constant population size model and the evaluation of the data was verified through a goodness-of-fit test between the observed and simulated data (Raggedness statistic r) [20].

2.3. Phylogenetic Analyses and Molecular Dating

To employ the best evolutionary model for each gene we used jModelTest 2.1.7 [21]. HYK + InvGamma model is selected by the Akaike Information Criterion (AIC) to the Cyt b gene. The molecular timing tree estimation was assessed based on the evolutionary rate of the Cyt b gene using its known evolutionary rate (2.1% substitutions/site/lineage/million years) [22], set to strict model and Yule tree, the length of chain set to 107, and sampling every 1000th tree in Beauti 1.8.0 [23]. The phylogenetic tree was performed in Beast 1.8.0 [23] and the stationarity and convergence of runs were checked using Tracer 1.5 [24]. The tree with the best likelihood was summarized using TreeAnnotator 1.8.0 [23]. Other species belonging to the Merops genus were also included in the molecular dating phylogenetic tree: M. albicollis, M. boehmi, M. bullockoides, M. gularis, M. hirundineus, M. leschenaulti, M. malimbicus, M. nubicus, M. nubicoides, M. orientalis, M. persicus, and M. viridis [10]. The analysis was conducted on the Cipres platform [25].

2.4. Genotyping and Microsatellite Analyses

A total of eleven polymorphic microsatellite loci were selected [26]. The corresponding primers were used to genotype 279 specimens of M. apiaster. The thermal cycling was done under the following conditions: 95 °C for 5 min, 29 cycles of 95 °C denaturation for 30 s, 51–61 °C annealing temperature for 1.5 min, 72 °C for 30 s of extension and 4 °C for 15 min. We used the Type-it PCR Kit for Microsatellites (Qiagen, Cat. No. 206243, Hilden, Germany), see Reference [3] for the reaction details.
PCR products were analyzed on an ABI gel electrophoresis sequencer conducted by GATC Biotech AG (Konstanz, Germany). The software PEAK ScannerTM version 2.0 was used to determine the allele sizes. The alleles were arranged by size in classes and plotted in the most probable class of allele size. The presence of null alleles and genotyping mistakes (amplification mistakes, recording of stutter peaks and short allele dominance) were verified using the software Microchecker version 2.2.3 (University of Hull, Hull, United Kingdom) [27]. Departures from the Hardy–Weinberg proportion and linkage equilibrium were verified for all loci using Genepop version 4.2 (Laboratiore de Genetique et Environment, Montpellier, France) [28], with significance estimated by Fisher exact, running the Markov chain algorithm with dememorization of 10,000.
The following parameters were calculated for all populations in Genetix version 4.05.2 (Université de Montpellier II, Montpellier, France) [29]: number of alleles (Na), allele richness (AR), observed heterozygosity (Ho), expected heterozygosity (He), gene diversity (HS), total gene diversity (HT). Fstat version 2.9.3 (UNIL, Lausanne, Switzerland) [30] was employed to determine the genetic polymorphism for each microsatellite locus (Na—number of alleles, Ar—average allelic richness, Ho—observed heterozygosity, Hs—gene diversity, HT—overall gene diversity, and Fst—genetic differentiation between populations). The levels of genetic differentiation between populations were determined using the Fst test [31] and the analysis of molecular variance (AMOVA) was conducted in Arlequin version 3.1 (Institute of Ecology and Evolution, University of Bern, Switzerland) [32]. The discriminant analysis of principal components was performed using R package adegenet; this approach uses the k-mean clustering algorithm according to the Bayesian Information Criterion to detect genetic clusters in the population [33].
The genetic structure was tested with a Bayesian clustering method and implemented in Structure version 2.3.4 using the admixture model without prior information about the populations [34]. MCMC parameters were fixed at 50,000 and 20,000, with run lengths of 106 iterations and K varying from 1 to 18 clusters. The calculations were repeated five times to confirm convergence among estimated settings. The most probable number of clusters (K) was based on the harmonic mean estimator [34] and on the second order rate of variation of the likelihood function referring to KK) [35].
The program Bottleneck [36] was used to test the heterozygosity excess in each population to identify population bottlenecks. We used the two-phase model of microsatellite evolution [37]. By means of GeneClass version 2.0 [38], we determined the direction of recent gene flow of M. apiaster between Europe and Africa by setting the parameters proposed by Paetkau et al. 2004 [39].

3. Results

3.1. Sequence Data

In total, we sequenced the mtDNA Cyt b and COI of 175 samples. We included samples of M. apiaster available in Genbank as well. Overall, we analyzed Cyt b of 299 samples. We found 48 polymorphic sites and 34 parsimony informative sites in Cyt b. The dataset shows substantial genetic variability. We identified 56 haplotypes, 22 of them with more than one sample. Two dominant haplotypes came from several localities in Europe, Asia, and Southern Africa: H1 with 77 samples from Greece, Spain, Macedonia, Germany, Turkey, Southern Africa and China; and H2 with 57 samples from Greece, Spain, Macedonia, Germany, Slovakia, Tunisia, Southern Africa and China (Figure 1, Table 1).
The overall haplotype network based on Cyt b yielded two main haplotypes (Figure 1) and was congruent with the dataset of COI (Table S2). Various unique haplotypes of low frequency derived from them and formed a star-like structure. The unique haplotypes were distinguished from the main haplotype by only one or two mutation steps. Overall, bee-eaters from all localities shared both main common haplotypes (H1 and H2), with the exception of bee-eaters from Turkey, which showed only haplotype H1, and those from Slovakia, Tunisia and Georgia, which had only haplotype H2 (Figure 1, Table 1). It is remarkable, that bee-eaters, which breed in Southern Africa, did not reveal unique haplotypes, but shared haplotypes with European populations, indicating that they probably derived from them.
Additionally, the haplotype network structure (Figure 1) revealed recent expansion events, which have been confirmed by the negative results of the demographic population parameters, test Fu’s Fs and Tajima’s D (Table 1). Data on genetic diversity, quantified by haplotype diversity (Hd) and nucleotide diversity (Pi), is summarized in Table 2. The overall values demonstrated high haplotype diversity 0.878 (Hd) and low nucleotide diversity 0.0043 (Pi). To test the hypothesis of population expansion in M. apiaster, the distribution of pairwise differences was calculated using the DnaSP version 5.0. Unimodal mismatch distribution and the goodness-of-fit of the observed data (Raggedness statistic r) revealed a model of sudden population expansion, subsequently to the last glacial bottleneck event (Figure 2). This finding was also supported by the Raggedness statistic with r = 0.0722 between observed and simulated data.

3.2. Microsatellite Data

Out of the eleven microsatellite loci, six loci produced well-defined alleles, which were in Hardy-Weinberg and linkage equilibrium (BE129, BE72, BE192, BE216, BE233, and BE231). Therefore, we restricted the following analysis to these six microsatellite loci. Additionally, we did not find evidence for large allele dropout, null alleles and stuttering signals using Microchecker 2.2.3. The diversity estimates differed between the microsatellite loci (Table 3) and between the populations analyzed (Table 4). The total number of alleles per locus ranged from 7 to 17 alleles and the average allelic richness (Ar) varied from 1.477 to 1.874. The average gene diversity (Hs) of all loci was 0.717 (range 0.456–0.867), whereas the overall gene diversity (HT) range was 0.458–0.885 (Table 3). Additionally, the lowest observed heterozygosity (Ho) and expected heterozygosity (He) were found for the population of Northern Africa, respectively 0.472 and 0.495, whereas the highest scores (Ho = 0.750) were observed in Scandinavia and Northern Asia (Table 4). The pairwise differences Fst varied from 0 to 3.6% (Table 5). Interestingly, the Fst values from the population of Southern Africa were the lowest. The AMOVA results demonstrated that most of the molecular variance was distributed within populations (92%) rather than among groups (2.5%) or among populations (1.2%) (Table 6).
Furthermore, the PCoA analysis grouped the populations into two clusters, however, samples from the distinct groups appear in all clades, evidencing a lack of genetic structure and indicating panmixia (Figure 3).
The estimated number of clusters obtained through the Bayesian cluster analysis, in a structure using the ΔK method, was six clusters. Overall, the different sampled bee-eater populations share genetic information; however, a slight differentiation was perceived in the populations from Western Asia, East Europe, Scandinavia and Southern Africa with K = 5 (Figure 4). No signal of a potential bottleneck was detected, based on the L shape of the distribution of allele frequencies observed in all the populations (Table S3).
The direction of the population movements based on the gene flow among the populations and estimated throughout the software GeneClass version 2.0 [38], detected the following routes: RLP/BW (western Germany) to Northern Asia, Saxony-Anhalt (eastern Germany) to West Asia and to RLP/BW, Northern Africa to East Europe, Southern Africa to Greece, Western Asia to East Europe and Greece to Western Asia.

3.3. Divergence Time Estimation

Molecular clock estimation, based on the evolutionary rate of Cyt b, revealed splitting events within the Meropidae family taking place during the Miocene/Pliocene/Pleistocene periods. The speciation of M. apiaster apparently took place in the mid-Pliocene with a highly supported sister relationship with M. malimbicus. The population establishment of the European bee-eater probably occurred during the late Pleistocene, about 0.6 Million years ago (Mya) (Figure 5).

4. Discussion

4.1. Historical Demography of the European Bee-eater

Our data reveal a weak phylogeographic structure between the populations of European bee-eaters, which suggests high gene flow between populations from different localities, resulting in a genetic panmixia. These results are in agreement with previous DNA studies from other migratory birds of the Western Palearctic [2,3,40], which revealed low phylogeographic structures. The lack of genetic structure observed for the populations of M. apiaster is likely to be explained by shared ancestral polymorphism due to the recent divergence in this group.
The low genetic differentiation of Palearctic birds can be explained by the frequent climatic oscillation during the Pleistocene, which repeatedly forced the birds to retreat to more suitable regions in the Mediterranean, subtropical and tropical Africa and Asia during glacial periods [8]. During warm periods, populations moved north again and differentiated locally. Migratory birds can present low genetic structure because they share reproductive sites; consequently, a reduction in the amount of exclusive of alleles within the populations occur. Thereby, populations from different localities in Europe likely mix with each other during the winter season and colonize other points after they return to Europe, which has probably resulted in the observed gene flow between birds from different populations [41]. It is also possible, that the lineages mixed earlier when they retracted to refugia.
Furthermore, some bird species have the ability to change the migratory pathway because of social or local environmental circumstances [42]. This plasticity was necessary for the adaptation of birds in different regions and different climatic conditions [43], however it has not yet been recognized in European bee-eaters.
The outcome of the demographic tests of Tajima’s D, Fu’s Fs, and Fu and Li’s F can explain the demographic forces that have defined the historical population dynamics of European bee-eaters. The negative values herein accessed reflect an excess of rare polymorphisms in all populations of M. apiaster, which is linked to population expansion and positive selection. Indeed, historical signs of population expansion northwards toward central Europe and Great Britain have been recorded for European bee-eaters, which were related to global warming and the appearance of new suitable habitats [13].
The slight genetic differentiation observed in the populations from Eastern Europe and Western Asia (Figure 4) might be associated with the distinct migration routes used by them. A similar pattern of the slight distinction between eastern and Western Europe was previously observed from microsatellite data [14]. European bee-eaters use two migration routes to Africa: Bee-eaters breeding in Western Europe move across the Straits of Gibraltar to overwintering areas in Western Africa, between Senegal and Nigeria, and Eastern European populations travel down to Southern Africa past the Eastern Mediterranean Sea and the Sahara Desert. These two migratory routes seem to drive the populations to differentiate genetically into two distinct groups [14].

4.2. Lack of Differentiation of the African Populations from Merops apiaster

The population, which breeds in the western half of Southern Africa and Southern Namibia and spends the non-breeding season in Malawi, was expected to show unique alleles and to form a distinct cluster in the structure analyses. Nonetheless, they in fact share haplotypes and belong to the same genetic pool as the Eurasian migratory populations. As the African populations share the same haplotypes and microsatellite variation with Eurasian populations, they might have derived from migrants which remained in Southern Africa for breeding.

4.3. Molecular Dating of the European Bee-eater

The Neoaves appeared in the Cretaceous but rapidly radiated in the early Tertiary, with intensive speciation processes during the Eocene. The origin of the Coraciiformes, to which the Meropidae belong, is dated to the Late Eocene circa 55 Mya [43]. Contributions of Woodruff 2010 [44] on the biogeography of Southeast Asia highlighted that landscape modifications of this region have resulted in multiple sea-level oscillations, successive biotic compressions, and expansions. Intensive diversification of taxa in the genus Merops was dated to the Pliocene and Pleistocene periods. Many bird species of the rainforests have also differentiated since the late Miocene and mid-Pliocene [45].
The evolution of the genus Merops partially agrees with the speciation process linked to dispersal routes for passerine birds proposed by Johansson et al., 2008 [46]. The ancestry of this family might have started in Asia from where bee-eaters dispersed in two directions during the Miocene: One route to Africa may have been through Asia and the second route over Madagascar to continental Africa.
In conclusion, molecular time estimations indicate that the radiation of bee-eaters took place mainly during the Pliocene. Populations of M. apiaster evolved about 0.6 Mya, during a period, which is characterized by climatic oscillations between warm and cold periods in the northern hemisphere. Although European bee-eaters breed from Western Europe to the Far East, they do not show phylogeographic differentiation. Haplotypes show a panmictic pattern, which probably resulted from population mixings in the winter quarters and in refugia during the last glaciation periods. European bee-eaters have spread and now breed in Southern Africa, but they remain migrants.

Supplementary Materials

The following are available online at https://www.mdpi.com/1424-2818/11/1/12/s1, Table S1: Sampling localities of Merops apiaster specimens used in this study and including also the other species belonging to the family Meropidae. GenBank accession numbers MH217684 - MH218120, Table S2: Genetic diversity and demographic population parameters based on the Cytochrome Oxidase I gene of different populations of Merops apiaster. Significant values p < 0.05 are displayed in bold, Table S3: Wilcoxon sign rank tests for heterozygosity excess in six populations of Merops apiaster based on six microsatellite loci and the allele frequency distribution. Three mutation models were tested: IAM-Infinite alleles model; SMM- Stepwise mutation model and TPM- Two-phase mutation model.

Author Contributions

C.C.d.M.M. and M.W. developed the ideas. C.C.d.M.M. conducted the laboratory work, bioinformatic analysis and wrote the manuscript. M.W. provided the laboratory supplies and contributed with important input for the improvement of the manuscript. A.B. and H.T.B. provided samples and reviewed the manuscript. E.W. and X.W. helped with laboratory support and reviewed the manuscript.

Funding

This research was partially funded by CAPES/Ciências sem fronteiras 1104-13/06. We acknowledge financial support by Deutsche Forschungsgemeinschaft within the funding programme Open Access Publishing, by the Baden-Württemberg Ministry of Science, Research and the Arts and by Ruprecht-Karls-Universität Heidelberg.

Acknowledgments

CCMM was supported by a PhD scholarship from CAPES/Ciências sem fronteiras 1104-13/06. We express gratitude to the staff and curators of the following bird collections for allowing us to use the tissues under their care for this study: Claire Spottiswood; FMNH, Chicago, IL, USA; LSUMNS, Baton Rouge, LA, USA; ZMUC, Copenhagen, Denmark; SMNH, Stockholm, Sweden; Zoological Museum of Bashkir University, Russia (ZMBU); Zoological Museum of Belarusian University, Azerbaijan. The laboratory procedures were conducted at the Institute of Pharmacy and Molecular Biotechnology, Heidelberg University. We also acknowledge Hedwig Sauer-Gürth for laboratory support and Alisa Samoylova for helping with laboratory support during the development of this project.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Nittinger, F.; Haring, E.; Pinsker, W.; Wink, M.; Gamauf, A. Out of Africa? Phylogenetic relationships between Falco biarmicus and the other hierofalcons (Aves: Falconidae). J. Zool. Syst. Evol. Res. 2005, 43, 321–331. [Google Scholar] [CrossRef]
  2. Shephard, J.M.; Ogden, R.; Tryjanowski, P.; Olsson, O.; Galbusera, P. Is population structure in the European white stork determined by flyway permeability rather than translocation history? Ecol. Evol. 2013, 3, 4881–4895. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Wang, E.; Van Wijk, R.E.; Santhosh, M.; Wink, M. Gene flow and genetic drift contribute to high genetic diversity with low phylogeographical structure in European hoopoes (Upupa epops). Mol. Phylogenet. Evol. 2017, 113, 113–125. [Google Scholar] [CrossRef] [PubMed]
  4. Cresswell, W. Migratory connectivity of Palaearctic-African migratory birds and their responses to environmental change: The serial residency hypothesis. Ibis 2014, 156, 493–510. [Google Scholar] [CrossRef]
  5. Irwin, D.E.; Irwin, J.H.; Smith, T.B. Genetic variation and seasonal migratory connectivity in Wilson’s warblers (Wilsonia pusilla): Species-level differences in nuclear DNA between western and eastern populations: Migratory connectivity and genetic divergence. Mol. Ecol. 2011, 20, 3102–3115. [Google Scholar] [CrossRef] [PubMed]
  6. Hewitt, G. The genetic legacy of the Quaternary ice ages. Nature 2000, 405, 907–913. [Google Scholar] [CrossRef] [PubMed]
  7. Newton, I.; Brockie, K. The Migration Ecology of Birds; Elsevier/Academic Press: Amsterdam, The Netherlands, 2008; ISBN 978-0-12-517367-4. [Google Scholar]
  8. Moreau, R.E. The main vicissitudes of the european avifauna since the Pliocene. Ibis 2008, 96, 411–431. [Google Scholar] [CrossRef]
  9. Fry, C.H. Family Meropidae (bee-eaters). In Handbook of the Birds of the World; Lynx Edicions: Barcelona, Spain, 2001; Volume 6, pp. 286–341. ISBN 84-87334-30-X. [Google Scholar]
  10. Marks, B.D.; Weckstein, J.D.; Moyle, R.G. Molecular phylogenetics of the bee-eaters (Aves: Meropidae) based on nuclear and mitochondrial DNA sequence data. Mol. Phylogenet. Evol. 2007, 45, 23–32. [Google Scholar] [CrossRef]
  11. Burt, B.D. Plumage-based phylogenetic analyses of the Merops. Mol. Phylogenet. Evol. 2004, 146, 481–492. [Google Scholar]
  12. Howard, R.; Moore, A.; Dickinson, E.C. The Howard and Moore Complete Checklist of the Birds of the World, Enlarged 3rd ed.; Princeton University Press: Princeton, NJ, USA, 2003; ISBN 978-0-691-11701-0. [Google Scholar]
  13. Fraser, P.A.; Rogers, M.J.; Calandrella, S.L. Report on scarce migrant birds in Britain in 2001. Br. Birds 2003, 96, 626–649. [Google Scholar]
  14. Ramos, R.; Song, G.; Navarro, J.; Zhang, R.; Symes, C.T.; Forero, M.G. Population genetic structure and long-distance dispersal of a recently expanding migratory bird. Mol. Phylogenet. Evol. 2016, 99, 194–203. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Sambrook, J.; Fritsch, E.F.; Maniatis, T. Molecular Cloning: A Laboratory Manual; Cold Spring Harbor Laboratory Press: Cold Spring Harbor, NY, USA, 1989; ISBN-13 978-0879693091. [Google Scholar]
  16. Sorenson, M.D.; Ast, J.C.; Dimcheff, D.E.; Yuri, T.; Mindell, D.P. Primers for a PCR-Based Approach to Mitochondrial Genome Sequencing in Birds and Other Vertebrates. Mol. Phylogenet. Evol. 1999, 12, 105–114. [Google Scholar] [CrossRef] [PubMed]
  17. Brumfield, R.T.; Edwards, S.V. Evolution into and out of the Andes: A Bayesian analysis of historical diversification in Thamnophilus Antshrikes. Evolution 2007, 61, 346–367. [Google Scholar] [CrossRef]
  18. Librado, P.; Rozas, J. DnaSP v5: A software for comprehensive analysis of DNA polymorphism data. Bioinformatics 2009, 25, 1451–1452. [Google Scholar] [CrossRef] [PubMed]
  19. Rogers, A.H.; Harpending, H. Population growth makes waves in the distribution of pairwise genetic differences. Mol. Biol. Evol. 1992, 9, 552–569. [Google Scholar] [CrossRef]
  20. Harpending, H.C. Signature of Ancient Population Growth in a Low-Resolution Mitochondrial DNA Mismatch Distribution. Hum. Biol. 1994, 66, 591–600. [Google Scholar]
  21. Posada, D. jModelTest: Phylogenetic Model Averaging. Mol. Biol. Evol. 2008, 25, 1253–1256. [Google Scholar] [CrossRef] [Green Version]
  22. Weir, J.T.; Schluter, D. Calibrating the avian molecular clock: Avian Molecular Clock. Mol. Ecol. 2008, 17, 2321–2328. [Google Scholar] [CrossRef]
  23. Drummond, A.J.; Suchard, M.A.; Xie, D.; Rambaut, A. Bayesian Phylogenetics with BEAUti and the BEAST 1.7. Mol. Biol. Evol. 2012, 29, 1969–1973. [Google Scholar] [CrossRef] [Green Version]
  24. Rambaut, A.; Drummond, A.J. Tracer v. 1.5; Institute of Evolutionary Biology, University of Edinburgh: Edinburgh, UK, 2014. [Google Scholar]
  25. Miller, M.A.; Pfeiffer, W.; Schwartz, T. Creating the CIPRES Science Gateway for inference of large phylogenetic trees. In Proceedings of the 2010 Gateway Computing Environments Workshop (GCE), New Orleans, LA, USA, 14 November 2010. [Google Scholar]
  26. Dasmahapatra, K.K.; Lessels, C.M.; Mateman, A.C.; Amos, W. Microsatellite loci in the European bee-eater, Merops apiaster. Mol. Ecol. 2004, 500–502. [Google Scholar] [CrossRef]
  27. Van Oosterhout, C.; Hutchinson, W.F.; Wills, D.P.M.; Shipley, P. micro-checker: Software for identifying and correcting genotyping errors in microsatellite data. Mol. Ecol. Notes 2004, 4, 535–538. [Google Scholar] [CrossRef]
  28. Rousset, F. Genepop’007: A complete re-implementation of the genepop software for Windows and Linux. Mol. Ecol. Resour. 2008, 8, 103–106. [Google Scholar] [CrossRef]
  29. Belkhir, K. GENETIX (v4. 04): A Windows Program for Population Genetic Analysis; Laboratorie Genome, Populations Interactions, UPR 5000 du CNRS; Universite Montpellier II: Montpellier, France, 1999. [Google Scholar]
  30. Goudet, J. FSTAT, A Program to Estimate and Test Gene Diversity and Fixation Indices (Version 2.9.3.2); Lausanne University: Lausanne, Switzerland, 2001. [Google Scholar]
  31. Weir, B.S.; Cockerham, C.C. Estimating F-statistics for the analysis of population structure. Evolution 1984, 38, 1358–1370. [Google Scholar] [CrossRef]
  32. Excoffier, L.; Laval, G.; Schneider, S. Arlequin (version 3.0): An integrated software package for population genetics data analysis. Evol. Bioinform. 2005, 1, 47–50. [Google Scholar] [CrossRef]
  33. Jombart, T. adegenet: A R package for the multivariate analysis of genetic markers. Bioinformatics 2008, 24, 1403–1405. [Google Scholar] [CrossRef]
  34. Pritchard, J.K.; Stephens, M.; Donnelly, P. Inference of population structure using multilocus genotype data. Genetics 2000, 155, 945–959. [Google Scholar]
  35. Evanno, G.; Regnaut, S.; Goudet, J. Detecting the number of clusters of individuals using the software structure: A simulation study. Mol. Ecol. 2005, 14, 2611–2620. [Google Scholar] [CrossRef]
  36. Piry, S.; Luikart, G.; Cornuet, J.M. BOTTLENECK: A computer program for detecting recent reductions in the effective population size using allele frequency data. J. Hered. 1999, 90, 502–503. [Google Scholar] [CrossRef]
  37. Di Rienzo, A.; Peterson, A.C.; Garza, J.C.; Valdes, A.M.; Slatkin, M.; Freimer, N.B. Mutational processes of simple-sequence repeat loci in human populations. Proc. Natl. Acad. Sci. USA 1994, 91, 3166–3170. [Google Scholar] [CrossRef]
  38. Piry, S.; Alapetite, A.; Cornuet, J.-M.; Paetkau, D.; Baudouin, L.; Estoup, A. GENECLASS2: A Software for Genetic Assignment and First-Generation Migrant Detection. J. Hered. 2004, 95, 536–539. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Paetkau, D.; Slade, R.; Burden, M.; Estoup, A. Genetic assignment methods for the direct, real-time estimation of migration rate: A simulation-based exploration of accuracy and power. Mol. Ecol. 2004, 13, 55–65. [Google Scholar] [CrossRef] [PubMed]
  40. Arbabi, T.; Gonzalez, J.; Witt, H.-H.; Klein, R.; Wink, M. Mitochondrial phylogeography of the Eurasian Reed Warbler Acrocephalus scirpaceus and the first genetic record of A. s. fuscus in Central Europe. Ibis 2014, 156, 799–811. [Google Scholar] [CrossRef]
  41. Dhanjal-Adams, K.L.; Bauer, S.; Emmenegger, T.; Hahn, S.; Lisovski, S.; Liechti, F. Spatiotemporal group dynamics in a long-distance migratory bird. Curr. Biol. 2018, 28, 2824–2830. [Google Scholar] [CrossRef] [PubMed]
  42. Guillemain, M.; Sadoul, N.; Simon, G. European flyway permeability and abmigration in Teal Anas crecca, an analysis based on ringing recoveries. Ibis 2005, 147, 688–696. [Google Scholar] [CrossRef]
  43. Ericson, P.G.; Anderson, C.L.; Britton, T.; Elzanowski, A.; Johansson, U.S.; Kallersjo, M.; Ohlson, J.I.; Parsons, T.J.; Zuccon, D.; Mayr, G. Diversification of Neoaves: Integration of molecular sequence data and fossils. Biol. Lett. 2006, 2, 543–547. [Google Scholar] [CrossRef] [PubMed]
  44. Woodruff, D.S. Biogeography and conservation in Southeast Asia: How 2.7 million years of repeated environmental fluctuations affect today’s patterns and the future of the remaining refugial-phase biodiversity. Biodivers. Conserv. 2010, 19, 919–941. [Google Scholar] [CrossRef]
  45. Fjeldså, J.; Bowie, R.C.K. New perspectives on the origin and diversification of Africa’s forest avifauna. Afr. J. Ecol. 2008, 46, 235–247. [Google Scholar] [CrossRef]
  46. Johansson, U.S.; Fjeldså, J.; Bowie, R.C.K. Phylogenetic relationships within Passerida (Aves: Passeriformes): A review and a new molecular phylogeny based on three nuclear intron markers. Mol. Phylogenet. Evol. 2008, 48, 858–876. [Google Scholar] [CrossRef]
Figure 1. Haplotype network of Merops apiaster using statistical parsimony based on the cytochrome b gene. Colors correspond to the specific origins.
Figure 1. Haplotype network of Merops apiaster using statistical parsimony based on the cytochrome b gene. Colors correspond to the specific origins.
Diversity 11 00012 g001
Figure 2. Mismatch distribution of frequencies of the pairwise number of differences among Merops apiaster populations. The dotted line represents the observed values in a stable population with a constant population size. The solid line represents the observed frequency of the pairwise mismatch distribution of mtDNA.
Figure 2. Mismatch distribution of frequencies of the pairwise number of differences among Merops apiaster populations. The dotted line represents the observed values in a stable population with a constant population size. The solid line represents the observed frequency of the pairwise mismatch distribution of mtDNA.
Diversity 11 00012 g002
Figure 3. Discriminant analysis of principal components based on six microsatellite loci of Merops apiaster. The dots correspond to the individuals while the ellipses represent the clusters. Colors represent each different cluster. RLP/BW—Rhineland-Palatinate, Baden-Württemberg and Saxony-Anhalt (Germany).
Figure 3. Discriminant analysis of principal components based on six microsatellite loci of Merops apiaster. The dots correspond to the individuals while the ellipses represent the clusters. Colors represent each different cluster. RLP/BW—Rhineland-Palatinate, Baden-Württemberg and Saxony-Anhalt (Germany).
Diversity 11 00012 g003
Figure 4. Graphic representation of Bayesian cluster analysis conducted in a structure based on six microsatellite loci performed with a K from 2 to 6. Population names are provided under the graphics. Each individual is represented by one vertical line and the colors point toward the probability of belonging to each of the genetic clusters.
Figure 4. Graphic representation of Bayesian cluster analysis conducted in a structure based on six microsatellite loci performed with a K from 2 to 6. Population names are provided under the graphics. Each individual is represented by one vertical line and the colors point toward the probability of belonging to each of the genetic clusters.
Diversity 11 00012 g004
Figure 5. Phylogenetic relationships and molecular dating of Merops apiaster inferred throughout Bayesian analysis using Cyt b sequences (in total 1059 bp). Posterior probabilities are placed at the nodes.
Figure 5. Phylogenetic relationships and molecular dating of Merops apiaster inferred throughout Bayesian analysis using Cyt b sequences (in total 1059 bp). Posterior probabilities are placed at the nodes.
Diversity 11 00012 g005
Table 1. Genetic diversity and demographic population parameters based on the cytochrome b gene of different populations of Merops apiaster. * p < 0.05.
Table 1. Genetic diversity and demographic population parameters based on the cytochrome b gene of different populations of Merops apiaster. * p < 0.05.
RegionLocalityNHsHdPiFu and Li’s FFu’s FsTajima’s D
Western EuropeRLP/BW149130.8790.0041−2.427 *−4.417 *−1.857 *
Saxony-Anhalt6717230.8860.0037−2.868 *−6.179 *−1.472
Greece3923250.9450.0036−2.249−18.44 *−1.751
Spain9632380.8860.0027−3.919 *−29.18 *−2.189 *
Eastern EuropeSlovakia2212150.8700.0045−1.866−5.870 *−1.473
Macedonia2014180.9580.0028−2.730 *−9.967 *−1.916 *
Western AsiaWestern Asia5350.7000.0030−0.5770.804−0.561
Central AsiaChina3216860.9270.0087−4.7645 *−2.060−2.525 *
Southern AfricaSouth Africa7550.9050.0032−0.392−1.654−0.330
All29656860.8780.0043−6.800 *−64.19 *−2.524 *
N—Number of samples, H—Number of Haplotypes, s—Number of polymorphic (segregating) sites, Hd—Haplotype diversity, Pi—Nucleotide diversity. RLP/BW—Rhineland-Palatinate, Baden-Württemberg (Germany).
Table 2. Haplotype table based on cytochrome b of Merops apiaster. Single haplotypes were not included. N = haplotype frequency.
Table 2. Haplotype table based on cytochrome b of Merops apiaster. Single haplotypes were not included. N = haplotype frequency.
HapVariable SitesNLocalities
H1TTTTTCCACTGGAATCTCTGGCTTCATCAACTGCCTTTGACC64China, Greece, Spain, Macedonia, Saxony-Anhalt, Baden-Württemberg, Rhineland-Palatinate, Turkey, Southern Africa.
H2.....T....................................32China, Greece, Spain, Slovakia, Macedonia, Saxony-Anhalt, Egypt, Tunisia, Baden-Württemberg, Rhineland-Palatinate.
H3.....T.........T..........................27Saxony-Anhalt, Spain, Macedonia, BW/RLP, Slovakia, Southern Africa.
H4....CT....................................17China, Greece, Spain, Saxony-Anhalt, Rhineland-Palatinate.
H5.....T........C...........................13China, Spain, Saxony-Anhalt, Greece.
H6.....T...........................A........14China, Greece, Spain, Southern Africa, Macedonia, Rhineland-Palatinate.
H7.....T.................C......T...........7China, Greece, Macedonia.
H8..................C.......................7Spain, Greece and Macedonia.
H9.....T.G..................................6Slovakia.
H10.....T..............................C.....5China and Greece.
H11..C..T....................................4Saxony-Anhalt.
H12.........................C................4Kazakhstan, Rhineland-Palatinate, Greece.
H13.....T.....................T........C.....4Spain and China.
H14C.........................................3Southern Africa.
H15.....T..............A.....................3Spain.
H16.....T..T...............T........A........3China and Spain.
H17..................................T.......3Macedonia and Greece.
H18...C.T....................................2China.
H19.....T.......................G............2Spain and Slovakia.
H20.....T........C................C..........2Spain and Greece.
H21C.........................C...............2Saxony-Anhalt and Macedonia.
H22.....T.........T............G.............2Saxony-Anhalt and Slovakia.
Table 3. Genetic diversity at six microsatellite loci in European bee-eaters.
Table 3. Genetic diversity at six microsatellite loci in European bee-eaters.
LocusNaArHoHsHTFST
BE129131.8740.7660.8550.8760.020
BE72171.8050.9040.8460.8620.030
BE19291.8290.6780.7840.8350.050
BE21671.4770.4410.4960.5130.024
BE23381.5220.4020.4560.4580.008
BE231161.8560.7930.8670.8850.041
Na—the observed number of alleles, Ar—average allelic richness, Ho—observed heterozygosity, Hs—gene diversity, HT—overall gene diversity, Fst—genetic differentiation between populations (estimated based on Weir and Cockerham 1984).
Table 4. Population size (N), the average number of alleles per locus (Na), observed heterozygosity (Ho) and expected heterozygosity (He) at six microsatellite loci in the European bee-eater.
Table 4. Population size (N), the average number of alleles per locus (Na), observed heterozygosity (Ho) and expected heterozygosity (He) at six microsatellite loci in the European bee-eater.
PopulationNNaHoHe
Spain42.5000.6660.479
RLP/BW (Germany)218.6660.6870.730
Saxony-Anhalt (Germany)15310.330.6990.716
Greece146.8330.6580.688
Scandinavia22.8330.7500.541
Eastern Europe327.1660.6240.599
Western Asia164.0900.6160.691
Northern Africa32.8330.4720.495
Northern Asia44.1660.7500.692
Sub-Saharan Africa45.8330.7140.719
RLP/BW—Rhineland-Palatinate and Baden-Württemberg (Germany), Scandinavia—Denmark and Sweden, Northern Africa—Egypt and Tunisia, Western Asia- Turkey, Georgia and Kuwait, Northern Asia—Russia and Kazakhstan, Sub-Saharan Africa—Southern Africa and Malawi.
Table 5. Pairwise Fst comparison matrix based on the microsatellite data of the European bee-eater.
Table 5. Pairwise Fst comparison matrix based on the microsatellite data of the European bee-eater.
SpainRLP/BWSaxonyGreeceScandinaviaEastern EuropeWestern AsiaNorthern AfricaNorthern AsiaSub-Saharan Africa
Spain-
RLP/BW0.0053-
Saxony-Anhalt0.00960.0003-
Greece0.00710.00290.0121-
Scandinavia0.00000.00850.00960.0007-
Eastern Europe0.02750.02280.02720.03050.0353-
Western Asia0.00740.00770.01030.02400.00740.0243-
Northern Africa−0.0514−0.0457−0.0141−0.0580.0366−0.1407−0.0630
Northern Asia0.0000−0.00530.0039−0.0580.00000.0274−0.0012−0.0102
Sub-Saharan Africa−0.0181−0.0010−0.00010.00290.00000.02110.0063−0.0574−0.0090-
RLP/BW—Rhineland-Palatinate and Baden-Württemberg (Germany), Scandinavia—Denmark and Sweden, Western Asia- Turkey, Georgia and Kuwait, Northern Africa—Egypt and Tunisia, Northern Asia—Russia and Kazakhstan, Sub-Saharan Africa—Southern Africa and Malawi. In bold p < 0.05.
Table 6. Analysis of molecular variance (AMOVA) for six microsatellite loci to test the variation among the following breeding populations of the bee-eaters: Group 1 (Spain, Germany (RLP, BW, Saxony-Anhalt), Greece); Group 2 (Denmark, Sweden); Group 3 (Slovakia); Group 4 (Turkey, Georgia, Kuwait); Group 5 (Egypt, Tunisia); Group 6 (Kazakhstan, Russia) and Group 7 (Southern Africa, Malawi).
Table 6. Analysis of molecular variance (AMOVA) for six microsatellite loci to test the variation among the following breeding populations of the bee-eaters: Group 1 (Spain, Germany (RLP, BW, Saxony-Anhalt), Greece); Group 2 (Denmark, Sweden); Group 3 (Slovakia); Group 4 (Turkey, Georgia, Kuwait); Group 5 (Egypt, Tunisia); Group 6 (Kazakhstan, Russia) and Group 7 (Southern Africa, Malawi).
Source of VariationSum of SquaresPercentage of VariationFSTp ValueAverage F-Statistics Over All Loci
FISFSTFIT
Among groups31.0642.5220.010410.0000.043830.03050.0733
Among populations within groups11.0521.248
Among individuals within populations535.8194.21769
Within populations512.00092.01154
FIS—variance among individuals, FIT—variance within-individuals, and FS—variance among populations.

Share and Cite

MDPI and ACS Style

Carneiro de Melo Moura, C.; Bastian, H.-V.; Bastian, A.; Wang, E.; Wang, X.; Wink, M. Pliocene Origin, Ice Ages and Postglacial Population Expansion Have Influenced a Panmictic Phylogeography of the European Bee-Eater Merops apiaster. Diversity 2019, 11, 12. https://doi.org/10.3390/d11010012

AMA Style

Carneiro de Melo Moura C, Bastian H-V, Bastian A, Wang E, Wang X, Wink M. Pliocene Origin, Ice Ages and Postglacial Population Expansion Have Influenced a Panmictic Phylogeography of the European Bee-Eater Merops apiaster. Diversity. 2019; 11(1):12. https://doi.org/10.3390/d11010012

Chicago/Turabian Style

Carneiro de Melo Moura, Carina, Hans-Valentin Bastian, Anita Bastian, Erjia Wang, Xiaojuan Wang, and Michael Wink. 2019. "Pliocene Origin, Ice Ages and Postglacial Population Expansion Have Influenced a Panmictic Phylogeography of the European Bee-Eater Merops apiaster" Diversity 11, no. 1: 12. https://doi.org/10.3390/d11010012

APA Style

Carneiro de Melo Moura, C., Bastian, H. -V., Bastian, A., Wang, E., Wang, X., & Wink, M. (2019). Pliocene Origin, Ice Ages and Postglacial Population Expansion Have Influenced a Panmictic Phylogeography of the European Bee-Eater Merops apiaster. Diversity, 11(1), 12. https://doi.org/10.3390/d11010012

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop