Unisexual rock lizard might be outcompeting its bisexual progenitors in the Caucasus

We compared the distributions, abundances and ecological requirements of parthenogenetic lizard Darevskia ‘dahli’ and its bisexual progenitors, D. portschinskii and D. mixta, in Georgia. We developed a regression model relating the species abundances with the distribution of climates. Darevskia portschinskii lives in warmer and drier climates than D. mixta; D. ‘dahli’ has the intermediate requirements. Temperature is more important than humidity for D. portschinskii, humidity is more important for D. mixta and both temperature and humidity are important for D. ‘dahli’. Suitable habitats of all three species overlap broadly; however, the observed ranges partly overlap only for D. ‘dahli’ and D. portschinskii. The observed abundance of each species, related to its predicted abundance, is lower at the sites with potential competitors. Darevskia ‘dahli’ occupies a higher proportion of the suitable habitats and has higher abundances than the progenitor species. Competition with D. ‘dahli’ is an important factor determining current distribution pattern of D. portschinskii and D. mixta. The parthenogen is a stronger competitor than the bisexual breeders and potential advantages of the bisexual reproduction remain unrealized in the given temporal and spatial scale. To explain domination of bisexually breeding lizards on the global scale, considering climate changes in geological timescale is necessary. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 101 , 447–460.     

Secondary contact followed by gene flow between divergent mitochondrial lineages of a widespread Neotropical songbird (Zonotrichia capensis)

Understanding how genetic and phenotypic differences that arise in geographically isolated populations influence the outcome of secondary contact advances our knowledge of speciation. In the present study, we investigate the secondary contact between divergent lineages of a widespread Neotropical songbird, the Rufous‐collared sparrow (Zonotrichia capensis). Zonotrichia capensis is morphologically and behaviourally diverse, and shows a pattern of lineage diversification produced by a Pleistocene expansion and colonization of South America from a probable Central American origin. Consistent with previous results, we find three lineages throughout the species range, showing between 1.5% and 2.5% divergence in mitochondrial control region sequences. These lineages come into secondary contact in the Dominican Republic, La Paz (Bolivia), and North‐eastern Argentina. We use DNA microsatellite data to study a broad secondary contact zone in North‐eastern Argentina, finding that Bayesian clustering analyses do not assign individuals to their respective mitochondrial lineages. Overall, we did not observe nuclear genetic discontinuities in the study area. We conclude that, if genetic, morphological, and/or cultural differences accumulated among lineages during isolation, they were insufficient to prevent gene flow after secondary contact. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 111 , 863–868.     

Nuclear DNA recapitulates the cryptic mitochondrial lineages of Lumbricus rubellus and suggests the existence of cryptic species in an ecotoxological soil sentinel

Mitochondrial DNA analysis has revealed two distinct phylogenetic lineages within the ecotoxological sentinel earthworm model Lumbricus rubellus Hoffmeister, 1843. The existence of these lineages could complicate ecotoxicological studies that use the species as a sentinel for soil contamination testing, as they may respond differently to contamination; however, as mitochondrial haplotypes are not always expected to segregate in the same way as chromosomal DNA in natural populations, we further investigated this issue by using nuclear DNA markers (microsatellites) to measure genetic diversity, differentiation, and gene flow in sympatric populations of the two L. rubellus lineages at two sites in South Wales. Our results show that sympatric populations of the two lineages are more genetically differentiated than geographically distant populations of the same lineage, and Bayesian clustering analysis revealed no evidence of gene flow between the lineages at either site. Additionally, DNA sequencing of these microsatellite loci uncovered substantial differentiation between lineages at homologous flanking regions. Overall our findings indicate a high degree of nuclear genetic differentiation between the two lineages of L. rubellus, implying reproductive isolation at the two study sites and therefore the potential existence of cryptic species. The existence of two cryptic taxa has major implications for the application of L. rubellus as an ecotoxicological sentinel. It may therefore be necessary to consider the lineages as separate taxa during future ecotoxicological studies. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2013, 110 , 780–795.     

Molecular structure of the allelic variants of (AAT)n microsatellite locus Du47D in the parthenogenetic species Darevskia unisexualis and bisexual parental species D. valentini and D. raddei

Microsatellite repeats are one of the most widespread elements of the eukaryotic genome, but are poorly studied in species with clonal reproduction. PCR analysis and DNA sequencing were used to study the molecular structure of the allelic variants of microsatellite locus Du47D in the parthenogenetic species Darevskia unisexualis and its evolutionary ancestors, bisexual species D. raddei and D. valentini, of the genus Darevskia (Lacerta saxicola complex). Sequencing showed that the allelic variants of the D. unisexualis Du47D locus and the alleles of its D. raddei and D. valentini orthologs have a perfect microsatellite cluster structure, differ in number of ATT monomeric units, and have certain species-specific combinations of nucleotide substitutions, deletions, and insertions in the microsatellite-flanking DNA sequences. The Du47D alleles that the parthenogenetic species inherited from D. valentini or from D. raddei were identified.     

Testing the ‘Pleistocene species pump’ in alpine habitats: lineage diversification of flightless ground beetles (Coleoptera: Carabidae: Nebria) in relation to altitudinal zonation

Alpine species often have similar demographic responses to Pleistocene climate changes, but exhibit different spatial patterns of genetic diversity. Using a comparative phylogeographical approach, we examined the factors influencing lineage formation in three alpine carabid beetles of the genus Nebria Latreille inhabiting the California Sierra Nevada. These flightless beetles differ in altitudinal zonation and habitat preferences, but overlap spatially, have limited dispersal capacities and share life history characteristics. Species distribution modelling predicted decreasing population connectivity in relation to increasing altitudinal preferences. Diversity patterns at the cytochrome oxidase subunit I gene revealed north–south genetic structure and recent population growth in all three species. The high‐elevation‐restricted species, Nebria ingens Horn, exhibited a deep phylogeographical split, morphological divergence and evidence of limited, unidirectional gene flow towards the south. This was supported by additional data from three nuclear genes and isolation with migration analysis. Nebria spatulata Van Dyke, inhabiting an intermediate altitudinal range, exhibited fixed morphological differences between northern and southern populations, but showed limited structure. The broadly distributed Nebria ovipennis LeConte showed less structure and lacked morphological variation. Diversification of these Nebria species supports the role of altitudinal zonation in lineage formation and is consistent with the Pleistocene species pump model. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, ?? , ??–??.     

Molecular phylogeny of sea snakes reveals a rapidly diverged adaptive radiation

Evolutionary relationships within and between the marine hydrophiine sea snake groups have been inferred primarily using morphological characters, and two major groups traditionally are recognized. The Aipysurus group comprises nine species in two genera, and the taxonomically chaotic Hydrophis group comprises as many as 40 species, of which 27 are generally allocated to the genus Hydrophis and 13 to ten additional genera. In addition to these two major groups are three putatively ‘primitive’ monotypic genera, Hydrelaps darwiniensis, Ephalophis greyi and Parahydrophis mertoni. The present study investigated the evolutionary relationships of 23 representative species of marine hydrophiines, comprising 15 species from the Hydrophis group, six species from the Aipysurus group, and H. darwiniensis and P. mertoni, to address two broad aims. First, the aim was to provide a robust phylogeny for sea snakes to test previous phylogenetic hypotheses based on morphology, and thus provide some taxonomic stability to the group. Second, there was interest in evaluating the hypothesis that the Hydrophis group might represent a rapidly diverged adaptive radiation. A large mitochondrial DNA data set based on the cytochrome b gene (1080 bp, 401 parsimony informative) and the 16S rRNA gene (510 bp, 57 parsimony informative) was assembled and these data were analysed using parsimony, maximum‐likelihood and Bayesian approaches. All analyses yielded virtually the same optimal tree, confirming that hydrophiine sea snakes comprise at least three lineages. The Aipysurus group formed a strongly supported and well‐resolved monophyletic clade. The Hydrophis group also formed a strongly supported clade; however, resolution among the genera and species was very poor. Hydrelaps darwiniensis and P. mertoni formed a sister clade to the Hydrophis lineage. Our phylogeny was used to test the validity of previous taxonomic and phylogenetic hypotheses, and to demonstrate that the genus Hydrophis is not monophyletic. Genetic diversity relative to phenotypic diversity is four to seven times greater in the Hydrophis lineage compared with the Aipysurus lineage. The topology of our phylogenetic hypothesis, combined with the levels of genetic divergence relative to morphological diversity, demonstrate that the Hydrophis lineage represents a rapidly diverged adaptive radiation. The data are consistent with the hypothesis that this adaptive radiation may be due to historical sea level fluctuations that have isolated populations and promoted speciation. © 2006 The Linnean Society of London, Biological Journal of the Linnean Society, 2006, 89 , 523–539.     

The wheat curl mite Aceria tosichella (Acari: Eriophyoidea) is a complex of cryptic lineages with divergent host ranges: evidence from molecular and plant bioassay data

Aceria tosichella (the wheat curl mite, WCM) is a global pest of wheat and other cereals, causing losses by direct damage, as well as the transmission of plant viruses. The mite is considered to have an unusually wide host range for an eriophyoid species. The present study tested the commonly held assumption that WCM is a single, highly polyphagous species by assessing the host range of genetically distinct lineages of WCM occurring in Poland on different host plants. Genotyping was performed by analyzing nucleotide sequence data from fragments of the mitochondrial cytochrome c oxidase subunit I (COI) and the nuclear D2 region of 28S rDNA. Mean between‐lineage distance estimated using COI data was found to be one order of magnitude greater than the within‐clade lineage and, in some cases, comparable to distances between WCM lineages and a congeneric outgroup species. Host acceptance was tested by quantifying population growth for different WCM mitochondrial (mt)DNA lineages when transferred from source host plants to test plants. These experiments revealed significant differences in host colonization ability between mtDNA lineages, ranging from highly polyphagous to more host‐specific. The present study reveals that WCM is composed of several discrete genetic lineages with divergent host‐acceptance and specificity traits. Genetic variation for host acceptance within A. tosichella s.l. may act as a reproductive barrier between these lineages, most of which had narrow host ranges. Two lineages appear to have high pest potential on cereals, whereas several others appear to specialize on wild grass species. We conclude that WCM is not a homogeneous species comprising polyphagous panmictic populations rather it is a complex of genetically distinct lineages with variable host ranges and therefore variable pest potential. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2013, 109 , 165–180.     

Genetic patterns in the Lathyrus pannonicus complex (Fabaceae) reflect ecological differentiation rather than biogeography and traditional subspecific division

The endangered European relict species complex Lathyrus pannonicus shows distinct morphological variation, reflected by the number of subspecies recognized, and complicated patterns of genetic variation. The traditionally recognized subspecies appear to possess different ecological preferences and disjunct distributions, particularly in the western range of the species. In this study, L. pannonicus was investigated by the correlation of distance matrices based on phytosociological, ecological, molecular and morphological data. Ecological characteristics of selected stands of L. pannonicus throughout Europe were assessed using ‘Ellenberg values’ of all the constituent taxa in phytosociological relevés. Genetic distances were calculated using recently developed methods to analyse high degrees of intra‐individual nuclear‐encoded internal transcribed spacer variability. We found that the remarkable genetic (and morphological) diversity in L. pannonicus could not be explained solely by the fragmentation of the distributional range. Instead, patterns of morphological and genetic differentiation were a reflection of the moisture regime in the sampled stands. Two major lineages could be identified: (1) a lineage adapted to dry conditions (Ellenberg indicator F‐value ≤ 3.5) and (2) a lineage preferring moist conditions (F‐value ≥ 4.5). Although both lineages occurred in close proximity in the Pannonian area, they appeared to be reproductively isolated in general. Further data are needed to determine whether these genetically and ecologically defined lineages, or ecospecies, within the L. pannonicus species complex can be formalized as (Linnaean) species or subspecies. © 2011 The Linnean Society of London, Botanical Journal of the Linnean Society, 2011, 165 , 402–421.     

Phylogeography of the mottled spinefoot Siganus fuscescens: Pleistocene divergence and limited genetic connectivity across the Philippine archipelago

Historical isolation during Pleistocene low sea level periods is thought to have contributed to divergence among marine basin populations across the Coral Triangle. In the Philippine archipelago, populations in the South China Sea, Sulu Sea–inland seas, and Philippine Sea‐Celebes Sea basins might have been partially isolated. Meanwhile, present‐day broadscale oceanographic circulation patterns suggest connectivity between these basins. To evaluate hypotheses regarding the influence of historical and contemporary factors on genetic structure, phylogeographic patterns based on mitochondrial control region sequences for a reef‐associated fish, Siganus fuscescens, were analysed. Three distinct lineages were recovered. One lineage was identified as the morphologically similar species Siganus canaliculatus, while two lineages are monophyletic with S. fuscescens. Clade divergence and demographic expansion in S. fuscescens occurred during the Pleistocene. A strong signal of latitudinal structure was detected (Φ_CT = 0.188), driven by marked differences in clade distribution: one clade is widely distributed (clade A), while a second clade (clade B) has a restricted northern distribution. Regional structure of clade A is consistent with the basin isolation hypothesis (Φ_CT = 0.040) and suggests isolation of the South China Sea (Φ_CT = 0.091). Fine‐scale structure was observed in the South China Sea and south Philippine Sea, while Sulu Sea and inland seas were unstructured. Genetic structure across multiple spatial scales (archipelagic, regional, and fine‐scale within basins) suggests the influence of vicariant barriers and contemporary limits to gene flow in S. fuscescens that may be influenced by oceanographic circulation, geographical distance between available habitats, and latitudinal temperature differences.     

Fuzzy species limits in Mediterranean gorgonians (Cnidaria,Octocorallia): inferences on speciation processes

The study of the interplay between speciation and hybridization is of primary importance in evolutionary biology. Octocorals are ecologically important species whose shallow phylogenetic relationships often remain to be studied. In the Mediterranean Sea, three congeneric octocorals can be observed in sympatry: Eunicella verrucosa, Eunicella cavolini and Eunicella singularis. They display morphological differences and E. singularis hosts photosynthetic Symbiodinium, contrary to the two other species. Two nuclear sequence markers were used to study speciation and gene flow between these species, through network analysis and Approximate Bayesian Computation (ABC). Shared sequences indicated the possibility of hybridization or incomplete lineage sorting. According to ABC, a scenario of gene flow through secondary contact was the best model to explain these results. At the intraspecific level, neither geographical nor ecological isolation corresponded to distinct genetic lineages in E. cavolini. These results are discussed in the light of the potential role of ecology and genetic incompatibilities in the persistence of species limits.     

Phylogeography of the parasitic fly Batrachomyia in the Wet Tropics of north-east Australia, and susceptibility of host frog lineages in a mosaic contact zone

Parasites could differentially impact intraspecific host lineages due to genetic, phenotypic, ecological, or behavioural differences between the lineages, or the development of reproductive isolation between them. Batrachomyia (Diptera: Chloropidae) are flies that exclusively parasitize Australian frogs, and in the Wet Tropics rainforest of north‐east Australia larvae are largely restricted to the green‐eyed tree frog Litoria genimaculata (Anura: Hylidae). This frog species consists of two highly divergent genetic lineages that overlap in two nearby, but independent, contact zones. At one contact zone there has been extensive phenotypic divergence and speciation between the lineages whereas, at the other contact relatively lower levels of phenotypic divergence and reproductive isolation suggest that speciation has not occurred. In the present study we tested: (1) whether the deep phylogeographic divergence between northern and southern host populations is mirrored by congruent genetic structuring in the parasite populations and (2) whether the host lineages are differentially impacted by parasitism. We found that the two divergent frog lineages are parasitized by a single lineage of Batrachomyia, which exhibits strikingly little phylogeographic structuring. We found a significant difference in Batrachomyia prevalence between the host lineages at mixed lineage sites in both contact zones, with the magnitude and direction of this effect being consistent in both. The pattern did not differ between the two contacts even though recent phenotypic divergence and speciation has occurred between the lineages at one contact but not the other. Taken together, this suggests a fundamental difference in susceptibility between the genetically divergent host lineages. Using weight relative to body length as a measure of body condition, we found no differential impact of parasitism on the body condition of each host lineage, and no evidence that parasitism impacts the body condition of the host in general. © 2007 The Linnean Society of London, Biological Journal of the Linnean Society, 2007, 92 , 593–603.     

Phylogeographical analysis of Ligia oceanica (Crustacea: Isopoda) reveals two deeply divergent mitochondrial lineages

Isopods of the species Ligia oceanica are typical inhabitants of the rocky intertidal of the northern European coastline. The aim of this study was to assess the genetic structure of this species using mitochondrial and nuclear sequence data. We analysed partial mitochondrial cytochrome c oxidase subunit I (CO1) and 16S rRNA gene sequence data of 161 specimens collected from ten sites ranging from Spain to Norway. For selected specimens, we also sequenced the hypervariable V7 expansion segment of the nuclear 18S rRNA gene as a supplementary marker. Furthermore, we studied the infection rate of all analysed specimens by the alphaproteobacterium Wolbachia. Our analyses revealed two deeply divergent mitochondrial lineages for Ligia oceanica that probably diverged in the late Pliocene to mid Pleistocene. One lineage comprised specimens from northern populations (‘lineage N’) and one primarily those from France and Spain (‘lineage S’). Distribution patterns of the haplotypes and the genetic distances between both lineages revealed two populations that diverged before the Last Glacial Maximum. Given that we found no homogenization of mitochondrial haplotypes, our present results also reject any influence of Wolbachia on the observed mtDNA variability. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 112 , 16–30.     

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.     

Species delimitation with gene flow: A methodological comparison and population genomics approach to elucidate cryptic species boundaries in Malaysian Torrent Frogs

Accurately delimiting species boundaries is a nontrivial undertaking that can have significant effects on downstream inferences. We compared the efficacy of commonly used species delimitation methods (SDMs) and a population genomics approach based on genomewide single‐nucleotide polymorphisms (SNPs) to assess lineage separation in the Malaysian Torrent Frog Complex currently recognized as a single species (Amolops larutensis). First, we used morphological, mitochondrial DNA and genomewide SNPs to identify putative species boundaries by implementing noncoalescent and coalescent‐based SDMs (mPTP, iBPP, BFD*). We then tested the validity of putative boundaries by estimating spatiotemporal gene flow (fastsimcoal2 , ABBA‐BABA) to assess the extent of genetic isolation among putative species. Our results show that the A. larutensis complex runs the gamut of the speciation continuum from highly divergent, genetically isolated lineages (mean F_st = 0.9) to differentiating populations involving recent gene flow (mean F_st = 0.05; N_m > 5). As expected, SDMs were effective at delimiting divergent lineages in the absence of gene flow but overestimated species in the presence of marked population structure and gene flow. However, using a population genomics approach and the concept of species as separately evolving metapopulation lineages as the only necessary property of a species, we were able to objectively elucidate cryptic species boundaries in the presence of past and present gene flow. This study does not discount the utility of SDMs but highlights the danger of violating model assumptions and the importance of carefully considering methods that appropriately fit the diversification history of a particular system.     

Allozyme Variation Patterns and Multiple Hybridization Origins: Clonal Variation Among Four Sibling Parthenogenetic Caucasian Rock Lizards

Allozyme electrophoresis of four sibling parthenogenetic Caucasian rock lizards Darevskia unisexualis, D.uzzelli, D.sapphirina, and D.bendimahiensis found seven clones and five variable loci. The data supported the hypothesis that D.raddei and D.valentini are the parental species of all four parthenogens. Variation patterns in Darevskia were summarized. Species that originated from a single F₁ typically consisted of one widespread clone with a few rare clones. Species with multiple origins displayed variation only slightly higher than species with a single origin. This is contrary to other genera of parthenogenetic lizards, in which cases massive clonal variations were observed. This revised version was published online in August 2006 with corrections to the Cover Date.     

Ongoing hybridization obscures phylogenetic relationships in the Drosophila subquinaria species complex

Inferring evolutionary relationships among recently diverged lineages is necessary to understand how isolating barriers produce independent lineages. Here, we investigate the phylogenetic relationships between three incompletely isolated and closely related mushroom‐feeding Drosophila species. These species form the Drosophila subquinaria species complex and consist of one Eurasian species (D. transversa) and two widespread North American species (D. subquinaria and D. recens) that are sympatric in central Canada. Although patterns of pre‐ and post‐mating isolation among these species are well characterized, previous work on their phylogenetic relationships is limited and conflicting. In this study, we generated a multi‐locus data set of 29 loci from across the genome sequenced in a population sample from each species, and then, we inferred species relationships and patterns of introgression. We find strong statistical support that D. subquinaria is paraphyletic, showing that samples from the geographic region sympatric with D. recens are most closely related to D. recens, whereas samples from the geographic region allopatric with D. recens are most closely related to D. transversa. We present several lines of evidence that both incomplete lineage sorting and gene flow are causing phylogenetic discordance. We suggest that ongoing gene flow primarily from D. recens into D. subquinaria in the sympatric part of their ranges causes phylogenetic uncertainty in the evolutionary history of these species. Our results highlight how population genetic data can be used to disentangle the sources of phylogenetic discordance among closely related species.     

Divergence,gene flow,and speciation in eight lineages of trans‐Beringian birds

Determining how genetic diversity is structured between populations that span the divergence continuum from populations to biological species is key to understanding the generation and maintenance of biodiversity. We investigated genetic divergence and gene flow in eight lineages of birds with a trans‐Beringian distribution, where Asian and North American populations have likely been split and reunited through multiple Pleistocene glacial cycles. Our study transects the speciation process, including eight pairwise comparisons in three orders (ducks, shorebirds and passerines) at population, subspecies and species levels. Using ultraconserved elements (UCEs), we found that these lineages represent conditions from slightly differentiated populations to full biological species. Although allopatric speciation is considered the predominant mode of divergence in birds, all of our best divergence models included gene flow, supporting speciation with gene flow as the predominant mode in Beringia. In our eight lineages, three were best described by a split‐migration model (divergence with gene flow), three best fit a secondary contact scenario (isolation followed by gene flow), and two showed support for both models. The lineages were not evenly distributed across a divergence space defined by gene flow (M) and differentiation (F_ST), instead forming two discontinuous groups: one with relatively shallow divergence, no fixed single nucleotide polymorphisms (SNPs), and high rates of gene flow between populations; and the second with relatively deeply divergent lineages, multiple fixed SNPs, and low gene flow. Our results highlight the important role that gene flow plays in avian divergence in Beringia.     

The modality of nine Triturus newt hybrid zones assessed with nuclear,mitochondrial and morphological data

The shape of hybrid zones provides insight into genetic isolation between the evolutionary lineages involved: the greater the limitation to gene flow, the farther along the unimodal to bimodal continuum. We study hybrid zone modality in a group of closely related species at a variety of levels of genetic divergence. We explore the degree and the geography of hybridization in 12 transects of Triturus newts, for nine of the ten secondary contact zones in Europe, using nuclear and mitochondrial encoded genetic data and morphological characteristics. Species status is verified by examining the correlation between the three marker systems in a geographical context. At 40 allozyme loci, two transects classify as unimodal, three as bimodal and two as intermediate, while the signal for hybridization is weak or absent in five transects. One zone studied in duplicate was classified as intermediate in one region and as bimodal in another region. mtDNA introgression is frequent and extends beyond nuclear introgression in two transects. Morphology provides additional evidence for hybridity, including one transect for which the signal of nuclear gene flow is weak. Compared to simulations allowing panmixia, the observed allozyme transitions at contact zones show a deficit of backcrossing to various degrees. Over all transects, there is a weak negative relationship between the level of hybridization and allozyme genetic distance for species pairs, consistent with Bateson–Dobzhansky–Muller effects. This observation, based upon highly comparable data for a single genus, supports conclusions derived from analyses over a wide variety of other taxa. © 2014 The Authors. Biological Journal of the Linnean Society published by John Wiley & Sons Ltd on behalf of The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 113 , 604–622.