Genetic diversity of the Andean tetraploid cultivated potato (Solanum tuberosum L. subsp. andigena Hawkes) evaluated by chloroplast and nuclear DNA markers.

Andigena potatoes (Solanum tuberosum L. subsp. andigena Hawkes) (2n = 4x = 48) are native farmer-selected important cultivars that form a primary gene pool of the common potato (Solanum tuberosum L. subsp. tuberosum). The genetic diversity of 185 Andigena accessions and 6 Chilean native potatoes (S. tuberosum subsp. tuberosum) was studied using chloroplast DNA (ctDNA) microsatellites and nuclear DNA (nDNA) restriction fragment length polymorphism (RFLP) markers. Andigena potatoes had 14 ctDNA haplotypes and showed higher variability in the central Andes, particularly in Bolivia, whereas those in the northern regions of the distribution area were remarkably uniform with A1 ctDNA and Chilean subsp. tuberosum with T ctDNA. Most of 123 clearly scored RFLP bands using 30 single-copy probes were randomly distributed throughout the distribution area and proved the same gene pool shared among these widely collected accessions. Nevertheless, the geographic trend of the nDNA differentiation from north to south along the Andes and the correlated differentiation between nDNA and ctDNA (r = 0.120) could also be revealed by canonical variates analysis. These results suggest that the genetic diversity in Andigena was brought about primarily from cultivated diploid species but considerably modified through sexual polyploidization and intervarietal and (or) introgressive hybridization and long-distance dispersal of seed tubers by humans.     

Global Population Genetic Structure and Male-Mediated Gene Flow in the Green Turtle (Chelonia Mydas): RFLP Analyses of Anonymous Nuclear Loci

We introduce an approach for the analysis of Mendelian polymorphisms in nuclear DNA (nDNA), using restriction fragment patterns from anonymous single-copy regions amplified by the polymerase chain reaction, and apply this method to the elucidation of population structure and gene flow in the endangered green turtle, Chelonia mydas. Seven anonymous clones isolated from a total cell DNA library were sequenced to generate primers for the amplification of nDNA fragments. Nine individuals were screened for restriction site polymorphisms at these seven loci, using 40 endonucleases. Two loci were monomorphic, while the remainder exhibited a total of nine polymorphic restriction sites and three size variants (reflecting 600-base pair (bp) and 20-bp deletions and a 20-bp insertion). A total of 256 turtle specimens from 15 nesting populations worldwide were then scored for these polymorphisms. Genotypic proportions within populations were in accord with Hardy-Weinberg expectations. Strong linkage disequilibrium observed among polymorphic sites within loci enabled multisite haplotype assignments. Estimates of the standardized variance in haplotype frequency among global collections (FST = 0.17), within the Atlantic-Mediterranean (FST = 0.13), and within the Indian-Pacific (FST = 0.13), revealed a moderate degree of population substructure. Although a previous study concluded that nesting populations appear to be highly structured with respect to female (mitochondrial DNA) lineages, estimates of Nm based on nDNA data from this study indicate moderate rates of male-mediated gene flow. A positive relationship between genetic similarity and geographic proximity suggests historical connections and/or contemporary gene flow between particular rookery populations, likely via matings on overlapping feeding grounds, migration corridors or nonnatal rookeries.     

Molecular research technologies in mitochondrial diseases: the microarray approach

Mitochondria are ubiquitous in eukaryotic cells where they generate much of the cellular energy by the process of oxidative phosphorylation (OXPHOS). The approximately 1500 genes of the mitochondrial genome are distributed between the cytoplasmic, maternally-inherited, mitochondrial DNA (mtDNA) which encodes 37 genes and the nuclear DNA (nDNA) which encompasses the remaining mitochondrial genes. The interplay between the mtDNA and nDNA encoded mitochondrial genes and their role in mitochondrial disorders is still largely unclear. One approach for elucidating the pathophysiology of mitochondrial diseases has been to look at changes in the expression of mtDNA and nDNA-encoded genes in response to specific mitochondrial genetic defects. Initial studies of gene expression changes in response to mtDNA defect employed blot technologies to analyze changes in the expression of individual genes one at a time. While Southern/Northern blot experiments confirmed the importance of nDNA-mtDNA interactions in the pathophysiology of mitochondrial myopathy, the methodology used limited the number of genes that could be analyzed from each patient. This barrier has been overcome, in part by the advent of DNA microarray technology. In DNA microarrays gene sequences or oligonucleotides homologous to gene sequences are arrayed on a solid support. The RNA from the subject is then isolated, the mRNA converted to cDNA and the cDNA labeled with a fluorescent probe. The labeled cDNA is hybridized on the microarray and the fluorescence bound to each array is then quantified. Recently, these technologies have been applied to mitochondrial disease patient tissues and the presence of coordinate changes in mitochondrial gene expression confirmed.     

Oxidative damage to mitochondrial DNA is inversely related to maximum life span in the heart and brain of mammals.

DNA damage is considered of paramount importance in aging. Among causes of this damage, free radical attack, particularly from mitochondrial origin, is receiving special attention. If oxidative damage to DNA is involved in aging, long-lived animals (which age slowly) should show lower levels of markers of this kind of damage than short-lived ones. However, this possibility has not heretofore been investigated. In this study, steady-state levels of 8-oxo-7, 8-dihydro-2'-deoxyguanosine (8-oxodG) referred to deoxyguanosine (dG) were measured by high performance liquid chromatography (HPLC) in the mitochondrial (mtDNA) and nuclear (nDNA) DNA from the heart of eight and the brain of six mammalian species ranging in maximum life span (MLSP) from 3.5 to 46 years. Exactly the same digestion of DNA to deoxynucleosides and HPLC protocols was used for mtDNA and nDNA. Significantly higher (three- to ninefold) 8-oxodG/dG values were found in mtDNA than in nDNA in all the species studied in both tissues. 8-oxodG/dG in nDNA did not correlate with MLSP across species either in the heart (r=-0.68; P<0.06) or brain (r = 0.53; P<0.27). However, 8-oxodG/dG in mtDNA was inversely correlated with MLSP both in heart (r=-0.92; P<0.001) and brain (r=-0.88; P<0.016) tissues following the power function y = a(.)x(b), where y is 8-oxodG/dG and x is the MLSP. This agrees with the consistent observation that mitochondrial free radical generation is also lower in long-lived than in short-lived species. The results obtained agree with the notion that oxygen radicals of mitochondrial origin oxidatively damage mtDNA in a way related to the aging rate of each species.-Barja, G., Herrero, A. Oxidative damage to mitochondrial DNA is inversely related to maximum life span in the heart and brain of mammals.     

Nuclear and chloroplast DNA differentiation in Andean potatoes.

Over 3500 accessions of Andean landraces have been known in potato, classified into 7 cultivated species ranging from 2x to 5x (Hawkes 1990). Chloroplast DNA (ctDNA), distinguished into T, W, C, S, and A types, showed extensive overlaps in their frequencies among cultivated species and between cultivated and putative ancestral wild species. In this study, 76 accessions of cultivated and 19 accessions of wild species were evaluated for ctDNA types and examined by ctDNA high-resolution markers (ctDNA microsatellites and H3 marker) and nuclear DNA restriction fragment length polymorphisms (RFLPs). ctDNA high-resolution markers identified 25 different ctDNA haplotypes. The S- and A-type ctDNAs were discriminated as unique haplotypes from 12 haplotypes having C-type ctDNA and T-type ctDNA from 10 haplotypes having W-type ctDNA. Differences among ctDNA types were strongly correlated with those of ctDNA high-resolution markers (r = 0.822). Differentiation between W-type ctDNA and C-, S-, and A-type ctDNAs was supported by nDNA RFLPs in most species except for those of recent or immediate hybrid origin. However, differentiation among C-, S-, and A-type ctDNAs was not clearly supported by nDNA RFLPs, suggesting that frequent genetic exchange occurred among them and (or) they shared the same gene pool owing to common ancestry.     

A Clinical,Neuropathological and Genetic Study of Homozygous A467T POLG-Related Mitochondrial Disease

Mutations in the nuclear gene POLG (encoding the catalytic subunit of DNA polymerase gamma) are an important cause of mitochondrial disease. The most common POLG mutation, A467T, appears to exhibit considerable phenotypic heterogeneity. The mechanism by which this single genetic defect results in such clinical diversity remains unclear. In this study we evaluate the clinical, neuropathological and mitochondrial genetic features of four unrelated patients with homozygous A467T mutations. One patient presented with the severe and lethal Alpers-Huttenlocher syndrome, which was confirmed on neuropathology, and was found to have a depletion of mitochondrial DNA (mtDNA). Of the remaining three patients, one presented with mitochondrial encephalomyopathy, lactic acidosis and stroke-like episodes (MELAS), one with a phenotype in the Myoclonic Epilepsy, Myopathy and Sensory Ataxia (MEMSA) spectrum and one with Sensory Ataxic Neuropathy, Dysarthria and Ophthalmoplegia (SANDO). All three had secondary accumulation of multiple mtDNA deletions. Complete sequence analysis of muscle mtDNA using the MitoChip resequencing chip in all four cases demonstrated significant variation in mtDNA, including a pathogenic MT-ND5 mutation in one patient. These data highlight the variable and overlapping clinical and neuropathological phenotypes and downstream molecular defects caused by the A467T mutation, which may result from factors such as the mtDNA genetic background, nuclear genetic modifiers and environmental stressors.     

Geographic Structure of Mitochondrial and Nuclear Gene Polymorphisms in Australian Green Turtle Populations and Male-Biased Gene Flow

The genetic structure of green turtle (Chelonia mydas) rookeries located around the Australian coast was assessed by (1) comparing the structure found within and among geographic regions, (2) comparing microsatellite loci vs. restriction fragment length polymorphism analyses of anonymous single copy nuclear DNA (ascnDNA) loci, and (3) comparing the structure found at nuclear DNA markers to that of previously analyzed mitochondrial (mtDNA) control region sequences. Significant genetic structure was observed over all regions at both sets of nuclear markers, though the microsatellite data provided greater resolution in identifying significant genetic differences in pairwise tests between regions. Inferences about population structure and migration rates from the microsatellite data varied depending on whether statistics were based on the stepwise mutation or infinite allele model, with the latter being more congruent with geography. Estimated rates of gene flow were generally higher than expected for nuclear DNA (nDNA) in comparison to mtDNA, and this difference was most pronounced in comparisons between the northern and southern Great Barrier Reef (GBR). The genetic data combined with results from physical tagging studies indicate that the lack of nuclear gene divergence through the GBR is likely due to the migration of sGBR turtles through the courtship area of the nGBR population, rather than male-biased dispersal. This example highlights the value of combining comparative studies of molecular variation with ecological data to infer population processes.     

Origin and phylogeography of African savannah elephants (<Emphasis Type="Italic">Loxodonta africana</Emphasis>) in Kruger and nearby parks in southern Africa

African savannah elephants (Loxodonta africana) occur in fragmented and isolated populations across southern Africa. Transfrontier conservation efforts aim at preventing the negative effects of population fragmentation by maintaining and restoring linkages between protected areas. We sought to identify genetic linkages by comparing the elephants in Kruger National Park (South Africa) to populations in nearby countries (Botswana, Mozambique, Zambia and Zimbabwe). We used a 446 base pair mitochondrial DNA (mtDNA) control region fragment (141 individuals) and 9 nuclear DNA (nDNA) microsatellite markers (69 individuals) to investigate phylogenetic relationships and gene flow among elephant populations. The mtDNA and nDNA phylogeographic patterns were incongruent, with mtDNA patterns likely reflecting the effects of ancient female migrations, with patterns persisting due to female philopatry, and nDNA patterns likely reflecting male-mediated dispersal. Kruger elephant heterozygosity and differentiation were examined, and were not consistent with genetic isolation, a depleted gene pool or a strong founder effect. Mitochondrial DNA geographic patterns suggested that the Kruger population was founded by elephants from areas both north and south of Kruger, or has been augmented through migration from more than one geographic source. We discuss our findings in light of the need for conservation initiatives that aim at maintaining or restoring connectivity among populations. Such initiatives may provide a sustainable, self-regulating management approach for elephants in southern Africa while maintaining genetic diversity within and gene flow between Kruger and nearby regions.     

Regulation of energy metabolism in human cells in aging and diabetes: FoF(1), mtDNA, UCP, and ROS

Recent advances in bioenergetics consist of discoveries related to rotational coupling in ATP synthase (FoF(1)), uncoupling proteins (UCP), reactive oxygen species (ROS) and mitochondrial DNA (mtDNA). As shown in cloned sheep, mammalian genomes are composed of both nuclear DNA (nDNA) and maternal mtDNA. Oxidative phosphorylation (oxphos) varies greatly depending on cellular activities, and is regulated by both gene expression and the electrochemical potential difference of H(+) (Delta muH(+)). The expression of both mtDNA (by mtTFA) and nDNA for oxphos and UCP (by NRFs, etc.) is coordinated by a factor called PGC-1. The Delta muH(+) rotates an axis in FoF(1) that is regulated by inhibitors and ATP-sensitive K(+)-channels. We cultured human rho(o) cells (cells without mtDNA) in synthetic media and elucidated relationships among mtDNA, nDNA, Delta muH(+), UCPs, ROS, and apoptosis. These cells lack oxphos-dependent ROS formation and survive under conditions of high O(2). Cells cultured in the absence of ROS scavengers have proliferated for 40 years. UCPs lower Delta muH(+) and prevent ROS formation and resulting apoptosis. These results were applied to diabetology and gerontology. The pancreatic rho(o) cells did not secrete insulin, and mtDNA mutations caused diabetes, owing to the deficient Delta muH(+). Insulin resistance was closely related to UCPs and other energy regulators. The resulting high-glucose environment caused glycation of proteins and ROS-mediated apoptosis in vascular cells involved in diabetic complications. Telomeres, oxphos, and ROS are determinants in cellular aging. Cell division and ROS shortened telomeres and accelerated aging. In aged cells, Delta muH(+) was reduced by the slow respiration, and this change induced apoptosis. Cybrids made from aged cytoplasts and rho(o) cells showed that both decreased expression of nDNA, and somatic mutations of mtDNA are involved in the slowing of respiration in aged cells.     

Mitochondrial DNA genetic diversity among four ethnic groups in Sierra Leone

Although there are numerous ethnic groups in Sierra Leone, the Mende and Temne together account for approximately 60% of the total population. To see if genetic differences could be observed among ethnic groups in Sierra Leone, the nucleotide sequence of the hypervariable 1 (HV1) region of mitochondrial DNA (mtDNA) was determined from samples of the two major ethnic groups, the Mende (n=59) and Temne (n=121), and of two minor ethnic groups, the Loko (n=29) and Limba (n=67). Among these 276 HV1 sequences, 164 individual haplotypes were observed. An analysis of molecular variance indicated that the distribution of these haplotypes within the Limba sample was significantly different from that of the other ethnic groups. No significant genetic variation was seen between the Mende, Temne, and Loko. These results indicate that distinguishing genetic differences can be observed among ethnic groups residing in historically close proximity to one another. Furthermore, we observed some mitochondrial DNA haplotypes that are common among the Sierra Leone ethnic groups but that have not been observed in other published studies of West African ethnic groups. Therefore, we may have evidence for mtDNA lineages that are unique to this region of West Africa.     

Genetic heterogeneity in familial dysbetalipoproteinemia. The E2(lys146----gln) variant results in a dominant mode of inheritance

As determined by isoelectric focusing, most patients with familial dysbetalipoproteinemia (FD) exhibit the homozygous apolipoprotein (apo) E2E2 phenotype. Only rarely does FD develop in the more common heterozygous phenotypes E3E2 or E4E2. In fact, only 1 to 4% of the E2E2 homozygotes will develop FD. We wondered whether this reduced penetrance of FD in E2E2 homozygotes could be due to additional heterogeneity in the APOE*2 allele. In the literature a number of different mutations causing an E2 isoelectric focusing variant have been described. To study the genetic heterogeneity of the APOE gene, hybridization of enzymatically amplified genomic DNA with mutation-specific oligonucleotide probes was applied. All FD patients (n = 40) with the E2E2 phenotype appeared to be homozygous for the common E2(arg158----cys) mutation. However, all three unrelated patients with the E3E2 phenotype exhibited the rare E2(lys146----gln) mutation due to an A----C substitution at nucleotide position 3,847 of the APOE gene. This mutation was not found among normolipidemic individuals with the E2E2 (n = 13) or E3E2 phenotype (n = 120) selected from a random population sample. Family studies of the three probands heterozygous for the E*2(lys146----gln) allele showed that this rare allele predisposes to FD with high penetrance. We conclude that FD is a genetically heterogeneous disease entity, displaying a recessive mode of inheritance with strongly reduced penetrance in case of the common E2(arg158----cys) variant and with a dominant mode of inheritance with high penetrance in case of the rare E2(lys146----gln) mutant. It should be noted that in this dominant form presymptomatic diagnosis is possible.     

Neoplastic transformation is associated with coordinate induction of nuclear and cytoplasmic oxidative phosphorylation genes

Neoplastic transformation was found to have a marked effect on the expression of nuclear DNA (nDNA)- and mitochondrial DNA (mtDNA)-encoded oxidative phosphorylation (OXPHOS) genes. Examining three pairs of human diploid fibroblasts and their SV 40-transformed counterparts revealed that mRNAs for the nuclear-encoded ATP synthase beta and the adenine nucleotide translocator (ANT) isoform 1 and 2 genes were markedly induced, whereas the mRNA for the ANT isoform 3 gene remained unchanged. The mRNA levels for the mtDNA-encoded 12 S rRNA, ND2, ATPase6+8, COIII, ND5+6, and Cytb genes were also increased, whereas the mtDNA number declined. Similar analysis of a cervical carcinoma (HeLa), fibrosarcoma (HT1080), and an Epstein-Barr virus (EBV)-transformed lymphoblastoid line (EBV-L) revealed that all three ANT isoforms were also expressed in these cells. Hence, changes in the expression of OXPHOS genes may be a common feature of transformed cells.     

Patients with idiopathic cardiomyopathy belong to the same mitochondrial DNA gene family of Parkinson's disease and mitochondrial encephalomyopathy

Comparison of total mitochondrial DNA sequences of patients with idiopathic (deleted or hypertrophic) cardiomyopathy with those of patients with Parkinson's disease and mitochondrial encephalomyopathies revealed distinct clustering of point mutations among patients. Furthermore, an inverse relation was found between the total number of base-substitution and life span of the patients. Among point mutations found in each patient, sequentially diverged six clusters consisting of 14, 10, 7, 1, 2, and 3 mutations, respectively, were detected. Five sub-clusters consisting of 2, 2, 11, 1, and 1 mutations, respectively, were detected. From each cluster, the patient's unique mutations were diverged with three types of the mutations specific for the disease. The divergence allowed construction of a phylogenetic tree which clearly indicated that patients with idiopathic cardiomyopathy belong to the same mitochondrial DNA gene family of Parkinson's disease and mitochondrial encephalomyopathies.     

Bioenergetics in human evolution and disease: implications for the origins of biological complexity and the missing genetic variation of common diseases

Two major inconsistencies exist in the current neo-Darwinian evolutionary theory that random chromosomal mutations acted on by natural selection generate new species. First, natural selection does not require the evolution of ever increasing complexity, yet this is the hallmark of biology. Second, human chromosomal DNA sequence variation is predominantly either neutral or deleterious and is insufficient to provide the variation required for speciation or for predilection to common diseases. Complexity is explained by the continuous flow of energy through the biosphere that drives the accumulation of nucleic acids and information. Information then encodes complex forms. In animals, energy flow is primarily mediated by mitochondria whose maternally inherited mitochondrial DNA (mtDNA) codes for key genes for energy metabolism. In mammals, the mtDNA has a very high mutation rate, but the deleterious mutations are removed by an ovarian selection system. Hence, new mutations that subtly alter energy metabolism are continuously introduced into the species, permitting adaptation to regional differences in energy environments. Therefore, the most phenotypically significant gene variants arise in the mtDNA, are regional, and permit animals to occupy peripheral energy environments where rarer nuclear DNA (nDNA) variants can accumulate, leading to speciation. The neutralist–selectionist debate is then a consequence of mammals having two different evolutionary strategies: a fast mtDNA strategy for intra-specific radiation and a slow nDNA strategy for speciation. Furthermore, the missing genetic variation for common human diseases is primarily mtDNA variation plus regional nDNA variants, both of which have been missed by large, inter-population association studies.     

Genetic counseling in Usher syndrome: linkage and mutational analysis of 10 Colombian families

Usher Syndrome (US), an autosomal recessive disease, is characterized by retinitis pigmentosa (RP), vestibular dysfunction, and congenital sensorineural deafness. There are three recognized clinical types of the disorder. In order to improve genetic counseling for affected families, we conducted linkage analysis and DNA sequencing in 10 Colombian families with confirmed diagnosis of US (4 type I and 6 type II). Seventy-five percent of the US1 families showed linkage to locus USH1B, while the remaining 25% showed linkage to loci USH1B and USH1C. Among families showing linkage to USH1B we found two different mutations in the MYO7A gene: IVS42-26insTTGAG in exon 43 (heterozygous state) and R634X (CGA-TGA) in exon 16 (homozygous state). All six US2 families showed linkage to locus USH2A. Of them, 4 had c.2299delG mutation (1 homozygote state and 3 heterozygous); in the remaining 2 we did not identify any pathologic DNA variant. USH2A individuals with a 2299delG mutation presented a typical and homogeneous retinal phenotype with bilateral severe hearing loss, except for one individual with a heterozygous 2299delG mutation, whose hearing loss was asymmetric, but more profound than in the other cases. The study of these families adds to the genotype-phenotype characterization of the different types and subtypes of US and facilitates genetic counseling in these families. We would like to emphasize the need to perform DNA studies as a prerequisite for genetic counseling in affected families.     

Localization of a mutant gene for cleft palate and ankyloglossia in an X-linked Icelandic family.

Common congenital malformations such as cleft lip and cleft palate are in most cases multifactorial in origin, involving both environmental and genetic components. Molecular biology techniques have enabled the successful chromosomal localization of many mutant genes from disorders that exhibit simple Mendelian segregation, whether autosomally dominant (e.g., Huntington's disease), autosomal recessive (e.g., cystic fibrosis), or X-linked (e.g., Duchenne muscular dystrophy). Studying the genetic aspect of multifactorial disorders is more complex. It requires a model family or families within which the common multifactorial phenotype is displayed as a single gene defect. Such a model has been recently exploited in the form of a large Icelandic family (over 280 members) exhibiting X-linked secondary cleft palate (CP) and ankyloglossia (A) (tongue-tied) as a single gene mutation. Using this family and the large bank of well-characterized DNA probes available for the human X chromosome, the gene for CP + A was localized by linkage analysis to Xq13-q21.1 (LOD score = 3.07, linked to anonymous probe DXYS1). Further fine mapping, using other X probes from this region (confirmed by analysis of DNA from a deletion cell-line) has placed the gene between markers DXYS12 and DXS17 (LOD score = 4.1) at Xq21.3-q22. The approximate distance between these two probes is 5 centimorgans (cM), equivalent to approximately 5 million base pairs. Now that the limits of genetic linkage have been fully tested and there are two markers flanking the defect locus, strategies are being pursued to clone the gene responsible.(ABSTRACT TRUNCATED AT 250 WORDS)     

Multiplex polymerase chain reaction for simultaneous quantitation of human nuclear, mitochondrial, and male Y-chromosome DNA: application in human identification

Human forensic casework requires sensitive quantitation of human nuclear (nDNA), mitochondrial (mtDNA), and male Y-chromosome DNA from complex biomaterials. Although many such systems are commercially available, no system is capable of simultaneously quantifying all three targets in a single reaction. Most available methods either are not multiplex compatible or lack human specificity. Here, we report the development of a comprehensive set of human-specific, target-specific multiplex polymerase chain reaction (PCR) assays for DNA quantitation. Using TaqMan-MGB probes, our duplex qPCR for nDNA/mtDNA had a linear quantitation range of 100 ng to 1 pg, and our triplex qPCR assay for nDNA/mtDNA/male Y DNA had a linear range of 100-0.1 ng. Human specificity was demonstrated by the accurate detection of 0.05 and 5% human DNA from a complex source of starting templates. Target specificity was confirmed by the lack of cross-amplification among targets. A high-throughput alternative for human gender determination was also developed by multiplexing the male Y primer/probe set with an X-chromosome-based system. Background cross-amplification with DNA templates derived from 14 other species was negligible aside from the male Y assay which produced spurious amplifications from other nonhuman primate templates. Mainstream application of these assays will undoubtedly benefit forensic genomics.     

Disturbed population genetics: suspected introgressive hybridization between two Mnais damselfly species (Odonata)

Mnais costalis and M. pruinosa are damselflies (Odonata: Calopterygidae) with low dispersal abilities, both during their aquatic stream-living immature stage and their flying adult stage. A previous nuclear DNA (nDNA) sequencing and morphology study showed that these two species are very closely related, and cohabit widely in western Japan. The two species, however, segregate microhabitats along a stream: M. costalis lives in the lower reaches, and M. pruinosa in the upper reaches. In this study, our analyses were based on mitochondrial DNA (mtDNA), which usually mutates faster and is more variable among individuals than nDNA, and which is inherited maternally. We found that most COI haplotypes were shared between the two species, and that for most study sites interspecific riverine genetic structures were not clarified by mtDNA analysis. Incongruent population genetic structures based on nDNA and mtDNA suggested hybridization and introgression of mtDNA between the two species.