The Japanese silver crucian carp, Carassius langsdorfii, is a naturally polyploid fish that reproduces by gynogenesis, resulting in a geographically clustered clonal population. To determine the amount of clonal diversity in this species, we used genetic markers to compare fish from several Japanese locations. Samples were collected from Lake Kasumigaura (Ibaraki Pref.), Lake Imba (Chiba Pref.), Lake Fukushimagata (Niigata Pref.), Lake Biwa (Shiga Pref.), Lake Koyama (Tottori Pref.), and the Niyodo and Monobe Rivers (Kochi Pref.). The genetic markers used were the microsatellite loci GF1*, GF17*, and GF29*, as well as isozyme (creatinekinase CK*) and mtDNA (D-loop region) RFLP. In the seven locations sampled, 61 clonal lines were identified by combining the genotypes of the three microsatellite loci. Populations at all locations consisted of six to 28 clonal lines. Fourteen common clonal lines were detected in two or more locations, and 13 of these were observed in Lake Biwa. We suggest that these clonal lines have been widely distributed in Japanese freshwaters through the intentional stocking of fish species such as “gengorobuna” (C. cuvieri) and “ayu” (Plecoglossus altivelis altivelis) obtained from Lake Biwa and its adjacent waters. The CK* genotypes and mtDNA haplotypes of fish from Lake Kasumigaura differed from those of other locations.
Unisexual fishes are usually polyploid, and they frequently reproduce by a non-bisexual reproductive mechanism such as gynogenesis or hybridogenesis (Dawley, 1989). The Japanese crucian carp is one such polyploid fish that reproduces by gynogenesis.
Studies of the various reproductive modes in unisexual fishes, including the cytological mechanisms that underlie them, occasionally make taxonomic revision necessary. It has until now been accepted that the Japanese silver crucian carp is represented by diploid (found only in western Japan), triploid and tetraploid forms (Kobayasi, 1971), the latter two forms reproducing gynogenetically. The diploid crucian carp “kinbuna”, Carassius buergeri subsp. was thought to be distributed only in north-east Japan (Nakamura, 1969), although a western Japanese diploid bisexual form, known as “Okinbuna”, C. buergeri subsp., (Ochiai et al. 1979), was found to have highly similar electrophoretic and morphological characteristics (Taniguchi and Sakata, 1977). The triploid silver crucian carp produces unreduced triploid eggs (Yamashita et al., 1993), meaning that the off-spring are usually genetically identical to the parent (Dong and Taniguchi, 1996; Umino et al., 1996). The triploid silver crucian carp is thus genetically independent from the diploid form (Ohara et al., 2000). On the basis of these findings, we consider the Japanese silver crucian carp, C. langsdorfii (Matsubara and Ochiai, 1965), is a naturally polyploidy fish that reproduces gynogenetically (Yamashita et al., 1993; Dong and Taniguchi, 1996).
Populations of some unisexual fishes have been shown to consist of several clonal lines (Moore, 1984). Wild Japanese silver crucian carp also consist of several genetically divergent clonal lines (Dong et al., 1996; Umino et al., 1997; Ohara et al., 1998). Large geographically-determined genetic divergences have also been reported (Taniguchi and Sakata, 1977; Shimizu et al., 1993); however, some clonal lines have been found to be present at distant locations in Japan (Ohara et al., 2000). A comprehensive study of the distribution of clones throughout Japan has not yet been carried out.
Microsatellite DNA polymorphism might be the most sensitive marker available for the identification of clonal lines by determining the zygotic condition and genotype of each locus (Ohara et al., 1998). Mitochondrial DNA analysis is not very useful for identifying clones, but it can be effective for demonstrating maternal origins (Murakami et al., 2001); therefore mtDNA analysis might be applicable in studies of clonal evolution. It is also known that creatinekinase (CK*) polymorphism is an effective marker for detecting genetic variation among different geographic locations in the Japanese silver crucian carp (Taniguchi and Sakata, 1977; Dong and Taniguchi, 1996).
In this study, we investigated the clonal component, the distribution of clones, of samples collected from six locations around Japan using microsatellite DNA polymorphism, RFLP of mitochondrial DNA and electrophoretic patterns of CK*.
MATERIALS AND METHODS
Sample collection and ploidy determination: We collected 74 fish samples from Lake Biwa in the Shiga Pref. (L. Biwa), 31 samples from Lake Fukushimagata in the Niigata Pref. (L. Fuku), 71 samples from Lake Kasumigaura in the Ibaraki Pref. (L. Kasumi), and 45 samples from Lake Imba in the Chiba Pref. (L. Imba) (Fig. 1). We also used previously published data to compare clonal components with each other (Ohara et al., 1999; 2000). From these earlier reports, we quote sixteen clonal lines (KOC-001∼−015, −017) in 237 samples from two rivers, the Niyodo and Monobe Rivers in the Kochi pref. (Kochi) (Ohara et al., 1999), and six clonal lines (KOY-001∼006) in 29 samples from Lake Koyama (L. Koyama) in the Tottori Pref. (Ohara et al., 2000). Because the distance between the two rivers in the Kochi prefecture was relatively close, we pooled the data from them, denoting them by the name Kochi. We determined the ploidy level of each sample by analyzing the major diameter of the erythrocyte (up to 15 μm) using the methods of Sezaki et al. (1977) and Onozato et al. (1983). As individuals of crucian carp verified as triploid were defined to be Japanese silver crucian carp, Carassius langsdorfii (Dong et al., 1997), we used only the Japanese silver crucian carp (triploid) to genetic analysis in this study.
Microsatellite DNA marker loci: Extraction of DNA was performed according to the methods of Takagi et al. (1997). Microsatellite primers, which were developed by Zheng et al. (1995), were used to detect the loci GF1*, GF17*, and GF29* in this experiment. The reverse primer had the 5′-end labeled with biotin. The PCR was programmed for seven cycles of 1 min at 94°C, 30 sec at 53°C, and 30 sec at 72°C, and 33 cycles of 30 sec at 90°C, 30 sec at 53°C, and 30 sec at 72°C. Following amplification, the PCR product was mixed with denaturing stop dye, heated at 95°C for 15 min, and electrophoresis was performed using a 6% denatured polyacrylamide gel. Chemiluminescence detection was carried out according to Perez-Enriquez et al. (1998). After electrophoresis, DNA was transferred to a nylon membrane by blotting, and the membrane was then dried and UV cross-linked. DNA on the membrane was detected using a PhototopeTM – Star Detection Kit (New England Biolabs). The sequence ladder obtained from the M13 was used as a size marker, and was prepared using a PhototopeTM CircumVent Kit (New England Biolabs). The nomenclature of loci and alleles follow that of Shaklee et al. (1990).
Identification of clonal lines and estimation of genetic diversity: When one or more individuals belonged to one combined genotype of the three microsatellite loci, we recognized it as an independent clonal line. In the two band types, one allele was recognized to be duplicated in one of the two bands (Ohara et al., 1999). The clonal diversity (D) was estimated by the formula, D = 1 - ∑ ni (ni - 1) / N (N - 1) (Moore, 1984), where ni is the number of individuals of a given clone and N is the total number of individuals. To estimate the number of different clonal lines at each location; the octave method (Ohtsuka and Tsuji, 1997) was used.
Isozymes: Muscle samples were preserved in a freezer at −20°C. For isozyme analysis, polymorphism of creatinekinase (CK; EC: 184.108.40.206) was detected by horizontal starch-gel electrophoresis. The designation of alleles and genotypes followed Dong and Taniguchi (1996).
mtDNA RFLP: Analysis of mtDNA RFLP followed Ohara et al. (1998). The region of approximately 2.1 kbp in length, containing the whole D-loop region of the mtDNA, was amplified by PCR. This region contains part of the cytochrome b gene and the 12SrRNA gene region. The primer sequences used were L15556 and H1067. The restriction endonucleases used in this study were Hinf I, Rsa I, Mbo I, and Taq I. A composite mtDNA haplotype, consisting of four letters and representing the fragment pattern generated by each of the endonucleases, was compiled for each individual. The haplo-type diversity h = 2n (1–∑xi2) / (2n - 1) was calculated (Nei, 1987). Differences among the haplotypes were determined using the program FreeTree (Pavlíek et al., 1999). Genetic distances (d) of haplotypes were computed from Nei-Li's coefficient of similarity (s), as d = 1 - s. The dendrogram was constructed by the UPGMA and bootstrap values (for 1000 resamplings) were computed for every node.
Microsatellite analysis: Combining the genotypes of the three microsatellite loci investigated, 61 clones were identified from the six locations (Tables 1 and 2). A total of 39 clonal lines were newly discovered. Fifteen unique clonal lines (designated as BIW-001 to -015) were detected in L. Biwa, eight (NII-001 to -008) in L. Fuku., ten (KAS-001 to -010) in L. Kasumi, and six (IMB-001 to -006) in L. Imba (Table 2). Fourteen common clonal lines were detected in two or more locations and designated as common clones COM-001 to -014. The common clonal lines KOC-001, -002, -004, -005, -006, -008, -010, -013, -014, -015, -017, and KOY-003, that had been recognized previously were renamed as COM-003 to -014, respectively (Table 1) (Ohara et al., 1999; Ohara et al., 2000). COM-005 and COM-007 were found in five locations, and COM-009 and COM-013 in four locations. The four clones had a common allele, *216 in GF17* and genotype *abc in *CK. Thirteen of the 14 common clonal lines were observed in L. Biwa, and 11 in the Kochi prefecture samples.
The genotypes in the three microsatellite loci, genotypes of CK* and haplotype of mtDNA of 14 common clonal lines in Japanese silver crucian carp collected from six locations
The genotypes in the three microsatellite loci, genotypes of CK* and haplotype of mtDNA of 47 unique clonal lines in Japanese silver crucian carp collected from six locations
The level of genetic variability observed in the microsatellite DNA is summarized in Table 3. The total numbers of clones differ by location, ranging from six in L. Koyama to 28 in L. Biwa. Genotype frequencies of triallelic and diallelic fish in all locations were almost the same, and homozygous genotypes were rarely found. The clonal diversity is distributed from 0.559 in L. Kasumi to 0.949 in L. Biwa (Table 3). The total number of clonal lines expected was 40.53 in L. Biwa, 17.33 in the Kochi pref. and 18.03 in L. Imba (Table 3); other locations were in disagreement with a log-normal distribution because of the insufficient number of samples. Isozymes: Three alleles, *a, *b, and *c, and six genotypes, *aab, *abb, *abc, *bbb, *bbc, and *bcc were found in CK* (Table 3). Individuals having identical genotypes in their microsatellite loci also possessed the same genotypes in CK*, and each of the six genotypes included many clones. The most common genotypes differed among L. Kasumi (*bbb), L. Imba (*aab) and the other three locations (*abc) (Table 3). Eleven of the 14 most common clones had genotype *abc.
Genetic diversity in three microsatellite loci, CK* and mtDNA RFLP in Japanese silver crucian carp collected from six locations
mtDNA RFLP: In the mtDNA RFLP analysis, the numbers of detected fragment patterns were five in Hinf I, three in Rsa I, six in Taq I and three in Mbo I (Table 3). The haplo-types were determined by combining these fragment patterns. The eight haplotype designations, #1 to #8, were the same as those described by Ohara et al. (1998). Five new haplotypes, #13 to #17, were found in this study, and their fragment patterns are shown in Table 3. The haplotypes #9 to #12 were already listed in a previous publication (Ohara et al., 2000). The individuals having identical genotypes in their microsatellite loci also have the same mtDNA haplo-type, and each haplotype except for #2, #5, and #7 contained many clones. The haplotype diversities in six locations ranged widely, from 0.067 to 0.787 (Table 3), being extremely low in L. Koyama. The haplotype frequencies differed by location: haplotypes #1 and #6 were found in all locations, and the frequency of haplotype #6 was highest in L. Koyama (0.966) and L. Fuku (0.581). The frequencies of haplotypes #13 and #16 were highest in L. Kasumi (0.634) and L. Imba (0.518), respectively. Haplotype #13 was considerably different from the other eleven haplotypes, and this branch had again a bootstrap value of 100% (Fig. 2).
Efficiency of microsatellite DNA in the identification of clonal lines
Sixty-one clonal lines were detected by the combined genotypes of microsatellite DNA markers. The validity of their clonal nature and the effectiveness of these markers for the identification of clones were shown by the fact that the individuals having identical microsatellite genotypes belonged to the same CK* genotypes and mtDNA haplo-types. The number of clonal lines observed in the Japanese silver crucian carp was rather high in comparison with other gynogenetic fishes, such as species of Poeciliopsis, Menidia, and Phoxinus, although the index of clonal diversity was similar to that of other unisexual species (Moore, 1984; Echelle et al., 1989; Elder and Schlosser, 1995). The observed numbers of clones in L. Biwa and L. Imba were 70% of the expected total number of clonal lines, whereas the observed and expected numbers of clones were almost identical in the Kochi prefecture. However, we used only three microsatellite loci in this study. If we were to use more loci, there is a possibility that we might identify many more clones.
Genetic divergence of Japanese silver crucian carp
In an earlier study (Taniguchi and Sakata, 1977), the genotypes *bbb of CK* was highly prevalent in L. Kasumi, whereas most of the other locations except for the Kanto district consisted of *abc genotypes. In this study, the *bbb genotype was again mostly detected in L. Kasumi, and *abc was mostly detected in western Japan. Shimizu et al. (1993) classified Japanese crucian carp into three groups, the “Pacific”, “Sea of Japan”, and “Western” types. These results support the hypothesis that Japanese silver crucian carp in the Kanto district are different from those of other locations.
It is believed that the Japanese silver crucian carp originated from hybridization between two unknown species (Taniguchi, 1974; Shimizu et al., 1993). Murakami et al. (2001) suggested that it consisted of two lineages with two different maternal ancestors. In this study, mtDNA haplotype #13 was divergent from the other haplotypes (Fig. 2). Because the four clonal lines with haplotype #13 did not have the allele *a of CK*, haplotype #13 must have originated from a different maternal ancestor.
Our findings also suggest that the Japanese silver crucian carp is at least diphyletic. The parental species have not been determined (Vrijenhoek et al., 1989; Shimizu et al., 1993). Allele *c of CK* has never been found in diploid crucian carp in Japan (Taniguchi and Sakata, 1977). Allele *B of AMY-2* (Shimizu et al., 1993) and the tandem repeated DNA sequence (Cal3nDr) (Murakami and Fujitani, 1997) are also specific to Japanese silver crucian carp. Therefore, it is possible that the silver crucian carp did not arise in Japan and that its roots might be traceable to other East Asian countries.
Historical distribution of the common clones
In this study, 14 common clonal lines were detected in two or more (up to five) locations. Echelle et al. (1989) reported that identical clonal lines of the unisexual genus Menidia were naturally distributed at several different locations along the shores of the Gulf of Mexico. In contrast, in the selfing hermaphroditic fish Rivulus marmoratus, identical clonal lines were not detected in two different years at the same location (Turner et al., 1992). There are at least two possible explanations for the presence in different locations of common clones of the Japanese silver crucian carp. The first is that common clones become naturally distributed over a long evolutionary period. The second possibility is that the common clones were distributed through anthropogenic means.
As to the second historical process, intentional stocking of freshwater fish species in the most recent 50 years could have accelerated the spread of distribution ranges. Fish species such as Cyprinus carpio, Hypomesus transpacificus, and Plecoglossus altivelis altivelis have been intentionally introduced to ponds and rivers in Japan to enhance fish resources over the last 50 years. This intentional stocking has caused spreads in the distributions of non-targeted species such as Rhodeus ocellatus ocellatus and Zacco platypus (Hosoya, 1993). In the genus Carassius, a strain called “kawachibuna” or “herabuna”, improved from “gengorobuna” (C. cuvieri), which was originally endemic to L. Biwa, has now been stocked throughout Japan (Hosoya, 1993). Eighteen million C. cuvieri individuals are released annually in Japan. We also found this species in the six locations described in this paper. In this study, 13 of the 14 common clonal lines were found in L. Biwa, and this result supports the notion that the common clones could have been accidentally introduced to locations through the stocking of C. cuvieri into irrigation ponds and dams throughout Japan. There is no clear evidence for the natural range expansion of common clones; however, they were spread with the stocking seeds of C. cuvieri from L. Biwa and its adjacent waters during the last 50 years.
In studies of the population structure of any Japanese freshwater fish species, it is necessary to take into account the possible effects of human-induced disturbance, both deliberate and accidental. The recognition in this study of both common and unique clones of Japanese crucian carp should provide insights that will be of undoubted value in further studies on the origins and mechanisms of clonal divergence of this species.
We would like to thank all the people who helped us collect fish samples: M. Maehata of the Lake Biwa Museum, N. Watanabe and K. Suzuki of the Ibaraki Prefectural Freshwater Fisheries Experimental Station, G. Yamashita of the Chiba Prefectural Freshwater Fisheries Experimental Station, K. Sitizyo of Tottori University, and S. Sato of the Niigata Prefectural Freshwater Fisheries Experimental Station. We thank Dr. M. J. Grygier of the Lake Biwa Museum and Dr. M. Ikeda of Tohoku University for their critical reading of the English manuscript.