Phylogeography of the greater horseshoe bat, Rhinolophus ferrumequinum: contrasting results from mitochondrial and microsatellite data

Phylogeographical studies are typically based on haplotype data, occasionally on nuclear markers such as microsatellites, but rarely combine both. This is unfortunate because the use of markers with contrasting modes of inheritance and rates of evolution might provide a more accurate and comprehensive understanding of a species' history. Here we present a detailed study of the phylogeography of the greater horseshoe bat, Rhinolophus ferrumequinum , using 1098 bp of the mitochondrial ND2 gene from 45 localities from across its Palaearctic range to infer population history. In addition, we re-analysed a large microsatellite data set available for this species and compared the results of both markers to infer population relationships and the historical processes influencing them. We show that mtDNA, the most popular marker in phylogeography studies, yielded a misleading result, and would have led us to conclude erroneously that a single expansion had taken place in Europe. Only by combining the mitochondrial and microsatellite data sets are we able to reconstruct the species' history and show two colonization events in Europe, one before the Last Glacial Maximum (LGM) and one after it. Combining markers also revealed the importance of Asia Minor as an ancient refugium for this species and a source population for the expansion of the greater horseshoe bat into Europe before the LGM.     

Understanding the phylogeographic patterns of European hedgehogs, Erinaceus concolor and E. europaeus using the MHC

The genome of the European hedgehog, Erinaceus concolor and E. europaeus, shows a strong signal of cycles of restriction to glacial refugia and postglacial expansion. Patterns of expansion, however, differ for mitochondrial DNA (mtDNA) and preliminary analysis of nuclear markers. In this study, we determine phylogeographic patterns in the hedgehog using two loci of the major histocompatibility complex (MHC), isolated for the first time in hedgehogs. These genes show long persistence times and high polymorphism in many species because of the actions of balancing selection. Among 84 individuals screened for variation, only two DQA alleles were identified in each species, but 10 DQB alleles were found in E. concolor and six in E. europaeus. A strong effect of demography on patterns of DQB variability is observed, with only weak evidence of balancing selection. While data from mtDNA clearly subdivide both species into monophyletic subgroups, the MHC data delineate only E. concolor into distinct subgroups, supporting the preliminary findings of other nuclear markers. Together with differences in variability, this suggests that the refugia history and/or expansion patterns of E. concolor and E. europaeus differ.     

Evolutionary patterns of the mitochondrial genome in the Moorish gecko, Tarentola mauritanica

A previous study on the evolutionary patterns of Tarentola mauritanica demonstrated that low levels of mitochondrial diversity observed in the European populations relative to nuclear markers were consistent with a selective sweep hypothesis. In order to unravel the mitochondrial evolutionary history in this European population and two other lineages of T. mauritanica (Iberian and North African clades), variation within 22 nearly complete mitogenomes was analyzed. Surprisingly, each clade seems to have a distinct evolutionary history; with both the European and Iberian clades presenting a decrease of polymorphism, which in the former is consistent with departure from neutrality of the mtDNA (positive or background selection), but in the latter seems to be the result of a bottleneck after a population expansion. The pattern exhibited by the North African clade seems to be a consequence of adaptation to certain mtDNA variants by positive selection.     

mtDNA variation predicts population size in humans and reveals a major Southern Asian chapter in human prehistory

The relative timing and size of regional human population growth following our expansion from Africa remain unknown. Human mitochondrial DNA (mtDNA) diversity carries a legacy of our population history. Given a set of sequences, we can use coalescent theory to estimate past population size through time and draw inferences about human population history. However, recent work has challenged the validity of using mtDNA diversity to infer species population sizes. Here we use Bayesian coalescent inference methods, together with a global data set of 357 human mtDNA coding-region sequences, to infer human population sizes through time across 8 major geographic regions. Our estimates of relative population sizes show remarkable concordance with the contemporary regional distribution of humans across Africa, Eurasia, and the Americas, indicating that mtDNA diversity is a good predictor of population size in humans. Plots of population size through time show slow growth in sub-Saharan Africa beginning 143-193 kya, followed by a rapid expansion into Eurasia after the emergence of the first non-African mtDNA lineages 50-70 kya. Outside Africa, the earliest and fastest growth is inferred in Southern Asia approximately 52 kya, followed by a succession of growth phases in Northern and Central Asia (approximately 49 kya), Australia (approximately 48 kya), Europe (approximately 42 kya), the Middle East and North Africa (approximately 40 kya), New Guinea (approximately 39 kya), the Americas (approximately 18 kya), and a second expansion in Europe (approximately 10-15 kya). Comparisons of relative regional population sizes through time suggest that between approximately 45 and 20 kya most of humanity lived in Southern Asia. These findings not only support the use of mtDNA data for estimating human population size but also provide a unique picture of human prehistory and demonstrate the importance of Southern Asia to our recent evolutionary past.     

Assessing inter‐drainage connections: patterns of genetic diversity in an Iberian cyprinid fish

The phylogeography of Iberian freshwater fish species strongly reflects the geomorphological history of Iberian water drainages and formation of mountain ranges. In particular, the Iberian Central Massif limits two ichthyogeographical districts based on fish distribution assemblages: the Northwest district, which comprises all river basins north of the Central massif, and the Central‐west district, which comprises the Tejo and Sado basins. Despite this scenario, the cyprinid Pseudochondrostoma polylepis occurs in both the Tejo and Mondego rivers (one either side of the Estrela Mountains, in the western end of the Central Massif). By analysing the population structure and demographic history of this straight‐mouth nase, based on both mtDNA and microsatellites, we have found evidence of a relatively recent common origin for both Mondego and Tejo populations. Shared haplotypes and lack of divergent lineages for mitochondrial marker, high levels of admixture and common ancestry at the nuclear level, together with signatures of asymmetrical gene flow from Tejo to Mondego suggest a Mondego colonization through a founder event with an origin in Tejo. We suggest that this colonization might have resulted from the occurrence of river connections/captures between Mondego and Tejo headwater tributaries in the recent past. With this study we reinforce that local‐scale geomorphological events and the adaptation/dispersal ability of a species to new environments can have a major impact on its demographic history, colonization processes, and ultimately its distribution. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2013, 109 , 656–669.     

Extensive sampling of polar bears (Ursus maritimus) in the Northwest Passage (Canadian Arctic Archipelago) reveals population differentiation across multiple spatial and temporal scales

As global warming accelerates the melting of Arctic sea ice, polar bears (Ursus maritimus) must adapt to a rapidly changing landscape. This process will necessarily alter the species distribution together with population dynamics and structure. Detailed knowledge of these changes is crucial to delineating conservation priorities. Here, we sampled 361 polar bears from across the center of the Canadian Arctic Archipelago spanning the Gulf of Boothia (GB) and M'Clintock Channel (MC). We use DNA microsatellites and mitochondrial control region sequences to quantify genetic differentiation, estimate gene flow, and infer population history. Two populations, roughly coincident with GB and MC, are significantly differentiated at both nuclear (F_ST = 0.01) and mitochondrial (Φ_ST = 0.47; F_ST = 0.29) loci, allowing Bayesian clustering analyses to assign individuals to either group. Our data imply that the causes of the mitochondrial and nuclear genetic patterns differ. Analysis of mtDNA reveals the matrilineal structure dates at least to the Holocene, and is common to individuals throughout the species’ range. These mtDNA differences probably reflect both genetic drift and historical colonization dynamics. In contrast, the differentiation inferred from microsatellites is only on the scale of hundreds of years, possibly reflecting contemporary impediments to gene flow. Taken together, our data suggest that gene flow is insufficient to homogenize the GB and MC populations and support the designation of GB and MC as separate polar bear conservation units. Our study also provide a striking example of how nuclear DNA and mtDNA capture different aspects of a species demographic history.     

Nuclear markers confirm taxonomic status and relationships among highly endangered and closely related right whale species

Right whales (genus: Eubalaena) are among the most endangered mammals, yet their taxonomy and phylogeny have been questioned. A phylogenetic hypothesis based on mitochondrial DNA (mtDNA) variation recently prompted a taxonomic revision, increasing the number of right whale species to three. We critically evaluated this hypothesis using sequence data from 13 nuclear DNA (nuDNA) loci as well as the mtDNA control region. Fixed diagnostic characters among the nuclear markers strongly support the hypothesis of three genetically distinct species, despite lack of any diagnostic morphological characters. A phylogenetics analysis of all data produced a strict consensus cladogram with strong support at nodes that define each right whale species as well as relationships among species. Results showed very little conflict among the individual partitions as well as congruence between the mtDNA and nuDNA datasets. These data clearly demonstrate the strength of using numerous independent genetic markers during a phylogenetics analysis of closely related species. In evaluating phylogenetic support contributed by individual loci, 11 of the 14 loci provided support for at least one of the nodes of interest to this study. Only a single marker (mtDNA control region) provided support at all four nodes. A study using any single nuclear marker would have failed to support the proposed phylogeny, and a strong phylogenetic hypothesis was only revealed by the simultaneous analysis of many nuclear loci. In addition, nu DNA and mtDNA data provided complementary levels of support at nodes of different evolutionary depth indicating that the combined use of mtDNA and nuDNA data is both practical and desirable.     

Comparative phylogeography and asymmetric hybridization between cryptic bat species

Cryptic speciation and hybridization are two key processes that affect the origin and maintenance of biodiversity and our ability to understand and estimate it. To determine how these two processes interact, we studied allopatric and sympatric colonies of two cryptic bat species (Eptesicus serotinus and Eptesicus isabellinus) with parapatric distribution in the Iberian Peninsula. These species are the main reservoir for the most commonly rabies virus found in bats in Europe: the European bat Lyssavirus type 1 (EBLV‐1). We used mtDNA and nuclear microsatellite markers to confirm the taxonomic status of both species and to show a more pronounced and geographically based genetic structure in E. isabellinus than in its sibling E. serotinus. Using approximate Bayesian computation (ABC), we inferred rapid range expansion in both species after the Last Glacial Maximum until reaching their present distributions. ABC analysis also supported interspecific differences in genetic diversity and structure, pointing to an earlier expansion of E. isabellinus northward. We found no evidence of mitochondrial introgression between species, but nuclear markers identified a male‐mediated ongoing asymmetric hybridization from E. isabellinus to E. serotinus (28% hybrids in E. serotinus and 5% in E. isabellinus) in the contact zone. Although none of the bats studied tested positive for Lyssavirus RNA, the asymmetric hybridization supports the potential for the recently suggested interspecific transmission of EBLV‐1 from E. isabellinus into E. serotinus.     

Phylogenetic and phylogeographic analysis of Iberian lynx populations

The Iberian lynx (Lynx pardinus), one of the world's most endangered cat species, is vulnerable due to habitat loss, increased fragmentation of populations, and precipitous demographic reductions. An understanding of Iberian lynx evolutionary history is necessary to develop rational management plans for the species. Our objectives were to assess Iberian lynx genetic diversity at three evolutionary timescales. First we analyzed mitochondrial DNA (mtDNA) sequence variation to position the Iberian lynx relative to other species of the genus LYNX: We then assessed the pattern of mtDNA variation of isolated populations across the Iberian Peninsula. Finally we estimated levels of gene flow between two of the most important remaining lynx populations (Do?ana National Park and the Sierra Morena Mountains) and characterized the extent of microsatellite locus variation in these populations. Phylogenetic analyses of 1613 bp of mtDNA sequence variation supports the hypothesis that the Iberian lynx, Eurasian lynx, and Canadian lynx diverged within a short time period around 1.53-1.68 million years ago, and that the Iberian lynx and Eurasian lynx are sister taxa. Relative to most other felid species, genetic variation in mtDNA genes and nuclear microsatellites were reduced in Iberian lynx, suggesting that they experienced a fairly severe demographic bottleneck. In addition, the effects of more recent reductions in gene flow and population size are being manifested in local patterns of molecular genetic variation. These data, combined with recent studies modeling the viability of Iberian lynx populations, should provide greater urgency for the development and implementation of rational in situ and ex situ conservation plans.     

Morphology and nuclear markers reveal extensive mitochondrial introgressions in the Iberian Wall Lizard species complex

Mitochondrial markers are still often used alone to identify evolutionary units, despite widespread evidence for processes such as incomplete lineage sorting or introgressive hybridization that may blur past population history. The combination of mitochondrial DNA data with other sources of information (morphology, nuclear genes) is a powerful tool to reveal when and why mitochondrial markers are potentially misleading. In this study, we evaluate the performance of mtDNA markers to unravel the evolutionary history of Spanish lizards from the Podarcis hispanicus species complex. We first uncover several cases of discordance between morphological and mitochondrial data in delimitation of taxa. To assess the origin of these discordances, we analysed the same populations using several independent nuclear loci. Both morphological and nuclear markers identified the same three evolutionary units in the region, while mitochondrial data revealed four deeply divergent lineages. We suggest here that the most likely scenario to explain this discordance is ancient mitochondrial introgression originating from a fourth evolutionary unit presently absent from the study area. Notably, this resulted in a complete replacement of the original lineage in a large part of the distribution of one of the taxa investigated. We discuss the potential evolutionary scenarios leading to this complete mitochondrial replacement and suggest why the previous studies have failed to recover the correct history of this species complex.     

Iberian red deer: paraphyletic nature at mtDNA but nuclear markers support its genetic identity

Red deer populations in the Iberian glacial refugium were the main source for postglacial recolonization and subspecific radiation in north‐western Europe. However, the phylogenetic history of Iberian red deer (Cervus elaphus hispanicus) and its relationships with northern European populations remain uncertain. Here, we study DNA sequences at the mitochondrial control region along with STR markers for over 680 specimens from all the main red deer populations in Spain and other west European areas. Our results from mitochondrial and genomic DNA show contrasting patterns, likely related to the nature of these types of DNA markers and their specific processes of change over time. The results, taken together, bring support to two distinct, cryptic maternal lineages for Iberian red deer that predated the last glacial maximum and that have maintained geographically well differentiated until present. Haplotype relationships show that only one of them contributed to the northern postglacial recolonization. However, allele frequencies of nuclear markers evidenced one main differentiation between Iberian and northern European subspecies although also supported the structure of both matrilines within Iberia. Thus, our findings reveal a paraphyletic nature for Iberian red deer but also its genetic identity and differentiation with respect to northern subspecies. Finally, we suggest that maintaining the singularity of Iberian red deer requires preventing not only restocking practices with red deer specimens belonging to other European populations but also translocations between both Iberian lineages.     

Microsatellite analysis supports mitochondrial phylogeography of the hellbender (Cryptobranchus alleganiensis)

We investigated genetic diversity of the hellbender (Cryptobranchus alleganiensis) throughout its range in the eastern US using nuclear markers and compared our results to a previously published mitochondrial analysis. A variety of nuclear markers, including protein-coding gene introns and microsatellites were tested but only microsatellites were variable enough for population level analysis. Microsatellite loci showed moderate among population sharing of alleles, in contrast to the reciprocal monophyly exhibited by mitochondrial DNA. However, analyses using F-statistics and Bayesian clustering algorithms showed considerable population subdivision and clustered hellbender populations into the same major groups as the mtDNA. The microsatellites combined with the mtDNA data suggest that gene flow is severely restricted or non-existent among eight major groups, and potentially among populations (rivers) within groups. The combined mtDNA and microsatellite data suggest that the currently recognized hellbender subspecies are paraphyletic. We suggest that the eight independent groups identified in our study should be managed as such, rather than basing conservation decisions on the two named subspecies of hellbender.     

Range instability leads to cytonuclear discordance in a morphologically cryptic ground squirrel species complex

The processes responsible for cytonuclear discordance frequently remain unclear. Here, we employed an exon capture data set and demographic methods to test hypotheses generated by species distribution models to examine how contrasting histories of range stability vs. fluctuation have caused cytonuclear concordance and discordance in ground squirrel lineages from the Otospermophilus beecheyi species complex. Previous studies in O. beecheyi revealed three morphologically cryptic and highly divergent mitochondrial DNA lineages (named the Northern, Central and Southern lineages based on geography) with only the Northern lineage exhibiting concordant divergence for nuclear genes. Here, we showed that these mtDNA lineages likely formed in allopatry during the Pleistocene, but responded differentially to climatic changes that occurred since the last interglacial (~120,000 years ago). We find that the Northern lineage maintained a stable range throughout this period, correlating with genetic distinctiveness among all genetic markers and low migration rates with the other lineages. In contrast, our results suggested that the Southern lineage expanded from Baja California Sur during the Late Pleistocene to overlap and potentially swamp a contracting Central lineage. High rates of intraspecific gene flow between Southern lineage individuals among expansion origin and expansion edge populations largely eroded Central ancestry from autosomal markers. However, male‐biased dispersal in this system preserved signals of this past hybridization and introgression event in matrilineal‐biased X‐chromosome and mtDNA markers. Our results highlight the importance of range stability in maintaining the persistence of phylogeographic lineages, whereas unstable range dynamics can increase the tendency for lineages to merge upon secondary contact.     

The role of mitochondrial introgression in illuminating the evolutionary history of Nearctic treefrogs

Inferring the evolutionary history of lineages often becomes difficult when gene histories are in conflict with each other. Introgression, for example, can cause DNA sequences from one species to be more similar to sequences of a different species and lead to incongruence amongst gene trees. However, incorporating congruent and incongruent locus‐specific phylogenetic estimates with the geographical distribution of lineages may provide valuable insight into evolutionary processes important to speciation. In this study, we investigated mitochondrial introgression within the Hyla eximia group to better understand its role in illuminating the evolutionary history and phylogeography of these treefrogs. We reconstructed and compared the matrilineal history of the Hyla eximia group with estimates of evolutionary history inferred from nuclear genes. We tested for introgression within the mitochondrial and nuclear genes using a posterior predictive checking approach. Reconstructions of the species tree based on the mitochondrial DNA (mtDNA) and nuclear DNA data were strongly discordant. Introgression between lineages was widespread in the mtDNA data set (145 occurrences amongst 11 of the 16 lineages), but uncommon in the nuclear genes (12 occurrences amongst four of the 16 lineages). Nonetheless, the geographical structuring of mtDNA within species provides valuable information on biogeographical areas, ancient areas of hybridization, and unique histories of lineages within the H. eximia group. These results suggest that the combination of nuclear, mitochondrial, and spatial information can provide a more complete picture of ‘how evolutionary history played out’, particularly in cases where mitochondrial introgression is known to occur. © 2014 The Linnean Society of London     

Patterns of nuclear and mitochondrial DNA variation in Iberian populations of <Emphasis Type="Italic">Emys orbicularis</Emphasis> (Emydidae): conservation implications

The European pond turtle (Emys orbicularis) is threatened and in decline in several regions of its natural range, due to habitat loss combined with population fragmentation. In this work, we have focused our efforts on studying the genetic diversity and structure of Iberian populations with a fine-scale sampling (254 turtles in 10 populations) and a representation from North Africa and Balearic island populations. Using both nuclear and mitochondrial markers (seven microsatellites, ∼1048 bp nDNA and ∼1500 bp mtDNA) we have carried out phylogenetic and demographic analyses. Our results show low values of genetic diversity at the mitochondrial level although our microsatellite dataset revealed relatively high levels of genetic variability with a latitudinal genetic trend decreasing from southern to northern populations. A moderate degree of genetic differentiation was estimated for Iberian populations (genetic distances, F _ST values and clusters in the Bayesian analysis). The results in this study combining mtDNA and nDNA, provide the most comprehensive population genetic data for E. orbicularis in the Iberian Peninsula. Our results suggest that Iberian populations within the Iberian–Moroccan lineage should be considered as a single subspecies with five management units, and emphasize the importance of habitat management rather than population reinforcement (i.e. captive breeding and reintroduction) in this long-lived species.     

Genealogy of the nuclear beta-fibrinogen locus in a highly structured lizard species: comparison with mtDNA and evidence for intragenic recombination in the hybrid zone

The study of nuclear genealogies in natural populations of nonmodel organisms is expected to provide novel insights into the evolutionary history of populations, especially when developed in the framework of well-established mtDNA phylogeographical scenarios. In the Iberian Peninsula, the endemic Schreiber's green lizard Lacerta schreiberi exhibits two highly divergent and allopatric mtDNA lineages that started to split during the late Pliocene. In this work, we performed a fine-scale analysis of the putative mtDNA contact zone together with a global analysis of the patterns of variation observed at the nuclear beta-fibrinogen intron 7 (beta-fibint7). Using a combination of DNA sequencing with single-strand conformational polymorphism (SSCP) analysis, we show that the observed genealogy at the beta-fibint7 locus reveals extensive admixture between two formerly isolated lizard populations while the two mtDNA lineages remain essentially allopatric. In addition, a private beta-fibint7 haplotype detected in the single population where both mtDNA lineages were found in sympatry is probably the result of intragenic recombination between the two more common and divergent beta-fibint7 haplotypes. Our results suggest that the progressive incorporation of nuclear genealogies in investigating the ancient demography and admixture dynamics of divergent genomes will be necessary to obtain a more comprehensive picture of the evolutionary history of organisms.     

线粒体DNA和人类进化

线粒体ＤＮＡ（ｍｔＤＮＡ）由于自身比较独特的遗传特性（母系遗传、缺乏重组和进化速率高）而被广泛地应用于人类群体的起源和演化研究。通过对其全序列的限制性酶切和Ｄ－环高变区序列数据的分析,ｍｔＤＮＡ较好地阐明了人类学中诸如现代人类起源、人群过去动态的估计以及单个人群的区域性微分化和人口历史学等问题。综述了近年来世界各人群ｍｔＤＮＡ的研究进展、研究方法的改进、ｍｔＤＮＡ与核基因标记结果的异同、ｍｔＤＮＡ     

Phylogeography of the false smooth snakes, Macroprotodon (Serpentes, Colubridae): mitochondrial DNA sequences show European populations arrived recently from Northwest Africa

Mitochondrial DNA (1075 bp: cytochrome b, 300 bp; 12S rRNA, 393 bp; and 16S rRNA, 382 bp) corroborates the monophyly of the genus Macroprotodon and of the species M. mauritanicus, M. abubakeri, and M. brevis. The subspecies M. brevis ibericus is also monophyletic. The mtDNA tree presented here indicates that M. cucullatus consists of at least two separate units and may possibly represent a primitive morphology rather than a species in its own right. However, this hypothesis is tentative since it is only reflects the history of a single evolutionary unit (mtDNA). A definitive understanding of the evolution of M. cucullatus will not be possible until informative nuclear markers are added to the mitochondrial data. Macroprotodon appears to have originated in the Maghreb region of NW Africa and speciated there around 4-5.5 million years ago around the end of the Miocene period, after which its three main lineages may each have expanded north into more mesic conditions. The group also spread eastwards into coastal areas of Libya quite recently and on to Egypt and Israel. Later still, M. b. ibericus from extreme north Morocco reached the Iberian Peninsula, and M. mauritanicus from Tunisia or Algeria colonised the Balearic Islands of Menorca and Mallorca. Both these range extensions may result from very recent natural colonisations or even from accidental human introduction. Recency of origin of Iberian and Balearic populations is indicated by uniformity of their mtDNA even across large distances, and its great similarity to that of populations in source regions. Isolated populations assigned to M. cucullatus in the Hoggar mountains (southern Algeria) and Western Sahara are probably relicts from quite recent periods of climatic amelioration in the North African desert.     

The incomplete natural history of mitochondria

Mitochondrial DNA (mtDNA) has been used to study molecular ecology and phylogeography for 25 years. Much important information has been gained in this way, but it is time to reflect on the biology of the mitochondrion itself and consider opportunities for evolutionary studies of the organelle itself and its ecology, biochemistry and physiology. This review has four sections. First, we review aspects of the natural history of mitochondria and their DNA to show that it is a unique molecule with specific characteristics that differ from nuclear DNA. We do not attempt to cover the plethora of differences between mitochondrial and nuclear DNA; rather we spotlight differences that can cause significant bias when inferring demographic properties of populations and/or the evolutionary history of species. We focus on recombination, effective population size and mutation rate. Second, we explore some of the difficulties in interpreting phylogeographical data from mtDNA data alone and suggest a broader use of multiple nuclear markers. We argue that mtDNA is not a sufficient marker for phylogeographical studies if the focus of the investigation is the species and not the organelle. We focus on the potential bias caused by introgression. Third, we show that it is not safe to assume a priori that mtDNA evolves as a strictly neutral marker because both direct and indirect selection influence mitochondria. We outline some of the statistical tests of neutrality that can, and should, be applied to mtDNA sequence data prior to making any global statements concerning the history of the organism. We conclude with a critical examination of the neglected biology of mitochondria and point out several surprising gaps in the state of our knowledge about this important organelle. Here we limelight mitochondrial ecology, sexually antagonistic selection, life-history evolution including ageing and disease, and the evolution of mitochondrial inheritance.     

Genetic Polymorphism of 11 Allozyme Loci in Populations of Wall Lizards (Podarcis sp.) from the Iberian Peninsula and North Africa

The taxonomy of Iberian and North African wall lizards (Podarcis sp.) has been controversial. Recently, morphological and mtDNA sequence data have provided new information on differentiation within these lacertids. To compare these results to those provided by nuclear markers, we investigated variation at 11 polymorphic protein loci using conventional electrophoresis and isoelectric focusing in 11 populations belonging to seven different mtDNA lineages. A total of 62 alleles were found. Populations belonging to the same mtDNA type presented high genetic similarity, whereas strong differentiation was observed between groups. These results are consistent with those previously obtained from morphological and mtDNA analysis and support the idea that Iberian and North African Podarcis are composed of several well-differentiated entities, some of which are already recognized as species, whereas others (belonging to the P. hispanica complex) clearly need taxonomic revision.