High Genetic Diversity among Breeding Red-Backed Shrikes Lanius collurio in the Western Palearctic

Revealing the genetic population structure in abundant avian species is crucial for understanding speciation, conservation, and evolutionary history. The Red-backed Shrike Lanius collurio, an iconic songbird renowned for impaling its prey, is widely distributed as a breeder across much of Europe, Asia Minor and western Asia. However, in recent decades, many populations have declined significantly, as a result of habitat loss, hunting along migration routes, decrease of arthropod food, and climate change e.g., severe droughts in Africa. Within this context, gene flow among different breeding populations becomes critical to ensure the survival of the species, but we still lack an overview on the genetic diversity of the species. In this paper, we analyzed the mitochondrial cytochrome b gene (mtDNA) and the cytochrome c oxidase subunit 1 gene (mtDNA) of 132 breeding Red-backed Shrikes from across the entire breeding range to address this knowledge gap. Our results revealed consistent genetic diversity and 76 haplotypes among the Eurasian populations. Birds are clustered in two major groups, with no clear geographical separation, as a direct consequence of Pleistocene glaciations and apparent lineage mixing in refugia. This has led to genetic panmixia.


Introduction
Quaternary cold periods have greatly shaped the evolution and distribution of present biota worldwide [1], and this is particularly evident in the northern temperate regions [2].During the Pleistocene era-starting 2.5 million years ago and ending just 11,000 years before present-the climate in northern temperate regions was characterized by dramatic oscillations, commonly known as "glaciations" [3,4], each of them lasting anywhere between 21,000 to 100,000 years [5].These severe glacial episodes forced biota to retreat to southern refugia [6,7], which implied admixing of populations.Subsequently, after the retreat of the ice cap, populations expanded and recolonized the northern latitudes [8,9].This complex process of post-glacial population expansion northward from the southern refugia occurred at least 30 times during the Pleistocene [3] and was in strong interdependence with the vegetation availability [10].As the flora and fauna recolonized northern habitats, populations from different refugia came in contact and merged, sometimes creating the so-called hybrid or suture zones [8,11].Based on morphological features alone, the admixture of population lineages and suture zones were evident in numerous species across the Palearctic [12].But it was only after the emergence of molecular techniques and their wide application in zoology that the history of post-glaciation colonization and its impact on the genetic structure of European biota really unfolded [13].
Across the animal kingdom, birds have a particular significance for phylogeographic and phylogenetic studies due to their high locomotive capacities, which stimulates both gene flow and speciation events but also, in many taxa, fidelity to the natal site, which slows down population admixture [14].Furthermore, due to their ubiquitous status and charismatic nature, birds have been at the forefront of genetic studies in zoology in recent decades [15].
The Shrike family (Aves, Laniidae) is currently comprised of 34 species, belonging to four genera: Corvinella, Urosteles, Eurocephalus, and Lanius [16,17].Several species, notably the ones from western Palearctic, are well studied, partly due to their charisma and habit of impaling their prey.Nonetheless, within the true shrikes-genus Lanius, only few taxa have been the focus of genetic studies: the Lesser Grey Shrike (L.minor) [18], the Great Grey Shrike (L.excubitor) [19], the Southern Grey Shrike (L.meridionalis) [20,21], and the Loggerhead Shrike (L.ludovicianus) from North America [22,23].The genetic structure of the Red-backed Shrike (L.collurio), which, despite its population decline in recent decades [24] still remains the most common shrike species across Europe [25,26], has not been yet elucidated.Furthermore, the Red-backed Shrike (hereafter RBS) is illegally trapped in the Mediterranean basin in great numbers, with as many as 30,000 birds/year in the Sinai Peninsula alone [27,28].As such, gene flow becomes important in ensuring the wellbeing of the species [29].To approach this knowledge gap, we analyzed the mitochondrial DNA of 132 breeding RBSs from 14 countries across the breeding range.Mitochondrial DNA (hereafter mtDNA) is an efficient marker for indicating the rate of evolution and speciation, haplotype diversity, and natal origin [30,31].For our study, we selected the cytochrome oxidase B (hereafter Cytb) and cytochrome c oxidase subunit 1 (hereafter CO1) mitochondrial genes, regarded as two of the most efficient molecular markers for avian population genetics based on mtDNA [32].However, it should be noted that mtDNA is only maternally inherited and non-recombinant, and therefore reflects solely the evolutionary history of females.Moreover, a number of studies [33][34][35] have challenged the advantages of mtDNA usage in phylogeny, but we argue that our results, albeit derived from just 1.551 base pairs, are in agreement with similar European studies on avian phylogeography (see Discussion for more details).The objectives of our research were (i) to evaluate the genetic diversity and population structure, in conjunction with (ii) phylogeography and evolutionary history of the RBS breeding in the western Palearctic.Based on our findings, we discuss the factors that shaped the current genetic structure of RBS across its breeding range and future implications for its conservation.

Sampling
We collected blood (30 µL) from breeding birds (n = 17) in Germany (Rhein-Neckar area) by puncturing the brachial vein, during the 2016 to 2017 breeding seasons.We further received blood, feathers, tissue, or buccal swabs from RBSs breeding in the following countries: Belgium, Bulgaria, Czech Republic, France, Hungary, Latvia, The Netherlands, Poland, Romania, Russia, Spain, Sweden, and Ukraine.For our collaborators with no previous experience in collecting buccal swabs, we provided a video protocol on how to sample buccal mucosa from live birds (see Video S1).All samples obtained from live birds were collected by accredited bird ringers, following national regulations.Blood and buccal swabs were conserved in EDTA buffer (0.1 M Tris, pH 7.4, 10% EDTA, 1% NaF, 0.1% Thymol) or in ethanol (70%).Feathers were stored in dry envelopes.Samples from Sweden were received from the Swedish National History Museum.Samples from Russia (except 91037 to 91048) and Ukraine were obtained from the Zoological Museum of Moscow University.For each sample we allocated an internal voucher at the collection of the Institute for Pharmacy and Molecular Biotechnology (IPMB).Full details of all samples included in our study are found in Table S1.

Analysis of Mitochondrial Markers, Genetic Diversity, and Divergence Time
We assembled the raw sequences and visually inspected and corrected the alignments in CodonCode Aligner 5.01 (CodonCode Corporation, Centerville Massachusetts, USA).We trimmed all Cytb sequences to a fixed length of 909 base pairs and the CO1 to 642 base pairs.We used MEGA 7.0 [44] to check that our sequences contained no stop codons and that the open reading frame was identical for both Cytb and CO1 sequences.We concatenated these two genes in Seaview 4.7 [45] and performed all further analyses with the concatenated sequence.Upon verification, we employed DnaSP 6 [46] to check for nucleotide and haplotype diversity indexes and to group all samples from a certain country in a single population.Consequently, we analyzed the *arp file in Arlequin 3.5.2.2 [47] to obtain F ST values.Via Arlequin software we also corrected the haplotype diversity values for the sample size.
To accurately identify the correct position of the nucleotide polymorphic sites in our trimmed sequences, we used two reference genes.For Cytb, we used a L. collurio complete Cytb sequence available in GenBank, accession number EF621589.1.However, for CO1 we could not find a complete sequence; we therefore used a L. isabellinus complete mitochondrial genome sample, GenBank accession number KP995437.1,from which we used the CO1 as a reference.We employed PopArt (Population Analysis with Reticulate Trees) [48] to visualize the haplotype network and geographic distribution across the western Palearctic.
To estimate evolutionary divergence time, we used the concatenated aligned sequences of our 132 RBS samples, to which we added 5 shrike species, with an identical alignment i.e., Cytb and CO1.We first retrieved from GenBank the complete mitochondrial genomes of the 5 species, after which we selected the CytB and CO1 and, in Seaview 4.7 we concatenated these two genes.From GenBank, we retrieved information for the following shrike species: L. schach (NC_030604.1),L. isabellinus (KP995437.1),L. cristatus (NC_028333.1),L. tephronotus (NC_028333.1), and L. sphenocercus (KU884610.1).As an outgroup, we used two corvid species, namely Corvus corax (KX245148.1)and Corvus frugilegus (Y18522.2),for which the alignment methodology was the same as mentioned above.We conducted the analyses in MEGA 7.0, using the Kimura 2-parameter model and modelled the rate variation among sites with a gamma distribution (shape parameter = 1).The analysis involved 139 nucleotide sequences i.e., 132 RBS, the 5 other Lanius species and the 2 Corvus species; each sequence had 1551 positions.
Finally, we performed all molecular dating with the programs BEAST 1.8.0 and BEAUTi 1.8.0 [49], using the CIPRES cyberinfrastructure for phylogenetic research [50].To establish the divergence time, we applied the substitution rate of 2.1% substitutions/site/ million years for Cytb [51] and we estimated the rate for CO1.To select the most suitable evolutionary model, we used the jModeltest 2.1.10software [52].We further applied the HKY nucleotide substitution model, with a strict clock and Yule speciation process, a 10 6 chain length and sampling every 1000th tree.We then employed TreeAnnotator 1.8.0 [49] to summarize the sample of trees and FigTree 1.4.3 for final visualization [53].

Results
Our 1551-bp concatenated alignment revealed 76 haplotypes for the 132 RBS samples, with 0.96 haplotype diversity, 0.009 nucleotide diversity (π), and 110 segregating sites (Table 1).Taken separately, Cytb was characterized by 63 polymorphic sites and CO1 by 47 sites (Table S2).For the concatenated genes, we found the highest haplotype diversity (i.e., 1) in Bulgaria, France, Latvia and Romania (Table 1).The birds from Germany, Poland, Hungary, and the Czech Republic were characterized by a haplotype diversity of 0.974, 0.95, 0.934, and 0.857, respectively.The most common haplotype, haplotype 4, was represented by 23 individuals and it was spread across the entire breeding area (Figure 3 and Table S3).It was followed by haplotype 7, with a total of eight birds in Bulgaria, Romania, Spain, and Sweden.Haplotype 3 was represented by seven birds in Belgium, Germany, Latvia, and Spain.The next ten common haplotypes were distributed across 31 birds, with no apparent pattern.The remaining 63 haplotypes were singletons i.e., 63 RBS in our dataset were characterized by a unique haplotype.The pairwise F ST values (Table 2) indicated low genetic differentiation between the majority of the 14 populations in our dataset.A certain level of differentiation was observed between shrikes from Czech Republic and The Netherlands, as well as between Russia and The Netherlands, and Russia and Ukraine.The value for Russia and Ukraine was most curious, due to their geographical position as neighboring countries and the absence of ecological barriers.Yet, the distance between the two-sampling locations was approximately 1160 km, which might explain the relatively low sharing of genetic material, in spite of their mobility and low philopatry rates.The haplotype network (Figures 1 and 2) indicated a two-clade structure, with higher diversity in the larger clade.The two clades were separated by 26 mutation steps.The geographic distribution of haplotypes illustrates a higher number of singletons in France, Latvia, Hungary, and Romania (Figure 3 and Table S3).Our phylogenetic tree (Figure 4, Appendix A, and Figure S1) and estimate of evolutionary divergence (Table 3) indicate a separation time between the two major clades of approximately 1.2 million years ago.Between L. collurio and its sister species-L.isabellinus, the divergence event is estimated at 1.9 million years before present.In comparison with L. cristatus, L. tephronotus, and L. schach, the RBS separated approximately 4.2 million years ago.In our analysis, L. sphenocercus scoreed the highest differentiation, separating itself from the other shrikes 6 million years ago.Regarding the outgroup, the two Corvus species, the divergence is estimated to have happened 12.3 million years before present.
Our phylogenetic tree (Figure 4, Appendix A, and Figure S1) and estimate of evolutionary divergence (Table 3) indicate a separation time between the two major clades of approximately 1.2 million years ago.Between L. collurio and its sister species-L.isabellinus, the divergence event is estimated at 1.9 million years before present.In comparison with L. cristatus, L. tephronotus, and L. schach, the RBS separated approximately 4.2 million years ago.In our analysis, L. sphenocercus scoreed the highest differentiation, separating itself from the other shrikes 6 million years ago.Regarding the outgroup, the two Corvus species, the divergence is estimated to have happened 12.3 million years before present.

Behind the Mask: Genetic Panmixia
Our mtDNA analysis provides the most complete view to date on the genetic structure of the RBSs breeding in the western Palearctic.The high number (n = 76), diversity index (0.96), and random distribution of haplotypes characterizing our 132 birds from 14 countries indicate significant genetic diversity and low evidence for population differentiation.This situation is not unique to the RBS, but actually a common phylogeographic status for many bird species currently distributed in Eurasia.Genetic admixture has been previously shown to characterize the Hoopoe (Upupa epops) population [54], the White-winged Snowfinch (Montifringila nivalis) [55], the Common Sandpiper (Actitis hypoleucos) [56], the Marsh Warbler (Acrocephalus palustris) [57], the Paddyfield Warbler (Acrocephalus agricola) [58], the Eurasian Reed Warbler (Acrocephalus scirpaceus) [59], the Western Capercaillie (Tetrao urogallus) [60], the Hazel Grouse (Tetrastes bonasia) [61], the Eurasian Blackcap (Sylvia atricapilla) [62], the European Bee-eater (Merops apiaster) [63], and the Eurasian Jay (Garrulus glandarius) [64], to name just a few.For the RBS, this genetic background is a clear result of Pleistocene climatic history (see next chapter), with different lineages coming together into refugia and, possibly, even bird movements between different refugia.We further speculate that, during the interglacial periods, lineages from different refugia came in contact at the suture zones, in the maximum population expansion times.Currently, this genetic diversity is under no evident selection pressure, due to their migratory behavior, which hinders sedentarism and population differentiation.Moreover, the RBS is characterized by low philopatry rates i.e., <10% [65,66] (plus Luís Reino, Franck Hollander, Piotr Tryjanowski, Marcin Tobółka, and Boris Nikolov pers.com.) and a wide continuous breeding range i.e., from Portugal until central Siberia [24], which further promotes gene flow among breeding populations.For European shrike species, random haplotype distribution has also been found in the Lesser Grey Shrike (L.minor) [18,67].We also refer the readers to the Southern Grey Shrike (L.meridionalis koenigi) [20], which, although endemic to the Canary Islands and subject of a different population history (see next chapter), is likewise characterized by a high haplotype diversity i.e., 0.815 versus 0.96 in L. collurio (this study).Similarly, in non-European shrikes, the endemic island Loggerhead Shrike (L.ludovicianus anthonyi) indicates a random distribution of haplotypes across its restricted distribution in the northern California Chanel Islands [23].In Africa, the Fiscal Shrike (L.collaris) is characterized by a haplotype diversity of 0.932, as shown by a study which found 26 haplotypes among 45 individuals [68].
In contrast, sedentary bird species in Europe, especially those with fragmented distribution, tend to show a better-defined genetic structure.To illustrate this, we mention the Black Grouse (Lyrurus tetrix), which overall shows specific haplotypes for each subpopulation in the Austrian Alps [69].Regarding the distinct haplotype clades structure of the RBS, among European bird taxa, the occurrence of a two clade structure with no geographical distinction has been previously observed in the Common Redstart (Phoenicurus phoenicurus) [70] and the Eurasian Collared Dove (Streptopelia decaocto) [71].For the two haplotype groups characterizing the Common Redstart, Johnsen et al. mention a high maximum intraspecific distance (%) i.e., 5.08 [42], whereas for the RBS the distance is 2.76.Additionally, in Eurasia, two haplotype clade structures have been observed in the Middle Spotted Woodpecker (Dendrocoptes medius), but with a clear differentiation between the European and Asian haplotypes [72].A comparable situation has been indicated in the Little Owl (Athene noctua), with two distinctive clades for western and eastern European haplotypes [73].Apparent western and eastern European haplotype clades have also been suggested for the Great Reed Warbler (Acrocephalus arundinaceus) [74].For the Tawny Owl (Strix aluco), one genetic study has revealed the existence of three separate haplotype clades, one for each of the three main refugia (Iberian Peninsula, Italian Peninsula, and the Balkans) [75].Similarly, on the North American continent, two distinct mitochondrial clades have been detected in the Northern Raven (Corvus corax) [76] and the Snow Goose (Anser caerulescens caerulescens) [77], and in both species a certain degree of geographic structuring of the clades is reported.Most commonly, high mtDNA divergence results from cryptic speciation, hybridization of closely related species or separation in isolated refugia during the glacial times.For the RBS, we suggest that the current situation is a result of the latter, as the molecular divergence time of the two clades coincides with the Pleistocene period.Nonetheless, because of limitations in our dataset and localities where samples were collected, we cannot exclude the possibility of introgressive hybridization with closely-related taxa or that we might have simply overlooked an extant lineage.Particularly for the RBS, a peculiar curiosity arises from the existence of its closest taxon, the Red-tailed Shrike (Lanius phoenicuroides).This species breeds at the western edge of the RBS areal and was previously considered a subspecies of the RBS, it is known to hybridize with the RBS and, its field identification provides great challenges due to morphological similarities.We believe the inclusion of Red-tailed Shrike samples in our RBS dataset might have offered crucial details regarding the two-clade structure, but at this stage we could not obtain a satisfactory number of samples.Nonetheless, we are confident that future analysis will shed light on this aspect.
Concerning the genetic diversity of country specific RBS populations in our study, in Belgium, Czech Republic, The Netherlands, and Russia there are less singletons and haplotype diversity is lowest (Figure 3, Table S3), which accounts for lower genetic diversity.Nonetheless, we mention that in our datasets per country, most samples were collected in the same region, simply because of collaborator's availability.Relating to the RBS high genetic diversity across the studied populations, we underline two main hypotheses to explain this pattern: (i) Because our dataset revealed higher haplotype diversity in Bulgaria, Romania, France, and Latvia, we believe the glacial refugia for the European shrikes was in the Balkans, as for many other biota taxa [13].The presence of France among the countries with the highest diversity in our dataset might also indicate a refugia in the Italian Peninsula, but unfortunately our data does not offer additional evidence in this direction and as such, it will be overconfident to make further statements.However, our diversity indexes per country generally imply that, the further away the RBS population expanded, the lower the genetic diversity [7,30].Additionally, all RBS from western Europe migrate to Africa via the Balkans [65,78,79], even the birds which breed in the Iberian Peninsula [80].The former information, correlated with our genetic results, indicates a clear population expansion pattern from the Balkans towards western Europe.(ii) Secondly, for Belgium and The Netherlands, this lower genetic diversity can be attributed to population crashes in the second half of last century, mainly as a response of habitat loss and intensive farming practices.In Belgium, for example, in the 1960's the population crashed to only 600 pairs from the original numbers of about 5000 pairs [81], but recent counts indicate a total population of 3700 pairs, distributed only in Wallonia (in Flanders there are only 1 to 5 RBS pairs) [82].Therefore, our results suggest that the RBS population in western Europe is more vulnerable, in terms of population conservation.Nonetheless, we mention Spain, where although the RBS is suffering a steady decline, with some local populations having dropped by 95% [83,84], genetic diversity is still high.Yet, from a species conservation perspective, the current prevailing threat for the RBS comes from the hunting of migratory birds in the Mediterranean countries [27,[85][86][87], where, recent data revealed that in Egypt alone more than 30,000 RBS are taken each autumn [28].
To conclude, we emphasize the high genetic diversity for the breeding populations included in our study, the characteristic two-clade haplotype structure, and the apparent gene flow, which is reassuring when taking into account the hunting of thousands of RBS in the Mediterranean basin, the alarming rate of habitat loss across Europe and the climatic changes occurring along the migratory pathways.

A Blast from the Past: Pleistocene Upbringing
Population structure is determined by the genotype distribution in space and time and is the result of past events as well as ongoing processes [88].Across the western Palearctic, the genetic population structure for majority of taxa has been shaped by the complex climatic shifts, vegetation composition, and refugium locations in the Quaternary [13].For many species of plants, amphibians, reptiles, birds, and mammals, this translates into genetic admixture, triggered by the numerous glaciations which forced populations to retreat to refugia, where the fusion of different groups of individuals took place; subsequently, genetic admixture was reinforced when breeding populations from different glacial refugia fused.Overall, for plants species and animals with limited mobility, some population differentiation has occurred, and individuals can be grouped in clades belonging to their current geographic location or clades belonging to certain refugia, with the highest diversity concentrated near the glacial refugia e.g., Balkans, as observed in our RBS study as well.The Black Alder (Alnus glutinosa) [89] and oak trees (genus Quercus) [90] have also spread to North Europe from the Balkan Peninsula.Schrimpf et al. revealed that South-east Europe is the hotspot for genetic diversity of the Noble Crayfish (Astacus astacus) [91].In the European Fire-bellied Toads (Bombina bombina & Bombina variegata), nuclear and mitochondrial phylogeography indicate region specific clades across their breeding range [92].For reptiles, the European Pond Turtle (Emys orbicularis) shows high intraspecific differentiation, with haplotypes limited to certain geographic areas [93].Furthermore, Whipsnakes also show population differentiation and it appears that they colonized Europe from the Balkan refugia [94,95].Wall Lizard genetic population structure likewise indicate the Balkans as refugia and diversity hotspot [96].Among mammals, the Brown Bear (Ursus arctos) haplotypes can be divided in western and eastern clades, but only Romanian bears belong to both clades [97], giving further evidence of the high genetic diversity which occurs in South-east Europe.
Whereas in the northern temperate regions, many taxa show genetic admixture caused by Pleistocene glaciations and recent expansion from the refugia, across the tropical areas, bird species, especially the sedentary ones, show a clearly differentiated structure, with haplotypes being safely assigned to certain regions.We mention the African Blue Tits (Cyanistes teneriffae) in the Canary Islands, which belong to six haplotype groups, one for each of the seven islands in the archipelago (Lanzarote having just one haplotype) [98].For the South American continent, a recent study revealed that Polioptila gnatcatchers have a remarkable phylogeographic differentiation, mainly attributed to the landscape changes which occurred during the late Pliocene [99].
In conclusion, the RBS mitochondrial phylogeography indicates genetic admixture among the sampled populations, a common feature for bird species inhabiting the western Palearctic.Moreover, for RBS, we recorded higher genetic diversity in South-east Europe, further consolidating the theory that, during the Pleistocene glaciations, the biggest refugia for this species was in the Balkan Peninsula.From this refugia, the species has subsequently expanded and colonized the whole of Europe.

Conclusions
To the best of our knowledge, this study is the first to shed light upon the genetic population structure of RBSs.Across the western Palearctic, breeding RBSs show high genetic diversity, with a total of 76 haplotypes and 0.96 haplotype diversity in our dataset from 132 birds in 14 countries.To our surprise, two major lineages are apparent, with no geographic distinction, similar to the situation of the Common Redstart and the Eurasian Collared Dove.The random geographic distribution of haplotypes, a general indicator for population admixture, is common and widespread among European biota.This situation originates in the Pleistocene era, whose climate was marked by dramatic oscillations between periods of extreme cold and extensive ice coverage and periods of relative warmth.These climatic rotations, which implied a 'back-and-forth' process of population expansion and subsequent retreat to refugia, represent the driving factor that shaped today's genetic admixture in RBS across their breeding range in the western Palearctic.In future research, we will try to elucidate the evolutionary past of the RBS by using genomic next generation sequencing analyses, which allow a much higher resolution than mtDNA sequences [100].

Diversity 2019 , 18 Figure 1 .
Figure 1.Minimum-spanning haplotype network based on the concatenated Cytb and CO1 genes.Circle size reflects frequency of the haplotypes.The small black circles indicate missing node haplotypes.Each small vertical line traversing the connecting line between two haplotypes corresponds to one mutation step.

Figure 1 .
Figure 1.Minimum-spanning haplotype network based on the concatenated Cytb and CO1 genes.Circle size reflects frequency of the haplotypes.The small black circles indicate missing node haplotypes.Each small vertical line traversing the connecting line between two haplotypes corresponds to one mutation step.

Figure 2 .
Figure 2. Geographic distribution of the two main haplotype clades, based on concatenated Cytb and CO1, as represented by our mtDNA data from 132 birds in 14 countries.On the map, the unit near each circle indicates number of samples for the respective population.

Figure 3 .
Figure 3. Geographic distribution of the 76 haplotypes based on concatenated Cytb and CO1, as represented by our mtDNA data from 132 birds in 14 countries.On the right-side legend panel, each shared haplotype is assigned with a different color, whereas all singleton haplotypes are marked with grey.On the map, the unit near each circle indicates number of samples for the respective population.

Figure 2 .
Figure 2. Geographic distribution of the two main haplotype clades, based on concatenated Cytb and CO1, as represented by our mtDNA data from 132 birds in 14 countries.On the map, the unit near each circle indicates number of samples for the respective population.

Figure 2 .
Figure 2. Geographic distribution of the two main haplotype clades, based on concatenated Cytb and CO1, as represented by our mtDNA data from 132 birds in 14 countries.On the map, the unit near each circle indicates number of samples for the respective population.

Figure 3 .
Figure 3. Geographic distribution of the 76 haplotypes based on concatenated Cytb and CO1, as represented by our mtDNA data from 132 birds in 14 countries.On the right-side legend panel, each shared haplotype is assigned with a different color, whereas all singleton haplotypes are marked with grey.On the map, the unit near each circle indicates number of samples for the respective population.

Figure 3 .
Figure 3. Geographic distribution of the 76 haplotypes based on concatenated Cytb and CO1, as represented by our mtDNA data from 132 birds in 14 countries.On the right-side legend panel, each shared haplotype is assigned with a different color, whereas all singleton haplotypes are marked with grey.On the map, the unit near each circle indicates number of samples for the respective population.

Table 1 .
Molecular diversity indexes for the 14 populations of breeding Red-backed Shrike (RBS), obtained from the concatenated genes cytochrome oxidase B (Cytb) and cytochrome c oxidase subunit 1 (CO1).

Table 2 .
Population pairwise F ST indicators.Values which indicate low genetic differentiation (<0.05) are typed in bold.The rest of the values indicate moderate genetic differentiation (0.05 to 0.25).Units with significant p-values (<0.05) are marked with "+" above the diagonal.

Table 3 .
Estimates of evolutionary divergence (base substitutions per site) over sequence pairs between the two main L. collurio clades, versus the other five Laniidae species and the Corvidae.