Open Access
How to translate text using browser tools
1 May 2002 Genetic Differentiation of Piscivorous Chub (genus Opsariichthys) in Japan, Korea and Russia
Toshio Okazaki, Sang-Rin Jeon, Tadao Kitagawa
Author Affiliations +
Abstract

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.

Fig. 1

Map of East Asia showing the sites where piscivorous chub were collected and the endemic range of this fish (shaded; Berg, 1948; Uchida, 1939; Wu et al., 1964; Nakamura, 1969). A) Lake Biwa (a;Mano, b;Sugoshi, c;Ado); B) L. Mikata; C) Somjin River (Yongdu-ri); D)Yongsan R. (Kumsong-ri); E) Kum R. (Changwi-ri); F) Han R.(a;Sewol-ri, b;Kahung-ri); G) Amur R (Troickoje).; H) L. Khanka.

i0289-0003-19-5-601-f01.gif

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

i0289-0003-19-5-601-t01.gif

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

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

i0289-0003-19-5-601-t02.gif

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

i0289-0003-19-5-601-t03.gif

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.

Fig. 2

UPGMA phenogram for the populations of the piscivorous chub of Japan and Korea based on genetic distances (Nei, 1978) using 43 loci.

i0289-0003-19-5-601-f02.gif

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

i0289-0003-19-5-601-t04.gif

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

i0289-0003-19-5-601-t05.gif

Fig. 3

NJ clustering based on restriction analysis (A) and NJ phenogram based on genetic distances estimated from the partial mitochondrial ND II gene sequences (B) of nine composite haplotypes of the piscivorous chub of Japan, Korea and Russia. Numbers at end of tree forks indicate haplotype number (shown in Table 5) and numbers at forks indicate bootstrap values (1,000 replicates).

i0289-0003-19-5-601-f03.gif

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.

Fig. 4

Variable nucleotide sites in haplotypes (shown in Table 5) from the aligned 496 bp nucleotide sequence of the ND II gene of the piscivorous chub, Zacco temmincki and Aphyocypris chinensis. Identity with the first sequence denoted by dots.

i0289-0003-19-5-601-f04.gif

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.

REFERENCES

1.

S. Anderson, A. T. Bankier, B. G. Barrell, M. H. L. de Bruijn, A. R. Coulson, J. Drouin, I. C. Eperon, D. P. Nierlich, B. A. Roe, F. Sanger, P. H. Schreier, A. J. H. Smith, R. Staden, and I. G. Young . 1981. Sequence and organization of the human mitochondrial genome. Nature 290:457–465. Google Scholar

2.

H. Aoyagi 1957. General Notes on the Freshwater Fishes of the Japanese Archipelago. Taishukan. Tokyo. in Japanese. Google Scholar

3.

T. Asahida, T. Kobayashi, K. Saitoh, and I. Nakayama . 1996. Tissue preservation and total DNA extraction from fish stored at ambient temperature using buffers containing high concentration of urea. Fish Sci 62:727–730. Google Scholar

4.

P. Banarescu 1968. Revision of the genera Zacco and Opsariichthys (Pisces, Cyprinidae). Vestn Cesko-Slov Spolecn Zool 32:305–311. Google Scholar

5.

L. S. Berg 1948. Freshwater Fishes of the U.S.S.R. and adjacent countries II. Zool Inst Acad Sci. English Transl Israel Prog Sci Transl Ltd. Jerusalem. Google Scholar

6.

L. Bernatchez and C. C. Wilson . 1998. Comparative phylogeography of Nearctic and Palearctic fishes. Mol Ecol 7:431–452. Google Scholar

7.

W. M. Brown, M. J. George, and A. C. Wilson . 1979. Rapid evolution of animal mitochondrial DNA. Proc Natl Acad Sci USA 76:1967–1971. Google Scholar

8.

Y. S. Chang, F. I. Huang, and T. B. Lo . 1994. The complete nucleotide sequence and gene organization of carp (Cyprinus carpio) mitochondrial genome. J Mol Evol 38:138–155. Google Scholar

9.

J. W. Clayton and D. N. Tretiak . 1972. Amine-citrate buffers for pH control in starch gel electrophoresis. J Fish Res Bd Canada 29:1169–1172. Google Scholar

10.

W. S. Grant, R. Bakkala, F. M. Utter, D. J. Teel, and T. Kobayashi . 1983. Biochemical genetic population structure of yellowfin sole, Limanda aspera, of the North Pacific Ocean and Bering Sea. Fish Bull 81:667–677. Google Scholar

11.

H. J. Hall and L. W. Nawrocki . 1995. A rapid method for detecting mitochondrial DNA variation in the brown trout, Salmo trutta. J Fish Biol 46:360–364. Google Scholar

12.

H. Harris and D. A. Hopkinson . 1976. Handbook of Enzyme Electrophoresis in Human Genetics. Am Elsevier Publ Co Inc. New York. Google Scholar

13.

K. Hosoya 1993. Cypriniformes. In “Fishes of Japan with Pictorial Keys to the Species”. Ed by T. Nakabo Tokai Univ Press. Tokyo. in Japanese. Google Scholar

14.

M. Kimura 1980. A simple method for estimating evolutionary rates of base substitution through comparative studies of nucleotide sequences. J Mol Evol 16:111–120. Google Scholar

15.

I. Koizumi 1977. Deep sea sediments and history of the Sea of Japan. Kagaku 47:45–51. in Japanese. Google Scholar

16.

S. Kumar, K. Tamura, and M. Nei . 1993. MEGA: Molecular Evolutionary Genetics Analysis, version 1.02. Pennsylvania State Univ. Univ Park. Google Scholar

17.

G. U. Lindberg 1972. Large-scale fluctuations of sea level in the Quaternary period: Hypothesis based on biogeographical evidence. Nauka, Leningrad. Japanese translation by T. Shinbori and F. Kanemitsu . 1981. Tokai Univ Press. Tokyo. Google Scholar

18.

C. L. Markert and I. Faulhaber . 1965. Lactate dehydrogenase isozyme patterns of fish. J Exp Zool 159:319–332. Google Scholar

19.

B. May, J. E. Wright, and M. Stoneking . 1979. Joint segregation of biochemical loci in Salmonidae: results from experiments with Salvelinus and review of the literature on other species. J Fish Res Bd Canada 36:1114–1128. Google Scholar

20.

A. Meyer 1994. DNA technology and phylogeny of fish. In “Genetics and Evolution of Aquatic Organisms”. Ed by A. R. Beaumont Chapman & Hall. London. Google Scholar

21.

N. Mizuno 1987. Formation of the Japanese freshwater fishes. In “Freshwater Fishes in Japan”. Ed by N. Mizuno and A. Goto . Tokai Univ Press. Tokyo. pp. 231–244. in Japanese. Google Scholar

22.

M. Nakamura 1969. Cyprinid fishes of Japan. Spec Publ Res Inst Nat Resources. Tokyo. in Japanese. Google Scholar

23.

M. Nei 1978. Estimation of average heterozygosity and genetic distance from a small number of individuals. Genetics 89:583–590. Google Scholar

24.

M. Nei and W. H. Li . 1979. Mathematical model for studying genetic variation in terms of restriction endonucleases. Proc Natl Acad Sci USA 76:5269–5273. Google Scholar

25.

K. Numachi, S. Nagahora, and M. Iwata . 1979. Genetic demonstration of hybrid between chum and pink salmon in the northwest Pacific. Rep Otsuchi Mar Res Cen Univ Tokyo 5:87–95. in Japanese. Google Scholar

26.

T. Okazaki, S. Jeon, and T. Kobayashi . 1993. Genetic variation and population structure in Korean piscivorous chub, Opsariichthys bidens from Korea. J Basic Sci Sang Myung Univ 7:1–10. in Korean. Google Scholar

27.

S. R. Palumbi, A. Martin, S. Romano, W. O. McMillan, L. Stice, and G. Grabowski . 1991. The Simple Fool's Guide to PCR, Version 2. Univ Hawaii Zool Dept. Honolulu. Google Scholar

28.

G. L. Ridgway, S. W. Sherburne, and R. D. Lewis . 1970. Polymorphism in the esterases of Atlantic herring. Trans Am Fish Soc 99:147–151. Google Scholar

29.

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

30.

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

31.

J. B. Shaklee, F. W. Allendorf, D. C. Morizot, and G. S. Whitt . 1990. Gene nomenclature for protein-coding loci in fish. Trans Am Fish Soc 119:2–15. Google Scholar

32.

P. H. Sneath and R. R. Sokal . 1973. Numerical Taxonomy. Freeman. San Francisco. Google Scholar

33.

R. R. Sokal and F. J. Rohlf . 1969. Biometry. Freeman. San Francisco. Google Scholar

34.

T. Sunaga 1980. Piscivorous chub. In “Freshwater Organisms in Japan”. Ed by T. Kawai, H. Kawanabe, and M. Mizuno . Tokai Univ Press. Tokyo. pp. 30–36. in Japanese. Google Scholar

35.

S. Tanaka 1970. Growth of a cyprinid fish “Hasu”, Opsariichthys uncirostris Temminck et Schlegel, in Lake Mikata. Bull Fac Education Toyama Univ 18:137–148. in Japanese. Google Scholar

36.

Y. Tomoda 1989. The transition of Lake Biwa. Aoki Shoten. Tokyo. in Japanese. Google Scholar

37.

K. Tsubokawa 1988. Zoogeographical consideration on the freshwater fish fauna of the Okayama Region, Japan. Bull Kurashiki Mus Nat Hist 3:1–30. in Japanese. Google Scholar

38.

K. Uchida 1939. The Fishes of Tyosen (Korea) I. Bull Fish Exp St, Government -General Tyosen. 6:1–458. in Japanese. Google Scholar

39.

H. Ujiie and M. Ichikura . 1977. Quaternary environment of the Sea of Japan. Kaiyo-kagaku 9:316–322. in Japanese. Google Scholar

40.

H. Wu, et al 1964. The cyprinid fishes of China I. Sci Tech Press. Shanghai. in Japanese. Google Scholar

41.

R. Zardoya and I. Doadrio . 1999. Molecular evidence on the evolutionary and biogeographical patterns of European cyprinids. J Mol Evol 49:227–237. Google Scholar
Toshio Okazaki, Sang-Rin Jeon, and Tadao Kitagawa "Genetic Differentiation of Piscivorous Chub (genus Opsariichthys) in Japan, Korea and Russia," Zoological Science 19(5), 601-610, (1 May 2002). https://doi.org/10.2108/zsj.19.601
Received: 18 April 2001; Accepted: 1 March 2002; Published: 1 May 2002
KEYWORDS
Asian mainland
genetic differentiation
Japan
mt-DNA
Opsariichthys
Back to Top