The genetic diversity and population genetic structure of thread-sail filefish,
(Temminck & Schlegel), were examined with a nucleotide sequence analysis of a 495bp fragment of the 5′-end of the cytochrome
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 F
, suggesting rapid, recent population growth from an ancestral population and sudden population expansion. The estimated pairwise fixation indices (F
) 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.
(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, the
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 of
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
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
), 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
along the Korean coast using sequence analysis of the mtDNA
gene to obtain basic data for fisheries resources management.
Materials and Methods
- Samples and DNA extraction
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 (
). 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).
Sampling sites and dates geographical coordinates and number of examined individuals (N) in fiveStephanolepis cirrhiferpopulations
Sampling sites and dates geographical coordinates and number of examined individuals (N) in five Stephanolepis cirrhifer populations
- Nucleotide sequence analysis
The polymerase chain reaction (PCR) was used to amplify the mtDNA
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),
Sampling locations of the five Stephanolepis cirrhifer populations analyzed in this study
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.
- Data analysis
Sequence data were aligned using DnaSP (ver.4.90.1; Rozas and Rozas, 1997) to examine the haplotypes of the
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 (
) 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 F
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 F
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 F
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).
Results and Discussion
Sequence analysis of a 495-bp fragment of the 5′-end of the
gene detected 17 variable nucleotide sites among 113 individuals, defining a total of 16 haplotypes (
). 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
; 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
Variable nucleotide sites in the mitochondrial DNACytbgene (495 bp) and the defined haplotypes ofStephanolepis cirrhiferDots indicate that the nucleotide is identical to that in the Hap1 sequence.
Variable nucleotide sites in the mitochondrial DNA Cytb gene (495 bp) and the defined haplotypes of Stephanolepis cirrhifer Dots indicate that the nucleotide is identical to that in the Hap1 sequence.
), suggesting that the observed haplotypes did not provide evidence of geographical association.The haplotype and nucleotide diversities are shown in
. The average haplotype diversity was 0.517 and the average nucleotide diversity was 0.00153,indicating that the
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 (
), 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
. These results suggest that
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 F
). Thus, the
population in Uljin may be subject con-
A single minimum spanning tree of the 16 mtDNA Cytb haplotypes of Stephanolepis cirrhifer(Table 2). Circle sizes reflect haplotype abundances.
Distribution of mtDNACytbhaplotypes among five populations ofStephanolepis cirrhifer
Distribution of mtDNA Cytb haplotypes among five populations of Stephanolepis cirrhifer
tinuously to high rates of gene flow from southern sources. The population structure of many marine
Measures of mtDNA diversity and neutrality indices calculated for each population ofStephanolepis cirrhifer*Significant differentiation (P<0.05).South Sea, with part of the TWC running into the Korea.
Measures of mtDNA diversity and neutrality indices calculated for each population of Stephanolepis cirrhifer *Significant differentiation (P<0.05).South Sea, with part of the TWC running into the Korea.
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
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
. With both tests used (Fu’s F
), the populations had negative values and deviated significantly from the neutral evolution model, except one. Fu’s F
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 (
). All of our results suggest sudden population expansion.
Our study suggests that the migration and reproductive behavior of
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 that
in the South Sea and East Sea Korean populations forms a single panmictic stock.Thus, the
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 of
, it is necessary to obtain more mtDNA sequence information by expanding the sampling areas around Korea and elsewhere in the Pacific Ocean.
Pairwise FSTestimates among the five populations ofStephanolepis cirrhifer“?” is no significant differentiation (P>0.05). Probability of differentiation with P value was calculated from 1000 replications
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
Results of the hierarchical analysis of molecular variance (AMOVA) based on mtDNACytbsequence data forStephanolepis cirrhifer
Results of the hierarchical analysis of molecular variance (AMOVA) based on mtDNA Cytb sequence data for Stephanolepis cirrhifer
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.
This study was supported by a research grant(project no. #20088033-1) from the Ministry of Land,Transport, and Maritime Affairs, of the Republic of Korea.
Molecular Markers Natural History and Evolution
Chapman & Hall
Phylogeography: the History and Formation of Species
Harvard University Press
Structural dynamics and ecology of flatfish populations
Neth J Sea Res
TCS: a computer program to estimate gene genealogies
Connectivity of marine populations:open or closed?
The role of long distance dispersal in replenishing marine populations
Gulf Caribb Rec
Kendall Jr AW
A comparison of larval fish assemblages in the temperate zone of the northeast Pacific and northwest Atlantic Oceans
Bull Mar Sc
Arlequin ver.3.0: an integrated software package for population genetics data analysis
Depleted marine resources: an approach to quantification based on the FAO capture database.FAO Fisheries Circular
Rome Fisheries and Aquaculture Organization. DepartmentFishery Information Data and Statistics Unit
Conservation and Evolution
Cambridge University Press
New York NY
Shallow population histories in deep evolutionary lineages of marine fishes: insights from sardines and anchovies and lessons for conservation
Mitochondrial cytochrome b analysis reveals two invasive lionfish species with strong founder effects in the western Atlantic
J Fish Biol
Influence of hydrographic features on larval fish distribution during the south-westerly monsoon in the waters of Taiwan western North Pacific Ocean
J Fish Biol
Study of the cytochrome b gene sequence in populations of Taiwan
J Forensic Sci
The current system in the Yellow and East China Seas
Population structure of the olive flounder (Paralichthys olivaceus) in Korea inferred from microsatellite marker analysis
J Fish Biology
Genetic population structure of the lionfishPterois miles(Scorpaenidae Pteroinae) in the Gulf of Aqaba and northern Red Sea
Population structure and speciation in tropical seas: global phylogeography of the sea urchin Diadema
The fishes of the Japanese Archipelago
Tokai University Press
REAP: An integrated environment for the manipulation and phylogenetic analysis of restriction data
Evolution of animal mitochondrial DNA: relevance for population biology and systematic
Ann Rev Ecol Syst
Molecular Evolutionary Genetics
Columbia University Press
Genetic divergence reproductive isolation and marine speciation
Annu Rev Ecol Syst
Population growth makes waves in the distribution of pairwise genetic differences
Mol Biol Evol
DnaSP version 2.0: a novel software package for extensive molecular population genetics analysis
Comput Appl Biosci
Molecular Cloning A Laboratory Manual
Cold Spring Harbor Press
Arlequin: a software for popula-tion genetics data analysis. User manual ver 2.000
Genetics and Bi-ometry Lab Dept.of Anthropology University of Geneva
Application of microsatellite markers to population genetics studies of Japanese flounderParalichthys olivaceus
The Japan Sea intermediate water; its characteristics and circulation
Pairwise comparisons of mitochondrial DNA sequences in stable and exponentially growing populations montgomery
Population genetic structure and larval dispersal potential of spottedtail gobySynechogobius ommaturusin the north-west Pacific
J Fish Biol
Wild stock structure of Girella punctata in Japan revealed shallow genetic differentiation but subtle substructure in subsidiary distributions
Separating the wheat from the chaff:patterns of genetic differentiation in high gene flow species
Fishes of Korea with Pictorial Key and Systematic List
Population genetic structure of crimson snapperLutjanus erythropterusin East Asia revealed by analysis of the mitochondrial control region
ICES J Mar Sci