Nuclear gene sequences resolve species phylogeny and mitochondrial introgression in Leptocarabus beetles showing trans-species polymorphisms

We studied the phylogenetic relationships among Japanese Leptocarabus ground beetles, which show extensive trans-species polymorphisms in mitochondrial gene genealogies. Simultaneous analysis of combined nuclear data with partial sequences from the long-wavelength rhodopsin, wingless, phosphoenolpyruvate carboxykinase, and 28S rRNA genes resolved the relationships among the five species, although separate analyses of these genes provided topologies with low resolution. For both the nuclear gene tree resulting from the combined data from four genes and a mitochondrial cytochrome oxidase subunit I (COI) gene tree, we applied a Bayesian divergence time estimation using a common calibration method to identify mitochondrial introgression events that occurred after speciation. Three mitochondrial lineages shared by two or three species were likely subject to introgression due to interspecific hybridization because the coalescent times for these lineages were much shorter than the corresponding speciation times estimated from nuclear gene sequences. We demonstrated that when species phylogeny is fully resolved with nuclear gene sequence data, comparative analysis of nuclear and mitochondrial gene trees can be used to infer introgressive hybridization events that might cause trans-species polymorphisms in mitochondrial gene trees.     

Genome‐wide single‐nucleotide polymorphism data reveal cryptic species within cryptic freshwater snail species—The case of the Ancylus fluviatilis species complex

DNA barcoding utilizes short standardized DNA sequences to identify species and is increasingly used in biodiversity assessments. The technique has unveiled an unforeseeably high number of morphologically cryptic species. However, if speciation has occurred relatively recently and rapidly, the use of single gene markers, and especially the exclusive use of mitochondrial markers, will presumably fail in delimitating species. Therefore, the true number of biological species might be even higher. One mechanism that can result in rapid speciation is hybridization of different species in combination with polyploidization, that is, allopolyploid speciation. In this study, we analyzed the population genetic structure of the polyploid freshwater snail Ancylus fluviatilis, for which allopolyploidization was postulated as a speciation mechanism. DNA barcoding has already revealed four cryptic species within A. fluviatilis (i.e., A. fluviatilis s. str., Ancylus sp. A–C), but early allozyme data even hint at the presence of additional cryptic lineages in Central Europe. We combined COI sequencing with high‐resolution genome‐wide SNP data (ddRAD data) to analyze the genetic structure of A. fluviatilis populations in a Central German low mountain range (Sauerland). The ddRAD data results indicate the presence of three cryptic species within A. fluviatilis s. str. occurring in sympatry and even syntopy, whereas mitochondrial sequence data only support the existence of one species, with shared haplotypes between species. Our study hence points to the limitations of DNA barcoding when dealing with organismal groups where speciation is assumed to have occurred rapidly, for example, through the process of allopolyploidization. We therefore emphasize that single marker DNA barcoding can underestimate the true species diversity and argue in strong favor of using genome‐wide data for species delimitation in such groups.     

Introgression and selection shaped the evolutionary history of sympatric sister‐species of coral reef fishes (genus: Haemulon)

Closely related marine species with large overlapping ranges provide opportunities to study mechanisms of speciation, particularly when there is evidence of gene flow between such lineages. Here, we focus on a case of hybridization between the sympatric sister‐species Haemulon maculicauda and H. flaviguttatum, using Sanger sequencing of mitochondrial and nuclear loci, as well as 2422 single nucleotide polymorphisms (SNPs) obtained via restriction site‐associated DNA sequencing (RADSeq). Mitochondrial markers revealed a shared haplotype for COI and low divergence for CytB and CR between the sister‐species. On the other hand, complete lineage sorting was observed at the nuclear loci and most of the SNPs. Under neutral expectations, the smaller effective population size of mtDNA should lead to fixation of mutations faster than nDNA. Thus, these results suggest that hybridization in the recent past (0.174–0.263 Ma) led to introgression of the mtDNA, with little effect on the nuclear genome. Analyses of the SNP data revealed 28 loci potentially under divergent selection between the two species. The combination of mtDNA introgression and limited nuclear DNA introgression provides a mechanism for the evolution of independent lineages despite recurrent hybridization events. This study adds to the growing body of research that exemplifies how genetic divergence can be maintained in the presence of gene flow between closely related species.     

Exploring species boundaries with multiple genetic loci using empirical data from non‐biting midges

Over the past decade, molecular approaches to species delimitation have seen rapid development. However, species delimitation based on a single locus, for example, DNA barcodes, can lead to inaccurate results in cases of recent speciation and incomplete lineage sorting. Here, we compare the performance of Automatic Barcode Gap Discovery (ABGD), Bayesian Poisson tree processes (PTP), networks, generalized mixed Yule coalescent (GMYC) and Bayesian phylogenetics and phylogeography (BPP) models to delineate cryptic species previously detected by DNA barcodes within Tanytarsus (Diptera: Chironomidae) non‐biting midges. We compare the results from analyses of one mitochondrial (cytochrome c oxidase subunit I [COI]) and three nuclear (alanyl‐tRNA synthetase 1 [AATS1], carbamoyl phosphate synthetase 1 [CAD1] and 6‐phosphogluconate dehydrogenase [PGD]) protein‐coding genes. Our results show that species delimitation based on multiple nuclear DNA markers is largely concordant with morphological variation and delimitations using a single locus, for example, the COI barcode. However, ABGD, GMYC, PTP and network models led to conflicting results based on a single locus and delineate species differently than morphology. Results from BPP analyses on multiple loci correspond best with current morphological species concept. In total, 10 lineages of the Tanytarsus curticornis species complex were uncovered. Excluding a Norwegian population of Tanytarsus brundini which might have undergone recent hybridization, this suggests six hitherto unrecognized species new to science. Five distinct species are well supported in the Tanytarsus heusdensis species complex, including two species new to science.     

Mito‐nuclear discordance across a recent contact zone for California voles

To examine the processes that maintain genetic diversity among closely related taxa, we investigated the dynamics of introgression across a contact zone between two lineages of California voles (Microtus californicus). We tested the prediction that introgression of nuclear loci would be greater than that for mitochondrial loci, assuming ongoing gene flow across the contact zone. We also predicted that genomic markers would show a mosaic pattern of differentiation across this zone, consistent with genomes that are semi‐permeable. Using mitochondrial cytochrome b sequences and genome‐wide loci developed via ddRAD‐seq, we analyzed genetic variation for 10 vole populations distributed along the central California coast; this transect included populations from within the distributions of both parental lineages as well as the putative contact zone. Our analyses revealed that (1) the two lineages examined are relatively young, having diverged ca. 8.5–54 kya, (2) voles from the contact zone in Santa Barbara County did not include F1 or early generation backcrossed individuals, and (3) there appeared to be little to no recurrent gene flow across the contact zone. Introgression patterns for mitochondrial and nuclear markers were not concordant; only mitochondrial markers revealed evidence of introgression, putatively due to historical hybridization. These differences in genetic signatures are intriguing given that the contact zone occurs in a region of continuous vole habitat, with no evidence of past or present physical barriers. Future studies that examine specific isolating mechanisms, such as microhabitat use and mate choice, will facilitate our understanding of how genetic boundaries are maintained in this system.     

Cryptic species and morphological plasticity in long-lived bivalves (Unionoida: Hyriidae) from inland Australia

Molecular (mitochondrial DNA, isozyme) and morphological diversity of freshwater mussels (Family Hyriidae) was examined at 21 sites encompassing four large river systems, across southwest Queensland, Australia. Evidence was found for two major morphological groups. One group, which occurred in every river system, closely matched a recognized species (Velesunio ambiguus) both morphologically and in a well-supported lineage within a mitochondrial phylogeny generated from partial cytochrome c oxidase subunit I (COI) sequences. The second group most closely matched Velesunio wilsonii in shell morphology but formed three deeply divergent mitochondrial DNA lineages. All four lineages occurred sympatrically in some areas and displayed corresponding fixed differences at nuclear allozyme loci, which suggests an absence of recent hybridization and the presence of separate species.     

Two Mitochondrial Barcodes for one Biological Species: The Case of European Kuhl's Pipistrelles (Chiroptera)

The Kuhl’s pipistrelle (Pipistrellus kuhlii) is a Western Palaearctic species of bat that exhibits several deeply divergent mitochondrial lineages across its range. These lineages could represent cryptic species or merely ancient polymorphism, but no nuclear markers have been studied so far to properly assess the taxonomic status of these lineages. We examined here two lineages occurring in Western Europe, and used both mitochondrial and nuclear markers to measure degrees of genetic isolation between bats carrying them. The sampling focused on an area of strict lineage sympatry in Switzerland but also included bats from further south, in North Africa. All individuals were barcoded for the COI gene to identify their mitochondrial lineages and five highly polymorphic microsatellite loci were used to cluster them according to their nuclear genotypes. Despite this low number of nuclear markers, all North African nuclear genotypes were grouped in a highly distinct subpopulation when compared with European samples sharing the same mitochondrial barcodes. The reverse situation prevailed in Switzerland where bats carrying distinct barcodes had similar nuclear genotypes. There was a weak east/west nuclear structure of populations, but this was independent of mitochondrial lineages as bats carrying either variant were completely admixed. Thus, the divergent mitochondrial barcodes present in Western Europe do not represent cryptic species, but are part of a single biological species. We argue that these distinct barcodes evolved in allopatry and came recently into secondary contact in an area of admixture north of the Alps. Historical records from this area and molecular dating support such a recent bipolar spatial expansion. These results also highlight the need for using appropriate markers before claiming the existence of cryptic species based on highly divergent barcodes.     

Introgressive hybridization in fishes: the biochemical evidence

Biochemical methods can detect variation at individual genetic loci, making possible the direct assessment of natural hybridization and introgression between fish populations. Protein electro-phoresis has been used to confirm and extend knowledge of many situations where species hybrids have been detected by morphological analyses. New cases of natural hybridization, including some at the subspecies level, have also been identified. Biochemical studies have provided the first conclusive evidence of natural post F₁ hybrids and of introgression between fish taxa. The strongest cases for introgression have used a combined analysis of nuclear protein genes and taxaspecific maternally inherited mitochondrial DNA variation. Information on the significance of introgression as a source of gene flow between taxa, particularly below the species level where sympatric subspecies and sibling species are involved, should expand in the future as the numbers and types of nuclear and mitochondrial DNA loci which can be assayed for variation increase. The full importance of introgressive hybridization in speciation may then be understood.     

Nuclear locus divergence at the early stages of speciation in the Orchard Oriole complex

As two lineages diverge from one another, mitochondrial DNA should evolve fixed differences more rapidly than nuclear DNA due to its smaller effective population size and faster mutation rate. As a consequence, molecular systematists have focused on the criteria of reciprocal monophyly in mitochondrial DNA for delimiting species boundaries. However, mitochondrial gene trees do not necessarily reflect the evolutionary history of the taxa in question, and even mitochondrial loci are not expected to be reciprocally monophyletic when the speciation event happened very recently. The goal of this study was to examine mitochondrial paraphyly within the Orchard Oriole complex, which is composed of Icterus spurius (Orchard Oriole) and Icterus fuertesi (Fuertes' Oriole). We increased the geographic sampling, added four nuclear loci, and used a range of population genetic and coalescent methods to examine the divergence between the taxa. With increased taxon sampling, we found evidence of clear structure between the taxa for mitochondrial DNA. However, nuclear loci showed little evidence of population structure, indicating a very recent divergence between I. spurius and I. fuertesi. Another goal was to examine the genetic variation within each taxon to look for evidence of a past founder event within the I. fuertesi lineage. Based on the high amounts of genetic variation for all nuclear loci, we found no evidence of such an event – thus, we found no support for the possible founding of I. fuertesi through a change in migratory behavior, followed by peripheral isolates speciation. Our results demonstrate that these two taxa are in the earliest stages of speciation, at a point when they have fixed differences in plumage color that are not reflected in monophyly of the mitochondrial or nuclear DNA markers in this study. This very recent divergence makes them ideal for continued studies of species boundaries and the earliest stages of speciation.     

ddRAD‐seq data reveal significant genome‐wide population structure and divergent genomic regions that distinguish the mallard and close relatives in North America

Recently evolved species typically share genetic variation across their genomes due to incomplete lineage sorting and/or ongoing gene flow. Given only subtle allele frequency differences at most loci and the expectation that divergent selection may affect only a tiny fraction of the genome, distinguishing closely related species based on multi‐locus data requires substantial genomic coverage. In this study, we used ddRAD‐seq to sample the genomes of five recently diverged, New World “mallards” (Anas spp.), a group of dabbling duck species characterized by diagnosable phenotypic differences but minimal genetic differentiation. With increased genomic sampling, we aimed to characterize population structure within this group and identify genomic regions that may have experienced divergent selection during speciation. We analyzed 3,017 autosomal ddRAD‐seq loci and 177 loci from the Z‐chromosome. In contrast to previous studies, the ddRAD‐seq data were sufficient to assign individuals to their respective species or subspecies and to generate estimates of gene flow in a phylogenetic framework. We find limited evidence of contemporary gene flow between the dichromatic mallard and several monochromatic taxa, but find evidence for historical gene flow between some monochromatic species pairs. We conclude that the overall genetic similarity of these taxa likely reflects retained ancestral polymorphism rather than recent and extensive gene flow. Thus, despite recurring cases of hybridization in this group, our results challenge the current dogma predicting the genetic extinction of the New World monochromatic dabbling ducks via introgressive hybridization with mallards. Moreover, ddRAD‐seq data were sufficient to identify previously unknown outlier regions across the Z‐chromosome and several autosomal chromosomes that may have been involved in the diversification of species in this recent radiation.     

The phylogeographic importance of the Strait of Gibraltar as a gene flow barrier in terrestrial arthropods: a case study with the scorpion Buthus occitanus as model organism

The phylogenetic relationship between Buthus occitanus populations across the Strait of Gibraltar was investigated using nuclear 18S/ITS-1 DNA sequences and mitochondrial 16S and COI DNA sequences. All analyses showed that the European samples are highly separated from North African samples, and also suggest the existence of three main groups within this species complex, i.e., an European, an Atlas (=Moroccan samples) and a Tell-Atlas group (=Tunisian samples). The European clade was subdivided into three distinct subclades. The application of a previous calibration of the molecular clock of another buthid species suggested that most of the detected mitochondrial DNA lineages including the European lineages are about three times older than the re-opening of the Gibraltar Strait, and consequently, that other and older vicariant events are responsible for the observed phylogeographic structure of this species complex. Concerning the Moroccan samples, a discordance between nuclear and mitochondrial gene markers was observed. The 18S/ITS-1 gene tree could not resolve the phylogenetic relationships among the Moroccan B. occitanus subspecies and the closely related species B. atlantis, whereas mitochondrial genes suggested the co-existence of several old phylogenetic lineages in Morocco. We hypothesized that this difference may be explained by male-biased gene flow and gene conversion at the tandemly repeated 18S/ITS-1 gene regions.     

Cryptic divergence and strong population structure in the colonial invertebrate Pycnoclavella communis (Ascidiacea) inferred from molecular data

We studied sequence variation in the mitochondrial gene cytochrome c oxidase subunit I (COI) for 135 individuals from eight Mediterranean populations of the colonial ascidian Pycnoclavella communis across most of its presently known range of distribution in the Mediterranean. Three haplotypes from Atlantic locations were also included in the study. Phylogenetic, phylogeographic and population genetic analyses were used to unravel the genetic variability within and between populations. The study revealed 32 haplotypes for COI, 29 of them grouped within two Mediterranean lineages of P. communis (mean nucleotide divergence between lineages was 8.55%). Phylogenetic and network analyses suggest the possible existence of cryptic species corresponding to these two lineages. Population genetic analyses were restricted to the five populations belonging to the main genetic lineage, and for these localities we compared the information gleaned from COI sequence data and from eight microsatellite loci. A high genetic divergence between populations was substantiated using both kinds of markers (COI, global Fst=0.343; microsatellite loci, global Fst=0.362). There were high numbers of private haplotypes (COI) and alleles (microsatellites) in the populations studied. Restricted gene flow and inbreeding occur in the present range of distribution of the species. Microsatellite loci showed a strong incidence of failed amplifications, which we attribute to the marked intraspecies variability that hampered the application of these highly specific markers. Our results show important genetic variability at all levels studied, from within populations to between basins, possibly coupled to speciation processes. This variability is attributable to restricted gene flow among populations due to short-distance dispersal of the larvae.     

Reduced introgression of the Y chromosome between subspecies of the European rabbit (Oryctolagus cuniculus) in the Iberian Peninsula

The role of the Y chromosome in speciation is unclear. Hybrid zones provide natural arenas for studying speciation, as differential introgression of markers may reveal selection acting against incompatibilities. Two subspecies of the European rabbit (Oryctolagus cuniculus) form a hybrid zone in the Iberian Peninsula. Previous work on mitochondrial DNA (mtDNA), Y- and X-linked loci revealed the existence of two divergent lineages in the rabbit genome and that these lineages are largely subspecies-specific for mtDNA and two X-linked loci. Here we investigated the geographic distribution of the two Y chromosome lineages by genotyping two diagnostic single nucleotide polymorphisms in a sample of 353 male rabbits representing both subspecies, and found that Y chromosome lineages are also largely subspecies-specific. We then sequenced three autosomal loci and discovered considerable variation in levels of differentiation at these loci. Finally, we compared estimates of population differentiation between rabbit subspecies at 26 markers and found a surprising bimodal distribution of F(ST)values. The vast majority of loci showed little or no differentiation between rabbit subspecies while a few loci, including the SRY gene, showed little or no introgression across the hybrid zone. Estimates of population differentiation for the Y chromosome were surprisingly high given that there is male-biased dispersal in rabbits. Taken together, these data indicate that there is a clear dichotomy in the rabbit genome and that some loci remain highly differentiated despite extensive gene flow following secondary contact.     

Hybrid origins of honeybees from italy (Apis mellifera ligustica) and sicily (A. m. sicula)

The genetic variability of honeybee populations Apis mellifera ligustica, in continental Italy, and of A. m. sicula, in Sicily, was investigated using nuclear (microsatellite) and mitochondrial markers. Six populations (236 individual bees) and 17 populations (664 colonies) were, respectively, analysed using eight microsatellite loci and DraI restriction fragment length polymorphism (RFLP) of the cytochrome oxidase I (COI)-cytochrome oxidase II (COII) region. Microsatellite loci globally confirmed the southeastern European heritage of both subspecies (evolutionary branch C). However, A. m. ligustica mitochondrial DNA (mtDNA) appeared to be a composite of the two European (M and C) lineages over most of the Italian peninsula, and only mitotypes from the African (A) lineage were found in A. m. sicula samples. This demonstrates a hybrid origin for both subspecies. For A. m. ligustica, the most widely exported subspecies, this hybrid origin has long been obscured by the fact that in the main area of queen production (from which most of the previous ligustica bee samples originated) the M mitochondrial lineage is absent, whereas it is present almost everywhere else in Italy. This presents a new view of the evolutionary history of European honeybees. For instance, the Iberian peninsula was considered as the unique refuge for the M branch during the quaternary ice periods. Our results show that the Apennine peninsula played a similar role. The differential distribution of nuclear and mitochondrial markers observed in Italy seems to be a general feature of introgressed honeybee populations. Presumably, it stems from the social nature of the species in which both genome compartments are differentially affected by the two (individual and colonial) reproduction levels.     

Molecular phylogeography of Apalone spinifera (Reptilia, Trionychidae)

The molecular evaluation of morphologically based species delimitations of many herpetofauna has improved the understanding of evolutionary processes and the rigor of conservation efforts. Previous evidence for a deep lineage divide between south-eastern + northern subspecies of the softshell turtle Apalone spinifera ( A. s. aspera  +  A. s. hartwegi  +  A. s. spinifera ) and western subspecies ( A. s. pallida  +  A. s. emoryi  +  A. s. guadalupensis ) was re-examined with a broader sampling using mitochondrial and nuclear markers. The south-eastern + northern clade and the western clade maintained mitochondrial reciprocal monophyly. We molecularly confirmed a geographical boundary between these two clades in Oklahoma, and developed a phylogeographical hypothesis that invokes stream capture. We evaluated whether these mitochondrial lineages represent distinct species by surveying these clades for divergence at the nuclear intron R35 and two nuclear genes, Cmos and recombination activating gene 1 (RAG-1). The nuclear loci showed no phylogenetic resolution and only the nuclear intron exhibited significant nearest neighbour statistics. Taken together, the nuclear data suggest that taxonomic elevation of the two mitochondrial clades would be currently unjustified.     

Molecular data and ecological niche modelling reveal the evolutionary history of the common and Iberian moles (Talpidae) in Europe

According to mitochondrial data, the common mole, Talpa europaea, is paraphyletic. This could be explained by either an ancient introgression of mtDNA from the Iberian blind mole T. occidentalis to T. europaea, or the existence of a differentiated taxonomic entity in northern Spain that needs to be described. In this study, we combined mitochondrial (Cytb) and nuclear (HDAC2) data to investigate these two alternative hypotheses. Based on both mitochondrial and nuclear data and an extensive geographical sampling (399 sequenced individuals), we show that the populations of T. europaea from Spain and south‐western France (south of the Loire River) are phylogenetically closer to T. occidentalis than to T. europaea. The Spanish–French lineage has some morphological characters resembling more to T. occidentalis (e.g. eyes) and others resembling more to T. europaea (external measurements, mesostyle of the first upper molar). It also seems to have several distinctive dental characters, suggesting that it should be recognized as a new species. Within the three lineages, we found a marked phylogeographical pattern, with several allopatric or parapatric lineages, dating from the Pleistocene. Our genetic data combined with species distribution models support the presence of several putative glacial refugia during glacial maxima for each species.     

DNA barcoding of bark and ambrosia beetles reveals excessive NUMTs and consistent east‐west divergence across Palearctic forests

A comprehensive DNA barcoding library is very useful for rapid identification and detection of invasive pest species. We tested the performance of species identification in the economically most damaging group of wood‐boring insects – the bark and ambrosia beetles – with particular focus on broad geographical sampling across the boreal Palearctic forests. Neighbour‐joining and Bayesian analyses of cytochrome oxidase I (COI) sequences from 151 species in 40 genera revealed high congruence between morphology‐based identification and sequence clusters. Inconsistencies with morphological identifications included the discovery of a likely cryptic Nearctic species of Dryocoetes autographus, the possible hybrid origin of shared mitochondrial haplotypes in Pityophthorus micrographus and P. pityographus, and a possible paraphyletic Xyleborinus saxeseni. The first record of Orthotomicus suturalis in North America was confirmed by DNA barcoding. The mitochondrial data also revealed consistent divergence across the Palearctic or Holarctic, confirmed in part by data from the large ribosomal subunit (28S). Some populations had considerable variation in the mitochondrial barcoding marker, but were invariant in the nuclear ribosomal marker. These findings must be viewed in light of the high number of nuclear insertions of mitochondrial DNA (NUMTs) detected in eight bark beetle species, suggesting the possible presence of additional cryptic NUMTs. The occurrence of paralogous COI copies, hybridization or cryptic speciation demands a stronger focus on data quality assessment in the construction of DNA barcoding databases.     

Genetic connectivity between north and south Mid-Atlantic Ridge chemosynthetic bivalves and their symbionts

Transform faults are geological structures that interrupt the continuity of mid-ocean ridges and can act as dispersal barriers for hydrothermal vent organisms. In the equatorial Atlantic Ocean, it has been hypothesized that long transform faults impede gene flow between the northern and the southern Mid-Atlantic Ridge (MAR) and disconnect a northern from a southern biogeographic province. To test if there is a barrier effect in the equatorial Atlantic, we examined phylogenetic relationships of chemosynthetic bivalves and their bacterial symbionts from the recently discovered southern MAR hydrothermal vents at 5°S and 9°S. We examined Bathymodiolus spp. mussels and Abyssogena southwardae clams using the mitochondrial cytochrome c oxidase subunit I (COI) gene as a phylogenetic marker for the hosts and the bacterial 16S rRNA gene as a marker for the symbionts. Bathymodiolus spp. from the two southern sites were genetically divergent from the northern MAR species B. azoricus and B. puteoserpentis but all four host lineages form a monophyletic group indicating that they radiated after divergence from their northern Atlantic sister group, the B. boomerang species complex. This suggests dispersal of Bathymodiolus species from north to south across the equatorial belt. 16S rRNA genealogies of chemoautotrophic and methanotrophic symbionts of Bathymodiolus spp. were inconsistent and did not match the host COI genealogy indicating disconnected biogeography patterns. The vesicomyid clam Abyssogena southwardae from 5°S shared an identical COI haplotype with A. southwardae from the Logatchev vent field on the northern MAR and their symbionts shared identical 16S phylotypes, suggesting gene flow across the Equator. Our results indicate genetic connectivity between the northern and southern MAR and suggest that a strict dispersal barrier does not exist.     

Humid forest refugia, speciation and secondary introgression between evolutionary lineages: differentiation in a Near Eastern brown frog, Rana macrocnemis

Brown frogs of the complex Rana macrocnemis demonstrate various degrees of differentiation between the two widespread forms, macrocnemis and camerani , in different parts of Anatolia, the Caucasus Isthmus and northern Iran. In order to reveal whether or not the forms represent monophyletic evolutionary lineages, we analysed the geographic distribution of characters of external morphology, alleles at three polymorphic allozyme loci and mitochondrial DNA haplotypes. All three character sets showed a highly congruent pattern in a limited area of humid forests in the south-west Caucasus. Frogs from this area (form macrocnemis) represent a monophyletic lineage which forms a narrow hybrid zone with populations of camerani inhabiting the southern Caucasus and north-east Turkey. However, in other parts of the group's range, hybridization between the two lineages resulted in the merger of forms and the formation of a clinal pattern of variations. The role of landscape-dependent diversifying selection in the evolution of the group is discussed. Conceptual difficulties in delimiting the borders of subspecies are pointed out.