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

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

Complex patterns of morphological and mtDNA variation in Lacerta perspicillata (Reptilia; Lacertidae)

Lacerta perspicillata is a north-west African lacertid lizard that shows considerable intraspecific variation, with three subspecies described on the basis of colour pattern and body size. Recent observations of a population containing two morphological forms and more than one deep genetic lineage, as well as an apparent lack of concordance between forms and genetic lineages, suggest that the complexity is greater than previously thought. To analyse and quantify this variation, we studied the variability within this species at two levels: (1) external morphology (multivariate analysis of scalation, body dimensions, and colour pattern) and (2) mtDNA (sequencing and single-strand confirmation polymorphism analysis). Fifty-two individuals were studied at Taza, northern Morocco. Two morphological groupings (ostensibly representing two previously described subspecies) and two deep mtDNA lineages were detected at this site, with complete correspondence between the two. This, together with an apparent lack of hybrids, would normally support respective full species recognition. However, analysis of 98 individuals from other populations demonstrated that the situation is highly complex with the same genetic lineages having reversed morphotypes in other areas, making such a designation difficult. Across the other studied populations, we found no support for any of the currently recognized subspecies. The lack of congruence between mtDNA lineages and morphometric patterns (in some cases) and the morphological similarity among lizards from different lineages suggest ecophenotypic convergence or multiple introgressive hybridization. The study highlights the tremendous complexity that may exist within a taxon and the inadequacy of older alpha-taxonomy based designations in describing it.  © 2007 The Linnean Society of London, Biological Journal of the Linnean Society , 2007, 90 , 479–490.     

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.     

Genetic structure and differentiation of Gentiana atuntsiensis W. W. Smith and G. striolata T. N. Ho (Gentianaceae) as revealed by ISSR markers

The Hengduan Mountains are the core region of the Himalaya hotspot, and are renowned for their high levels of endemism. Gentiana atuntsiensis and G. striolata are two closely related and morphologically similar species endemic to this region. In this study, inter-simple sequence repeat (ISSR) markers were used to investigate the genetic structure within these two species, as well as the differentiation between them. An analysis of molecular variance-derived estimate demonstrated only 13.5% of genetic differentiation between the two species. Considering their adjacent distribution patterns, low genetic divergence, and clear clustering into two groups, it is hypothesized that the two species arose from a rapid and recent speciation event induced mainly by geographical isolation. A relatively high level of genetic diversity was revealed in each species (Shannon's index of diversity: H _sp = 0.324 and 0.391; H _pop = 0.225 and 0.274; for G. atuntsiensis and G. striolata , respectively). Most of the genetic variation was partitioned within populations ( Ø _ST = 0.232 and 0.226 in G. atuntsiensis and G. striolata , respectively). The large population sizes, outcrossing breeding system, and small, light seeds that disperse widely may explain the genetic structure in both species.  © 2007 The Linnean Society of London. Botanical Journal of the Linnean Society , 2007, 154 , 225–232.     

Falling apart and merging: diversification of slender salamanders (Plethodontidae: Batrachoseps) in the American West

The plethodontid genus Batrachoseps , the slender salamanders, is the most diverse clade of salamanders in western North America, but it has posed taxonomic difficulties because it contains many morphologically cryptic species. A segment of the mitochondrial DNA gene cytochrome b was studied for 278 individuals densely sampled from throughout the range of all 18 described species and several undescribed species. Phylogenetic analyses of the mtDNA data identify six major clades, one corresponding to the subgenus Plethopsis and five within a monophyletic subgenus Batrachoseps. All major clades and most species within these clades display strong phylogeographic structuring. Comparisons of mtDNA and allozyme data show that several allozymically cohesive groups are not monophyletic with respect to mtDNA. We suggest that this phenomenon results from fragmentation of populations, divergence in allopatry, and then recontact and gradual merging of units caused predominantly by male-mediated gene flow. The mtDNA offers evidence that populations were once more isolated than they are now, while the patterns of allozyme variation reflect recent and current interactions among populations. The complex patterns of morphological, allozymic and mtDNA variation associated with the constantly changing geological landscape give insight into the nature of processes responsible for species formation in Batrachoseps .  © 2002 The Linnean Society of London. Biological Journal of the Linnean Society , 2002, 76 , 361–391.     

Infraspecific variation and phylogeography of the high-polyploid Iberian endemic Anthoxanthum amarum Brot. (Poaceae; Pooideae) assessed by random amplified polymorphic DNA markers (RAPDs) and morphology

High morphological and ecological diversity has been observed in the north-west Iberian Peninsula endemic Anthoxanthum amarum Brot., for which two different morphotypes (northern and southern) have been described on the basis of qualitative traits and geographical origin. In the present study, a combined molecular and morphological method was applied to ten populations of this species with the following aims: (1) to assess whether the variation observed was taxonomically meaningful; (2) to assess the influence of the environment on the variation in the morphological characters; and (3) to track the potential phylogeographical information deduced from our random amplified polymorphic DNA marker (RAPD) data in order to draw inferences about the past history of this species in the north-west Iberian Peninsula. To achieve these aims, 26 macromorphological characters were recorded in 279 specimens, and 77 RAPD phenotypes were identified in the 79 plants studied. The association analyses performed using the morphological and molecular data showed that no clear separation existed between the morphs, and a strong correlation between qualitative characters and the environment was detected. Moreover, both the multivariate analyses and the assignment test based on RAPD data revealed that the genetic variation was hierarchically structured, and three genetically distinct groups could be identified. Two of these clusters might correspond to different expansion routes proposed in the literature for different plant species in the north-west Iberian Peninsula.  © 2007 The Linnean Society of London, Botanical Journal of the Linnean Society , 2007, 155 , 179–192.     

Biogeographical connections between the Maghreb and the Mediterranean peninsulas of southern Europe

The glacial–interglacial cycles have caused severe range modifications of species' distributions. In Europe, thermophilic species had to retreat into geographically distinct southern refugia during glaciations. This process produced strong genetic imprints, which are still detectable by the present pattern of genetic differentiation and the distribution of regional diversity. To reveal the biogeographical imprints in the western Mediterranean, we analysed 26 populations of the butterfly Maniola jurtina spread over large areas of its European and North African distribution range. The samples were analysed using allozyme electrophoresis. We detected three genetic groups, divided into Western Europe, Central/Eastern Europe, and Italy with the Maghreb. The North African samples randomly cluster within the Italian samples. Even the population sampled in Morocco is genetically closely related to these samples and not to the geographically neighbouring Iberian ones. Parameters of genetic diversity showed similar values over the whole study area. The observed genetic pattern reflects possible glacial refugia in Europe located in the Iberian Peninsula and the Balkans. For North Africa and Italy, our data reveal a colonization of Africa originating from Italy.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 98 , 693–703.     

Spatio-temporal genetic variation in sympatric and allopatric Mediterranean Cicada species (Hemiptera, Cicadidae)

Cicada barbara (Stål) and Cicada orni L. are two Mediterranean cicadas, very similar in morphology, that produce distinct acoustic mating signals and that have partially overlapping distribution ranges in the Iberian Peninsula, occurring in sympatry in several locations. In the present study, six microsatellite loci were analysed in C. barbara , four of which were also analysed in C. orni . Geographical and temporal genetic variation in these species was studied. No evidence of hybridization was found, enabling us to infer that the isolating barriers between these species are efficient. Partitioning of geographic variation in each species, revealed the following patterns: Iberian Peninsula and Northwestern African populations of C. barbara showed higher differentiation between than within each region, supporting C. barbara subspecific taxonomy ( C. barbara lusitanica in the Iberian Peninsula and C. barbara barbara in Northwestern Africa) and highlighting isolation coincident with the presence of physical barriers to gene-flow; differentiation between populations of C. orni from both sides of the Pyrenees was very low, and this mountain range may not constitute a significant barrier for the dispersal of this species; Greek populations of C. orni were found to be highly differentiated from Western European populations; and no pattern of isolation-by-distance was found in either species within the Iberian Peninsula.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 96 , 249–265.     

Phylogenetic analyses of mtDNA sequences reveal three cryptic lineages in the widespread neotropical frog Leptodactylus fuscus (Schneider, 1799) (Anura, Leptodactylidae)

Leptodactylus fuscus is a neotropical frog ranging from Panamá to Argentina, to the east of the Andes mountains, and also inhabiting Margarita, Trinidad, and the Tobago islands. We performed phylogenetic analyses of 12S rRNA, 16S rRNA, tRNA-Leu, and ND1 mitochondrial (mt) DNA sequences from specimens collected across the geographic distribution of L. fuscus to examine two alternative hypotheses: (i) L. fuscus is a single, widely distributed species, or (ii) L. fuscus is a species complex. We tested statistically for geographic association and partitioning of genetic variation among mtDNA clades. The mtDNA data supported the hypothesis of several cryptic species within L. fuscus. Unlinked mtDNA and nuclear markers supported independently the distinctness of a 'northern' phylogenetic unit. In addition, the mtDNA data divided the southern populations into two clades that showed no sister relationship to each other, consistent with high differentiation and lack of gene flow among southern populations as suggested by allozyme data. Concordance between mtDNA and allozyme patterns suggests that cryptic speciation has occurred in L. fuscus without morphological or call differentiation. This study illustrates a case in which lineage splitting during the speciation process took place without divergence in reproductive isolation mechanisms (e.g. advertisement call in frogs), contrary to expectations predicted using a biological species framework.  © 2006 The Linnean Society of London, Biological Journal of the Linnean Society, 2006, 87 , 325–341. No claim to original US government works.     

Genetic structure of the critically endangered endemic Cochlearia polonica (Brassicaceae): efficiency of the last-chance transplantation

Cochlearia polonica , a narrow endemic of southern Poland, is one of the rarest and most endangered species of the European flora. All natural populations are extinct and the species has survived in only one transplanted population derived from 14 original individuals. Using AFLPs, the genetic variation and spatial structure of this population were analysed approximately 30 years after transplantation. The incidence of polymorphic AFLP bands (30.46%) is low compared with data from a natural population of another Cochlearia species, C. tatrae . Principal co-ordinates and spatial autocorrelation analyses demonstrated the presence of significant genetic structure. It is recommended that conservation efforts on C. polonica should preserve the complete population area, because local extinctions may lead to a loss of genetic information. The presence of genetic structure should also be taken into account during the sampling of material (plants or seeds) for ex situ conservation measures.  © 2007 The Linnean Society of London, Botanical Journal of the Linnean Society , 2007, 155 , 527–532.     

Microsatellite DNA polymorphism confirms reproductive isolation and reveals differences in population genetic structure of cryptic pipistrelle bat species

Previous studies have indicated that the common European pipistrelle bat ( Pipistrellus pipistrellus ) comprises two cryptic species, P. pipistrellus and Pipistrellus pygmaeus , which differ in echolocation call frequency and mitochondrial DNA sequence. However, levels of divergence based on nuclear markers have not been examined, and hence the potential for male-mediated gene flow between the species cannot be discounted. Moreover, little is known about population structure and migration patterns in either species. Here, we describe the use of microsatellites to investigate nuclear DNA differentiation between, and the pattern of population genetic structure within, the two cryptic pipistrelle species. In total, 1300 individuals from 82 maternity colonies were sampled across the British Isles and Continental Europe. We show, using multivariate analyses, that colonies of the same species are generally genetically more similar to each other than to those from the other species regardless of geographical location. Our findings support the hypothesis that the species are reproductively isolated. Significant patterns of genetic isolation by distance were identified in both species, indicating that mating may occur before any long-distance autumnal migration. The presence of a sea channel does not confer higher levels of genetic differentiation among colonies over and above distance alone in either species. Differences in genetic population structure were identified between the species, with P. pipistrellus showing a wider range of levels of genetic differentiation among colonies and a stronger relationship between genetic and geographical distance than P. pygmaeus . Differences in dispersal, mating behaviour, colony size and/or postglacial colonization patterns could contribute to the differences observed.  © 2007 The Linnean Society of London, Biological Journal of the Linnean Society , 2007, 90 , 539–550.     

Current and historic gene flow of the sand goby Pomatoschistus minutus on the European Continental Shelf and in the Mediterranean Sea

Phylogeographical patterns of the sand goby Pomatoschistus minutus (Gobiidae, Teleostei) were studied by means of sequence and single-stranded conformational polymorphism analysis of a 283-bp fragment of the cytochrome b locus of the mtDNA. A total of 228 individuals sampled at 13 sites throughout the species's distributional range revealed a moderate level of diversity and a low but significant level of overall genetic differentiation at all but one site. The goby sample from the Adriatic Sea differed in sequence by approximately 10% from the Atlantic P. minutus and is thought to belong to a cryptic species of the genus Pomatoschistus . Limited genetic differentiation with a weak pattern of isolation-by-distance was recorded throughout the distributional range of the typical P. minutus . Phylogeographical analysis suggested a contiguous range expansion in the Atlantic and Baltic basins during the Eemian and evidence for a glacial refugium in the southern North Sea during the Weichselian. In P. minutus from the western Mediterranean Sea a high number of endemic haplotypes as well as the most common Atlantic haplotype were recorded in appreciable frequencies. This might be explained by secondary contact between different mitochondrial lineages, which evolved in allopatry.   © 2004 The Linnean Society of London, Biological Journal of the Linnean Society , 2004, 83 , 561–576.     

Species-richness of the Anopheles annulipes complex (Diptera: Culicidae) revealed by tree and model-based allozyme clustering analyses

The Australasian Anopheles annulipes complex contains at least ten sibling species, some of which are important vectors of myxomatosis in rabbits. We aimed to establish how many species occurred among specimens from 61 sites throughout Australia, scored for 32 putative allozyme loci. We compared the number of species predicted from tree-based clustering of operational taxonomic units (OTUs) with that from a novel model-based Bayesian clustering approach for individual genotypes. We rejected the hypothesis of conspecificity of OTUs if they differed by at least 20% fixed differences and 0.300 Nei's standard genetic distance D . According to these criteria, 18–25 species occur, making this the most species-rich anopheline complex known to date. A conservative estimate from the Bayesian analysis was 15–20 species. There was large overlap in the assignment of individuals to clusters inferred from the Bayesian and tree-based analyses. The genetic clustering of northern and southern distributed species and an apparent cline in alleles of the enzyme glucose phosphate isomerase suggest that a latitude-dependent factor, such as temperature, may have played a role in speciation and the subsequent distribution of species. Ecological niche modelling of clusters predicted that none occur in New Guinea, emphasizing that additional, as yet unsampled, species may occur.  © 2007 The Linnean Society of London, Biological Journal of the Linnean Society , 2007, 91 , 523–539.     

Phylogeography and subspecies taxonomy of dunlins (Calidris alpina) in western Palearctic analysed by DNA microsatellites and amplified fragment length polymorphism markers

In a circumpolar wader, the dunlin ( Calidris alpina ), there are 11 named subspecies, but only five mitochondrial DNA (mtDNA) lineages have been found. In the present study, we investigated the genetic structure of dunlins in western Palearctic (from East Greenland to Taimyr peninsula) using DNA microsatellites and amplified fragment length polymorphism (AFLP) markers that may detect more recent differentiation than mtDNA. In this region, we consider four described subspecies; alpina , schinzii , arctica and centralis , together comprising two mtDNA lineages. We analyse seven polymorphic microsatellite loci and 91 AFLP markers in 287 and 152 unrelated individuals, respectively, originating from 17 populations. Neither microsatellites nor AFLPs reveal distinct groups that correspond to currently recognized subspecies. There is a clear pattern of isolation by distance in microsatellites. Our results do not contradict the former mtDNA results that there are two phylogenetic lineages (approximately corresponding to schinzii and centralis ) that have met and formed a cline ( alpina ). We find no difference between schinzii and arctica (East Greenland). We conclude that, given the lack of distinct groups and the gradual changes in microsatellite allele frequencies, these markers provide little genetic support for the dunlin subspecies taxonomy in the western Palearctic.  © 2007 The Linnean Society of London, Biological Journal of the Linnean Society , 2007, 92 , 713–726.     

Species diversity and endemism in the Daphnia of Argentina: a genetic investigation

Although the temperate regions of South America are known to have a diverse daphniid fauna, there has been no genetic evaluation of the existing taxonomic system or of the affinities between the North and South American faunas. The present study analyses mitochondrial DNA sequences and allozyme variation to investigate species diversity in 176 Daphnia populations from Argentina. This work established the presence of at least 15 species in Argentina, six of which are either undescribed or are currently misidentified and two of which represent range extensions of North American taxa. Eleven of the Argentine species appear endemic to South America, while the remaining four also occur in North America. In the latter cases, the close genetic similarity between populations from North and South America indicates the recent exchange of propagules between the continents. While biological interactions and habitat availability have undoubtedly contributed to the observed species distributions, chance dispersal has apparently played a dominant role in structuring large-scale biogeographical patterns in this genus and probably in other passively-dispersed organisms.  © 2004 The Linnean Society of London, Zoological Journal of the Linnean Society , 2004, 140 , 171−205.     

Genetic differentiation among the Maculinea species (Lepidoptera: Lycaenidae) in eastern Central Europe

The present study aimed to analyse the level of genetic variation in the eastern Central European (Slovenia, Hungary, and Romania/Transylvania) populations of the Large Blues ( Maculinea ) to analyse the pattern of differentiation both between and within the species. One objective was to compare the level of differentiation between the two disputed species ( Maculinea alcon and Maculinea rebeli ) with that among the other species. Imagos were collected from 23 localities in eastern Central Europe in 2002. Enzyme polymorphism was analysed using polyacrylamide gel electrophoresis. Fourteen enzyme loci were studied in all samples. In the analysis of the data, F -statistics and Nei's genetic distances were calculated and a dendrogram (unweighted pair group method with arithmetic mean) was constructed on the basis of the distance matrix. A multivariate analysis of variance was performed to study the pattern of genetic differentiation among the samples. Principal component analysis analysis was also carried out using the allele frequency data of the samples. Our results indicated that the large blues are generally less polymorphic than other European lycaenid butterflies studied. At the same time, the level of genetic differentiation was high, even among local populations within the species. A low level of genetic variation within the populations coupled with strong differentiation among them implies the effect of genetic drift. Strong genetic differentiation of four Maculinea species ( M. alcon , Maculinea teleius, Maculinea nausithous , and Maculinea arion ) was confirmed. Significant differentiation was not found between M. alcon and M. rebeli .  © 2007 The Linnean Society of London, Biological Journal of the Linnean Society , 2007, 91 , 11–21.     

Phenotypic divergence and reproductive isolation between sympatric forms of Japanese threespine sticklebacks

The threespine stickleback ( Gasterosteus aculeatus ) species complex is well suited for identifying the types of phenotypic divergence and isolating barriers that contribute to reproductive isolation at early stages of speciation. In the present study, we characterize the patterns of genetic and phenotypic divergence as well as the types of isolating barriers that are present between two sympatric pairs of threespine sticklebacks in Hokkaido, Japan. One sympatric pair consists of an anadromous and a resident freshwater form and shows divergence in body size between the forms, despite the lack of genetic differentiation between them. The second sympatric pair consists of two anadromous forms, which originated before the last glacial period and are currently reproductively isolated. These two anadromous forms have diverged in many morphological traits as well as in their reproductive behaviours. Both sexual isolation and hybrid male sterility contribute to reproductive isolation between the anadromous species pair. We discuss the shared and unique aspects of phenotypic divergence and reproductive isolation in the Japanese sympatric pairs compared with postglacial stickleback species pairs. Further studies of these divergent species pairs will provide a deeper understanding of the mechanisms of speciation in sticklebacks.  © 2007 The Linnean Society of London, Biological Journal of the Linnean Society , 2007, 91 , 671–685.     

Microscale variation in alpine grasslands: AFLPs reveal a high level of genotypic diversity in Primula minima (Primulaceae)

The microscale variation and spatial genetic structure of the alpine plant species Primula minima L. was analysed using AFLPs. AFLP analysis based on three primer combinations and 123 fragments revealed no identical genotypes among the 86 studied samples from a 300 × 300 cm plot. Variation within the study plot was high: Nei's gene diversity was 0.22, Shannon's information index 0.33 and the percentage of polymorphic bands was 60.9. Cluster analysis revealed four main groups of genetically similar individuals and mapping these individuals resulted in a clear spatial pattern, with samples from the same group often located close together. The observed microscale structure was corroborated using a Mantel test, which revealed significant correlation of genetic and spatial distances, and by the results of a spatial autocorrelation analysis that indicated a high level of similarity between adjacent samples. An analysis of molecular variance revealed clear differentiation (18%) between the spatial groups. Overall gene flow within the plot was 1.11 and ranged from 0.33 between the spatially most distant groups to 2.33 between directly neighbouring groups. The extraordinary level of diversity detected in this study indicates an unexpectedly strong relevance of reproduction by seed for the species P. minima in alpine grasslands. The strong microscale variation suggests, however, that there is limited dispersal of seeds. Clonal reproduction is of subsidiary importance to sexual reproduction and seems to occur only over very small distances.  © 2007 The Linnean Society of London, Botanical Journal of the Linnean Society , 2007, 155 , 549–556.     

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.     

Complex biogeographical distribution of genetic variation within Podarcis wall lizards across the Strait of Gibraltar

Abstract Aim To examine the effect of a known geological barrier on genetic variation within a wall lizard species complex. Location The Iberian Peninsula and North Africa. Methods Sequencing of partial 12S rRNA and cytochrome b mtDNA. Results The current distribution of genetic variability is not related to the opening of the Strait of Gibraltar. Conclusions Podarcis hispanica in North Africa is probably a species complex. The Strait of Gibraltar should not be used as a known barrier to gene flow in other land based organisms without careful sampling to test for multiple crossings since its formation.