Pattern of Morphological Diversification in the Leptocarabus Ground Beetles (Coleoptera: Carabidae) as Deduced from Mitochondrial ND5 Gene and Nuclear 28S rDNA Sequences

Choong-Gon Kim*, Hong-Zhang Zhou{dagger}, Yûki Imura{ddagger}, Osamu Tominaga§, Zhi-Hui Su2,* and Syozo Osawa*

*JT Biohistory Research Hall, Osaka, Japan; and
{dagger}Institute of Zoology, Chinese Academy of Sciences, Beijing, China;
{ddagger}Department of Gynecology, Tôkyû General Hospital, Tokyo, Japan; and
§Shibatsuji-Cho, Nara, Japan


    Abstract
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results and Discussion
 Acknowledgements
 References
 
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. yokoae/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.


    Introduction
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results and Discussion
 Acknowledgements
 References
 
Leptocarabus is a well-known higher taxon in the subtribe Carabina of the family Carabidae. Although often ranked as one of the subgenera of the grand genus Carabus, we treat it as a distinct genus for convenience. The genus contains 20 species (Imura and Mizusawa 1996Citation ; Kwon and Lee 1984Citation ). All of the Leptocarabus species are hindwingless and can move only by walking. They are distributed throughout the northern and eastern Eurasian continent including the Korean Peninsula, the Japanese Islands, Sakhalin, and the Kurils. Most of the species are large and elongate in shape and black to brown in color, while some others are characterized by a smaller size and metallic tint to the body (fig. 1 ). The Leptocarabus species have been well studied taxonomically and are classified into several subgenera (see below) (Nakane 1961, 1962Citation ; Ishikawa 1972, 1992Citation ; Kwon and Lee 1984Citation ; Imura and Mizusawa 1996Citation ). According to Imura and Mizusawa (1996)Citation , they are classified into five species groups (hereinafter called "types" for convenience), Procerulus (P) (=subgenus Leptocarabus [s. str.]), Semiopacus (S) (=subgenus Adelocarabus), Koreanus (K) (=subgenus Weolseocarabus), Canaliculatus (C) (=subgenus Aulonocarabus), and Truncaticollis (T) (no proposed subgenus), based on the morphological characters (table 1 ). The P type contains Leptocarabus marcilhaci, Leptocarabus yokoae, and Leptocarabus sp. from central China and Leptocarabus procerulus, Leptocarabus kumagaii, Leptocarabus hiurai, and Leptocarabus kyushuensis from Japan. This type is said to be closely related to the K type based on morphological characters (Ishikawa 1972, 1989Citation ). The S type contains Leptocarabus semiopacus, Leptocarabus seishinensis, a few other species mainly from Korea, and Leptocarabus arboreus from Japan and Sakhalin (Kwon and Lee 1984Citation ; Ishikawa 1992;Citation Imura and Mizusawa 1996Citation ). The K type is composed of a single species, Leptocarabus koreanus, which inhabits the Korean Peninsula (Kwon and Lee 1984Citation ; Imura and Kezuka 1992Citation ). The C type species are widely distributed in the eastern Eurasian continent, Sakhalin, Hokkaido, and the Kurils. The T type contains four species which are widely distributed from the northeastern part of the Eurasian continent to Alaska (Imura and Mizusawa 1996Citation ).



View larger version (43K):
[in this window]
[in a new window]
 
Fig. 1.—Leptocarabus ground beetles. 1—L. marcilhaci; 2—L. yokoae; 3—L. kyushuensis kyushuensis; 4—L. kyushuensis cerberus; 5—L. hiurai; 6—L. procerulus procerulus; 7—L. procerulus miyakei; 8—L. kumagaii; 9—L. semiopacus; 10—L. seishinensis seishinensis; 11—L. arboreus arboreus; 12—L. arboreus gracillimus; 13—L. koreanus coreanicus; 14—L. kurilensis rausuanus; 15—L. canaliculatus jankowskiellus; 16—L. truncaticollis. Scale bar: 10 mm.

 

View this table:
[in this window]
[in a new window]
 
table 1 Morphological Characters of the Five Species Groups (Types)

 
The taxonomic and phylogenetic studies for the genus Leptocarabus based on morphology need to be reexamined with more explicit procedures, such as molecular evolution studies, because discrepancies have often been reported between morphology and molecular studies for some other carabid group (Su et al. 1996b, 1998Citation ). In the present paper, phylogenetic relationships have been estimated by analyzing a large part of the ND5 gene sequences and a reasonable length of the nuclear 28S rDNA of 52 examples covering the representative Leptocarabus species and geographic races from their known localities. The molecular phylogeny uncovered an interesting evolutionary history of Leptocarabus that cannot be reached by cladistic analysis that depends on morphology alone.


    Materials and Methods
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results and Discussion
 Acknowledgements
 References
 
Isolation of DNA
Names and localities of the specimens used are shown in table 2 and figure 2 . To prevent degradation of DNA, the specimens were immediately killed in 95% ethanol and stored in ethanol until use. A single adult individual of each species was used for DNA extraction. Total DNA was extracted from thorax muscle (10–25 mg) using the QIAamp tissue kit (QIAGEN GmbH, Germany).


View this table:
[in this window]
[in a new window]
 
table 2 List of the Leptocarabus Specimens Used in this Study

 


View larger version (34K):
[in this window]
[in a new window]
 
Fig. 2.—Specimen locations of the Leptocarabus species used in the present study. (ô) S type species; () P type species; ({diamondsuit}) K type species; ({blacksquare}) C type species; (*) T type species; ({blacktriangleup}) outgroup species. Numerals correspond to the locality numbers in table 2 and figures 3–5 . Approximate distribution boundaries for respective lineages are shown by dotted lines.

 
PCR Amplification and DNA Sequencing
Altogether, 52 specimens from the Eurasian continent and the Japanese Islands were analyzed. Total DNA was used as a template for amplification of ND5 DNA and 28S rDNA fragments by the polymerase chain reaction (PCR) (Saiki et al. 1988Citation ). The 1,084-bp sequence which contains a 1,069-bp 3' region of the ND5 gene, 8 bp of noncoding sequence, and a 7-bp 5' terminus of the Phe-tRNA gene was amplified by the ver 1.04 primer (5'-GTC ATA CTC TAA ATA TAA GCT A-3') and the ver 1.06 primer (5'-CCT GTT TCT GCT TTA GTT CA-3') (Su et al. 1996aCitation ). The 852–881-bp sequence of the 5' region of 28S rDNA was amplified by the 28S-01 primer (5'-GAC TAC CCC CTG AAT TTA AGC AT-3') and the 28SR-01 primer (5'-GAC TCC TTG GTC CGT GTT TCA AG-3'), which were designed from the 28S rDNA sequences of Drosophila melanogaster (Tautz et al. 1988Citation ) and Xenopus laevis (Clark et al. 1984Citation ). Mostly, the two primers used for PCR were sufficient to read 1,084 bp of the ND5 sequences and 852–881 bp of the 28S rDNA sequences. In some cases, four internal primers were used for the ND5 gene: forward primers Ezo-2 (5'-TTC ATC TTT TAA CTC ATG CA-3') and RC4-4 (5'-GAT CAA GGT TGA AAT GAA T-3'), and reverse primers AO-3 (5'-ATA TTC ATT TCA ACC TTG ATC-3') and RCE-2 (5'-TGC ATG AGT TAA AAG ATG AA-3') (Su et al. 1998Citation ). PCR amplifications were carried out in a 100-µl mixture containing 10 mM Tris-HCl (pH 8.3), 50 mM KCl, 1.5 mM MgCl2, 0.01% gelatin, 0.2 mM of each dNTP, 100 pM of each primer, and 2.5 U of Taq polymerase (TaKaRa EX Taq, Takara). Each reaction mixture was overlaid by one drop of mineral oil. PCR was performed for 35 cycles of denaturation at 94°C for 1 min, primer annealing at 50°C for 1 min, and extension at 70°C for 2 min. The final single cycle was performed under the same conditions but with an extension step at 70°C for 7 min using a DNA Thermal Cycler 480 (Perkin Elmer). The PCR product was purified with QIAquick PCR purification kit (QIAGEN GmbH).

Direct sequencing was performed with an automated ABI PRISM 377 DNA sequencer using the dideoxy chain termination method (Sanger, Nicklen, and Coulson 1977Citation ). The reaction mixture for cycle sequencing consisted of 6 µl of dRhodamine terminator cycle sequencing ready reaction with AmpliTaq DNA Polymerase, FS (Applied Biosystems, Foster City, Calif.), 0.1–0.3 pmol/µl of DNA, 2.4 µl (1 pmol/µl) of sequencing primer, and distilled water to a total volume of 15 µl. The cycle-sequencing conditions were 25 cycles at 96°C for 10 s, 50°C for 5 s, and 60°C for 4 min, followed by an indefinite hold at 4°C using a GeneAmp PCR system 9600 (Perkin Elmer). The DNA product was cleaned with Centri-Sep spin columns (Applied Biosystems) and vacuum-dried before applying. The nucleotide sequence data reported in this paper will appear in the DDBJ, EMBL, and GenBank nucleotide sequence databases with the accession numbers shown in table 2 .

Phylogenetic Analysis
The ND5 and the 28S rDNA sequences were aligned and compared using the multiple-alignment program CLUSTAL W (Thompson, Higgins, and Gibson 1994Citation ) and DNASIS, version 3.7 (Hitachi Software Engineering, Japan). The evolutionary distances (D) were computed by Kimura’s (1980)Citation two-parameter method, and the phylogenetic trees were constructed by the neighbor-joining (NJ) method (Saitou and Nei 1987Citation ) and the unweighted pair grouping method with arithmetic means (UPGMA). All of these processes were performed with the DNA sequence analysis package SINCA, version 3.0 (Fujitsu System Engineering, Japan). A maximum-parsimony (MP) tree was also constructed with PHYLIP, version 3.5 (Felsenstein 1993Citation ). Bootstrap analysis was performed for all the trees (Felsenstein 1985Citation ) based on 500 resamplings. The gene sequences of two Calosomina species were used as an outgroup.

Dating
For setting the timescale, a 0.01D unit corresponding to 3.6 Myr was used (Su et al. 1998Citation ; Osawa et al. 1999;Citation revised by Su et al. [unpublished data]).


    Results and Discussion
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results and Discussion
 Acknowledgements
 References
 
Sequence Divergence and G+C Content of DNA
The maximum sequence divergences of the ND5 gene were 3.81% for the Japanese species and 11.24% for the species on the Eurasian continent, indicating that the Japanese species are much more closely related to each other than are the species on the continent. The maximum sequence divergences of the 28S rDNA were 0.81% for the Japanese species and 2.76% for the continental species (table 3 ). The 28S rDNA evolved about one fourth as fast as the ND5 gene, so a meaningful resolution between allied species was not possible on the 28S rDNA phylogenetic tree. There were no insertions/deletions and no length variations for the ND5 gene throughout the species examined. A few deletions were found in some 28S rDNA sequences. The G+C content was nearly constant (21 ± 0.5%) for the ND5 gene, and it was 58 ± 0.8% for 28S rDNA in the Leptocarabus species that were analyzed.


View this table:
[in this window]
[in a new window]
 
table 3 Pairwise Sequence Divergence (by Kimura’s [1980] Two-Parameter Method) for ND5 DNA and 28S rDNA of the Representative Leptocarabus Species that Were Examined in this Study

 
Relation Between Phylogeny and Morphology
figures 3 and 4 show the NJ phylogenetic trees of the ND5 gene and the nuclear 28S rDNA, respectively. The UPGMA and MP trees gave essentially the same topology (not shown). Three major lineages are recognized in both the ND5 and 28S rDNA trees. Lineage 1 is composed of three Chinese species, L. marcilhaci, L. yokoae, and Leptocarabus sp., all belonging to the Procerulus (P) type (for classification of the type = species group, see Introduction). Lineage 2 includes six species from the northeastern Eurasian continent including the Korean Peninsula and Hokkaido, Japan: L. koreanus (Koreanus [K] type), L. truncaticollis (Truncaticollis [T] type), L. seishinensis and L. semiopacus (both Semiopacus [S] type), and L. canaliculatus and L. kurilensis (both Canaliculatus [C] type). All the Japanese species belong to lineage 3, which contains L. kyushuensis, L. hiurai, L. procerulus, and L. kumagaii (all P type) and L. arboreus (S type). Of special interest is that L. truncaticollis is very different from other Leptocarabus species in appearance (see fig. 1 ), and yet it belongs to lineage 2 together with the S and C type species. Thus, clustering of the species examined is closely related to their geographic distribution and does not correlate with morphologically defined type. The members of the P type appear in both lineage 1 and lineage 3; those of the S type appear in lineages 2 and 3, which are, respectively, related to the species belonging to the other types in the same lineage on the phylogenetic trees.



View larger version (40K):
[in this window]
[in a new window]
 
Fig. 3.—Neighbor-joining phylogenetic tree of the ND5 gene of the genus Leptocarabus. Throughout figures 3–5 , distance (D) denotes Kimura’s (1980) two-parameter evolutionary distance. The value at the node represents the percentage bootstrap confidence level based on 500 resamplings. The tree was outgroup-rooted using the ND5 gene sequences of two species of the Calosomina, Calosoma inquisitor and Campalita chinense. The symbol of type (=species group) by morphology (Imura and Mizusawa 1996Citation ) is indicated on the right. For explanation of the types, see the text.

 
Interpretations of the Results
There may be several ways to interpret the phylogenetic trees obtained in this study. One is that the P type is the ancestral form of the genus, and the appearance of the P type in lineages 1 and 3 is simply explained by shared ancestral characters. This is possible, but the ND5 tree suggests an almost simultaneous diversification of the three lineages. On the 28S rDNA tree, lineage 1 forms the most basal branch of the genus. Since the base substitution rate of 28S rDNA is much slower than that of the ND5 gene, the molecular clock by the ND5 gene may be more accurate than that by 28S rDNA. The primary purpose of the use of 28S rDNA is simply to show the existence of the three lineages, each containing the same members as the ND5 tree. Thus, we tentatively explain the parallel appearance of the P and S types in lineages 1 and 3.

The possibility of lineage sorting of ancestral polymorphism is now considered. If the ancestral species or population is assumed to contain both P and S (and the other) types, the current character distribution may be explained by random lineage sorting. For the reason mentioned above, there is at present no way to specify the ancestral form. Suppose that the P type is the ancestral form. The distribution range of the P type in lineage 1 is exclusively distributed in central China, where no other types have been discovered in reasonably extensive expeditions by a number of entomologists. The distribution range of lineage 1 does not overlap with that of the other types (see fig. 2 ). There is no evidence showing the presence of polymorphism in a single population. Of course, this does not exclude its absence in the past. The ancestral polymorphism–random lineage sorting possibility may be theoretically pleasing, but it cannot be verified.

One might argue that the ND5 phylogeny would be brought about by introgression of the mitochondrial gene via hybridization. This possibility can be largely ruled out by the overall congruence of the 28S rDNA and the ND5 phylogenies. Furthermore, the distributional isolation of the three lineages, which would have occurred long time ago, is not consistent with the introgression hypothesis.

One plausible interpretation for the appearance of the same type in different lineages would be that parallel evolution took place for the P type in lineages 1 and 3 and for the S type in lineages 2 and 3. If the P type is the ancestral form, parallelism should be considered only for the S type. In addition, the intermingled occurrence of different types in lineages 2 and 3 suggests the morphological transformation from one type to another within the respective lineages. The situation resembles that of Ohomopterus (Su et al. 1996bCitation ), in which taxonomically the "same species" or the members belonging to the same species group (=type) appear in more than two different places on the ND5 tree. Su et al. (1996bCitation ) proposed that parallel evolution took place in different lineages, possibly through discontinuous morphological transformation called "type-switching."

Diversification Within the Respective Lineages
The branching of the three lineages mentioned above seems to have started within a short time, as judged from the ND5 tree. No ancestral lineage can be defined (see above). Within lineage 2, to which the northern and eastern Eurasian species belong, initiation of diversification into various species is much older than that in lineage 3, to which all of the Japanese species belong. Leptocarabus koreanus, L. truncaticollis, L. seishinensis, and L. semiopacus, respectively, form a well-defined cluster. The L. canaliculatus specimens from various localities are clustered together, and L. kurilensis from Hokkaido is included in this cluster, although L. canaliculatus and L. kurilensis are usually regarded as two distinct species by most morphologists. This suggests that L. kurilensis branched off from L. canaliculatus, with accompanying morphological changes. Thus, taxonomy based on morphology at the species level is consistent with the ND5 and 28S rDNA phylogeny in lineage 2, except that L. kurilensis was not separated from L. canaliculatus. In lineage 3 (fig. 5 ), two sublineages are recognized: the L. kyushuensis sublineage (sublineage 1), which is further divided into two clades containing, respectively, the habitants in Kyushu and those in Honshu (the Chugoku district), and sublineage 2, which consists of the habitants in Hokkaido, Honshu, Shikoku, and Kyushu. Within sublineage 2, four clades are recognized: clade a (L. procerulus miyakei in Kyushu), clade b (L. hiurai in Shikoku), clade c (L. arboreus in Hokkaido), and clade d (the habitants in Honshu [L. procerulus procerulus, L. kumagaii, and L. arboreus]). Thus, clustering is more or less linked to geographical distribution. Especially noteworthy are the inhabitants of Honshu (sublineage 2, clade d). As noted above, in clade d there occur three morphological species, L. procerulus (P type), L. kumagaii (P type), and L. arboreus (S type). These three species are clearly separable morphologically, and yet the evolutionary distances among them are either null or very small (less than 0.5% difference), such that they are unresolved on the ND5 phylogenetic tree. This would represent "evolution in action"; i.e., these (sub)species recently started their speciation, in which "type-switching" would be involved. Thus, besides intermingled occurrence of the species of different types in one sublineage, taxonomically the "same species" appear in different sublineages or clades. For example, L. procerulus appears in clade d (subspecies procerulus) and clade a (subspecies miyakei), and L. arboreus appears in clade c (subspecies arboreus and pararboreus) and clade d (many subspecies). The "species" mentioned above are most probably paraphyletic as judged from the ND5 phylogenetic tree; i.e., each "species" arose in parallel in the different pyhlogenetic lines with minor morphological differences and can be recognized as subspecies or local races. Although the 28S rDNA tree is mostly consistent with the ND5 tree, some discrepancies between them do exist, especially in lineage 3. In the 28S rDNA phylogeny, L. hiurai form a clade with L. kyushuensis (Sublineage 1), while in the ND5 phylogeny, L. hiurai clusters with the other members of sublineage 2. In the 28S rDNA tree, L. procerulus procerulus of Honshu forms a clade with L. procerulus miyakei of Kyushu, whereas in the ND5 tree, L. procerulus procerulus of Honshu clusters with the other Honshu members, L. kumagaii and L. arboreus. As was noted above, the maximum D among the Japanese species or races in the 28S rDNA tree is less than 0.005 (0.81%), in contrast to that in the ND5 gene (maximum D = 11.42%). Therefore, it is not possible to decide whether the apparent discrepancies between the ND5 and the 28S rDNA phylogenies are really meaningful or whether or not this situation can be explained by the introgression hypothesis.



View larger version (28K):
[in this window]
[in a new window]
 
Fig. 5.—Neighbor-joining phylogenetic tree of the ND5 gene of the Japanese Leptocarabus.

 
Origin
The origin of Leptocarabus may be traced back on a phylogenetic tree of the genus Carabus (s. lat.) covering the representative species from all over the world (Imura et al. 1998Citation ). According to Imura et al. (1998)Citation , Leptocarabus started to radiate slightly after the "big bang," i.e., a large-scale radiation into the major subgenera or divisions of Carabus (s. lat.). The "big bang" and the radiation of Leptocarabus, which presumably occurred somewhere in the eastern Eurasian continent, were estimated to have taken place 45–40 and 35–30 MYA, respectively. The phylogenetic tree in figure 3 suggests that the ancestral form of Leptocarabus was divided, presumably by geographic barriers, into three lineages, followed by species diversification in each lineage.

The origin of the Japanese Leptocarabus (lineage 3) is now considered. As noted above, the Japanese species started to diversify much later (ca. 12–10 MYA) than the continental species (ca. 28–25 MYA). We tentatively assume that the ancestor of the Japanese Leptocarabus inhabited the ancient Japanese area in the eastern periphery of the continent. Upon its split from the continent ca. 15 MYA, followed by the archipelago formation by an extensive submergence, the Japanese Leptocarabus ancestor was isolated on some island(s). Upon subsequent upheaval of the Japanese Islands (9–6 MYA), the ancestor spread all over Japan, resulting in various species (even different types) and subspecies upon isolation by various barriers such as straits, tectonic lines, rivers, mountains, etc. (For geohistory of the Japanese Islands, see Su et al. [1998Citation ]). It has been generally believed that all the Japanese Leptocarabus immigrated from the Korean Peninsula via land-bridges in the glacial era (<2 MYA) (Ishikawa 1989Citation ). Indeed, morphologies including genitalia of the S type Japanese Leptocarabus are very similar to those of L. seishinensis or L. semiopacus from the Korean Peninsula, and the P type Japanese Leptocarabus species resemble L. koreanus from Korea, or L. yokoae, L. marchilhaci, and Leptocarabus sp. from central China. In particular, L. kyushuensis from Japan is morphologically quite similar to L. yokoae. However, lineage 3 includes solely the Japanese species and none from the Korean Peninsula or China, suggesting that there exist no direct sister relationships between the Japanese and the Korean (and Chinese) species. Moreover, the emergence of the Japanese species is much older (12–10 MYA) than previously thought (<2 MYA). Thus, the present phylogenetic analyses and dating would reject the Korean origin hypothesis, although it is possible, but not very likely, that such a Korean species, if they existed, became extinct or have not been discovered.



View larger version (35K):
[in this window]
[in a new window]
 
Fig. 4.—Neighbor-joining phylogenetic tree of nuclear 28S rDNA of the genus Leptocarabus.

 

    Acknowledgements
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results and Discussion
 Acknowledgements
 References
 
We thank the following colleagues for invaluable specimens used in this study: M. Arai, M. Anzai, Y. Aoki, S. Imasaka, K. Kato, J.-G. Kim, H. Kinugasa, T. Kobayashi, H. Koike, T. Kosaka, M. Kunikane, K. Maekawa, K. Masunaga, S. V. Murzin, H. Matsuzawa, Y. Nagahata, S. Nakamura, S. Nomura, T. S. Okada, M. Okamoto, Y. Okajima, D. V. Obydov, J.-C. Paik, A. Plutenko, K. Seya, M. Shiraishi, R. Shimoyama, T. Takano, S. Uchida, K. Yahiro, and N. Yasuda. We are grateful to Hideko Tanaka for skillful technical assistance. This study was performed as a part of the joint project of JT Biohistory Research Hall, Japan, and the Chinese Academy of Sciences, Beijing, China. H.-Z.Z. is supported by grants KZ952-SI-106 and STZ97-1-03 from the Chinese Academy of Sciences and grant 39893360 from the National Natural Science Foundation of China.


    Footnotes
 
Naruya Saitou, Reviewing Editor

1 Keywords: Leptocarabus ground beetles, mitochondrial ND5 gene, nuclear 28S rDNA, phylogeny, parallel evolution, Back

2 Address for correspondence and reprints: Zhi-Hui Su, JT Biohistory Research Hall, 1-1 Murasaki-Cho, Takatsuki, Osaka 569-1125, Japan. E-mail: su.zhihui{at}ims.brh.co.jp Back


    References
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results and Discussion
 Acknowledgements
 References
 

    Clark, C. G., B. W. Tague, V. C. Ware, and S. A. Gerbi. 1984. Xenopus laevis 28S ribosomal RNA: a secondary structure model and its evolutionary and functional implications. Nucleic Acids Res. 12:6197–6220.[Abstract]

    Felsenstein, J. 1985. Confidence limits on phylogenies: an approach using the bootstrap. Evolution 39:783–791.

    ———. 1993. PHYLIP (phylogeny inference package). Version 3.5c. Distributed by the author, Department of Genetics, University of Washington, Seattle.

    Imura, Y., and H. Kezuka. 1992. Geographical and individual variation of carabid beetles in the species of the subtribe Carabina (3), carabid beetles of the southern part of the Korean Peninsula. Illustrations of selected insects in the world. Series B (Coleoptera), No. 3. Mushi-Sha, Tokyo.

    Imura, Y., C.-G. Kim, Z.-H. Su, and S. Osawa. 1998. An attempt at the higher classification of the Carabina (Coleoptera, Carabidae) based on morphology and molecular phylogeny, with special reference to Apotomopterus, Limnocarabus and Euleptocarabus. Elytra Tokyo 26:17–35.

    Imura, Y., and K. Mizusawa. 1996. The Carabus of the world. Mushi-Sha, Tokyo.

    Ishikawa, R. 1972. Studies on Leptocarabus and its related subgenera of the genus Carabus L. (Coleoptera: Carabidae). Bull. Natl. Sci. Mus. Tokyo 15:19–27.

    ———. 1989. The Japanese Carabina: geographical distribution and speciation within an archipelago (Coleoptera: Carabidae). Nature and Culture No.1. The University Museum, The University of Tokyo.

    ———. 1992. Taxonomic studies on Leptocarabus (Adelocarabus) arboreus (Lewis) (Coleoptera, Carabidae). TMU Bull. Nat. Hist. 1:1–40.

    Kimura, M. 1980. A simple method for estimating evolutionary rate of base substitutions through comparative studies of nucleotide sequences. J. Mol. Evol. 16:111–120.[ISI][Medline]

    Kwon, Y. J., and S.-M. Lee. 1984. Classification of the subfamily Carabinae from Korea (Coleoptera: Carabidae). Insecta Koreana, Series 4, Editorial Committee of Insecta Koreana, Seoul, Korea.

    Nakane, T. 1961. New or little known Coleoptera from Japan and its adjacent regions. XV. Fragm. Coleopt. 1:1–6.

    ———. 1962. Insecta Japonica, Series 2, Part 3. Coleoptera: Carabidae (1). Hokuryukan, Tokyo.

    Osawa, S., Z.-H. Su, C.-G. Kim, M. Okamoto, O. Tominaga, and Y. Imura. 1999. Evolution of the carabid ground beetles. Adv. Biophys. 36:65–106.[ISI][Medline]

    Saiki, R. K., D. H. Gelfand, S. Stoffl, S. T. Scharf, R. Higuchi, G. T. Horn, K. B. Mullis, and H. A. Erlich. 1988. Primer-directed enzymatic amplification of DNA with a thermostable DNA polymerase. Science 239:487–491.

    Saitou, N., and M. Nei. 1987. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol. Biol. Evol. 4:406–425.[Abstract]

    Sanger, F., S. Nicklen, and A. R. Coulson. 1977. DNA sequencing with chain-terminating inhibitors. Proc. Natl. Acad. Sci. USA 74:5463–5468.

    Su, Z.-H., T. Ohama, T. S. Okada, K. Nakamura, R. Ishikawa, and S. Osawa. 1996a. Phylogenetic relationships and evolution of the Japanese Carabinae ground beetles based on mitochondrial ND5 gene sequences. J. Mol. Evol. 42:124–129.

    Su, Z.-H., O. Tominaga, T. Ohama, E. Kajiwara, R. Ishikawa, T. S. Okada, K. Nakamura, and S. Osawa. 1996b. Parallel evolution in radiation of Ohomopterus ground beetles inferred from mitochondrial ND5 gene sequences. J. Mol. Evol. 43:662–671.

    Su, Z.-H., O. Tominaga, M. Okamoto, and S. Osawa. 1998. Origin and diversification of hindwingless Damaster ground beetles within the Japanese Islands as deduced from mitochondrial ND5 gene sequences (Coleoptera, Carabidae). Mol. Biol. Evol. 15:1026–1039.[Abstract]

    Tautz, D., J. M. Hancock, D. A. Webb, C. Tautz, and G. A. Dover. 1988. Complete sequences of the rRNA gene of Drosophila melanogaster. Mol. Biol. Evol. 5:366–376.[Abstract]

    Thompson, J. D., D. G. Higgins, and T. J. Gibson. 1994. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, positions-specific gap penalties and weight matrix choice. Nucleic Acids Res. 22:4673–4680.[Abstract]

Accepted for publication September 28, 1999.





This Article
Abstract
FREE Full Text (PDF)
Alert me when this article is cited
Alert me if a correction is posted
Services
Email this article to a friend
Similar articles in this journal
Similar articles in ISI Web of Science
Similar articles in PubMed
Alert me to new issues of the journal
Add to My Personal Archive
Download to citation manager
Search for citing articles in:
ISI Web of Science (10)
Request Permissions
Google Scholar
Articles by Kim, C.-G.
Articles by Osawa, S.
PubMed
PubMed Citation
Articles by Kim, C.-G.
Articles by Osawa, S.