Genetic Diversity of Thread-sail Filefish Stephanolepis cirrhifer Populations in Korean Coastal Waters Inferred from Mitochondrial DNA Sequence Analysis
- Author: Yoon Moongeun, Jung Ju-Yeon, Nam Yoon Kwon, Kim Dong Soo
- Organization: Yoon Moongeun; Jung Ju-Yeon; Nam Yoon Kwon; Kim Dong Soo
- Publish: Fisheries and aquatic sciences Volume 14, Issue1, p16~21, 31 March 2011
-
ABSTRACT
The genetic diversity and population genetic structure of thread-sail filefish,
Stephanolepis cirrhifer (Temminck & Schlegel), were examined with a nucleotide sequence analysis of a 495bp fragment of the 5′-end of the cytochromeb gene in 113 fish collected from five populations from the south and east coasts of the Korean Peninsula. Seventeen variable nucleotide sites and 16 haplotypes were defined. The observed haplotypes had a shallow haplotype genealogy and no geographical association. Most of the populations had high haplotype diversity and low nucleotide diversity, and significant negative values for Fu’s FS, suggesting rapid, recent population growth from an ancestral population and sudden population expansion. The estimated pairwise fixation indices (FST) indicate that substantial gene flow occurs among these populations.Thread-sail filefish in the South Sea of Korea and East Sea Korean populations forms a single panmictic population. Thus, thread-sail filefish in these areas should be treated as one management unit.
-
KEYWORD
Cytochrome b gene , Management unit , Population structure , Sea currents , Stephanolepis cirrhifer
-
Thread-sail filefish
Stephanolepis cirrhifer (Temminck & Schlegel) is an economically important fish species that is widely distributed in the western Pacific Ocean (Masuda et al., 1984; Youn, 2002). It is mostly captured in the coastal areas of the Korean Peninsula, south Japan, and China in the East China Sea, mostly in the Kuroshio Current from the tropical Pacific Ocean. In 1985, theS. cirrhifer catch exceeded 250,000 tonnes in the Northwest Pacific,but this decreased rapidly to less than 350 tonnes in 2002 (FAO, 2004). Stocks of this species have been decreasing as a result of overfishing and changes in its environment. This situation makes it urgent to evaluate the biological and genetic characteristics ofS. cirrhifer populations to maintain genetic variation,allowing the sustainable harvest of this species.Genetic variation is important for the long-term survival of natural populations because it confers the ability to adapt to environmental changes, thereby increasing the fitness of the population (Frankel and Soule 1981). Indeed, the structure of fish populations has attracted considerable interest because of its importance for the management of fisheries (Bailey,1997; Zhang et al., 2006). The estimation of population genetic structure using molecular markers is now commonly used to understand the dynamics of natural populations, especially resource recovery, and to monitor populations for fisheries management and conservation purposes. The population genetic structures of marine fishes are influenced by their larval behavior and dispersal patterns, which are determined by oceanographic features, such as sea currents, hydrological conditions, and physical barriers (Doyle et al., 1993; Umino et al., 2009; Hsieh et al., 2010; Kim et al., 2010). Most marine species have limited intraspecific population structures and high levels of gene flow influenced by sea currents(Palumbi, 1994; Sekino and Hara, 2001; Kochzius and Blohm, 2005; Song et al., 2010). However, little is known about the genetic structure of
S. cirrhifer in the Pacific Ocean.Nucleotide sequence information from maternally inherited mitochondria DNA (mtDNA) genes is useful for investigating inter- and intraspecific genetic variation among closely related taxa, species,or regional populations (Moritz et al., 1987: Avise,1994). The cytochrome
b gene (Cytb ), which contains both conserved and variable sites, is a useful marker for assessing the genetic structure of fish populations(Grant and Bowen, 1998; Hamner et al., 2007; Hwa et al., 2009).This study examined the genetic variation and population structure of
S. cirrhifer along the Korean coast using sequence analysis of the mtDNACytb gene to obtain basic data for fisheries resources management.Stephanolepis cirrhifer samples were collected from five sites located on two coastlines (the south (3 sites) and east (2 sites) coasts) of Korea in 2009 and 2010 (Table 1, Fig. 1). Collected samples were stored at -20°C or kept in 100% ethanol at room temperature until use. Genomic DNA was extracted using a conventional SDS/proteinase K method, followed by organic extraction and ethanol precipitation(Sambrook and Russell, 2001). Purified DNA was dried at room temperature, and then dissolved in TE buffer (10 mM Tris-HCl, 1 mM EDTA, pH 8.0).> Nucleotide sequence analysis
The polymerase chain reaction (PCR) was used to amplify the mtDNA
Cytb gene with the newly designed primers StepCytbF (5′-CCGGAACTAATG ACTTGAAA-3′) and StepCytbR (5′-AGCATTGTCT ACCGAGAAGCC-3′) based on the complete mitogenomic sequence of S. cirrhifer available in GenBank (NC_003177). PCR amplification was performed with a DNA Engine thermocycler (MJ Research, Tokyo, Japan) in 20-μL reactions containing 1-2 μg of genomic DNA, 2 μM of each primer, 0.25 mM of each dNTP, 1 unit of TaKaRa LA Taq DNA polymerase (TaKaRa Shuzo, Shiga, Japan),and 2 μL of 10× LA Taq reaction buffer (TaKaRa Shuzo). PCR conditions consisted of preheating at 94°C for 5 min, followed by 35 cycles of 94°C for 30 s, 55°C for 30 s, and 72°C for 30 s, with a final 5 min extension at 72°C. Size of the PCR product was verified by 1.5% agarose gel electrophoresis after ethidium bromide staining. The PCR product was purified with the AccuPrep PCR Purification Kit(Bioneer, Daejeon, Korea). After cycle sequencing with the ABI PRISM BigDye Terminator v3.1 Cycle Sequencing Ready Reaction Kit (Applied Biosystems,Foster City, CA, USA), the purified PCR product was sequenced directly on an ABI 3730xl DNA Analyzer(Applied Biosystems) with the same primers.
Sequence data were aligned using DnaSP (ver.4.90.1; Rozas and Rozas, 1997) to examine the haplotypes of the
Cytb gene. A parsimony network connecting the observed haplotypes was plotted with TCS (ver. 1.21; Clement et al., 2000) to resolve their genealogy. The haplotype diversity (h ) and nucleotide diversity (π) within the populations were estimated according to Nei (1987), based on Kimura’s twoparameter distance method, using K and DA in the program REAP (McElroy et al., 1993).Pairwise population FST values were calculated to estimate the genetic differentiation between populations, according to Slatkin and Hudson (1991), using the program Arlequin (ver. 3.1; Excoffier et al., 2005).The significance of each FST value was tested using 10,000 random permutations. Analysis of molecular variance (AMOVA) was used to test for population structure with Arlequin (ver. 3.1). Neutral expectation and historic demographic expansions were investigated by examining Fu’s FS and Tajima’s
D and mismatch distributions with the sudden expansion model (Rogers and Harpending 1992). Both the neutrality test and mismatch distribution analysis were performed in Arlequin (ver. 2000; Schneider et al., 2000).Sequence analysis of a 495-bp fragment of the 5′-end of the
Cytb gene detected 17 variable nucleotide sites among 113 individuals, defining a total of 16 haplotypes (Table 2). All of the substitutions at each site were biallelic, suggesting the occurrence of a single base substitution between sequences and no saturation of the substitutions, except at one site that was triallelic. A network diagram of the 16 observed haplotypes is shown in Fig. 2; the circle sizes reflect the haplotype abundances in the study areas. The focal haplotype, Hap1, was abundant, whereas the other haplotypes were rare and most likely derived from the focal haplotype, suggesting a shallow haplotype genealogy. The haplotype network had a star-like shape. The focal haplotype was common in all populations, while the genealogically related haplotypes were scattered throughout all populations(Table 3), suggesting that the observed haplotypes did not provide evidence of geographical association.The haplotype and nucleotide diversities are shown in Table 4. The average haplotype diversity was 0.517 and the average nucleotide diversity was 0.00153,indicating that the
S. cirrhifer populations examined have high haplotype and low nucleotide diversities.Thus, a shallow haplotype genealogy and the star-like shape pattern, and high level of haplotype diversity,but low nucleotide diversity within populations suggest the recent appearance and rapid population growth from an ancestral population with a small effective population size (Slatkin and Hudson, 1991;Rogers and Harpending, 1992; Avise, 2000).Fixation indices, based on mitochondrial haplotypes,showed no significant genetic differentiation between any population pair (Table 5), perhaps as a result of the high gene flow. Results of AMOVA were also not significant for the South Sea or East Sea Korean populations, suggesting that no population structure existed throughout the examined range of
S.cirrhifer . These results suggest thatS. cirrhifer in the South Sea and East Sea forms a single panmictic population. Indeed, large numbers shared the focal haplotype, suggesting that they arose from one large breeding population. Low genetic diversity in the Uljin population, from the northern-most East Sea, in comparison with the southern populations was observed in both the haplotype and nucleotide diversities. Indeed, this population has no genetic differentiation from southern populations, inferred from pairwise FST values (Table 5). Thus, theS.cirrhifer population in Uljin may be subject con-tinuously to high rates of gene flow from southern sources. The population structure of many marine
fishes often depends on the active migration of adults and passive dispersal of larvae (Waples, 1998;Cowen et al., 2000, 2003). The major current along the Korean coast is the Tsushima Warm Current(TWC), which branches off the Kuroshio Current in the East China Sea, with part of the TWC running into the Yellow Sea and the main part entering the East Sea along the Korean Peninsula (Senjyu, 1999;Ichikawa and Beardsley, 2002). Thus, our results suggest that the TWC transports
S. cirrhifer larvae in the East China Sea to the East Sea on the Korean coast.Neutrality indices calculated for each population and the entire sample are shown in Table 4. With both tests used (Fu’s FS and Tajima’s
D ), the populations had negative values and deviated significantly from the neutral evolution model, except one. Fu’s FS and Tajima’sD for the entire population were -15.07895(P<0.0001) and -2.12964 (P<0.0001), respectively.The mismatch distributions of all populations were unimodal (Fig. 3). All of our results suggest sudden population expansion.Our study suggests that the migration and reproductive behavior of
S. cirrhifer are similar to those of other marine fish species among which haplotypes with recently diverged nucleotides are distributed via the long-range dispersal of larvae among populations or regions (Sekino and Hara,2001; Kim et al., 2010). Our findings support the notion thatS. cirrhifer in the South Sea and East Sea Korean populations forms a single panmictic stock.Thus, theS. cirrhifer in the South Sea and East Sea Korean populations should be treated as one management unit. To have a comprehensive understanding of the population structure ofS. cirrhifer , it is necessary to obtain more mtDNA sequence information by expanding the sampling areas around Korea and elsewhere in the Pacific Ocean.-
[Table 1.] Sampling sites and dates geographical coordinates and number of examined individuals (N) in five Stephanolepis cirrhifer populations
-
[Fig. 1.] Sampling locations of the five Stephanolepis cirrhifer populations analyzed in this study
-
[Table 2.] Variable nucleotide sites in the mitochondrial DNA Cytb gene (495 bp) and the defined haplotypes of Stephanolepis cirrhifer
-
[Fig. 2.] A single minimum spanning tree of the 16 mtDNA Cytb haplotypes of Stephanolepis cirrhifer(Table 2). Circle sizes reflect haplotype abundances.
-
[Table 3.] Distribution of mtDNA Cytb haplotypes among five populations of Stephanolepis cirrhifer
-
[Table 4.] Measures of mtDNA diversity and neutrality indices calculated for each population of Stephanolepis cirrhifer
-
[Table 5.] Pairwise FST estimates among the five populations of Stephanolepis cirrhifer “?” is no significant differentiation (P>0.05). Probability of differentiation with P value was calculated from 1000 replications
-
[Table 6.] Results of the hierarchical analysis of molecular variance (AMOVA) based on mtDNA Cytb sequence data for Stephanolepis cirrhifer
-
[Fig. 3.] Mismatch distribution of Stephanolepis cirrhifer population in the south and east coast Korean peninsula. Bar graph: observed mismatch distribution; line graph: expected distribution for sudden expansion model.