Genetic admixture between the Iberian endemic lizards Podarcis bocagei and Podarcis carbonelli: evidence for limited natural hybridization and a bimodal hybrid zone

When recently diverged taxa come into contact, the extent of introgression between them is related to the degree of differentiation that they have achieved. Studying contact zones is therefore essential to understand if differentiated taxa are reproductively isolated and, ultimately, if they are likely to remain distinct. Recent work on Iberian and North African wall lizards ( Podarcis ) has documented the existence of multiple evolutionary units, diagnosable both by genetic markers and morphology, but suggests that gene flow between distinct forms has occurred. Therefore, we were interested in evaluating how species boundaries are maintained in the areas where they meet. In this work, we study the contact zone between Podarcis bocagei and Podarcis carbonelli . We sampled a transect including the only locality where these two lizards are known to occur in syntopy and analysed a battery of 15 unlinked nuclear genetic markers and mitochondrial DNA. We also conducted a preliminary analysis of morphology and fertility. Using model-based clustering approaches, we show that the two species hybridize in the population where they have direct contact, but evidences of introgression are low for nearby populations. Although a significant number of individuals show evidence of admixture, this hybrid zone is clearly bimodal, suggesting strong barriers to gene flow, of which the putative nature are discussed. Interestingly, morphological analyses do not support the existence of intermediate forms among individuals that are admixed genetically. Taken together, these results constitute further evidence validating P. bocagei and P. carbonelli as distinct species.     

Contrasting patterns of population subdivision and historical demography in three western Mediterranean lizard species inferred from mitochondrial DNA variation

Pleistocene climatic oscillations were a major force shaping genetic variability in many taxa. We analyse the relative effects of the ice ages across a latitudinal gradient in the Western Mediterranean region, testing two main predictions: (i) species with historical distributions in northern latitudes should have experienced greater loss of suitable habitat, resulting in higher extinction of historical lineages than species distributed in southern latitudes, where the effects of the ice ages were not as drastic. This would be reflected in the observation of lower diversity and number of differentiated lineages in northern areas. (ii) a signature of demographic expansion following the climate amelioration should be obvious in northern species, whereas in the south evidence of long-term effective population size stability should be observed. We used as models three species of wall lizards (Podarcis bocagei, Podarcis carbonelli and Podarcis vaucheri) that replace each other along the study area. We investigated the patterns of mitochondrial DNA diversity and subdivision and obtained demographic parameter estimates for each species. Our results suggest that P. bocagei, the northernmost species, bears low genetic diversity, a shallow coalescent history and marks of a demographic expansion. In contrast, P. vaucheri, the species with a southernmost distribution, shows deeper coalescence events, complex geographical substructure and no evidence for population growth. The species with an intermediate distribution, P. carbonelli, shows average levels of diversity, substructure and population growth. Taken together, these results conform to our main predictions and are explained by a differential influence of the ice ages on distinct latitudes.     

Phylogeography of the Italian wall lizard, Podarcis sicula, as revealed by mitochondrial DNA sequences

In a phylogeographical survey of the Italian wall lizard, Podarcis sicula, DNA sequence variation along an 887-bp segment of the cytochrome b gene was examined in 96 specimens from 86 localities covering the distribution range of the species. In addition, parts of the 12S rRNA and 16S rRNA genes from 12 selected specimens as representatives of more divergent cytochrome b haploclades were sequenced (together about 950 bp). Six phylogeographical main groups were found, three representing samples of the nominate subspecies Podarcis sicula sicula and closely related subspecies and the other three comprising Podarcis sicula campestris as well as all subspecies described from northern and eastern Adriatic islands. In southern Italy a population group with morphological characters of P. s. sicula but with the mitochondrial DNA features of P. s. campestris was detected indicating a probably recent hybridization zone. The present distribution patterns were interpreted as the consequence of natural events like retreats to glacial refuges and postglacial area expansions, but also as the results of multiple introductions by man.     

Phylogenetic relationships and biogeography of Podarcis species from the Balkan Peninsula, by bayesian and maximum likelihood analyses of mitochondrial DNA sequences

Wall lizards of the genus Podarcis (Sauria, Lacertidae) comprise 17 currently recognized species in southern Europe, where they are the predominant nonavian reptile group. The taxonomy of Podarcis is complex and unstable. Based on DNA sequence data, the species of Podarcis falls into four main groups that have substantial geographic coherence (Western island group, southwestern group, Italian group, and Balkan Peninsula group). The Balkan Peninsula species are divided into two subgroups: the subgroup of P. taurica (P. taurica, P. milensis, P. gaigeae, and perhaps P. melisellensis), and the subgroup of P. erhardii (P. erhardii and P. peloponnesiaca). In the present study, the question of phylogenetic relationships among the species of Podarcis encountered in the Balkan Peninsula was addressed using partial mtDNA sequences for cytochrome b (cyt b) and 16S rRNA (16S). The data support the monophyly of Podarcis and suggest that there are three phylogenetic clades: the clade A (P. taurica, P. gaigeae, P. milensis, and P. melisellensis); the clade B (P. erhardii and P. peloponnesiaca), and the clade C (P. muralis and P. sicula). By examining intraspecific relationships it was found that extant populations of P. erhardii are paraphyletic. Furthermore, subspecies previously defined on the basis of morphological characteristics do not correspond to different molecular phylogenetic clades, suggesting that their status should be reconsidered. The distinct geographic distribution of the major clades of the phylogenetic tree and its topology suggest a spatial and temporal sequence of phylogenetic separations that coincide with some major paleogeographic separations during the geological history of the Aegean Sea. The results stress the need for a reconsideration of the evolutionary history of Balkan Podarcis species and help overcome difficulties that classical taxonomy has encountered at both the species and subspecies level.     

Phylogeography of the Marsh Fritillary Euphydryas aurinia (Lepidoptera: Nymphalidae) in the UK

We sampled extant (and extinct) populations of Euphydryas aurinia to examine the phylogeography of the species in the UK. We were interested in whether the genetic structuring of populations reflects anything other than a single recent post-glacial colonization event. Four hundred base pairs of the mitochondrial Cytochrome b gene were sequenced from individuals from populations throughout the UK, as well as populations in France and Portugal and seven polymorphic allozyme loci resolved. The mean number of allozyme alleles per locus was 4 and isolation by distance was shown not to be a factor in the geographic structuring of genetic diversity either with or without the inclusion of the French data (Mantel statistic Z = 0.015 and 0.112 respectively, P >0.5). Cytochrome b nucleotide diversity (average number of nucleotide differences per site between two sequences) was low overall (0.003, n = 63) but mean cytochrome b gene diversity over all populations was 0.77. The presence of the Portuguese and French haplotype in Scottish populations indicates that the Iberian peninsula was likely to be one glacial refugium for E. aurinia populations. The pattern of mitochondrial DNA found in the UK could be interpreted in one of two ways: (1) two separate colonization events or (2) a single slow colonization event. Allele frequency distributions followed a similar geographical pattern as mitochondrial DNA haplotypes. An AMOVA assigned just 2.68% of allozyme genetic variation to the grouping tested, providing more support for the single colonization event theory.     

Comparing patterns of nuclear and mitochondrial divergence in a cryptic species complex: the case of Iberian and North African wall lizards (Podarcis, Lacertidae)

Combining different sources of information is essential for a complete understanding of the process of genetic differentiation between species. The Iberian and North African wall lizard ( Podarcis ) species complex has been the object of several studies regarding morphological and mitochondrial DNA variation but, so far, no large-scale survey of nuclear variation within this group has been accomplished. In this study, ten polymorphic allozyme loci were studied in 569 individuals collected across the Iberian Peninsula and North Africa. The obtained data were analysed using both conventional population genetic tools and recent Bayesian model-based clustering methods. Our results show that there are several well-differentiated entities corroborating the major splits observed in mtDNA analyses. These groups correspond not only to the fully recognized species Podarcis bocagei , Podarcis carbonelli , and Podarcis vaucheri but also to multiple forms within the polytypic Podarcis hispanica , all of which have a similar level of differentiation to that observed between the acknowledged species. However, relationships between forms are weakly supported both by population and individual clustering methods, suggesting a scenario of a rapid diversification that contrasts to the clear bifurcating model assumed from previous mtDNA analyses. Individual multilocus analyses report few individuals misassigned or apparently admixed, some of which are most likely explained by the persistence of high levels of ancestral polymorphism. Other admixed individuals, however, are probably the result of limited levels of gene flow between forms.  © 2007 The Linnean Society of London, Biological Journal of the Linnean Society , 2007, 91 , 121–133.     

Contrasting phylogeographical patterns between mainland and island taxa of the Pinus luchuensis complex

Species whose geographical distribution encompasses both mainland and island populations provide an ideal system for examining isolation and genetic divergence. In this study, paternally transmitted chloroplast DNA (cpDNA) and maternally transmitted mitochondrial DNA (mtDNA) were used to estimate population structure and phylogeography of Pinus luchuensis, a species found in eastern China (ssp. hwangshanensis), Taiwan (ssp. taiwanensis), and the Ryukyu Archipelago (ssp. luchuensis). Gene genealogies of both mtDNA and cpDNA reveal two major lineages. Molecular dating indicates that these lineages diverged before the colonization of P. luchuensis subspecies in Taiwan and the Ryukyu Archipelago. Both mtDNA and cpDNA show a lack of correspondence between molecular phylogeny and subspecies designation. Phylogeographical analysis suggests that paraphyly of the subspecies is the result of recent divergence rather than secondary contacts. In spite of the short divergence history of P. luchuensis on islands, the island populations show the same degree of genetic divergence as mainland populations. Low levels of genetic diversity in the mainland ssp. hwangshanensis suggest demographic bottlenecks. In contrast, the high heterogeneity of genetic composition for island populations is likely to be associated with a history of multiple colonization from the mainland. The spatial apportionment of organelle DNA polymorphisms is consistent with a pattern of stepwise colonization on island populations.     

Vicariance divergence and gene flow among islet populations of an endemic lizard

Allopatry and allopatric speciation can arise through two different mechanisms: vicariance or colonization through dispersal. Distinguishing between these different allopatric mechanisms is difficult and one of the major challenges in biogeographical research. Here, we address whether allopatric isolation in an endemic island lizard is the result of vicariance or dispersal. We estimated the amount and direction of gene flow during the divergence of isolated islet populations and subspecies of the endemic Skyros wall lizard Podarcis gaigeae, a phenotypically variable species that inhabits a major island and small islets in the Greek archipelago. We applied isolation-with-migration models to estimate population divergence times, population sizes and gene flow between islet-mainland population pairs. Divergence times were significantly correlated with independently estimated geological divergence times. This correlation strongly supports a vicariance scenario where islet populations have sequentially become isolated from the major island. We did not find evidence for significant gene flow within P. g. gaigeae. However, gene-flow estimates from the islet to the mainland populations were positively affected by islet area and negatively by distance between the islet and mainland. We also found evidence for gene flow from one subspecies (P. g. weigandi) into another (P. g. gaigeae), but not in the other direction. Ongoing gene flow between the subspecies suggests that even in this geographically allopatric scenario with the sea posing a strong barrier to dispersal, divergence with some gene flow is still feasible.     

Biometrical analysis reveals major differences between the two subspecies of the European rabbit

The climatic oscillations that have occurred in the last few million years have strongly affected species distribution ranges. Highly divergent genetic lineages arose, some of which correspond to recognized subspecies that currently occupy small geographical areas. Understanding the implications of the genetic differences between these subspecies is crucial for proper conservation of Evolutionarily Significant Units. We use the two European rabbit subspecies, Oryctolagus cuniculus cuniculus and O. c. algirus, in the Iberian Peninsula as a model to investigate the repercussions at the biometric level of their largely recognized genetic differentiation. To accomplish this we analysed the ear and hind foot length, and the body mass of 999 adult rabbits from 27 locations across the distribution range of both subspecies in their native range, the Iberian Peninsula. Our results show biometric differences between the two subspecies, also explained by geographical location and sex, O. c. algirus being lighter and having shorter ear and hind foot lengths. We examine these findings under an evolutionary framework, and discuss their implications for current conservation efforts. Future research should focus on the ecological implications of these biometric differences, namely potential different habitat use and anti‐predatory strategies in the species' native range.     

Molecular phylogeny and biogeography of the wall-lizard Podarcis erhardii (Squamata: Lacertidae)

Erhard's wall lizard, Podarcis erhardii (Sauria: Lacertidae), is highly diversified in Greece and especially in the southern Aegean region. Out of the 28 recognized subspecies, 27 are found in Greece from the North Sporades island-complex in the North Aegean (grossly south of the 39th parallel) to the island of Crete in the South. The species exhibits great morphological and ecological plasticity and inhabits many different habitats from rocky islets and sandy shores to mountaintops as high as 2000m. By examining intraspecific variability at a segment of the mitochondrial gene cytochrome b we have found that that extant populations of P. erhardii are paraphyletic. Furthermore, we have found that subspecies previously defined on the basis of morphological characteristics do not correspond to different molecular phylogenetic clades, so that their status should be reconsidered. The DNA based biogeographical and phylogenetic history of Podarcis in Southern Greece is congruent with available paleogeographic data of the region, which supports the view that DNA sequences may be a useful tool for the study of palaeogeography.     

Genetic diversity of flat-headed vole (Alticola strelzowi (Kastschenko, 1899)) inferred from cytochrome b variation

Variation of the cytochrome b gene fragment was examined in 27 flat-headed voles Alticola strelzowi from different parts of the species range. A total of 15 haplotypes were described, while the species as a whole was characterized by low levels of genetic differentiation and polymorphism. The haplotypes fell into three haplogroups, one of which corresponded to the subspecies A. s. strelzowi, and the other two, to A. s. desertorum. Based on different index values, the level of genetic polymorphism in the later subspecies was considered to be higher than in the first one. Phylogeographic analysis suggested post-glacial dispersal of flat-headed voles from a single refugium located in Western Altai. Using different techniques, relatively recent colonization of the Central Altai territory was demonstrated (subspecies A. s. strelzowi), which determined low level of genetic variation in this territory.     

Differences in chemical signals may explain species recognition between an island lizard, Podarcis atrata, and related mainland lizards, P. hispanica

Chemical signals can be the basis of interspecific recognition and speciation in many animals. The Columbretes Islands wall lizard, Podarcis atrata is very close genetically to the mainland Iberian wall lizard Podarcis hispanica. However, a previous study suggested that chemosensory interspecific recognition would avoid reproductive interactions and hybridization between these two species. These results suggested that chemicals used in intraspecific communication might differ in composition and/or proportions between these two species. In this paper, we used gas chromatography–mass spectrometry (GC–MS) to characterize the chemical composition of the lipophilic fraction from femoral gland secretions of male P. atrata and P. hispanica. The analysis showed that chemicals found in femoral secretions varied in composition and proportions between species and between populations. Seven steroids and two unidentified waxy esters, were exclusive of P. atrata lizards from the islands. In contrast, nine steroids and other six compounds were only found in mainland P. hispanica. There were also differences in proportions of shared compounds between species. Moreover, all these differences were higher between P. atrata and P. hispanica than between any population of P. hispanica. Chemical differences might be consequence of genetic differences, but they could also be explained by adaptation to different habitats with different climatic conditions or diet resources. Compounds that are specific of each species, or differences in the pattern of compounds, could explain species recognition. Therefore, these results of chemical composition and previous studies of chemosensory recognition reinforce the fact that the genetic differences between P. hispanica and P. atrata may result in an effective reproductive isolation between these two taxa.     

Strong genetic differentiation due to multiple founder events during a recent range expansion of an introduced wall lizard population

Biological invasions represent ideal systems for the study of evolutionary processes associated with colonization events. It has been hypothesized that the genetic diversity is generally decreasing from the centre of the range to the margins due to multiple founder events. Invasive populations offer the opportunity to test this hypothesis at a fine spatial and temporal scale. We analysed the genetic structure of a large expanding non-native population of the Common Wall Lizard (Podarcis muralis) in Passau (Germany) using thirteen microsatellite loci. We analyzed the genetic structure and levels of admixture across a transect reflecting the expansion process and tested for a loss of genetic diversity and an increase of genetic differentiation from the centre to the invasion front. Our results demonstrate that significant genetic population structure can emerge rapidly at a small spatial scale. We found a trend for an increase in genetic differentiation and a decrease in genetic diversity from the invasion centre to the expanding range margin, suggesting that genetic drift is the major factor causing this pattern. The correlation between genetic diversity and average genetic differentiation was significant among sites. We hypothesize that the territoriality of P. muralis generates sufficient rates of noncontiguous and stratified dispersal from longer established sites to maintain significant genetic diversity at the invasion front. Simultaneously, territoriality might restrict the colonization success of migrants at established sites, so that in combination with founder events a strong differentiation arises.     

Colonization and divergence: phylogeography and population genetics of the Atlantic coast beach mice

Barrier island taxa provide an opportunity to investigate recent evolutionary processes, such as colonization and isolation of recently diverged taxa, and provide important insights into understanding contemporary diversity and the assessment of conservation units. Using rapidly evolving genetic markers (mitochondrial DNA and microsatellites), we studied the Atlantic coast beach mouse subspecies (Peromyscus polionotus decoloratus, P. p. niveiventris, and P. p. phasma). Our data indicate that each of the extant coastal subspecies (P. p. niveiventris and P. p. phasma) is comprised of unique haplotypes indicative of their isolation, while the extinct subspecies, P. p. decoloratus, contain a single haplotype, which was shared with P. p. phasma. Moreover, all the coastal haplotypes originate from a single mainland haplotype found in central Florida, USA. The microsatellite data indicated high levels of genetic structure among our sampled populations. Additionally, these data group the populations into three distinct genetic clusters, with each of the extant coastal subspecies belonging to their own cluster and the mainland individuals forming a separate cluster. The extant Atlantic coast beach mice are on separate evolutionary trajectories, thus representative of separate taxonomic units. Therefore, the data support that two extant subspecies on the Florida Atlantic coast fit the Distinct Population Segment designation and should be managed and conserved as two separate independent units.     

Phylogeography of Balkan wall lizard (Podarcis taurica) and its relatives inferred from mitochondrial DNA sequences

Wall lizards of the genus Podarcis (Sauria, Lacertidae) comprise 17 currently recognized species in southern Europe, where they are the predominant reptile group. The taxonomy of Podarcis is complex and unstable. Based on DNA sequence data the species of Podarcis falls into four main groups that have substantial geographical conherence (western island group, southwestern group, Italian group and Balkan group). The Balkan species are divided in two subgroups: the subgroup of Podarcis taurica (P. taurica, P. milensis, P. gaigeae and perhaps P. melisellensis), and the subgroup of Podarcis erhardii (P. erhardii and P. peloponnesiaca). We addressed the question of phylogenetic relations among the species of the P. taurica subgroup encountered in Greece, as they can be inferred from partial mtDNA (cyt b and 16S) sequences. Our data support the monophyly of P. taurica subgroup and suggest that P. gaigeae, P. milensis and P. melisellensis form a clade, which thereinafter connects to P. taurica. Within the previous clade, P. gaigeae is more closely related to P. milensis than to P. melisellensis. However, the specimens of P. taurica were subdivided in two different groups. The first one includes the specimens from northeastern Greece, and the other group includes the specimens from the rest of continental Greece and Ionian islands. Because the molecular clock of the cyt b and 16 rRNA genes was not rejected in our model test, it is possible to estimate times of speciation events. Based on the splitting of the island of Crete from Peloponnisos [c. 5 million years ago (Ma)], the evolutionary rate for the cyt b is 1.55% per million years (Myr) and for the 16S rRNA is 0.46% per Myr. These results suggest that the evolutionary history of P. taurica in Greece is more complex than a single evolutionary invasion. The data analysed, stress the need for a reconsideration of the evolutionary history of Greek Podarcis species and help overcome difficulties that classical taxonomy has encountered at both the species level.     

新疆草兔的种群遗传结构和亚种分化

新疆草兔 (Lepus capensis) 的群体遗传结构至今无系统的研究报道,亚种水平的分类也长期存在争议.该文测定了形态分类上的新疆草兔3个亚种共87个个体的线粒体DNA (mtDNA)控制区(control region,D-Loop)592 bp的序列,经分析发现148个多态性位点,共定义了44个单倍型.新疆草兔的单倍型多样度(h,0.977 ± 0.005)和核苷酸多样度(π,0.064 ± 0.031)都较高,显示了较高的遗传多样性.分子变异分析(AMOVA)结果显示,4个地理群体间的显著分化可能是由地理隔离造成的.群体遗传结构分析显示,新疆草兔包含4个进化枝,并且每个进化枝都对应特定的分布区域,显示了明显的系统地理结构.该研究的结果支持形态分类上草兔西域亚种(L.c.lehmanni)的分类地位; 但中亚亚种(L.c.centrasiaticus)被分为两个独立的进化枝,提示可能存在两个亚种; 帕米尔亚种(L.c.pamirensis)与其他亚种间的遗传距离在13%以上,提示其可能已达到种的分化水平.     

Genetic structure of kestrel populations and colonization of the Cape Verde archipelago

Genetic diversity and population structure were studied in eight populations of the kestrel Falco tinnunculus to identify the genetic consequences of spatial distribution and to infer the colonization patterns of the Cape Verde archipelago. We studied genetic differentiation and gene flow among seven island populations and one mainland population using nine microsatellite loci. Within the archipelago, differentiation was strong and genetic diversity and heterozygosity were low but variable among populations. Two subspecies F. tinnunculus neglectus on the northwestern islands and F. tinnunculus alexandri on all the other islands were identified as genetically distinct units. F. t. alexandri could be further separated into two groups on eastern and southern islands. Populations are probably founded by birds originating from the mainland. Immigration is more likely to the eastern and southern populations, whereas the northwestern islands with the lowest genetic diversity and highest differentiation are likely to exhibit fewer founding events by immigrants. The number of founding events on each island may depend not only on geographical distance to neighbouring populations, but also on directional immigration due to the northeastern trade winds. This may explain differences in genetic differentiation and diversity between populations and subspecies and may enable allopatric speciation.     

Variation in morphology, gait characteristics and speed of locomotion in two populations of lizards

Locomotor behaviour varies between two subspecies of the Spanish wall lizard Podarcis hispanica. One subspecies inhabits the Columbretes islands, the other lives on the Spanish mainland. Size standardized voluntary speeds (as measured in unrestrained laboratory conditions) are lower in the island population (P. h. atrata) than in the mainland population (P. h. hispanica). Maximal running performance (when chased) is much higher in the mainland population than in the island population. High speed video recordings show that subspecies differ in gait characteristics: individuals from the mainland modulate running velocity primarily by modifying stride length, individuals from the island primarily by altering stride frequency. P. h. hispanica's strategy for modulating speed probably allows this mainland subspecies to attain higher maximal speeds than the island subspecies P. h. atrata. Theoretical considerations suggest that at high speeds, P. h. hispanica's running style is energetically more favourable, but this hypothesis awaits experimental verification. We suggest that the differences in locomotion efficiency between the subspecies result from differences in predation pressure between the mainland and the island. The mainland study site has a higher predator diversity and offers less hiding opportunities to the lizards.     

Head shape allometry and proximate causes of head sexual dimorphism in Podarcis lizards: joining linear and geometric morphometrics

Podarcis bocagei and Podarcis carbonelli are two lacertid species endemic to the western Iberian Peninsula, and both show head size and shape sexual dimorphism. We studied immature and adult head sexual dimorphism and analysed ontogenetic trajectories of head traits with body and head size, aiming to shed light on the proximate mechanisms involved. Immatures were much less dimorphic than adults, but geometric morphometric techniques revealed that head shape sexual differences are already present at this stage. Males and females differed in allometry of all head characters with body size, with males showing a disproportionate increase of head size and dimensions. On the other hand, head dimensions and head shape changed with increasing head size following similar trends in both sexes, possibly indicating developmental restrictions. Consequently, adult sexual dimorphism for head characters in these species is the result of both shape differences in the immature stage and hypermetric growth of the head in relation to body size in males.  © 2008 The Linnean Society of London, Biological Journal of the Linnean Society , 2008, 93 , 111–124.