Advanced
Cytochrome oxidase subunit I (COI) DNA sequence divergence between two cryptic species of Oryzias in South Korea
Cytochrome oxidase subunit I (COI) DNA sequence divergence between two cryptic species of Oryzias in South Korea
Journal of Ecology and Environment. 2013. Sep, 36(3): 159-166
Copyright ©2013, The Ecological Society of Korea
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 : March 11, 2013
  • Accepted : July 12, 2013
  • Published : September 27, 2013
Download
PDF
e-PUB
PubReader
PPT
Export by style
Article
Author
Metrics
Cited by
TagCloud
About the Authors
Dong-Su In
Biological Resource Center, Kongju National University, Gongju, 314-701, Korea
Eun-Sook Choi
Department of Biology Education, Kongju National University, Gongju, 314-701, Korea
Ju-Duk Yoon
Biological Resource Center, Kongju National University, Gongju, 314-701, Korea
Jeong-Hui Kim
Department of Biology Education, Kongju National University, Gongju, 314-701, Korea
Jun-Il Min
Department of Animation, Kongju National University, Gongju, 314-701, Korea
Seung-Ho Baek
Department of Biology Education, Kongju National University, Gongju, 314-701, Korea
Min-Ho Jang
Department of Biology Education, Kongju National University, Gongju, 314-701, Korea
jangmino@kongju.ac.kr

Abstract
Oryzias latipes and Oryzias sinensis are indigenous species found in Japan, China, and other East Asian countries, including Korea. Based on morphological differences, the species have been classified distinctly. However, the range of morphological characters such as the number of gill rakers, vertebrae, and spots on the lateral body overlaps and is too vague for clear identification, so their classification based on their morphological characteristics remains uncertain. In this study, the mitochondrial cytochrome oxidase subunit I (COI) gene, which is used for DNA barcoding, was applied to clarify interspecific variation of O. latipes and O. sinensis . Intraspecific genetic diversity was calculated to identify correlations with geographic distributions. We studied two species collected from 55 locations in Korea. All individuals carried a 679-base pair gene without deletion or insertion. Between species, 525 base pairs of the gene were shared. The Kimura two parameter (K2P) distance of O. latipes and O. sinensis was 0.41% and 1.39%, respectively. Mean divergence within genera was 23.5%. Therefore, the species were clearly different. The distance between O. latipes and O. sinensis was 14.0%, which is the closest within genera. Interestingly O. latipes from the Japanese and Korean group represented 16.5% distant. These results were derived from geohistorical and anthropogenic environmental factors. The O. latipes haplotypes were joined in only one group, but O. sinensis was divided into two groups, one is found in the Han River and upper Geum River watershed; the other is found in the remaining South Korean watersheds. Further studies will address the causes for geographic speciation of O. sinensis haplotypes.
Keywords
INTRODUCTION
Morphological identification, the most common method for classification, has been employed for many years. Studies involve measuring and comparing the phenotypes of target species, sometimes by microscopic examination. This method is normal and sound when it applies to well-known taxa; however, it is not always efficient for species identification, unless the phenotype is obvious. To overcome this problem, a taxonomic system based on DNA analysis was established in the early 2000s (Tautz et al. 2002, 2003).
Mitochondrial DNA cytochrome oxidase subunit I (mtDNA COI) has been proposed for used in animal bioidentification because the COI gene could differentiate nearly all animal species with only a single gene sequence (Hebert et al. 2003, 2005). In the early 2000s, various taxonomic studies were carried out using the COI gene
Lager Image
Oryzias collection sites in South Korea.
(Hebert et al. 2003, 2004a, 2004b, 2005, Smith et al. 2006, 2007, Bucklin et al. 2010). The universal primer for the COI gene is very stable and attaches at the 5′ end of the animal gene (Folmer et al. 1994, Zhang and Hewitt 1997). The COI provides profound insights into phylogenetic studies more than other mitochondrial genes such as cytochrome b (Simmons and Weller 2001), because changes in the amino acid sequence occur more slowly (Lynch and Jarrell 1993). Moreover, mtDNA extraction and sequencing are easy because the small multi-copy mitochondrial genome is located in each cell. Inter-specific variation is high and intra-specific variation is low, making species classification clears (Jin 2012).
Oryzias latipes and Oryzias sinensis (genus Oryzias ) are indigenous species distributed in South Korea, China, Japan, and East Asia. Both species are shoaling fish that stay near the water surface and favor static, shallow waters such as reservoirs, swamps, and small streams (Kim and Park 2002). They were classified by morphological analysis as different species in the early 1990s; the China-western Korea group is designated O. sinensis and the eastern Korea group is O. latipes (Kim and Lee 1992, Kim and Kim 1993). However, the range of classification keys, number of gill rakers, vertebrae, and spots on the lateral body overlap and are too vague for clear identification, so debate regarding their classification continues.
In this study, the mitochondrial COI region used for DNA barcoding was applied to clarify the interspecific variation of O. latipes and O. sinensis . Intraspecific genetic diversity was also calculated to identify the correlation between geographic distributions.
MATERIALS AND MEHTODS
- Sample collection and identification of genus Oryzias
Sampling was conducted from 2009 to 2011 and sites were selected based on literature and government reports (Appendix 1). Other sites were chosen based on the possibility of Oryzias habitation, even in the absence of records of their presence. Samples were collected at 55 sites, most of which were located in the mid-lower parts of streams ( Fig. 1 , Appendix 1). At each site, three specimens were chosen for preservation in 99.9% ethanol; other specimens were released at the site of capture. After sample moved to the laboratory, each specimen was identified, labeled for individual identification, and measured total length (mm), standard length (mm) and body weight (g). Spots on the lateral body were used for identification (Park et al. 2006); if spots were distributed on the caudal peduncle, the specimen was classified as O. latipes and if not, it was classified as O. sinensis .
- DNA analysis
We analyzed 165 specimens from 55 sites. DNA extracts were prepared from muscle tissue (about 30 mg) using SolGent™ Genomic DNA Prep Kit (Solgent, Daejoen, South Korea). To amplify the COI gene from mitochondrial DNA, PCR was performed with primers FishF2 (5′-TCAACYAATCAYAAAGATAT-3′) and FishR2 (5′-ACTTCYGGGTGRCCRAARAA-3′) (Ward et al., 2005). The 50 μl PCR mixture included 5 μl 10× reaction buffer, 1 μl each primer (10 pmole/μl), 4 μl dNTP (2.2 mM each), 0.25 μl TaKaRa Ex Taq™ (Takara, Shiga, Japan) and 4 μl template DNA (50 ng/μl). Amplifications were performed using a C1000 gradient thermal cycler (Bio-RAD, Foster City, CA, USA). Cycling included an initial step of 3 min at 95℃ followed by 40 cycles of 0.5 min at 95℃, 0.5 min at 55℃, and 1 min at 72℃, followed by elongation for 7 min at 72℃. PCR products were analyzed by 1% agarose gel electrophoresis, and then purified with a Solg ™ PCR Purification Kit (SolGent, Daejoen, South Korea). Bidirectional sequencing was performed with the FishF2 and FishR2 primers and an 3730xl DNA Analyzer (Applied Biosystems ® , Foster City, CA, USA) with the BigDye ® Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems ® ). Contigs of DNA sequences were generated in CodonCode aligner ver. 3.7.1 (CodonCode Corporation, Dedham, MA, USA). Sequences were aligned with BioEdit ver. 7.0.0 by using the ClustalW Multiple alignment method. Aligned sequences of genus Oryzias and references (6 other Oryzias species and Hyporhamphus sajori as the outgroup) from GenBank (gi208341979; gi253960407; gi314909980; gi253960477; gi253960463; gi253960435; gi253960393; gi254939398; gi253960421) were compared using MEGA5 (Tamura et al. 2011). Neighbor-joining (NJ) was used for phylogenetic analysis. To assess statistical reliability, an NJ tree was generated by bootstrapping with 1000 replications (Felsenstein 1993) and Kimura two parameter (K2P) distances were determined to provide a graphic representation of the patterning of divergence between species. Sequence divergences were calculated using the K2P distance model (Kimura 1980). All Oryzias sequences include references calculated with NETWORK ver. 4.6.11 using the median-joining network option for intraspecific relationship (Bandelt et al. 1999) to identify relationships.
RESULTS
- Analysis of the mtDNA COI gene
The COI sequence length was 679-bp and there were no insertions or deletions in any specimen. The COI gene in O. latipes included 629 conserved bp and 23 variable bp. Base compositions were: A, 22.8%; T, 19.5%; G, 30.9%; C, 26.8%. In O. sinensis , the COI gene included 620 conserved and 32 variable bp. Base composition was slightly different in O. latipes : A, 23.7%; T, 17.9%; G, 31.9%; C, 26.5%. The sequence variation in O. sinensis was relatively greater than in O. latipes . The overlapping sequence in both species was 525-bp long.
The K2P distance of O. latipes was 0.41%, and the distance for O. sinensis was 3 times larger, at 1.39%. Mean divergence within genera was 23.5%. Therefore, the species were obviously different. The distance of O. latipes and O. sinensis was 14.0%, which is the closest within genera. Conversely, O. latipes and O. marmoratus showed 26.5% distance, so the farthest within genera. Moreover, interestingly O. latipes from Japanese group and Korean group represented 16.5% of distance ( Table 1 ).
- Molecular phylogeny of genus Oryzias
The molecular phylogeny results indicate O. latipes and O. sinensis in South Korea were clearly branched under 100 bootstrapping (BP) value. O. latipes formed a single group including the GenBank reference gi208341979. The O. sinensis group with GenBank reference gi314909980 was divided into two groups. The BP value of the O. sinensis groups was 99. O. latipes in Japan comprise a different branch with a BP value < 89 ( Fig. 2 ). A total of 19 O. latipes haplotypes from 30 sites were identified. The estimated number of mutations in the shortest tree within torso was 3, and the estimated number of mutations in the shortest tree was 25. There was only one network of O. latipes in South Korea. In contrast, 14 O. sinensis haplotypes were
Lager Image
Molecular phylogeny of Oryzias species using the Neighbor-joining method in MEGA5.
Matrix of K2P (Kimura 2 Parameter) for 9Oryziasspecies.Hyporhamphus sajoriis an outgroup
Lager Image
Matrix of K2P (Kimura 2 Parameter) for 9 Oryzias species. Hyporhamphus sajori is an outgroup
Lager Image
Intraspecific relationship between haplotypes in Korea constructed by median-joining networks with Oryzias latipes and O. sinensis. Circles indicateeach haplotype and geographical location; bold lines represent mutated regions.
identified across 25 sites; the estimated number of mutations in the shortest tree was 32. Thus, the O. sinensis group is divided into two subclades by the Charyeong mountain chain; subclade A is from the Han River watershed to the northern part of the Geum River watershed and subclade B is from the southern part of the Geum River watershed to the Yeongsan and Nakdong River watershed ( Fig. 3 ).
DISCUSSION
- COI region for species identification for Oryzias species
O. latipes and O. sinensis are considered different species; they have a different number of chromosomes (Park et al. 2006) and are reproductively incompatible (Kim and Kim 1993). Recent results from mitochondrial DNA cytochrome b and allozyme confirmed the biogeographic and reproductive isolation of these species (Sakaizumi and Jeon 1987, Matsuda et al. 1997, Takehana et al. 2004). These results are consistent with karyotype analyses ( O. latipes : eastern Korea grout, 2n = 48 and O. sinensis : China- western Korea group, 2n = 46) (Uwa 1986, Kim and Moon 1987, Uwa and Jeon 1987). Haplotype analysis of the mtDNA control region (Kang et al. 2005) revealed a similar distribution of O. sinensis .
In this study, we clearly classified two species using the mitochondrial COI region. Thus, the COI region could be used to supplement current taxonomic keys that are largely focused on morphological characteristics. In many cases, morphological classification generates identification errors due to vague criteria. Therefore, species classification using DNA markers such as the COI region (DNA barcodes) is an appropriate alternative.
NJ trees clearly differentiated O. latipes and O. sinensis independently. The COI sequences provided strong support for the monophyletic origin of O. latipes and O. sinensis . BP of each species was the same for 100. Thus, the species are completely different. The K2P distance of O. sinensis was higher than that of O. latipes , indicating the intraspecific variation of O. sinensis is greater. Two clusters occurred within O. sinensis , and K2P distance was 2.17%. We suspect there is more variation within O. sinensis . The early study conducted by Min (1997) indicated a similar result, suggesting the genetic variation of O. sinensis groups was about 2.5 times higher than that of O. latipes . Genetic variations are normally accumulated through many factors such as stability, group and habitat size, etc. (Soule 1972, 1976, Valentine 1976, Stonecking et al. 1981). Thus, we suspect differences in habitat or environmental conditions exist.
- Geographic distribution of genus Oryzias
Oryzias latipes is also distributed in China and Japan, inhabiting streams that flow into the eastern part of Korea and the islands part of southwestern Korea. In contrast, O. sinensis is distributed in the streams and islands of western Korea and China (Kim and Moon 1987, Kim and Lee 1992). Through allozyme studies of O. latipes in Japan, northern and southern groups were distinguished and a relationship was identified between the China and southern groups (Sakaizumi 1986). An isozyme study indicated the Korean and Japanese O. latipes groups are genetically different species (Min 1997). Such distribution could be generated by geohistorical and environmental factors. O. sinensis flowed into western Korea from southern China when the Korean peninsula was connected via the paleo- Hwangho (Yellow River) during the regression period of the Pliocene. O. latipes originated from the paleo-Amur River (Kim et al. 2005). Geographic characteristics of the Korean peninsula also support the distinct distribution of genus Oryzias . The Taebaek Mountains divide South Korea vertically into eastern and western regions, serving as a physical barrier. Stream morphology and fish species significantly differ between these regions. Thus, although similar species flowed from different sources into the Korean peninsula, they have maintained their original genetic identities because of the geographic characteristics of the Korean peninsula. In this study, O. latipes intruded into the western part of Korea. We believe this is due to anthropogenic causes, rather than natural displacement. Unintentional human releases are a significant cause for fish displacement. Also this newly found O. sinensis inhabitation sites (N3, YS 12 and YS 13). As we demonstrated, O. latipes from Japan and Korea are incorporated into different clades and the K2P distance was 16.50%. Thus, the phylogeny of these groups should be reconsidered for classification. In addition, overall consideration of the genetic relationship between species in Korea, China, and Japan is required for establishing genetic divergence in Asian countries.
Collection sites information of genusOryzias(H: Han River Watershed, G: Geum River Watershed, YS: Yeongsan-Seomjin River Watershed, N: Nakdong River Watershed)
Lager Image
Collection sites information of genus Oryzias (H: Han River Watershed, G: Geum River Watershed, YS: Yeongsan-Seomjin River Watershed, N: Nakdong River Watershed)
Acknowledgements
This work was supported by the research grant of the Kongju National University in 2012 (Project No.: 2011-0896).
References
Bandelt HJ , Forster P , Rohl A 1999 Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol 16 37 - 48    DOI : 10.1093/oxfordjournals.molbev.a026036
Bucklin A , Hopcroft RR , Kosobokova KN , Nigro LM , Ortman BD , Jennings RM , Sweetman CJ 2010 DNA barcoding of Arctic Ocean holozooplankton for species identification and recognition. Deep-Sea Res Pt II 57 40 - 48    DOI : 10.1016/j.dsr2.2009.08.005
Felsenstein J , Kishino H 1993 Is there something wrong with the bootstrap on phylogeny? A reply to Hills and Bull. Syst Biol 42 193 - 200
Folmer O , Black M , Hoeh W , Lutz R , Vrijenhoek R 1994 DNA primers for amplification of mitochondrial cytochrg ome c oxidase subunit I from diverse metazoan invertebrates. Mol Mar Biol Biotech 3 294 - 299
Hebert PDN , Cywinska A , Ball SL , deWaard JR 2003 Biological identifications through DNA barcodes. Proc R Soc Lond B Bio Sci 270 313 - 322    DOI : 10.1098/rspb.2002.2218
Hebert PDN , Gregory TR 2005 The promise of DNA barcoding for taxonomy. Syst Biol 54 852 - 859    DOI : 10.1080/10635150500354886
Hebert PDN , Penton EH , Burns JM , Janzen DH , Hallwachs W 2004a Ten species in one: DNA barcoding reveals cryptic species in the neotropical skipper butterfly Astraptes fulgerator. Proc Natl Acad Sci USA 101 14812 - 14817    DOI : 10.1073/pnas.0406166101
Hebert PDN , Stoeckle MY , Zemlak TS , Francis CM 2004b Identification of birds through DNA barcodes. PLoS Biol 2 1657 - 1663    DOI : 10.1371/journal.pbio.0020312
Jin SD , PhD Dissertation 2012 Phylogenetic analysis of birds in korea using mtDNA variation. Chungnam National University Daejeon, South Korea PhD Dissertation
Kang TW , Lee EH , Kim MS , Paik SG , Kim S , Kim CB 2005 Molecular phylogeny and geography of Korean medaka fish (Oryzias latipes). Mol Cells 20 151 - 156
Kim IS , Choi Y , Lee CR , Lee YJ , Kim BJ , Kim JH 2005 The Great Illustrated Fish Book in Korea. Kyohak Publishing Seoul (in Korean)
Kim IS , Lee EH 1992 New record of ricefish, Oryzias latipes sinensis (Pisces, Oryzidae) from Korea. Korean J Syst Zool (in Korean) 8 177 - 182
Kim IS , Moon KC 1987 The karyotype of a ricefish, Oryzias latipes from southern Korea. Kor J Zool (in Korean with English abstract) 30 379 - 386
Kim IS , Kim SY 1993 Interspecific hybridization between Oryzias sinensis and O. latipes from Korea. Kor J Ichthyol (in Korean with English abstract) 5 113 - 121
Kim IS , Park JY 2002 Freshwater Fishes of Korea. Kyo-hak Publishing Seoul (in Korean)
Kimura M 1980 A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol 16 111 - 120    DOI : 10.1007/BF01731581
Lynch M , Jarrell PE 1993 A method for calibrating molecular clocks and its application to animal mitochondrial DNA. Genetics 135 1197 - 1208
Min MS 1997 Genetic variation and speciation of 2 species of genus Oryzias (Pisces, Adrianichthyidae) in Korea. Korean J Syst Zool (in Korean) 13 9 - 20
Matsuda M , Yamagishi M , Sakizumi M , Jeon SR 1997 Mitochondrial DNA variation in the Korean wild population of medaka, Oryzias latipes. Kor J Limnol 30 119 - 128
Park MK , Kim IS , Park JY 2006 Morphological characters and geographical distributions of Korean ricefish (Pisces, Adrianichthyidae), Oryzias. Kor J Ichthyol 18 163 - 169
Sakaizumi M 1986 Genetic divergence in wild populations of medaka Oryzias latipes (Pisces: Oryziatidae) from Japan and China. Genetica 69 119 - 125    DOI : 10.1007/BF00115131
Sakaizumi M , Jeon SR 1987 Two divergent groups in the wild populations of medaka Oryzias latipes (Pisces: Oryziatidae) in Korea. Kor J Limnol 20 13 - 20
Simmons RB , Weller SJ 2001 Utility and evolution of cytochrome b in insects. Mol Phylogenet Evol 20 196 - 210    DOI : 10.1006/mpev.2001.0958
Smith MA , Wood DM , Janzen DH , Hallwachs W , Hebert PDN 2007 DNA barcodes affirm that 16 species of apparently generalist tropical parasitoid flies (Diptera, Tachinidae) are not all generalists. Proc Natl Acad Sci USA 104 4967 - 4972    DOI : 10.1073/pnas.0700050104
Smith MA , Woodley NE , Janzen DH , Hallwachs W , Hebert PDN 2006 DNA barcodes reveal cryptic host-specificity within the presumed polyphagous members of a genus of parasitoid flies (Diptera: Tachinidae). Proc Natl Acad Sci USA 103 3657 - 3662    DOI : 10.1073/pnas.0511318103
Soule M 1972 variation problem: the gene flow variation hypothesis. Taxon 20 37 - 50    DOI : 10.2307/1218532
Soule M 1976 Allozyme variation: its determinants in a space and time. In: Molecular Evolution (Ayala FJ, ed). Sinauer Associates Sunderland 60 - 77
Stoneking M , Wagner DJ , Hildebrand AC 1981 Genetic evidence suggesting subspecific differences between northern and southern populations of brook trout (Salvelinus fontinalis). Copeia 4 810 - 819    DOI : 10.2307/1444182
Takehana Y , Jeon SR , Sakaizumi M 2004 Genetic structure of Korean wild populations of the medaka Oryzias latipes inferred from Allozymic variation. Zool Sci 21 977 - 988    DOI : 10.2108/zsj.21.977
Tamura K , Peterson N , Peterson D , Stecher G , Nei M , Kumar S 2011 MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol 28 2731 - 2739    DOI : 10.1093/molbev/msr121
Tautz D , Arctander P , Minelli A , Thomas RH , Vogler AP 2002 DNA points the way ahead in taxonomy. Nature 418 479 -
Tautz D , Arctander P , Minelli A , Thomas RH , Vogler AP 2003 A plea for DNA taxonomy. Trends Ecol Evol 18 70 - 74    DOI : 10.1016/S0169-5347(02)00041-1
Uwa H 1986 Karyotype evolution and geographical distribution in the ricefish, genus Oryzias (Oryziidae). In: Indo-Pacific Fish Biology: Proceedings of the Second International Conference on Indo-Pacific Fishes (Uyeno T, Arai T, Taniuchi T, Matsuura K, eds). Ichthyological Society of Japan Tokyo 867 - 876
Uwa H , Jeon SR 1987 Karyotypes in two divergent groups of a ricefish, Oryzias latipes, from Korea. Kor J Limnol 20 139 - 147
Valentine JW 1976 Genetic strategies of adaptation. In: Molecular Evolution (Ayala FJ, ed). Sinauer Associates Sunderland, MA 78 - 94
Ward RD , Zemlak TS , Innes BH , Last PR , Hebert PDN 2005 DNA barcoding Australia’s fish species. Philos T Roy Soc B 360 1847 - 1857    DOI : 10.1098/rstb.2005.1716
Zhang DX , Hewitt GM 1997 Assessment of the universality and utility of a set of conserved mitochondrial COI primers in insects. Insect Mol Biol 6 143 - 150    DOI : 10.1111/j.1365-2583.1997.tb00082.x