We investigated the genetic structure of Korean and Japanese ayu
populations by examining 669 individuals from 14 populations using three microsatellite loci. Genetic variation did not differ significantly among the populations examined in terms of allelic number and heterozygosity. Korean populations were genetically close to each other, implying that persistent gene flow has occurred in these populations. This suggests that eastern populations in Korea form a single large population and all of the Korean populations are distinct from the Japanese populations. Pairwise population
estimates, principal component analyses, and a neighbor-joining tree showed that genetic separation between the southern and pooled eastern coast populations was probably influenced by restricted gene flow. Hierarchical analysis of molecular variance (AMOVA) revealed a weak but significant genetic structure among three ayu groups (eastern and southern coasts of Korea and the Japan coast), and no genetic variation within groups. The estimated genetic population structure and potential applications of microsatellite markers may aid in the proper management of ayu populations.
, is widely distributed in Ko-rea and Japan and is an ecologically important inland fish (Han et al., 2003). Two different ecological forms of ayu, an amphi-dromous form that normally migrates between rivers and the sea and a landlocked form, have different life histories. Both forms exist in Japan, whereas only the amphidromous form is found in Korea (Iguchi et al., 1999; Ikeda and Taniguchi 2002). Recently, the number of ayu returning from the sea has declined, possibly because of environmental degradation in rivers caused by industrial and uncontrolled development. Therefore, knowledge of wild populations is essential for ef-fective natural resource management and the conservation of native aquatic biodiversity (Ryman et al., 1995). Genetic variation is important for the long-term survival of natural populations because it confers the ability to adapt to environ-mental changes, thereby increasing fitness (Frankel and Soulé, 1981). A lack of genetic variation caused by inbreeding can be detrimental to fitness. Estimates of genetic variation in wild populations, monitored using appropriate molecular markers, are important for preventing undesirable changes in produc-tion. Hence, the biological and genetic characteristics of ayu populations should be evaluated to maintain genetic variation.
To date, isozymes have been widely used as markers in studies of ayu population genetics. Taniguchi et al. (1983) studied genetic variability and differentiation among amphi-dromous, landlocked, and hatchery populations of ayu in Ja-pan. Seki and Taniguchi (1985) and Nishida (1985) studied genetic divergence among amphidromous ayu populations in Japan. Although the variability of isozymes is beneficial for population genetic analyses of ayu, isozyme analysis re-quires careful collection and handling of tissues (Park et al., 1993). Furthermore, the resolution of isozymes is most effec-tive at regional levels (Han et al., 2003). Han et al. (2003) found substantial gene flow that was sufficient to genetically homogenize 11 natural Korean ayu populations. In addition, a limited number of studies have applied genetic analyses to Korean amphidromous ayu populations. Seki et al. (1988) and Sawashi et al. (1998) showed genetic divergence between Ko-rean and Japanese ayu populations, but they only studied four populations in Korea.
Microsatellites are highly polymorphic nuclear loci that have been used successfully in studies of population genet-ics, pedigree analysis, parentage assignment, and linkage mapping. Among the many types of DNA markers, micro-satellites are particularly useful because they are evenly dis-tributed in genomes, have a codominant Mendelian manner of inheritance, and are easily genotyped via PCR. Takagi et al. (1999) demonstrated the great potential of microsatellites as indicators of genetic variability and divergence among ayu populations, finding higher levels of polymorphism than were obtained during previously isozyme analyses.
The present study investigated genetic variation and popu-lation structure in natural populations of
collected from Korea and Japan by analyzing microsatellite loci.
Materials and Methods
- Fish samples
Samples of ayu were collected from 10 rivers located in eastern and southern Korea in 1998 (
). Ayu samples were also collected from the Namdae River and the Wangpi River in 1997. The Kochi River and the Biwa River populations in Japan were studied previously by Takagi et al. (1999) and were compared with the Korean populations. Wild fish were caught at a single location at each site within a few days using a pot, frozen with dry ice, and stored at -20℃ until use.
- DNA extraction and microsatellite genotyping
For each ayu sample, DNA was extracted from a fin-clip following a slight modification of the methods described by Taggart et al. (1992). Fin tissue was placed in 700 μL TNES-Urea (10 mM Tris-HCl pH 7.5, 1.5 M NaCl, 10 mM EDTA, 0.5% sodium dodecyl sulfate, and 4 M Urea) and 5 μL of proteinase K (50 μg/μL final concentration). The mixture was then shaken gently and incubated overnight at 37℃. DNA was purified by successive extractions with phenol : chloroform : isoamylalchol (25:24:1) and chloroform-:-isoamylachol (24:1), respectively. DNA was precipitated with 3 M sodium acetate trihydrate and a double volume of 99% cold ethanol. The precipitate was decanted, washed with 70% ethanol, and air-dried. The DNA pellet was resuspended in 100 μL TE buf-fer (10 mM Tris-HCl, 1 mM EDTA pH 7.2) and stored at 4℃ prior to PCR analysis.
Sampling locations of 14 ayu populations analyzed in this study (see Table 1 for site names).
Sample sites sample number and sampling date of ayu in the present study
Sample sites sample number and sampling date of ayu in the present study
For the microsatellite analysis, three primers,
), were screened using the annealing temper-atures and PCR cycles described by Takagi et al. (1997). The forward primer from each primer set was 5-fluorescent labeled with one of three dyes: 6-FAM, HEX, or NED (PE Applied Biosystems, Foster City, CA, USA . PCR amplification of six microsatellite loci was conducted using an RTC 200 instrument (MJ Research, Wiltham, MA, USA in 10 mL of solution con-taining 10-50 ng DNA, 1× ExTaq buffer, 0.2 mM dNTPs, 10 pmol of each primer, and 0.25 U Taq DNA polymerase (Takara, Ohtsu, Japan . The amplification protocol included an initial de-naturation for 11 min at 95℃ followed by 35 cycles of 1 min at 94℃, 1 min at the optimal annealing temperature (the anneal-ing temperature for each locus is listed in
), and 1 min at 72℃, with a final extension step of 5 min at 72℃. The sizes of fluorescence-labeled allele fragments were measured on an ABI PRISM 3130XL automated sequencer, followed by analy-sis with GeneMapper version 3.7 (Applied Biosystems).
- Data analyses
The genetic diversity of each location was estimated by the number of alleles per locus and observed (Ho) and expected (He) heterozygosities, which were calculated using FSTAT ver-sion 2.9.3 (Goudet, 2001) and GENEPOP version 1.2 (Raymond and Rousset, 1995). The inbreeding coefficient, F
, was calcu-lated in an analysis of variance framework following Weir and Cockerham (1984) using GENEPOP. Departure from Hardy-Weinberg equilibrium (HWE) and linkage disequilibrium were calculated using GENEPOP version 1.2 (Raymond and Rousset, 1995). Tests for the occurrence of null alleles were performed with MICRO-CHECKER version 2.2.3 (Van Oosterhout et al., 2004). Pairwise
values were used to estimate genetic differ-entiation between population pairs according to Slatkin (1995) using FSTAT. An analysis of molecular variance (AMOVA) was employed to define the grouping of genetic variation in hierarchi-cal arrangements using Arlequin version 3.11 (Excoffier et al., 2005). Genetic relationships among populations were assessed by principal component analysis (PCA) based on the covariance matrix of gene frequencies using PCA-GEN 1.2 (Goudet, 1999).
Nucleotide sequence of 3 microsatellite PCR primers repeat motif and amplification condition in Korean and Japanese populations*Core repeat motif cloned Plecoglossus altivelis,†PCR annealing temperature were optimized for Plecoglossus altivelis,‡Estimated size of the PCR fragment when compared with M13 sequence fragment of known length.
Nucleotide sequence of 3 microsatellite PCR primers repeat motif and amplification condition in Korean and Japanese populations *Core repeat motif cloned Plecoglossus altivelis, †PCR annealing temperature were optimized for Plecoglossus altivelis, ‡ Estimated size of the PCR fragment when compared with M13 sequence fragment of known length.
In addition, after constructing a genetic distance matrix based on a set of gene frequencies in different populations, which was estimated according to Nei (1972), a NJ tree was constructed for each replicated genetic distance matrix via bootstrapping of 1000 replications using NEIGHBOR in the PHYLIP version 3.5 software package (Felsenstein, 1993).
Results and Discussion
We examined 669 individuals from 14 ayu populations in Korea and Japan using three microsatellite DNA markers. Al-lele size in base pairs (S), the total number of alleles (A
), and observed (H
) and expected (H
) heterozygosities for the three loci (
) for each population are shown in
. The number of alleles per locus for
re-vealed high levels of polymorphism, ranging from 10 to 19 and from 10 to 17, respectively. Average Ho and He values ranged from 0.579 in Nakpoong to 0.765 in Wangpi, respectively, and no linkage disequilibrium was found in the Korean and Japanese populations. These results suggest that all of the microsatellite loci were polymorphic, with differences being detected in the number of alleles and observed heterozygosity in the examined Korean ayu populations. Four of the 12 Korean populations and one Japanese population (landlocked) showed significant devia-tion from the observed allele frequencies for HWE, suggesting that null alleles were present at some loci, as determined by MICRO-CHECKER.
estimates in ayu based on the microsatellites are given in
. Korean populations had high genetic distances when compared with landlocked (Biwa) populations. In addi-tion, high genetic distance was observed between landlocked (Biwa) and amphidromous (Kochi) populations in Japan. The observed results suggest there was low or restricted gene flow between the 13 amphidromous populations (12 Korean popula-tions and one Japanese population) and the landlocked ayu pop-ulation. This genetic differentiation between the two ecological forms of ayu in Japan has been detected in allozymes (Seki et al., 1985 ), mitochondrial DNA (Iguchi et al., 1999), and microsatel-lite markers (Takagi et al., 1999). Multi-locus pairwise estimates
Genetic variabilities at 3 loci of microsatellite DNA in Korean and Japanese ayuSize in base pair of alleles (S), total number of alleles (AT), Observed (HO) and expected (HE) heterozygosities, probability value estimates regarding deviation from Hardy-Weinberg equilibrium (P) and inbreeding coefficient (FIS).Departure from Hardy-Weinberg equilibrium: NS, not significant;*P< 0.05;**P< 0.01;***P< 0.001.
Genetic variabilities at 3 loci of microsatellite DNA in Korean and Japanese ayu Size in base pair of alleles (S), total number of alleles (AT), Observed (HO) and expected (HE) heterozygosities, probability value estimates regarding deviation from Hardy-Weinberg equilibrium (P) and inbreeding coefficient (FIS). Departure from Hardy-Weinberg equilibrium: NS, not significant; *P < 0.05; **P < 0.01; ***P < 0.001.
FSTvalues between samples (below diagonal) and probability of differentiation with P value in FSTestimate (above diagonal)+, significant; -, not significant inFST. Significance was tested at the 5%.
FST values between samples (below diagonal) and probability of differentiation with P value in FST estimate (above diagonal) +, significant; -, not significant in FST. Significance was tested at the 5%.
showed that differences between Korean and Japanese amphidromous populations were significant, although they had low
values (ranging from 0.019 to 0.048). These findings were also evident in the PCA scatter plots (
). The first two axes together explained 78% of the total genetic variation. The first (PC1) and second (PC2) axes explained 66% (
= 0.05) and 12% (
< 0.005) of the variance among the populations, respec-tively, and Japanese and Korean populations were separated. This implies a distinct geographical difference between Korea and Japan. Genetic differentiation between Korean and Japanese ayu populations was confirmed by Seki et al. (1988). Iguchi et al. (1999) found that although a phylogenetic tree with one Korean and six Japanese ayu populations showed few associations with geographic locations, there was a larger extent of net nucleotide substitutions between Korean and Japanese samples.
The lack of differentiation (
) and the PCA indicated (
) that there were no geographical trends among the popu-lations on the East Sea coast of Korea. Indeed, matrices of lin-earized genetic distance (
) and geographical distance (
) in kilometers between samples within the eastern Korean popula-tions were compared statistically to assess the effect of isola-tion by distance. The extent of isolation by distance could not be inferred from scatter plots of
. This analysis found no significant correlations with distance for populations in east-ern Korea (data not shown), suggesting there was a single large population. This is thought to represent the influence of a certain
Principal component analysis plotting the relationships between Korean (●) and Japanese (○) ayu populations.
Analysis of molecular variance (AMOVA) based on microsatellite DNA variation in amphidromous ayuThe percentage of variance (%), probability estimated from permutation (P), and theF-statistics (φ) are given at each hierarchical level. Numbers in paren-thesis correspond to populations sampled inFig. 1.
Analysis of molecular variance (AMOVA) based on microsatellite DNA variation in amphidromous ayu The percentage of variance (%), probability estimated from permutation (P), and the F-statistics (φ) are given at each hierarchical level. Numbers in paren-thesis correspond to populations sampled in Fig. 1.
level of gene flow through migration. The lack of regional equi-librium with isolation by distance within the Korean populations suggests that they may still remain in an unstable condition be-cause the equi-librium pattern with isolation by distance should require a sufficiently long period of time to achieve a stable condition (Hutchison and Templeton, 1999). However, pairwise population
estimates between the Jook population (South Sea coast) and all other populations were weak (
values for all pairs were lower than 0.05), but there was substantial dif-ferentiation in microsatellite variation with significant
values (< 0.05) compared to values for all other population pairs (
), suggesting that significant population subdivisions existed at small spatial scales.
The genetic structure of ayu populations in Korea and Japan was estimated by AMOVA (
). The variation within three groups (East Sea coast, South Sea coast, and Japan coast) was 2.10% (
< 0.05), suggesting the possibility of sub-structure among the populations. Analysis of variation within the three groups found no genetic variation within the groups (0.40%;
= 0.111); 97.5% of the total variation was due to differences within populations (
< 0.001). The hierarchical pattern of genetic differentiation among groups of ayu on the Korean and Japanese coasts indicates there were weak but historical patterns of isolation and restriction on gene flow. Pairwise population
estimates and AMOVA results were consistent with previous results based on isozyme data (Han et al., 2003), although few associations with geographic locations in amphidromous ayu populations were found in the present study. These results suggest that microsatellite loci can provide a powerful method for revealing genetic variation, with increased accuracy and resolution compared with isozyme markers.
The populations examined in this study were clustered us-ing the neighbor-joining (NJ) method (
). All of the Korean populations were separated from the two Japanese populations. The Jook River population was also separated from the cluster of the other Korean populations. Our findings support the pat-tern of genetic structure in wild ayu populations in Korea and Japan that has been revealed by genetic analyses (pairwise pop-ulation
estimates and AMOVA).
In conclusion, our results suggest that the ayu populations on the East Sea coast of Korea form a single population, and all of
Neighbor-joining tree of ayu populations. Nodal values for bootstrap support over 50% of the 1000 replicated trees.
the Korean populations are distinct from the Japanese popula-tions. The observed importance of genetic variation and genetic structure will provide a means for defining evolutionary and conservation units for the management and sustainable use of ayu resources. Further analyses of other genetic markers, such as maternally inherited mitochondrial DNA, and studies that in-clude more populations from other areas of Japan would help identify gene flow patterns in ayu.
This research was supported a grant from the NFRDI (RP-2011-BT-007).
Arlequin (version 3.0): an integrated software package for population genetics data analysis.
PHYLIP (Phylogeny Interference Package), Ver-sion 3.5c
Department of Genomic Sciences University of Washington
Seattle WA US.
Conservation and Evolution.
Cam-bridge University Press
New York NY US.
PCAGEN, Version 1.2.
Department of Ecology and Evolution University of Lausanne
FSTAT A Program to Estimate and Test Gene Diversities and Fixation Indices (Version 2.9.3).
Department of Ecology and Evolution University of Lausanne
Genetic variations of natural and hatchery populations of Korean ayu (Plecoglossus altivelis) by iso-zyme markers.
Correlation of pairwise ge-netic and geographic distance measures: inferring the relative in-fluences of gene flow and drift on the distribution of genetic vari-ability.
DOI : 10.2307/2640449
Genetic variation and geographic population structure of amphidromous ayu Plecoglossus altivelis as examined by mitochondrial DNA se-quencing.
Genetic variation and divergence in populations of ayu Plecoglossus altivelis including endangered subspecies inferred from PCR-RFLP analysis of the mitochondrial DNA D-loop region.
DOI : 10.1046/j.1444-2906.2002.00384.x
Substantial genetic differentiation in ayu Plecoglos-sus altivelis of the Japan and Ryukyu Islands.
Bull Jpn Sco Sci Fish
DOI : 10.2331/suisan.51.1269
Low lev-els of intraspecific variation in the mitochondrial DNA of chum salmon (Oncorhynchus keta).
Mol Mar Biol Biotechnol
GENEPOP (version 1.2): popula-tion genetics software for exact tests and ecumenicism.
Protection of intraspecific biodi-versity of exploited fishes.
Rev Fish Biol Fish
DOI : 10.1007/BF01103814
Distribution and genetic characteristics of the ayu on islands in the Tsushima Current area.
Jpn J Ichthyol
Genetic divergence among local popula-tions of ayu Plecoglossus altivelis in southwestern Japan.
Rep Usa Mar Boil Inst Kochi Univ
Genetic divergence among natural populations of ayu from Japan and Korea.
Nippon Suisan Gakkaishi
DOI : 10.2331/suisan.54.559
A measure of population subdivision based on micro-satellite allele frequencies.
Isolation and characterization of microsatellite loci from red sea bream Pagrus major and detection in closely related species.
Microsatellite DNA poly-morphism to reveal genetic divergence in ayu Plecoglossus altivelis.
Genetic variability and differ-entiation of amphidromous landlocked and hatchery populations of ayu Plecoglossus altivelis.
Bull Jpn Soc Sci Fish
DOI : 10.2331/suisan.49.1655
Van Oosterhout C
MICRO-CHECKER: software for identifying and correcting ge-notyping errors in microsatellite data.
Mol Ecol Notes
DOI : 10.1111/j.1471-8286.2004.00684.x
Estimating F-statistics for analysis of population structure.
DOI : 10.2307/2408641