The Study of Genetic Diversity and Population Structure of the Korean Fleshy Shrimp, Fenneropenaeus chinensis, Using Newly Developed Microsatellite Markers
The Study of Genetic Diversity and Population Structure of the Korean Fleshy Shrimp, Fenneropenaeus chinensis, Using Newly Developed Microsatellite Markers
Journal of Life Science. 2015. Dec, 25(12): 1347-1353
Copyright © 2015, Korean Society of Life Science
This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License ( which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
  • Received : September 21, 2015
  • Accepted : November 27, 2015
  • Published : December 30, 2015
Export by style
Cited by
About the Authors
은하 신
희정 공
보혜 남
영옥 김
봉석 김
동균 김
철민 안
형택 정
우진 김

The fleshy shrimp, Fenneropenaeus chinensis , is the family of Penaeidae and one of the most economically important marine culture species in Korea. However, its genetic characteristics have never been studied. In this study, a total of 240 wild F. chinensis individuals were collected from four locations as follows: Narodo (NRD, n = 60), Beopseongpo (BSP, n = 60), Chaesukpo (CSP, n = 60), and Cheonsuman (CSM, n = 60). Genetic variability and the relationships among four wild F. chinensis populations were analyzed using 13 newly developed microsatellite loci. Relatively high levels of genetic variability (mean allelic richness = 16.87; mean heterozygosity = 0.845) were found among localities. Among the 52 population loci, 13 showed significant deviation from the Hardy–Weinberg equilibrium. Neighbor-joining, principal coordinate, and molecular variance analyses revealed the presence of three subpopulations (NRD, CSM, BSP and CSP), which was consistent with clustering based on genetic distance. The mean observed heterozygosity values of the NRD, CSM, BSP, and CSP populations were 0.724, 0.821, 0.814, and 0.785 over all loci, respectively. These genetic variability and differentiation results of the four wild populations can be applied for future genetic improvement using selective breeding and to design suitable management guidelines for Korean F. chinensis culture.
The fleshy shrimp, Fenneropenaeus chinensis , is an economically important species in the family Penaeidae. F. chinensis is distributed from the west coast of the Korean Peninsula to the east coast of northern China [3] and is characterized by long-distance migrations of up to 2,000 km between spawning and feeding sites [18] . The F. chinensis South Korean aquaculture industry started in 1960, and production increased very quickly to reach 1,533 metric tons in 1997. However, a white spot syndrome virus (WSSV) disease outbreak decreased production dramatically to about 998 metric tons in 1998 [13] . Disease outbreaks cause mass mortality among cultured F. chinensis worldwide, particularly in Asian countries. To solve this problem, F. chinensis was replaced with Litopenaeus vannamei as the major cultured species in South Korea in 2003 [11] , but damage to shrimp farming associated with WSSV has been increasing in recent years. The West Sea Mariculture Research Center (Taean, Korea) conducts a breeding program to produce shrimp strains that are more resistant to WSSV. If a large number of cultured shrimp are released from aquaculture facilities, they could alter the genetic composition of wild populations by either displacing them or interbreeding. This could reduce the population’s ability to adapt to new environments. Accordingly, basic knowledge about the geographic distribution, genetic diversity, and population differences of F. chinensis is important. It is now widely recognized that this information can be obtained through recently developed molecular genetics techniques. However, few reports have been published about F. chinensis population genetics in South Korea.
Because microsatellite markers and simple sequence repeats have high levels of polymorphism, co-dominant inheritance, genome-wide distributions, and high reproducibility, they have been applied widely in population genetics, genetic linkage map, genetic diversity, and phylogenetic studies [10 , 16] . However, only a few microsatellite markers can be used for each shrimp species because of poor conservation of the microsatellite flanking sequences among different crustacean species, particularly in F. chinensis [15 , 27] . Several microsatellite markers have been isolated from F. chinensis [5 , 9 , 24] . Unfortunately, despite the commercial importance of this shrimp species in Korea, studies describing its genetic background are scarce. In the present study, we assessed genetic diversity within and among wild Korean F. chinensis populations and examined the genetic structure among these populations using new microsatellite DNA markers.
Materials and Methods
- Sample collection and DNA preparation
A total of 240 wild F. chinensis individuals were collected from four locations as follows: Narodo (NRD, n=60), Beopseongpo (BSP, n=60), Chaesukpo (CSP, n=60), and Cheonsuman (CSM, n=60) ( Fig. 1 and Table 1 ). All muscle tissue samples were stored in 100% ethanol prior to DNA extraction. The tissue was homogenized in lysis buffer (MFX-2000; Toyobo, Osaka, Japan) containing 20 mg/ml proteinase K. Total DNA was isolated using the MagExtractor MFX-6100 automated DNA extraction system (Toyobo Co., Tokyo, Japan). The extracted genomic DNA was quantified using a Nanodrop ND-1000 spectrophotometer (Thermo Fisher Scientific, Barrington, IL, USA) and stored at −20℃ until use.
PPT Slide
Lager Image
Sampling sites (◆) for the four wild Fenneropenaeus chinensis populations: Narodo (NRD); Beopseongpo (BSP); Chaesukpo (CSP); and Cheonsuman (CSM).
Geographic area, location, abbreviation, locality, collection date, and sample size collected
PPT Slide
Lager Image
NRD, Narodo; BSP, Beopseongpo; CSP, Chaesukpo; CSM, Cheonsuman
- Polymerase chain reaction (PCR) and microsatellite genotyping
Thirteen microsatellite loci that were amplified in F. chinensis were used in this study. Detailed information about the primers is presented in Table 2 . Each PCR contained three primer sets that were labeled at the 5‘ end of the forward primer with 6-FAM, HEX, or NED dyes (Applied Biosystems, Foster City, CA, USA). PCR amplification was carried out in a 10-μL reaction mixture containing 0.25 U Ex-Taq DNA polymerase (TaKaRa Biomedical Inc., Shiga, Japan), 1× PCR buffer, 0.2 mM dNTP mix, 10 pM of each primer, and 100 ng template DNA, using the PTC 200 DNA Engine (MJ Research, Waltham, MA, USA). The PCR conditions were: 5 min at 95℃, followed by 35 cycles of 30 sec at 95℃, 45 sec at 58℃, and 45 sec at 72℃, with a final extension of 10 min at 72°C. Microsatellite polymorphisms were screened using an ABI PRISM 3130 XL automated DNA sequencer (Applied Biosystems), and alleles were designated according to PCR product size relative to molecular size markers (GENESCAN 400 HD [ROX]; Applied Biosystems).
The 13Fenneropenaeus chinensismicrosatellite loci used in this study, with primer sequences, repeat motifs, and annealing temperatures (Ta) for polymerase chain reaction amplification
PPT Slide
Lager Image
The 13 Fenneropenaeus chinensis microsatellite loci used in this study, with primer sequences, repeat motifs, and annealing temperatures (Ta) for polymerase chain reaction amplification
- Statistical analysis
The number of samples ( n ), expected heterozygosity ( He ), and observed heterozygosity ( Ho ) were calculated using the Arlequin software package (ver. 3.0; [8] ). Tests for allelic richness ( AR ), number of alleles per locus ( Na ), and deviations from Hardy–Weinberg equilibrium (HWE) were estimated using GENEPOP ver. 4.0 (http://kimura.univmontp2. fr/~rousset/Genepop.htm), and the adjusted P-values for both analyses were obtained using a sequential Bonferroni test for multiple comparisons. MICRO-CHECKER 2.2.3 was used to test for the presence of null alleles. FST values (1,000 permutations; [25] ) were calculated using Arlequin 3.0. The population structure patterns were further investigated using a model-based Bayesian clustering procedure in STRUCTURE (ver. 2.3; [21] ), which assigns individuals to K populations based on their multilocus genotype. STRUCTURE was run for K = 2 using a burn-in length of 50,000 and a run of 50,000 steps. The overall inbreeding coefficient (F IS ; [25] ) was also estimated using GENEPOP ver. 4.0. Analysis of molecular variance (AMOVA; [7] ) was used to test for population structure with Arlequin 3.0. A principal coordinate analysis (PCoA) was performed using GENALEX (ver. 6.0; [19] ). A neighbor-joining (NJ) tree was constructed based on the chord distances (D CE ; [2] ) to evaluate the genetic population relationships using the POPULATION program (ver. 1.2.30). Bootstrap values were calculated using 1,000 replicates.
Results and Discussion
- Genetic variability
Microsatellite markers have been used in various shrimp species, and the pedigree of mixed populations can be determined using a few microsatellite markers [14] . However, little information is available on the use of microsatellite markers in F. chinensis [27] for understanding population structure or genetic diversity [15] . In this study, F. chinensis , a commercially and recreationally valuable species, were collected from the west coast of the Korean Peninsula (four separate collections in 2006) to illuminate its population structure and genetic diversity.
The genetic variability indices for the four populations are summarized in Table 3 . Allelic richness per locus ( AR ) ranged from 3 to 31 across all populations. Mean AR values were in the following order: NRD (16.2), CSM (17.5), BSP (16.8), and CSP (16.9). The total number of alleles per locus ( Na ) ranged from 3 to 31. KFcg959 had the highest AR and Na values in the NRD population. Mean Na ranged from 16.2 to 17.5. The mean expected and observed heterozygosities ranged from 0.835 to 0.861 and from 0.724 to 0.821, respectively. Mean heterozygosity was highest in CSM ( He = 0.861, Ho = 0.821), followed by BSP ( He = 0.843, Ho = 0.814), CSP ( He = 0.839, Ho = 0.785), and NRD ( He = 0.835, Ho = 0.724). The mean He value was higher than the mean Ho value in all populations. Among the 52 population-locus cases (4 populations × 13 loci), 13 cases showed significant deviations ( p <0.01). All populations departed from HWE at the KFcg614 locus. The NRD population departed from HWE at 6 of the 13 microsatellite loci, and the CSM, BSP, and CSP populations departed from HWE at three, two, and two microsatellite loci, respectively. The F IS value estimated for all populations at the KFcg614 locus was significantly different from zero ( p <0.05); all populations were in an excess heterozygosity condition for the KFc673 and KFc959 loci in NRD, the KFc42, KFc438, KFc657, and KFc959 loci in CSM; the KFc37, KFc42, KFc438, KFc568, and KFc658 loci in BSP; and the KFc37 and KFc658 loci in CSP.
Allelic variability at 13 microsatellite loci surveyed in the fourFenneropenaeus chinensispopulations
PPT Slide
Lager Image
Allelic richness (AR), number of alleles per locus (Na), expected heterozygosity (He), observed heterozygosity (Ho), inbreeding coefficient (FIS), and probability of significant deviation from Hardy–Weinberg equilibrium (P) are given for each population and locus. *P < 0.05; **P < 0.01. Calculations assume that individuals with one microsatellite band are homozygous for the allele. NRD, Narodo; BSP, Beopseongpo; CSP, Chaesukpo; CSM, Cheonsuman
Our results reveal that the newly developed microsatellite markers were a powerful approach to monitoring genetic diversity between the four geographically different wild Korean F. chinensis populations investigated. Heterozygosity is an important measure of population diversity at the genetic level. The mean observed heterozygosity values of the NRD, CSM, BSP, and CSP populations were 0.724, 0.821, 0.814, and 0.785 over all loci, respectively. These values were lower than the mean expected heterozygosity for the four populations. The overall heterozygosity in this study was high and differed from the reported mean heterozygosity, which is known to be low (7.3%), within crustacean populations as a whole [12] .
In the present study, 18 of the 52 population-locus cases deviated significantly from HWE after applying the Bonferroni correction. However, the KFc37, KFc41, and KFc42 loci in the NRD; the KFc614 locus in the CSM; and the KFc614 locus in the CSP population also deviated from HWE without excess heterozygosity. These deviations may have been caused by selection, population mixing, nonrandom mating, presence of null alleles, or the limited sample size used in our analysis [20 , 22] . The presence of null alleles, for example, is a classical source affecting the accuracy of microsatellite loci during parentage assignment [1] and null allele frequencies > 5% can compromise pedigree estimates [17] .
- Genetic differentiation between the four populations
FST is the proportion of total genetic diversity that separates groups, and values range from 0 to 1. If there is no population substructure (i.e., no stable groups), FST will approach 0. An FST range of 0.00-0.05 indicates little genetic differentiation [26] . Significant pairwise FST values ( p <0.05) were observed between the NRD and CSM, NRD and BSP, NRD and CSP, CSM and BSP, and CSM and CSP populations ( Table 4 ). In our study, pairwise FST tests detected low levels of genetic differentiation among the populations, particularly between the CSP and BSP populations, indicating that a geographical barrier was not effectively maintaining genetic integrity among the populations in the four areas.
PairwiseFSTestimates (below the diagonal) and Cavalli-Sforza and Edwards’s chord distance (DCE, above the diagonal) between the fourFenneropenaeus chinensispopulations
PPT Slide
Lager Image
Significance was tested at the 5% level with the Bonferroni correction applied for multiple tests. *p<0.05. Pairwise FST (Weir and Cockerham, 1984) and DCE distance (Cavalli-Sforza and Edwards, 1967) are measures of genetic differentiation and genetic distance between populations, respectively. NRD, Narodo; BSP, Beopseongpo; CSP, Chaesukpo; CSM, Cheonsuman
Genetic distances (D CE ) were also calculated for all possible population pairs. The D CE measure ranged from 0.0380 to 0.3582. The smallest estimate for D CE was between BSP and CSP (0.0380), whereas the highest estimate was between NRD and CSP (0.3582); the genetic distances between NRD and CSM, NRD and BSP, CSM and BSP, and CSM and CSP were 0.1864, 0.3030, 0.0991, and 0.1216, respectively ( Table 4 ). The result of the genetic distances among the F. chinensis populations was further confirmed by the findings of PCoA ( Fig. 2 ). The NJ tree constructed based on D CE , indicated that the four populations were allocated into three major groups ( Fig. 3 ); that is, one group included the CSM population, one group included the BSP and CSP populations, and one group included the NRD population. The PCoA produced a result similar to that of the cluster analysis ( Fig. 3 ) and showed clear separation of the four populations into three clusters.
PPT Slide
Lager Image
Principal coordinates analysis to determine the relationships among the wild Fenneropenaeus chinensis population based on chord distance (Cavalli-Sforza and Edwards, 1967)
PPT Slide
Lager Image
Neighbor-joining tree illustrating the relationships between the four Fenneropenaeus chinensis populations based on the Cavalli-Sforza and Edwards chord distance.
The results of hierarchical AMOVA tests to form putative groups estimated from each NJ tree topology provided additional evidence to support previous findings in F. chinensis populations. The NJ tree topology included a group for the NRD population and another group for the CSM, BSP, and CSP populations. In this case, the fixation index among groups ( FCT ) was not significant ( FCT =0.0147, P=0.2494), whereas it was significant among populations within groups ( FSC =0.0099, P=0.0000). The NJ tree topology defined a group that included the NRD, BSP, and CSP populations and another group that included the CSM population. In this case, the fixation index among groups ( FCT ) was not significant ( FCT =0.0010, P=0.4943), whereas it was significant among populations within groups (F SC =0.0167, P=0.0000). These results suggest that the populations remained structured within at least one group in each test case, as the FSC values estimated in both test cases were significant. Therefore, we separated them into three groups. The first group included the NRD population, the second included the CSM population, and the third included the BSP and CSP populations. These results show that the fixation index among groups ( FCT =0.0220, P=0.1643), and among populations within groups (F SC =−0.0012, P=0.8603), were not significant ( Table 5 ). The NJ tree topology and AMOVA results suggest that the four populations could be assigned to a group that included the NRD population, a group that included the CSM population, and a group that included the BSP and CSP populations ( Table 5 ). An evaluation of the evolutionary relationships among the four wild populations showed that BSP and CSP had the highest degree of genetic identity, followed by CSM and BSP, CSM and CSP, NRD and CSM, and NRD and BSP, whereas NRD and CSP were most distantly related. A previous study demonstrated the most likely explanation for this result is that wild populations of Chinese F. chinensis , which are extensively distributed in the Yellow and Bohai Seas, are two independent populations [15] . These two populations share the same wintering ground that lies in the mid-depth waters of the Yellow Sea [4] . An independent population of F. chinensis has been found off the west and south coasts of the Korean Peninsula [18] ; thus, further attempts are necessary to evaluate the genetic population relationships in this species. The goal of this study was to examine the importance of conservation and further genetic improvements in wild Korean F. chinensis . Information on genetic variation and differentiation in these four wild populations can be applied for future genetic improvement thorough selective breeding and to design suitable management guidelines for Korean F. chinensis .
Hierarchical analysis of molecular variance results for groups with several combinations of populations
PPT Slide
Lager Image
Probability values associated with F statistics are shown in parentheses. F statistic significantly greater than zero. * p<0.05. NRD, Narodo; BSP, Beopseongpo; CSP, Chaesukpo; CSM, Cheonsuman
The authors would like to acknowledge the help provided by Dr. Litticia Bryant at Queensland University of Technology in regard to constructive comments on the manuscript. This work was supported by a grant from the National Fisheries Research and Development Institute (R2015038).
Castro J. , Bouza C. , Presa P. , Pino-Querido A. , Riaza A. , Ferreiro I. , Sanchez L. , Martines P. 2004 Potential sources of error in parentage assessment of turbot (Scophthalmus maximus) using microsatellite loci Aquaculture 242 119 - 135    DOI : 10.1016/j.aquaculture.2004.09.010
Cavalli-Sforza L. L. , Edwards A. W. 1967 Phylogenetic analysis. Models and estimation procedures Am. J. Hum. Genet. 19 233 - 257
Deng J. Y. , Kang Y. D. , Jiang Y. W. 1983 A summary of surveys of the Penaeid shrimp spawning ground in the Bohai Bay Mar. Fish. Res. 5 17 - 23
Deng J. Y. , Ye C. C. , Liu Y. C. , Deng J. Y. 1990 its biology and management Ocean Press Beijing, China Penaeus chinensis in the Bohai and Yellow seas 48 - 91
Dong S. , Kong J. , Zhang T. , Meng X. , Wang R. 2006 Parentage determination of Chinese shrimp (Fenneropenaeus chinensis) based on microsatellite DNA markers Aquaculture 258 283 - 288    DOI : 10.1016/j.aquaculture.2006.04.044
Dong S. , Kong J. , Meng X. , Zhang Q. , Zhang T. , Wang R. 2008 Microsatellite DNA markers associated with resistance to WSSV in Penaeus (Fenneropenaeus) chinensis Aquaculture 282 138 - 141    DOI : 10.1016/j.aquaculture.2008.05.037
Excoffier L. , Smouse P. E. , Quattro J. M. 1992 Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data Genetics 131 479 - 491
Excoffier L. , Laval G. , Schneider S. 2005 Arlequin (version 3.0): an integrated software package for population genetics data analysis Evol. Bioinform 1 47 - 50
Gao H. , Kong J. , Yan B. , Yu F. , Luan S. , Cai S. 2008 Twelve new microsatellite markers for the Chinese shrimp Fenneropenaeus chinensis Mol. Ecol. Resour. 8 325 - 327    DOI : 10.1111/j.1471-8286.2007.01943.x
Goldstein D. G. , Pollock D. D. 1997 Launching microsatellites: a review of mutation processes and methods of phylogenetic inference J. Hered. 88 335 - 342    DOI : 10.1093/oxfordjournals.jhered.a023114
Jang I. K. , Jun J. C. , Jo G. J. , Cho Y. R. , Seo H. C. , Kim B. L. , Kim J. S. 2007 Polyculture of fleshy shrimp Fenneropenaeus chinensis and white shrimp Litopenaeus vannamei with river puffer Takifugu obscurus in shrimp ponds J. Aquacul. 20 278 - 288
Hedgecock D. , Tracey M. L. , Nelson K. , Abele L. G. 1982 the Biology of Crustacea Academic Press New York Genetics 284 - 403
Heo M. S. 2005 Characterization of the white spot syndrome baculovirus (WSBV) infection in fresh shrimp, Penaeus chinensis, cultured in Korea J. Life Sci. 15 248 - 252    DOI : 10.5352/JLS.2005.15.2.248
Herbinger C. M. , Doyle R. W. , Pitman E. R. , Paquet D. , Mesa K. T. , Morris D. B. , Wright J. M. , Cook D. 1995 DNA fingerprint based analysis of parental and maternal effects on offspring growth and survival in communally reared rainbow trout Aquaculture 137 245 - 256    DOI : 10.1016/0044-8486(95)01109-9
Liu P. , Meng X. H. , Kong J. , He Y. Y. , Wang Q. Y. 2006 Polymorphic analysis of microsatellite DNA in wild populations of Chinese shrimp (Fenneropenaeus chinensis) Aquac. Res. 37 556 - 562    DOI : 10.1111/j.1365-2109.2006.01461.x
Liu Z. J. , Cordes J. F. 2004 DNA marker technologies and their applications in aquaculture genetics Aquaculture 238 1 - 37    DOI : 10.1016/j.aquaculture.2004.05.027
Marshall T. C. , Slate J. , Kruuk L. E. B. , Pemberton J. M. 1998 Statistical confidence for likelihood-based pater- nity inference in natural populations Mol. Ecol. 7 639 - 655    DOI : 10.1046/j.1365-294x.1998.00374.x
Meng X. H. , Wang Q. Y. , Jang I. K. , Liu P. , Kong J. 2009 Genetic differentiation in seven geographic populations of the fleshy shrimp Penaeus (Fenneropenaeus) chinensis based on microsatellite DNA Aquaculture 287 46 - 51    DOI : 10.1016/j.aquaculture.2008.10.030
Peakall R. , Smouse P. E. 2006 GENALEX 6; genetic analysis in excel. Population genetic software for research and teaching Mol. Ecol. Notes 6 288 - 295    DOI : 10.1111/j.1471-8286.2005.01155.x
Pemberton J. M. , Slate J. , Bancroft D. R. , Barrett J. A. 1995 Non amplifying alleles at microsatellite loci: a caution for parentage and population studies Mol. Ecol. 4 249 - 252    DOI : 10.1111/j.1365-294X.1995.tb00214.x
Pritchard J. K. , Wen X. , Falush D. 2007 Documentation for structure software: version 2.2 University of Chicago Chicago 1 - 36
Rousset F. , Raymond M. 1995 Testing heterozygote excess and deficiency Genetics 140 1413 - 1419
Wang H. , Li F. , Xiang J. , Zhang C. , Yu K. 2008 Microsatellite-centromere distances and microsatellite diversity in different ploidy classes of Chinese shrimp (Fenneropenaeus Chinensis) Genetica 132 43 - 50
Wang X. , Guo X. , Zhang Y. , Meng X. , Qiu X. , Liu S. , Zhang T. 2009 Development of polymorphic EST-derived SSR markers for the shrimp, Fenneropenaeus chinensis Conserv. Genet. 10 1455 - 1457    DOI : 10.1007/s10592-008-9759-y
Weir B. S. , Cockerham C. C. 1984 Estimating F-statistics for the analysis of population structure Evolution 38 1358 - 1370    DOI : 10.2307/2408641
Wright S. 1998 Evolution and the genetics of population, variability within and among natural populations The University of Chicago Press Chicago
Xu P. , Zhou L. H. , Xiang J. H. 2001 Isolating microsatellite DNA of Chinese shrimp Penaeus chinensis Oceanol. Limnol. Sin. 32 255 - 259