Shallow Population Genetic Structures of Thread-sail Filefish (Stephanolepis cirrhifer) Populations from Korean Coastal Waters

Article information

Asian-Australas J Anim Sci. 2012;25(2):170-176
Publication date (electronic) : 2012 February 01
doi : https://doi.org/10.5713/ajas.2011.11122
*Corresponding Author : D. S. Kim. Tel: +82-51-629-5914, Fax: +82-51-629-5908, E-mail: dongskim@pknu.ac.kr
Institute of Marine Living Modified Organisms, Pukyong National University, Busan 608-737, Korea
Received 2011 April 28; Accepted 2011 June 11.

Abstract

Genetic diversities, population genetic structures and demographic histories of the thread-sail filefish Stephanolepis cirrhifer were investigated by nucleotide sequencing of 336 base pairs of the mitochondrial DNA (mtDNA) control region in 111 individuals collected from six populations in Korean coastal waters. A total of 70 haplotypes were defined by 58 variable nucleotide sites. The neighbor-joining tree of the 70 haplotypes was shallow and did not provide evidence of geographical associations. Expansion of S. cirrhifer populations began approximate 51,000 to 102,000 years before present, correlating with the period of sea level rise since the late Pleistocene glacial maximum. High levels of haplotype diversities (0.974±0.029 to 1.000±0.076) and nucleotide diversities (0.014 to 0.019), and low levels of genetic differentiation among populations inferred from pairwise population FST values (−0.007 to 0.107), support an expansion of the S. cirrhifer population. Hierarchical analysis of molecular variance (AMOVA) revealed weak but significant genetic structures among three groups (FCT = 0.028, p<0.05), and no genetic variation within groups (0.53%; FSC = 0.005, p = 0.23). These results may help establish appropriate fishery management strategies for stocks of S. cirrhifer and related species.

INTRODUCTION

Filefish (family Monacanthidae) include 95 species widely distributed in both temperate and tropical seas (Nelson, 1994). The thread-sail filefish, Stephanolepis cirrhifer (Temminck and Schlegel), is an economically important species for the Korean fishery and aquaculture industries. S. cirrhifer aquaculture has received particular attention because of the fish’s high per weight market price and rapid growth rates (market size is reached in only one year) (Miyajima et al., 2011). Although, the species is also widely distributed in the western Pacific Ocean, occurring in the East Sea, the Yellow Sea, and the East China Sea (Masuda et al., 1984; Shao et al., 1990; Nelson, 1994; Ni and Kwok, 1999; Kim et al., 2005) most commercial catches are from the southern coastal area of the Korean Peninsula and the Kuroshio Current of the tropical Pacific Ocean.

In 1985, the S. cirrhifer catch in the Northwest Pacific exceeded 250,000 tonnes, but it had decreased to less than 950 tonnes by 2002 (FAO, 2004). Stocks of filefish species are decreasing as a result of overfishing and environmental degradation. Therefore, stock enhancement programs are critically needed to increase severely depleted fishery stocks. Hatchery-based release of filefish is one approach to revitalizing stocks and yields. However, hatchery stocks are susceptible to genetic degradation and reduction of genetic diversity: release of hatchery stocks can threaten the integrity of wild populations by changing the genetic composition and diversity of the gene pool (Allendorf and Ryman, 1987). Thus, evaluation of the biological and genetic characteristics of natural S. cirrhifer populations requires urgent attention to maintain genetic diversity levels and promote sustainable harvests of natural populations.

The primary aim of many historical and evolutionary studies of marine species is to ascertain oceanographic conditions, climatic factors, and vicariant events disrupting gene flow and genetic diversity, in relation to biogeographic patterns of demographic expansion and contraction. Genetic diversity is important for the long-term persistence and survival of natural populations because it confers the ability to adapt to changes in environmental conditions (Frankel and Soulé, 1981). The spatial distributions of population genetic structures are influenced by the effects of gene flow, natural selection, genetic drift and mutation (Slatkin, 1987), which may be affected by both current environmental conditions (Weersing and Toonen, 2009) and historical biogeographical processes (Avise, 2000). Assessments of the genetic structure of fish populations can provide invaluable information for successful conservation or management of exploited species (Bailey, 1997; Zang et al., 2006).

Information on nucleotide sequences in maternally inherited mitochondrial DNA (mtDNA) genes is useful for investigating interspecific and intraspecific genetic diversity among closely related taxa, species, or populations (Wilson et al., 1985; Moritz et al., 1987; Avise, 1994). The intraspecific population genetic structures of S. cirrhifer population are still poorly understood. A recent study by Yoon et al. (2011) found no differences in the nucleotide sequences of the mtDNA cytochrome b gene (Cyt b) between S. cirrhifer populations off the southern and eastern coasts of the Korean Peninsula. They suggest that substantial gene flow occurs within and between populations of S. cirrhifer, and that western Pacific S. cirrhifer belong to a single panmictic population; however, this study has not investigated the genetic structure of populations from the west coast of Korea.

In the mtDNA genome, the sequence variability of the non-coding control region is often higher than the variability in coding regions; thus, sequencing of the non-coding region has been recommended for investigating genetic diversity at the intraspecific level (Moritz et al., 1987). Up until now, the mtDNA control region has been generally used to study genetic diversities, population structures, and intraspecific phylogenesis of fish (Brown et al., 1993; Stepien and Faber, 1998; Sato et al., 2001; Guarnieo et al., 2002).

In the present study, we examine the genetic diversity and population structure of S. cirrhifer along the Korean coast using sequence analysis of the 5′ end of the mtDNA control region to obtain basic data for fisheries resources management.

MATERIALS AND METHODS

Sampling and sequencing

Stephanolepis cirrhifer specimens were collected in 2009 and 2010 from six sites located on three Korean seas: the Yellow Sea (one site), South Sea of Korea (South Sea) (three sites), and East Sea (two sites) (Table 1). All the samples examined, except the Yellow Sea, for the present mtDNA (control region gene) analysis had been used previously for the mtDNA analysis (Yoon et al., 2011). Specimens were stored at −80°C or preserved in ethanol at room temperature until DNA extraction. DNA was extracted from approximately 70 mg of each specimen with a QIAGEN Blood and Cell Culture DNA Midi Kit (Qiagen, Germany) according to the manufacturer’s instructions.

Sampling sites, dates, geographical coordinates, number of individuals examined, and measures of mtDNA diversity in six Stephanolepis cirrhifer populations

The polymerase chain reaction (PCR) was used to amplify the mtDNA control region with the primers StepCR-F (5′-CTAGCTCCCAAGCTAGGATT-3′) and StepCR-R (5′-TGGTGAGCCACGTATTGCAA-3′), newly designed on the basis of the complete mitogenomic sequence of S. cirrhifer, available in the GenBank genetic sequence data bank (NC_003177). PCR amplification was performed with a DNA Engine thermocycler (MJ Research, Tokyo, Japan) in 20-μl reaction tubes containing 100 ng 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). The PCR conditions consisted of preheating at 94°C for 5 min, followed by 35 cycles of 94°C for 30 s, 58°C for 30 s, and 72°C for 30 s, with a final 7 min extension at 72°C.

The size of the PCR product was verified by 1.0% 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). The primers used for sequencing were the same as those used for PCR amplification.

Data analysis

The sequence fragments obtained in this study were aligned with GENETIX-WIN ver. 4.0.1 (Software Development Co., Ltd, Tokyo, Japan) to identify sequence variants. The integrated software package DnaSP v. 4.90.1 (Rozas and Rozas, 1997) was used to determine haplotypes. Haplotype diversity (h) and nucleotide diversity (π) were estimated according to Nei (1987), based on Kimura’s two-parameter distance method, using the K and DA in the REAP program (McElroy et al., 1993). Genetic relationships among haplotypes were reconstructed using the neighbor-joining method, generated with the Seqboot, Neighbor, and Consensus options in PHYLIP v. 3.6 (Felsenstein, 1993). A bootstrap analysis of 1,000 replicates evaluated support for phylogenetic relationships after construction of a genetic distance matrix based on nucleotide divergences between haplotypes, estimated according to Nei (1987) and the Kimura two-parameter model (Kimura, 1980).

Pairwise population FST values were calculated to estimate genetic differentiation between populations (according to Slatkin, 1991) using ARLEQUIN v. 3.1 (Excoffier et al., 2005). The statistical significance of FST values was tested using 10,000 random permutations. The variance component was evaluated at each hierarchical level using analysis of molecular variance (AMOVA) (Excoffier et al., 2005) to assess genetic divergences within and among groups, using ARLEQUIN v. 3.1.

Neutral expectations and historic demographic expansions were investigated by examining Fu’s FS and Tajima’s D mismatch distributions with the sudden expansion model (Rogers and Harpending, 1992). A goodness-of-fit test was used to test the validity of the sudden expansion model using a parametric bootstrap approach based on the sum of squared deviations (SSD) to compare the observed and the estimated mismatch distributions (Schneider and Excoffier, 1999). Both the neutrality test and the mismatch distribution analysis were performed in ARLEQUIN v. 3.1 (Excoffier et al., 2005).

The mutation rate of the S. cirrhifer control region gene over the estimated time since expansion has not been determined. The evolutionary rate for the control region seems to vary among major taxonomic groups of marine fishes. For the Japanese sea bass (Lateolabrax japonicas) and spotted sea bass (Lateolabrax maculates), the sequence divergence rate of the control region is 6% per million years (Liu et al., 2006). In Nibea albiflora, the sequence divergence rate in the control region is 5% to 10% per million years (Han et al., 2008). In the present study, we calculated a sequence divergence rate of 5% to 10% per million years for the control region of S. cirrhifer.

RESULTS AND DISCUSSION

Genetic diversity

Analysis of a 336 base pair fragment of the 5′-end of the mtDNA control region in 104 individuals of S. cirrhifer from six populations defined a total of 70 haplotypes. The neighbor-joining tree of the 70 haplotypes was shallow and did not provide evidence for geographical associations (Figure 1). This suggests a signature of population expansion and/or high gene flow among populations.

Figure 1

Neighbor-joining tree for 70 control region haplotypes of Stephanolepis cirrhifer. Bootstrap supports of >50% in 1,000 replicates are shown.

Haplotype diversities (h) and nucleotide diversities (π) were high for all populations, ranging from 0.974±0.029 (JEJ) to 1.000±0.076 (WID), and 0.014 (JEJ, YES, POH) to 0.019 (ULJ), respectively (Table 1). Grant and Bowen (1998) suggest four basic scenarios for population demographic histories of marine fishes with different measures of haplotype and nucleotide diversities. They suggest that high haplotype and nucleotide diversities may be attributed to a long evolutionary history in a large population; examples include the Japanese anchovy (h = 0.91, π = 0.01) and the Atlantic bluefish (h = 0.70, π = 0.01). This condition might also be observed in re-mixed populations of individuals from historically separated populations (Avise, 2000). Therefore, genetic diversity values of S. cirrhifer populations might reflect strong dispersal capacities among ancestral populations with a long evolutionary history in Korean waters.

High genetic diversities within and between populations provide a potential genetic resource for future adaptation (Savolainen and Kuittinen, 2000; Hurt and Phillip, 2004). Recent research has examined genetic diversity in other fishes to determine the present status of populations in this area (e.g., Liu et al., 2006; Shui et al., 2009; Xiao et al., 2009).

Population genetic structure

Estimation of population structures using molecular markers is a powerful approach to understanding the dynamics of natural populations. Pairwise population FST estimates are generally low (Table 2), suggesting little genetic differentiation between population pairs, perhaps attributable to high gene flow. High gene flow is characteristic of marine fishes with passively dispersed planktonic larvae and adult migration abilities; their dispersals may be strongly affected by sea currents, and population substructures are limited (Han et al., 2008; Xiao et al., 2009; Kim et al., 2010). The Tsushima Warm Current (TWC) diverges from the Kuroshio Current, with the main flow entering the East Sea west of the Korean Peninsula and a subsidiary flow entering the Yellow Sea (Senjyu, 1999; Ichikawa and Beardsley, 2002). Thus, the TWC may transport larvae into the southern part of the East Sea, near Pohang and Uljin on the Korean coast, and to the Yellow Sea, near Wi-do. Previous studies based on mtDNA Cyt b nucleotide variation have suggested that the South Sea and East Sea S. cirrhifer populations comprise a single population (Yoon et al., 2011). Our mtDNA control region results also showed that estimates of the population differentiation seem to provide similar patterns of genetic differentiation among the South Sea and East Sea populations. However, pairwise population FST estimates between Wi-do and all other populations are relatively large (0.051 to 0.107), with significant p values (<0.05), compared to values for all other population pairs (−0.007 to 0.043) (Table 2), suggesting that significant population subdivisions exist at small spatial scales.

Pairwise FST (below diagonal) and p (above diagonal) values among populations of Stephanolepis cirrhifer

The genetic structure of S. cirrhifer populations around the Korean Peninsula was estimated by AMOVA (Table 3). Overall genetic variation among six populations (Analysis 1) was 2.34% (FCT = 0.02, p<0.05), suggesting the possibility of substructures among populations. Analysis of variation within three groups (Yellow Sea, South Sea, and East Sea; Analysis 2) indicated no genetic variation within groups (0.53%; FSC = 0.005, p = 0.23), 96.7% of the total variation due to differences within populations (FST = 0.033, p<0.05), and 2.77% of the total variation due to differences among groups (FCT = 0.028, p<0.05). The hierarchical pattern of genetic differentiation among groups of S. cirrhifer in Korean waters indicates weak but historical patterns of isolation and restriction of gene flow between groups. These results are consistent with previous results based on the mtDNA Cyt b sequence data (Yoon et al., 2011), suggesting that the sequence variability of the mtDNA control region of S. cirrhifer could be expected to provide a powerful means with an increased accuracy and resolution to reveal genetic variation as compared with the those of the Cyt b gene. However, differences between the three groups likely arose during isolation into refugia during glacial periods. Rare migrations of strays from the various refugia may have occurred during glaciations. Pleistocene glaciations (2.5 Ma to 10,000 years ago) influenced not only historical demographic patterns (glacial population extinctions and interglacial re-colonizations) but also the demographics of contemporary populations of phylogroups of marine species (Seeb and Crane, 1999; Kitano et al., 2007).

Results of the hierarchical analysis of molecular variance (AMOVA) based on mtDNA control region sequence data for Stephanolepis cirrhifer

At present, the East Sea is a semi-enclosed marginal sea, which was geographically isolated from the Pacific Ocean during the Pleistocene glacial period (Nishimura, 1974), by a large land bridge that extended from the Yellow Sea to Taiwan (Kimura, 1996, 2000). The present population structure of S. cirrhifer populations may represent the influence of this vicariant barrier. Populations of S. cirrhifer may have been isolated in the Yellow Sea, the South Sea, and the East Sea during times of glacial maxima.

Historic demography

Neutral expectations and historical demographic expansions of S. cirrhifer were investigated using Fu’s FS and Tajima’s D tests, and mismatch distributions with the sudden expansion model. Fu’s FS (−25.427; p = 0.000) and Tajima’s D (−1.708; p = 0.020) showed significant deviations from the neutral evolution model. The mismatch distributions of S. cirrhifer in Korean coastal waters appeared to be unimodal (Figure 2). Significant differences for the sums of the square deviations (pSSD<0.05) between the observed and simulated mismatch distributions indicate that the population is at equilibrium (i.e., a non-expansion phase). The goodness-of-fit test (pSSD = 0.657) did not reject the null hypothesis of sudden population expansion for S. cirrhifer populations. Thus, both the neutrality tests and the mismatch distributions indicate a recent expansion of S. cirrhifer populations.

Figure 2

Mismatch distribution constructed using pairwise differences among the mtDNA haplotypes of Stephanolepis cirrhifer. The bars are observed pairwise differences values and the open circles and solid line are the expected mismatch distributions under the sudden expansion model.

Using τ (3.436), the expansion was estimated to have occurred 51 to 102 Ma before present, which is consistent with sea level rise since the late Pleistocene glacial maximum. Rising sea level removed barriers to migration and dispersal during the period of postglacial warming, establishing potential gene flow between previously isolated populations. High levels of genetic diversity and low levels of genetic differentiation determined in the present study also support population expansion since the late Pleistocene.

In conclusion, this study found high levels of gene flow among populations of S. cirrhifer, weak but significant genetic differentiation of populations on the west coast of Korea (Yellow Sea) and pooled populations in southern and eastern Korean coastal waters. Incorporation of additional population in the Yellow Sea thus provided an increased resolution in the geographic differentiation among three groups of populations. Moreover, the genetic structures of S. cirrhifer populations indicate the presence of three geographic population groups (Yellow Sea, South Sea, and East Sea), each of which may require separate fisheries management strategies.

ACKNOWLEDGMENTS

This work was supported by a research grant (project no. 20088033-1) from the Ministry of Land, Transport, and Maritime Affairs, of the Republic of Korea.

References

Allendorf FW, Ryman N. 1987. Population genetics and fishery management. Genetic Management of Hatchery Stocks In : Ryman N, Utter F, eds. University of Washington Press. Seattle, WA, USA: p. 141–159.
Avise JC. 1994. Molecular markers, natural history and evolution Chapman & Hall. New York, USA:
Avise JC. 2000. Phylogeography: The History and Formation of Species Harvard University Press. Cambridge, MA, USA:
Bailey KM. 1997;Structural dynamics and ecology of flatfish populations. Neth J Sea Res 37:269–280.
Brown JR, Beckenbach AT, Smith MJ. 1993;Intraspecific DNA sequence variation of the mitochondrial control region of white sturgeon (Acipenser transmontanus). Mol Biol Evol 10:326–341.
Excoffier L, Laval G, Schneider S. 2005;Arlequin ver. 3.0: an integrated software package for population genetics data analysis. Evol Bioinform 1:47–50.
FAO. 2004. Depleted marine resources: an approach to quantification based on the FAO capture database. FAO Fisheries Circular, No. 1011 p. 18–19. Rome: Fisheries and Aquaculture Organization. Department, Fishery Information, Data and Statistics Unit.
Felsenstein J. 1993. PHYLIP (Phylogeny Interference Package) Version 3.5c Department of Genetics, Washington State University. Seattle, USA:
Frankel OH, Soulé ME. 1981. Conservation and evolution Cambridge University Press. New York, USA:
Grant WS, Bowen BW. 1998;Shallow population histories in deep evolutionary lineages of marine fishes: insights from sardines and anchovies and lessons for conservation. J Hered 89:415–426.
Guarnieo I, Franzellitti S, Ungaro N, Tommasini S, Piccinetti C, Tinti F. 2002;Control region haplotype variation in the central Mediterranean common sole indicates geographical isolation and population structuring in Italian stock. J Fish Biol 60:1459–1474.
Han ZQ, Gao TX, Yanagimoto T, Sakurai Y. 2008;Genetic population structure of Nibea albiflora in Yellow Sea and East China Sea. Fish Sci 74:544–552.
Hurt C, Phillip H. 2004;Conservation genetics in aquatic species: general approaches and case studies in fishes and spring snails of arid lands. Aquat Sci 66:402–413.
Ichikawa H, Beardsley RC. 2002;The current system in the Yellow and East China Seas. J Oceanog 58:77–92.
Kim IS, Choi Y, Lee CL, Lee YJ, Kim BJ, Kim JH. 2005. Illustrated book of Korean fishes Kyo-Hak Pub Co. Seoul: (in Korean).
Kim WJ, Kim KK, Han HS, Nam BH, Kim YO, Kong HJ, Noh JK, Yoon M. 2010;Population structure of the olive flounder (Paralichthys olivaceus) in Korea inferred from microsatellite marker analysis. J Fish Biol 76:1958–1971.
Kimura M. 1980;A simple method for estimating evolutionary rates of base substitutions that compare studies of nucleotide sequences. J Mol Evol 16:111–120.
Kimura M. 1996;Quaternary paleogeography of the Ryukyu. Arc J Geol 105:259–285.
Kimura M. 2000;Paleography of the Ryukyu Islands. Tropics 10:5–24.
Kitano J, Mori S, Peichel CL. 2007;Phenotypic divergence and reproductive isolation between sympatric forms of Japanese threespine sticklebacks. Biol J Linn Soc 91:671–685.
Liu JX, Gao TX, Yokogawa K, Zhang YP. 2006;Differential population structuring and demographic history of two closely related fish species, Japanese sea bass (Lateolabrax japonicus) and spotted sea bass (Lateolabrax maculatus) in Northwestern Pacific. Mol Phylogenet Evol 39:799–811.
Masuda H, Amaoka K, Araga C, Uyeno T, Yoshino T. 1984. The fishes of the Japanese Archipelago Tokai University Press. Tokyo, Japan:
McElroy D, Moran P, Bermingham E, Kornfield I. 1993;REAP: an integrated environment for the manipulation and phylogenetic analysis of restriction data. J Hered 83:157–158.
Miyajima Y, Masuda R, Kurihara A, Kamata R, Yamashita Y, Takeuchi T. 2011;Juveniles of thread-sail filefish, Stephanolepis cirrhifer, can survive and grow by feeding on moon jellyfish Aurelia aurita. Fish Sci 77:41–48.
Moritz C, Dowling TW, Brown WM. 1987;Evolution of animal mitochondrial DNA: relevance for population biology and systematic. Ann Rev Ecol Syst 18:269–292.
Nei M. 1987. Molecular evolutionary genetics Columbia University Press. New York, USA:
Nelson JS. 1994. Fishes of the World 3rd Ednth ed. John Wiley & Sons Inc. New York, USA:
Ni IH, Kwok KY. 1999;Marine fish fauna in Hong Kong waters. Zool Stud 38:130–152.
Nishimura S. 1974. Formation of the sea of Japan Tsukiji Shokan. Tokyo, Japan:
Rogers AR, Harpending H. 1992;Population growth makes waves in the distribution of pairwise genetic differences. Mol Biol Evol 9:552–569.
Rozas J, Rozas R. 1997;DnaSP version 2.0: a novel software package for extensive molecular population genetics analysis. Comput Appl Biosci 13:307–311.
Sato S, Ando J, Ando H, Urawa S, Urano A, Abe S. 2001;Genetic variation among Japanese populations of chum salmon inferred from the nucleotide sequences of the mitochondrial DNA control region. Zoolog Sci 18:99–106.
Savolainen O, Kuittinen H. 2000. Small population processes. Forest Conservation Genetics: Principles and Practice In : Young A, Boshier D, Boyle T, eds. CABI Publishing. Wallingford, UK: p. 91–100.
Seeb LW, Crane P. 1999;Allozymes and mitochondrial DNA discriminate Asian and North American populations of chum salmon in mixed-stock fisheries along the south coast of the Alaska Peninsula. Trans Am Fish Soc 128:88–103.
Senjyu T. 1999;The Japan Sea intermediate water; its characteristics and circulation. J Oceanogr 55:111–122.
Schneider S, Excoffier L. 1999;Estimation of past demographic parameters from the distribution of pairwise differences when mutation rates vary among sites: application to human mitochondrial DNA. Genetics 152:1079–1089.
Shao KT, Lin CP, Ho LT, Lin PL. 1990;Study on the fish communities from northern and southern waters of Taiwan by analyzing the impingement data. J. Fish. Soc Taiwan 17:73–90.
Slatkin M, Hudson RR. 1991;Pairwise comparisons of mitochondrial DNA sequences in stable and exponentially growing populations. Genetics 129:555–562.
Stepien CA, Faber JE. 1998;Population genetic structure, phylogeography, and spawning philopatry in walleye (Stizostedion vitreum) from mtDNA control region sequences. Mol Ecol 7:1757–1769.
Shui BN, Han ZQ, Gao TX, Miao ZQ, Yanagimoto T. 2009;Mitochondrial DNA variation in the East China Sea and Yellow Sea populations of Japanese Spanish mackerel Scomberomorus niphonius. Fish Sci 75:593–600.
Weersing K, Toonen RJ. 2009;Population genetics, larval dispersal, and connectivity in marine systems. Mar Ecol Prog Ser 393:1–12.
Wilson AC, Cann RL, Carr SM, George M, Gyllensten UB, Helm-Btchowski K, Higuchi RG, Palumbi SR, Prager EM, Sage RD, Stone-King M. 1985;Mitochondrial DNA and two perspectives on evolutionary genetics. Biol J Linn Soc 26:375–400.
Xiao Y, Zhang Y, Gao T, Yanagimoto T, Yabe M, Sakurai Y. 2009;Genetic diversity in the mtDNA control region and population structure in the small yellow croaker Larimichthys polyactis. Environ Biol Fishes 85:303–314.
Yoon M, Jung JY, Nam YK, Kim DS. 2011;Genetic diversity of thread-sail filefish (Stephanolepis cirrhifer) populations in Korean coastal waters inferred from mitochondrial DNA sequence analysis. Fish Aqua Sci 14:16–21.
Zhang J, Cai Z, Huang L. 2006;Population genetic structure of crimson snapper Lutjanus erythropterus in East Asia, revealed by analysis of the mitochondrial control region. ICES J Mar Sci 63:693–704.

Article information Continued

Figure 1

Neighbor-joining tree for 70 control region haplotypes of Stephanolepis cirrhifer. Bootstrap supports of >50% in 1,000 replicates are shown.

Figure 2

Mismatch distribution constructed using pairwise differences among the mtDNA haplotypes of Stephanolepis cirrhifer. The bars are observed pairwise differences values and the open circles and solid line are the expected mismatch distributions under the sudden expansion model.

Table 1

Sampling sites, dates, geographical coordinates, number of individuals examined, and measures of mtDNA diversity in six Stephanolepis cirrhifer populations

Sampling site Abbreviation Date of collection Geographical co-ordinates N No. of haplotypes h π

Latitude Longitude
Yellow Sea
 Wi-do WID Oct. 2010 35°35′33.14″N 126°16′23.65″E 7 7 1.000±0.076 0.018
South Sea
 Jeju JEJ Aug. 2010 33°03′14.54″N 126°32′39.67″E 18 15 0.974±0.029 0.014
 Yeosu YES Aug. 2010 34°41′22.99″N 127°39′59.62″E 23 21 0.992±0.015 0.014
 Geoje GUJ July. 2010 34°43′15.11″N 128°41′07.46″E 24 18 0.967±0.024 0.015
East Sea
 Pohang POH Oct. 2010 36°02′50.30″N 129°24′31.75″E 18 15 0.980±0.024 0.014
 Uljin ULJ Oct. 2009 36°58′04.20″N 129°26′09.39″E 21 18 0.986±0.019 0.019

N: Sample size, h: Haplotype diversity, π: Nucleotide diversity.

Table 2

Pairwise FST (below diagonal) and p (above diagonal) values among populations of Stephanolepis cirrhifer

Population WID JEJ YES GUJ POH ULJ
WID 0.05 0.01 0.02 0.01 0.01
JEJ 0.051 0.78 0.13 0.02 0.08
YES 0.074 −0.014 0.15 0.08 0.35
GUJ 0.103 0.020 0.016 0.23 0.56
POH 0.107 0.043 0.027 0.011 0.38
ULJ 0.089 0.025 0.004 −0.007 0.002

Table 3

Results of the hierarchical analysis of molecular variance (AMOVA) based on mtDNA control region sequence data for Stephanolepis cirrhifer

Source of variation df Sum of squares Variance components Percentage of variation F index p
Analysis 1. One gene pool
 Among populations 5 18.321 0.06127 Va 2.34 FCT = 0.02 <0.05
 Within populations 105 267.931 2.55172 Vb 97.66
Analysis 2. Three gene pools (Yellow Sea, South Sea, East Sea)
 Among groups 2 9.797 0.07311 Va 2.77 FCT = 0.028 <0.05
 Within groups 3 8.524 0.01393 Vb 0.53 FSC = 0.005 0.23
 Within populations 105 267.931 2.55172 Vc 96.70 FST = 0.033 <0.05