Development & Reproduction
Korean Society of Developmental Biology
Short communication

Genetic Distances Within-Population and Between-Population of Tonguesole, Cynoglossus spp. Identified by PCR Technique

Jong-Man Yoon
Corresponding Author : Jong-Man Yoon. Dr. Prof. Dept. of Aquatic Life Medicine, College of Ocean Science and Technology, Kunsan National University, Gunsan 54150, Korea. Tel.: +82-63-469-1887, E-mail: jmyoon@kunsan.ac.kr

© Copyright 2019 The Korean Society of Developmental Biology. This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Received: Jul 03, 2019 ; Revised: Aug 14, 2019 ; Accepted: Sep 05, 2019

Published Online: Sep 30, 2019

ABSTRACT

The higher fragment sizes (>2,100 bp) are not observed in the two C. spp. populations. The six oligonucleotides primers OPA-11, OPB-09, OPB-14, OPB-20, OPC-14, and OPC-18 were used to generate the unique shared loci to each tonguesole population and shared loci by the two tonguesole populations. The hierarchical polar dendrogram indicates two main clusters: Gunsan (GUNSAN 01–GUNSAN 11) and the Atlantic (ATLANTIC 12–ATLANTIC 22) from two geographic populations of tonguesoles. The shortest genetic distance displaying significant molecular difference was between individuals’ GUNSAN no. 02–GUNSAN no. 01 (genetic distance=0.038). In the long run, individual no. 02 of the ATLANTIC tonguesole was most distantly related to GUNSAN no. 06 (genetic distance=0.958). These results demonstrate that the Gunsan tonguesole population is genetically different from the Atlantic tonguesole population. The potential of PCR analysis to identify diagnostic markers for the identification of two tonguesole populations has been demonstrated. As a rule, using various oligonucleotides primers, this PCR method has been applied to identify polymorphic/specific markers particular to species and geographical population, as well as genetic diversity/polymorphism in diverse species of organisms.

Keywords: Cynoglossus spp.; Euclidean genetic distance; Genetic cluster; Pairwise population matrix; Polar dendrogram

INTRODUCTION

Cynoglossus spp. is, biologically warm water teleost species, belonging to family Cynoglossidae, order Pleuronectiformes, class Actinopterygii, widely distributed on the coast of the Yellow Sea, southern sea in the Korean Peninsula, as well as in several areas in East China and the Bo Hai of China. Under the natural ecosystem, the tonguesoles inhabit in the flats consisting of a lot of mud, silt and slime in the coast where the clams, crabs and the other shellfishes are abundant. Tonguesole is one of the most popular marine products and economically important aquacultural species in Korea because of their taste and nourishing value, and Koreans consume them plentifully. The tonguesoles live buried in sandy mud at a number of depths. However, in spite of their economic and scientific consequences, a little information currently exists regarding the physiological and ecological levels of tonguesole species in Korea. Nevertheless, this tonguesole, which is distinguished important reproductively (Lee et al., 2000; Kang et al., 2012; Lim et al., 2016), morphologically (Choi, 1996; Gu et al., 2009; Kwun & Kim, 2010), as well as ecologically (Choi et al., 1995; Baeck & Huh, 2004) are a few genetically and/or molecular-biologically (Kwun & Kim, 2010; Noh et al., 2016) studied like other finfishes. There are more wishes to understand the genetic traits and conformation of this fish population in order to evaluate precisely the genetic consequences.

The specific markers particular to species, genus or geographical population have all been employed in the identification and discrimination of individuals, species and populations, hybrid parentages and genetic diagnostics (Callejas & Ochando, 1998; Huang et al., 2000; Kim et al., 2000; Yoon & Kim, 2004). Thus, the applications of PCR techniques to aquaculture had been to identify genetic distances among a few of fish species and/or invertebrates apart from geographic sites (Partis & Wells, 1996; Callejas & Ochando, 1998; McCormack et al., 2000; Kwun & Kim, 2010; Song & Yoon, 2013; Oh & Yoon, 2014; Yoon, 2015; Noh et al., 2016). In the present study, to elucidate the genetic distances and differences between geographical tonguesole populations, the author performed a clustering analysis of two tonguesole populations collected from two sites of Gunsan and the Atlantic.

MATERIALS AND METHODS

Muscle tissues were collected separately from Gunsan and the Atlantic tonguesole populations, respectively. These tonguesole muscles were collected in sterile tubes, immediately placed on ice, and stored at −80°C until necessary. PCR analysis was performed on the muscle extracts from 22 individuals, using six oligonucleotides primers. DNA extraction should be performed according to the separation and extraction methods as described previously (Song & Yoon, 2013). 600 μL of chloroform was added to the mixture and then inverted (no phenol). After washing several washings, samples of muscle tissues were placed into 10 mL test tubes, to which 3 volumes of lysis buffer (155 mM NH4Cl, 10 mM KHCO3, 1mM EDTA) was added, and the mixture tubes were lightly inverted. The precipitates obtained were diffused with lysis buffer (10 mM Tris-HCl, pH 8.0, 10 mM EDTA, 100 mM NaCl, 0.5% SDS). Samples were added 15 μL proteinase K solutions (10 mg/mL). After incubation, there was added 300 μL of 6 M NaCl and softly pipetted for a few of min. The DNA pellets were incubation-dried for 2 hrs, held at −80°C until analysis, and then dissolved in the TE buffer (10 mM Tris-HCl, pH 8.0; 1 mM EDTA). Ice-cold 70% ethanol was added, and then the samples were centrifuged at 19,621 g for 5 minutes to extract DNA from the lysates. The concentration of the extracted genomic DNA was measured by the optical density at 260 nm, with a spectrophotometer (Beckman Coulter, Buckinghamshire, UK). The DNA pellets were then incubation-dried for more than 12 hours, maintained at −80°C until needed and then dissolved in the distilled water.

The oligonucleotides primers were purchased from Operon Technologies, USA. The OPA-11 (5’-CAATCGCCGT-3’), OPB-09 (5’-TGGGGGACTC-3’), OPB-14 (5’-TCCGCTCTGG-3’), OPB-20 (5’-GGACCCTTAC-3’), OPC-14 (5’-TGCGTGCTTG-3’) and OPC-18 (5’-TGAGTGGGTG-3’) were shown to yield the unique shared loci to each population and the number of shared loci by the two populations which could be clearly scored. The author used these oligonucleotides primers to determine the genetic variations, DNA polymorphisms, and similarity of the tonguesole from Gunsan and the Atlantic. PCR process was performed using a Programmable DNA Thermal Cyclers (MJ Research Inc., Waltham, MA, USA). DNA amplification was performed in 25 μL samples, which contained 10 ng of template DNA, 20 μL of premix (Bioneer Corp., Daejeon, Korea), and 1 unit of primer. Amplification products were generated via electrophoresis on 1.4% agarose (Bioneer Corp., Daejeon, Korea) gel containing TBE (90 mM Tris, pH 8.5; 90 mM borate; 2.5 mM EDTA). The 100 bp ladder marker (Bioneer Corp., Daejeon, Korea) was utilized as a DNA molecular weight marker. The electrophoresed agarose gels were illuminated by ultraviolet rays, and photographed using a photoman direct copy system (PECA Products, Beloit, WI, USA).

Similarity matrix including Bandsharing (BS) values were calculated according to the presence/absence of amplified products at specific positions in the same gel from the DNA profiles. The degree of variability was calculated by use of the Dice coefficient (F), which is given by the formula: F=2 nab / (na+nb), where nab is the number of bands shared between the samples a and b, na is the total number of bands for sample a and nb is the total number of bands for sample b (Jeffreys & Morton, 1987; Yoke-Kqueen & Radu, 2006). The average of within-population similarity was calculated by pairwise comparison between individuals within a population. The relatedness between different individuals in the tonguesole populations was generated according to the bandsharing values and similarity matrix. A hierarchical polar dendrogram was assembled using similarity matrices to yield a cluster tree, which was supported by the Systat version 10 (SPSS Inc., Chicago, IL, USA). The Systat software was also used to analyze genetic distances, Euclidean distances within and between tonguesole populations, means, standard errors, and t-test scores.

RESULTS AND DISCUSSION

Here, the six oligonucleotides primers OPA-11, OPB-09, OPB-14, OPB-20, OPC-14, and OPC-18 were used to generate the unique shared loci to each population and shared loci by the two populations. The complexity of the banding patterns varied dramatically between the primers from the two zones. The size of the DNA fragments also varied widely, from 100 to 2,050 bp (Fig. 1). The fragment numbers in each size interval have been computed from the integrated fragments attained with every six oligonucleotides primers. The higher fragment sizes (>2,100 bp) are not observed in the two C. spp. populations. Eighty fragments ranging from 200 bp to 2,200 bp were unambiguously counted in the black tiger shrimp (Penaeus monodon) (Tassanakajon et al., 1998). The DNA fragments obtained using the four primers ranged from 100 to 2,300 bp in the brittle star (Amphiura filiformis) (McCormack et al., 2000). The bandsharing value between individuals no. 01 and no. 02 was 0.962, which was the highest value identified within the Gunsan tonguesole population. The bandsharing value between individual’s no. 11 of the Gunsan tonguesole and no. 20 of the Atlantic tonguesole was 0.204, which was the lowest observed. The 22 unique shared loci to each population generated by OPA-11 oligonucleotides primer were in Gunsan tonguesole population, as summarized in Table 2. Excitingly, the oligonucleotides primer OPC-14 generated 44 unique shared loci to each population, which were identifying each species in the Atlantic population (Table 2). Especially, the oligonucleotides primer OPC-18 generated 33 unique shared loci to each population, which were classifying each population in the Atlantic tonguesole. Unusually, the primer OPB-20 detected 22 shared loci by the two populations, which were identical in all samples, as summarized in Table 2. The oligonucleotides primer OPC-18 also detected 66 shared loci by the two populations, major and/or minor fragments of 400 bp, 800 bp, and 1,600 bp in all samples.

dr-23-3-297-g1
Fig. 1. Distribution of fragment sizes of Gunsan and the Atlantic population of Cynoglossus spp.. Solid grey lines: Gunsan population. Solid black lines: the Atlantic population. The fragment numbers in each size interval have been computed from the integrated fragments attained with every six oligonucleotides primers. The higher fragment sizes (>2,100 bp) are not observed in the two C. spp. populations.
Download Original Figure
Table 1. Similarity matrix, including bandsharing values, calculated using Nei and Li’sindex, of the similarity of tonguesoles, Cynoglossus spp. from Gunsan and the Atlantic, respectively
Bandsharing values
from Gunsan from the Atlantic
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22
1 - 0.962 0.800 0.945 0.806 0.793 0.904 0.783 0.833 0.860 0.570 0.397 0.388 0.402 0.470 0.422 0.386 0.350 0.366 0.292 0.423 0.376
2 - 0.812 0.890 0.895 0.862 0.946 0.821 0.824 0.921 0.550 0.380 0.395 0.387 0.539 0.464 0.428 0.411 0.430 0.324 0.886 0.394
3 - 0.856 0.872 0.787 0.833 0.779 0.759 0.796 0.589 0.913 0.905 0.446 0.473 0.927 0.852 0.773 0.782 0.545 0.908 0.788
4 - 0.843 0.821 0.959 0.839 0.861 0.915 0.598 0.792 0.773 0.802 0.481 0.911 0.421 0.767 0.798 0.604 0.913 0.819
5 - 0.863 0.821 0.826 0.830 0.870 0.556 0.379 0.373 0.386 0.382 0.407 0.371 0.333 0.363 0.236 0.364 0.326
6 - 0.826 0.804 0.819 0.775 0.638 0.382 0.323 0.384 0.459 0.422 0.785 0.710 0.399 0.324 0.260 0.389
7 - 0.840 0.844 0.928 0.568 0.395 0.358 0.378 0.493 0.398 0.361 0.342 0.361 0.286 0.390 0.350
8 - 0.795 0.900 0.596 0.314 0.311 0.297 0.421 0.361 0.324 0.324 0.319 0.247 0.319 0.313
9 - 0.890 0.383 0.406 0.374 0.395 0.426 0.412 0.370 0.356 0.370 0.261 0.335 0.338
10 - 0.696 0.538 0.510 0.520 0.502 0.442 0.401 0.362 0.396 0.310 0.575 0.453
11 - 0.408 0.405 0.363 0.368 0.332 0.322 0.274 0.353 0.204 0.272 0.225
12 - 0.876 0.844 0.730 0.845 0.829 0.837 0.850 0.725 0.811 0.849
13 - 0.779 0.763 0.712 0.688 0.721 0.818 0.718 0.720 0.588
14 - 0.777 0.735 0.661 0.670 0.708 0.653 0.747 0.676
15 - 0.750 0.693 0.706 0.807 0.658 0.811 0.785
16 - 0.878 0.848 0.883 0.771 0.793 0.877
17 - 0.888 0.805 0.680 0.808 0.781
18 - 0.866 0.703 0.714 0.838
19 - 0.639 0.739 0.842
20 - 0.597 0.681
21 - 0.893
22 -
Download Excel Table
Table 2. The number of unique shared loci to each population and number of shared loci by the two populations generated by DNA analysis using 6 oligonucleotides primers in of tonguesole, Cynoglossus spp. from Gunsan and the Atlantic, respectively
Item No. of unique shared loci to each population No. of shared loci by the two populations
Primer \ Population Gunsan Atlantic Two populations
OPA-11 22 33 0
OPB-09 11 11 0
OPB-14 11 11 0
OPB-20 11 11 22
OPC-14 22 44 0
OPC-18 11 33 66
Total no. 88 143 88
Average no. per primer 14.67 23.83 14.67
Download Excel Table
Table 3. Multiple calculations of average bandsharing values (mean±SE) between two tonguesole populations were generated in line with the bandsharing values and similarity matrix
Population Gunsan Atlantic
Gunsan 0.799±0.016a 0.452±0.016c
Atlantic - 0.765±0.011b

Each value is a result of three different experiments.

a-c Values with different superscript are significantly different, p<0.05.

Download Excel Table

Polymorphisms are determined by the banding patterns of primer-amplified products at specific positions (Tassanakajon et al., 1998; Nozaki et al., 2000; Yoon & Kim, 2001). Several researchers have studied the sizes of DNA fragments in the PCR profiles of five species of Eastern Pacific abalone (genus Haliotis) (Muchmore et al., 1998), black tiger shrimp (Penaeus monodon) (Tassanakajon et al., 1998), the brittle star (Amphiura filiformis) (McCormack et al., 2000) and wild shrimp populations (Yoon & Kim, 2003a). The specific primer proved useful in the identification of individuals and/or populations, resulting from variations in DNA polymorphisms among individuals/populations (Liu et al., 1998; Yoon & Kim, 2003b; Kim et al., 2004; Yoon & Kim, 2004; Kwun & Kim, 2010; Noh et al., 2016). The PCR method, using several primers, was applied to the identification of three endemic Spanish barbel species: Barbus bocagei, B. graellsii and B. sclateri (Callejas & Ochando, 1998). Results indicated that Barbus bocagei and B. graellsii were more closely related to each other than they were to B. sclateri. In this study, PCR analysis has revealed a significant genetic distance between two population/species couples. PCR enabled to detect the existence of population/species discrimination and genetic variation in two geographic populations of tonguesole. The hierarchical polar dendrogram indicates two main cluster 1 (GUNSAN 01–GUNSAN 11) and cluster 2 (ATLANTIC 12–ATLANTIC 22) from two geographic tonguesole populations.

The shortest genetic distance displaying significant molecular difference was between individuals’ GUNSAN no. 02 and GUNSAN no. 01 (genetic distance=0.038). In the long run, individual no. 02 of the ATLANTIC tonguesole was most distantly related to GUNSAN no. 06 (genetic distance=0.958). These results demonstrate that the Gunsan tonguesole population is genetically different from the Atlantic tonguesole population. Phylogenetic relationships among 5 Haliotis species and one hybrid were conducted by calculation of the distance coefficient and construction of a phylogenetic tree based on PCR data (Kim et al., 2000). These branched off into two clusters: cluster I was formed by H. discus hannai, H. discus, H. gigantea, H. sieboldii, and the hybrid, which was subsequently re-divided into two sub-clusters. As mentioned above, the potential of PCR analysis to identify diagnostic markers for the identification of two tonguesole populations has been demonstrated. In teleost, a phylogenetic tree was constructed using UPGMA cluster analysis based on a total of 3,744 distinguishable fragments in gynogenetic clones from the silver crucian carp, Carassius auratus gibelio Block (Zhou et al., 2000). The inter-population similarity indices and genetic distance values indicated that the Jamuna-Padma population pair of the Indian major carp (Catla catla) was genetically closer than the Halda-Jamuna and the Halda-Padma population pairs, which agreed with the geographical distances between them (Islam et al., 2005). The two phylogenetic trees resulting from the neighbor-joining and parsimony analyses showed the same topology in distinguishing the Mullidae species (Mamuris et al., 1999). This confirms that the method is a suitable tool for DNA comparisons, both within and between individuals, species, and populations. As a rule, using various oligonucleotides primers, this PCR method has been applied to identify specific markers specific to species and geographical population, as well as genetic diversity/polymorphism in diverse species of organisms (McCormack et al., 2000; Yoon & Kim, 2003a; Kim et al., 2004; Kwun & Kim, 2010; Noh et al., 2016).

dr-23-3-297-g2
Fig. 2. Hierarchical polar dendrogram of genetic distances, obtained from two geographic populations of tonguesole, Cynoglossus spp. from Gunsan (GUNSAN 01–GUNSAN 11) and the Atlantic (ATLANTIC 12–ATLANTIC 22), respectively. The relatedness between different individuals in the tonguesoles was generated according to the bandsharing values and similarity matrix.
Download Original Figure

CONFLICT OF INTERESTS

The authors declare no potential conflict of interest.

ACKNOWLEDGEMENTS

The author would like to thank the referees who aided with thorough and insightful rectification. The author also states my appreciation to my laboratory colleagues and undergraduate students for their help in sample collection, and their help with the PCR techniques during this process.

REFERENCES

1.

Baeck GW, Huh SH. 2004; Age and growth of red tongue-sole (Cynoglossus joyneri) in the southern sea of Korea. J Korean Fish Socs. 37:307-311

2.

Callejas C, Ochando MD. 1998; Identification of Spanish barbel species using the RAPD technique. J Fish Biol. 53:208-215

3.

Choi Y. 1996; Two foreign species of tongue soles (Cyno-glossidae). J Korean Fish Soc. 29:267-269.

4.

Choi Y, Kim IS, Ryu BS, Chung EY, Park JY. 1995; Ecology of Cynoglossus joyneri Guenther from the western coast of Korea. Korean J Ichthyol. 7:56-63.

5.

Gu SO, Hwang HJ, Kim JK. 2009; Morphological descrip-tion of Cynoglossus joyneri larvae collected from the Yellow Sea. Korean J Fish Aquat Sci. 42:725-728

6.

Huang BX, Peakall R, Hanna PJ. 2000; Analysis of genetic structure of blacklip abalone (Haliotis rubra) populations using RAPD minisatellite and microsatellite markers. Mar Biol. 136:207-216

7.

Islam MS, Ahmed ASI, Azam MS, Alam MS. 2005; Genetic analysis of three river populations of Catla catla (HAMILTON) using randomly amplified polymorphic DNA markers. Asian-Australas J Anim Sci. 18:453-457

8.

Jeffreys AJ, Morton DB. 1987; DNA fingerprints of dogs and cats. Anim Genet. 18:1-15

9.

Kang HW, Lim HK, Kang DY, Han HS, Do YH, Park JS. 2012; Maturation and spawning of the female tongue sole Cynoglossus semilaevis in the West Coast of Korea. Dev Reprod. 16:87-93.

10.

Kim JY, Park CY, Yoon JM. 2004; Genetic differences and DNA polymorphism in oyster (Crassostrea spp.) analysed by RAPD-PCR. Korean J Genet. 26:123-134.

11.

Kim SK, Jung YH, Han SH, Oh YS, Ko MH, Oh MY. 2000; Phylogenetic relationship among Haliotis spp. distributed in Korea by the RAPD analysis. Korean J Genet. 14:179-185.

12.

Koh MC, Lim CH, Chua SB, Chew ST, Phang STW. 1998; Random amplified polymorphic DNA (RAPD) finger-prints for identification of red meat animal species. Meat Sci. 48:275-285

13.

Kwun HJ, Kim JK. 2010; Validation of morphology-based identification of two Cynoglossidae larvae using mito-chondrial DNA. Korean J Fish Aquat Sci. 43:482-488

14.

Lee JS, Kim SY, Ma KH, Huh SH. 2000; Ovarian develop-ment and reproductive cycle of the female red tongue sole Cynoglossus joyneri (Teleostei: Cynoglossidae). J Korean Fish Soc. 33:554-558.

15.

Lim HK, Park JY, Kang HW. 2016; Effect of luteinising hormone releasing hormone analogue and human cho-rionic gonadotropin on spermiation in the tongue sole Cynoglossus semilaevis. J Fish Mar Sci Educ. 28:596-602

16.

Liu Z, Li P, Argue BJ, Dunham RA. 1998; Inheritance of RAPD markers in channel catfish (Ictalurus punctatus), blue catfish (I. furcatus), and their F1, F2 and backcross hybrids. Anim Genet. 29:58-62

17.

Mamuris Z, Stamatis C, Bani M, Triantaphyllidis C. 1999; Taxonomic relationships between four species of the Mullidae family revealed by three genetic methods: Allozymes, random amplified polymorphic DNA and mitochondrial DNA. J Fish Biol. 55:572-587

18.

McCormack GP, Powell R, Keegan BF. 2000; Comparative analysis of two populations of the brittle star Amphiura filiformis (Echinodermata: Ophiuroidea) with different life history strategies using RAPD markers. Mar Biotechnol. 2:100-106.

19.

Muchmore ME, Moy GW, Swanson WJ, Vacquier VD. 1998; Direct sequencing of genomic DNA for charac-terization of a satellite DNA in five species of Eastern Pacific abalone. Mol Mar Biol Biotechnol. 7:1-6.

20.

Noh ES, Kang HS, An CM, Park JY, Kim EM, Kang JH. 2016; Rapid and specific identification of genus Cynoglossus by multiplex PCR assays using species-specific derived from the COI region. J Life Sci. 26:1007-1014

21.

Nozaki T, Mishiba K, Mii M, Koba T. 2000; Construction of synteny groups of Brassica alboglabra by RAPD markers and detection of chromosome aberrations and distorted transmission under the genetic background of B. campestris. Theor Appl Genet. 101:538-546

22.

Tassanakajon A, Pongsomboon S, Jarayabhand P, Klinbunga S, Boonsaeng VV. 1998; Genetic structure in wild popu-lations of black tiger shrimp (Penaeus monodon) using randomly amplified polymorphic DNA analysis. J Mar Biotech. 6:249-254.

23.

Yoke-Kqueen C, Radu S. 2006; Random amplified poly-morphic DNA analysis of genetically modified organisms. J Biotechnol. 127:161-166

24.

Yoon JM. 2009; Geographic variations in four freshwater crab (Eriocheir sinensis) populations throughout its distribution range. Dev Reprod. 13:97-103.

25.

Yoon JM. 2015; Geographical variations and genetic dis-tances of three Saxidomus purpuratus populations ascertained by PCR analysis. Dev Reprod. 19:259-264

26.

Yoon JM. 2018; Genetic variations of intra- and between-razor clam Solen corneus population identified by PCR analysis. Dev Reprod. 22:193-198

27.

Yoon JM, Kim GW. 2003a; Genetic differences between cultured and wild penaeid shrimp (Penaeus chinensis) populations analysed by RAPD-PCR. Korean J Genet. 25:21-32.

28.

Yoon JM, Kim JY. 2004; Genetic differences within and between populations of Korean catfish (S. asotus) and bullhead (P. fulvidraco) analysed by RAPD-PCR. Asian-Australas J Anim Sci. 17:1053-1061

29.

Yoon JM, Kim YH. 2003b; Wide marsh clam (Corbicula spp.) populations from three sites analysed by RAPD-PCR-AGE. Bull Electrochem. 19:337-348.

30.

Zhou L, Wang Y, Gui JF. 2000; Analysis of genetic heterogeneity among five gynogenetic clones of silver crucian carp Carassius auratus gibelio Block, based on detection of RAPD molecular markers. Cytogenet Cell Genet. 88:133-139