The Japanese spined loach, Cobitis takatsuensis, has some unique morphological and ecological features among Japanese Cobitis species. Mitochondrial DNA analyses were conducted to investigate the magnitude of intraspecific differentiation and phylogenetic relationships among Japanese congeners of C. takatsuensis. PCR-RFLP (polymerase chain reaction-restriction fragment length polymorphism) analysis of the ND1 region with 16 restriction enzymes was carried out for thirteen C. takatsuensis populations throughout the species' distributional range. Populations in each river system possessed their own haplotypes, with marked genetic differentiation between the populations from Chugoku and Kyushu (ChugokuKyushu group) and those from Shikoku (Shikoku group). The two allopatric groups also showed different color pattern. Subsequently, sequencing analysis of part (725 bp) of the cytochrome b gene was carried out for C. takatsuensis and six other closely-related Japanese species. The phylogenetic tree indicated the presence of two major mtDNA lineages within Japanese Cobitis. It was noteworthy that the Chugoku-Kyushu and Shikoku groups of C. takatsuensis were included in separate mtDNA major lineages, and each group was closely related to other species. It is inferred that the distinct mtDNA relationship between the two allopatric C. takatsuensis groups is a result of the parallel evolution or mtDNA introgression, rather than divergence by geographic isolations.
The Japanese spined loach, Cobitis takatsuensis (Mizuno, 1970), is small species inhabiting only stony bottoms in mountain streams of western part of Japan (Kimizuka et al., 1982; Shimizu and Mizuno, 1994). It is rare species listed as ‘Threatened I’B rank in the 1997 Red List of Threatened Animals of Japan (Environment Agency of Japan, 1997). This species is morphologically characterized by thin caudal peduncle, incomplete sexual dimorphism, and unique color pattern (Mizuno 1970, 1975). These ecological and morphological features are peculiar among Japanese Cobitis species.
In addition, C. takatsuensis exhibits considerable morphological variation among local populations (Shimizu and Mizuno, 1994; Mizuno, 1995). In particular, marked differences in color pattern exist between populations from Chugoku District (western Honshu) and Kyushu Island, and those from Shikoku Island. The former were roughly represented by striated, and the latter by spotted types (Shimizu and Mizuno, 1994; Fig. 1). Although they have been regarded as intraspecific geographical variations, their taxonomic status remains obscure.
In this study, we conducted restriction fragment length polymorphism (RFLP) and sequence analyses on some parts of the mtDNA to investigate the magnitude of genetic differentiation between the two types of C. takatsuensis and their phylogenetic relationships among Japanese congeners, C. biwae, C. matsubarai, and three undescribed species of Cobitis (C. sp. 1, C. sp. 2, and C. sp. 3), which taxonomically follow as Hosoya (1993) (Table 1). Recently, Saitoh et al. (2000) reported the presence of three distinct mtDNA lineages within C. biwae, corresponding to geographical proximity (eastern part of Honshu, central-western part of Honshu, and Pacific slope of Shikoku). We considered them as distinct groups in C. biwae here.
Distribution and chromosome number of Cobitis and Niwaella species in Japan, according to Hosoya (1993).
MATERIALS AND METHODS
Specimens were caught by hand net during 1995 and 1999. Total of forty-seven individuals of C. takatsuensis were obtained from thirteen locations in twelve river systems (Fig. 2, Table 2). In addition, nine individuals of other species (C. biwae, C. matsubarai, C. sp. 1, C. sp. 2 subsp. 2, C. sp. 3, and Niwaella delicata) were included for inter-specific phylogenetic analysis. Three groups of C. biwae (sensu Saitoh et al. 2000), are herein temporarily referred to as Eastern, Western (consisted of diploid and tetraploid races (Kobayashi, 1976)), and Kochi groups, respectively. Accordingly, four individuals of C. biwae representing the three groups and a tetraploid race were included (Table 2). The identification of diploid or tetraploid specimens was determined from blood-smear preparations (Sezaki and Kobayashi, 1978). The remaining five individuals represented four other Cobitis species and a sister genus species Niwaella delicata as an outgroup. Saitoh et al. (2000) reported close genetic relationships among subspecies within C. sp. 2, and also that karyotype races within C. matsubarai. For simplification in the analysis, a single subspecies (subsp. 2) and a race (2n=94, inferred from collecting locality) each represented (Table 2). Samples were frozen immediately following collection, and stored at −20°C until analysis. Voucher specimens were deposited in the collection of the Fisheries Research Laboratory of Mie University (FRLM: catalogue numbers were listed in Table 2).
Sampling localities, sample sizes, and catalogue numbers. Locality numbers correspond to those in Fig. 2.
Total DNA was extracted from approximately 100 mg of frozen muscle tissue, as described by Asahida et al. (1996).
PCR-RFLP for intraspecific analysis
Restriction analysis was conducted for all specimens of C. takatsuensis.
A segment of about 2.0 kbp, containing a complete NADH dehydrogenase subunit 1 (ND1) and part of the 16S ribosomal RNA (16SrRNA) gene of mtDNA, was amplified by the polymerase chain reaction (PCR). A pair of primers (forward: 5′ - ACC CCG CCT GTT TAC CAA AAA CAT - 3′ and reverse: 5′ - GGT ATG AGC CCG ATA GCT TA - 3′) described by Hall and Nawrocki (1995) was used. However, because these were not effective for fish from Shikoku, another pair of primers, modified for carp (forward: 5′ - GCC TCG CCT GTT TAC CAA AAA CAT- 3′ and reverse: 5′ - GGT ATG GGC CCG AAA GCT TT - 3′), were used. PCR consisted of 30 cycles of 1 min. at 94 °C, 1 min. at 50°C and 2 min. at 72°C. Subsequently, PCR products were digested with sixteen restriction enzymes that specifically recognize four or five base nucleotides (Table 3; from New England Biolabs, Beverly, MA; and Takara Shuzo, Kyoto), following the protocols recommended by the manufacturers. The restricted fragments were separated by horizontal electrophoresis in 3% agarose gel. Digested fragments were visualized and photographed on an ultraviolet transilluminator after ethidium bromide staining. Nucleotide sequence divergences between mtDNA haplotypes were estimated according to Nei and Li (1979), using the D program from REAP (MacElroy et al., 1992). Phylogenetic trees were generated by the UPGMA method using the Neighbor program from PHYLIP ver. 3.572 (Felsenstein, 1996).
Fragment patterns and approximate molecular size (in base pairs; bp) generated by sixteen restriction enzymes.
Sequencing for interspecific phylogenetic analysis
Sequence analysis was performed on mtDNA fragments comprising part of the cytochrome b gene. Seven individuals of C. takatsuensis, representing each RFLP cluster, and the nine individuals of closely-related species (Table 2) were sequenced. A pair of primers (forward: 5′ - TGA CTT GAA RAA CCA YCG TTG - 3′ and reverse: 5′ - RGC RAA KAR RAA RTA YCA TTC - 3′) described by Palumbi et al., (1991) were used. PCR conditions for the sequencing analysis were similar to those for RFLP analysis, with the annealing temperature set at 54°C. Amplified DNA was purified with a QIA quick spin column (Qiagen, Germany). Direct sequencing of purified products was achieved using the Ready Reaction Dye Terminator Cycle Sequencing Kit (Applied Biosystems) following the manufacturer's protocol. Primers used were the same as those for PCR. Sequences were analyzed on an Applied Biosystems Model 377 automated DNA sequencer. Each DNA strand was determined the sequence two or three times. DNA sequence data were edited using DNASIS programs (Hitachi Software Engineer. Co. Ltd.). The nucleotide sequence data reported here are available from DDBJ, EMBL and GeneBank accession numbers AB039337– AB039352.
Three different methods were conducted to infer the phylogenetic relationships. A distance matrix was calculated based on Kimura's two-parameter method (Kimura, 1980) and clustered by the neighbor-joining method (Saitou and Nei, 1987) using MEGA ver. 1.0 (Kumar et al., 1993), and a maximum likelihood analysis performed using DNAML in PHYLIP ver. 3.572 (Felsenstein, 1996). A maximum parsimony analysis was performed with the Heuristic Search program from PAUP ver. 3.1.1 (Swofford, 1993). The robustness of each phylogeny was assessed by bootstrap analyses consisting of 1,000 replicates (Felsenstein, 1985).
PCR-RFLP for intraspecific analysis
The digested mtDNA fragment patterns are listed alphabetically for each of the 16 enzymes (Table 3). The total length of the digested DNA fragments varied in some cases, since the detection of small-sized fragments (< 70 bp) was not possible with the method used. The composite genotypes (haplotypes) detected among C. takatsuensis populations and the location of each are given in Table 4. All 16 enzymes produced variable restriction patterns. Restriction patterns detected in the populations from Chugoku and Kyushu were so different from those from Shikoku as to have hardly any fragments in common. Therefore, length-difference method (Nei and Li, 1979) cannot be applied for the analysis between them (Nei, 1987). Accordingly, we referred to them as Chugoku-Kyushu and Shikoku groups, respectively. The sequence divergences of the two groups were calculated separately.
Composite PCR-RFLP haplotypes for 13 populations of C. takatsuensis.
Within the Chugoku-Kyushu group, ten of the sixteen enzymes were polymorphic. A total of eight haplotypes (C1-7 and K1) were observed among the 31 fish analyzed (Table 4). Populations in each river system were fixed for their own haplotype (Table 4). The resulted UPGMA tree (Fig. 3. A) revealed some sub-clusters corresponding to geographical proximity. The estimated range of pairwise sequence divergence between haplotypes varied from 0.13 to 1.33%.
Within the Shikoku group, twelve of the sixteen enzymes were polymorphic. A total of six haplotypes (S1-6) were observed among the 16 fish analyzed (Table 4). The estimated range of pairwise sequence divergence among all haplotypes varied from 0.24 to 2.14%. Each population possessed it's own haplotypes (Table 4). In the UPGMA tree (Fig. 3. B), haplotypes detected in the same populations (S3 and S4, S5 and S6) were clustered.
Sequencing for interspecific phylogenetic analysis
A 725 base pair region of the cytochrome b gene was successfully sequenced for all specimens. No insertions or deletions were observed.
The NJ tree (Fig. 4) constructed from all three codon positions (725 bp) on the cytochrome b sequence indicated two major lineages of Cobitis supported by high bootstrap value (98% and 99%, respectively). Nucleotide divergence between the two major lineages averaged 17.21%. It was noteworthy that the haplotypes of C. takatsuensis from the ChugokuKyushu (C2, C4, C7, and K1) and Shikoku (S1, S3, and S6) groups were included in separate major lineages. Similarly, haplotypes of C. biwae from the Western group (W2 and W4) were included in different major lineages to those from the Eastern (E) and Kochi (K) groups. Cobitis takatsuensis from the Shikoku group (S1, S3, and S6) were clustered very closely with C. biwae from the Western group (W2 and W4). On the other hand, C. takatsuensis from the Chugoku-Kyushu group (C2, C4, C7, and K1) was associated with C. biwae from the Eastern group (E) with a high bootstrap value (95%).
Maximum likelihood analysis was conducted for the two individuals of C. takatsuensis (C7 and S1) representing each group, and the nine individuals of other species which were analyzed by the NJ analysis. The transition / transversion ratio was set at 4.6, following the method of Oohara et al. (1999). Although a change in the transition / transversion ratio value to accommodate a wide range (from 1 to 10) was attempted, the topology of the tree was not entirely changed. The topology of the ML tree (Fig. 5. A) was concordant with that constructed by the NJ method (Fig. 4). The two major lineages of Cobitis were completely (100%) supported by bootstrap replications. The each close relationship between C. takatsuensis (S1) and C. biwae (W2 and W4), and C. takatsuensis (C7) and C. biwae (E) were also strongly supported (100% and 96%, respectively).
The nucleotide substitutions at the third position of codon between each taxon were counted from 10 to 106 of the 242 sites, indicating saturation between distantly-related taxa. There was a possibility that multiple substitutions could not be adequately corrected in maximum parsimony analysis (Nei, 1987). Therefore, only first and second codon positions data set was used for MP analysis. The latter produced a topology concordant with those obtained by the NJ (Fig. 4) and ML (Fig. 5. A) methods, with the exception of an ambiguity in the branching order deep within the lineage I (Fig. 5. B).
The present result indicated the marked mtDNA divergence between C. takatsuensis populations from ChugokuKyushu and those from Shikoku, corresponding to color pattern differentiation proposed by Shimizu and Mizuno, (1994) (Fig. 1). RFLP haplotype divergences within the ChugokuKyushu and Shikoku groups were 0.13–1.33 % and 0.24– 2.14%, respectively. Since present RFLP data were estimated from a part of mtDNA, we could not directly compare to other population studies that have used whole mtDNA. Nevertheless, these values seem to be relatively high compared with other intraspecific divergence in freshwater species summarized by Bernatchez and Wilson (1998) (the maximum intraspecific divergence value < 1%, in the half of 42 species). Moreover, no shared haplotypes among river systems were observed. Within the each group, such considerable differentiations among populations must be attributed to the isolations for their exclusive mountain stream habitats (Kimizuka et al., 1982; Shimizu and Mizuno, 1994). On the other hand, it is unlikely that the differentiation between the two groups (19.24% sequence divergence in cytochrome b gene sequence) has resulted from such local isolations. If a conventional 2% divergence per Myr (Brown et al. 1979) is applied to this data, the two groups of C. takatsuensis could have diverged about 9–10 million years ago. This time estimation implies that the divergence occurred much long before the geographical isolation among western Chugoku, northern Kyushu and Shikoku Regions (the isolation time is shorter than 20,000 years, after Kaizuka, et al. 1995).
The present phylogenetic results provided some hints for distinct mtDNA compositions between the two allopatric groups of C. takatsuensis. The phylogenetic trees revealed the presence of two major mtDNA lineages within Japanese Cobitis. It is noteworthy that the two divergent groups within C. takatsuensis (Chugoku-Kyushu and Shikoku) were included in separate major Cobitis lineages, as well as three distinct groups of C. biwae (Eastern, Western, and Kochi groups). A simple, general explanation for the distinctive mtDNA composition within one species is possible random lineage sorting of ancestral polymorphic mtDNA (Billington and Herbert, 1991). In the present case, however, the divergence time of the two lineages (approximately 8–10 million years ago) calculated from the conventional vertebrate mtDNA clock (2% divergence per Myr; Brown et al. 1979) is much too long to have survived ancestral mtDNA polymorphisms (Avise, 1994). Accordingly, this explanation should be rejected.
Two alternative possible hypotheses exist satisfying the distinct mtDNA composition between the two groups. The first is that the two C. takatsuensis groups may have evolved independently from distinct lineages (Fig. 6. A). If this is correct, the two groups of C. takatsuensis must be regarded as two distinct species.
The second is that mtDNA introgression from another species to one of the C. takatsuensis groups may have occurred. The most parsimonious scenario under the introgression hypothesis is that mtDNA exchange between ancestors of Chugoku-Kyushu group C. takatsuensis and Western group C. biwae (W2 and W4) may have occurred (Fig. 6. B). In this scenario, two groups of C. takatsuensis can be considered as having evolved monophyletically. On the other hand, since C. matsubarai and C. sp. 1 were inferred as allotetraploid species whose maternal origins are C. biwae (Saitoh et al. 2000), it is reasonable that C. biwae is also considered as a monophyletic group. In fish species, interspecific mtDNA introgression, including a reciprocal case as in the present hypothesis, have been reported in several species (Avise and Saunders, 1984; Billington and Herbert, 1991; Dowling and Hoeh,1991; Dowling and DeMarais, 1993; Bernatchez et al., 1995, etc.). Moreover, evolutions originated from interspecific hybridization, being highly suggestive the possibility of mtDNA introgressions, had reported among several Cobitis species (Vasil'ev et al., 1989; Kim and Lee, 1990; Sezaki et al., 1994; Ráb and Slavík 1996; Saitoh et al. 2000). Published karyological data (Table 5) indicates that C. takatsuensis exhibits little karyotypic variations even between the two divergent groups (Kimizuka et al., 1982). Although C. biwae shows some variations, their karyotypes are characterized by unusual compositions among Cobitis species (high proportion of meta-centric chromosomes; Ueno et al., 1980; Ueno, 1981). The above seems to support the monophyletic evolutions of C. takatsuensis and C. biwae lineages, respectively.
Brief comparison of published data of karyotypes of C. takatsuensis and C. biwae.
However, it is impossible to distinguish whether parallel evolution or mtDNA introgression (or a combination of both) is responsible for the distinct mtDNA compositions of the two groups of C. takatsuensis, based solely upon mtDNA data. The two above hypotheses are presently the subject of ongoing nucleic DNA analyses.
We thank I. Oohara (National Research Institute of Fisheries Science), K. Kawamura (National Research Institute of Aquaculture), S. Kimura (Fisheries Research Laboratory of Mie University), and K. Saitoh (Tohoku National Fisheries Research Institute) for variable advice and assistance. We are grateful to G. S. Hardy (Thames, New Zealand) and two anonymous reviewers for their critically reading and improving the manuscript. We also thank Y. Fujioka (Shiga Prefectural Fisheries Experimental Station), I. Tanaka (Masuda City, Shimane Prefecture), D. Inoue (Kitakyushu High School) and N. Uechi (Mie University) for providing specimens used in the present study.
- 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
- J. C. Avise and N. C. Saunders . 1984. Hybridization and introgression among species of sunfish (Lepomis): analysis by mitochondrial DNA and allozyme markers. Genetics 108:237–255. Google Scholar
- J. C. Avise 1994. Molecular markers, natural history, and evolution. Chopman and Hall. New York. Google Scholar
- L. Bernatchez, H. Glémet, C. C. Wilson, and R. G. Danzmann . 1995. Introgression and fixation of Arctic char (Salvelinus alpinus) mitochondrial genome in an allopatric population of brook trout (Salvelinus fontinails). Can J Fish Aquat Sci 52:179–185. Google Scholar
- L. Bernatchez and C. C. Wilson . 1998. Comparative phylogeography of Nearctic and Palearctic fishes. Mol Ecol 7:431–452. Google Scholar
- N. Billington and P. D. N. Herbert . 1991. Mitochondrial DNA diversity in fishes and its implications for introductions. Can J Fish Aquat Sci 48:Suppl 180–94. Google Scholar
- W. M. Brown, M. J. George, and A. C. Willson . 1979. Rapid evolution of animal mitochondrial DNA. Proc Natl Acad Sci USA 76:1967–1971. Google Scholar
- T. E. Dowling and W. R. Hoeh . 1991. The extent of introgression outside the contact zone between Notropis cornutus and Notropis chrysocephalus (Teleostei: Cyprinidae). Evolution 45:944–956. Google Scholar
- T. E. Dowling and B. D. DeMarais . 1993. Evolutionary significance of introgressive hybridization in cyprinid fishes. Nature 362:444–446. Google Scholar
- Environment Agency of Japan 1997. Red list of threatened animals in Japan. Wildlife Protection Division, Nature Conservation Bureau, Environment Agency of Japan. Tokyo. (in Japanese). Google Scholar
- J. Felsenstein 1985. Confidence limits on phylogenies: an approach using the bootstrap. Evolution 39:783–791. Google Scholar
- J. Felsenstein 1996. PHYLIP (phylogeny inference package), version 3.572. University of Washington. Seattle. Google Scholar
- 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
- K. Hosoya 1993. Cobitidae. In “Fishes of Japan with Pictorial Keys to The Species”. Ed by T. Nakabo Tokai Univ Press. Tokyo. 231–234. (in Japanese). Google Scholar
- S. Kaizuka, Y. Naruse, Y. Ohta, and K. Koike . 1995. Nature in Japan 4: Plains and coasts of Japan. New ed Iwanami-Shoten. Tokyo. 248. (in Japanese). Google Scholar
- I. S. Kim and J. H. Lee . 1990. Diploid-triploid complex of the spined loach Cobitis sinensis and C. longicorpus (Pisces, Cobitidae). Korea J Ichthyol 2:203–210. Google Scholar
- Y. Kimizuka, H. Kobayashi, and N. Mizuno . 1982. Geographic distributions and karyotypes of Cobitis takatsuensis and Niwaella delicata (Cobitididae). Jpn J Ichthyol 29:305–310. Google Scholar
- M. Kimura 1980. A simple method for estimating evolutionary rate of base substitutions through comparative studies of nucleotide sequences. J Mol Evol 16:111–120. Google Scholar
- H. Kobayashi 1976. Comparative study of karyotypes in the small and large races of spinous loaches (Cobitis biwae). Zool Mag (Japan) 85:84–87. (in Japanese with English abstract). Google Scholar
- S. Kumar, K. Tamura, and M. Nei . 1993. MEGA: Molecular Evolutionary Genetics Analysis Version 1.01. Pennsylvania State Univ. Pennsylvania. Google Scholar
- D. MacElroy, P. Moran, E. Bermingham, and I. Konfield . 1992. REAP: an integrated environment for the manipulation and phylogenetic analysis of restriction data. J Heredity 83:157–158. Google Scholar
- N. Mizuno 1970. A new species of cobitid fish from Japan. Mem Ehime Univ Sci Ser B (Biol) 6:133–143. Google Scholar
- N. Mizuno 1975. On the notes of Cobitis takatsuensis. The Freshwater Fishes 1:80–82. (in Japanese). Google Scholar
- N. Mizuno 1995. Cobitis takatsuensis. In “Freshwater Fishes of Japan”. Ed by H. Kawanabe and N. Mizuno . Yamatokeikokusha. Tokyo. 394. (in Japanese). Google Scholar
- M. Nei and W. H. Li . 1979. Mathematical model for studying genetic variation in terms of restriction endnucleases. Proc Natl Acad Sci USA 76:5269–5273. Google Scholar
- M. Nei 1987. Molecular Evolutionary Genetics. Columbia University Press. New York. Google Scholar
- I. Oohara, K. Sawano, T. Okazaki, and T. Kobayashi . 1999. Reexamination of the molecular phylogeny of the masu salmon in the genus Oncorhynchus. Proc 5th Ind-Pac Fish Conf Noumea 1997Ed by B. Seret and J. Y. Sire . Soc Fr Ichthyol. Paris. 417–426. Google Scholar
- S. R. Palumbi, A. Martin, S. Romano, W. O. MacMillan, L. Stice, and G. Grabowski . 1991. The Simple Fool's Guide to PCR, Version 2. Univ Hawaii Zool Dept. Honolulu. Google Scholar
- P. Ráb and O. Slavík . 1996. Diploid-triploid-tetraploid complex of the spined loach, genus Cobitis in PŠovka Creek: the first evidence of the new species of Cobitis in the ichthyofauna of the Czech Republic. Acta Univ Carol (Biol) 39:201–214. Google Scholar
- K. Saitoh, T. Kobayashi, R. Ueshima, and K. I. Numachi . 2000. Analyses of mitochondrial and satellite DNAs on spined loaches of the genus Cobitis from Japan have revealed relationships among populations of three diploid-tetraploid complexes. Folia Zool 49:Suppl 19–16. Google Scholar
- 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
- K. Sezaki and H. Kobayashi . 1978. Comparison of erythrocytic size between diploid and tetraploid in spinous loach, Cobitis biwae. Bull Jpn Soc Sci Fish 44:851–854. Google Scholar
- K. Sezaki, S. Watabe, Y. Ochiai, and K. Hashimoto . 1994. Biochemical genetic evidence for a hybrid origin of spined loach, Cobitis taenia taenia, in Japan. J Fish Biol 44:683–691. Google Scholar
- T. Shimizu and N. Mizuno . 1994. Cobitis takatsuensis. In. “Report for The Project of Conservation of the Rare Aquatic Species”. Ed by. Japan Fisheries Resources Conservation Ass. Tokyo. 378–385. (in Japanese). Google Scholar
- D. L. Swofford 1993. PAUP: phylogenetic analysis using parsimony, version 3.1.1. Computer program distributed by the Illinois Natural History Survey. Champaign, Illinois. Google Scholar
- K. Ueno, S. Iwai, and Y. Ojima . 1980. Karyotypes and geographical distribution in the genus Cobitis (Cobitidae). Bull Jpn Soc Sci Fish 46:9–18. (in Japanese with English abstract). Google Scholar
- K. Ueno 1981. Karyotypes of cobitid fishes, mainly concerning about chromosomal polymorphism and polyploidy. Marine Science 13:60–70. (in Japanese). Google Scholar
- V. P. Vasil'ev, K. D. Vasil'eva, and A. G. Oshinov . 1989. Evolution of a diploidtriploid-tetraploid complex in fishes of the genus Cobitis (Pisces, Cobitidae). In “Evolution and Ecology of Unisexual Vertebrates”. Ed by R. M. Dawley and J. P. Bogart . New York State Museum. New York. 153–169. Google Scholar