Next Article in Journal
Enabling Anticancer Therapeutics by Nanoparticle Carriers: The Delivery of Paclitaxel
Previous Article in Journal
Novel Recovery of Nano-Structured Ceria (CeO2) from Ce(III)-Benzoxazine Dimer Complexes via Thermal Decomposition
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genetic Diversity and Differentiation of the Orange-Spotted Grouper (Epinephelus coioides) Between and Within Cultured Stocks and Wild Populations Inferred from Microsatellite DNA Analysis

State Key Laboratory of Biocontrol, Institute of Aquatic Economic Animals and the Guangdong Province Key Laboratory for Aquatic Economic Animals, School of Life Sciences, Sun Yat-Sen University, Guangzhou 510275, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2011, 12(7), 4378-4394; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms12074378
Submission received: 26 May 2011 / Revised: 17 June 2011 / Accepted: 28 June 2011 / Published: 6 July 2011
(This article belongs to the Section Biochemistry)

Abstract

:
In the present study, we employed microsatellite DNA markers to analyze the genetic diversity and differentiation between and within cultured stocks and wild populations of the orange-spotted grouper originating from the South China Sea and Southeast Asia. Compared to wild populations, genetic changes including reduced genetic diversity and significant differentiation have taken place in cultured grouper stocks, as shown by allele richness and heterozygosity studies, pairwise Fst, structure, molecular variance analysis, as well as multidimensional scaling analysis. Although two geographically adjacent orange-spotted grouper populations in China showed negligible genetic divergence, significant population differentiation was observed in wild grouper populations distributed in a wide geographical area from China, through Malaysia to Indonesia. However, the Mantel test rejected the isolation-by-distance model of genetic structure, which indicated the genetic differentiation among the populations could result from the co-effects of various factors, such as historical dispersal, local environment, ocean currents, river flows and island blocks. Our results demonstrated that microsatellite markers could be suitable not only for genetic monitoring cultured stocks but also for revealing the population structuring of wild orange-spotted grouper populations. Meanwhile, our study provided important information for breeding programs, management of cultured stocks and conservation of wild populations of the orange-spotted grouper.

1. Introduction

The orange-spotted grouper (Epinephelus coioides), an economically high valued marine food fish species belongs to the subfamily Epinephelinae (family, Serranidae) [1], which inhabits in a large area from eastern Africa, south to at least Durban (South Africa), east to the western Pacific from the Ryukyu Islands to Australia, and eastward to Palau and Fiji [2]. However, due to overfishing and habitat destruction, the wild populations of the orange-spotted grouper have declined greatly in recent years and this species has been classified as nearly threatened [3]. Over the last decades, various attempts have been made to conserve this grouper species. The most important achievement is successful practices in aquaculture for the orange-spotted grouper, which can alleviate the fishing pressure on wild populations [4]. According to the FAO fishery statistics, the global aquaculture production of the orange-spotted grouper has dramatically increased nearly 40-fold between 1999 and 2008 [5]. At present, it has become one of the grouper species most commonly cultured on commercial scales in Asian-Pacific region and a major food fish in Hong Kong live fish markets in China [6].
However, aquaculture practices are likely to reduce the genetic diversity by founder effects, inbreeding and random genetic drift in hatchery-reared stocks, which may cause loss of disease resistance and reduce environmental adaptability, thus limit the genetic potential for selective breeding [79]. For this reason, some research reports focused on genetic monitoring of the differences between hatchery stocks and wild populations in various fish species including Turbot (Scophthalmus maximus), Atlantic salmon (Salmo salar), Sea bream (Pagrus major), Japanese flounder (Paralichthys olivaceus), Arctic charr (Salvelinus alpinus), common carp (Cyprinus carpio) using DNA fingerprinting techniques such as microsatellite and mitochondrial markers [1015]. Unfortunately, despite the importance of orange-spotted grouper for commercial aquaculture in Southeast Asia, literature about its genetic background is scarce and mainly focuses on the genetic differentiation among geographical populations [16,17], therefore the genetic difference between cultured stocks and wild populations is still unknown.
Marine fishes are considered to be genetically less differentiated compared to freshwater fishes due to high gene flow in open marine environment [18]. The dispersal potential of planktonic eggs and larvae and the absence of geographical barriers between different populations all contribute low genetic differentiation in marine fishes. Nevertheless, the study on population genetic structure of marine fishes is crucial and can provide important information for conservation and management of the genetically different populations. Meanwhile, such study can also suggest essential theoretical and practical guidance for improved breeding programs, as the genetically differentiated populations have the potential to be used for hybridization to improve the genetic quality of cultured fish [19,20].
We have recently isolated polymorphic microsatellite markers and shown that these markers have the potential to characterize the genetic diversity and analyze the population structure of the orange-spotted grouper [21]. In the present study, we employed microsatellite DNA markers to compare the genetic diversity and differentiation between and within cultured stocks and wild populations of the orange-spotted grouper originating from the South China Sea and Southeast Asia including Malaysia and Indonesia. Our goals were to evaluate the putative genetic changes in cultured stocks as a result of founder effect, random genetic drift and inbreeding during aquaculture practices compared to wild populations, and examine the population differentiation within samples distributed in a wide area from China, through Malaysia to Indonesia in order to facilitate the breeding programs of the orange-spotted grouper.

2. Materials and Methods

2.1. Sample Collection and DNA Extraction

Our samples were mainly obtained from the South China Sea, from which almost all the cultured and wild orange-spotted grouper in China mainland presently originated, due to its ideal subtropical climate. Sampling locations are shown in Figure 1. In this study, all samples from China, including three cultured stocks (HZ, SZ and ZJ) and two wild populations (DYB and HNI) were acquired within this region. In detail, HZ stock (n = 26) was the first generation progeny of a broodstock reared in a government-sponsored hatchery (Guangdong Daya Bay Fisheries Development Center, Huizhou City, China). Both SZ (n = 15) and ZJ (n = 34) stocks were also the first generation offspring from private hatcheries located in Shenzhen and Zhanjiang Cities, respectively. These hatcheries were small- or medium-sized with a limited number of orange-spotted grouper. Both wild caught and captive cultured sources of the broodstock were investigated. In detail, large-sized groupers in the wild were collected during the spawning season and treated as broodstock directly. Alternatively, captured juveniles were cultured in captivity until sexual maturation and then used as broodstock. However, information about the exact proportion of these two different sources of individuals in each broodstock and management practices on cultured stocks were not available. Two wild populations of China, DYB (n = 50) and HNI (n = 40), were caught in the South China Sea near Daya Bay and Hainan Island, respectively. Furthermore, another two wild populations collected in the Southeast Asia were used as ‘outgroup’ for comparison with the Chinese samples, including the Sandakan population of Malaysia (SDK, n = 25) and Tarakan population of Indonesia (TRK, n = 11). All samples of the wild orange-spotted grouper were caught by hook and line in 2009.
Caudal fin or muscle tissue samples of the orange-spotted grouper were clipped and immediately preserved in 90% ethanol until DNA extraction was performed. Genomic DNA was isolated according to the phenol-chloroform extraction protocol by Hoelzel and Green [22].

2.2. Microsatellite Genotyping

A total of 11 microsatellite loci were selected for genotyping fish individuals from the seven samples (Table 1). Forward primers were 5′-labeled with a fluorescent dye HEX or 6-FAM. PCR amplification was performed on PTC-200 thermal cycler (MJ Research) in a 20-μL volume containing 0.5 μm of each primer, 0.2 mM dNTP, 1.5 mM MgCl2, 1XPCR buffer, 1U Taq DNA polymerase (Fermentas), and 20ng template genomic DNA. PCR program was as follows: 5 min at 94 °C followed by 30 cycles of 30 s at 94 °C, 30 s at the annealing temperature (Table 1), and 30 s at 72 °C with a final extension of 5 min at 72 °C. PCR products were separated on an ABI PRISM 3730 DNA automated sequencer (Applied Biosystems). Fragment size was measured according to the ROX-500 standard using GeneMapper (Applied Biosystems). Genotype data was exported to Excel tables for analysis. MicroChecker [28] was used to check the genotyping errors and possible presence of null alleles.

2.3. Statistical Analysis

Intra-population genetic diversity was assessed by computing the number of alleles per locus (A), Allelic richness (Ar), observed heterozygosity (Ho) and expected heterozygosity (He) using the program FSTAT 2.9.3 [29]. Since allele number is influenced by size of different samples, we used allele richness for comparison in this study [30]. Differences in genetic diversity parameters between cultured stocks and wild populations were tested using nonparametric analysis (Mann-Whitney U test; [31]). The software GENEPOP 4.0 computer package [32] was used to test linkage disequilibrium between pairwise loci and departure from Hardy-weinberg equilibrium (HWE) across all loci. The Markov chain method was employed to calculate an unbiased estimate of the p-value in order to test deficiency or excess of heterozygote with the following parameters (dememorization = 1000, batches = 500, and iterations per batches = 1000). Sequential Bonferroni correction was used to adjust the significance of HWE and linkage disequilibrium tests [33]. The program BOTTLENECK version 1.2.02 [34] was used under the infinite allele model (IAM), stepwise-mutation model (SMM) and two-phase model (TPM) with 1000 iterations to check the existence of bottleneck inferred by heterozygosity excess in the seven samples. The significance was tested using Wilcoxon sign rank test. Following the recommendations of Kuikart & Cornuet [34], if significant results were obtained under all these three models, populations can be concluded as having experienced bottleneck. Inbreeding coefficients (ƒ) were also calculated using the program GDA [35].
Inter-population genetic differentiation was estimated with pairwise Fst values and significance tests of pairwise Fst was computed by a permutation with 10,000 replicates using ARLEQUIN 3.11 [36]. Furthermore, the genetic relationship among samples was visualized by MDS (Multidimensional scaling) based on pairwise Fst values using SPSS 13.0. Analysis of molecular variance (AMOVA) was also calculated using ARLEQUIN 3.11 [36]. AMOVA was used to partition genetic variance hierarchically between the wild group and cultured group. However, AMOVA was not used to partition wild groups in different geographical regions (China, Malaysia and Indonesia). As two of these three groups consisted of a single population, this would influence the analysis in a discriminatory way. The correlation between genetic distance and geographical distance was evaluated using the Mantel test by the program TFPGA version 1.3 [37]. We also analyzed the population structure of all samples using the program Structure 2.2.3 [38]. This program was used to infer the number of putative clusters (K) and assign individuals into corresponding clusters. We performed this analysis under admixture model and using 105 iterations after a 105 burn-in length with K ranging from 1 to 7. The most likely K value was inferred by calculating ΔK using the method of Evanno et al. [39].

3. Results

3.1. Polymorphisms of Microsatellites

Microchecker analysis showed no evidence for null alleles and allele stuttering. Examination of the Hardy-Weinberg equilibrium (HWE) showed that two loci (CA-6 and Ec_157) [26,21] exhibited non conformity to HWE in many samples. Therefore, these two loci were excluded for subsequent analysis. Among 63 locus population cases, 11 cases significantly deviated from HWE. After sequential Bonferroni correction, only three cases showed significant (Table 2). Examination of linkage disequilibrium using Fisher’s exact test showed no disequilibrium in these microsatellite loci. A total of 165 alleles were detected across nine microsatellite loci, all of which were polymorphic. The number of alleles per locus ranged from 4.286 (GAA-1) to 17.143 (D496) with an average value of 9.651. The expected heterozygosity (He) at each locus varied from 0.460 (Mbo061) to 0.916 (D496) with a mean value of 0.716.

3.2. Intra-Population Genetic Diversity

For the Chinese samples, a statistically significant (P < 0.05) reduction of allele richness (Ar) was observed in cultured stocks of orange-spotted grouper (Ar between 5.176 and 6.674, mean = 5.821) compared to the wild populations (Ar between 7.176 and 7.454, mean = 7.311) (Table 2). The mean expected heterozygosity (He) in the two groups (Cultured and Wild) were 0.689 and 0.748, respectively. However, genetic diversity in terms of heterozygosity was not markedly reduced (P > 0.05) in contrast to the significant reduction of Ar. The Inbreeding coefficient (ƒ) value in the cultured stocks did not show significant differences with that of the wild populations of the Chinese samples (P > 0.05) (Table 2). On the other hand, for the wild populations originating from different geographical areas, the highest level of genetic diversity was detected in Malaysian population SDK (Ar = 8.421, He = 0.812), and the lowest was in Indonesian population TRK (Ar = 6.367, He = 0.641) (Table 2). Bottleneck analysis showed no evidence of recent bottleneck for any of the seven samples except for the HZ stock, which showed heterozygosity excess under IAM of microsatellites (P = 0.037). The results of bottleneck analysis most likely indicate that all these seven samples have not experienced recent bottleneck.

3.3. Inter-Population Genetic Differentiation

Pairwise Fst analysis showed significant genetic structure among all the seven samples except two comparisons where the differentiation was not statistically significant (Table 3). The genetic differentiation between the cultured stocks and the wild populations was statistically significant (Fst = 0.024, P < 0.05), although the differentiation between cultured stock SZ and wild population DYB was not significant (P > 0.05, Table 3). For the two ‘outgroup’ samples, Malaysia SDK and Indonesia TRK not only differentiated with each other (Fst = 0.075, P < 0.05), but also significantly differentiated (0.020 < Fst < 0.049, P < 0.05) with the other samples. The largest differentiation was between the wild population TRK and the cultured stock HZ (Fst = 0.129, P < 0.05), whereas the smallest divergence was presented between the two Chinese wild populations DYB and HNI (Fst = 0.002, P > 0.05). The genetic relationships suggested by pairwise Fst values were clearly visualized in the multidimensional scaling analysis. MDS results (Stress = 0.048; RSQ = 0.985) showed the cultured stocks except for SZ markedly differentiated from the wild populations (Figure 2) and also supported the ‘outgroup’ samples (SDK and TRK) differentiated from the other samples. In simulations of the Bayesian clustering method with the program Structure, the mean L(K) and ΔK all suggested three clusters as the most likely population structure, as discontinuity was presented at this K/ΔK value (Figure 3). In the results, the cultured stocks except for SZ and wild populations showed much different plotting, and SDK and TRK also showed markedly different plotting to the other samples and to one another (Figure 4). This pattern of structuring was similar to that of MDS analysis. The Mantel test showed no significant relationship between genetic distance and geographical distance among the four wild populations (r = −0.342, Z = 224.016, upper tail P = 0.779 and lower tail P = 0.221), rejecting the isolation-by-distance model of genetic structure. AMOVA results showed that statistically significant variation occurred among populations within groups (2.987 %, P < 0.05) and between cultured and wild groups of the orange-spotted grouper in China (1.381 %, P < 0.05; Table 4). This result was consistent with that of Fst analysis.

4. Discussion

4.1. Microsatellites Polymorphism

This is the first study of genetic diversity and differentiation of the orange-spotted grouper between cultured stocks and wild populations using microsatellite markers. In this study, we found that the number of alleles per locus varied between 4.286 and 17.143, and the expected heterozygosity varied from 0.460 to 0.916 (Table 2). The polymorphism level of these microsatellite loci in the orange-spotted grouper was similar to that of some other marine fish species, such as cod (Gadus morhua) [40], sea bream [41], king fish (Seriola lalandi) [42] and Asian seabass (Lates calcarifer) [19], suggesting that these polymorphic microsatellites were sufficient to reveal intraspecific diversity of the orange-spotted grouper. Also, the polymorphism level detected in this study is much higher than that in previous studies that employ microsatellite analysis of the orange-spotted grouper, for two reasons [16,43]. Firstly, we employed more microsatellite loci and each of these loci had a high polymorphism level. Secondly, we used an automatic DNA sequencer to isolate PCR products and genotype, thus could find more alleles than the conventional silver-staining method [44].

4.2. Genetic Variability Between Cultured Stocks and Wild Populations

Several studies have shown that aquaculture practices reduce genetic variability in hatchery-reared stocks of various fish species [1315]. As high genetic variation was considered to be related to the adaptive fitness in changing environments, losing genetic variation was detrimental to the domestication process of cultured stocks [45]. In the present study, significantly reduced genetic diversity in terms of Ar was also observed in cultured stocks of the orange-spotted grouper originating from the South China Sea, when compared to the wild populations of the same origin (Table 2). However, the reduction of He in cultured stocks was not statistically significant, though it seemed to have declined (Table 2). This non significant reduction in heterozygosity in cultured stocks in contrast to wild populations was also reported in other studies [10,11,46]. Allele richness is a more sensitive measure of genetic perturbations than heterozygosity [7,47]. Loss of rare alleles in populations can greatly influence Ar. However, this loss has little effect on heterozygosity. This has been backed up by previous theoretical and empirical studies [11,48,49]. In our study, the decline of genetic variation in cultured stocks could be caused mainly by founder effects. Founder effects could have occurred in breeding programs by using broodstock with a small number of individuals, which led to the loss of genetic diversity in cultured stocks. The founding of a stock with broodstock comprised of a few wild individuals is prone to have great effect on genetic diversity, especially on allelic diversity [50]. It is a fact that most of the hatcheries of the orange-spotted grouper in the region of southern China are small- or medium-sized with inadequate parental fish. This is mainly because the wild spawners, as the major source of broodstock in hatchery, were difficult to catch because of overfishing and habitat destruction. The decline of genetic variation in cultured stocks may be caused by the increased effect of genetic drift resulting from using a small number of parental individuals [11]. Besides founder effects, the genetic changes in cultured stocks compared to wild populations are likely due to artificial and natural selection existing in the culture environment [51,46]. However this selection could not be reflected using neutral markers as microsatellites in our study.
With respect to genetic variation changes of the cultured stocks relative to wild populations, we should not neglect the stock enhancement programs for some grouper species in China. In recent years, the government-sponsored hatcheries of the orange-spotted grouper have released a number of hatchery-reared fingerlings into the South China Sea [52]. However, the fingerlings’ genetic characteristics and adaptability to a wild environment were never assessed. If a large number of cultured groupers with reduced genetic variation compared to the wild populations were released into the South China Sea, the genetic composition of the wild populations in this area is likely to be effected. Therefore, these programs should be carried out with caution.
Genetic analysis involving allele richness and heterozygosity showed that genetic diversity among the orange-spotted grouper, dusky grouper (Epinephelus marginatus) and Hawaiian grouper (Epinephelus quernus) was comparable [53,26]. However, compared to migratory fishes, such as Atlantic cod [40], red sea bream [41] and king fish [37], these epinephelus fishes are less diverse. This might indicate that populations of the orange-spotted grouper are small and isolated [16]. The wild populations in China (HNI and DYB) had an intermediate value of genetic variability compared to the two ‘outgroup’ populations, being lower than the Malaysian SDK population but higher than Indonesian TRK population (Table 2). Based on this result, we conclude that the present wild population of the orange-spotted grouper in the South China Sea still maintain a normal level of genetic diversity and the population’s genetic variability is relatively unaffected, despite the population decline during the past decades owing to overfishing and habitat destruction. However, such conclusions must be made with caution, because the genetic diversity information of the wild grouper populations of China prior to population decline was not available. Consequently, a firm demonstration of this conclusion must be made according to further genetic analysis of historical and contemporary samples before and after population decline. Also, a reasonable sampling strategy encompassing sample size and distribution is needed in order to obtain more credible analytical results. Nevertheless, our results of the genetic variability between cultured stocks and wild populations of the orange-spotted grouper provided important information for aquaculture practices of this species, which was never studied before. Future studies on genetic diversity of this species should include more cultured samples and wild populations in other regions around the South China Sea to acquire a more comprehensive picture of the genetic variation of this species.

4.3. Genetic Differentiation Between and Within Cultured Stocks and Wild Populations

Pairwise Fst, Structure and AMOVA analysis revealed significant differentiation between the cultured stocks and wild populations (Table 3, Figure 4, Table 4). Significant differentiation between cultured stocks and wild populations was also observed in other fishes, such as salmon (Salmo salar) [11], grass carp (Ctenopharyngodon idella) [54,55] and Asian seabass [19]. The genetic divergence between cultured stocks and wild populations of the orange-spotted grouper may be caused by artificial selection, founder effects and random genetic drift in the cultured stocks or in the breeding program. However, artificial selection could have a minor effect on such differentiation, as the cultured stocks in our study have not been selected extensively. Founder effects should not be ruled out, because the broodstock in the hatcheries were comprised of only a few wild individuals as described above. As no evidence of bottleneck was detected in the cultured stocks, these founder effects could have been insignificant. On the other hand, the effect of random genetic drift could be an important reason for the significant differentiation. It was suggested that genetic differentiation was more influenced by a random genetic drift in a small sample than a big one [56]. The random genetic drift may also be caused by unequal sex ratio or differential reproductive contributions of the brooders in the hatcheries [57]. The orange-spotted grouper is such a species with unequal sex ratio due to late sex reversal.
Regarding the differentiation among wild populations, the insignificant differentiation between two Chinese wild populations (HNI and DYB) indicated that these two samples should be considered as one population (Table 3). This result is very consistent with a previous study of the orange-spotted grouper that no obvious difference (Fst = 0.011) was observed between two wild populations collected in different locations (Taiwan and Guangdong) of the South China Sea [17]. However, beyond the range of the South China Sea, Chinese wild grouper populations were significantly divergent from the two ‘outgroup’ populations, SDK (Malaysia) and TRK (Indonesia), according to Fst and Structure analysis (Table 3, Figure 4). Significant population structuring in a wide range of geographical areas was also observed in other grouper species, such as humpback grouper (Cromileptes altivelis) [58], dusky grouper [53], Rock grouper (Epinephelus adscensionis) [59]. This differentiation might be caused by many factors, such as biological factors including larvae dispersing ability, site fidelity, ocean currents and so on. Although the historical dispersal dynamics is a determinant key in spatially structured populations [60], it is beyond the scope of this study. The elucidation of all these factors needs a more detailed phylogeographical study.
The detection of differentiation between and within cultured stocks and wild populations in our study provided much valuable information for the breeding program of this species. The significant genetic differentiation among these wild populations might be used to produce hybrid vigour, although further experiments are required. At the same time, our results supplemented the previous population genetics studies, which showed important implications for the conservation of this species.

5. Conclusions

To our knowledge, this study was the first genetic monitoring research for cultured stocks and wild populations of the orange-spotted grouper. Our work demonstrated that genetic changes, including reduced genetic diversity and significant differentiation, have taken place in cultured grouper stocks compared to the wild populations due to founder effects and random genetic drifts during aquaculture practices. This result provided important information for ongoing breeding programs and stock enhancement programs. We also found considerable population differentiation among the orange-spotted grouper populations distributed in a wide geographical area from China, through Malaysia to Indonesia, although negligible genetic differentiation was observed between Chinese populations. This could result from the co-effects of various factors, but not merely a matter of isolation-by-distance. The genetic differences among geographical populations could provide more choice for the selective breeding work of the orange-spotted grouper.

Acknowledgments

We would like to thank Chewswee Kang (Malaysia), Yongzhong Liufu (China), Haifa Zhang (China) and a number of anonymous grouper farmers for collecting samples and Guili Zhu (Duke-NUS Graduate Medical School Singapore) for comments on the manuscript. This study was supported by the National Key Technology R & D program of China (No. 2007BAD29B03, 2007BAD29B01), as well as the National Natural Science Foundation (No. 31001112) and the National High Technology R & D program of China (No. 2006AA10A414).

References

  1. Nelson, JS. Fishes of the World; John Wiley and Sons: New York, NY, USA, 1994. [Google Scholar]
  2. Heemstra, PC; Randall, JE. FAO Fisheries Synopsis: Groupers of the World (Family Serranidae, Subfamily Epinephelinae); Food and Agriculture Organization of the United Nations: Rome, Italy, 1993; Volume 16, pp. 1–382. [Google Scholar]
  3. IUCN red list of threatened species; IUCN: Cambridge, UK, 2011. Available online: http://www.iucnredlist.org accessed on 28 June 2011.
  4. Tupper, M; Sheriff, N. Capture-based aquaculture of groupers. In FAO Fisheries Technical Paper; FAO: Rome, Italy, 2008; pp. 217–253. [Google Scholar]
  5. FAO, Cultured Aquatic Species Information Programme, Epinephelus coioides; FAO Fisheries and Aquaculture Department: Rome, Italy, 2010. Available online: http://www.fao.org accessed on 28 June 2011.
  6. Froese, R; Pauly, D. FishBase, version (06/2011); FishBase: Paris, France, 2011. Available online: http://www.fishbase.org accessed on 29 June 2011.
  7. Allendorf, FW; Phelps, SR. Loss of genetic variation in a hatchery stock of cutthroat trout. Trans. Am. Fish. Soc 1980, 109, 537–543. [Google Scholar]
  8. Lacy, RC. Loss of genetic diversity from managed populations: Interacting effects of drift, mutation, selection, and population subdivision. Conserv. Biol 1987, 2, 143–158. [Google Scholar]
  9. Lind, CE; Evans, BS; Knauer, J; Taylor, JJ; Jerry, DR. Decreased genetic diversity and a reduced effective population size in cultured silver-lipped pearl oysters (Pinctada maxima). Aquaculture 2009, 286, 12–19. [Google Scholar]
  10. Coughlan, JP; Imsland, AK; Galvin, PT; Fitzgerald, RD; Naevdal, G; Cross, TF. Microsatellite DNA variation in wild populations and farmed strains of turbot from Ireland and Norway: A preliminary study. J. Fish. Biol 1998, 52, 916–922. [Google Scholar]
  11. Norris, AT; Bradley, DG; Cunningham, EP. Microsatellite genetic variation between and within farmed and wild Atlantic salmon (Salmo salar) populations. Aquaculture 1999, 180, 247–264. [Google Scholar]
  12. Perez-Enriquez, R; Takagi, M; Taniguchi, N. Genetic variability and pedigree tracing of a hatchery-reared stock of red sea bream (Pagrus major) used for stock enhancement, based on microsatellite DNA markers. Aquaculture 1999, 173, 413–423. [Google Scholar]
  13. Sekino, M; Hara, M; Taniguchi, N. Loss of microsatellite and mitochondrial DNA variation in hatchery strains of Japanese flounder Paralichthys olivaceus. Aquaculture 2002, 213, 101–122. [Google Scholar]
  14. Lundrigan, TA; Reist, JD; Ferguson, MM. Microsatellite genetic variation within and among Arctic charr (Salvelinus alpinus) from aquaculture and natural populations in North America. Aquaculture 2005, 244, 63–75. [Google Scholar]
  15. Kohlmann, K; Kersten, P; Flajshans, M. Microsatellite-based genetic variability and differentiation of domesticated, wild and feral common carp (Cyprinus carpio L.) populations. Aquaculture 2005, 247, 253–266. [Google Scholar]
  16. Antoro, S; Na-Nakorn, U; Koedprang, W. Study of genetic diversity of orange-spotted grouper, Epinephelus coioides, from Thailand and Indonesia using microsatellite markers. Mar. Biotechnol 2006, 8, 17–26. [Google Scholar]
  17. Wang, J; Guo, F; Ding, S; Wang, J. Genetic diversity of different geographical stocks of Epinephelus coioides by microsatellite DNA. Mar. Sci 2009, 11, 60–64. [Google Scholar]
  18. Mccusker, MR; Bentzen, P. Historical influences dominate the population genetic structure of a sedentary marine fish, Atlantic wolffish (Anarhichas lupus), across the North Atlantic Ocean. Mol. Ecol 2010, 19, 4228–4241. [Google Scholar]
  19. Yue, GH; Zhu, ZY; Lo, LC; Wang, CM; Lin, G; Feng, F; Pang, HY; Li, J; Gong, P; Liu, HM; Tan, J; Chou, R; Lim, H; Orban, L. Genetic variation and population structure of Asian seabass (Lates calcarifer) in the Asia-Pacific region. Aquaculture 2009, 293, 22–28. [Google Scholar]
  20. Meng, X; Wang, Q; Jang, I; Liu, P; Kong, J. Genetic differentiation in seven geographic populations of the fleshy shrimp Penaeus (Fenneropenaeus) chinensis based on microsatellite DNA. Aquaculture 2009, 287, 46–51. [Google Scholar]
  21. Molecular Ecology Resources Primer Development Consortium. Permanent genetic resources added to molecular ecology resources database, 1 August 2009–30 September 2009. Mol. Ecol. Resour 2010, 10, 232–236.
  22. Hoelzel, AR; Green, A. Analysis of population-level variation by sequencing PCR-amplified DNA. In Molecular Genetic Analysis of Populations: A Practical Approach; IRL Press: Oxford, UK, 1992; pp. 159–187. [Google Scholar]
  23. Dong, Q; Liu, C; Guo, Y; Liu, L; Xu, T. Isolation and characterization of polymorphic microsatellite markers in Epinephelus awoara. J. Fish. China 2007, 31, 841–847. [Google Scholar]
  24. Zatcoff, MS; Ball, AO; Chapman, RW. Characterization of polymorphic microsatellite loci from black grouper, Mycteroperca bonaci (Teleostei: Serranidae). Mol. Ecol. Notes 2002, 2, 217–219. [Google Scholar]
  25. Zhu, ZY; Lo, LC; Lin, G; Xu, YX; Yue, GH. Isolation and characterization of polymorphic microsatellites from red coral grouper (Plectropomus maculatus). Mol. Eco. Notes 2005, 5, 579–581. [Google Scholar]
  26. Rivera, M; Graham, GC; Roderick, GK. Isolation and characterization of nine microsatellite loci from the Hawaiian grouper Epinephelus quernus (Serranidae) for population genetic analyses. Mar. Biotechnol 2003, 5, 126–129. [Google Scholar]
  27. Ramirez, MA; Patricia-Acevedo, J; Planas, S; Carlin, JL; Funk, SM; McMillan, WO. New microsatellite resources for groupers (Serranidae). Mol. Ecol. Notes 2006, 6, 813–817. [Google Scholar]
  28. Van Oosterhout, C; Hutchinson, WF; Wills, D; Shipley, P. MICRO-CHECKER: Software for identifying and correcting genotyping errors in microsatellite data. Mol. Ecol. Notes 2004, 4, 535–538. [Google Scholar]
  29. Hered, GJ. FSTAT (Version 1.2): A computer program to calculate F-statistics. J. Hered 1995, 86, 485–486. [Google Scholar]
  30. Mousadik, AE; Petit, RJ. High level of genetic differentiation for allelic richness among populations of the argan tree [Argania spinosa (L.) Skeels] endemic to Morocco. Theor. Appl. Genet 1996, 92, 832–839. [Google Scholar]
  31. Sokal, RR; Rohlf, FJ. Biometry: The Principles and Practice of Statistics in Biological Research, 3rd ed; W.H. Freeman and Co: New York, NY, USA, 1995; pp. 427–434. [Google Scholar]
  32. Rousset, F. GENEPOP'007: A complete re-implementation of the GENEPOP software for Windows and Linux. Mol. Ecol. Resour 2008, 8, 103–106. [Google Scholar]
  33. Rice, WR. Analyzing tables of statistical tests. Evolution 1989, 43, 223–225. [Google Scholar]
  34. Cornuet, JM; Luikart, G. Description and power analysis of two tests for detecting recent population bottlenecks from allele frequency data. Genetics 1996, 144, 2001–2014. [Google Scholar]
  35. Lewis, P; Zaykin, D. Genetic data analysis: Computer program for the analysis of allelic data, Version 1.1; Paul Lewis Phylogenetics: Storrs, CT, USA, 2002. Available online: http://hydrodictyon.eeb.uconn.edu/people/plewis/software.php accessed on 21 June 2011.
  36. Excoffier, L; Laval, G; Schneider, S. Arlequin ver. 3.0: An integrated software package for population genetics data analysis. Evol. Bioinform 2005, 1, 47–50. [Google Scholar]
  37. Miller, MP. Tools for Population Genetic Analyses (TFPGA) 1.3: A Windows Program for the Analysis of Allozyme and Molecular Population Genetic Data; Utah State University: Logan, UT, USA, 1997. [Google Scholar]
  38. Hubisz, MJ; Falush, D; Stephens, M; Pritchard, JK. Inferring weak population structure with the assistance of sample group information. Mol. Ecol. Resour 2009, 9, 1322–1332. [Google Scholar]
  39. Evanno, G; Regnaut, S; Goudet, J. Detecting the number of clusters of individuals using the software Structure a simulation study. Mol. Ecol 2005, 14, 2611–2620. [Google Scholar]
  40. Bentzen, P; Taggart, CT; Ruzzante, DE; Cook, D. Microsatellite polymorphism and the population structure of Atlantic cod (Gadus morhua) in the Northwest Atlantic. Can. J. Fish. Aquat. Sci 1996, 53, 2706–2721. [Google Scholar]
  41. Takagi, M; Taniguchi, N; Cook, D; Doyle, RW. Isolation and characterization of microsatellite loci rom red sea bream Pagrus major and detection in closely related species. Fish. Sci 1997, 63, 199–204. [Google Scholar]
  42. Nugroho, E; Ferrel, DJ; Smith, P; Taniguchi, N. Genetic divergence of kingfish from Japan, Australia and New Zealand inferred by microsatellite DNA and mitochondrial DNA control region markers. Fish. Sci 2001, 67, 843–850. [Google Scholar]
  43. Koedprang, W; Na-Nakorn, U; Nakajima, M; Taniguchi, N. Evaluation of genetic diversity of eight grouper species Epinephelus spp. based on microsatellite variations. Fish. Sci 2007, 73, 227–236. [Google Scholar]
  44. Sun, XW; Liu, CY; Liang, LQ. Isolation of microsatellite with enrichment by magnetic beads in grass carp Ctenopharyngodon idella. J. Fish. China 2005, 29, 286–482. [Google Scholar]
  45. Koljonen, ML; Tähtinen, J; Säisä, M; Koskiniemi, J. Maintenance of genetic diversity of Atlantic salmon (Salmo salar) by captive breeding programmes and the geographic distribution of microsatellite variation. Aquaculture 2002, 212, 69–92. [Google Scholar]
  46. Mjolnerod, IB; Refseth, UH; Karlsen, E; Balstad, T; Jakobsen, KS; Hindar, K. Genetic differences between two wild and one farmed population of Atlantic salmon (Salmo salar) revealed by three classes of genetic markers. Hereditas 1997, 127, 239–248. [Google Scholar]
  47. Allendorf, FW. Genetic drift and the loss of alleles versus heterozygosity. Zoo Biol 1986, 5, 181–190. [Google Scholar]
  48. Comps, B; Gomory, D; Letouzey, J; Thiebaut, B; Petit, RJ. Diverging trends between heterozygosity and allelic richness during postglacial colonization in the European beech. Genetics 2001, 157, 389–397. [Google Scholar]
  49. Skaala, O. Genetic population structure of Norwegian brown trout. J. Fish Biol 1992, 41, 631–646. [Google Scholar]
  50. Allendorf, FW; Ryman, N. Genetic management of hatchery stocks. In Population Genetics and Fishery Management; Ryman, N, Utter, F, Eds.; Washington Sea Grant Program, University of Washington Press: Seattle, WA, USA, 1987; pp. 141–159. [Google Scholar]
  51. Eknath, AE; Doyle, RW. Indirect selection for growth and life-history traits in Indian carp aquaculture: I. Effects of broodstock management. Aquaculture 1985, 49, 73–84. [Google Scholar]
  52. Oceanic and Fisheries Administrator of Guangdong Province, Guangzhou, China. 2010. Available online: http://www.gdofa.gov.cn accessed on 28 June 2011.
  53. De Innocentiis, S; Sola, L; Cataudella, S; Bentzen, P. Allozyme and microsatellite loci provide discordant estimates of population differentiation in the endangered dusky grouper (Epinephelus marginatus) within the Mediterranean Sea. Mol. Ecol 2001, 10, 2163–2175. [Google Scholar]
  54. Zhang, DC; Yu, LN; FYL. A study on genetic diversity of wild and cultured grass carp populations using RAPD. Freshwater Fish 2004, 4, 5–7. [Google Scholar]
  55. Liu, F; Xia, J; Bai, Z; Fu, J; Li, J; Yue, GH. High genetic diversity and substantial population differentiation in grass carp (Ctenopharyngodon idella) revealed by microsatellite analysis. Aquaculture 2009, 297, 51–56. [Google Scholar]
  56. Wade, MJ. Extinction and recolonization: Their effects on the genetic differentiation of local population. Evolution 1988, 42, 995–1005. [Google Scholar]
  57. Li, J; Wang, G; Bai, Z. Genetic variability in four wild and two farmed stocks of the Chinese freshwater pearl mussel (Hyriopsis cumingii) estimated by microsatellite DNA Markers. Aquaculture 2009, 287, 286–291. [Google Scholar]
  58. Sugama, K; Tridjoko, H; Budi, S; Cholik, F. Genetic variation and population structure in the humpback grouper, Cromileptes altivelis, throughout its range in Indonesian waters. Indo. Fish. Res. J 1999, 1, 32–38. [Google Scholar]
  59. Carlin, JL; Robertson, DR; Bowen, BW. Ancient divergences and recent connections in two tropical Atlantic reef fishes Epinephelus adscensionis and Rypticus saponaceous (Percoidei: Serranidae). Mar. Biol 2003, 143, 1057–1069. [Google Scholar]
  60. Hastings, A; Harrison, S. Metapopulation dynamics and genetics. Annu. Rev. Ecol. Syst 1994, 25, 167–188. [Google Scholar]
Figure 1. Sampling locations of the orange-spotted grouper in China and the Southeast Asia. HZ: cultured stock in Daya Bay, Huizhou; SZ: cultured stock in Shenzhen; ZJ: cultured stock in Zhanjiang; DYB: wild population in Daya Bay; HNI: wild population in Hainan Island; SDK: wild population in Sandakan, Malaysia; TRK: wild population in Tarakan, Indonesia.
Figure 1. Sampling locations of the orange-spotted grouper in China and the Southeast Asia. HZ: cultured stock in Daya Bay, Huizhou; SZ: cultured stock in Shenzhen; ZJ: cultured stock in Zhanjiang; DYB: wild population in Daya Bay; HNI: wild population in Hainan Island; SDK: wild population in Sandakan, Malaysia; TRK: wild population in Tarakan, Indonesia.
Ijms 12 04378f1
Figure 2. Multi-dimensional scaling plot of four wild populations and three cultured stocks of the orange-spotted grouper for genetic distribution based on pairwise Fst values.
Figure 2. Multi-dimensional scaling plot of four wild populations and three cultured stocks of the orange-spotted grouper for genetic distribution based on pairwise Fst values.
Ijms 12 04378f2
Figure 3. The log probability of data [L(K) ± SD] over 20 runs for each K and ΔK values, where the highest level of structure is suggested to be the true number of clusters.
Figure 3. The log probability of data [L(K) ± SD] over 20 runs for each K and ΔK values, where the highest level of structure is suggested to be the true number of clusters.
Ijms 12 04378f3
Figure 4. Bayesian analysis of the genetic structre based on nine microsatellite loci. Each individual is represented by a vertical line, which is coloured according to the assigned groups at estimated K = 3 (see Figure 3).
Figure 4. Bayesian analysis of the genetic structre based on nine microsatellite loci. Each individual is represented by a vertical line, which is coloured according to the assigned groups at estimated K = 3 (see Figure 3).
Ijms 12 04378f4
Table 1. Sequence of 11 pairs of microsatellite primers.
Table 1. Sequence of 11 pairs of microsatellite primers.
LocusGenBank no.Primer sequence (5′-3′)Tann (°C)Reference
D496DQ914905F: TTACTGGCAGCAATGGAC50[23]
R: GATGTATGACTACGAATGG
Mbo061AY063512F: TGAAGAATGTCAGATATTTTGTGGTG53[24]
R: TCCCAAGAGTGTTGAAGTTATAGG
Pm12AY688385F: AGAAAAAGCTCCACAACACAACAA55[25]
R: GAGCCCCAGTCCCAAATATTG
CA-2AF539606F: GACTTGATTCAGCAAAATAAAGATG55[26]
R: AGAGACGGTGCCAGTAAATGAA
CA-6AF539608F: GTGTTGCTGGGGTTACTAATGAAG50[26]
R: TTAGACACATTGTCACGATGGTCC
GAA-1AF539612F: GGAGTGTTAAATATGCCCACCA60[26]
R: CAGAAATCGCTGAGGACAAGAG
RH_GATA_003DQ223790F: GGGCAATTTGGTTCTTCACA57[27]
R: TGTCAATGCCACAGGATACA
Ec_122GQ267997F: CATTCCTTAAAGTATTCTGTG55[21]
R: CCACAGCCAGTCTAGGTATTC
Ec_154GQ429007F: AGCTGCTCAACAGGTTGTGTT56[21]
R: CAAGTTCCATATGTGCTCTGACA
Ec_157GQ429008F: TGGAACAAGTTGGCATGGTA56[21]
R: CAAATACAACACCCTAGATTTT
Ec_158GQ429009F: TGAGAGACAGTGGAGCACAAA56[21]
R: CGTGGTTACATTCTACCCCCTA
Tann, annealing temperature.
Table 2. Summary statistics of microsatellites in orange-spotted grouper.
Table 2. Summary statistics of microsatellites in orange-spotted grouper.
LocusHZ
n = 26
SZ
n = 15
ZJ
n = 34
DYB
n = 50
HNI
n = 40
TRK
n = 11
SDK
n = 25
Mean
D496A1312172325131717.143
Ar9.41710.6919.03112.64712.76412.26011.54211.193
He0.9000.9260.8310.9450.9450.9310.9340.916
ƒ0.0630.0620.0530.0520.0450.0470.038
P0.1850.9170.2700.4180.4500.0910.061
Mbo061A4561011566.714
Ar2.6923.9923.9744.9305.5084.9004.2474.320
He0.2460.4050.3430.4590.5180.5840.6670.460
ƒ−0.100−0.095−0.097−0.087−0.081−0.087−0.148
P1.0001.0001.0000.8970.9481.0000.000*
Pm12A68888797.714
Ar4.6347.2855.4696.5286.5136.8106.6746.273
He0.7270.8710.6120.8260.7840.8310.8110.780
ƒ0.0610.0360.0470.042−0.0110.0270.027
P0.2270.3950.1730.3050.0010.3510.036
CA-2A4659103187.857
Ar3.9735.5653.4704.8945.8422.81810.6275.313
He0.7320.7350.4850.6070.6720.1770.8860.613
ƒ0.0210.0430.0440.0200.0340.025−0.009
P0.5220.9220.1350.6570.5531.0000.001
GAA-1A36736234.286
Ar2.9994.8964.4602.9543.7042.0002.9793.427
He0.6570.5980.6120.5340.6060.2470.5490.543
ƒ−0.008−0.0060.0060.0190.0460.0140.016
P0.4200.1490.1010.9000.1751.0000.842
RH_GATA_003A758119577.429
Ar5.1884.5264.8035.9086.2775.0005.4625.309
He0.6980.5400.6480.6720.7310.6160.7430.664
ƒ−0.020−0.0180.017−0.027−0.032−0.0340.003
P0.4170.5100.3020.1310.0130.1310.478
Ec_122A1011111517111613.000
Ar7.4519.2977.9249.5079.70711.52810.4759.413
He0.8080.8800.8470.8870.8760.9440.8830.875
ƒ−0.006−0.032−0.019−0.050−0.044−0.039−0.054
P0.0340.3950.0460.9530.7490.6550.062
Ec_154A9710201872113.143
Ar6.9216.1186.0219.8729.7717.00012.8428.364
He0.7940.7380.7400.8880.8960.7260.9440.818
ƒ−0.087−0.091−0.065−0.169−0.138−0.109−0.117
P0.0140.0530.000*0.0940.7920.9760.873
Ec_158A49713125179.571
Ar3.3077.6955.3677.2657.0034.99110.9446.653
He0.6080.8740.7580.8210.7930.7140.8870.779
ƒ0.0400.0380.0520.0120.0140.0380.047
P0.0960.7480.0350.0780.1250.6980.000*
MeanA6.6677.6678.77812.44412.8896.44412.667
Ar5.1766.6745.6137.1677.4546.3678.421
He0.6860.7300.6530.7380.7580.6410.812
A, number of alleles; Ar, allele richness; He, expected heterozygosity; ƒ, inbreeding coefficient; P, Hardy-Weinberg equilibrium P value;
*significant departure from HWE after sequential Bonferroni correction.
Table 3. Matrix of pairwise Fst values (below diagonal) and P value (above diagonal) among four wild populations and three cultured stocks of the orange-spotted grouper based on nine microsatellite loci (Correction for multiple comparison: P < 0.0024).
Table 3. Matrix of pairwise Fst values (below diagonal) and P value (above diagonal) among four wild populations and three cultured stocks of the orange-spotted grouper based on nine microsatellite loci (Correction for multiple comparison: P < 0.0024).
PopHZSZZJDYBHNITRKSDK
HZ00.0000.0000.0000.0000.0000.000
SZ0.05300.0000.0630.0000.0000.000
ZJ0.0440.03900.0000.0000.0000.000
DYB0.0430.0100.03200.2790.0000.000
HNI0.0380.0300.0390.00200.0000.000
TRK0.1290.0700.1120.0310.04900.000
SDK0.0470.0410.0730.0320.0200.0750
Table 4. Analysis of molecular variances (AMOVA) of microsatellites between the wild group and cultured group of the orange-spotted grouper in China.
Table 4. Analysis of molecular variances (AMOVA) of microsatellites between the wild group and cultured group of the orange-spotted grouper in China.
Source of variationSum of squaresVariance componentsPercentage variationP value
Among groups18.4480.0471.3810.014
Among populations with groups43.4860.1022.9870.030
Among individuals within populations623.593−0.027−0.801−0.008
Within individuals657.0003.28496.4330.036
Total1342.5273.406100.000

Share and Cite

MDPI and ACS Style

Wang, L.; Meng, Z.; Liu, X.; Zhang, Y.; Lin, H. Genetic Diversity and Differentiation of the Orange-Spotted Grouper (Epinephelus coioides) Between and Within Cultured Stocks and Wild Populations Inferred from Microsatellite DNA Analysis. Int. J. Mol. Sci. 2011, 12, 4378-4394. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms12074378

AMA Style

Wang L, Meng Z, Liu X, Zhang Y, Lin H. Genetic Diversity and Differentiation of the Orange-Spotted Grouper (Epinephelus coioides) Between and Within Cultured Stocks and Wild Populations Inferred from Microsatellite DNA Analysis. International Journal of Molecular Sciences. 2011; 12(7):4378-4394. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms12074378

Chicago/Turabian Style

Wang, Le, Zining Meng, Xiaochun Liu, Yong Zhang, and Haoran Lin. 2011. "Genetic Diversity and Differentiation of the Orange-Spotted Grouper (Epinephelus coioides) Between and Within Cultured Stocks and Wild Populations Inferred from Microsatellite DNA Analysis" International Journal of Molecular Sciences 12, no. 7: 4378-4394. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms12074378

Article Metrics

Back to TopTop