Cytochrome oxidase subunit I (COI) DNA sequence divergence between tow cryptic species of Oryzias in South Korea
- DOI : 10.5141/ecoenv.2013.159
- Author: In Dong-Su, Choi Eun-Sook, Yoon Ju-Duk, Kim Jeong-Hui, Min Jun-Il, Baek Seung-Ho, Jang Min-Ho
- Organization: In Dong-Su; Choi Eun-Sook; Yoon Ju-Duk; Kim Jeong-Hui; Min Jun-Il; Baek Seung-Ho; Jang Min-Ho
- Publish: Journal of Ecology and Environment Volume 36, Issue3, p159~166, 27 Sep 2013
-
ABSTRACT
Oryzias latipes andOryzias 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 ofO. latipes andO. 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 ofO. latipes andO. sinensis was 0.41% and 1.39%, respectively. Mean divergence within genera was 23.5%. Therefore, the species were clearly different. The distance betweenO. latipes andO. sinensis was 14.0%, which is the closest within genera. InterestinglyO. latipes from the Japanese and Korean group represented 16.5% distant. These results were derived from geohistorical and anthropogenic environmental factors. TheO. latipes haplotypes were joined in only one group, butO. 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 ofO. sinensis haplotypes.
-
KEYWORD
cytochrome oxidase subunit I , geographic distribution , K2P distance , Oryzias latipes , Oryzias sinensis
-
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
(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 andOryzias sinensis (genusOryzias ) 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 designatedO. sinensis and the eastern Korea group isO. 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 andO. sinensis . Intraspecific genetic diversity was also calculated to identify the correlation between geographic distributions.> 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 asO. latipes and if not, it was classified asO. sinensis .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 genusOryzias and references (6 otherOryzias species andHyporhamphus 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). AllOryzias 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.> 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%. InO. sinensis , the COI gene included 620 conserved and 32 variable bp. Base composition was slightly different inO. latipes : A, 23.7%; T, 17.9%; G, 31.9%; C, 26.5%. The sequence variation inO. sinensis was relatively greater than inO. latipes . The overlapping sequence in both species was 525-bp long.The K2P distance of
O. latipes was 0.41%, and the distance forO. sinensis was 3 times larger, at 1.39%. Mean divergence within genera was 23.5%. Therefore, the species were obviously different. The distance ofO. latipes andO. sinensis was 14.0%, which is the closest within genera. Conversely,O. latipes andO. marmoratus showed 26.5% distance, so the farthest within genera. Moreover, interestinglyO. 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 andO. sinensis in South Korea were clearly branched under 100 bootstrapping (BP) value.O. latipes formed a single group including the GenBank reference gi208341979. TheO. sinensis group with GenBank reference gi314909980 was divided into two groups. The BP value of theO. sinensis groups was 99.O. latipes in Japan comprise a different branch with a BP value < 89 (Fig. 2). A total of 19O. 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 ofO. latipes in South Korea. In contrast, 14O. sinensis haplotypes wereidentified 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).> COI region for species identification for Oryzias species
O. latipes andO. 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 andO. 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 ofO. 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 andO. sinensis independently. The COI sequences provided strong support for the monophyletic origin ofO. latipes andO. sinensis . BP of each species was the same for 100. Thus, the species are completely different. The K2P distance ofO. sinensis was higher than that ofO. latipes , indicating the intraspecific variation ofO. sinensis is greater. Two clusters occurred withinO. sinensis , and K2P distance was 2.17%. We suspect there is more variation withinO. sinensis . The early study conducted by Min (1997) indicated a similar result, suggesting the genetic variation ofO. sinensis groups was about 2.5 times higher than that ofO. 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 ofO. 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 JapaneseO. 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 genusOryzias . 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 foundO. 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.-
[Fig. 1.] Oryzias collection sites in South Korea.
-
[Fig. 2.] Molecular phylogeny of Oryzias species using the Neighbor-joining method in MEGA5.
-
[Table 1.] Matrix of K2P (Kimura 2 Parameter) for 9 Oryzias species. Hyporhamphus sajori is an outgroup
-
[Fig. 3.] Intraspecific relationship between haplotypes in Korea constructed by median-joining networks with Oryzias latipes and O. sinensis. Circles indicate each haplotype and geographical location; bold lines represent mutated regions.
-
[Appendix 1_1] Collection sites information of genus Oryzias (H: Han River Watershed, G: Geum River Watershed, YS: Yeongsan-Seomjin River Watershed, N: Nakdong River Watershed)