Examining genetic structure in a bogus yucca moth: a sequential approach to phylogeography

Biological and biochemical diversity in different biotypes of spotted stem borer, Chilo partellus (Swinhoe) in India

Because of variation in incidence and severity of damage by Chilo partellus (Swinhoe) in different geographical regions, it is difficult to identify stable sources of resistance against this pest. Therefore, the present studies were undertaken on biological attributes (damage in resistant and susceptible genotypes, survival and development) and biochemical profiles (amino acids and lipophilic compound) of C. partellus populations from eight geographical regions to understand it's population structure in India. There was a significant variation in biological attributes and biochemical profiles of C. partellus populations from different geographical regions. Based on virulence and biological attributes, similarity index placed the C. partellus populations in five groups. Likewise, lipophilic and amino acid profiling also placed the C. partellus populations in five groups. However, the different clusters based on biological and biochemical attributes did not include populations from the same regions. Similarity index based on virulence, biological attributes, and amino acids and lipophilic profiles placed the C. partellus populations in six groups. The C. partellus populations from Hisar, Hyderabad, Parbhani and Coimbatore were distinct from each other, indicating that there are four biotypes of C. partellus in India. The results suggested that sorghum and maize genotypes need to be tested against these four populations to identify stable sources of resistance. However, there is a need for further studies to establish the restriction in gene flow through molecular approaches across geographical regions to establish the distinctiveness of different biotypes of C. partellus in India.

Phylogeography and Population Structure of the Mediterranean Corn Borer, Sesamia nonagrioides (Lepidoptera: Noctuidae), Across Its Geographic Range

Sesamia nonagrioides (Lefèbvre) (Lepidoptera: Noctuidae), is a widespread insect pest in Africa, the Middle East, and Europe. However, its pest status varies across its distribution range. It is a major pest of maize in Europe and of sugarcane in Iran. In Africa, it is a major pest of maize in West Africa but not considered as a pest in East Africa. Recent surveys conducted in 2015 recorded S. nonagrioides to be a major pest of sugarcane in Ethiopia and reported the species for the first time in Botswana, outside its known geographic range. The genetic relationship of these records with the previously recorded population of S. nonagrioides was investigated using the cytochrome oxidase subunit I region of the mitochondrial genome. In total, 113 individuals across the geographic range of the species were analyzed and 63 haplotypes were identified. Phylogenetic analysis separated the populations into two clades with no distinct geographic distribution pattern. The genetic differentiation was also not associated with host plants and geographic distances. Results of the molecular analysis revealed the long-time establishment of S. nonagrioides population in Botswana and identified the newly recorded sugarcane population from Ethiopia as part of the wild host population in the country. The phylogeographic patterns observed among population of S. nonagrioides have probably been shaped by Pleistocene's climatic oscillations and geographic range expansions from different refugia with secondary contact and admixture. Possible reasons for the host-plant expansion by the Ethiopian population are discussed.

What happened in the South American Gran Chaco? Diversification of the endemic frog genus Lepidobatrachus Budgett, 1899 (Anura: Ceratophryidae)

The Chaco is one the most neglected and least studied regions of the world. This highly-seasonal semiarid biome is an extensive continuous plain without any geographic barrier, and in spite of its high species diversity, the events and processes responsible have never been assessed. Miocene marine introgressions and Pleistocene glaciations have been mentioned as putative drivers of diversification for some groups of vertebrates in adjacent biomes of southern South America. Here we used multilocus data (one mitochondrial and six nuclear loci) from the three species of the endemic frog genus Lepidobatrachus (Lepidobatrachus asper, Lepidobatrachus laevis, and Lepidobatrachus llanensis) to determine if any of the historical events suggested as drivers of vertebrate diversification in southern South America are related to the diversification of the genus and if the Chaco is indeed a biome without barriers. Using fossil calibration in a coalescent framework we estimated that the genus diversified in the second half of the Miocene, coinciding with marine introgressions. Genetic patterns and historical demography suggest an important role of old archs and cratons as refuges during floods. In one species of the genus, L. llanensis, genetic structure reveals some breaks along the landscape, the main one of which corresponds to an area of the central Chaco that may act as a climatic barrier. Additionally, we found differential effects of the main Chacoan rivers on species of Lepidobatrachus that could be related to the time of persistence of populations in the areas influenced by these rivers.

Specialization in the yucca-yucca moth obligate pollination mutualism: A role for antagonism?

PREMISE OF THE STUDY

Specialized brood pollination systems involve both mutualism and antagonism in the overall interaction and have led to diversification in both plants and insects. Although largely known for mutualism, the role of the antagonistic side of the interaction in these systems has been overlooked. Specialization may be driven by plant defenses to feeding by the insect larvae that consume and kill developing plant ovules. The interaction among yuccas and yucca moths is cited as a classic example of the importance of mutualism in specialization and diversification. Pollinators moths are very host specific, but whether this specificity is due to adult pollination ability or larval feeding ability is unclear. Here, I test the potential role of antagonism in driving specialization among yuccas and yucca moths.

METHODS

I examined the ability of the most-polyphagous yucca moth pollinator, Tegeticula yuccasella, to pollinate and develop on five Yucca species used across its range. Yucca species endemic to the Great Plains and Texas were transplanted to a common garden in Syracuse, New York and exposed to the local pollinator moth population over 3 years.

KEY RESULTS

Local moths visited all but one of the Yucca species, but had drastically lower rates of successful larval development on non-natal Yucca species in comparison to the local host species.

CONCLUSION

Specialization in many brood pollination systems may be strongly influenced by the antagonistic rather than the mutualistic side of the overall interaction, suggesting that antagonistic coevolution is a possible source of diversification.

The roles of geography and founder effects in promoting host-associated differentiation in the generalist bogus yucca moth Prodoxus decipiens

There is ample evidence that host shifts in plant-feeding insects have been instrumental in generating the enormous diversity of insects. Changes in host use can cause host-associated differentiation (HAD) among populations that may lead to reproductive isolation and eventual speciation. The importance of geography in facilitating this process remains controversial. We examined the geographic context of HAD in the wide-ranging generalist yucca moth Prodoxus decipiens. Previous work demonstrated HAD among sympatric moth populations feeding on two different Yucca species occurring on the barrier islands of North Carolina, USA. We assessed the genetic structure of P. decipiens across its entire geographic and host range to determine whether HAD is widespread in this generalist herbivore. Population genetic analyses of microsatellite and mtDNA sequence data across the entire range showed genetic structuring with respect to host use and geography. In particular, genetic differentiation was relatively strong between mainland populations and those on the barrier islands of North Carolina. Finer scale analyses, however, among sympatric populations using different host plant species only showed significant clustering based on host use for populations on the barrier islands. Mainland populations did not form population clusters based on host plant use. Reduced genetic diversity in the barrier island populations, especially on the derived host, suggests that founder effects may have been instrumental in facilitating HAD. In general, results suggest that the interplay of local adaptation, geography and demography can determine the tempo of HAD. We argue that future studies should include comprehensive surveys across a wide range of environmental and geographic conditions to elucidate the contribution of various processes to HAD.

Testing for coevolutionary diversification: linking pattern with process

Coevolutionary diversification is cited as a major mechanism driving the evolution of diversity, particularly in plants and insects. However, tests of coevolutionary diversification have focused on elucidating macroevolutionary patterns rather than the processes giving rise to such patterns. Hence, there is weak evidence that coevolution promotes diversification. This is in part due to a lack of understanding about the mechanisms by which coevolution can cause speciation and the difficulty of integrating results across micro- and macroevolutionary scales. In this review, we highlight potential mechanisms of coevolutionary diversification, outline approaches to examine this process across temporal scales, and propose a set of minimal requirements for demonstrating coevolutionary diversification. Our aim is to stimulate research that tests more rigorously for coevolutionary diversification.

Species delimitation and morphological divergence in the scorpion Centruroides vittatus (Say, 1821): insights from phylogeography

Scorpion systematics and taxonomy have recently shown a need for revision, partially due to insights from molecular techniques. Scorpion taxonomy has been difficult with morphological characters as disagreement exists among researchers with character choice for adequate species delimitation in taxonomic studies. Within the family Buthidae, species identification and delimitation is particularly difficult due to the morphological similarity among species and extensive intraspecific morphological diversity. The genus Centruroides in the western hemisphere is a prime example of the difficulty in untangling the taxonomic complexity within buthid scorpions. In this paper, we present phylogeographic, Ecological Niche Modeling, and morphometric analyses to further understand how population diversification may have produced morphological diversity in Centruroides vittatus (Say, 1821). We show that C. vittatus populations in the Big Bend and Trans-Pecos region of Texas, USA are phylogeographically distinct and may predate the Last Glacial Maximum (LGM). In addition, we suggest the extended isolation of Big Bend region populations may have created the C. vittatus variant once known as C. pantheriensis.

Cytogeography of the Humifusa clade of Opuntia s.s. Mill. 1754 (Cactaceae, Opuntioideae, Opuntieae): correlations with pleistocene refugia and morphological traits in a polyploid complex

Ploidy has been well studied and used extensively in the genus Opuntia to determine species boundaries, detect evidence of hybridization, and infer evolutionary patterns. We carried out chromosome counts for all members of the Humifusa clade to ascertain whether geographic patterns are associated with differences in ploidy. We then related chromosomal data to observed morphological variability, polyploid formation, and consequently the evolutionary history of the clade. We counted chromosomes of 277 individuals from throughout the ranges of taxa included within the Humifusa clade, with emphasis placed on the widely distributed species, Opuntia humifusa (Raf.) Raf., 1820 s.l. and Opuntia macrorhiza Engelm., 1850 s.l. We also compiled previous counts made for species in the clade along with our new counts to plot geographic distributions of the polyploid and diploid taxa. A phylogeny using nuclear ribosomal ITS sequence data was reconstructed to determine whether ploidal variation is consistent with cladogenesis. We discovered that diploids of the Humifusa clade are restricted to the southeastern United States (U.S.), eastern Texas, and southeastern New Mexico. Polyploid members of the clade, however, are much more widely distributed, occurring as far north as the upper midwestern U.S. (e.g., Michigan, Minnesota, Wisconsin). Morphological differentiation, although sometimes cryptic, is commonly observed among diploid and polyploid cytotypes, and such morphological distinctions may be useful in diagnosing possible cryptic species. Certain polyploid populations of Opuntia humifusa s.l. and Opuntia macrorhiza s.l., however, exhibit introgressive morphological characters, complicating species delineations. Phylogenetically, the Humifusa clade forms two subclades that are distributed, respectively, in the southeastern U.S. (including all southeastern U.S. diploids, polyploid Opuntia abjecta Small, 1923, and polyploid Opuntia pusilla (Haw.) Haw., 1812) and the southwestern U.S. (including all southwestern U.S. diploids and polyploids). In addition, tetraploid Opuntia humifusa s.l., which occurs primarily in the eastern U.S., is resolved in the southwestern diploid clade instead of with the southeastern diploid clade that includes diploid Opuntia humifusa s.l. Our results not only provide evidence for the polyphyletic nature of Opuntia humifusa and Opuntia macrorhiza, suggesting that each of these represents more than one species, but also demonstrate the high frequency of polyploidy in the Humifusa clade and the major role that genome duplication has played in the diversification of this lineage of Opuntia s.s. Our data also suggest that the southeastern and southwestern U.S. may represent glacial refugia for diploid members of this clade and that the clade as a whole should be considered a mature polyploid species complex. Widespread polyploids are likely derivatives of secondary contact among southeastern and southwestern diploid taxa as a result of the expansion and contraction of suitable habitat during the Pleistocene following glacial and interglacial events.

Inference of population history by coupling exploratory and model-driven phylogeographic analyses

Understanding the nature, timing and geographic context of historical events and population processes that shaped the spatial distribution of genetic diversity is critical for addressing questions relating to speciation, selection, and applied conservation management. Cladistic analysis of gene trees has been central to phylogeography, but when coupled with approaches that make use of different components of the information carried by DNA sequences and their frequencies, the strength and resolution of these inferences can be improved. However, assessing concordance of inferences drawn using different analytical methods or genetic datasets, and integrating their outcomes, can be challenging. Here we overview the strengths and limitations of different types of genetic data, analysis methods, and approaches to historical inference. We then turn our attention to the potentially synergistic interactions among widely-used and emerging phylogeographic analyses, and discuss some of the ways that spatial and temporal concordance among inferences can be assessed. We close this review with a brief summary and outlook on future research directions.

Host-associated divergence and incipient speciation in the yucca moth Prodoxus coloradensis (Lepidoptera: Prodoxidae) on three species of host plants

A wide range of evolutionary processes have been implicated in the diversification of yuccas and yucca moths, which exhibit ecological relationships that extend from obligate plant-pollinator mutualisms to commensalist herbivory. Prodoxus coloradensis (Lepidoptera: Prodoxidae) is a yucca moth, which feeds on the flowering stalks of three Yucca species as larvae, but does not provide pollination service. To test for evidence of host-associated speciation, we examined the genetic structure of P. coloradensis using mitochondrial (cytochrome oxidase I) and nuclear (elongation factor 1 alpha) DNA sequence data. Multilocus coalescent simulations indicate that moths on different host plant species are characterized by recent divergence and low levels of effective migration, with large effective population sizes and considerable retention of shared ancestral polymorphism. Although geographical distance explains a proportion of the mitochondrial and nuclear DNA variation among moths on different species of Yucca, the effect of host specificity on genetic distance remains significant after accounting for spatial isolation. The results of this study indicate that differentiation within P. coloradensis is consistent with the evolution of incipient species affiliated with different host plants, potentially influenced by sex-biased dispersal and female philopatry.

Population genetic analysis of an Eastern U.S. wheat powdery mildew population reveals geographic subdivision and recent common ancestry with U.K. and Israeli populations

The structure of the U.S. wheat powdery mildew population (Blumeria graminis f. sp. tritici) has not been previously investigated, and the global evolutionary history of B. graminis f. sp. tritici is largely unknown. After gathering 141 single-ascosporic B. graminis f. sp. tritici isolates from 10 eastern U.S. locations, 34 isolates from the United Kingdom, and 28 isolates from Israel, we analyzed pathogen population structure using presumptively neutral markers. DNA was extracted from conidia, primers for 12 "housekeeping" genes were designed, and amplicons were examined for polymorphism. Four genes were found to contain a total of 12 single-nucleotide polymorphisms in the U.S. population and were also analyzed in the U.K. and Israeli populations. In total, 25 haplotypes were inferred from the four concatenated genes, with 2 haplotypes comprising over 70% of the U.S. population. Using Hudson's tests and analysis of molecular variance, we found the wheat mildew isolates subdivided into four groups corresponding to distinct regions: the mid-Atlantic United States, the southern United States, the United Kingdom, and Israel. Genotypic diversity was greatest in samples from the United Kingdom, Israel, Virginia, and Kingston, NC. Using rarefaction, a procedure that compensates for differing sample sizes when estimating population richness and diversity, we found that cooler locations with greater conduciveness to regular powdery mildew epidemics had the greatest haplotype richness. Our results suggest that the eastern U.S. B. graminis f. sp. tritici population is young, descended recently from Old World populations with isolation and genetic drift, and is currently subdivided into northern and southern subpopulations.

Evolutionary history of the endangered fish Zoogoneticus quitzeoensis (Bean, 1898) (Cyprinodontiformes: Goodeidae) using a sequential approach to phylogeography based on mitochondrial and nuclear DNA data

Background

Tectonic, volcanic and climatic events that produce changes in hydrographic systems are the main causes of diversification and speciation of freshwater fishes. Elucidate the evolutionary history of freshwater fishes permits to infer theories on the biotic and geological evolution of a region, which can further be applied to understand processes of population divergence, speciation and for conservation purposes. The freshwater ecosystems in Central Mexico are characterized by their genesis dynamism, destruction, and compartmentalization induced by intense geologic activity and climatic changes since the early Miocene. The endangered goodeid Zoogoneticus quitzeoensis is widely distributed across Central México, thus making it a good model for phylogeographic analyses in this area.

Results

We addressed the phylogeography, evolutionary history and genetic structure of populations of Z. quitzeoensis through a sequential approach, based on both microsatellite and mitochondrial cytochrome b sequences. Most haplotypes were private to particular locations. All the populations analysed showed a remarkable number of haplotypes. The level of gene diversity within populations was H¯_d = 0.987 (0.714 – 1.00). However, in general the nucleotide diversity was low, π = 0.0173 (0.0015 – 0.0049). Significant genetic structure was found among populations at the mitochondrial and nuclear level (Φ_ST = 0.836 and F_ST = 0.262, respectively). We distinguished two well-defined mitochondrial lineages that were separated ca. 3.3 million years ago (Mya). The time since expansion was ca. 1.5 × 10⁶ years ago for Lineage I and ca. 860,000 years ago for Lineage II. Also, genetic patterns of differentiation, between and within lineages, are described at different historical timescales.

Conclusion

Our mtDNA data indicates that the evolution of the different genetic groups is more related to ancient geological and climatic events (Middle Pliocene, ca. 3.3 Mya) than to the current hydrographic configuration of the basins. In general, mitochondrial and nuclear data supported the same relationships between populations, with the exception of some reduced populations in highly polluted basins (Lower Lerma River), where the effects of genetic drift are suggested by the different analyses at the nuclear and mitochondrial level. Further, our findings are of special interest for the conservation of this endangered species.

Phylogeography of the tree lizard,Urosaurus ornatus: responses of populations to past climate change

Isolation due to both geological barriers and range contractions during the Pleistocene glacial maxima has been an important cause of diversification of arid-adapted species in the North American deserts. Tree lizards, Urosaurus ornatus, are distributed across much of the southwestern arid regions and can tolerate a wide range of environments. Thus, they may have avoided large-scale shifts in distribution caused by Pleistocene climate change and any subsequent evolutionary impacts. Cytochrome b sequences were sampled from U. ornatus across the northern part of their range to test if current structure of these populations resulted from post-Pleistocene range expansion and habitat fragmentation, or prior geological isolation. Phylogenetic analyses found geographical structuring of populations consistent with a model of long-term geographical isolation corresponding to each of the desert regions. The two post-Pleistocene hypotheses were not well supported as estimated times of divergence predated the retreat of the last continental ice sheet. Populations in different regions were impacted by different processes. Southern populations of U. ornatus appear to have remained largely independent of more derived northern and eastern populations during Pleistocene climate change, while populations in regions containing more derived populations showed evidence of more recent range expansion (Colorado Plateau). As populations of U. ornatus attest to, the complex and dynamic history of the southwestern USA has left a deep-rooted and multifaceted imprint on genetic and phylogeographical structure of the species living there.

Phylogeography and population genetics of the maize stalk borer Busseola fusca (Lepidoptera, Noctuidae) in sub-Saharan Africa

The population genetics and phylogeography of African phytophagous insects have received little attention. Some, such as the maize stalk borer Busseola fusca, display significant geographic differences in ecological preferences that may be congruent with patterns of molecular variation. To test this, we collected 307 individuals of this species from maize and cultivated sorghum at 52 localities in West, Central and East Africa during the growing season. For all collected individuals, we sequenced a fragment of the mitochondrial cytochrome b. We tested hypotheses concerning the history and demographic structure of this species. Phylogenetic analyses and nested clade phylogeographic analyses (NCPA) separated the populations into three mitochondrial clades, one from West Africa, and two--Kenya I and Kenya II--from East and Central Africa. The similar nucleotide divergence between clades and nucleotide diversity within clades suggest that they became isolated at about the same time in three different refuges in sub-Saharan Africa and have similar demographic histories. The results of mismatch distribution analyses were consistent with the demographic expansion of these clades. Analysis of molecular variance (amova) indicated a high level of geographic differentiation at different hierarchical levels. NCPA suggested that the observed distribution of haplotypes at several hierarchical levels within the three major clades is best accounted for by restricted gene flow with isolation by distance. The domestication of sorghum and the introduction of maize in Africa had no visible effect on the geographic patterns observed in the B. fusca mitochondrial genome.

Cryptic Neogene vicariance and Quaternary dispersal of the red-spotted toad (Bufo punctatus): insights on the evolution of North American warm desert biotas

We define the geographical distributions of mitochondrial DNA (mtDNA) lineages embedded within a broadly distributed, arid-dwelling toad, Bufo punctatus. These patterns were evaluated as they relate to hypothesized vicariant events leading to the formation of desert biotas within western North America. We assessed mtDNA sequence variation among 191 samples from 82 sites located throughout much of the species' range. Parsimony-based haplotype networks of major identified lineages were used in nested clade analysis (NCA) to further elucidate and evaluate shallow phylogeographic patterns potentially associated with Quaternary (Pleistocene-Holocene) vicariance and dispersal. Phylogenetic analyses provided strong support for three monophyletic lineages (clades) within B. punctatus. The geographical distributions of the clades showed little overlap and corresponded to the general boundaries of the Peninsular Desert, and two continental desert regions, Eastern (Chihuahuan Desert-Colorado Plateau) and Western (Mojave-Sonoran deserts), geographically separated along the Rocky Mountains and Sierra Madre Occidental. The observed divergence levels and congruence with postulated events in earth history implicate a late Neogene (latest Miocene-early Pliocene) time frame for separation of the major mtDNA lineages. Evaluation of nucleotide and haplotype diversity and interpretations from NCA reveal that populations on the Colorado Plateau resulted from a recent, likely post-Pleistocene, range expansion from the Chihuahuan Desert. Dispersal across historical barriers separating major continental clades appear to be recent, resulting in secondary contacts in at least two areas. Given the observed contact between major clades, we speculated as to why the observed deep phylogeographic structure has not been eroded during the multiple previous interglacials of the Pleistocene.

Age, phylogeography and population structure of the microendemic banded spring snail, Mexipyrgus churinceanus

Recent theoretical and empirical studies of phylogeography and population structure indicate that many processes influence intraspecific evolutionary history. The present study represents the first examination of various forces influencing the spatial and temporal patterns of sequence variation in the freshwater Mexican banded spring snail, Mexipyrgus churinceanus. This snail occurs in one of the most critically endangered centres of freshwater endemism, the desert ecosystem of Cuatro Ciénegas. From cytochrome b mtDNA sequence variation, there is strong evidence of long-term isolation of three regions, suggesting that these regions represent evolutionarily distinct lineages. Molecular clock estimates of clade age indicate a time to most recent common ancestor of approximately 2.5 million years ago (Ma). The three regions differ considerably in the historical and demographic forces affecting population structure. The western populations have extremely low mtDNA diversity consistent with a severe bottleneck dating to 50,000 years before present (bp). The nearby Rio Mesquites drainage is characterized by fragmentation events, restricted gene flow with isolation by distance, and higher levels of mtDNA polymorphism. These patterns are consistent with the long-term stability of this drainage along with habitat heterogeneity and brooding contributing to population isolation and restricted gene flow. Southeastern populations show evidence of range expansion and a strong influence of genetic drift. Migration rates between drainages indicate very little gene flow between drainages except for asymmetric migration from the Rio Mesquites into both western and southeastern drainages.

Phylogeographical history of the sponge Crambe crambe (Porifera, Poecilosclerida): range expansion and recent invasion of the Macaronesian islands from the Mediterranean Sea

We studied sequence variation in the nuclear ribosomal internal transcribed spacers (ITS-1 and ITS-2) in 111 individuals from 11 populations/localities of the sponge Crambe crambe across the core species range in the western Mediterranean Sea and Atlantic Ocean. We report the first confirmed instance of intragenomic variation in sponges. Phylogeographical, nested clade and population genetic analyses were used to elucidate the species' evolutionary history. The study revealed highly structured populations affected by restricted gene flow and isolation-by-distance. A contiguous range expansion in the whole distribution area of the sponge was inferred. Phylogenetic analyses indicate a recent origin of most sequence types that could be explained by a recent origin of the species or a by recent bottleneck event in the studied area. A recent expansion of the distribution range to the Macaronesian region from the Mediterranean Sea was also detected, suggesting that C. crambe was recently introduced from the Mediterranean Sea to the Atlantic Ocean via human-mediated transport, and that the pattern observed is not the result of a natural biogeographical relationship between these zones.

Phylogeography, genetic structure and diversity in the endangered bearded vulture (Gypaetus barbatus, L) as revealed by mitochondrial DNA

Bearded vulture populations in the Western Palearctic have experienced a severe decline during the last two centuries that has led to the near extinction of the species in Europe. In this study we analyse the sequence variation at the mitochondrial control region throughout the species range to infer its recent evolutionary history and to evaluate the current genetic status of the species. This study became possible through the extensive use of museum specimens to study populations now extinct. Phylogenetic analysis revealed the existence of two divergent mitochondrial lineages, lineage A occurring mainly in Western European populations and lineage B in African, Eastern European and Central Asian populations. The relative frequencies of haplotypes belonging to each lineage in the different populations show a steep East-West clinal distribution with maximal mixture of the two lineages in the Alps and Greece populations. A genealogical signature for population growth was found for lineage B, but not for lineage A; futhermore the Clade B haplotypes in western populations and clade A haplo-types in eastern populations are recently derived, as revealed by their peripheral location in median-joining haplotype networks. This phylogeographical pattern suggests allopatric differentiation of the two lineages in separate Mediterranean and African or Asian glacial refugia, followed by range expansion from the latter leading to two secondary contact suture zones in Central Europe and North Africa. High levels of among-population differentiation were observed, although these were not correlated with geographical distance. Due to the marked genetic structure, extinction of Central European populations in the last century re-sulted in the loss of a major portion of the genetic diversity of the species. We also found direct evidence for the effect of drift altering the genetic composition of the remnant Pyrenean population after the demographic bottleneck of the last century. Our results argue for the management of the species as a single population, given the apparent ecological exchangeability of extant stocks, and support the ongoing reintroduction of mixed ancestry birds in the Alps and planned reintroductions in Southern Spain.