Hybridization and speciation in the carabid beetles of the subgenusOhomopterus (Coleoptera, Carabidae, genusCarabus)

Evolutionary history and diversity in the ball roller beetle Canthon cyanellus

To understand the evolutionary history of species, it is necessary to know the mechanisms for reproductive isolation, divergence-time between populations, and the relative action of the evolutionary forces (e.g., mutation, genetic drift, gene flow) within and between populations of the same, or closely related species. Although Canthon is one of the more diverse genera of neotropical beetles, insufficient research has been done to comprehend the divergent patterns that explain its speciation process. The absence of diagnostic morphological characters and the wide geographic variation of qualitative traits in Scarabaeinae obscures species delimitation, genealogical limits between populations, and its taxonomy. Canthon cyanellus is one of the best-known species in ecological and evolutionary aspects. It is a widely distributed species in the tropical forests of America. Also, the current deforestation has facilitated its incursion into open areas. Individuals from different populations have similar morphological characters but show wide variation in body color throughout their distribution, which makes it difficult to delimit the subspecies that comprise it. Recently, studies have been carried out to elucidate the pre-and postzygotic isolation mechanisms between populations and the historical biogeographical processes favoring cladogenesis events during the Pleistocene. Morphological variation of the male genitalia does not correspond to the phylogeographic structure. However, the morphological differences in one of the pieces of the endophallic sclerites have allowed a preliminary delimitation of some genetically differentiated clades. Finally, we consider that the joint analysis of traditional morphological taxonomy and phylogeography is important to understand the speciation process in the C. cyanellus complex.

Reproductive isolation via divergent genital morphology due to cascade reinforcement in Ohomopterus ground beetles

Secondary contact between incipient species and selection against maladaptive hybridization can drive reinforcement between populations in contact and result in reproductive character displacement (RCD). Resultant divergence in mating traits within a species may generate downstream reproductive isolation between populations with displaced and non-displaced traits, referred to as the cascade reinforcement hypothesis. We examined this hypothesis using three allopatric populations of the ground beetle Carabus maiyasanus with a genital lock-and-key system. This species shows RCD in male and female genital morphologies in populations in contact with the sister species C. iwawakianus. In a reciprocal mating experiment using three allopatric populations with differences in male and female genital sizes, insemination failure increased as the difference in genital size increased. Based on the reproductive isolation index, insemination failure was the major postmating-prezygotic isolation barrier, at least in one population pair with comparable total isolation to those of other species pairs. By contrast, there was only incomplete premating isolation among populations. These results suggest that RCD in genital morphologies drives incipient allopatric speciation, supporting the cascade reinforcement hypothesis. These findings provide insight into the roles of interspecific interactions and subsequent trait diversification in speciation processes.

Relationships between reproductive character displacement in genital morphology and the population-level cost of interspecific mating: implications for the Templeton effect

Natural selection against maladaptive interspecific reproductive interactions may cause greater divergence in mating traits between sympatric populations than between allopatric populations in a pair of species, known as reproductive character displacement (RCD) which is evidence for the lock-and-key hypothesis of genital evolution. However, the relative importance of various processes contributing to RCD in genital morphology (e.g. reinforcement, reproductive interference, and population filtering or the Templeton effect) is not clear. Here, we examined hypotheses for RCD in genital morphology, with a special focus on the Templeton effect (which predicts that only highly differentiated populations can exist in sympatry). We examined population-level fitness costs in interspecific mating between Carabus maiyasanus and Carabus iwawakianus with RCD in genital morphology. A mating experiment using populations with various degrees of RCD in genital morphology showed no evidence for consistently lower interspecific mating costs in C. maiyasanus populations in contact with displacement in genital morphology than in remote populations, contrary to the predictions of the Templeton effect. Alternatively, interspecific mating costs varied among populations. Observed relationships between the sizes of genital parts concerning isolation and interspecific mating costs across populations suggested that population-level fitness costs do not necessarily decrease during the process leading to RCD. Our results provide insight into ecological and evolutionary processes during secondary contact in closely related species.

Local ecological divergence of two closely related stag beetles based on genetic, morphological, and environmental analyses

The process of phenotypic adaptation to the environments is widely recognized. However, comprehensive studies integrating phylogenetic, phenotypic, and ecological approaches to assess this process are scarce. Our study aims to assess whether local adaptation may explain intraspecific differentiation by quantifying multidimensional differences among populations in closely related lucanid species, Platycerus delicatulus and Platycerus kawadai, which are endemic saproxylic beetles in Japan. First, we determined intraspecific analysis units based on nuclear and mitochondrial gene analyses of Platycerus delicatulus and Platycerus kawadai under sympatric and allopatric conditions. Then, we compared differences in morphology and environmental niche between populations (analysis units) within species. We examined the relationship between morphology and environmental niche via geographic distance. P. kawadai was subdivided into the “No introgression” and “Introgression” populations based on mitochondrial COI gene – nuclear ITS region discordance. P. delicatulus was subdivided into “Allopatric” and “Sympatric” populations. Body length differed significantly among the populations of each species. For P. delicatulus, character displacement was suggested. For P. kawadai, the morphological difference was likely caused by geographic distance or genetic divergence rather than environmental differences. The finding showed that the observed mitochondrial–nuclear discordance is likely due to historical mitochondrial introgression following a range of expansion. Our results show that morphological variation among populations of P. delicatulus and P. kawadai reflects an ecological adaptation process based on interspecific interactions, geographic distance, or genetic divergence. Our results will deepen understanding of ecological specialization processes across the distribution and adaptation of species in natural systems.

The Syngameon Enigma

Despite their evolutionary relevance, multispecies networks or syngameons are rarely reported in the literature. Discovering how syngameons form and how they are maintained can give insight into processes such as adaptive radiations, island colonizations, and the creation of new hybrid lineages. Understanding these complex hybridization networks is even more pressing with anthropogenic climate change, as syngameons may have unique synergistic properties that will allow participating species to persist. The formation of a syngameon is not insurmountable, as several ways for a syngameon to form have been proposed, depending mostly on the magnitude and frequency of gene flow events, as well as the relatedness of its participants. Episodic hybridization with small amounts of introgression may keep syngameons stable and protect their participants from any detrimental effects of gene flow. As genomic sequencing becomes cheaper and more species are included in studies, the number of known syngameons is expected to increase. Syngameons must be considered in conservation efforts as the extinction of one participating species may have detrimental effects on the survival of all other species in the network.

Topology of syngameons

Syngameons are sets of species linked by interspecific hybridization. Common observations regarding the structure of syngameons are that hybridization propensity is not uniform across species and that patterns of hybridization are dominated by a few species. I use computer simulations to test these claims in naturally occurring syngameons selected from the literature and from personal observation. Natural syngameons, especially those involving plants, typically exhibit nonrandom structure: The first three order statistics for the number of hybrid partners and the variance in the number of hybrid partners are larger than chance alone would predict. The structure of two insect syngameons examined is not significantly different from random. To test a hypothesis that variation in hybridization propensity across species in natural syngameons is simply an artifact of hybridization opportunity, I examine the structure of four artificial syngameons (fertility relationships) produced by full diallel crosses. Three of four artificial syngameons exhibit nonrandom structure, as the observed variation in number of successful crosses is larger than chance alone would predict. In general, there are no significant results involving the order statistics. Finally, I discuss biogeographic, ecological, and phylogenetic hypotheses for variation in hybridization propensity across species in natural syngameons.

A Natural Population Derived from Species Hybridizationin the Drosophila ananassae Species Complexon Penang Island, Malaysia

We surveyed natural population of the Drosophila ananassae species complex on Penang Island, Malaysia. Analyses of phenotypic traits, chromosome arrangements, molecular markers, and reproductive isolation suggest the existence of two species: D. ananassae and D. cf. parapallidosa. Molecular marker analysis indicates that D. cf. parapallidosa carries chromosome Y and 4 introgressions from D. ananassae. Thus, D. cf. parapallidosa seems to be a hybrid descendant that recently originated from a natural D. parapallidosa♀× D. ananassae♂ cross. Furthermore, D. cf. parapallidosa behaves differently from authentic D. parapallidosa with respect to its reproductive isolation from D. ananassae. Premating isolation is usually seen in only the D. ananassae♀× D. parapallidosa♂ cross, but we observed it in crosses of both directions between D. ananassae and D. cf. parapallidosa. In addition, hybrid males from the D. ananassae♀× D. parapallidosa♂ cross are usually sterile, but they were fertile when D. ananassae♀ were mated with D. cf. parapallidosa ♂. We attempted an artificial reconstruction of the hybrid species to simulate the evolutionary process(es) that produced D. cf. parapallidosa. This is a rare case of natural hybrid population in Drosophila and may be a useful system for elucidating speciation with gene flow.

A DNA barcode library for ground beetles (Insecta, Coleoptera, Carabidae) of Germany: The genus Bembidion Latreille, 1802 and allied taxa

As molecular identification method, DNA barcoding based on partial cytochrome c oxidase subunit 1 (COI) sequences has been proven to be a useful tool for species determination in many insect taxa including ground beetles. In this study we tested the effectiveness of DNA barcodes to discriminate species of the ground beetle genus Bembidion and some closely related taxa of Germany. DNA barcodes were obtained from 819 individuals and 78 species, including sequences from previous studies as well as more than 300 new generated DNA barcodes. We found a 1:1 correspondence between BIN and traditionally recognized species for 69 species (89%). Low interspecific distances with maximum pairwise K2P values below 2.2% were found for three species pairs, including two species pairs with haplotype sharing (Bembidion atrocaeruleum/Bembidion varicolor and Bembidion guttula/Bembidion mannerheimii). In contrast to this, deep intraspecific sequence divergences with distinct lineages were revealed for two species (Bembidion geniculatum/Ocys harpaloides). Our study emphasizes the use of DNA barcodes for the identification of the analyzed ground beetles species and represents an important step in building-up a comprehensive barcode library for the Carabidae in Germany and Central Europe as well.

Ecological differentiation and habitat unsuitability maintaining a ground beetle hybrid zone

Exogenous selection via interactions between organisms and environments may influence the dynamics of hybrid zones between species in multiple ways. Two major models of a hybrid zone allowed us to hypothesize that environmental conditions influence hybrid zone dynamics in two ways. In the first model, an environmental gradient determines the mosaic distribution at the boundary between ecologically differentiated species (mosaic hybrid zone model). In the second model, a patch of unsuitable habitat traps a hybrid zone between species whose hybrids are unfit (tension zone model). To test these, we examined the environmental factors influencing the spatial structure of a hybrid zone between the ground beetles Carabus maiyasanus and C. iwawakianus using GIS‐based quantification of environmental factors and a statistical comparison of species distribution models (SDMs). We determined that both of the hypothetical processes can be important in the hybrid zone. We detected interspecific differences in the environmental factors in presence localities and their relative contribution in SDMs. SDMs were not identical between species even within contact areas, but tended to be similar within the range of each species. These results suggest an association between environments and species, and provide evidence that ecological differentiation between species plays a role in the maintenance of the hybrid zone. Contact areas were characterized by a relatively high temperature, low precipitation, and high topological wetness. Thus, the contact areas were regarded as being located in an unsuitable habitat with a drier climate, where those populations are likely to occur in patches with limited precipitation concentrated. A comparison of spatial scales suggests that exogenous selection via environmental factors may be weaker than endogenous selection via genitalic incompatibility.

Explosive radiation of the genus Schizopera on a small subterranean island in Western Australia (Copepoda:Harpacticoida): unravelling the cases of cryptic speciation, size differentiation and multiple invasions

A previously unsurveyed calcrete aquifer in the Yilgarn region of Western Australia revealed an unprecedented diversity of copepods, representing 67% of that previously recorded in this whole region. Especially diverse was the genus Schizopera, with up to four morphospecies per bore and a significant size difference between them. Aims of this study were to: (1) survey the extent of this diversity using morphological and molecular tools; (2) derive a molecular phylogeny based on COI; and (3) investigate whether high diversity is a result of an explosive radiation, repeated colonisations, or both, size differentiation is a result of parallel evolution or different phylogeny, and whether Schizopera is a recent invasion in inland waters. More than 300 samples were analysed and the COI fragment successfully amplified by PCR from 43 specimens. Seven species and one subspecies are described as new, and three possible cryptic species were detected. Reconstructed phylogenies reveal that both explosive radiation and multiple colonisations are responsible for this richness, and that Schizopera is probably a recent invasion in these habitats. No evidence for parallel evolution was found, interspecific size differentiation being a result of different phylogeny. Sister species have parapatric distributions and show niche partitioning in the area of overlap.

Resource partitioning or reproductive isolation: the ecological role of body size differences among closely related species in sympatry

1. Body size differences among coexisting related species are common, but the actual effect of these differences in mitigating interspecific interactions, such as resource competition and reproductive interference, is poorly understood. 2. Local assemblages of the ground beetle genus Carabus (subgenus Ohomopterus) typically consist of two or more species of varying sizes. Through foraging and mating experiments using four Ohomopterus species in parapatry and sympatry, we examined whether interspecific body size differences are effective in partitioning food resources or reducing reproductive interference. 3. Because larval Ohomopterus feed exclusively on earthworms, body size differences may be related to partitioning earthworms of different sizes. However, larvae did not exhibit differences in selectivity or attack success on earthworms of different sizes based on larval body size, indicating little possibility of partitioning food by body size. 4. In contrast, interspecific mating behaviours, such as mate recognition, mounting, and copulation, were hindered when body size differences were large; copulation was frequently accomplished between parapatric species with smaller body size differences. 5. These results suggest that body size differences between species effectively reduce reproductive interference, rather than resource competition. Although body size differences in coexisting closely related species have been considered to function in resource partitioning, they may function primarily in reproductive isolation and thereby facilitate coexistence of species.

Reproductive interference determines persistence and exclusion in species interactions

1. Reproductive interference is a negative interspecific sexual interaction that adversely affects the fitness of males and females during reproductive process. Theoretical studies suggest that because reproductive interference is characterized by positive frequency dependence it is far more likely to cause species exclusion than the density dependence of resource competition. However, the respective contributions of resource competition and reproductive interference to species exclusion, which have been frequently observed in many competition studies, remain unclear. 2. We show that reproductive interference is a far more critical cause of species exclusion than resource competition in the competition between Callosobruchus bean weevil species. In competition experiments over several generations, we manipulated the initial relative abundance of the adzuki bean beetle, Callosobruchus chinensis, and the southern cowpea beetle, Callosobruchus maculatus. When the initial adult ratio of C. chinensis : C. maculatus were 6 : 2 and 4 : 4, C. chinensis excluded C. maculatus. However, when C. maculatus was four times more abundant than C. chinensis at the start, we observed the opposite outcome. 3. A behavioural experiment using adults of the two species revealed asymmetric reproductive interference. The fecundity and longevity of C. maculatus females, but not those of C. chinensis females, decreased when the females were kept with heterospecific males. Fecundities of females of both species decreased as the number of heterospecific males increased. In contrast, resource competition at the larval stage resulted in higher survival of C. maculatus than of C. chinensis. 4. These results suggest that the positive frequency-dependent effect of reproductive interference resulted in species exclusion, depending on the initial population ratio of the two species, and the asymmetry of the interference resulted in C. chinensis being dominant in this study, as in previous studies. Classical competition studies should be reviewed in light of this evidence for reproductive interference.

Diversification in a fluctuating island setting: rapid radiation of Ohomopterus ground beetles in the Japanese Islands

The Japanese Islands have been largely isolated from the East Asian mainland since the Early Pleistocene, allowing the diversification of endemic lineages. Here, we explore speciation rates and historical biogeography of the ground beetles of the subgenus Ohomopterus (genus Carabus) based on nuclear and mitochondrial gene sequences. Ohomopterus diverged into 15 species during the Pleistocene. The speciation rate was 1.92 Ma(-1) and was particularly fast (2.37 Ma(-1)) in a group with highly divergent genitalia. Speciation occurred almost solely within Honshu, the largest island with complex geography. Species diversity is highest in central Honshu, where closely related species occur parapatrically and different-sized species co-occur. Range expansion of some species in the past has resulted in such species assemblages. Introgressive hybridization, at least for mitochondrial DNA, has occurred repeatedly between species in contact, but has not greatly disturbed species distinctness. Small-island populations of some species were separated from main-island populations only after the last glacial (or the last interglacial) period, indicating that island isolation had little role in speciation. Thus, the speciation and formation of the Ohomopterus assemblage occurred despite frequent opportunities for secondary contact and hybridization and the lack of persistent isolation. This radiation was achieved without substantial ecological differentiation, but with marked differentiation in mechanical agents of reproductive isolation (body size and genital morphology).

Mechanical barriers to introgressive hybridization revealed by mitochondrial introgression patterns in Ohomopterus ground beetle assemblages

To reveal the role of diverged body size and genital morphology in reproductive isolation among closely related species, we examined patterns of, and factors limiting, introgressive hybridization between sympatric Ohomopterus ground beetles in central Japan using mitochondrial NADH dehydrogenase subunit 5 (ND5) gene sequences. We sampled 17 local assemblages that consisted of two to five species and estimated levels of interspecific gene flow using the genetic distance, D(A), and maximum-likelihood estimates of gene flow. Sharing of haplotypes or haplotype lineages was detected between six of seven species that occurred in the study areas, indicating mitochondrial introgression. The intensity and direction of mitochondrial gene flow were variable among species pairs. To determine the factors affecting introgression patterns, we tested the relationships between interspecific D(A) and five independent variables: difference in body size, difference in genital size, phylogenetic relatedness (nuclear gene sequence divergence), habitat difference, and species richness of the assemblage. Body and genital size differences contributed significantly to preventing gene flow. Thus, mechanical isolation mechanisms reduce the chance of introgressive hybridization between closely related species. Our results highlight the role of morphological divergence in speciation and assemblage formation processes through mechanical isolation.

Utility of Nuclear Allele Networks for the Analysis of Closely Related Species in the Genus Carabus, Subgenus Ohomopterus

Nuclear DNA sequence data for diploid organisms are potentially a rich source of phylogenetic information for disentangling the evolutionary relationships of closely related organisms, but present special phylogenetic problems owing to difficulties arising from heterozygosity and recombination. We analyzed allelic relationships for two nuclear gene regions (phosphoenolpyruvate carboxykinase and elongation factor-1a), along with a mitochondrial gene region (NADH dehydrogenase subunit 5), for an assemblage of closely related species of carabid beetles (Carabus subgenus Ohomopterus). We used a network approach to examine whether the nuclear gene sequences provide substantial phylogenetic information on species relationships and evolutionary history. The mitochondrial gene genealogy strongly contradicted the morphological species boundary as a result of introgression of heterospecific mitochondria. Two nuclear gene regions showed high allelic diversity within species, and this diversity was partially attributable to recombination between various alleles and high variability in the intron region. Shared nuclear alleles among species were rare and were considered to represent shared ancestral polymorphism. Despite the presence of recombination, nuclear allelic networks recovered species monophyly more often and presented genetic differentiation patterns (low to high) among species more clearly. Overall, nuclear gene networks provide clear evidence for separate biological species and information on the phylogenetic relationships among closely related carabid beetles.

Reconstructing species phylogeny of the carabid beetles Ohomopterus using multiple nuclear DNA sequences: heterogeneous information content and the performance of simultaneous analyses

We attempted a phylogenetic reconstruction for the carabid subgenus Ohomopterus (genus Carabus), a notable case of radiation with mitochondrial introgression across species. Sequence data from five nuclear single copy loci were used, including wingless (Wg), phosphoenolpyruvate carboxykinase (PepCK), cytochrome c (Cytc), elongation factor-1alpha (EF-1alpha), and an anonymous single copy locus (Carab1). Sequences of Cytc, EF-1alpha, and Carab1 included intron or intron-like parts with length variation. The analysis of individual loci resulted in low resolution of the phylogenetic relationships, and the monophyly of several morphologically recognized species for which multiple specimens were analyzed was not revealed. Several specimens were heterozygous, with non-monophyletic alleles observed in three of the five loci at which alleles in heterozygotes were separated. In a simultaneous analysis of the five loci with ambiguously aligned parts eliminated and heterozygotic sites treated as missing, the resulting tree was well resolved, but the branch support was generally weak because of conflicting phylogenetic signals from different loci. We also attempted to incorporate allelic sequence data plus the ambiguously aligned parts in the analysis, by using all possible combinations of alleles from different loci in heterozygotic individuals, but the resultant tree was not supported more strongly. Nonetheless, these simultaneous analyses provided support for the monophyly of several species and species groups, and revealed the basic evolutionary trend of OHOMOPTERUS: initial widespread groups with simpler genitalia and the origination of exaggerated genitalia in a derived clade. This study exemplifies problems inherent in the phylogenetic reconstruction of closely related organisms where low levels of variation limit the information content from each locus, while heterozygosity, different phylogenetic history of multiple loci, and alignment ambiguity further hamper phylogenetic reconstruction unless several loci converge on a uniform signal.

Sympatric speciation by sexual conflict

It is well established that sexual conflict can drive an endless coevolutionary chase between the sexes potentially leading to genetic divergence of isolated populations and allopatric speciation. We present a simple mathematical model that shows that sexual conflict over mating rate can result in two other general regimes. First, rather than "running away" from males, females can diversify genetically into separate groups, effectively "trapping" the males in the middle at a state characterized by reduced mating success. Female diversification brings coevolutionary chase to the end. Second, under certain conditions, males respond to female diversification by diversifying themselves. This response results in the formation of reproductively isolated clusters of genotypes that emerge sympatrically.

Extensive trans-species mitochondrial polymorphisms in the carabid beetles Carabus subgenus Ohomopterus caused by repeated introgressive hybridization

To study the potential importance of introgressive hybridization to the evolutionary diversification of a carabid beetle lineage, we studied intraspecific and trans-species polymorphisms in the mitochondrial NADH dehydrogenase subunit 5 (ND5) gene sequence (1083 bp) in four species of the subgenus Ohomopterus (genus Carabus) in central and eastern Honshu, Japan. Of the four species, C. insulicola is parapatric with the other three, and can hybridize naturally with at least two. This species possesses two haplotypes of remote lineages. We classified ND5 haplotypes using polymerase chain reaction-restriction fragment length polymorphism with TaqI endonuclease for 524 specimens, and sequenced 143 samples. Analysis revealed that each species was polyphyletic in its mitochondrial DNA phylogeny, representing a marked case of trans-species polymorphism. Recent one-way introgression of mitochondria from C. arrowianus nakamurai to C. insulicola, and from C. insulicola to C. esakii, was inferred from the frequency of identical sequences between these species and from direct evidence of hybridization in their contact zones. Other intraspecific polymorphisms in the four species may be due to undetected introgressive hybridization (e.g. C. insulicola to C. maiyasanus) or from stochastic lineage sorting of ancestral polymorphisms. This beetle group has a genital lock-and-key system, with species-specific or subspecies-specific genital morphology that may act as a barrier to hybridization. However, our results demonstrate that introgressive hybridization has occurred multiple times, at least for mitochondria, despite differences among, and stability within, morphological characters that distinguish local populations. Thus, hybridization and introgression could have been key processes in the evolutionary diversification of Ohomopterus.