A comparison of allozyme variation, restriction fragment length polymorphisms in the mitochondrial DNA and partial sequences of the ND II gene (496 bp) was made among two lacustrine populations of the piscivorous chub (Opsariichthys uncirostris uncirostris) in Japan, four fluvial populations in Korea, one lacustrine population in Russia and one specimen from the Amur River (O. u. amurensis). All analyses indicate that the Japanese populations of piscivorous chub are separable from the Asian mainland populations of Korea and Russia. The latter populations were further divided into the Korean and Russian fish. Although opinions are divided on the phylogenetic position of the population in Lake Mikata, Japan, which shows unique morphological traits intermediate between those of the population in Lake Biwa and the mainland populations, the current analysis indicates a closer relationship to the population in Lake Biwa.
INTRODUCTION
The piscivorous chub, a member of the cyprinid genus Opsariichthys, is widely distributed in East Asia. Its natural range extends from the Amur basin to the rivers that enter the Yellow Sea on the Korea Peninsula and from the Yangtze basin in China to the northern parts of Vietnam and Laos including Hainan Island (Berg, 1948; Uchida, 1939; Wu et al., 1964). The fish is also found in central Japan (Nakamura, 1969; Fig. 1). Banarescu (1968) recognized four subspecies in Opsariichthys uncirostris on the basis of the number of scales; O. u. uncirostris in Japan; O. u. amurensis in the Amur drainage, parts of North China and Korea; O. u. bidens in Central and South China; O. u. hainanensis in Hainan Island. Banarescu (1968) also notified that the Korean specimens are intergrades between O. u. uncirostris and O. u. amurensis.
The population in Lake Biwa is unique insofar as these fish have a larger number of scales on the lateral line than those in the other piscivorous populations (Nakamura, 1969). The fish distributed in the Amur basin, Korea and China have a small number of scales (Banarescu, 1968), while the number of scales in the population in Lake Mikata is intermediate between those in Lake Biwa and those on the Asian mainland (Tanaka, 1970). The above differences have led to some hypotheses as to the origin of the Japanese forms and their routes of dispersal from Korea. They include a sea route via the Sea of Japan and an overland route (Tsubokawa, 1988; Tomoda, 1989).
This study examines the genetic differentiation of the piscivorous chub that Banarescu (1968) suggested to be closely related. The samples were collected in Japan, Korea and Russia and examined using the data of allele-frequency, mitochondrial(mt) DNA restriction fragment length polymorphisms and mtDNA sequences, which could also elucidate the dispersal route to Japan.
MATERIALS AND METHODS
Specimens
Fish were caught with traps and hand nets in each locality from 1990 to 1996. Specific locations are shown in Fig. 1 and sample sizes are listed in Tables 2, 3 and 5. Opsariichthys u. uncirostris were collected in the Hasu River as it entered Lake Mikata, and O. u. amurensis were collected in the Somjin R., the Yongsan R., the Kum R., the Han R., the Amur R. and Lake Khanka. Fish collected respectively at three sites on Lake Biwa and at two sites on the Han River were pooled for analyses. Whole fish were frozen or preserved in ethanol immediately after collection until they were processed for analysis. Only analysis of mtDNA was performed in the case of the specimens from the Amur River and Lake Khanka, which were preserved in ethanol.
The dark chub, Zacco temmincki (collected in the Onga River, Fukuoka Prefecture) and the venus fish, Aphyocypris chinensis (collected in the Chikugo River, Fukuoka Prefecture), comprising subfamily Danioninae along with the piscivorous chub (Hosoya, 1993), were included as an outgroup for restriction and sequence analyses.
Number of lateral scales were counted from the first scale in contact with the pectoral girdle to the structural base of the caudal fin.
Allozyme analysis
Analysis on horizontal starch gels was performed as described by May et al. (1979) and Harris and Hopkinson (1976). Liver, skeletal muscle, eye ball and heart tissues were used, and thirty protein systems were examined (Table 1). The nomenclature of loci and genes followed that of Shaklee et al. (1990). If multiple loci code for a protein, the loci are numbered starting with the one which encodes the proteins with the least anodal migration. The most common allele at a given locus of fish collected in Lake Biwa was designated *100, and other alleles were assigned numbers according to the mobility of their homomeric protein products relative to that of the *100 allele.
Table 1
Electrophoretic analysis of allozymes from piscivorous chub collected in Japan and Korea
The genetic basis of banding variants was inferred by application of the following criteria: (1) the banding patterns had to be consistent with the known molecular structure of the protein; and (2) when a gene was expressed in more than one tissue, variant phenotypes had to be parallel among tissues (Grant et al., 1983).
Restriction analysis
Total DNA was extracted from approximately 500 mg of the white muscle tissue of each fish by the method of Asahida et al. (1996). Fragments of about 1.5 kilo base pairs (kbp) in length from the mitochondrial gene for tRNAPro to that for 12S rRNA, including the control region, were amplified by PCR (Saiki et al., 1988) from whole genomic DNA using primers PRO-L(CTACCTCCAACTCCCAAAGC) and 12SAR-H (ATAGTGGGGTATCTAATCCCAGTT; Palumbi et al., 1991). PCR was conducted for 30 cycles in a PJ-9600 apparatus (Perkin Elmer Cetus, CA) at 94°C for 0.5 min., 55°C for 0.5 min., and 72°C for 2.0 min. All products of PCR were confirmed to be of equal length by electrophoresis and were subjected directly to digestion with restriction endonucleases.
The products of PCR were digested according to recommendations of the manufacturer of enzymes. We used enzymes that recognize four or five nucleotides (Table 4; from New England Biolabs, Beverly, MA; and Takara Shuzo, Kyoto). Restriction fragments were separated by electrophoresis on horizontal 3.0% agar-ose gels. The bands of DNA were visualized for photography by staining with ethidium bromide and exposure to UV light. The restriction analysis was conducted to supplement sequence analysis.
Sequence analysis
Nine individuals representing each haplotype resulting from restriction fragment length polymorphism (RFLP) analysis, were selected for sequence analysis. A partial region of the mitochondrial ND II gene (496 bp following tRNAMet) which is considered one of the most variable region among protein coding genes in mtDNA (Meyer, 1994) was amplified by PCR using the following primer pairs: (5′-TAAGCTATCGGGCTCATACC-3′) and (5′-TTTGAAGGCTYTTGGTCTAA-3′). These were designed based on the reported primer (Hall and Nawrocki, 1995) and the sequences of mitochondrial DNA from carp (Chang et al., 1994) and human (Anderson et al., 1981). PCR was conducted for 30 cycles in the abovementioned apparatus at 94°C for 1.0 min., 52°C for 1.0 min., and 72°C for 2.0 min. Amplified DNA was resolved by electrophoresis on 2% agar-ose gels, stained with ethidium bromide and purified with a QIA quick spin column (Qiagen, Germany).
Sequence data were obtained using the dye-primer Taq cycle-sequencing reaction and an automated DNA sequencer (Applied Biosystem 377A). Sequences were determined three times. The nucleotide sequence data reported in this paper are available from DDBJ, EMBL and GenBank accession numbers AB050615- 050619.
Data analysis
Nei's unbiased genetic distance (Nei, 1978) values were estimated from allele-frequency data and a phenogram was generated by the unweighted pair-group method (UPGMA; Sneath and Sokal, 1973). Since an outgroup was not included in the analysis, we used this procedure.
Nucleotide sequence divergence was estimated by comparing the electrophoretic patterns of the two DNA sequences in question (length-difference method; Nei and Li, 1979) digested by a restriction endonuclease. DNA sequence data were processed using Higgins parameter in DNASIS programs (Hitachi Software Engineer. Co. Ltd. Japan). Genetic distances for sequence analysis were obtained using Kimura's two parameter method (Kimura, 1980). Phenograms were generated by neighbor-joining analysis (NJ; Saitou and Nei, 1987) and the reliability of each interior branch was tested by 1,000 bootstrap replications both using the MEGA 1.02 computer package (Kumar et al., 1993)
RESULT
Numbers of lateral scales
The distribution of the number of scales on the various specimens is shown in Table 2. The range for fish from Lake Biwa did not overlap that of fish from the Asian mainland. The average number of scales on fish from Lake Mikata was intermediate between that on fish from Lake Biwa and that from the Asian mainland, as described by Tanaka (1970). Significant differences (Welch's t test; P<0.01) were observed between Lake Mikata population and Lake Biwa population and between Lake Mikata and Asian mainland populations.
Table 2
Distribution of the number of scales on the lateral line for each population of the piscivorous chub
Allozyme analysis
Among 44 presumable loci surveyed, 14 loci seemed polymorphic in the piscivorous chub populations when we applied the criterion that a locus should be considered polymorphic in populations in which the frequency of the most common allele is less than 0.99. Though the PGM-1 locus was highly polymorphic, it was not routinely scored because of poor resolution. At all polymorphic loci, deviation of the observed number from the expected number was not significant (χ2 test on the assumption of a Hardy-Weinberg equilibrium).
Allelic frequencies at polymorphic and diagnostic loci for the collected specimens are listed in Table 3. The allelic frequencies of the Korean populations are taken from a previous report (Okazaki et al., 1993). The alleles in Japanese fish were replaced by alleles different from those in the Korean samples at three loci. By contrast, the genetic diversity within the Japanese populations was quite small. A significant difference in frequencies of alleles was barely observed at the Ldh-1 locus (G test, P<0.01; Sokal and Rohlf, 1995) and the Ah-2 locus (G test, P<0.05) when we compared populations from Lake Biwa and Lake Mikata. Two-allele polymorphism at the Ldh-1 locus was only observed in the population from Lake Biwa. In case of the Korean fish some characteristic traits were observed in the examined populations. The population from the Kum River was characterized by the unique allele Ah-2-160 which was found at high frequency. Populations from the Somjin and Yongsan Rivers were polymorphic for the Ada locus but no variation was seen in the remaining populations. The estimates of expected average heterozygosity (revised by reference to sample size; Nei, 1978) ranged from 1.2% to 3.2% (Table 3).
Table 3
Observed allelic frequencies at polymorphic and diagnostic loci for populations of piscivorous chub from Japan and Korea
A phenogram (UPGMA), generated on the basis of Nei's unbiased genetic distances (Fig. 2), indicates that the Japanese populations of piscivorous chub are separable from the Korean populations.
Restriction analysis
Seven out of 14 restriction endonucleases yielded multiple banding patterns (restriction morphs; Table 4). Samples from Japan and the Asian mainland were characterized by different morphs produced by Alu I, Hae III, ScrF I and Taq I. The specimens from the Amur River and Lake Khanka, both consisting a single river system, yielded unique morphs upon digestion by Hinf I, ScrF I and Taq I.
Table 4
Restriction enzymes used to examine piscivorous chub mt-DNA and their resulting fragment patterns
Table 5 shows 9 composite genotypes (haplotypes) detected among 53 piscivorous chub examined and their localities. The populations of Lake Biwa, Lake Mikata, the Kum River and the Hun River had their own haplotypes. The specimens examined from the Amur River and Lake Khanka shared a same haplotype. A phenogram drawn by the NJ procedure (Fig. 3-A) revealed two definite clusters, from Japan and the Asian mainland, as was also found after the analysis of allozymes (Fig. 2). The mainland cluster is clearly divided into the Korean samples and the specimens from the Amur River and Lake Khanka.
Table 5
Composite haplotypes for restriction enzyme polymorphisms detected in piscivorous chub,Zacco temmincki andAphyocypris chinensis, showing numbers and locales of fish with each haplotype
Sequence analyses
Four hundred and ninety-six base pairs of the nucleotide sequences following tRNAMet were determined from the 9 specimens of the piscivorous chub representing each RFLP haplotype and two specimens of the outgroup species. Fig. 4 summarizes the variable nucleotide positions.
A constructed phenogram by the NJ procedure (Fig. 3-B) as well as those obtained in the allozyme and in the restriction analyses revealed that the Japanese populations are separable from the Asian mainland populations of Korea and Russia. This phenogram also indicates that the Amur River specimen is separated from the Korean samples. The above is supported by high bootstrap values. The difference was observed in the branching order among haplotype-6, -7 and -8 between the restriction and sequence analyses. The probability of these relationships, however, did not receive high bootstrap support. In spite of the slight difference observed in the branching order among the Korean fish, the topology of UPGMA phenogram was almost same with that of NJ phenogram (data are not shown).
DISCUSSION
Differentiation between fish from Japan and the Asian mainland
Results obtained by analysis of both allozymes and mtDNA indicated that the Japanese populations of piscivorous chub are separable from the Asian mainland populations of Korea and Russia. The latter populations were further divided into the Russian and Korean fish. Although the topology of the Korean fish differed slightly among results of the three analyses, the phenograms that we obtained were almost identical (Figs. 2 and 3).
It is noteworthy that the mtDNA restriction analysis revealed that almost all populations examined had their own haplotypes. This result indicates that analysis of restriction site polymorphism, combining PCR and digestion of DNA, is a powerful tool for discrimination among river populations of this fish.
Originality of the Lake Mikata fish
In spite of the unique feature in the number of scales, the current analysis indicates a genetic similarity between Opsariichthys u. uncirostris (Banarescu, 1968) in Lake Mikata and that in Lake Biwa. Until recently, the endemic range of the piscivorous chub in Japan was confined to Lake Biwa, Lake Mikata and their adjoining rivers (Nakamura, 1969). Nowadays, however, this species is found in many other rivers and lakes as a consequence of transplantations from Lake Biwa incidental to transplantation of the Ayu fish, Plecoglossus altivelis, from Lake Biwa (Sunaga, 1980). Based on the above, Tomoda (1989) supposed that the piscivorous chub in Lake Mikata population was originated in Lake Biwa. However, the unique haplotypes observed in this lake (Table 5), which are clearly distinguishable from those observed in Lake Biwa, allows us to eliminate the possibility that the Lake Mikata population was derived from Lake Biwa population by recent transplantation. This is also supported by the absence of the Ldh-1-(-65) allele in the Lake Mikata population, which appeared exclusively in the Lake Biwa population at high frequency (Table 3).
Dispersal of piscivorous chub to Japan
It is considered that some organisms that are now native to Japan have close relationship to those on the Asian mainland (Lindberg, 1972). Freshwater fishes are believed to have penetrated Japan when Japan and the Asian Continent were joined by a land bridge during certain glacial periods (Lindberg, 1972). Since the Japanese freshwater fishes with the closest affinity to those on the Korean Peninsula are generally distributed in western Japan, it is presumed that the most of them came from the Korean Peninsula (Aoyagi, 1957; Mizuno, 1987).
The pattern of distribution of the piscivorous chub is quite unique, in that their natural distribution does not include western Japan (Nakamura, 1969), although populations have been established after they were transplanted to these areas (Sunaga, 1980). Tsubokawa (1988) and Tomoda (1989) hypothesized that the piscivorous chub never inhabited western Japan in the past and penetrated the areas around Lake Mikata directly from Korea through the Sea of Japan when it was a putative freshwater lake. However, geographic and paleontologic data indicate that the Sea of Japan was never a freshwater lake at least after the late Miocene (about 15 million years (MY) ago; Koizumi, 1977; Ujiie and Ichikura, 1977). The presumed timing of the divergence of the piscivorous chub between Japan and Korea, as deduced from the present sequence analysis (6.9–7.3% between them), is the first half of the Pliocene (4 MY ago) if we simply apply the vertebrate mtDNA clock (2% per MY (Brown et al., 1979). If we apply the figure of 0.76% per MY (Zardoya and Doadrio, 1999), which is the cyto-chrome b region of mtDNA clock examined for a number of cyprinid fishes, the divergence dates back to 10 MY ago. Considering a much faster molecular clock in ND II region used in this study compared to the cytochrome b region (Meyer, 1994), the timing of divergence may not go back as early as the period when the Sea of Japan was a putative freshwater lake though some errors are unavoidable in this method. The above indicates little possibility of a dispersal route from Korea through the Sea of Japan.
In tracing the dispersal route of the piscivorous chub, it is probably more reasonable to apply the accepted routes of other cyprinid fishes than to set up a specific hypothesis for this particular fish (Mizuno, 1987; Tomoda, 1989). It is possible to consider the piscivorous chub that penetrated western Japan survived in Lake Biwa and Lake Mikata alone through the sea level fluctuation. Similar supposition was proposed for some freshwater fishes in western Japan by Lindberg (1972).
Complete picture will emerge as O. u. hainanensis and O. u. bidens are included for further analysis in the future.
Acknowledgments
We thank Drs. N.Okada, M.Saneyoshi, V.Kharitonov, A. Brykov, V. Tsiger, H. Hosoya, K. Kawamura and T. Tabe for assistance in collecting the specimens used in this study. Grateful acknowledgment is to Dr. I. Oohara for his technical advice. We are also grateful to Dr. T. Seki for his interest and encouragement.