Development and Validation of Single Nucleotide Polymorphism (SNP) Markers from an Expressed Sequence Tag (EST) Database in Olive Flounder (Paralichthys olivaceus)

Jung Eun Kim, Young Mee Lee, Jeong-Ho Lee†, Jae Koo Noh, Hyun Chul Kim, Choul-Ji Park, Jong-Won Park, Kyung-Kil Kim
Author Information & Copylight
Corresponding Author : Jeong-Ho Lee, Genetics and Breeding Research Center, NFRDI, Geoje 656-842, Korea. Tel : +82-55-639-5811, Fax : +82-55-639-5809, E-mail :

Copyright © 2014 © Copyright an Official Journal of the Korean Society of Developmental Biology. All Rights Reserve. 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: Nov 28, 2014 ; Revised: Dec 5, 2014 ; Accepted: Dec 8, 2014


To successful molecular breeding, identification and functional characterization of breeding related genes and development of molecular breeding techniques using DNA markers are essential. Although the development of a useful marker is difficult in the aspect of time, cost and effort, many markers are being developed to be used in molecular breeding and developed markers have been used in many fields. Single nucleotide polymorphisms (SNPs) markers were widely used for genomic research and breeding, but has hardly been validated for screening functional genes in olive flounder. We identified single nucleotide polymorphisms (SNPs) from expressed sequence tag (EST) database in olive flounder; out of a total 4,327 ESTs, 693 contigs and 514 SNPs were detected in total EST, and these substitutions include 297 transitions and 217 transversions. As a result, 144 SNP markers were developed on the basis of 514 SNP to selection of useful gene region, and then applied to each of eight wild and culture olive flounder (total 16 samples). In our experimental result, only 32 markers had detected polymorphism in sample, also identified 21 transitions and 11 transversions, whereas indel was not detected in polymorphic SNPs. Heterozygosity of wild and cultured olive flounder using the 32 SNP markers is 0.34 and 0.29, respectively. In conclusion, we identified SNP and polymorphism in olive flounder using newly designed marker, it supports that developed markers are suitable for SNP detection and diversity analysis in olive flounder. The outcome of this study can be basic data for researches for immunity gene and characteristic with SNP.


The olive flounder (Paralichthys olivaceus) is an economically important fish in Korea and Japan, and total yields of aquacultured olive flounder increased and reached to 40,805 tons in 2011 (FAO) in Korea. Recently, breeding research is focused only culture techniques and development of research also insufficient to further increase production. Higher breeding research including the development of molecular markers is essential to produce of good breeds and high value fisheries and to enhance our competitiveness in fisheries industry. Identification and functional characterization of economic traits (or breeding traits) related genes, the development of molecular breeding techniques using DNA markers to molecular breeding is meaningful (Hubert et al., 2009; Yoon, 2006). Moreover, application of these genes in the aquaculture (molecular selective breeding) has attracted worldwide attention.

Recently, large-scale DNA sequencing such as BAC library projects and EST projects have provided usable data for a great variety of species (Hayes et al., 2005). In particular, the sequence analysis was significantly increased in economic worth of fish. This information can be used to increase the accuracy of selection breeding research for economical traits, thereby increasing the rate of genetic gain and production efficiency.

The molecular marker is an important element for breeding research, and there are various methods to used for the research. The kind of molecular techniques are RAPD (random-amplified polymorphic DNA), RFLP (restriction fragment length polymorphism), AFLP (amplified fragment length polymorphism), MS (microsatellites) and SNP (single nucleotide polymorphisms) (Liu & Cordes, 2004). Among possible molecular markers, microsatellites have been most useful for population genotyping technology, but is not suitable to high-quality data about influence gene expression and the functions. Unlike MS, single nucleotide polymerphisms (SNPs) are the most abundant type of variation in DNA sequence among individuals in a population (Fahrenkrug et al., 2002). Moreover, these markers appear to have low mutation rates (a mutation occurs every 2 × 108 nucleotides) and hence are extremely stable than other genetic markers (The International SNP Map Working Group, 2001). Although, SNPs have been recently proposed as an alternative to microsatellites for genotyping individuals (Vignal et al., 2002), have a disadvantage of time and cost-consuming when a certain selection mechanism does not exist. Whereas, if developed SNPs from of expressed sequence tags (ESTs) can be particularly valuable, also it is possible to reduce the time and cost for developing a useful SNP marker. In aquaculture species, SNPs are especially important if they cause differences in economic traits, or are linked to the mutations that do so.

For some aquaculture species, a large number of SNPs was developed from expressed sequence tag (EST) database as below; catfish (He et al., 2003), salmon (Hayes et al., 2007; Campbell & Narum, 2011; Everett et al., 2011), sea bream (Hayes et al., 2007; Campbell & Narum, 2011; Everett et al., 2011), Atlantic cod (Hubert, 2009), and turbot (Vera et al., 2011). But, SNP marker was hardly developed and validated in olive flounder even economically important fish species in Korea (He et al., 2011; He et al., 2012).

Wild and cultured olive flounder have different growth rate and resistance to disease because of differences in the growth environment. For example, wild olive flounder has generally small size and slow growth but strong immunity, whereas cultured have generally rapid growth and heavy weight but weak immunity. Each feature is associated with a SNP, and therefore SNP analysis using wild and culture samples can be a help to find regions involved superior traits. In this study, we develop useful marker to detect the single nucleotide polymorphism (SNP) from EST in olive flounder and then apply to sample of wild and cultured. Moreover, we were identified substitution of SNP and population diversity analysis using developed marker.


1. EST database and candidate SNP detection
1) ESTs database assemble

In this study used polymorphic sites were identified from expressed sequence tags (ESTs) derived from P. olivaceus. Using to the previous studied EST data, bioinformatic analysis was conducted to determine gene identities using Gene Master Software (Ensoltek). ESTs were assembled in clusters of contiguous sequences (contig) using ICAtools program (Parsons, 1995). Gene annotation procedures and homology searches of the sequenced ESTs have been locally done by BLASTX for amino acid similarity comparisons. Matches with the Expect value (E) less than 10xe-4 were considered to be significant. All ESTs that were not identified as orthologs of known genes were designated as unknown EST clones.

2) Candidate SNP detection

Total of 755 SNP genes were selected from chromosome 24. The selection of SNP gene criteria is as below; the frequency of minor allele is 30 percent, each gene is related with immune and growth, and evenly exists on chromosome 24. As a result, 144 markers were selected, developed the primer for detecting of specific gene. Primers of SNPs were designed using Primer 3 software (Rozen and Skaletsky, 2000) adjusting length products to 100-300 bp (Table 1).

Table 1. Sequences of primers used for the PCR
Marker name Forward primer Reverse primer Amplicon size
Download Excel Table
2. Validation of SNP in olive flounder
1) Sampling, amplification and analysis

Wild and cultured each of eight olive flounder were randomly obtained from East Sea and Jeju Island in Korea, respectively. Genomic DNA was extracted from pectoral fin with the standard method of phenol-chloroform (He et al., 2008) and a total of 144 primers were designed.

Amplification was carried out in the following conditions; in 20 uL volume including 2~4 uL template DNA (100~ 150 ng), 1× Gold Taq buffer (Amplified Biosystems), 2.5 mM of MgCl2, 0.2 mM of dNTP (Solgent), 10 pmol of both forward and reverse primers, 1.25 U of Gold taq DNA polymerase (Amplified biosystems) and 1 uL dimethyl sulfoxide (Sigma). Thermal cycling was conducted on a Applied Biosystems (Model #: 9902) as following conditions: initial denaturation at 95°C for 2 min, 35 cycles of denaturation at 95°C for 30 s, annealing temperature at 58°C for 30 s, and extension at 72°C for 30 s. Final extension step was performed at 72°C for 10 min. Amplification of each SNP was checked on 2% agarose gel using electrophoresis. PCR product was purified using DNA purification Kit, sequenced in forward direction with the Big Dye termination reaction chemistry kit (Nucleogen), and analyzed with an ABI 3730xl Genetic Analyzer (Applied Biosystems).

Sequences were aligned with the CLUSTAL W algorithm using MEGA 5.05 software. Heterozygosity was estimated using Arlequin software 3.5 and polymorphic information content (PIC) value using Power Marker V.3.25 software.


1. Identification of SNPs from EST database

The olive flounder EST-derived database contains 4,327 ESTs and identified 693 contigs and 514 SNP were detected. The average of SNP frequency was 1 per 216 bp of contig sequences. These substitutions included 297 transitions (C/T, 155; A/G, 142) and 217 transversions (A/T, 56; A/C, 57; T/G, 62; C/G, 42), with C/T rate (155, 52.8%) has captured most of substitutions and C/G being the least common (42, 19.4%) alteration observed (Fig. 1). A transition/ transversion ratio (ts/tv) was 1.369, SNP substitutions result of development marker shows transitions rate higher than transversions rate (Table 2).

Table 2. Summary statistics of SNP identification from olive flounder EST resources
Parametrs EST resources SNPs
Total number of ESTs searched 4,327
Number of contigs 693
Total SNPs detected 514
Total number of transitions 297
C/T 155
A/G 142
Total number of transversions 217
A/T 56
A/C 57
T/G 62
C/G 42
Total number of indels 0
Transition/transversion 1.369
Download Excel Table
Fig. 1. Substitution types detected in the EST database of olive flounder.
Download Origianl Figure

We were applied the designed primers on 16 sample (wild 8, cultured 8), and were performed PCR for amplifycation, with the result that 144 putative SNPs were selected for validation. Out of the 144 primer pairs, 25 primer did not amplify any product (17.4%) and 36 primer failed to sequencing (25.0%). Cause of failed sequencing is nonspecific bands (8, 5.6%), low PCR expression (15, 10.4%) and low alignment efficiency (13, 9.0%). But 83 markers (57.6%) were successfully analyzed in sequencing analysis. 51 out of 83 were confirmed monomorphic SNPs from the sample (wild 8, cultured 8), only 32 markers had detected polymorphic. 32 polymorphism SNP were observed in the variety genes region involved in immunity and protein synthesis (Table 4). The SNPs polymorphic are shown a total of 21 transitions (C/T and A/G), 11 transversions (A/T, A/C, T/G and C/G), and indel was not detected in polymorphic SNPs (Fig. 2). With C/T being the commonest (12) and C/G the least common (1) alteration observed. Transition/transversion ratio (ts/tv) was 1.910, SNP substitutions result of development marker shows transitions rate higher than transversions rate (Table 3). These results correspond with the EST resource SNP’s ones.

Table 4. Characteristics of the 32 polymorphic SNPs validated from the olive flounder EST database in this study
Primer name EST No. Gene identity E value
PO_02 98 Apolipoprotein AI precursor 2.00E-92
PO_03 38 Parvalbumin 4.00E-52
PO_23 23 Gelsolin, like 1 0.00E+00
PO_25 20 Ahsg protein 2.00E-41
PO_32 32 Trypsinogen 2 1.00E-132
PO_34 16 MHC class IIb antigen 1.00E-114
PO_36 15 Ferritin heavy subunit 4.00E-93
PO_44 13 Transferrin 1.00E-140
PO_46 13 V-Fos transformation effector 1.00E-140
PO_47 12 Family-2 cystatin 1.00E-27
PO_50 6 Testis enhanced gene transcript-like protein 2.00E-78
PO_52 11 Elastase 2 precursor 1.00E-148
PO_57 10 Elastase 1 precursor 1.00E-139
PO_59 9 Alpha-2-HS-glycoprotein 3.00E-60
PO_61 9 Apolipoprotein E1 5.00E-86
PO_62 9 Triosephosphateisomerase B 1.00E-127
PO_67 7 Alpha-1-antiproteinase-like protein 1.00E-129
PO_79 6 Glutathione-S-transferase 1.00E-107
PO_86 5 Ornithinedecarboxylaseantizyme large isoform 2.00E-86
PO_91 5 Troponin C, fast skeletal 2.00E-71
PO_92 4 Natural killer enhancing factor 1.00E-111
PO_109 3 Glyceraldehyde 3-phosphate dehydrogenase 3.00E-60
PO_110 3 Pyruvate kinase 4.00E-68
PO_116 3 Glycerol-3-phosphate dehydrogenase 1 (soluble) 1.00E-118
PO_121 69 Cytochrome c oxidase subunit I 0.00E+00
PO_122 65 Cytochrome c oxidase subunit II 1.00E-105
PO_126 7 Lactate dehydrogenase-A; NAD+:lactate oxidoreductase; LDH-A 1.00E-166
PO_130 7 Nephrosin 1.00E-109
PO_136 5 Proteasome subunit C10-11 1.00E-110
PO_172 2 Ovel protein similar to type I cytokeratin,enveloping layer 2.00E-71
PO_176 2 Ornithine decarboxylase antizyme small isoform 2.00E-74
PO_179 2 NADH dehydrogenase subunit 2 2.00E-64

PO_121 marker had the highest value of Heat 0.50 while Primer PO_126 marker had the lowest value of 0.12. The estimation of polymorphic information content (PIC) by marker ranged from 0.11 to 0.38 (Table 5).

Download Excel Table
Table 3. Summary statistics of SNP identification from new designed marker
Parametrs Detected SNPs
Total designed primer pairs 144
Successful assays 83
Failed assays 61
Detected SNP substitutions from sample 32
Undetected SNP substitutions from sample 51
Total number of transitions 21
C/T 12
A/G 9
Total number of transversions 11
A/T 2
A/C 2
T/G 6
C/G 1
Total number of indels 0
Transition/transversion 1.910
Download Excel Table
Fig. 2. Substitution types detected using the new designed marker in sample.
Download Origianl Figure
2. SNP genotype and allele frequency

The calculation on frequency of genotyping and allele of 32 SNP markers, there were differences in allele frequency between wild and cultured olive flounder. In particular, Primer PO_03 (Allele: TT, 100%), Primer PO_25 (Allele: TT, 100%), Primer PO_91 (Allele: GG, 100%), Primer PO_116 (Allele: CC, 100%) marker had 100 percent one allele type in wild sample while cultured sample had various types. Meanwhile, Primer PO_36 (Allele: TT, 100%), Primer PO_50 (Allele: CC, 100%), Primer PO_79 (Allele: CC, 100%), Primer PO_122 (Allele: AA, 100%), Primer PO_179 (Allele: AA, 100%) had one allele type in cultured sample while had various alleles in wild sample. The analysis on heterozygosity (He) by sample showed 0.34 in wild sample and 0.29 in cultured sample, respecttively. The range of heterozygosity by marker was identified from 0.12 to 0.50. Primer PO_61, Primer PO_62, Primer PO_121 marker had the highest value of Heat 0.50 while Primer PO_126 marker had the lowest value of 0.12. The estimation of polymorphic information content (PIC) by marker ranged from 0.11 to 0.38 (Table 5).

Table 5. Frequencies of genotype and allele for the SNPs in olive flounder
Download Excel Table


Genomic regions of immune and protein synthesisrelated genes is directly related to life of organisms so that its importance is well known (He et al., 2011; He et al., 2012). In addition, single base substitution of specific genomic region can be a positive or negative a factor for survival, growth and adaptation to the environment. Because these usefulness, the SNP research has been actively proceeding in various fish species; but hasn’t yet been reported on olive flounder. In our experiments, discovery of SNP based on EST database in olive flounder and developed a valuable marker for applications to various genomic analyses. The SNPs polymorphism show that a total of 21 transitions (C/T and A/G), and 11 transversions (A/T, A/C, T/G and C/G) while transitions are higher than transversions. This result is similar with other previous research (Vera et al., 2011). These result suggest that transitions are more easily substituted than transversion, because purine-purine (transitions) substitution is easy and stable than purine-pyrimidine (transversions) substitution.

Moreover, 32 SNPs genomic region of identified polymorphism are related to immune response and protein synthesis, functional changes in the these proteins causes the phenotypic change. The single nucleotide substitution of the specific region might influences on immunity and survival of organisms. In particular, glutathione S-transferase, nephrosin, elastases factors, MHC II alpha beta antigens and apolipoprotein involved in innate immunity (Dijkstra et al., 2013; Tsai ea al., 2004; Vera et al., 2011). In addition, ferritin has an important role for iron storage and detoxification, the storage of iron is closely related to the immune response to bacterial infection protected. Finally, transferrin could be related to the inflammation and enhanced resistance to oxidative stress by parasite infections (Paramá et al., 2007). Also, SNP detected genes in olive flounder show a high similarity with immune-related genes in mammals (Vera et al., 2011). The innate immune response is an important and highly developed defence mechanism against pathogens in fish. Examples of innate immunity include anatomic barriers, mechanical removal of pathogens, bacterial antagonism, pattern-recognition receptors, antigen-nonspecific defense compounds, complement pathways, phagocytosis, and inflammation fish (Ellis et al., 2001). The research outcome can be used as basis data for a study to identify the relationship between the nucleotide substitution in the immune-related genes specific region and immune response.

The allele frequency analysis using newly developed SNP marker shows that a variety of alleles is present in wild olive flounder samples, but does not shows the diversity in a few cultured samples. It shows the necessity to use different markers for the wild and cultured sample for genome analysis. Primer PO_36, Primer PO_50, Primer PO_79, Primer PO_122, Primer PO_179 marker are useful to detect SNP for wild olive flounder samples, while Primer PO_03, Primer PO_25, Primer PO_91, Primer PO_116 are suitable for cultured samples. In addition, this specific marker is important to distinguish between wild and culture olive flounder sample.

The diversity analysis using developed SNP marker in our research shows the heterozygosity value is higher in the wild sample of 0.34 than the cultured sample of 0.29, meaning higher diversity of the wild than the cultured. The cultured sample has got smaller diversity due to continuous inbreeding in the group without any inflow of new species (new mother group), resulting in lower diversity than the wild (Xing et al., 2014). Result of our heterozygosity is a correlation with diversity study outcome of olive flounder using microsatellite (Jung et al., 2009). It supports the appropriateness and usefulness of newly developed SNP marker to apply not only to SNP detection but also to other diversity analysis.

In conclusion, we was identified the substitution using the new marker from sample, it supports that development marker on this study is suitable for SNP detection in olive flounder. In addition, the marker is an appropriate tool for diversity analysis. The outcome of this study will be used in the selection of breeding with economic traits and will facilitate the molecular breeding programs of olive flounder.


This work was supported by a grant from the National Fisheries Research and Development Institute (NFRDI), contribution number RP-2014-BT-017.



Campbell NR, Narum SR. Development of 54 novel single-nucleotide polymorphism (SNP) assays for sockeye and coho salmon and assessment of available SNPs to differentiate stocks within the Columbia River. Mol Ecol Resour. 2011; 11:20-30.


Cenadelli S, Maran V, Bongioni G, Fusetti L, Parma P, Aleandri R. Identification of nuclear SNPs in gilthead seabream. J Fish Biol. 2007; 70:399-405.


Dijkstra JM, Grimholt U, Leong J, Koop BF, Hashimoto K. Comprehensive analysis of MHC class II genes in teleost fish genomes reveals dispensability of the peptide-loading DM system in a large part of vertebrates. BMC Evol Biol. 2013; 260:26-13.


Ellis AE. Innate host defence mechanisms of fish against viruses and bacteria. Dev Comp Immunol. 2001; 25:827-839.


Everett MV, Grau ED, Seeb JE. Short reads and non-model species exploring the complexities of next generation sequence assembly and SNP discovery in the absence of a reference genome. Mol Ecol Resour. 2011; 11:93-108.


Fahrenkrug SC, Freking BA, Smith TP, Rohrer GA, Keele JW. Single nucleotide polymorphism (SNP) discovery in porcine expressed genes. Anim Genet. 2002; 33:186-195.


Hayes B, Laerdahl JK, Lien S, Moen T, Berg P, Hindar K, Davidson WS, Koop BF, Adzhubei A and Hoyheim B. An extensive resource of single nucleotide polymorphism markers associated with Atlantic salmon (Salmo salar) expressed sequences. Aquaculture. 2007; 265:82-90.


Hayes B, Sonesson AK, Gjerde B. Evaluation of three strategies using DNA markers for traceability in aquaculture species. Aquaculture. 2005; 250:70-81.


He C, Chen L, Simmons M, Li P, Kim S, Liu ZJ. Putative SNP discovery in interspecific hybrids of catfish by comparative EST analysis. Anim Genet. 2003; 34:445-448.


He F, HS Wen SL Dong, B Shi CF Chen. Identification of single nucleotide polymorphism cytochrome P450-c19a and its relation to reproductive traits in Japanese flounder (Paralichthys olivaceus). Aquaculture. 2008; 279:177-181.


He F, Ma RQ, Wen HS, Li JF, Mu WJ, Liu M, Zhang YQ, Hu J, Qun L. Polymorphysims of CYP17-I gene in the exons were associated with the reproducetive endocrine of Japanese flounder (Paralichthys olivaceus). J Anim Sci. 2012; 6:794-799.


He F, Wen HS, Li JF, Yu DH, Ma RQ, Shi D, Mu WJ, Zhang YQ, Hu J, Liu M, Han GW, Zhang JN, Wang QQ, Yuan UR, Liu Q. Single nucleotide polymorphisms of the GnRHR gene associated with reproductive traits of Japanese flounder (Paralichthys olivaceus). J Anim Sci. 2011; 4:463-470.


Hohenlohe PA, Amish SJ, Catchen JM, Allendorf FW, Luikart G. Next-generation RAD sequencing identifies thousands of SNPs for assessing hybridization between rainbow and westslope cutthroat trout. Mol Ecol Resour. 2011; 11:117-122.


Hubert S, Bussey JT, Higgins B, Curtis BA, Bowman S. Development of single nucleotide polymorphism markers for Atlantic cod (Gadus morhua) using expressed sequences. Aquaculture. 2009; 296:7-14.


Jung DS, Noh JG, Myung JI, Lee JH, Kim HC, Park CJ, Min BH, Ha DS, Jeoun CY. Genetic variability comparison of wild populations and cultured stocks of flounder Paralichthys olivaceus based on microsatellited DNA Markers. Kor J Ichthyol. 2009; 21:221-226.


Liu ZJ, Cordes JF. DNA marker technologies and their applications in aquaculture genetics. Aquaculture. 2004; 238:1-37.


Paramá A, Castro R, Arranz JA, Sanmartín ML, Lamas J, Leiro J. Scuticociliate cysteine proteinases modulate turbot leucocyte functions. Fish Shellfish Immunol. 2007; 23:945-956.


Parsons JD. Improved tools for DNA comparison and clustering. Comput Appl Biosci. 1995; 11:603-613.


Rozen S, Skaletsky H. Primer3 on the WWW for general users and for biologist programmers. Methods in Mol Biol. 2000; 132:365-386.


The International SNP Map Working Group A map of human genome sequence variation containing 142 million single nucleotide polymorphisms. Nature. 2001; 409:928-933.


Tsai PL, Chen CH, Huang CJ, Chou CM, Chang GD. Purification and cloning of an endogenous protein inhibitor of carp nephrosin, an astacin metalloproteinase. J Biol Chem. 2004; 279:19-32.


Vera M, Antonio J, Dios A, Millán A, Pardo AG, Bouza C, Hermida M, Fernández C Herr n, Molina-Luzón M, Paulino M. Validation of single nucleotide polymorphism (SNP) markers from an immune Expressed Sequence Tag (EST) turbot, Scophthalmus maximus. Aquaculture. 2011; 313:31-41.


Vignal A, Milan D, San Cristobal M, Eggen A. A review on SNP and other types of molecular markers and their use in animal genetics. Genet Sel Evol. 2002; 34:275-305.


Xing K, Gao ML, Li HJ. Genetic differentiation between natural and hatchery populations of Manila clam (Ruditapes philippinarum) based on microsatellite markers. Genet Mol Res. 2014; 13:237-245.


Yoon JM. Genetic differences and variations in freshwater crab (Eriocheir sinensis) and swimming crab (Portunus trituberculatus). Dev Repord. 2006; 6:19-32.