Haplotype diversity of the nematode Pristionchus pacificus on Réunion in the Indian Ocean suggests multiple independent invasions

Pristionchus pacificus has been established as a nematode model system in evolutionary developmental biology and evolutionary ecology. Field studies in North and South America, Asia, Africa and Europe indicated that nematodes of the genus Pristionchus live in association with scarab beetles. Here, we describe the first account of soil‐ and beetle‐associated nematodes on an island setting by investigating the island of Réunion in the Indian Ocean. Réunion has high numbers of endemic insects and is one among several attractive islands for biodiversity studies. Being of volcanic origin, Réunion is 2–3 million years old, making it the youngest of the Mascareigne islands. We show that beetle‐ and soil‐derived nematodes on Réunion are nearly exclusively hermaphroditic, suggesting that selfing is favoured over gonochorism (outcrossing) during island colonization. Among members of four nematode genera observed on Réunion, Pristionchus pacificus was the most prevalent species. A total of 76 isolates, in association with five different scarab beetles, has been obtained for this cosmopolitan nematode. A detailed mitochondrial haplotype analysis indicates that the Réunion isolates of P. pacificus cover all four worldwide clades of the species. This extraordinary haplotype diversity suggests multiple independent invasions, most likely in association with different scarab beetles. Together, we establish Réunion as a case study for nematode island biogeography, in which the analysis of nematode population genetics and population dynamics can provide insight into evolutionary and ecological processes. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 100 , 170–179.     

Multi locus analysis of Pristionchus pacificus on La Réunion Island reveals an evolutionary history shaped by multiple introductions, constrained dispersal events and rare out-crossing

Pristionchus pacificus, recently established as a model organism in evolutionary biology, is a cosmopolitan nematode that has a necromenic association with scarab beetles. The diverse array of host beetle species and habitat types occupied by P. pacificus make it a good model for investigating local adaptation to novel environments. Presence of P. pacificus on La Réunion Island, a young volcanic island with a dynamic geological history and a wide variety of ecozones, facilitates such investigation in an island biogeographic setting. Microsatellite data from 20 markers and 223 strains and mitochondrial sequence data from 272 strains reveal rich genetic diversity among La Réunion P. pacificus isolates, shaped by differentially timed introductions from diverse sources and in association with different beetle species. Distinctions between volcanic zones and between arid western and wet eastern climatic zones have likely limited westward dispersal of recently colonized lineages and maintained a genetic distinction between eastern and western clades. The highly selfing lifestyle of P. pacificus contributes to the strong fine‐scale population structure detected, with each beetle host harbouring strongly differentiated assemblages of strains. Periodic out‐crossing generates admixture between genetically diverse lineages, creating a diverse array of allelic combinations likely to increase the evolutionary potential of the species and facilitate adaptation to local environments and beetle hosts.     

The nematode Pristionchus pacificus as a model system for integrative studies in evolutionary biology

Comprehensive studies of evolution have historically been hampered by the division among disciplines. Now, as biology moves towards an ‘‐omics’ era, it is more important than ever to tackle the evolution of function and form by considering all those research areas involved in the regulation of phenotypes. Here, we review recent attempts to establish the nematode Pristionchus pacificus as a model organism that allows integrative studies of development and evo‐devo, with ecology and population genetics. Originally developed for comparative study with the nematode Caenorhabditis elegans, P. pacificus provided insight into developmental pathways including dauer formation, vulva and gonad development, chemosensation, innate immunity and neurobiology. Its subsequent discovery across a wide geographic distribution in association with scarab beetles enabled its evaluation in a biogeographic context. Development of an evolutionary field station on La Réunion Island, where P. pacificus is present in high abundance across a number of widespread habitat types, allows examination of the microfacets of evolution – processes of natural selection, adaptation and drift among populations can now be examined in this island setting. The combination of laboratory‐based functional studies with fieldwork in P. pacificus has the long‐term prospective to provide both proximate (mechanistic) and ultimate (evolutionary and ecological) causation and might therefore help to overcome the long‐term divide between major areas in biology.     

Natural variation in chemosensation: lessons from an island nematode

All organisms must interact with their environment, responding in behavioral, chemical, and other ways to various stimuli throughout their life cycles. Characterizing traits that directly represent an organism's ability to sense and react to their environment provides useful insight into the evolution of life‐history strategies. One such trait for the nematode Pristionchus pacificus, chemosensation, is involved in navigation to beetle hosts. Essential for the survival of the nematode, chemosensory behavior may be subject to variation as nematodes discriminate among chemical cues to complete their life cycle. We examine this hypothesis using natural isolates of P. pacificus from La Réunion Island. We select strains from a variety of La Réunion beetle hosts and geographic locations and examine their chemoattraction response toward organic compounds, beetle washes, and live beetles. We find that nematodes show significant differences in their response to various chemicals and are able to chemotax to live beetles in a novel assay. Further, strains can discriminate among different cues, showing more similar responses toward beetle washes than to organic compounds in cluster analyses. However, we find that variance in chemoattraction response is not significantly associated with temperature, location, or beetle host. Rather, strains show a more concerted response toward compounds they most likely directly encounter in the wild. We suggest that divergence in odor‐guided behavior in P. pacificus may therefore have an important ecological component.     

Zur verlorenen Taubenwelt der Insel Réunion – Neue Erkenntnisse über ausgerottete Vertreter der Ordnung Columbiformes auf der Inselgruppe der Maskarenen

Since the beginning of the influence of the fauna of the Mascarene Islands as a result of activities by Europeans, various members of Columbiformes have become extinct on Réunion Island. The erection of an obviously involved species of the former family Raphidae proved to be erroneous, indeed, but nevertheless it was assumed that there were up to four different endemic members of Columbidae that vanished from the said island. Inquiries by the help of both historical and modern literature showed that on Réunion Island overall “only” three indigenous pigeon species were wiped out: Nesoenas duboisi, Alectroenas sp., and one more form of the former genus that because of the lack of sufficient subfossil bone remains cannot be definitely identified.     

Morphological and plumage colour variation in the Réunion grey white‐eye (Aves: Zosterops borbonicus): assessing the role of selection

The Réunion grey white‐eye (Zosterops borbonicus), a small passerine endemic to the island of Réunion (Mascarene archipelago), constitutes an extraordinary case of phenotypic variation within a bird species, with conspicuous plumage colour differentiation at a microgeographical scale. To understand whether natural selection could explain such variability, we compared patterns of variation in morphological and plumage colour traits within and among populations. To quantify morphological variation, we used measurements obtained by Frank Gill in the 1960s from 239 individuals collected in 60 localities distributed over the entire island of Réunion. To quantify colour variation, we measured the reflectance spectra of plumage patches of 50 males from a subset of Gill's specimens belonging to the five recognized plumage colour variants and used a visual model to project these colours in an avian‐appropriate, tetrachromatic, colour space. We found that variants occupy different regions of the avian colour space and that between‐variant differences for most plumage patches could be discriminated by the birds. Differences in morphology were also detected, but these were, in general, smaller than colour differences. Overall, we found that variation in both plumage colour and morphology among variants is greater than would be expected if genetic drift alone was responsible for phenotypic divergence. As the plumage colour variants correspond to four geographical forms, our results suggest that phenotypic evolution in the Réunion grey white‐eye is at least partly explained by divergent selection in different habitats or regions. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2015, 114 , 459–473.     

Environmental Variables Explain Genetic Structure in a Beetle-Associated Nematode

The distribution of a species is a complex expression of its ecological and evolutionary history and integrating population genetic, environmental, and ecological data can provide new insights into the effects of the environment on the population structure of species. Previous work demonstrated strong patterns of genetic differentiation in natural populations of the hermaphroditic nematode Pristionchus pacificus in its La Réunion Island habitat, but gave no clear understanding of the role of the environment in structuring this variation. Here, we present what is to our knowledge the first study to statistically evaluate the role of the environment in shaping the structure and distribution of nematode populations. We test the hypothesis that genetic structure in P. pacificus is influenced by environmental variables, by combining population genetic analyses of microsatellite data from 18 populations and 370 strains, with multivariate statistics on environmental data, and species distribution modelling. We assess and quantify the relative importance of environmental factors (geographic distance, altitude, temperature, precipitation, and beetle host) on genetic variation among populations. Despite the fact that geographic populations of P. pacificus comprise vast genetic diversity sourced from multiple ancestral lineages, we find strong evidence for local associations between environment and genetic variation. Further, we show that significantly more genetic variation in P. pacificus populations is explained by environmental variation than by geographic distances. This supports a strong role for environmental heterogeneity vs. genetic drift in the divergence of populations, which we suggest may be influenced by adaptive forces.     

Spectrum of CFTR mutations on Réunion Island: impact on neonatal screening

The large heterogeneity in the cystic fibrosis (CF) gene is the main difficulty for genotype characterization. Numerous studies have reported considerable variations in frequencies of CF transmembrane conductance regulator (CFTR) mutations in different populations, such as African, Asian, or European populations. To completely characterize the spectrum of mutations in the CFTR gene in the Réunion Island population, we screened 228 CF chromosomes using denaturing high-pressure liquid chromatography and denaturing gradient gel electrophoresis following by direct sequencing. We identified 27 mutations, accounting for 93% of CF chromosomes. They included three novel mutations (M1T, 3121-3C-->G, and L1324P), which are described in this paper. The detection of such a high proportion of Réunion Island CFTR mutations is important for improving neonatal screening of CF on Réunion Island.     

The geographic scale of diversification on islands: genetic and morphological divergence at a very small spatial scale in the Mascarene grey white-eye (Aves: <Emphasis Type="Italic">Zosterops borbonicus</Emphasis>)

Background  

Oceanic islands provide unique scenarios for studying the roles of geography and ecology in driving population divergence and speciation. Assessing the relative importance of selective and neutral factors in driving population divergence is central to understanding how such divergence may lead to speciation in small oceanic islands, where opportunities for gene flow and population mixing are potentially high. Here we report a case of genetic and morphological structure in the Mascarene grey white-eye (Zosterops borbonicus) a species that shows a striking, geographically structured plumage polymorphism on the topographically and ecologically complex island of Réunion, yet is monotypic on the relatively uniform neighbouring island of Mauritius.     

Demographic history of a recent invasion of house mice on the isolated Island of Gough

Island populations provide natural laboratories for studying key contributors to evolutionary change, including natural selection, population size and the colonization of new environments. The demographic histories of island populations can be reconstructed from patterns of genetic diversity. House mice (Mus musculus) inhabit islands throughout the globe, making them an attractive system for studying island colonization from a genetic perspective. Gough Island, in the central South Atlantic Ocean, is one of the remotest islands in the world. House mice were introduced to Gough Island by sealers during the 19th century and display unusual phenotypes, including exceptionally large body size and carnivorous feeding behaviour. We describe genetic variation in Gough Island mice using mitochondrial sequences, nuclear sequences and microsatellites. Phylogenetic analysis of mitochondrial sequences suggested that Gough Island mice belong to Mus musculus domesticus, with the maternal lineage possibly originating in England or France. Cluster analyses of microsatellites revealed genetic membership for Gough Island mice in multiple coastal populations in Western Europe, suggesting admixed ancestry. Gough Island mice showed substantial reductions in mitochondrial and nuclear sequence variation and weak reductions in microsatellite diversity compared with Western European populations, consistent with a population bottleneck. Approximate Bayesian computation (ABC) estimated that mice recently colonized Gough Island (~100 years ago) and experienced a 98% reduction in population size followed by a rapid expansion. Our results indicate that the unusual phenotypes of Gough Island mice evolved rapidly, positioning these mice as useful models for understanding rapid phenotypic evolution.     

The role of selection and historical factors in driving population differentiation along an elevational gradient in an island bird

Adaptation to local environmental conditions and the range dynamics of populations can influence evolutionary divergence along environmental gradients. Thus, it is important to investigate patterns of both phenotypic and genetic variations among populations to reveal the respective roles of these two types of factors in driving population differentiation. Here, we test for evidence of phenotypic and genetic structure across populations of a passerine bird (Zosterops borbonicus) distributed along a steep elevational gradient on the island of Réunion. Using 11 microsatellite loci screened in 401 individuals from 18 localities distributed along the gradient, we found that genetic differentiation occurred at two spatial levels: (i) between two main population groups corresponding to highland and lowland areas, respectively, and (ii) within each of these two groups. In contrast, several morphological traits varied gradually along the gradient. Comparison of neutral genetic differentiation (F_ST) and phenotypic differentiation (P_ST) showed that P_ST largely exceeds F_ST at several morphological traits, which is consistent with a role for local adaptation in driving morphological divergence along the gradient. Overall, our results revealed an area of secondary contact midway up the gradient between two major, cryptic, population groups likely diverged in allopatry. Remarkably, local adaptation has shaped phenotypic differentiation irrespective of population history, resulting in different patterns of variation along the elevational gradient. Our findings underscore the importance of understanding both historical and selective factors when trying to explain variation along environmental gradients.     

Importance of host oviposition pattern and plant size for the selection of Trichogramma strains to control the diamondback moth

The aim of this work was to select a candidate strain of Trichogrammatidae (Hymenoptera: Chalcidoidea) to control the diamondback moth, Plutella xylostella L. (Lepidoptera: Plutellidae), an important pest of cabbage in Europe. The parasitic efficiency of Trichogramma chilonis Ishii from Japan and Réunion Island, Trichogramma evanescens Westwood from Egypt, Trichogramma ostriniae Pang & Chen from Japan, and Trichogramma semblidis (Aurivillius) from France was studied with sentinel eggs in greenhouse‐grown cauliflower. The percentage of parasitized eggs was measured 1 m from the release point and 3 days after release. Two factors were studied: host‐egg density (three eggs per stem vs. 50 eggs per stem) and plant size (50–100 cm high plants vs. 100–150 cm high plants). Trichogramma evanescens from Egypt achieved the best parasitism. Parasitism efficiency of T. ostriniae from Japan and T. chilonis from Réunion Island was not influenced by egg density. In contrast, parasitism efficiencies of the three other strains were higher at the lower egg density than at the higher egg density. Parasitism efficiency was lower when the plants were bigger, with the exception of T. chilonis from Réunion Island, for which the level of parasitism was not affected by plant size. Thus, T. chilonis from Réunion Island appeared to be the most suitable Trichogramma strain for use in controlling P. xylostella, regardless of the growth stage of the crop. However, T. ostriniae from Japan may also be more suitable for use on cauliflower crops at the early stage of crop development. If potential problems related to the introduction of exotic species are taken into consideration, it is possible to consider the native T. evanescens.     

Phylogenetic divergence of island biotas: Molecular dates,extinction, and “relict” lineages

Island formation is a key driver of biological evolution, and several studies have used geological ages of islands to calibrate rates of DNA change. However, many islands are home to “relict” lineages whose divergence apparently pre‐dates island age. The geologically dynamic New Zealand (NZ) archipelago sits upon the ancient, largely submerged continent Zealandia, and the origin and age of its distinctive biota have long been contentious. While some researchers have interpreted NZ's biota as equivalent to that of a post‐Oligocene island, a recent review of genetic studies identified a sizeable proportion of pre‐Oligocene “relict” lineages, concluding that much of the biota survived an incomplete drowning event. Here, we assemble comparable genetic divergence data sets for two recently formed South Pacific archipelagos (Lord Howe; Chatham Islands) and demonstrate similarly substantial proportions of relict lineages. Similar to the NZ biota, our island reviews provide surprisingly little evidence for major genetic divergence “pulses” associated with island emergence. The dominance of Quaternary divergence estimates in all three biotas may highlight the importance of rapid biological turnover and new arrivals in response to recent climatic and/or geological disturbance and change. We provide a schematic model to help account for discrepancies between expected versus observed divergence‐date distributions for island biotas, incorporating the effects of both molecular dating error and lineage extinction. We conclude that oceanic islands can represent both evolutionary “cradles” and “museums” and that the presence of apparently archaic island lineages does not preclude dispersal origins.     

Late Pleistocene lineage divergence among populations of Neolitsea sericea (Lauraceae) across a deep sea‐barrier in the Ryukyu Islands

Aim Our goals were: (1) to assess the population genetic structure and demographic divergence history of a bird‐dispersed tree, Neolitsea sericea, endemic to East Asian land‐bridge islands; and (2) to interpret the results in the light of controversies over the dating and configurations of land bridges through the Japanese Ryukyu Island Arc. Location Japan–Ryukyu–Taiwan Island Arc and Chinese/Korean offshore islands. Methods We applied 10 nuclear microsatellites (nSSRs) and one chloroplast (cp) DNA sequence marker (psbA–trnH intergenic spacer) to 31 populations (397 and 326 individuals, respectively) from throughout the species’ range to infer current patterns of genetic diversity and structure, and pollen‐to‐seed migration ratios (r). A coalescent‐based isolation‐with‐migration (IM) model was applied to the combined nSSR/cpDNA data set to estimate lineage divergence times and population demographic parameters. Results The geographic structure of nSSRs and the distribution of most cpDNA haplotypes revealed two distinct lineages located in areas north and south of the ‘Tokara Gap’, a narrow (c. 37 km wide) but deep (> 1000 m) sea‐strait between the northern and central Ryukyus. Based on the IM analyses, we (1) dated the divergence of these northern and southern lineages to c. 0.07 Ma (90% highest posterior density interval: 0.02–0.38 Ma); (2) estimated a slightly smaller effective population size for the northern compared to the southern lineage; and (3) recovered only trivial signals of post‐divergence gene flow between them. Main conclusions The estimated divergence time for northern and southern lineages is consistent with geological evidence for the existence of land connections in the Tokara region during cold stages of the latest Pleistocene; it is thus incompatible with an ‘ancient sea‐barrier hypothesis’ for the Ryukyu Arc, where we would have expected much older divergences related to the initial formation of the Tokara and Kerama tectonic straits during the Pliocene. Multiple factors are likely to have had a role in the divergence of N. sericea, including not only land‐bridge submergence, but also island configuration, and/or constraints on adaptation along a latitudinal temperature gradient.     

Nematode biphasic ‘boom and bust’ dynamics are dependent on host bacterial load while linking dauer and mouth-form polyphenisms

Cross-kingdom interactions involve dynamic processes that shape terrestrial ecosystems and represent striking examples of co-evolution. The multifaceted relationships of entomopathogenic nematodes with their insect hosts and symbiotic bacteria are well-studied cases of co-evolution and pathogenicity. In contrast, microbial interactions in soil after the natural death of insects and other invertebrates are minimally understood. In particular, the turnover and succession of nematodes and bacteria during insect decay have not been well documented - although it represents a rich ecological niche with multiple species interactions. Here, we utilize developmentally plastic nematode Pristionchus pacificus and its associated scarab beetles as models. On La Réunion Island, we collected rhinoceros beetle Oryctes borbonicus, induced death, and placed carcasses in cages both on the island and in a mock-natural environment in the laboratory controlling for high spatial and temporal resolution. Investigating nematode population density and dispersal dynamics, we were able to connect two imperative plasticities, dauer and mouth form. We observed a biphasic ‘boom and bust’ dispersal dynamic of dauer larvae that corresponds to bacterial load on carcasses but not bacterial type. Strikingly, all post-dauer adults have the predatory mouth form, demonstrating novel intricate interactions on decaying insect hosts. Thus, ecologically relevant survival strategies incorporate critical plastic traits.     

Reproductive compatibility and genetic and morphometric variability among populations of the predatory mite,Amblyseius largoensis (Acari: Phytoseiidae), from Indian Ocean Islands and the Americas

The red palm mite (RPM), Raoiella indica Hirst (Acari: Tenuipalpidae), is an invasive phytophagous mite that was recently introduced into The Americas. The predatory mite Amblyseius largoensis Muma (Acari: Phytoseiidae) has been the only natural enemy consistently found in association with RPM. This study aimed to determine if A. largoensis populations from the Indian Ocean Islands (La Réunion and Mauritius) and the Americas (Brazil, Trinidad and Tobago and the USA) consist a taxonomic unit or a group of cryptic species. First, the morphological variability among the A. largoensis populations from these areas was evaluated through morphometric analyses of 36 morphological traits. Then, their genetic variability and phylogenetic relationships were assessed based on two target DNA fragments: the nuclear Internal Transcribed Spacer and the mithochondrial 12S rRNA. Finally, reproductive compatibility of the populations from La Réunion and Roraima, Brazil was evaluated. Morphometric differences between the A. largoensis specimens from La Réunion Island and the Americas were observed, most of them on the length of the setae. Molecular analysis indicated that the A. largoensis populations from the Indian Ocean Islands and the Americas belong to the same taxonomic entity, although to two well defined genetic groups. Crossings involving the A. largoensis populations from La Réunion Island and Roraima, Brazil revealed complete reproductive compatibility between these populations. Information on the morphometric and genetic variability among studied A. largoensis populations can be further exploited in future studies to follow colonization of Indian Ocean Islands populations in the Americas, in the case of field releases.     

Microsatellite DNA markers for the Stable Fly,Stomoxys calcitrans (Diptera: Muscidae)

Eight polymorphic microsatellite markers have been isolated from a microsatellite‐enriched DNA library from the Stable Fly, Stomoxys calcitrans. These loci exhibited four to 15 alleles per locus and an expected heterozygosity ranging from 0.19 to 0.84 in the three populations studied from La Réunion island (Indian Ocean). They should therefore be valuable for studying genetic diversity and population structure. Cross‐species amplification of these eight loci in the closely allied species Stomoxys niger niger was successful for four of these loci, one of which was monomorphic and one of which strongly departed from Hardy–Weinberg expectations.     

Do island plant populations really have lower genetic variation than mainland populations? Effects of selection and distribution range on genetic diversity estimates

Ecological and evolutionary studies largely assume that island populations display low levels of neutral genetic variation. However, this notion has only been formally tested in a few cases involving plant taxa, and the confounding effect of selection on genetic diversity (GD) estimates based on putatively neutral markers has typically been overlooked. Here, we generated nuclear microsatellite and plastid DNA sequence data in Periploca laevigata, a plant taxon with an island–mainland distribution area, to (i) investigate whether selection affects GD estimates of populations across contrasting habitats; and (ii) test the long‐standing idea that island populations have lower GD than their mainland counterparts. Plastid data showed that colonization of the Canary Islands promoted strong lineage divergence within P. laevigata, which was accompanied by selective sweeps at several nuclear microsatellite loci. Inclusion of loci affected by strong divergent selection produced a significant downward bias in the GD estimates of the mainland lineage, but such underestimates were substantial (>14%) only when more than one loci under selection were included in the computations. When loci affected by selection were removed, we did not find evidence that insular Periploca populations have less GD than their mainland counterparts. The analysis of data obtained from a comprehensive literature survey reinforced this result, as overall comparisons of GD estimates between island and mainland populations were not significant across plant taxa (N = 66), with the only exception of island endemics with narrow distributions. This study suggests that identification and removal of markers potentially affected by selection should be routinely implemented in estimates of GD, particularly if different lineages are compared. Furthermore, it provides compelling evidence that the expectation of low GD cannot be generalized to island plant populations.     

Pleistocene glacial cycles drove lineage diversification and fusion in the Yosemite toad (Anaxyrus canorus)

Species endemic to alpine environments can evolve via steep ecological selection gradients between lowland and upland environments. Additionally, many alpine environments have faced repeated glacial episodes over the past two million years, fracturing these endemics into isolated populations. In this “glacial pulse” model of alpine diversification, cycles of allopatry and ecologically divergent glacial refugia play a role in generating biodiversity, including novel admixed (“fused”) lineages. We tested for patterns of glacial pulse lineage diversification in the Yosemite toad (Anaxyrus [Bufo] canorus), an alpine endemic tied to glacially influenced meadow environments. Using double‐digest RADseq on populations densely sampled from a portion of the species range, we identified nine distinct lineages with divergence times ranging from 18 to 724 thousand years ago (ka), coinciding with multiple Sierra Nevada glacial events. Three lineages have admixed origins, and demographic models suggest these fused lineages have persisted throughout past glacial cycles. Directionality indices supported the hypothesis that some lineages recolonized Yosemite from east of the ice sheet, whereas other lineages remained in western refugia. Finally, refugial niche reconstructions suggest that low‐ and high‐elevation lineages have convergently adapted to similar climatic niches. Our results suggest glacial cycles and refugia may be important crucibles of adaptive diversity across deep evolutionary time.     

Multilocus genetic diversity and historical biogeography of the endemic wall lizard from Ibiza and Formentera,Podarcis pityusensis (Squamata: Lacertidae)

Two monophyletic sister species of wall lizards inhabit the two main groups of Balearic Islands: Podarcis lilfordi from islets and small islands around Mallorca and Menorca and Podarcis pityusensis from Ibiza, Formentera and associated islets. Genetic diversity within the endangered P. lilfordi has been well characterized, but P. pityusensis has not been studied in depth. Here, 2430 bp of mtDNA and 15 microsatellite loci were analysed from P. pityusensis populations from across its natural range. Two main genetic groupings were identified, although geographical structuring differed slightly between the mtDNA and the nuclear loci. In general, individuals from islets/islands adjacent to the main island of Ibiza were genetically distinct from those from Formentera and the associated Freus islands for both mtDNA and the nuclear loci. However, most individuals from the island of Ibiza were grouped with neighbouring islets/islands for nuclear loci, but with Formentera and Freus islands for the mitochondrial locus. A time‐calibrated Bayesian tree was constructed for the principal mitochondrial lineages within the Balearics, using the multispecies coalescent model, and provided statistical support for divergence of the two main P. pityusensis lineages 0.111–0.295 Ma. This suggests a mid‐late Pleistocene intraspecific divergence, compared with an early Pleistocene divergence in P. lilfordi, and postdates some major increases in sea level between 0.4 and 0.6 Ma, which may have flooded Formentera. The program IMa2 provided a posterior divergence time of 0.089–0.221 Ma, which was similar to the multispecies coalescent tree estimate. More significantly, it indicated low but asymmetric effective gene copy migration rates, with higher migration from Formentera to Ibiza populations. Our findings suggest that much of the present‐day diversity may have originated from a late Pleistocene colonization of one island group from the other, followed by allopatric divergence of these populations. Subsequent gene flow between these insular groups seems likely to be explained by recent human introductions. Two evolutionary significant units can be defined for P. pityusensis but these units would need to exclude the populations that have been the subjects of recent admixture.