A Natural Hybrid of Intergeneric Mating between a Female Pungtungia herzi and a Male Pseudorasbora parva (Cypriniformes: Cyprinidae)
A Natural Hybrid of Intergeneric Mating between a Female Pungtungia herzi and a Male Pseudorasbora parva (Cypriniformes: Cyprinidae)
Fisheries and aquatic sciences. 2015. Mar, 18(1): 99-107
Copyright © 2015, The Korean Society of Fisheries and Aquatic Science
This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial Licens ( permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
  • Received : September 11, 2014
  • Accepted : November 17, 2014
  • Published : March 30, 2015
Export by style
Cited by
About the Authors
Keun-Yong, Kim
Marine Eco-technology Institute, Co., Ltd, Busan 608-830, Korea
Myeong-Hun, Ko
Division of EcoScience, Ewha Womans University, Seoul 120-750, Korea
Sung Jang, Cho
Boryeong Freshwater Ecology Center, Boryeong 355-923, Korea
Woo-Jin, Kim
Biotechnology Research Center, National Fisheries Research and Development Institute, Busan 619-902, Korea
Min Ho, Son
Marine Eco-technology Institute, Co., Ltd, Busan 608-830, Korea
In-Chul, Bang
Department of Life Science and Biotechnology, Soonchunhyang University, Asan 336-745, Korea

A natural hybrid of a probable intergeneric mating between the striped shiner Pungtungia herzi and the stone morocco Pseudorasbora parva (Cypriniformes: Cyprinidae) was captured in the Geumho River, a tributary of the Nakdong River basin in Korea. Morphological characters and DNA sequences were analyzed to verify its hybrid state and identify the parentage of its parent species. The hybrid exhibited a phenotypic intermediacy between the two parent species in the number of vertebrae and the mouth shape. Out of 1,488 base pair (bp) positions of the nuclear recombination activating gene 1 gene ( rag1 ), which has a biparental mode of inheritance, 41-bp substitutions were detected between the two parent species, whereas an electropherogram of the hybrid displayed polymorphic double peaks at all of the base positions, along with one additional one, strongly indicating its hybrid state. Meanwhile, sequence comparison of the mitochondrial cytochrome b gene ( mt-cyb ) (1,140 bp), which has a maternal mode of inheritance, showed only 5–22-bp differences (97.6–99.5% identities) between the hybrid and Pu. herzi , but as many as 158–168-bp differences (85.2–86.1% identities) between the hybrid and Ps. parva , clearly indicating Pu. herzi as the maternal species. Thus, combined analyses of independent data sets (i.e., morphology and nuclear and mitochondrial DNA sequences) offered convincing evidence for the hybrid state of a naturally occurring hybrid resulting from intergeneric mating between a female Pu. herzi and a male Ps. parva .
As in plants, hybridization and introgression play important roles in the evolution of fish ( Hubbs, 1955 ; Dowling and DeMarais, 1993 ; Scribner et al., 2001 ). They represent the driving force behind the evolution of new lineages with the potential to rapidly adapt to new habitats ( DeMarais et al., 1992 ; Gerber et al., 2001 ; Grant and Grant, 1992 ). However, hybridization can also lead to extinction of rare indigenous species and presents challenges for the formulation of a biological species concept and definition, the reconstruction of phylogenetic trees and the practice for biological conservation ( Grant and Grant, 1992 ; Rhymer and Simberloff, 1996 ; Allendorf et al., 2001 ). Natural hybridization occurs across diverse piscine taxa ( Verspoor and Hammar, 1991 ) and is quite common in cyprinid fishes ( Hubbs, 1955 ; Gerber et al., 2001 ; Scribner et al., 2001 ; Ünver et al., 2008 ). Hybridization between fish species frequently occurs in nature resulting from the intrinsic properties of the species, external fertilization and weakly developed reproductive isolating mechanisms, as well as natural or man-made causes, such as disparity in the number of parent species, habitat disturbance and fragmentation, the introduction of non-native species, and climate change ( Hubbs, 1955 ; Rhymer and Simberloff, 1996 ; Allendorf et al., 2001 ; Scribner et al., 2001 ; Šorić, 2004 ; Ünver et al., 2008 ; Kwan et al., 2014 ).
Morphological analysis is a conventional approach to identification of naturally occurring hybrids because of the intermediate phenotypes. However, morphology alone is not always sufficient to determine the hybrid status and to trace the parent species ( Avise and Saunders, 1984 ; Allendorf et al., 2001 ; Scribner et al., 2001 ). Analyses of genetic markers with different evolutionary histories provide convincing evidence for hybrid status and can help determine the extent and direction of hybridization (e.g., Avise and Saunders, 1984 ; Dowling et al., 1989 ; Avise et al., 1990 ; Albert et al., 2006 ; Chelomina et al., 2008 ). In practice, the combined application of nuclear DNA (nDNA) and mitochondrial DNA (mtDNA) sequences has proven useful in determination of hybridization in fish species (e.g., Sonnenberg et al., 2007 ; Lee et al., 2009 ; Yoon et al., 2009 ). These techniques are simple and easy to implement as the result of recent advances in molecular techniques, such as polymerase chain reaction (PCR), DNA sequencing, etc.
The mitochondrial genome is haploid, located in the cytoplasm and is maternally inherited without recombination. Given these distinct properties, analyses of mtDNA sequences provide an opportunity to determine the direction of gene flow in a hybrid cross ( Avise and Saunders, 1984 ; Moritz et al., 1987 ; DeMarais et al., 1992 ; Rhymer and Simberloff, 1996 ). However, the information provided by mtDNA alone is limited due to its strict maternal mode of inheritance. Thus, to determine the directionality of species hybridization, mtDNA analysis is often used in conjunction with analysis of a nuclear gene, which is inherited biparentally or in a simple Mendelian fashion. Hybrid state is reliably identifiable in monomorphic (homozygous) versus polymorphic (heterozygous) states; i.e., each pure parent species possesses a clear DNA sequence with monomorphic base positions, whereas their hybrids exhibit polymorphisms different to the parent ( Scribner et al., 2001 ; Sonnenberg et al., 2007 ).
The striped shiner Pungtungia herzi Herzenstein and the stone morocco (or the false dace) Pseudorasbora parva Temminck and Schlegel (Cypriniformes: Cyprinidae) are widely distributed in East Asia and coexist in the river basins of Korea ( Kim and Park, 2002 ; Kim et al., 2005 ). A natural hybrid of a probable mating between Pu. herzi and Ps. parva was captured during a survey of ichthyofauna in Korea. Morphological characters as well as nDNA and mtDNA sequences were analyzed to verify the hybrid state and determine the parentage.
Materials and Methods
- Collection of fish specimens
A single natural hybrid ( n = 1; standard length: 64.7 mm) between Pu. herzi and Ps. parva was captured with a spoon net (mesh size: 4 × 4 mm) on 21 April 2010 in the Geumho River, a tributary of the Nakdong River basin in Korea ( Fig. 1 ). Specimens of its supposed parent species, Pu. herzi ( n = 16) and Ps. parva ( n = 11) were also sampled from the same tributary. The specimens were deposited in the freshwater fish collection of the Department of Marine Biotechnology of Soonchunhyang University (SUC; Asan, Korea). Detailed sampling information is provided in Table 1 .
PPT Slide
Lager Image
Specimens of a natural hybrid and its supposed parent species, Pungtungia herzi and Pseudorasbora parva (left) and fragments of their electropherograms showing the recombination activating gene 1 gene (rag1) (right). A) Pu. herzi (♀), B) the natural hybrid and C) Ps. parva (♂). The fragments of the electropherograms show differences in base positions for the two parent species and double peaks at the respective positions for the natural hybrid (indicated by inverted arrowheads).
Information on a natural hybrid and its supposed parent species,Pungtungia herziandPseudorasbora parva
PPT Slide
Lager Image
Information on a natural hybrid and its supposed parent species, Pungtungia herzi and Pseudorasbora parva
- Morphological measurements
A total of 25 morphological characters were used in the calculation of the hybrid index (Hindex) to measure the degree of phenotypic intermediacy of the hybrid in comparison to its parent species. The number of fin rays and vertebrae were counted using a soft X-ray photographer (Hitex HA-80; Hitex, Osaka, Japan). Total of six meristic and 19 morphometric characters were counted or measured as described by Hubbs and Lagler (1964) ( Table 2 ). The mean of Pu. herzi was set to zero, and the mean of Ps. parva to 100 to calculate the Hindex, which was computed using the following formula ( Nikoljukin, 1972 ):
PPT Slide
Lager Image
where H is the numerical value of a character of the hybrid, and M1 and M2 are those of the same character in Pu. herzi and Ps. parva , respectively.
Comparison of six meristic and 19 morphometric characters of a natural hybrid and its supposed parent species,Pungtungia herziandPseudorasbora parva, captured in the Geumho River, a tributary of the Nakdong River basin of Korea. Hybrid index (Hindex) values were calculated for each character state based on the mean values of the two parent species
PPT Slide
Lager Image
Comparison of six meristic and 19 morphometric characters of a natural hybrid and its supposed parent species, Pungtungia herzi and Pseudorasbora parva, captured in the Geumho River, a tributary of the Nakdong River basin of Korea. Hybrid index (Hindex) values were calculated for each character state based on the mean values of the two parent species
- Analyses of DNA sequences
A section of the pectoral or anal fin was excised from each specimen to extract genomic DNA (gDNA). The tissue was incubated with TNES-Urea buffer (10 mM Tris-HCl, pH 8.0; 125 mM NaCl; 10 mM EDTA, pH 8.0; 1% SDS; 8 M urea) ( Asahida et al., 1996 ) containing 100-μg proteinase K (Sigma-Aldrich, St. Louis, MO, USA) at 37℃ for one day, followed by separation of the aqueous phase with a phenol : chloroform : isoamyl alcohol (25 : 24 : 1) solution, and ethanol precipitation according to Sambrook and Russell (2001) . The extracted gDNA was resuspended in 1× TE buffer (10 mM Tris-HCl, pH 8.0; 1 mM EDTA, pH 8.0). Its quantity and quality were checked using a spectrophotometer, NanoDrop 1000 (Thermo Fisher Scientific, Wilmington, DE, USA) and by electrophoresis in a 0.7% agarose gel after staining with GelRed Nucleic Acid Gel Stain (Biotium, Hayward, CA, USA).
To amplify a complete fragment of the mitochondrial cytochrome b gene ( mt-cyb ) and a partial fragment of the nuclear recombination activating protein gene 1 gene ( rag1 ), a PCR was carried out in a 20-μL reaction volume using AccuPower ® PCR Premix (Bioneer, Daejeon, Korea), including 50-ng gDNA and 0.2 μM forward and reverse primers [trnE-14317f (5′-GAYTTGAAGAACCAYCGTTG-3′) and trnP-15601r (5′-ATTHTGGCTTTGGGAGYCAG-3′) for the mt-cyb , newly designed in this study, and RAG1-1495f3 (5′-CAGTAYCAYAAGATGTACCG-3′) and RAG1-3067r (5′-TTGTGAGCYTCCATRAACTT-3′) for the rag1 ( Kim and Bang, 2010 )]. The PCR was conducted using the following thermal cycling profile in a DNA Engine DYAD Peltier Thermal Cycler (MJ Research, Waltham, MA, USA): an initial denaturation at 94℃ for 3 min, 35 cycles of denaturation at 94℃ for 30 s, annealing at 55℃ for 30 s and elongation at 72℃ for 1 min. The reaction was completed with a final elongation at 72℃ for 7 min. The PCR product was purified with the AccuPrep ® PCR Purification Kit (Bioneer).
The purified PCR product of the rag1 of the hybrid fish was cloned into the pGEM ® -T Easy Vector (Promega, Madison, WI, USA) and transformed into competent cells ( Escherichia coli XL1-Blue MRF’; Stratagene, La Jolla, CA, USA) according to Chung et al. (1989) . A total of 12 white E. coli colonies were randomly picked, and their plasmid DNAs were extracted using the alkaline lysis method ( Sambrook and Russell, 2001 ).
After cycle sequencing using the ABI PRISM ® BigDye Terminator v3.1 Cycle Sequencing Ready Reaction Kit (Applied Biosystems, Foster City, CA, USA), the purified PCR products and plasmid DNAs were sequenced by Macrogen (Seoul, Korea), a commercial company, on an ABI 3730xl DNA Analyzer (Applied Biosystems) using both sets of PCR primers. The sequences analyzed in this study were deposited in GenBank under accession numbers KP053615–KP053620.
- Phylogenetic analysis
DNA sequences of the rag1 and mt-cyb of the natural hybrid and its parents ( Pu. herzi and Ps. parva ) analyzed in this study, and those of closely related gobionine species retrieved from GenBank ( ) were aligned with ClustalW in BioEdit 7.0.5 ( Hall, 1999 ). Sarcocheilichthys variegatus microoculus (GenBank accession numbers EU711157 and AB054124, respectively) and Sarcocheilichthys parvus (EU409625 and EF193422, respectively) were used as outgroups to reconstruct the phylogenetic trees. The nucleotide sequence matrices for the rag1 and mt-cyb consisted of 1,488 and 1,140 base pair (bp), respectively.
Phylogenetic analyses were carried out in PAUP* 4.0b10 ( Swofford, 2002 ). A neighbour-joining (NJ) tree was reconstructed with the Kimura 2-parameter model, and a maximum-parsimony (MP) tree with the heuristic search option with the random addition of sequences (10 replicates) and tree-bisection-reconnection (TBR) branch swapping. The robustness of the tree topologies was evaluated by bootstrap analyses for NJ and MP analyses with 1,000 pseudoreplicates.
- Morphology
The morphological data for six meristic and 19 morphometric characters along with their Hindex values are presented in Table 2 . A Hindex of 50 indicates an exact phenotypic intermediacy between the two parent species. The hybrid was considered to resemble Pu. herzi when the Hindex was < 30, and Ps. parva when the Hindex was > 70; it was considered to be intermediate between the two species when the Hindex was between 30 and 70 ( Šorić, 2004 ).
Of the meristic characters, the number of dorsal and ventral fin rays did not differ among the hybrid and parent species, but the number of anal fin rays of the hybrid and Pu. herzi (iii6) were different than those of Ps. parva (ii6). The scales above the lateral line of the hybrid were more similar to Ps. parva (Hindex = 71), whereas the lateral line scales were more similar to Pu. herzi (Hindex = 26), and the number of vertebrae was intermediate (Hindex = 39). Of the morphometric characters, the anal base length of the hybrid was closer to Pu. herzi (Hindex = 9), while the predorsal length, caudal peduncle depth, snout length and interorbital width were closer to Ps. parva (Hindex = 86–99). Four characters, including the prepectoral length, preanal length, ventral-anal length and eye diameter of the hybrid, were intermediate (Hindex = 39–53). The remaining nine morphometric characters deviated from both parent species.
Mouth shape is an important morphological character that distinguishes two parent species ( Mori, 1935 ). Pungtungia herzi has a thick, horseshoe-shaped mouth with a pair of barbels and its upper lip is longer than the lower lip, while Ps. parva has a short, thin upper lip in comparison to the lower lip and lacks barbels ( Fig. 2 ). The natural hybrid possessed thin and identical upper and lower lips lacking barbels.
PPT Slide
Lager Image
The shape of the mouth of a natural hybrid and its supposed parent species, Pungtungia herzi (♀) and Pseudorasbora parva (♂). A) Pu. herzi (♀), B) a natural hybrid and C) Ps. parva (♂).
- Sequence analyses of therag1andmt-cyb
Partial DNA sequences of the single-copy nuclear rag1 of the natural hybrid and its supposed parent species were analyzed to verify the hybrid state. Both parent species were monomorphic; the electropherograms consisted of unique single peaks at all 1,488-bp positions with 41-bp differences ( Fig. 3 ), showing 97.2% overall identity. Meanwhile, the natural hybrid was polymorphic; its electropherogram displayed double peaks at all 41-bp positions, along with one additional peak at position 711. The hybrid showed an identical ratio of double peaks at the polymorphic base positions where the parent species differed ( Fig. 1 ).
PPT Slide
Lager Image
Nucleotide sequence alignment of the nuclear recombination activating gene 1 gene (rag1) of a natural hybrid and its supposed parent species, Pungtungia herzi (♀) and Pseudorasbora parva (♂). The polymorphic base positions are indicated according to IUPAC codes.
Twelve randomly chosen PCR clones of the rag1 were sequenced to confirm the hybrid status at the polymorphic base positions. The results confirmed the presence of two distinct haplotypes. Seven of the PCR clones corresponded to Pu. herzi , and five to Ps. parva . The consensus sequences of the haplotypes showed 99.9–100% and 99.1–100% identities to those of the parent species, Pu. herzi and Ps. parva , respectively. In the phylogenetic tree, each haplotype consistently clustered with one of the parents, with the highest bootstrap supports ( Fig. 4A ). Thus, the hybrid possessed two heterogeneous copies of the rag1 , one from each parent species.
PPT Slide
Lager Image
Phylogenetic trees inferred from the nuclear recombination activating gene 1 gene (rag1) (A) and mitochondrial cytochrome b gene (mt-cyb) (B) of a natural hybrid and its supposed parent species, Pungtungia herzi (♀) and Pseudorasbora parva (♂). Bootstrap values above 50% for the neighbour-joining and maximum-parsimony trees, respectively, are shown at each branch node. Species analyzed in this study are in bold faced font. Two haplotypes of rag1 of the hybrid were deduced from twelve randomly chosen PCR clones.
Analyses of the complete nucleotide sequences of the mitochondrial mt-cyb of the natural hybrid and its supposed parents revealed the direction of hybridization. Pungtungia herzi and Ps. parva displayed 162–180-bp differences out of 1,140-bp positions, showing 83.7–85.7% identities. There were only 5–22-bp differences (97.6–99.5% identities) between the hybrid and P. herzi , whereas there were as many as 158–168-bp differences (85.2–86.1% identities) between the hybrid and Ps. parva . In the phylogenetic tree, the hybrid consistently clustered with Pu. herzi with the highest bootstrap support, showing the closest affiliation with the specimen collected in this study from the Nakdong River basin with 82% and 97% bootstrap supports in the NJ and MP trees, respectively ( Fig. 4B ). However, it showed a clear phylogenetic distinction from Ps. parva .
A natural hybrid of a probable intergeneric mating between Pu. herzi and Ps. parva was captured in the Geumho River, a tributary of the Nakdong River basin in Korea, and its hybrid state and parentage were verified, based on morphological characters and genetic markers. The putative hybrid displayed intermediate traits in some morphometric and meristic characters, although polarity was only observed in one parent or the other, and vigor from the parent species was also evident.
The nuclear rag1 indicated mixed inheritance in the hybrid with respect to its two parent species. Specifically, the hybrid had polymorphic base positions in its DNA sequence, one of which was inherited from Pu. herzi and the other from Ps. parva , in contrast to the base monomorphic positions in the parent species. In addition, the mitochondrial mt-cyb of the hybrid showed the highest sequence identity to Pu. herzi , implying that this was the maternal species. Thus, the genetic markers strongly suggested the maternity of Pu. herzi and paternity of Ps. parva for the natural hybrid. In addition, the hybrid was considered a F 1 hybrid, based on the identical ratio of double peaks at the polymorphic base positions, at which the parent species differed. The single-copy rag1 used in this study ( Kim and Bang, 2010 ) was apparently advantageous over a multi-copy ribosomal RNA gene ( Sonnenberg et al., 2007 ) for determining hybrid status, because the former is inherited in a simple Mendelian fashion.
Pungtungia herzi and Ps. parva are sympatric, but have different habitat preferences and spawning grounds. The former lives on stones or rocks and spawns in host nests (brood parasitism) or crevices between stones and rocks ( Yamane et al., 2009 ), whereas the later inhabits aquatic vegetation and spawns on various substrata such as plants, stones and mollusk shells ( Gozlan et al., 2010 ; Witkowski, 2011 ). In addition, Ps. parva exhibits aggressive paternal courtship ( Konishi and Takata, 2004 ). These ecological differences and behavioural characteristics appear to contribute to a pre-mating reproductive barrier. However, the spawning periods of Pu. herzi and Ps. parva overlap from April to May in Korea ( Kim and Park, 2002 ; Kim et al., 2005 ), which may enable mating. In addition, they show close taxonomic ( Banarescu and Nalbant, 1973 ; Kim, 1984 ; Kim and Kang, 1989 ) and molecular phylogenetic ( Yang et al., 2006 ; Liu et al., 2010 ) relationships among species of the Gobioninae, and have the same number of chromosomes ( Ojima et al., 1972 ). Furthermore, Pu. herzi was reported to engage in intergeneric mating with Pseudopungtungia nigra ( Kim et al., 1991 ), and Ps. parva with congeneric Pseudorasbora pumila pumila ( Koga and Goto, 2005 ). Taken together, these properties of the two parent species may increase the probability of intergeneric hybridization as presented in this study. Unfortunately, no other hybrid specimens were found in the region where the hybrid in this study was captured (S.J. Cho, personal observation). Thus, at present, the occurrence of natural hybrids between Pu. herzi and Ps. parva appears to be a rare phenomenon resulting from the accidental contact of gametes.
We are grateful to Dr. G-S Kim of Korea Institute of Ocean freshScience and Technology for his valuable technical assistance in molecular works and Prof. JY Park of Chonbuk National University for providing a soft X-ray photographer.
Albert V , Jónsson B , Bernatchez L 2006 Natural hybrids in Atlantic eels (Anguilla anguilla, A. rostrata): evidence for successful reproduction and fluctuating abundance in space and time Mol Ecol 15 1903 - 1916    DOI : 10.1111/j.1365-294X.2006.02917.x
Allendorf FW , Leary RF , Spruell P , Wenburg JK 2001 The problems with hybrids: setting conservation guidelines Trends Ecol Evol 16 613 - 622    DOI : 10.1016/S0169-5347(01)02290-X
Asahida T , Kobayashi T , Saitoh K , Nakayama I 1996 Tissue preservation and total DNA extraction from fish stored at ambient temperature using buffers containing high concentration of urea Fish Sci 62 727 - 730    DOI : 10.2331/suisan.62.727
Avise JC , Saunders NC 1984 Hybridization and introgression among species of sunfish (Lepomis): analysis by mitochondrial DNA and allozyme markers Genetics 108 237 - 255
Avise JC , Nelson WS , Arnold J , Koehn RK , Williams GC , Thorsteinsson V 1990 The evolutionary genetic status of Icelandic eels Evolution 44 1254 - 1262    DOI : 10.2307/2409286
Banarescu P , Nalbant TT 1973 Pisces, Teleostei, Cyprinidae (Gobioninae). Das Tierreich Lieferung 93 Walter de Gruyter Berlin
Chelomina GN , Rozhkovan KV , Ivanov SA 2008 Discrimination of interspecific hybrids in natural populations of Amur sturgeon fish by means of multilocus RAPD-PCR markers Cytol Genet 42 342 - 350    DOI : 10.3103/S0095452708050083
Chung CT , Niemela SL , Miller RH 1989 One-step preparation of competent Escherichia coli: transformation and storage of bacterial cells in the same solution Proc Natl Acad Sci USA 86 2172 - 2175    DOI : 10.1073/pnas.86.7.2172
DeMarais BD , Dowling TE , Douglas ME , Minckley WL , Marsh PC 1992 Origin of Gila seminuda (Teleostei: Cyprinidae) through introgressive hybridization: implications for evolution and conservation Proc Natl Acad Sci USA 89 2747 - 2751    DOI : 10.1073/pnas.89.7.2747
Dowling TE , DeMarais BD 1993 Evolutionary significance of introgressive hybridization in cyprinid fishes Nature 362 444 - 446    DOI : 10.1038/362444a0
Dowling TE , Smith GR , Brown WM 1989 Reproductive isolation and introgression between Notropis cornutus and Notropis chrysocephalus (Family Cyprinidae): comparison of morphology, allozymes, and mitochondrial DNA Evolution 43 620 - 634    DOI : 10.2307/2409064
Gerber AS , Tibbets CA , Dowling TE 2001 The role of introgressive hybridization in the evolution of the Gila robusta complex (Teleostei: Cyprinidae) Evolution 55 2028 - 2039    DOI : 10.1111/j.0014-3820.2001.tb01319.x
Gozlan RE , Andreou D , Asaeda T , Beyer K , Bouhadad R , Burnard D , Caiola N , Cakic P 2010 Pan-continental invasion of Pseudorasbora parva: towards a better understanding of freshwater fish invasions Fish Fish 11 315 - 340    DOI : 10.1111/j.1467-2979.2010.00361.x
Grant PR , Grant BR 1992 Hybridization of bird species Science 256 193 - 197    DOI : 10.1126/science.256.5054.193
Hall TA 1999 BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT Nucl Acid Symp Ser 41 95 - 98
Hubbs CL 1955 Hybridization between fish species in nature Syst Zool 4 1 - 20    DOI : 10.2307/2411933
Hubbs CL , Lagler KF 1964 Fishes of the Great Lakes Region The University of Michigan Press Ann Arbor
Kwan Y-S , Ko M-H , Won Y-J 2014 Genomic replacement of native Cobitis lutheri with introduced C. tetralineata through a hybrid system following the artificial connection of river systems Ecol Evol 4 1451 - 1465    DOI : 10.1002/ece3.1027
Kim I-S 1984 The taxonomic study of gudgeons of the subfamily Gobioninae (Cyprinidae) in Korea Bull Korean Fish Soc 17 436 - 448
Kim K-Y , Bang I-C 2010 Molecular phylogenetic position of Abbottina springeri (Cypriniformes; Cyprinidae) based on nucleotide sequences of RAG1 gene Korean J Ichthyol 22 273 - 278
Kim I-S , Kang E-J 1989 Comparative study on the urohyal of the subfamily Gobioninae of Korea Korean J Ichthyol 1 24 - 34
Kim I-S , Park J-Y 2002 Freshwater Fishes of Korea Kyohak Publishing Seoul, Korea
Kim I-S , Choi Y , Shim J-H 1991 An occurrence of intergeneric hybrid cross, Pungtungia herzi × Pseudopungtungia nigra from the Ungcheon River, Korea Korean J Ichthyol 3 42 - 47
Kim I-S , Choi Y , Lee C-L , Lee Y-J , Kim B-J , Kim J-H 2005 Illustrated Book of Korean Fishes Kyohak Publishing Seoul, Korea
Koga K , Goto A 2005 Genetic structures of allopatric and sympatric populations in Pseudorasbora pumila pumila and Pseudorasbora parva Ichthyol Res 52 243 - 250    DOI : 10.1007/s10228-005-0280-8
Konishi M , Takata K 2004 Size-dependent male-male competition for a spawning substrate between Pseudorasbora parva and Pseudorasbora pumila Ichthyol Res 51 184 - 187    DOI : 10.1007/s10228-003-0203-5
Lee I-R , Yang H , Kim J-H , Kim K-Y , Bang I-C 2009 Identification of a natural hybrid between the striped spine loach Cobitis tetralineata and the king spine loach Iksookimia longicorpa by analyzing mitochondrial COI and nuclear RAG1 sequences Korean J Ichthyol 21 287 - 290
Liu H , Yang J , Tang Q 2010 Estimated evolutionary tempo of East Asian gobionid fishes (Teleostei: Cyprinidae) from mitochondrial DNA sequence data Chinese Sci Bull 55 1501 - 1510    DOI : 10.1007/s11434-010-3159-7
Mori T 1935 Descriptions of two new genera and seven new species of Cyprinidae from Chosen Annot Zool Japon 15 161 - 181
Moritz C , Dowling TE , Brown WM 1987 Evolution of animal mitochondrial DNA: relevance for population biology and systematics Annu Rev Ecol Syst 18 269 - 292    DOI : 10.1146/
Nikoljukin MJ 1972 Distant Hybridization in Acipenseridae and Teleostei Moskova
Ojima Y , Hayashi M , Ueno K 1972 Cytogenetic studies in lower vertebrates. X. Karyotype and DNA studies in 15 species of Japanese Cyprinidae Jpn J Genet 47 431 - 440    DOI : 10.1266/jjg.47.431
Rhymer JM , Simberloff D 1996 Extinction by hybridization and introgression Annu Rev Ecol Syst 27 83 - 109    DOI : 10.1146/annurev.ecolsys.27.1.83
Sambrook J , Russell DW 2001 Molecular Cloning: A Laboratory Manual Third Edition Cold Spring Harbor Laboratory Press New York, USA
Scribner KT , Page KS , Bartron ML 2001 Hybridization in freshwater fishes: a review of case studies and cytonuclear methods of biological inference Rev Fish Biol Fish 10 293 - 323
Sonnenberg R , Nolte AW , Tautz D 2007 An evaluation of LSU rDNA D1-D2 sequences for their use in species identification Front Zool 4 6 -    DOI : 10.1186/1742-9994-4-6
Šorić VM 2004 A natural hybrid of Leuciscus cephalus and Alburnus alburnus (Pisces, Cyprinidae) from the Ibar River (Western Serbia) Arch Biol Sci Belgrade 56 23 - 32    DOI : 10.2298/ABS0402023S
Swofford DL 2002 PAUP*: Phylogenetic Analysis Using Parsimony (*and Other Methods). Ver. 4 Sinauer Associates Sunderland
Ünver B , Tatlidil H , Erk’akan F 2008 Biometrical features of intergeneric hybrid between Leuciscus cephalus (L.) and Chalcalburnus chalcoides (G.) (Osteichthyes-Cyprinidae) distributed in Lake Tödürge (Sivas-Turkey) Turk J Fish Aquat Sci 8 207 - 213
Verspoor E , Hammar J 1991 Introgressive hybridization in fishes: the biochemical evidence J Fish Biol 39 (Suppl A) 309 - 334    DOI : 10.1111/j.1095-8649.1991.tb05094.x
Witkowski A 2011 NOBANIS – Invasive Alien Species Fact Sheet – Pseudorasbora parva Retrieved from on January 30, 2014
Yamane H , Watanabe K , Nagata Y 2009 Flexibility of reproductive tactics and their consequences in the brood parasitic fish Pungtungia herzi (Teleostei: Cyprinidae) J Fish Biol 75 563 - 574    DOI : 10.1111/j.1095-8649.2009.02306.x
Yang J , He S , Freyhof J , Witte K , Liu H 2006 The phylogenetic relationships of the Gobioninae (Teleostei: Cyprinidae) inferred from mitochondrial cytochrome b gene sequences Hydrobiology 553 255 - 266    DOI : 10.1007/s10750-005-1301-3
Yoon E-Y , Lee I-R , Park SY , Kang E-J , Kim E-O , Yang S-K , Nam YK , Bang I-C 2009 Genetic identification of hybrids between Rhodeus uyekii and R. notatus by sequence analysis of RAG-1 gene J Aquacult 22 79 - 82