Phylogenetic value of jaw elements of lacertid lizards (Squamata: Lacertoidea): a case study with Oligocene material from France

Several extinct species are known from the family Lacertidae, but due to poor preservation, many of them are based on single bones. Here, we compare phylogenetic signals of disarticulated premaxillae, maxillae and dentaries of lacertids from four French Oligocene localities (Coderet, La Colombière, Roqueprune 2, Mas de Got B). We identified five morphotypes among the premaxillae, six among the maxillae, and ten among the dentaries. These morphotypes were scored as individual taxa per locality into three separate character matrices with the same 246 characters, one matrix for each jaw element. Subsequently, the phylogenetic position of the morphotypes was tested using maximum parsimony. The consensus trees with the dentaries and the maxillae found a large polytomy including all taxa except the outgroup taxon Gekko gecko. The consensus tree with the premaxillae showed a considerably more resolved topology but found all morphotype taxa outside Lacertidae. In a second step, we compared the constitution of our three datasets and the morphotype taxa. Our results suggest that a combination of convergent characters and missing data led to the outgroup position of the premaxilla morphotype taxa. The poor resolution of the maxillae strict consensus is likely a consequence of their fragmentary preservation. For the dentaries, a high amount of missing data due to the high number of morphotype taxa most likely caused the poor tree resolution. Indeed, tests with fewer morphotypes found tree resolutions comparable to the premaxilla data. When linking the morphotypes, five possible lacertid "species" were found. Comparison with already known French Oligocene lacertid species points to a slightly higher species richness of Lacertidae at that time than known before. Reliable species classification based on phylogeny only seems possible when combining the jaw elements or in association with other cranial and postcranial material, putting some doubt on species identifications based on single bones.

Fossil-calibrated time tree of Podarcis wall lizards provides limited support for biogeographic calibration models

Podarcis wall lizards are endemic to the Mediterranean Basin where they represent the predominant reptile group. Despite being extensively used as model organisms in evolutionary and ecological studies their phylogeny and historical biogeography are still incompletely resolved. Moreover, molecular clock calibrations used in wall lizard phylogeography are based on the assumption of vicariant speciation triggered by the abrupt Mediterranean Sea level rise at the end of the Messinian salinity crisis (MSC). However, the validity of this biogeographic calibration remains untested. In this study we inferred a robust time tree based on multilocus data and fossil calibrations using both gene concatenation and species-tree approaches and including models with gene-flow. We found five deeply divergent, geographically coherent, and well-supported clades comprising species from i) Iberian Peninsula and North Africa; ii) Western Mediterranean islands, iii) Sicilian and Maltese islands; and iv-v) Balkan region and Aegean islands. The mitochondrial tree shows some inconsistencies with the species tree that warrant future investigation. Diversification of main clades is estimated in a short time frame during the Middle Miocene and might have been associated with a period of global climate cooling with the establishment of a marked climatic zonation in Europe. Cladogenetic events within the main clades are scattered throughout the time tree, from the Late Miocene to the Early Pleistocene, suggesting that speciation events in wall lizards reflect a complex interplay between regional topography, climate and geological history rather than a shared major climatic or paleogeographic event. Our absolute time estimates, as well as a relative dating approach, demonstrate that the assumption of a causal link between sea-level rise at the end of the MSC and the diversification of many island endemics is not justified. This study reinforces the notion that multiple dispersal and vicariant events, at different time frames, are required to explain current allopatric distributions and to account for the historical assembly of Mediterranean biota, and cautions against the use of biogeographic calibrations based on the assumption of vicariance.

Ancient divergence, a crisis of salt and another of ice shaped the evolution of the west Mediterranean butterfly Euchloe tagis

The Mediterranean region is an extremely complex hotspot where, since the Miocene, extensive geological, habitat and climate changes have taken place, alternating between warm and cold periods. These phenomena have taken a toll on the genetic composition of species, and surviving lineages have often adapted locally and diverged to the point of (complete) speciation. To study these phenomena, in this study we used one of the most enigmatic butterflies, the Portuguese dappled white, Euchloe tagis, a west Mediterranean endemic with fragmented, morphologically differentiated populations whose status have long been disputed. Even its affiliations with other Anthocharidini are largely unresolved. We used mitochondrial and nuclear markers under a phylogenetic and phylogeographical framework to evaluate its placement among relatives and population differentiation, reconstructing its evolutionary history. We found that this species had a Miocene origin ~15 Mya and was nearest to Euchloe s.s. and Elphinstonia. Its populations showed high genetic diversity but all coalesced to 5.3 Mya, when European and all but one African population diverged. Our multiple findings concerning the evolution of E. tagis through a changing, narrow habitable area might provide a more general perspective on how species survive within this hotspot of paramount importance.

Dispersal and local persistence shape the genetic structure of a widespread Neotropical plant species with a patchy distribution

BACKGROUND AND AIMS

Isolated populations constitute an ideal laboratory to study the consequences of intraspecific divergence, because intrinsic incompatibilities are more likely to accumulate under reduced gene flow. Here, we use a widespread bromeliad with a patchy distribution, Pitcairnia lanuginosa, as a model to infer processes driving Neotropical diversification and, thus, to improve our understanding of the origin and evolutionary dynamics of biodiversity in this highly speciose region.

METHODS

We assessed the timing of lineage divergence, genetic structural patterns and historical demography of P. lanuginosa, based on microsatellites, and plastid and nuclear sequence data sets using coalescent analyses and an Approximate Bayesian Computation framework. Additionally, we used species distribution models (SDMs) to independently estimate potential changes in habitat suitability.

KEY RESULTS

Despite morphological uniformity, plastid and nuclear DNA data revealed two distinct P. lanuginosa lineages that probably diverged through dispersal from the Cerrado to the Central Andean Yungas, following the final uplift of the Andes, and passed through long-term isolation with no evidence of migration. Microsatellite data indicate low genetic diversity and high levels of inbreeding within populations, and restricted gene flow among populations, which are likely to be a consequence of bottlenecks (or founder events), and high selfing rates promoting population persistence in isolation. SDMs showed a slight expansion of the suitable range for P. lanuginosa lineages during the Last Glacial Maximum, although molecular data revealed a signature of older divergence. Pleistocene climatic oscillations thus seem to have played only a minor role in the diversification of P. lanuginosa, which probably persisted through adverse conditions in riparian forests.

CONCLUSIONS

Our results imply drift as a major force shaping the evolution of P. lanuginosa, and suggest that dispersal events have a prominent role in connecting Neotropical open and forest biomes.

Spotlight on islands: on the origin and diversification of an ancient lineage of the Italian wall lizard Podarcis siculus in the western Pontine Islands

Groups of proximate continental islands may conceal more tangled phylogeographic patterns than oceanic archipelagos as a consequence of repeated sea level changes, which allow populations to experience gene flow during periods of low sea level stands and isolation by vicariant mechanisms during periods of high sea level stands. Here, we describe for the first time an ancient and diverging lineage of the Italian wall lizard Podarcis siculus from the western Pontine Islands. We used nuclear and mitochondrial DNA sequences of 156 individuals with the aim of unraveling their phylogenetic position, while microsatellite loci were used to test several a priori insular biogeographic models of migration with empirical data. Our results suggest that the western Pontine populations colonized the islands early during their Pliocene volcanic formation, while populations from the eastern Pontine Islands seem to have been introduced recently. The inter-island genetic makeup indicates an important role of historical migration, probably due to glacial land bridges connecting islands followed by a recent vicariant mechanism of isolation. Moreover, the most supported migration model predicted higher gene flow among islands which are geographically arranged in parallel. Considering the threatened status of small insular endemic populations, we suggest this new evolutionarily independent unit be given priority in conservation efforts.

The Alboran volcanic-arc modulated the Messinian faunal exchange and salinity crisis

What process triggered the Mediterranean Sea restriction remains debated since the discovery of the Messinian Salinity Crisis (MSC). Recent hypotheses infer that the MSC initiated after the closure of the Atlantic-Mediterranean Betic and Rifean corridors, being modulated through restriction at the Gibraltar Strait. These hypotheses however, do not integrate contemporaneous speciation patterns of the faunal exchange between Iberia and Africa and several geological features like the evaporite distribution. Exchange of terrestrial biota occurred before, during and after the MSC, and speciation models support an exchange path across the East Alborán basin (EAB) located a few hundreds of km east of the Gibraltar Strait. Yet, a structure explaining jointly geological and biological observations has remained undiscovered. We present new seismic data showing the velocity structure of a well-differentiated 14–17-km thick volcanic arc in the EAB. Isostatic considerations support that the arc-crust buoyancy created an archipelago leading to a filter bridge across the EAB. Sub-aerial erosional unconformities and onlap relationships support that the arc was active between ~10–6 Ma. Progressive arc build-up leading to an archipelago and its later subsidence can explain the extended exchange of terrestrial biota between Iberia and Africa (~7–3 Ma), and agrees with patterns of biota speciation and terrestrial fossil distribution before the MSC (10–6.2 Ma). In this scenario, the West Alboran Basin (WAB) could then be the long-postulated open-marine refuge for the Mediterranean taxa that repopulated the Mediterranean after the MSC, connected to the deep restricted Mediterranean basin through a sill at the Alboran volcanic arc archipelago.

Comparative phylogeography of amphibians and reptiles in Algeria suggests common causes for the east-west phylogeographic breaks in the Maghreb

A series of phylogeographic studies in the Maghreb identified a repeated pattern of deep genetic divergence between an eastern (Tunisia) and western (Morocco) lineage for several taxa but lack of sampling in Algeria made it difficult to know if the range limits between the eastern and western lineages were shared among taxa or not. To address this question, we designed a comparative phylogeographic study using 8 reptile and 3 amphibian species with wide distribution in the Maghreb as models. We selected species where previous studies had identified an East-West phylogeographic divide and collected sampled in Algeria to 1) examine whether the simple East-West divergence pattern still holds after filling the sampling gap in Algeria or if more complex diversity patterns emerge; 2) if the E-W pattern still holds, test whether the limits between the E and W clades are shared between species, suggesting that common historical process caused the E-W divergences; 3) if E-W limits are shared between species, use information on the age of the divergence to identify possible geological or climatic events that could have triggered these E-W differentiations. We found that the E-W pattern was generally maintained after additional sampling in Algeria and identified two common disjunction areas, one around the Algeria-Morocco border, the other one in Kabylia (central Algeria), suggesting that common historical mechanisms caused the E-W divergences in the Maghreb. Our estimates for the times to most common recent ancestors to the E and W clades span a wide range between the Messinian salinity crisis and the Plio-Pleistocene limit (except for one older split), suggesting different origins for the initial divergences and subsequent preservation of the E and W lineages in common climatic refugia in the west and the east of the Maghreb.

The role of ecological factors in determining phylogeographic and population genetic structure of two sympatric island skinks (Plestiodon kishinouyei and P. stimpsonii)

We conducted comparative phylogeographic and population genetic analyses of Plestiodon kishinouyei and P. stimpsonii, two sympatric skinks endemic to islands in the southern Ryukyus, to explore different factors that have influenced population structure. Previous phylogenetic studies using partial mitochondrial DNA indicate similar divergence times from their respective closest relatives, suggesting that differences in population structure are driven by intrinsic attributes of either species rather than the common set of extrinsic factors that both presumably have been exposed to throughout their history. In this study, analysis of mtDNA sequences and microsatellite polymorphism demonstrate contrasting patterns of phylogeography and population structure: P. kishinouyei exhibits a lower genetic variability and lower genetic differentiation among islands than P. stimpsonii, consistent with recent population expansion. However, historical demographic analyses indicate that the relatively high genetic uniformity in P. kishinouyei is not attributable to recent expansion. We detected significant isolation-by-distance patterns among P. kishinouyei populations on the land bridge islands, but not among P. stimpsonii populations occurring on those same islands. Our results suggest that P. kishinouyei populations have maintained gene flows across islands until recently, probably via ephemeral Quaternary land bridges. The lower genetic variability in P. kishinouyei may also indicate smaller effective population sizes on average than that of P. stimpsonii. We interpret these differences as a consequence of ecological divergence between the two species, primarily in trophic level and habitat preference.

Phylogeny, phylogeography, and evolution in the Mediterranean region: News from a freshwater mussel (Potomida, Unionida)

The Potomida genus (Bivalvia, Unionida) has a Circum-Mediterranean distribution and like other freshwater mussel species, its populations have suffered dramatic declines. Although this genus is currently considered as monotypic, it has a long history of taxonomic revisions and presently many aspects of its systematics and evolutionary history are unclear. We sampled a total of 323 individuals from 39 different sites across the Potomida genus distribution, and sequenced two mitochondrial (16S rDNA and Cytochrome c Oxidase Subunit I) and one nuclear (28S rDNA) genes to clarify its phylogeny and phylogeographic history. Our results show that the genus includes two well-supported clades, one comprising solely the western Mediterranean species Potomida littoralis, and the other including two eastern Mediterranean species, the Greek endemic P. acarnanica and the Anatolian and Middle Eastern P. semirugata. We suggest that Potomida started radiating during the upper Miocene, and that both vicariance and dispersal events shaped the diversification and distribution of the genus along the Mediterranean region. P. littoralis is further divided in two mitochondrial lineages, one restricted to Europe and the other occurring mostly in North Africa. Moreover, some European basins present both lineages in sympatry. The conservation status of the three recognized species should be reevaluated, particularly P. acarnanica, since it is restricted to two Greek river basins presenting a high risk of extinction. Overall, our results clarify some important gaps in knowledge concerning the phylogeny, phylogeography and evolution of the Potomida genus in the Mediterranean region with important taxonomical, ecological and conservational implications.

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.

Changing geographic ranges of ticks and tick-borne pathogens: drivers, mechanisms and consequences for pathogen diversity

The geographic ranges of ticks and tick-borne pathogens are changing due to global and local environmental (including climatic) changes. In this review we explore current knowledge of the drivers for changes in the ranges of ticks and tick-borne pathogen species and strains via effects on their basic reproduction number (R₀), and the mechanisms of dispersal that allow ticks and tick-borne pathogens to invade suitable environments. Using the expanding geographic distribution of the vectors and agent of Lyme disease as an example we then investigate what could be expected of the diversity of tick-borne pathogens during the process of range expansion, and compare this with what is currently being observed. Lastly we explore how historic population and range expansions and contractions could be reflected in the phylogeography of ticks and tick-borne pathogens seen in recent years, and conclude that combined study of currently changing tick and tick-borne pathogen ranges and diversity, with phylogeographic analysis, may help us better predict future patterns of invasion and diversity.

Phylogeography and diversification history of the day-gecko genus Phelsuma in the Seychelles islands

BACKGROUND

Lying in a shallow continental shelf cyclically affected by oscillating sea levels since the Miocene, the Seychelles islands are particularly interesting for evolutionary studies. Recent molecular studies are generating an emerging picture of the origin of its biota, yet very little is known regarding their phylogeographic structure or on the factors promoting diversification within the archipelago. Here we aimed to obtain a detailed depiction of the genetic structure and evolution of one of the most widespread vertebrate groups in the archipelago: the day-geckos of the genus Phelsuma. In parallel, we aimed to infer divergence times between species and subspecies, testing a long-standing hypothesis that argues for different time since sympatry between species as the cause of their different morphological differentiation across the archipelago.

RESULTS

Molecular data corroborated the existence of two main lineages, corresponding to the two currently recognized species. Divergences between species likely date back to the Mio-Pliocene, while more recent, Pleistocenic, divergences are suggested within each species. Populations from outer islands share mtDNA haplotypes with inner island populations, suggesting very recent dispersals (or introductions). We found no evidence of current gene flow between species, but results pointed to the possibility of gene flow between (now allopatric) subspecies. Time estimates suggest a synchronous divergence within each species (between island groups).

CONCLUSIONS

The geographic patterns of genetic variation agree with previous taxonomic subdivisions within each species and the origin of outer islands populations is clearly tracked. The similar intraspecific divergence time estimates obtained suggest that the differential body-size differentiation between species within each group of islands may be driven by factors other than character displacement proportional to time since sympatry, as previously suggested. These factors could include different habitats/resources available within each island group, niche differentiation and/or character displacement. We also bring again into consideration the hypothesis of body size being influenced by the distribution of native vegetation and social systems within this group, although it remains to be tested. Our results highlight not only the necessity of clarifying the role of ecology and interspecific interactions in this group's morphological diversification and community assemblage, but also the importance of co-evolutionary mechanisms and their importance for appropriate conservation of island biodiversity. Further, we provide a detailed description of the phylogeographic structure of these taxa across these islands, which still remain poorly characterized in this respect.

Challenges and pitfalls in the characterization of anonymous outlier AFLP markers in non-model species: lessons from an ocellated lizard genome scan

In the last few years, dozens of studies have documented the detection of loci influenced by selection from genome scans in a wide range of non-model species. Many of those studies used amplified fragment length polymorphism (AFLP) markers, which became popular for being easily applicable to any organism. However, because they are anonymous markers, AFLPs impose many challenges for their isolation and identification. Most recent AFLP genome scans used capillary electrophoresis (CE), which adds even more obstacles to the isolation of bands with a specific size for sequencing. These caveats might explain the extremely low number of studies that moved from the detection of outlier AFLP markers to their actual isolation and characterization. We document our efforts to characterize a set of outlier AFLP markers from a previous genome scan with CE in ocellated lizards (Lacerta lepida). Seven outliers were successfully isolated, cloned and sequenced. Their sequences are noncoding and show internal indels or polymorphic repetitive elements (microsatellites). Three outliers were converted into codominant markers by using specific internal primers to sequence and screen population variability from undigested DNA. Amplification in closely related lizard species was also achieved, revealing remarkable interspecific conservation in outlier loci sequences. We stress the importance of following up AFLP genome scans to validate selection signatures of outlier loci, but also report the main challenges and pitfalls that may be faced during the process.

Evolution and radiation in the scorpion Buthus elmoutaouakili Lourenco and Qi 2006 (Scorpiones: Buthidae) at the foothills of the Atlas Mountains (North Africa)

When low dispersal ability of an organism meets geographical barriers, the evolution of inter- and intraspecific differentiation is often facilitated. In the Atlas massif of North Africa, the genus Buthus splits into several species and diverges into numerous genetic lineages, often following the orographic structures of mountain systems. Such high mountain ranges often act as barriers for species with reduced mobility even on small spatial scales. To study the effect of orographic structures on organisms with low dispersal ability, we collected 61 individuals of the scorpion species Buthus elmoutaouakili at 18 locations around the southwestern foothills of the High Atlas and Antiatlas and in the Sousse valley (western Morocco). We analyzed intraspecific differentiation patterns within this geographically restricted area of about 100 × 50 km using 452 bp of the cytochrome oxidase I mitochondrial gene. We detected 5 distinct genetic lineages. In a second analysis, we added 61 previously published sequences from Buthus species from Europe and North Africa. Using a molecular clock approach, we detected old splits (4-5 Ma) separating the samples from 1) the western High Atlas and north of these mountains, 2) the Sousse valley and adjoining mountain areas, and 3) the southwestern Antiatlas. Further differentiation happened in the first 2 geographical groups about 3 Ma. Thus, the divergence time estimates based on a Bayesian approach support the onset of differentiation into these main clades along the Pliocene (5-2.3 Ma) when climatic oscillations started and a constant global cooling preceded the glacial-interglacial cycles of the Pleistocene. Further genetic splits into parapatric groups are detectable for the Sousse valley main group in the early Pleistocene. The climatic oscillations of the Pliocene and early Pleistocene might have caused repeated range shifts, expansions, and retractions leading to repeated vicariance, hereby producing the hierarchical structure of genetic differentiation in B. elmoutaouakili. A taxonomic revision, including morphological and molecular data, is needed to assess the status of each of these Buthus scorpion lineages.

A molecular phylogeny of the eastern group of ocellated lizard genus Timon (Sauria: Lacertidae) based on mitochondrial and nuclear DNA sequences

Timon, a small genus of lacertid lizards, includes four species distributed in two separate ranges in the western and eastern part of the Mediterranean Basin. Phylogenetic relationships between the two groups have not been resolved, and the taxonomic situation of the two subspecies of the eastern representative of the genus, Timon princeps, is not clear. To address these questions, partial DNA sequences of two nuclear (β-fibrinogen intron 7 and C-mos) and three mitochondrial (cytochrome b, 12S rRNA and 16S rRNA) genes were analyzed. Based on the high genetic distance between the two subspecies of T. princeps we promote their taxonomic status to full species, Timon princeps and Timon kurdistanicus. Divergence time estimates based on other lacertid species suggest that the separation of the green (Lacerta) and ocellated (Timon) lizards took place around 12 My ago, and that the Eastern group underwent speciation around 4-5 my ago, perhaps associated with the uplifting of the Zagros mountains. As expected given this ancient divergence and complex paleogeography, considerable levels of genetic diversity are recovered within both taxa, with geographically close individuals showing very divergent haplotypes.

Molecular phylogenetics and historical biogeography of the west-palearctic common toads (Bufo bufo species complex)

In most pan-Eurasiatic species complexes, two phenomena have been traditionally considered key processes of their cladogenesis and biogeography. First, it is hypothesized that the origin and development of the Central Asian Deserts generated a biogeographic barrier that fragmented past continuous distributions in Eastern and Western domains. Second, Pleistocene glaciations have been proposed as the main process driving the regional diversification within each of these domains. The European common toad and its closest relatives provide an interesting opportunity to examine the relative contributions of these paleogeographic and paleoclimatic events to the phylogeny and biogeography of a widespread Eurasiatic group. We investigate this issue by applying a multiproxy approach combining information from molecular phylogenies, a multiple correspondence analysis of allozyme data and species distribution models. Our study includes 304 specimens from 164 populations, covering most of the distributional range of the Bufo bufo species complex in the Western Palearctic. The phylogenies (ML and Bayesian analyses) were based on a total of 1988 bp of mitochondrial DNA encompassing three genes (tRNAval, 16S and ND1). A dataset with 173 species of the family Bufonidae was assembled to estimate the separation of the two pan-Eurasiatic species complexes of Bufo and to date the main biogeographic events within the Bufo bufo species complex. The allozyme study included sixteen protein systems, corresponding to 21 presumptive loci. Finally, the distribution models were based on maximum entropy. Our distribution models show that Eastern and Western species complexes are greatly isolated by the Central Asian Deserts, and our dating estimates place this divergence during the Middle Miocene, a moment in which different sources of evidence document a major upturn of the aridification rate of Central Asia. This climate-driven process likely separated the Eastern and Western species. At the level of the Western Palearctic, our dating estimates place most of the deepest phylogenetic structure before the Pleistocene, indicating that Pleistocene glaciations did not have a major role in splitting the major lineages. At a shallow level, the glacial dynamics contributed unevenly to the genetic structuring of populations, with a strong influence in the European-Caucasian populations, and a more relaxed effect in the Iberian 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.

Population genetics, taxonomy, phylogeny and evolution of Borrelia burgdorferi sensu lato

In order to understand the population structure and dynamics of bacterial microorganisms, typing systems that accurately reflect the phylogenetic and evolutionary relationship of the agents are required. Over the past 15 years multilocus sequence typing schemes have replaced single locus approaches, giving novel insights into phylogenetic and evolutionary relationships of many bacterial species and facilitating taxonomy. Since 2004, several schemes using multiple loci have been developed to better understand the taxonomy, phylogeny and evolution of Lyme borreliosis spirochetes and in this paper we have reviewed and summarized the progress that has been made for this important group of vector-borne zoonotic bacteria.

Association of Mc1r variants with ecologically relevant phenotypes in the European ocellated lizard, Lacerta lepida

A comprehensive knowledge on the genetic basis of coloration is crucial to understand how new colour phenotypes arise and how they contribute to the emergence of new species. Variation in melanocortin-1 receptor (Mc1r), a gene that has been reported as a target for repeated evolution in a wide range of vertebrate taxa, was assessed in European ocellated lizards (Lacerta lepida) to search for associations with melanin-based colour phenotypes. Lacerta lepida subspecies' distribution is associated with the three major bio-climatic regions in the Iberian Peninsula. A nonconserved and derived substitution (T162I) was associated with the L. l. nevadensis phenotype (prevalence of brown scales). Another substitution (S172C) was associated with the presence of black scales in both L. l. lepida and L. l. iberica, but no mutations were found to be associated with the higher proportion of black in L. l. iberica. Extensive genotyping of Mc1r along the contact zone between L. l. nevadensis and L. l. lepida revealed low gene flow (only two hybrids detected). The implications of these findings are discussed in the context of previous knowledge about the evolutionary history of ocellated lizards.

Phylogeography and demographic history of Lacerta lepida in the Iberian Peninsula: multiple refugia, range expansions and secondary contact zones

BACKGROUND

The Iberian Peninsula is recognized as an important refugial area for species survival and diversification during the climatic cycles of the Quaternary. Recent phylogeographic studies have revealed Iberia as a complex of multiple refugia. However, most of these studies have focused either on species with narrow distributions within the region or species groups that, although widely distributed, generally have a genetic structure that relates to pre-Quaternary cladogenetic events. In this study we undertake a detailed phylogeographic analysis of the lizard species, Lacerta lepida, whose distribution encompasses the entire Iberian Peninsula. We attempt to identify refugial areas, recolonization routes, zones of secondary contact and date demographic events within this species.

RESULTS

Results support the existence of 6 evolutionary lineages (phylogroups) with a strong association between genetic variation and geography, suggesting a history of allopatric divergence in different refugia. Diversification within phylogroups is concordant with the onset of the Pleistocene climatic oscillations. The southern regions of several phylogroups show a high incidence of ancestral alleles in contrast with high incidence of recently derived alleles in northern regions. All phylogroups show signs of recent demographic and spatial expansions. We have further identified several zones of secondary contact, with divergent mitochondrial haplotypes occurring in narrow zones of sympatry.

CONCLUSIONS

The concordant patterns of spatial and demographic expansions detected within phylogroups, together with the high incidence of ancestral haplotypes in southern regions of several phylogroups, suggests a pattern of contraction of populations into southern refugia during adverse climatic conditions from which subsequent northern expansions occurred. This study supports the emergent pattern of multiple refugia within Iberia but adds to it by identifying a pattern of refugia coincident with the southern distribution limits of individual evolutionary lineages. These areas are important in terms of long-term species persistence and therefore important areas for conservation.

Multiple approaches to detect outliers in a genome scan for selection in ocellated lizards (Lacerta lepida) along an environmental gradient

Identification of loci with adaptive importance is a key step to understand the speciation process in natural populations, because those loci are responsible for phenotypic variation that affects fitness in different environments. We conducted an AFLP genome scan in populations of ocellated lizards (Lacerta lepida) to search for candidate loci influenced by selection along an environmental gradient in the Iberian Peninsula. This gradient is strongly influenced by climatic variables, and two subspecies can be recognized at the opposite extremes: L. lepida iberica in the northwest and L. lepida nevadensis in the southeast. Both subspecies show substantial morphological differences that may be involved in their local adaptation to the climatic extremes. To investigate how the use of a particular outlier detection method can influence the results, a frequentist method, DFDIST, and a Bayesian method, BayeScan, were used to search for outliers influenced by selection. Additionally, the spatial analysis method was used to test for associations of AFLP marker band frequencies with 54 climatic variables by logistic regression. Results obtained with each method highlight differences in their sensitivity. DFDIST and BayeScan detected a similar proportion of outliers (3-4%), but only a few loci were simultaneously detected by both methods. Several loci detected as outliers were also associated with temperature, insolation or precipitation according to spatial analysis method. These results are in accordance with reported data in the literature about morphological and life-history variation of L. lepida subspecies along the environmental gradient.

Phylogeographic patterns of genetic diversity in eastern Mediterranean water frogs have been determined by geological processes and climate change in the Late Cenozoic

AIM: Our aims were to assess the phylogeographic patterns of genetic diversity in eastern Mediterranean water frogs and to estimate divergence times using different geological scenarios. We related divergence times to past geological events and discuss the relevance of our data for the systematics of eastern Mediterranean water frogs. LOCATION: The eastern Mediterranean region. METHODS: Genetic diversity and divergence were calculated using sequences of two protein-coding mitochondrial (mt) genes: ND2 (1038 bp, 119 sequences) and ND3 (340 bp, 612 sequences). Divergence times were estimated in a Bayesian framework under four geological scenarios representing alternative possible geological histories for the eastern Mediterranean. We then compared the different scenarios using Bayes factors and additional geological data. RESULTS: Extensive genetic diversity in mtDNA divides eastern Mediterranean water frogs into six main haplogroups (MHG). Three MHGs were identified on the Anatolian mainland; the most widespread MHG with the highest diversity is distributed from western Anatolia to the northern shore of the Caspian Sea, including the type locality of Pelophylax ridibundus. The other two Anatolian MHGs are restricted to south-eastern Turkey, occupying localities west and east of the Amanos mountain range. One of the remaining three MHGs is restricted to Cyprus; a second to the Levant; the third was found in the distribution area of European lake frogs (P. ridibundus group), including the Balkans. MAIN CONCLUSIONS: Based on geological evidence and estimates of genetic divergence we hypothesize that the water frogs of Cyprus have been isolated from the Anatolian mainland populations since the end of the Messinian salinity crisis (MSC), i.e. since c. 5.5-5.3 Ma, while our divergence time estimates indicate that the isolation of Crete from the mainland populations (Peloponnese, Anatolia) most likely pre-dates the MSC. The observed rates of divergence imply a time window of c. 1.6-1.1 million years for diversification of the largest Anatolian MHG; divergence between the two other Anatolian MHGs may have begun about 3.0 Ma, apparently as a result of uplift of the Amanos Mountains. Our mtDNA data suggest that the Anatolian water frogs and frogs from Cyprus represent several undescribed species.

Isolation by environmental distance in mobile marine species: molecular ecology of franciscana dolphins at their southern range

The assessment of population structure is a valuable tool for studying the ecology of endangered species and drafting conservation strategies. As we enhance our understanding about the structuring of natural populations, it becomes important that we also understand the processes behind these patterns. However, there are few rigorous assessments of the influence of environmental factors on genetic patterns in mobile marine species. Given their dispersal capabilities and localized habitat preferences, coastal cetaceans are adequate study species for evaluating environmental effects on marine population structure. The franciscana dolphin, a rare coastal cetacean endemic to the Western South Atlantic, was studied to examine these issues. We analysed genetic data from the mitochondrial DNA and 12 microsatellite markers for 275 franciscana samples utilizing frequency-based, maximum-likelihood and Bayesian algorithms to assess population structure and migration patterns. This information was combined with 10 years of remote sensing environmental data (chlorophyll concentration, water turbidity and surface temperature). Our analyses show the occurrence of genetically isolated populations within Argentina, in areas that are environmentally distinct. Combined evidence of genetic and environmental structure suggests that isolation by distance and a process here termed isolation by environmental distance can explain the observed correlations. Our approach elucidated important ecological and conservation aspects of franciscana dolphins, and has the potential to increase our understanding of ecological processes influencing genetic patterns in other marine species.

Spatial genetic structure of Taxus baccata L. in the western Mediterranean Basin: past and present limits to gene movement over a broad geographic scale

English yew (Taxus baccata L., Taxaceae), a Tertiary relict, provides a seminal example of a widespread albeit locally endangered (often close to extinction) tree species. In order to gain detailed insights into the evolutionary dynamics of the species on a broad geographical scale, over 1000 trees from 91 populations of English yew in the western Mediterranean were analyzed using seven nuclear microsatellite markers. Our results revealed contrasting patterns of genetic structure at different spatial scales: genetic variation was highly structured at the local scale, while only a low proportion of the observed variation was attributed to regional differences. We also found a geographic gradient of decreasing diversity and increasing population divergence from northwest (central Europe and northern Iberian Peninsula) to southeast (Mediterranean Iberia and North Africa). The patterns revealed in this study probably reflect the combined effects of Quaternary climatic changes and recent impact of human activities, and potentially also more ancient events dating back to the Tertiary. Both climatic and anthropogenic factors seem to have conducted to a long history of population isolation, which may have contributed significantly to enhance population divergence through restricted gene flow and genetic drift.

Historical biogeography of the land snail Cornu aspersum: a new scenario inferred from haplotype distribution in the Western Mediterranean basin

Background

Despite its key location between the rest of the continent and Europe, research on the phylogeography of north African species remains very limited compared to European and North American taxa. The Mediterranean land mollusc Cornu aspersum (= Helix aspersa) is part of the few species widely sampled in north Africa for biogeographical analysis. It then provides an excellent biological model to understand phylogeographical patterns across the Mediterranean basin, and to evaluate hypotheses of population differentiation. We investigated here the phylogeography of this land snail to reassess the evolutionary scenario we previously considered for explaining its scattered distribution in the western Mediterranean, and to help to resolve the question of the direction of its range expansion (from north Africa to Europe or vice versa). By analysing simultaneously individuals from 73 sites sampled in its putative native range, the present work provides the first broad-scale screening of mitochondrial variation (cyt b and 16S rRNA genes) of C. aspersum.

Results

Phylogeographical structure mirrored previous patterns inferred from anatomy and nuclear data, since all haplotypes could be ascribed to a B (West) or a C (East) lineage. Alternative migration models tested confirmed that C. aspersum most likely spread from north Africa to Europe. In addition to Kabylia in Algeria, which would have been successively a centre of dispersal and a zone of secondary contacts, we identified an area in Galicia where genetically distinct west and east type populations would have regained contact.

Conclusions

Vicariant and dispersal processes are reviewed and discussed in the light of signatures left in the geographical distribution of the genetic variation. In referring to Mediterranean taxa which show similar phylogeographical patterns, we proposed a parsimonious scenario to account for the "east-west" genetic splitting and the northward expansion of the western (B) clade which roughly involves (i) the dispersal of ancestral (eastern) types through Oligocene terranes in the Western Mediterranean (ii) the Tell Atlas orogenesis as gene flow barrier between future west and east populations, (iii) the impact of recurrent climatic fluctuations from mid-Pliocene to the last ice age, (iv) the loss of the eastern lineage during Pleistocene northwards expansion phases.