Parallel evolution in radiation of Ohomopterus ground beetles inferred from mitochondrial ND5 gene sequences

Molecular phylogeny and diversification timing of the Nemouridae family (Insecta, Plecoptera) in the Japanese Archipelago

The generation of the high species diversity of insects in Japan was profoundly influenced by the formation of the Japanese Archipelago. We explored the species diversification and biogeographical history of the Nemouridae Billberg, 1820 family in the Japanese Archipelago using mitochondrial DNA and nuclear DNA markers. We collected 49 species among four genera: Indonemoura Baumann, 1975; Protonemura Kempny, 1898; Amphinemura, Ris 1902 and Nemoura Latreille, 1796 in Japan, China, South Korea and North America. We estimated their divergence times-based on three molecular clock node calibrations-using Bayesian phylogeography approaches. Our results suggested that Japanese Archipelago formation events resulted in diversification events in the middle of the Cretaceous (<120 Ma), speciation in the Paleogene (<50 Ma) and intra-species diversification segregated into eastern and western Japan of the Fossa Magna region at late Neogene (20 Ma). The Indonemoura samples were genetically separated into two clades-that of Mainland China and that of Japan. The Japanese clade clustered with the Nemouridae species from North America, suggesting the possibility of a colonisation event prior to the formation of the Japanese Archipelago. We believe that our results enhanced the understanding both of the origin of the species and of local species distribution in the Japanese Archipelago.

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).

Phylogeography and introgressive hybridization of the ground beetle Carabus yamato in Japan based on mitochondrial gene sequences

To study the phylogeography of the ground beetle Carabus yamato in Japan, we compared 1,020-bp sequences of the mitochondrial NADH dehydrogenase subunit 5 (ND5) gene from 373 specimens from 37 localities with those of three parapatric species (C. albrechti, C. kimurai, and C. japonicus) that might share mitochondrial lineages with C. yamato through introgressive hybridization. We found 81 haplotypes from C. yamato. Of these, 17 haplotypes were considered to be of an introgressed lineage from C. albrechti, based on the phylogeny and geographic distribution. In addition, one haplotype of C. kimurai was likely an introgressant from C. yamato. Putative introgression events among the four species were restricted to these two directional cases. We analyzed the phylogeography of C. yamato using nested clade phylogeographical analysis and population genetic parameters. The mitochondrial lineages of C. yamato were estimated to have diverged no more than approximately 1.12 million years ago, implying that the estimated historical events occurred after the Early Pleistocene. Carabus yamato was inferred to have experienced a contraction of its distribution range, followed by recent range expansion. Populations in the western and eastern regions, segregated by Ise Bay and the Nobi and Okazaki Plains, diverged in the mitochondrial clades. The northern and most western populations possessed one clade only (except an introgressed lineage), whereas eastern and some southwestern populations possessed several diverged clades, which were considered to be ancestral; these populations may have been associated with refugia during glacial periods.

Rapid diversification of male genitalia and mating strategies in Ohomopterus ground beetles

We analysed evolutionary diversification and covariation in male genitalia and four mating traits related to sexual selection, i.e. testis size, spermatophore size, copulation duration and post-copulatory guarding duration, in Ohomopterus ground beetles using phylogenetically independent contrasts. Male genital size and mating duration have evolved more rapidly than body size and the other traits studied. Male genital size was negatively correlated with copulation duration, suggesting that elongated male genitalia may enable decreased time investment in a single copulation because it is more effective at facilitating spermatophore deposition. Male genital size was positively correlated with spermatophore size, suggesting coevolution between offensive and defensive male mating tactics because the elongated male genitalia may be advantageous in displacement of rivals' plug-like spermatophores, and decreased mating duration may intensify sperm competition. Thus, the remarkable diversity of male genitalia in Ohomopterus may have been facilitated by the interplay between inter- and intrasexual selection processes.

Establishment of hybrid-derived offspring populations in the Ohomopterus ground beetles through unidirectional hybridization

An approach to deduce the mechanism of stabilization of the hybrid-derived populations in the Ohomopterus ground beetles has been made by comparative studies on the phylogenetic trees of the mitochondrial and nuclear DNA. A phylogenetic tree based on the internal transcribed spacer (ITS) of nuclear ribosomal gene roughly reflects the relations of morphological species group, while mitochondrial NADH dehydrogenase subunit 5 (ND5) gene shows a considerable different topology on the tree; there exist several geographically-linked lineages, most of which consist of more than one species. These results suggest that the replacement of mitochondria has occurred widely in the Ohomopterus species. In most cases, hybridization is unidirectional, i.e., the species A (♂) hybridized with another species B (♀) and not vice versa, with accompanied replacement of mitochondria of A by those of B. The results also suggest that partial or complete occupation of the distribution territory by a hybrid-derived morphological species. The morphological appearance of the resultant hybrid-derivatives are recognized as that of the original species A. Emergence of a morphological new species from a hybrid-derived population has been exemplified.

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.

Phylogeny and evolution of Digitulati ground beetles (Coleoptera, Carabidae) inferred from mitochondrial ND5 gene sequences

Genealogical trees have been constructed using mitochondrial ND5 gene sequences of 87 specimens consisting of 32 species which have been believed to belong to the division Digitulati (one of the lineages of the subtribe Carabina) of the world. There have been recognized six lineages, which are well separated from each other. Each lineage contains the following genus: (1) the lineage A: Ohomopterus from Japan; (2) the lineage B: Isiocarabus from eastern Eurasian Continent; (3) the lineage C: Carabus from China which are further subdivided into three sublineages; (4) the lineage D: Carabus from USA; (5) the lineage E: Carabus from the Eurasian Continent, Japan and North America; and (6) the lineage F: Eucarabus from the Eurasian Continent. Additionally, the genus Acrocarabus which had been treated as a constituent of the division Archicarabomorphi has been recognized to be the 7th lineage of the division Digitulati from the ND5 genealogical analysis as well as morphology. These lineages are assumed to have radiated within a short period and are largely linked to their geographic distribution.

Pattern of phylogenetic diversification of the Cychrini ground beetles in the world as deduced mainly from sequence comparisons of the mitochondrial genes

The phylogenetic position of the tribe Cychrini within the subfamily Carabinae (the family Carabidae) was estimated by comparing the nucleotide sequences of the mitochondrial NADH dehydrogenase subunit 5 (ND5) gene and the nuclear 28S ribosomal DNA (rDNA). The phylogenetic trees suggest that the Cychrini would most probably be the oldest line within the Carabinae. Phylogenetic trees were constructed by comparing the mitochondrial cytochrome C oxidase subunit I (COI) gene sequences from 33 species of the Cychrini from various localities that include the whole distribution ranges of the representative species within all the known genera in the world. The trees suggest that the Cychrini members radiated into a number of phylogenetic lineages within a short period, starting about 44 million years ago (MYA). Most of the phylogenetic lineages or sublineages are geographically linked, each consisting of a single or only a few species without scarce morphological differentiation in spite of their long evolutionary histories (silent or near-silent evolution [see Adv. Biophys. 36 (1999) 65; J. Mol. Evol. 53 (2001) 517]). The fact suggests that the geographic isolation per se did not bring about conspicuous morphological differentiation. The phylogenetic lineages of the Cychrini well correspond to the taxonomically defined genera and the subgenera.

Mode of morphological differentiation in the Latitarsi-ground beetles (Coleoptera, Carabidae) of the world inferred from a phylogenetic tree of mitochondrial ND5 gene sequences

The Latitarsi is one large division of the subtribe Carabina (subfamily Carabinae, family Carabidae), and has been considered as a discrete morphological group consisting of 17 genera. The phylogenetic relationships and evolutionary pattern of the Latitarsi ground beetles have been investigated by analyzing mitochondrial NADH dehydrogenase subunit 5 (ND5) gene sequences. The phylogenetic tree suggests that the Latitarsi members do not form a single cluster, i.e., not monophyletic and at least 16 lineages belonging to the so-called Latitarsi emerged at about the same time of the Carabina radiation together with the members of other divisions. This suggests that these lineages (A, B, C, H, L, N, O, P, Q, R, S, T, U, V, W and X in Fig. 2a) may be treated each as a phylogenetically distinct division equivalent to other divisions. The group with bootstrap value of more than 80 percent has been considered as a single lineage (division) with two exceptions, V and X. The independency of each lineage has been assumed by the traditional morphology as well as a single clustering on the trees constructed by independent methods, unchanged topology by replacement of outgroups, etc. Generally speaking, the members in a single lineage are geographically linked. Many phylogenetic lineages are composed of a single or only a few species without conspicuous morphological differentiation. In contrast to such a "silent morphological evolution", a remarkable morphological differentiation occasionally took place in several lineages.

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.

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.

Rapid evolution in the Nebria gregaria group (Coleoptera: Carabidae) and the paleogeography of the Queen Charlotte Islands

Morphological differentiation in the ground beetles of the Nebria gregaria group, found on the Queen Charlotte Islands, has been used as support for the glacial refugium proposed for the northwest coast of North America. Two members of this species group, N. charlottae and N. louiseae, are restricted to cobble beaches in this archipelago. A third, N. haida, is found only in alpine regions of the archipelago and the adjacent mainland. The remaining two species of the gregaria group, N. lituyae and N. gregaria, show highly restricted distributions in the mountains of the Alaska panhandle and on the beaches of the Aleutian Islands, respectively. To determine the relationships of the five species, we conducted phylogenetic analyses on nucleotide sequence data obtained from five regions of the mitochondrial DNA. In total, 1835 bp were analyzed. The results suggest that one species, N. lituyae, does not belong in the gregaria group, and that only seven mutations separated the two most divergent of the four remaining species. We also conducted random amplified polymorphic DNA fingerprinting analyses on genomic DNA extracted from the five species. Analyses of genetic diversity revealed a lack of molecular differentiation among the Queen Charlotte species, suggesting that these populations may be postglacial in origin and that together N. gregaria, N. charlottae, N. louiseae, and N. haida might represent local variations of a single species. These results are consistent with conclusions derived for the morphological and genetical differentiation among Gasterosteus populations in the archipelago.

Differentiation within the genus Leptocarabus (excl. L. kurilensis) in the Japanese Islands as deduced from mitochondrial ND5 gene sequences (Coleoptera, Carabidae)

The phylogenetic trees have been constructed for the mitochondrial ND5 gene sequences from the Japanese Leptocarabus ground beetles, which contain 101 specimens collected from nearly the complete distribution ranges of them consisting of five morphological species, i.e., Leptocarabus procerulus, L. kumagaii, L. hiurai, L. kyushuensis and L. arboreus. On the trees, there are recognized two major lineages, each of which is further divided into two or more sublineages. The phylogenetic lines are geographically linked. Two or more species occur in a single lineage, and the same species appear in different lines. We suggest that transformation from one type of morphology to another took place in parallel in various periods of evolution of the Japanese Leptocarabus. From the phylogenetic tree and the dating from the nucleotide substitution rate and the geohistorical data it is inferred that the ancestry of all the Japanese Leptocarabus species was derived from a protoform of L. kyushuensis inhabited the ancient Japan area, followed by separation into two lineages after split of the Japanese Islands from the Eurasian Continent. They then propagated distribution to occupy their own habitat ranges, during which the morphological transformation took place in some lineages.

The current state of insect molecular systematics: a thriving Tower of Babel

Insect molecular systematics has undergone remarkable recent growth. Advances in methods of data generation and analysis have led to the accumulation of large amounts of DNA sequence data from most major insect groups. In addition to reviewing theoretical and methodological advances, we have compiled information on the taxa and regions sequenced from all available phylogenetic studies of insects. It is evident that investigators have not usually coordinated their efforts. The genes and regions that have been sequenced differ substantially among studies and the whole of our efforts is thus little greater than the sum of its parts. The cytochrome oxidase I, 16S, 18S, and elongation factor-1 alpha genes have been widely used and are informative across a broad range of divergences in insects. We advocate their use as standards for insect phylogenetics. Insect molecular systematics has complemented and enhanced the value of morphological and ecological data, making substantial contributions to evolutionary biology in the process. A more coordinated approach focused on gathering homologous sequence data will greatly facilitate such efforts.

Pattern of morphological diversification in the Leptocarabus ground beetles (Coleoptera: Carabidae) as deduced from mitochondrial ND5 gene and nuclear 28S rDNA sequences

Most of the mitochondrial NADH dehydrogenase subunit 5 (ND5) gene and a part of nuclear 28S ribosomal RNA gene were sequenced for 14 species of ground beetles belonging to the genus Leptocarabus. In both the ND5 and the 28S rDNA phylogenetic trees of Leptocarabus, three major lineages were recognized: (1) L. marcilhaci/L. yokoael/Leptocarabus sp. from China, (2) L. koreanus/L. truncaticollis/L. seishinensis/L. semiopacus/L. canaliculatus/L. kurilensis from the northern Eurasian continent including Korea and Hokkaido, Japan, and (3) all of the Japanese species except L. kurilensis. Clustering of the species in the trees is largely linked to their geographic distribution and does not correlate with morphological characters. The species belonging to different species groups are clustered in the same lineages, and those in the same species group are scattered among the different lineages. One of the possible interpretations of the present results would be that morphological transformations independently took place in the different lineages, sometimes with accompanying parallel morphological evolution, resulting in the occurrence of the morphological species belonging to the same species group (= type) in the different lineages.

Molecular evidence for parallel evolution of adaptive syndromes in fig-breeding Lissocephala (Drosophilidae)

Afrotropical Lissocephala (Drosophilidae) breed strictly in syconia (figs) of Ficus (Moraceae) and have accordingly evolved specific features including modified female and eggshell morphologies, ovipositing, larval foraging, and mating behaviors. These various traits may exist as two or three alternative states. Each species displays a specific suite of traits so closely coordinated with one another that alternative states of the overall suites of traits can be seen as "adaptive syndromes." Three clear-cut adaptive syndromes can be recognized while two taxonomic lineages (juncta and sanu species groups) are traditionally accepted on the basis of male terminalia. A crucial evolutionary question results from the consideration that ecological clusters and taxonomic groups have conflicting compositions: the three syndromes are found in the juncta group while two of them occur in the sanu group. To resolve this conflict, we present molecular data which provide a robust phylogeny: mitochondrial DNA (12S + 16S ribosomal DNA and cytochrome b) sequence data are in agreement with one another regardless of the algorithm used. All molecular data consistently support male terminalia dichotomy. Such a level of consistency unambiguously indicates that parallel evolution of adaptive syndromes occurred. Thus, homoplasy may affect morphological and behavioral traits concomitantly when these are involved in a network of functional relationships.