Hybridization and the inheritance of female colour polymorphism in two ischnurid damselflies (Odonata: Coenagrionidae)

Female‐limited polychromatism is frequent in many species of Odonata. Ischnura elegans has three colour morphs: one male‐like coloured (androchrome) and two additional gynochrome brown morphs (infuscans and rufescens‐obsoleta morphs). A total of 19 progenies obtained from once‐mated females were reared in the laboratory in three generations. Results indicate that the colour morphs are controlled by the same genetic system as previously described for I. graellsii, i.e. an autosomal locus with female‐limited expression and with three alleles with a hierarchy of dominance (p^a > pⁱ > p°). Five interspecific crossings between female I. graellsii and male I. elegans, five crossings between hybrid females and male I. elegans and one crossing between female I. graellsii and a hybrid male further confirmed that the genetic system is the same in both species. A survey of morph frequencies in north‐west Spain revealed that I. elegans shows high variability in androchrome frequency (4–91%) between nearby populations, whereas in I. graellsii androchromes never are the majority morph (5–40%). The highest androchrome frequency in I. graellsii was found in populations closest to a locality where both species have hybridized, and that now has the highest androchrome frequency of I. elegans. We hypothesize that I. elegans genes have been incorporated into the genome of I. graellsii resulting in increased androchrome frequency in the latter species. Low androchrome frequency in I. elegans seems also related to the influence of I. graellsii genes. Therefore, we suggest that hybridization between both taxa is contributing to the temporal maintenance of contrasting androchrome frequencies in nearby populations. © 2005 The Linnean Society of London, Biological Journal of the Linnean Society, 2005, 85 , 471–481.     

Environmental (in)dependence of a hybrid zone: Insights from molecular markers and ecological niche modeling in a hybrid zone of Origanum (Lamiaceae) on the island of Crete

The role of environment and the relative significance of endogenous versus exogenous selection in shaping hybrid zones have been crucial issues in the studies of hybridization. Recent advances in ecological niche modeling (ENM) offer new methodological tools, especially in combination with the genotyping of individuals in the hybrid zone. Here, we study the hybrid zone between the widely known spices Origanum onites and Origanum vulgare ssp. hirtum in Crete. We analyze the genetic structure of both parental taxa and their hybrid Origanum × intercendens using AFLP markers on 15 sympatric and 12 allopatric populations and employ ecological niche modeling and niche similarity tests to study their niche patterns. We complement these analyses with seed viability measurements. Our study revealed that the hybridizing taxa O. onites and O. vulgare ssp. hirtum and the resulting genotypic classes showed geographical and environmental niche similarities based on the predictions of ENMs and the subsequent similarity tests. The occurrence of the hybrid zone is not directly dependent on environmental factors which favor the fitness of the hybrid compared to the parental taxa, but rather on aspects such as historical factors and management practices, which may contribute to the localization and maintenance of the contact zone between parental species. Our results suggest that if a minimum required niche differentiation between genotypic classes is not achieved, environmental dependence might not have a prominent role on the outcome of the hybridization.     

Gene flow,divergent selection and resistance to introgression in two species of morning glories (Ipomoea)

Gene flow is thought to impede genetic divergence and speciation by homogenizing genomes. Recent theory and research suggest that sufficiently strong divergent selection can overpower gene flow, leading to loci that are highly differentiated compared to others. However, there are also alternative explanations for this pattern. Independent evidence that loci in highly differentiated regions are under divergent selection would allow these explanations to be distinguished, but such evidence is scarce. Here, we present multiple lines of evidence that many of the highly divergent SNPs in a pair of sister morning glory species, Ipomoea cordatotriloba and I. lacunosa, are the result of divergent selection in the face of gene flow. We analysed a SNP data set across the genome to assess the amount of gene flow, resistance to introgression and patterns of selection on loci resistant to introgression. We show that differentiation between the two species is much lower in sympatry than in allopatry, consistent with interspecific gene flow in sympatry. Gene flow appears to be substantially greater from I. lacunosa to I. cordatotriloba than in the reverse direction, resulting in sympatric and allopatric I. cordatotriloba being substantially more different than sympatric and allopatric I. lacunosa. Many SNPs highly differentiated in allopatry have experienced divergent selection, and, despite gene flow in sympatry, resist homogenization in sympatry. Finally, five out of eight floral and inflorescence characteristics measured exhibit asymmetric convergence in sympatry. Consistent with the pattern of gene flow, I. cordatotriloba traits become much more like those of I. lacunosa than the reverse. Our investigation reveals the complex interplay between selection and gene flow that can occur during the early stages of speciation.     

A biogeographic genetic approach for testing the role of reinforcement: the case of Drosophila pseudoobscura and D. persimilis

Abstract.— The role of reinforcement in speciation can be explained by two distinct models. In model I, two diverged populations hybridize and produce fertile hybrids that successfully backcross (hybridization with gene flow). In model II, two populations hybridize but succeeding backcrosses are unproductive (hybridization without gene flow). Using Drosophila persimilis and D. pseudoobscura , we have tested model I by comparing the extent of heterospecific introgression in sympatric versus allopatric populations. We show that certain expectations of this particular model of reinforcement, which is based on hybridization and gene flow between divergent populations after secondary contact, are not realized in these two species. The evidence consists of the similarity of genetic distances as well as proportions of unique/rare alleles between sympatric and allopatric heterospecific populations and a negative correlation between genetic distance and geographical distance between heterospecific populations, which suggests ecological differentiation. This approach in quantifying differential gene flow has important consequences to studies that compare sympatric and allopatric isolation using genetic distance. Following model I, one would expect a pattern of higher prezygotic isolation in sympatric species compared to allopatric species of the same genetic distance simply as a result of an underestimation of genetic distance due to introgression between sympatric populations. We suggest more parsimonious explanations such as reinforcement without genetic exchange (model II) and ecological differentiation, which require high levels of preexisting reproductive isolation between populations.     

Rapid evolution of prezygotic barriers in non‐territorial damselflies

A central question in evolutionary biology concerns the accumulation of reproductive barriers during speciation. However, separating the reproductive barriers that have led to speciation from those that have secondarily accumulated (i.e. after initial divergence) is a widely recognized problem. Ideal candidate species for overcoming this problem are young species, where time for additional barriers to accrue has been limited. In the present study, we add to previous studies investigating the strength of reproductive barriers between the parapatric damselflies Ischnura elegans and Ischnura graellsii by quantifying seven prezygotic barriers between the allopatric pairs of I. elegans and Ischnura genei, as well as I. graellsii and I. genei. Specifically, we measured four premating (temporal, sexual, mechanical _I, and mechanical _II) and three postmating (oviposition success, fecundity, and fertility) barriers using experimental approaches and, for first time, we investigated the mechanisms causing mechanical isolation, which is the strongest reproductive barrier in ischnurans. The findings of the present study support the notion that premating barriers are generally strong and contribute significantly to total reproductive isolation in young lineages (65–98%), although they never solely lead to complete isolation. Asymmetry was generally stronger in premating than in postmating barriers, and was driven mostly through asymmetry in mechanical isolation, which is caused by morphological divergence of secondary sexual appendages. We found that barriers act multiplicatively in all species combinations tested, with the exception of sexual isolation, which was not detected. Our results are consistent with a recent allopatric speciation scenario driven by differences in male anal appendages, either impeding copulation or affecting female preferences. Taken together, the results from this and previous studies in diverse odonate genera suggest that premating barriers have evolved rapidly in ischnuran damselflies and, although reproductive isolation in ischnurans is more commonly the result of several barriers acting together, morphological divergence of secondary sexual appendages appears to be a common factor facilitating premating isolation in this group. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 113 , 485–496.     

Comparative phylogeography and asymmetric hybridization between cryptic bat species

Cryptic speciation and hybridization are two key processes that affect the origin and maintenance of biodiversity and our ability to understand and estimate it. To determine how these two processes interact, we studied allopatric and sympatric colonies of two cryptic bat species (Eptesicus serotinus and Eptesicus isabellinus) with parapatric distribution in the Iberian Peninsula. These species are the main reservoir for the most commonly rabies virus found in bats in Europe: the European bat Lyssavirus type 1 (EBLV‐1). We used mtDNA and nuclear microsatellite markers to confirm the taxonomic status of both species and to show a more pronounced and geographically based genetic structure in E. isabellinus than in its sibling E. serotinus. Using approximate Bayesian computation (ABC), we inferred rapid range expansion in both species after the Last Glacial Maximum until reaching their present distributions. ABC analysis also supported interspecific differences in genetic diversity and structure, pointing to an earlier expansion of E. isabellinus northward. We found no evidence of mitochondrial introgression between species, but nuclear markers identified a male‐mediated ongoing asymmetric hybridization from E. isabellinus to E. serotinus (28% hybrids in E. serotinus and 5% in E. isabellinus) in the contact zone. Although none of the bats studied tested positive for Lyssavirus RNA, the asymmetric hybridization supports the potential for the recently suggested interspecific transmission of EBLV‐1 from E. isabellinus into E. serotinus.     

Phylogeny and biogeography of the American live oaks (Quercus subsection Virentes): a genomic and population genetics approach

The nature and timing of evolution of niche differentiation among closely related species remains an important question in ecology and evolution. The American live oak clade, Virentes, which spans the unglaciated temperate and tropical regions of North America and Mesoamerica, provides an instructive system in which to examine speciation and niche evolution. We generated a fossil‐calibrated phylogeny of Virentes using RADseq data to estimate divergence times and used nuclear microsatellites, chloroplast sequences and an intron region of nitrate reductase (NIA‐i3) to examine genetic diversity within species, rates of gene flow among species and ancestral population size of disjunct sister species. Transitions in functional and morphological traits associated with ecological and climatic niche axes were examined across the phylogeny. We found the Virentes to be monophyletic with three subclades, including a southwest clade, a southeastern US clade and a Central American/Cuban clade. Despite high leaf morphological variation within species and transpecific chloroplast haplotypes, RADseq and nuclear SSR data showed genetic coherence of species. We estimated a crown date for Virentes of 11 Ma and implicated the formation of the Sea of Cortés in a speciation event ~5 Ma. Tree height at maturity, associated with fire tolerance, differs among the sympatric species, while freezing tolerance appears to have diverged repeatedly across the tropical–temperate divide. Sympatric species thus show evidence of ecological niche differentiation but share climatic niches, while allopatric and parapatric species conserve ecological niches, but diverge in climatic niches. The mode of speciation and/or degree of co‐occurrence may thus influence which niche axis plants diverge along.     

Spatial and host related genomic variation in partially sympatric cactophagous moth species

Surveys of patterns of genetic variation in natural sympatric and allopatric populations of recently diverged species are necessary to understand the processes driving intra- and interspecific diversification. The South American moths Cactoblastis cactorum, Cactoblastis doddi and Cactoblastis bucyrus are specialized in the use of cacti as host plants. These species have partially different geographic ranges and differ in patterns of host plant use. However, there are areas that overlap, particularly, in northwestern Argentina, where they are sympatric. Using a combination of genome-wide SNPs and mitochondrial data we assessed intra and interspecific genetic variation and investigated the relative roles of geography and host plants on genetic divergence. We also searched for genetic footprints of hybridization between species. We identified three well delimited species and detected signs of hybridization in the area of sympatry. Our results supported a hypothetical scenario of allopatric speciation in the generalist C. cactorum and genetic interchange during secondary geographic contact with the pair of specialists C. bucyrus and C. doddi that probably speciated sympatrically. In both cases, adaptation to new host plants probably played an important role in speciation. The results also suggested the interplay of geography and host plant use as drivers of divergence and limiting gene flow at intra and interspecific levels.     

Extensive hybridization and past introgression between divergent lineages in a quasi‐clonal hermaphroditic fish: Ramifications for species concepts and taxonomy

Extreme inbreeding is expected to reduce the incidence of hybridization, serving as a prezygotic barrier. Mangrove rivulus is a small killifish that reproduces predominantly by self‐fertilization, producing highly homozygous lines throughout its geographic range. The Bahamas and Caribbean are inhabited by two highly diverged phylogeographic lineages of mangrove rivulus, Kryptolebias marmoratus and a ‘Central clade’ closely related to K. hermaphroditus from Brazil. The two lineages are largely allopatric, but recently were found in syntopy on San Salvador, Bahamas, where a single hybrid was reported. To better characterize the degree of hybridization and the possibility of secondary introgression, here we conducted a detailed genetic analysis of the contact zone on San Salvador. Two mixed populations were identified, one of which contained sexually mature hybrids. The distribution of heterozygosity at diagnostic microsatellite loci in hybrids showed that one of these hybrids was an immediate offspring from the K. marmoratus x Central clade cross, whereas the remaining five hybrids were products of reproduction by self‐fertilization for 1–3 generations following the initial cross. Two hybrids had mitochondrial haplotypes of K. marmoratus and the remaining four hybrids had a haplotype of the Central clade, indicating that crosses go in both directions. In hybrids, alleles of parental lineages were represented in equal proportions suggesting lack of recent backcrossing to either of the parental lineages. However, sympatric populations of two lineages were less diverged than allopatric populations, consistent with introgression. Results are discussed in terms of applicability of the biological species concept for isogenic, effectively clonal, organisms.     

Maintenance of species boundaries in a Neotropical radiation of Begonia

A major goal of evolutionary biology is to determine the mechanisms generating biodiversity. In Begonia, one of the largest plant genera (1900+ species), it has been postulated that the high number of endemic species is a by‐product of low gene flow among populations, which predisposes the group to speciation. However, this model of divergence requires that reproductive barriers accumulate rapidly among diverging species that overlap in their geographic ranges, otherwise speciation will be opposed by homogenizing gene flow in zones of secondary contact. Here, we test the outcomes of secondary contact in Begonia by genotyping multiple sympatric sites with 12 nuclear and seven plastid loci. We show that three sites of secondary contact between B. heracleifolia and B. nelumbiifolia are highly structured, mostly containing parental genotypes, with few F1 hybrids. A sympatric site between B. heracleifolia and B. sericoneura contains a higher proportion of F1s, but little evidence of introgression. The lack of later‐generation hybrids contrasts with that documented in many other plant taxa, where introgression is extensive. Our results, in conjunction with previous genetic work, show that Begonia demonstrate properties making them exceptionally prone to speciation, at multiple stages along the divergence continuum. Not only are populations weakly connected by gene flow, promoting allopatric speciation, but species often show strong reproductive barriers in secondary contact. Whether similar mechanisms contribute to diversification in other large genera remains to be tested.     

Evidence of ecological niche shift in Rhododendron ponticum (L.) in Britain: Hybridization as a possible cause of rapid niche expansion

Biological invasions threaten global biodiversity and natural resources. Anticipating future invasions is central to strategies for combating the spread of invasive species. Ecological niche models are thus increasingly used to predict potential distribution of invasive species. In this study, we compare ecological niches of Rhododendron ponticum in its native (Iberian Peninsula) and invasive (Britain) ranges. Here, we test the conservation of ecological niche between invasive and native populations of R. ponticum using principal component analysis, niche dynamics analysis, and MaxEnt‐based reciprocal niche modeling. We show that niche overlap between native and invasive populations is very low, leading us to the conclusion that the two niches are not equivalent and are dissimilar. We conclude that R. ponticum occupies novel environmental conditions in Britain. However, the evidence of niche shift presented in this study should be treated with caution because of nonanalogue climatic conditions between native and invasive ranges and a small population size in the native range. We then frame our results in the context of contradicting genetic evidence on possible hybridization of this invasive species in Britain. We argue that the existing contradictory studies on whether hybridization caused niche shift in R. ponticum are not sufficient to prove or disprove this hypothesis. However, we present a series of theoretical arguments which indicate that hybridization is a likely cause of the observed niche expansion of R. ponticum in Britain.     

Maintaining their genetic distance: Little evidence for introgression between widely hybridizing species of Geum with contrasting mating systems

Within the plant kingdom, many genera contain sister lineages with contrasting outcrossing and inbreeding mating systems that are known to hybridize. The evolutionary fate of these sister lineages is likely to be influenced by the extent to which they exchange genes. We measured gene flow between outcrossing Geum rivale and selfing Geum urbanum, sister species that hybridize in contemporary populations. We generated and used a draft genome of G. urbanum to develop dd‐RAD data scorable in both species. Coalescent analysis of RAD data from allopatric populations indicated that the species diverged 2–3 Mya, and that historical gene flow between them was extremely low (1 migrant every 25 generations). Comparison of genetic divergence between species in sympatry and allopatry, together with an analysis of allele frequencies in potential parental and hybrid populations, provided no evidence of contemporary introgression in sympatric populations. Cluster‐ and species‐specific marker analyses revealed that, apart from four early‐generation hybrids, individuals in sympatric populations fell into two genetically distinct groups that corresponded exactly to their morphological species classification with maximum individual admixture estimates of only 1–3%. However, we did observe joint segregation of four putatively introgressed SNPs across two scaffolds in the G. urbanum population that was associated with significant morphological variation, interpreted as tentative evidence for rare, recent interspecific gene flow. Overall, our results indicate that despite the presence of hybrids in contemporary populations, genetic exchange between G. rivale and G. urbanum has been extremely limited throughout their evolutionary history.     

Unexpected cryptic species diversity in the widespread coral Seriatopora hystrix masks spatial‐genetic patterns of connectivity

Mounting evidence of cryptic species in a wide range of taxa highlights the need for careful analyses of population genetic data sets to unravel within‐species diversity from potential interspecies relationships. Here, we use microsatellite loci and hierarchical clustering analysis to investigate cryptic diversity in sympatric and allopatric (separated by 450 km) populations of the widespread coral Seriatopora hystrix on the Great Barrier Reef. Structure analyses delimited unique genetic clusters that were confirmed by phylogenetic and extensive population‐level analyses. Each of four sympatric yet distinct genetic clusters detected within S. hystrix demonstrated greater genetic cohesion across regional scales than between genetic clusters within regions (<10 km). Moreover, the magnitude of genetic differentiation between different clusters (>0.620 G”_ST) was similar to the difference between S. hystrix clusters and the congener S. caliendrum (mean G”_ST 0.720). Multiple lines of evidence, including differences in habitat specificity, mitochondrial identity, Symbiodinium associations and morphology, corroborate the nuclear genetic evidence that these distinct clusters constitute different species. Hierarchical clustering analysis combined with more traditional population genetic methods provides a powerful approach for delimiting species and should be regularly applied to ensure that ecological and evolutionary patterns interpreted for single species are not confounded by the presence of cryptic species.     

Speciation processes in putative island endemic sister bat species: false impressions from mitochondrial DNA and microsatellite data

Cases of geographically restricted co‐occurring sister taxa are rare and may point to potential divergence with gene flow. The two bat species Murina gracilis and Murina recondita are both endemic to Taiwan and are putative sister species. To test for nonallopatric divergence and gene flow in these taxa, we generated sequences using Sanger and next‐generation sequencing, and combined these with microsatellite data for coalescent‐based analyses. MtDNA phylogenies supported the reciprocally monophyletic sister relationship between M. gracilis and M. recondita; however, clustering of microsatellite genotypes revealed several cases of species admixture suggesting possible introgression. Sequencing of microsatellite flanking regions revealed that admixture signatures stemmed from microsatellite allele homoplasy rather than recent introgressive hybridization, and also uncovered an unexpected sister relationship between M. recondita and the continental species Murina eleryi, to the exclusion of M. gracilis. To dissect the basis of these conflicts between ncDNA and mtDNA, we analysed sequences from 10 anonymous ncDNA loci with *beast and isolation‐with‐migration and found two distinct clades of M. eleryi, one of which was sister to M. recondita. We conclude that Taiwan was colonized by the ancestor of M. gracilis first, followed by the ancestor of M. recondita after a period of allopatric divergence. After colonization, the mitochondrial genome of M. recondita was replaced by that of the resident M. gracilis. This study illustrates how apparent signatures of sympatric divergence can arise from complex histories of allopatric divergence, colonization and hybridization, thus highlighting the need for rigorous analyses to distinguish between such scenarios.     

Novel trophic niches drive variable progress towards ecological speciation within an adaptive radiation of pupfishes

Adaptive radiation is recognized by a rapid burst of phenotypic, ecological and species diversification. However, it is unknown whether different species within an adaptive radiation evolve reproductive isolation at different rates. We compared patterns of genetic differentiation between nascent species within an adaptive radiation of Cyprinodon pupfishes using genotyping by sequencing. Similar to classic adaptive radiations, this clade exhibits rapid morphological diversification rates and two species are novel trophic specialists, a scale‐eater and hard‐shelled prey specialist (durophage), yet the radiation is <10 000 years old. Both specialists and an abundant generalist species all coexist in the benthic zone of lakes on San Salvador Island, Bahamas. Based on 13 912 single‐nucleotide polymorphisms (SNPs), we found consistent differences in genetic differentiation between each specialist species and the generalist across seven lakes. The scale‐eater showed the greatest genetic differentiation and clustered by species across lakes, whereas durophage populations often clustered with sympatric generalist populations, consistent with parallel speciation across lakes. However, we found strong evidence of admixture between durophage populations in different lakes, supporting a single origin of this species and genome‐wide introgression with sympatric generalist populations. We conclude that the scale‐eater is further along the speciation‐with‐gene‐flow continuum than the durophage and suggest that different adaptive landscapes underlying these two niche environments drive variable progress towards speciation within the same habitat. Our previous measurements of fitness surfaces in these lakes support this conclusion: the scale‐eating fitness peak may be more distant than the durophage peak on the complex adaptive landscape driving adaptive radiation.     

Genomic signatures of divergent selection and speciation patterns in a ‘natural experiment’, the young parallel radiations of Nicaraguan crater lake cichlid fishes

Divergent selection is the main driving force in sympatric ecological speciation and may also play a strong role in divergence between allopatric populations. Characterizing the genome‐wide impact of divergent selection often constitutes a first step in unravelling the genetic bases underlying adaptation and ecological speciation. The Midas cichlid fish (Amphilophus citrinellus) species complex in Nicaragua is a powerful system for studying evolutionary processes. Independent colonizations of isolated young crater lakes by Midas cichlid populations from the older and great lakes of Nicaragua resulted in the repeated evolution of adaptive radiations by intralacustrine sympatric speciation. In this study we performed genome scans on two repeated radiations of crater lake species and their great lake source populations (1030 polymorphic AFLPs, n ～ 30 individuals per species). We detected regions under divergent selection (0.3% in the crater lake Xiloá flock and 1.7% in the older crater lake Apoyo radiation) that might be responsible for the sympatric diversifications. We find no evidence that the same genomic regions have been involved in the repeated evolution of parallel adaptations across crater lake flocks. However, there is some genetic parallelism apparent (seven out of 51 crater lake to great lake outlier loci are shared; 13.7%) that is associated with the allopatric divergence of both crater lake flocks. Interestingly, our results suggest that the number of outlier loci involved in sympatric and allopatric divergence increases over time. A phylogeny based on the AFLP data clearly supports the monophyly of both crater lake species flocks and indicates a parallel branching order with a primary split along the limnetic‐benthic axis in both radiations.     

Evidence for ecological speciation via a host shift in the holly leaf miner,Phytomyza glabricola (Diptera: Agromyzidae)

Evolutionary radiations have been well documented in plants and insects, and natural selection may often underly these radiations. If radiations are adaptive, the diversity of species could be due to ecological speciation in these lineages. Agromyzid flies exhibit patterns of repeated host‐associated radiations. We investigated whether host‐associated population divergence and evidence of divergent selection exist in the leaf miner Phytomyza glabricola on its sympatric host plants, the holly species, Ilex coriacea and I. glabra. Using AFLPs and nuclear sequence data, we found substantial genetic divergence between host‐associated populations of these flies throughout their geographic range. Genome scans using the AFLP data identified 13 loci under divergent selection, consistent with processes of ecological speciation. EF‐1α data suggest that I. glabra is the original host of P. glabricola and that I. coriacea is the novel host, but the AFLP data are ambiguous with regard to directionality of the host shift.     

Directional mitochondrial introgression and character displacement due to reproductive interference in two closely related Pterostichus ground beetle species

Reproductive interference due to interspecific hybridization can lead to character displacement among related species with overlapping ranges. However, no studies have examined which reproductive traits are most important in reducing reproductive interference. We conducted molecular analyses of two nuclear genes (28S and Wingless) and a mitochondrial gene (COI) from two closely related ground beetle species, Pterostichus thunbergi and Pterostichus habui (Coleoptera: Carabidae), with overlapping distributions. In addition, we examined four reproductive traits (body size, organ morphologies of intromittent and non‐intromittent male genital organs, and female reproductive period) in sympatric and allopatric habitats. We compared male genital morphology using geometric morphometric analysis. The species determined by morphology were classified into separate groups based on the phylogenetic tree constructed by the nuclear gene (Wingless). However, according to the mitochondrial genes examined, P. thunbergi was not monophyletic, whereas at the sympatric sites, these species formed a monophyletic clade. This incongruence suggests that interspecific hybridization and subsequent mitochondrial introgression from P. habui to P. thunbergi have occurred. Concerning genital morphology, both of the intromittent and nonintromittent organs of P. thunbergi differed more from P. habui at the sympatric sites than between allopatric sites, suggesting reproductive character displacement. Pterostichus thunbergi, which likely arrived in P. habui habitat in small numbers, would have experienced stronger selection pressures than P. habui.