Phylogenetic affiliation of ancient and contemporary humans inferred from mitochondrial DNA.

Nucleotide sequence analysis of the major non-coding region of human mitochondrial DNA (mtDNA) from three major races was extended with data from 27 contemporary Mongoloids (20 from southeast Asia, seven from America) and 11 Ancient Japanese bones (five from Jomon Age; 3000-6000 years BP, six from the early modern Ainu; 200-300 years BP). In both cases, the sequence was determined directly from the polymerase chain reaction products. Based on a comparison of the 482 base pair sequences from a total of 128 contemporary humans, the nucleotide diversity is estimated to be 1.46%, which is three times higher than the corresponding value estimated from restriction-enzyme analysis of the whole mtDNA genome. The phylogenetic tree revealed that all lineages are classified into at least five clusters designated as C1-C5. C1 consists exclusively of Africans, and most Asians and Europeans formed C2, C3, C5 and C4, respectively. Phylogenetic analysis also indicated that part of the Asians, including the Japanese, subsequently diverged from the majority of Africans, and that Asians can therefore be separated into two distinct groups. Native Americans, however, appeared only in C3 and C5, suggesting that the size of the founder population was not so large during the peopling of American. Nucleotide sequences derived from ancient bones in a highly polymorphic region were also compared with those of contemporary humans. The nucleotide diversity among the 139 sequences in the region was estimated to be 2.26%. A group of ancient Japanese, including both Jomon peoples and the Ainu, showed a close phylogenetic affiliation with one group of contemporary Japanese and southeast Asians.(ABSTRACT TRUNCATED AT 250 WORDS)     

Molecular phylogeny of macaques: implications of nucleotide sequences from an 896-base pair region of mitochondrial DNA

We determined the nucleotide sequences of an 896-base pair region of mitochondrial DNA (mtDNA) from 20 primates representing 13 species of macaques, a baboon, and a patas. We compared these sequences and the homologous sequences from four macaques and a human against each other and deduced the phylogenetic relationships of macaques. The results from the phylogenetic analyses revealed five groups among the macaques: (1) Barbary macaque, (2) two species of Sulawesi macaques, (3) Japanese, rhesus, Taiwanese, crab-eating, and stump-tailed macaques, (4) toque, pig-tailed, and lion-tailed macaques, and (5) Assamese and bonnet macaques. The phylogenetic position of Tibetan macaque remains ambiguous as to whether it belongs to the fourth or fifth group. Phylogenetic trees revealed that Barbary macaque diverged first from the other Asian macaques. Subsequently, the four groups of Asian macaques diverged from one another in a relatively short period of time. Within each group, most of the species diverged in a relatively short period of time following the divergence of the groups. Assuming that the Asian macaques diverged from the outgroup Barbary macaque three million years ago (MYA), the divergence times among groups of Asian macaques were estimated at 2.1-2.5 MYA and within groups at 1.4- 2.2 MYA. The intraspecific nucleotide diversity observed among three rhesus macaques was so large that they did not form a monophyletic cluster in the phylogenetic trees. Instead, one of them formed a cluster with Japanese and Taiwanese macaques, whereas the other two formed a separate cluster. This implies that either polymorphisms of mtDNA sequences that existed before the divergence of these three species (ca. 700,000 years ago) have been retained in rhesus macaques or introgression has occurred among the three species.      

Genetic relationship amongst the major non-coding regions of mitochondrial DNAs in wild boars and several breeds of domesticated pigs

We completed phylogenetic analysis of the major non-coding region of the mitochondrial DNA (mtDNA) from 159 animals of eight Euro-American and six East Asian domesticated pig breeds and 164 Japanese and five European wild boars. A total of 62 mtDNA haplotypes were detected. Alignment of these regions revealed nucleotide variations (including gaps) at 73 positions, including 58 sites with transition nucleotide substitutions, and two transversion substitutions. Phylogenetic analysis of the sequences could not organize domestic pig breeds into discrete clusters. In addition, many of the haplotypes found in members of diverged clustering groups were found primarily in Euro-American pig breeds, indicating extensive introgression of Asian domestic pigs into European breeds. Furthermore, phylogenetic analysis allocated the DNA sequences of non-coding regions into two different groups, and the deepest branchpoint of this porcine phylogeny corresponded to 86 000-136 000 years before present. This time of divergence would predate the historical period when the pig is thought to have been domesticated from the wild boar.     

Reduced-median-network analysis of complete mitochondrial DNA coding-region sequences for the major African, Asian, and European haplogroups

The evolution of the human mitochondrial genome is characterized by the emergence of ethnically distinct lineages or haplogroups. Nine European, seven Asian (including Native American), and three African mitochondrial DNA (mtDNA) haplogroups have been identified previously on the basis of the presence or absence of a relatively small number of restriction-enzyme recognition sites or on the basis of nucleotide sequences of the D-loop region. We have used reduced-median-network approaches to analyze 560 complete European, Asian, and African mtDNA coding-region sequences from unrelated individuals to develop a more complete understanding of sequence diversity both within and between haplogroups. A total of 497 haplogroup-associated polymorphisms were identified, 323 (65%) of which were associated with one haplogroup and 174 (35%) of which were associated with two or more haplogroups. Approximately one-half of these polymorphisms are reported for the first time here. Our results confirm and substantially extend the phylogenetic relationships among mitochondrial genomes described elsewhere from the major human ethnic groups. Another important result is that there were numerous instances both of parallel mutations at the same site and of reversion (i.e., homoplasy). It is likely that homoplasy in the coding region will confound evolutionary analysis of small sequence sets. By a linkage-disequilibrium approach, additional evidence for the absence of human mtDNA recombination is presented here.     

Intra- and interbreed genetic variations of mitochondrial DNA major non-coding regions in Japanese native dog breeds [Canis familiaris)

Mitochondrial DNA (mtDNA) major non-coding regions were amplified from 73 dogs of eight Japanese native dog breeds and from 21 dogs of 16 non-Japanese dog breeds by the polymerase chain reaction and their DNA sequences were determined. A total of 51 nucleotide positions within the non-coding region (969–972 base pairs) showed nucleotide variations of which 48 were caused by transition. These nucleotide substitutions were abundant in the region proximate to tRNA^Pro. In addition to the nucleotide substitutions, the dog mtDNA D-loop sequences had a heteroplasmic repetitive sequence (TACACGTÀCG) involving size variation. The DNA sequences of the non-coding region were classified into four different groups by phylogenetic analysis and the deepest branchpoints of this dog phylogeny was calculated to about 100 000 years before the present. Phylogenetic analysis showed that Japanese native dog breeds could not be clearly delimited as distinct breeds. Many haplotypes found in members of some clustering groups were seen in each dog breed, and interbreed nucleotide differences between Japanese dog breeds were almost the same as the intrabreed nucleotide diversities.     

Human Mitochondrial DNA Variation and Evolution: Analysis of Nucleotide Sequences from Seven Individuals

We have analyzed nucleotide sequence variation in an approximately 900-base pair region of the human mitochondrial DNA molecule encompassing the heavy strand origin of replication and the D-loop. Our analysis has focused on nucleotide sequences available from seven humans. Average nucleotide diversity among the sequences is 1.7%, several-fold higher than estimates from restriction endonuclease site variation in mtDNA from these individuals and previously reported for other humans. This disparity is consistent with the rapidly evolving nature of this noncoding region. However, several instances of convergent or parallel gain and loss of restriction sites due to multiple substitutions were observed. In addition, other results suggest that restriction site (as well as pairwise sequence) comparisons may underestimate the total number of substitutions that have occurred since the divergence of two mtDNA sequences from a common ancestral sequence, even at low levels of divergence. This emphasizes the importance of recognizing the large standard errors associated with estimates of sequence variability, particularly when constructing phylogenies among closely related sequences. Analysis of the observed number and direction of substitutions revealed several significant biases, most notably a strand dependence of substitution type and a 32-fold bias favoring transitions over transversions. The results also revealed a significantly nonrandom distribution of nucleotide substitutions and sequence length variation. Significantly more multiple substitutions were observed than expected for these closely related sequences under the assumption of uniform rates of substitution. The bias for transitions has resulted in predominantly convergent or parallel changes among the observed multiple substitutions. There is no convincing evidence that recombination has contributed to the mtDNA sequence diversity we have observed.     

Mitochondrial 16S rRNA sequence diversity of hominoids

We determined nucleotide sequences of the 16S rRNA gene of mitochondrial DNA (mtDNA) (about 1.6 kb) for 35 chimpanzee, 13 bonobo, 10 gorilla, 16 orangutan, and 23 gibbon individuals. We compared those data with published sequences and estimated nucleotide diversity for each species. All the ape species showed higher diversity than human. We also constructed phylogenetic trees and networks. The two orangutan subspecies were clearly separated from each other, and Sumatran orangutans showed much higher nucleotide diversity than Bornean orangutans. Some gibbon species did not form monophyletic clusters, and variation within species was not much different from that among species in the subgenus Hylobates.     

中国狼不同地理种群遗传多样性及系统进化分析

为了解中国狼不同地理种群遗传多样性及系统发育情况,从中国境内狼的主要分布区青海、新疆、内蒙古和吉林4个地区采集样品,用分子生物学技术手段成功地获得44个个体线粒体DNA控制区第一高变区(HVRⅠ)序列和40个线粒体Cyt b部分序列。线粒体控制区HVRⅠ共检测到51个变异位点,位点变异率为8.76%;线粒体Cyt b部分序列发现31个变异位点,位点变异率为5.33%,未见插入及缺失现象,变异类型全部为碱基置换。共定义了16个线粒体HVRⅠ单倍型,其中吉林与内蒙种群存在共享单倍型,估计这两地间种群亲缘关系较近。4个地理种群中新疆种群拥有较高的遗传多样性(0.94)。中国狼种群总体平均核苷酸多态性为2.27%,与世界其他国家地区相比,中国狼种群拥有相对较高的遗传多样性。通过线粒体HVRⅠ单倍型构建的系统进化树可以看出,中国狼在进化上分为2大支,其中位于青藏高原的青海种群独立为一支,推测其可能长期作为独立种群进化。基于青海种群与新疆,内蒙种群的线粒体Cyt b遗传距离,推测中国狼2个世系可能在更新世冰川时期青藏高原受地质作用急速隆起后出现分歧,分歧时间大约在1.1 MY前。     

Low diversity and biased substitution patterns in the mitochondrial DNA control region of sperm whales: implications for estimates of time since common ancestry

The mitochondrial DNA (mtDNA) control region was sequenced in 37 sperm whales from a large part of the global range of the species. Nucleotide diversity was several-fold lower than that reported for control regions of abundant and outbred mammals, but similar to that for populations known to have experienced bottlenecks. Relative neck tests did not suggest that the low diversity is due to a lower substitution rate in sperm whale mtDNA. Rather, it is more likely that demographic factors have reduced diversity. The pattern of nucleotide substitutions was examined by cladistic methods, facilitated by the apparent monophyly of lineages from the Southern Hemisphere, as defined by a single base pair deletion. Substitutions were nonrandom in nature, confined to a few "hot spots," and parallel substitutions constituted a majority of the inferred changes. The substitution pattern fitted a negative binomial distribution better than a Poisson distribution, and the bias in number of substitutions among sites was considerably higher than previously reported for the mtDNA control region of any species. A novel method of estimating time since common ancestry was developed, which utilizes the transition/transversion ratio R and the number of substitutions inferred from a parsimony analysis. Using this method, we estimated the age of sperm whale mtDNA diversity to be about 6,000-25,000 years, and when the uncertainty of R was accounted for, a range of about 1,000- 100,000 years was obtained.      

中国17个黄牛品种mtDNA变异特征与多态性分析

中国黄牛品种资源丰富,尚有28个地方固有品种.为了进一步深入了解这些宝贵遗传资源,本研究通过mtDNA变异特征与多态性分析揭示这些来自中国不同地域地方黄牛的母系起源与分子系统学特征.在17个品种84个体的mtDNA D-loop全序列中,一共检测到了102个核苷酸替代突变位.由此定义的53个单倍型被类聚为2个明显的单倍群：普通牛和瘤牛.mtDNA D-loop全序列变异的第一个特征是转换发生的频率远高于颠换;第二个特征是缺失与替代突变共存;第三个特征是缺失突变率比较高.所有D-loop全序列的核苷酸多样性和单倍型多样性分别为0.026 78±0.000 50和0.919±0.027.普通牛D-loop单倍型在北方牛种群中占有优势(80%～100%),而瘤牛单倍型在南方牛种群中占有优势(42.9%～100%),2种不同单倍型在中原牛种群中的分布也存在差异.2种不同单倍型在中国不同地域17个黄牛品种中的差异性分布揭示出了瘤牛mtDNA基因在中国黄牛中自南而北、由高到低的流动模式,这种基因流动模式的形成可能是由历史事件、地理隔离以及气候环境差异等造成的.     

Intraspecific DNA sequence variation of the mitochondrial control region of white sturgeon (Acipenser transmontanus)

Intraspecific sequence variation in the D-loop region of mtDNA in white sturgeon (Acipenser transmontanus), a relict North American fish species, was examined in 27 individuals from populations of the Columbia and Fraser rivers. Thirty-three varied nucleotide positions were present in a 462-nucleotide D-loop sequence, amplified using the polymerase chain reaction. Bootstrapped neighbor-joining and maximum- parsimony trees of sequences from 19 haplotypes suggest that the two populations have recently diverged. This is consistent with the hypothesis that the Columbia River, a Pleistocene refugium habitat, was the source of founders for the Fraser River after the last glacial recession. On the basis of a divergence time of 10-12 thousand years ago, the estimated substitution rate of the white sturgeon D-loop region is 1.1-1.3 x 10(-7) nucleotides/site/year, which is comparable to rates for hypervariable sequences in the human D-loop region. Furthermore, the ratio of mean percent nucleotide differences in the D- loop (2.27%) to that in whole mtDNA (0.54%, as estimated from restriction-enzyme data) is 4.3, which is similar to the fourfold-to- fivefold-higher substitution rate estimated for the human D-loop. The high nucleotide substitution rate of the hypervariable region indicates that the vertebrate D-loop has potential as a genetic marker in molecular population studies.      

Sequence polymorphism in a novel noncoding region of Pacific oyster mitochondrial DNA

Nucleotide sequence polymorphism in a 641-bp novel major noncoding region of mitochondrial DNA (mtDNA-NC) of the Pacific oyster Crassostrea gigas was analysed for 29 cultured individuals within the Goseong population. A total of 30 variable sites were detected, and the relative frequency of nucleotide alteration was determined to be 4.68. Alterations were mostly single nucleotide substitutions. Transition, transversion, both transition and transversion, and both transversion and nucleotide deletion were observed at 18, 9, 2 and 1 sites, respectively. Among 29 specimens, 22 haplotypes were identified, and pairwise genetic diversity of haplotypes was calculated to be 0.988 from multiple sequence substitutions using the two-parameter model. A phylogenetic tree, obtained for haplotypes by the neighbor-joining method, showed a single cluster of linkages. The cluster comprised 11 haplotypes associating with 14 specimens, while the other 11 haplotypes associating with 15 specimens were scattered. This mtDNA-NC presenting a high nucleotide sequence polymorphism is a potential mtDNA control region. It therefore can serve as a genetic marker for intraspecies phylogenetic analysis of the Pacific oyster and is more useful than the less polymorphic mtDNA coding genes.     

Model of amino acid substitution in proteins encoded by mitochondrial DNA

Mitochondrial DNA (mtDNA) sequences are widely used for inferring the phylogenetic relationships among species. Clearly, the assumed model of nucleotide or amino acid substitution used should be as realistic as possible. Dependence among neighboring nucleotides in a codon complicates modeling of nucleotide substitutions in protein-encoding genes. It seems preferable to model amino acid substitution rather than nucleotide substitution. Therefore, we present a transition probability matrix of the general reversible Markov model of amino acid substitution for mtDNA-encoded proteins. The matrix is estimated by the maximum likelihood (ML) method from the complete sequence data of mtDNA from 20 vertebrate species. This matrix represents the substitution pattern of the mtDNA-encoded proteins and shows some differences from the matrix estimated from the nuclear-encoded proteins. The use of this matrix would be recommended in inferring trees from mtDNA-encoded protein sequences by the ML method. Received: 3 May 1995 / Accepted: 31 October 1995     

Secondary structure of mitochondrial 12S rRNA among fish and its phylogenetic applications.

The complete 12S ribosomal RNA(rRNA) sequences from 23 gobioid species and nine diverse assortments of other fish species were employed to establish a core secondary structure model for fish 12S rRNA. Of the 43 stems recognized, 41 were supported by at least some compensatory evidence among vertebrates. The rates of nucleotide substitution were lower in stems than in loops. This may produce less phylogenetic information in stems when recently diverged taxa are compared. An analysis of compensatory substitution shows that the percentage of covariation is 68%, and the weighting factor for phylogenetic analyses to account for the dependence of mutations should be 0.66. Different stem-loop weighting schemes applied to the analyses of phylogenetic relationships of the Gobioidei indicate that down-weighting paired regions because of nonindependence could not improve the present phylogenetic analysis. A biased nucleotide composition (adenine% [A%] > thymine% [T%], cytosine% [C%] > guanine% [G%]) in the loop regions was also observed in the mammalian counterpart. The excess of A and C in the loop regions may be because of the asymmetric mechanism of mtDNA replication, which leads to the spontaneous deamination of C and A. This process may also be responsible for a transition-transversion bias and the patterns of nucleotide substitutions in both stems and loops.     

Sequence diversity within the reovirus S2 gene: reovirus genes reassort in nature, and their termini are predicted to form a panhandle motif.

To better understand genetic diversity within mammalian reoviruses, we determined S2 nucleotide and deduced sigma 2 amino acid sequences of nine reovirus strains and compared these sequences with those of prototype strains of the three reovirus serotypes. The S2 gene and sigma 2 protein are highly conserved among the four type 1, one type 2, and seven type 3 strains studied. Phylogenetic analyses based on S2 nucleotide sequences of the 12 reovirus strains indicate that diversity within the S2 gene is independent of viral serotype. Additionally, we found marked topological differences between phylogenetic trees generated from S1 and S2 gene nucleotide sequences of the seven type 3 strains. These results demonstrate that reovirus S1 and S2 genes have distinct evolutionary histories, thus providing phylogenetic evidence for lateral transfer of reovirus genes in nature. When variability among the 12 sigma 2-encoding S2 nucleotide sequences was analyzed at synonymous positions, we found that approximately 60 nucleotides at the 5' terminus and 30 nucleotides at the 3' terminus were markedly conserved in comparison with other sigma 2-encoding regions of S2. Predictions of RNA secondary structures indicate that the more conserved S2 sequences participate in the formation of an extended region of duplex RNA interrupted by a pair of stem-loops. Among the 12 deduced sigma 2 amino acid sequences examined, substitutions were observed at only 11% of amino acid positions. This finding suggests that constraints on the structure or function of sigma 2, perhaps in part because of its location in the virion core, have limited sequence diversity within this protein.     

Phylogenetic relationships within Hevea brasiliensis as deduced from a polymorphic mitochondrial DNA region

We have cloned a 4.5-kb mtDNA fragment showing a high RFLP polymorphism between various Hevea genotypes. Subcloning and sequencing of a 1.4-kb segment of this clone allowed us to design PCR amplification primers to isolate homologous mtDNA segments of about 0.9 kb from 23 representative genotypes of Hevea. Complete sequences from 4 genotypes showed between 6.7% and 20.2% of nucleotide diversity, suggesting the presence of a hypervariable, or hotspot, region. A sequence of 345 nucleotides within this region was determined for the 23 genotypes. The phylogenetic relationships inferred from the sequence comparison are in general agreement with the results obtained from mtDNA RFLP analysis, indicating that this polymorphic mtDNA region is a useful molecular marker for phylogenetic analysis within Hevea.     

Phylogenetic relationships amongTetrahymena species determined using the polymerase chain reaction

Summary The species of theTetrahymena pyriformis complex present a conundrum with regard to their highly conservative morphology and widely divergent molecular characteristics. We have investigated the phylogenetic relationships among these species using the nucleotide sequences from the histone H3II/H4II region of the genome. This region includes portions of the two histone coding sequences, as well as the intergenic region. The DNA sequences of these regions were amplified by the polymerase chain reaction (PCR) and the sequence of each was determined. Nucleotide substitutions and insertions/deletions within this set of sequences were compared to determine the phylogenetic relationships among the species of the complex. These data yield phylogenetic trees with identical topologies when different tree-building routines are used, indicating that the data are very robust.Glaucoma chattoni was used as an outgroup to root the trees for this analysis. The genome organization ofG. chattoni and the divergence of its histone H3II/H4II region sequence relative to those of the complex clearly indicate that this species has diverged considerably from the complex. These results show that PCR amplification analysis is feasible over considerable evolutionary distances. However, DNA-DNA hybridization may be more useful than sequence analysis in resolving the relationships among the closely related species in the complex.     

Genetic variation and phylogenetic relationship between three serow species of the genus Capricornis based on the complete mitochondrial DNA control region sequences

The molecular evidence of phylogenetic status regarding the Formosan serow (Capricornis swinhoei) is not robust and little is known about the genetic diversity of the Sumatran serow (Capricornis sumatraensis), which partly is due to the hardness in sample collection. Here we determined the sequences of the complete mitochondrial DNA control region (1,014 bp) of 19 Sumatran-serow individuals. Nine new haplotypes were defined based on 78 variable sites. Combined analysis with other 32 haplotypes downloaded from the public database, including 1 Sumatran-serow, 11 Formosan-serow and 20 Japanese-serow (Capricornis crispus) haplotypes, a relatively high level of nucleotide diversity was first observed in Sumatran serow (π = 0.0249). By comparative analysis with structural consensus sequences from other mammals, we have identified central, left and right domains and depicted the putative functional structure, including extend termination associated sequences and conserve sequence blocks, in mtDNA control region. The alignment of mtDNA control region revealed that both Sumatran and Japanese serow have two tandem repeats (TRs), but three TRs in Formosan serow. Phylogenetic analyses revealed that the Formosan serow is distinct species with the Japanese serow, but a sister group with the Sumatran serow. The divergence time estimated among three serow species revealed that Pleistocene climate changes and the uplift of Qinghai-Tibetan plateau might play an important role in the genetic differentiation of the serows. These results mainly provide the convinced evidence on the genetic relationship between three serow species.     

Phylogenetic relationships among Japanese, rhesus, Formosan, and crab- eating monkeys, inferred from restriction-enzyme analysis of mitochondrial DNAs [published erratum appears in Mol Biol Evol 1988 Sep;5(5):601]

Mitochondrial DNA (mtDNA) polymorphisms in four species of macaques, i.e., Japanese monkey (Macaca fuscata), rhesus monkey (M. mulatta), Formosan monkey (M. cyclopis), and crab-eating monkey (M. fascicularis), were analyzed to study phylogenetic relationships. When 17 restriction enzymes of 6-bp recognition were used, 42-49 sites were observed in the samples. The estimated number of nucleotide substitutions per site among Japanese, rhesus, and Formosan monkeys ranges from 0.0318 to 0.0396, and that between the crab-eating monkey and the other monkeys from 0.0577 to 0.0653. These findings suggest that the crab-eating monkey diverged from the other three approximately 1.5-3.0 Myr before the present (Mybp) and that the Japanese, rhesus, and Formosan monkeys diverged approximately 0.9-1.8 Mybp, although the branching order cannot be determined conclusively.      

Evidence of a secondary interspecific mitochondrial DNA introgression in the pond loach Misgurnus dabryanus (Teleostei: Cobitidae) population introduced in Japan

The pond loach Misgurnus dabryanus is a freshwater fish with a distribution range spanning the eastern part of the Asian continent, the Korean Peninsula, and Taiwan. The pond loach was transplanted to the Japanese archipelago through the co-inclusion with dojo loach (Misgurnus anguillicaudatus species complex) populations, which were imported live from China for food materials, and it is currently distributed widely across Japan. A previous mitochondrial DNA (mtDNA) analysis revealed that a pond loach population in Ehime Prefecture (Shikoku Island, Japan) included two highly diverged mtDNA groups (Groups I and II). To examine the origin of these two distinct forms of mtDNA within the Japanese pond loach population, we performed phylogenetic analyses using sequences based on the mtDNA of cytochrome oxidase b (cyt b) and the nuclear DNA recombination activating gene 1 (RAG-1). We also conducted a random amplified polymorphic DNA (RAPD) analysis to examine the establishment of reproductive isolation between sympatric pond loaches with two different mtDNA groups. Our mtDNA phylogenetic results indicated that the two diverged pond loach mtDNA sequences showed polyphyletic relationships among Misgurnus species and its related genus Cobitis. In contrast, there were no clear divergence in nuclear DNA among the pond loaches irrespective of their mtDNA groups, and they all formed monomorphic clades in the phylogenetic relationships among the species. The discrepancy between the mtDNA and nuclear DNA genes support that the existence of two diverged forms of DNA within the pond loach population could be attributed to past mtDNA introgressions from other species rather than convergent evolution. Previous mtDNA phylogenetic studies among Cobitidae revealed that the dojo loach also consisted of two genetically diverged polyphyletic clades: an original Misgurnus mtDNA and an introgressed mtDNA from Cobitis species. In our mtDNA result, the Group II haplotype of the pond loach was included in the mtDNA from the introgressed dojo loach. This suggested that the Group II haplotype was derived from introgressed dojo loach mtDNA. The close relationships between the introgressed dojo loach and the pond loach mtDNA indicated that this secondary introgression had recently occurred via hybridization in a recent artificial aquaculture or transportation process. Common RAG-1 alleles and RAPD bands were shared between the sympatric pond loaches with original and introgressed mtDNAs. This indicates that the introgressed mtDNA haplotype is included as one of the polymorphic genotypes within the pond loach populations, and does not represent existence of different cryptic species.