Molecular diversity and phylogeography of Andricus curtisii (Hymenoptera,Cynipidae) in Turkey

A 433-bp fragment of the mitochondrial cytochrome b gene and the entire nuclear ITS2 region were sequenced in an oak gall wasp species, Andricus curtisii, which represents 21 populations across its range in Turkey. Forty cytochrome b haplotypes and 63 ITS2 alleles were discovered among the 161 individuals studied. Remarkable genetic variation was detected in the species, with an average haplotype diversity (h) of 0.72 and 0.76 and an average nucleotide diversity (π) of 0.06 and 0.01 for the cyt b gene and the ITS region, respectively. Most populations were genetically distinct, possessing unique cytochrome b haplotypes. Maximum likelihood and Bayesian phylogenetic analyses and application of the haplotype network to both datasets revealed that the most basally located haplotypes/alleles were from the eastern part of Turkey. The detection of older sequences in the phylogenetic trees of the easterly located populations implied an eastern origin of the species. The application of a conventional insect mitochondrial DNA clock to the dataset suggested a splitting of ingroup haplotypes from the outgroup lineages predating the Pleistocene epoch and the formation of two major haplogroups that coincided with the oscillation of the early Pleistocene glaciations.     

Genetic diversity and population connectivity of the Asian green mussel Perna viridis in South China Sea,inferred from mitochondria DNA markers

The Asian green mussel Perna viridis is ecologically and economically important in the coastal regions of China. In order to characterize the genetic diversity and population connectivity of P. viridis in South China Sea, a 664 bp region of mitochondrial COI gene and a 293 bp region of 16S rRNA gene were sequenced and analyzed for 78 and 92 individuals from four populations in South China Sea, respectively. A total of 15 haplotypes were defined by 14 variable nucleotide sites in COI gene, and 7 haplotypes by 6 variable nucleotide sites in 16S rRNA gene. High haplotype diversity and low nucleotide diversity were observed in COI gene, while moderate haplotype diversity and low nucleotide diversity were observed in 16S rRNA gene. Pairwise F_ST values of COI gene were all negative and those of 16S rRNA gene ranged from −0.01409 to 0.10289. The results showed that no significant genetic divergence (or shallow genetic structure) and high levels of population connectivity among the four populations of P. viridis in South China Sea.     

Molecular insights into the genetic and haplotype diversity among four populations of Catla catla from Madhya Pradesh revealed through mtDNA cyto b gene sequences

In the present investigation, the genetic structure of four populations of Catla catla, sequences of mitochondrial gene, cytochrome b (cyto b) from four populations were sequenced and analyzed. The sequences of mitochondrial regions revealed high haplotype diversity and low nucleotide diversity. The lowest 249 polymorphic sites and 0.00 parsimony informative sites were detected in populations of Fish Federation Pond (CCFFB) whereas highest 330 polymorphic sites and 56 parsimony informative sites were detected in populations of Narmada River (CCNRH) in the cyto b gene sequences in Catla catla populations. The twelve different haplotypes were detected among the four populations studied, lowest population specific haplotype as 2.00 was observed in Fish Federation Pond (CCFFB) and highest was in Population of Narmada River and Tighra reservoir. Sequencing of cyto b gene revealed 12 number of haplotypes (h) with haplotype (gene) diversity (Hd) 0.8736 and nucleotide diversity (π) 0.6474. These data clearly indicated that, feral/wild population showing highest values of polymorphisms, parsimony, haplotype diversity showing good, healthy habitat is lotic water (Narmada River) and lentic water body (Tighra reservoir). The results also concluded that the partial cyto b is polymorphic and can be a potential marker to determine ecological habitat based genetic differentiation among the populations.     

Phylogeography of the Golden Jackal (Canis aureus) in India

The golden jackal (Canis aureus) is one of the most common and widely distributed carnivores in India but phylogeographic studies on the species have been limited across its range. Recent studies have observed absence of mitochondrial (mt) DNA diversity in European populations while some North African populations of golden jackal were found to carry gray wolf (Canis lupus lupaster) mtDNA lineages. In the present study, we sequenced 440 basepairs (bp) of control region (CR) and 412 bp of cytochrome b (cyt b) gene of mtDNA from 62 golden jackals sampled from India (n = 55), Israel (n = 2) and Bulgaria (n = 5), to obtain a total of eighteen haplotypes, comprising sixteen from India and one each from Israel and Bulgaria. Except for three previously described haplotypes represented by one cyt b and one CR haplotype both from India, and one CR haplotype from Bulgaria, all haplotypes identified in this study are new. Genetic diversity was high in golden jackals compared to that reported for other canids in India. Unlike the paraphyletic status of African conspecifics with the gray wolf, the Indian (and other Eurasian) golden jackal clustered in a distinct but shallow monophyletic clade, displaying no evidence of admixture with sympatric and related gray wolf and domestic dog clades in the region. Phylogeographic analyses indicated no clear pattern of genetic structuring of the golden jackal haplotypes and the median joining network revealed a star-shaped polytomy indicative of recent expansion of the species from India. Indian haplotypes were observed to be interior and thus ancestral compared to haplotypes from Europe and Israel, which were peripheral and hence more derived. Molecular tests for demographic expansion confirmed a recent event of expansion of golden jackals in the Indian subcontinent, which can be traced back ~ 37,000 years ago during the late Pleistocene. Our results suggest that golden jackals have had a potentially longer evolutionary history in India than in other parts of the world, although further sampling from Africa, the Middle East and south-east Asia is needed to test this hypothesis.     

Population genetic structure in wild and hatchery populations of white cloud mountain minnow (Tanichthys albonubes): Recommendations for conservation

We used mitochondrial (mt) cytochrome b gene (cyt b) to compare the genetic variability in three hatchery broodstocks of white cloud mountain minnow with the variability in six wild populations sampled in two river drainages. A total of 43 haplotypes in 102 specimens were observed, with no haplotype shared between wild and hatchery populations. The nucleotide diversity of the wild samples (0.048) was significantly higher than that of the hatchery ones (0.007), but the haplotype diversity was almost similar between them. Two major phylogenetic haplotype groups were revealed and estimated to diverge about 6.531 myr (million years) ago. Significant genetic differentiation was revealed between wild and hatchery populations as well as among nine sampled populations, suggesting at chance effect during the founding process for the hatchery population and a subsequent genetic drift. According to the network, the connection between wild and hatchery populations indicates that present hatchery populations originated from single wild population. We suggested that two regions (Pearl River system and Lu River) identified by reciprocal mtDNA monophyly and SAMOVA should be regarded as three different ESUs and two different MUs in South China, respectively.     

祁连山区特有物种黑虎耳草的居群遗传多样性研究

利用叶绿体基因(rbcL和trnS-G)及核糖体DNA内转录间隔区(ITS),对祁连山区狭域分布的特有物种黑虎耳草8个居群(115个个体)进行遗传多样性研究,旨在揭示黑虎耳草的居群遗传结构及其历史进化过程。结果表明:(1)所有取样个体共检测到4个cpDNA单倍型和9个ITS单倍型,其中祁连山东南部的居群固定较多的单倍型和特有单倍型,而西北部居群只固定少数几个广泛分布的单倍型,且遗传多样性普遍较低。(2)基于cpDNA数据和ITS数据的分子变异分析(AMOVA)显示,遗传变异主要来源于居群内。(3)基于cpDNA数据的中性检验表明,Tajima’s D(-1.012 30,P> 0.05)和Fu&Li’s D*(-2.066 77,P> 0.05)均为负值,均不显著;歧点分布分析结果显示,黑虎耳草居群经历过近期的扩张事件。根据物种现有遗传分布格局推测,黑虎耳草在第四纪冰期时退缩到祁连山东南部的边缘避难所,间冰期或冰期后回迁到祁连山西北部地区,在回迁过程中由于奠基者效应导致祁连山西北部的居群仅固定少数广泛分布的单倍型,并呈现出较低的遗传多样性;由于居群较小且相互隔离,该...     

Genetic diversity and population structure of Schizopygopsis younghusbandi Regan in the Yarlung Tsangpo River inferred from mitochondrial DNA sequence analysis

The genetic diversity and population structure of Schizopygopsis younghusbandi from six sites in the middle reach of the Yarlung Tsangpo River were examined based on mitochondrial DNA cytochrome b (Cyt b) gene and control region (D-loop) sequences. 1141 bp of complete Cyt b sequences and 737 bp of partial D-loop sequences for 153 individuals of S. younghusbandi were obtained by using PCR amplification and sequencing. The results showed that S. younghusbandi populations had high haplotype diversity and low nucleotide diversity, and the population genetic diversity of this species was at a low level. Analysis of molecular variance revealed that most genetic variation occurred within populations, and that genetic differentiation was at low or moderate levels. The network of haplotypes based on Cyt b gene showed that there were two groups within the examined populations. Neutrality tests and mismatch distribution analysis suggested that this species might have undergone a population expansion and the expansion time was estimated to be 0.25–0.46 Ma BP. All these results would be crucial to establish scientific strategies for effective conservation and sustainable exploitation of wild resources of S. younghusbandi.     

Genetic variation in the mitochondrial 16S rRNA gene of the American dog tick,Dermacentor variabilis (Acari: Ixodidae)

Genetic variation in the mitochondrial (mt) 16S ribosomal RNA (rRNA) gene was examined for the American dog tick, Dermacentor variabilis (Say, 1821). Nine different haplotypes were detected among 369 adult D. variabilis collected from four localities in Canada. There were eight variable nucleotide positions in the 404 bp sequence alignment. Individuals of haplotype 1 occurred at frequency of >75% at all localities. Five haplotypes were detected at only one of the four localities. High haplotype diversity and low nucleotide diversity, combined with significantly negative F_s values for ticks at three localities, suggest a recent population expansion. Genetic differences were found between populations at different localities, but a Mantel regression analysis revealed no association between genetic differences and geographical distances. There was also no association between tick haplotype and the prevalence of the bacterium, Rickettsia montanensis Weiss and Moulder, 1984, in D. variabilis among localities or on opposite sides of Blackstrap Lake (Saskatchewan). The 16S rDNA haplotypes from Canadian populations of D. variabilis formed a clade with those from the eastern and central U.S.A., to the exclusion of D. variabilis from geographically isolated populations in the western U.S.A. Although sample sizes for D. variabilis in the eastern U.S.A. are small, there may be genetic divergence between populations in Canada and those in the eastern U.S.A., which may have implications for studies on the pathogenic agents transmitted by D. variabilis to its hosts.     

Genetic structure and diversity of the Diachasmimorpha longicaudata species complex in Thailand: SSCP analysis of mitochondrial 16S rDNA and COI DNA sequences

Genetic variation among 20 populations of Diachasmimorpha longicaudata (Ashmead) (Hymenoptera: Braconidae) in Thailand was investigated using single strand conformation polymorphism (SSCP) analysis of mitochondrial DNA sequences. From a total 641 individual parasitoids, seven distinct haplotypes containing a total of 32 polymorphic sites were observed from cytochrome oxidase subunit I (COI) sequences along with five distinct haplotypes containing a total 16 polymorphic sites from 16S rDNA sequences. Values obtained through pairwise F_ST comparisons and analysis of molecular variance (AMOVA) indicated significant genetic differentiation among D. longicaudata populations in Thailand. Congruent relationships showing separation of these populations into three groups were obtained from Neighbor joining and Bayesian phylogenetic tree analyses along with the use of haplotype networks. This SSCP analysis of populations of the D. longicaudata species complex is the first report using molecular population genetic methods to analyze the structure of this parasitoid species in Thailand. This may provide useful information for release of parasitoid strains to maximize their benefit in biological control programs.     

Chloroplast DNA Microsatellites Reveal Contrasting Phylogeographic Structure in Mahogany (Swietenia macrophylla King, Meliaceae) from Amazonia and Central America

Big-leaf mahogany (Swietenia macrophylla King) is one of the most valuable and overharvested timber trees of tropical America. In order to better characterize geographic patterns of genetic variation, we performed a phylogeographic analysis of S. macrophylla based on six polymorphic chloroplast genome simple sequence repeat loci (cpSSRs) analyzed in 16 populations (N?=?245 individuals) distributed across Central America and the Brazilian Amazon. Of the 31 total cpDNA haplotypes identified, 16 occurred in Central America and 15 in Amazonia with no single haplotype shared between the two regions. Populations from Central America showed moderate differentiation (F _ST ?=?0.36) while within population genetic diversity was generally high (mean Nei’s H _E ?=?0.639). In contrast, the Amazonian populations were strongly differentiated (F _ST ?=?0.91) and contained relatively low genetic diversity (mean H _E ?=?0.176), except for one highly diverse population (H _E ?=?0.925) from eastern Amazonia. Spatial analysis of molecular variance (SAMOVA) identified a single Central American phylogroup and four Amazonian phylogroups, indicating stronger phylogeographic structure within Amazonia. The results demonstrate distinctive regional patterns of S. macrophylla differentiation, and the first evidence of a strong phylogeographic break between Central American and South American mahogany populations. We suggest that the frequent occurrence of hurricanes in Central America, the differences in the glacial histories and in the duration and intensity of anthropogenic disturbance during the late Holocene may have played important roles in the geographic structuring of cpDNA lineages in the two regions. The high private haplotype diversity in Brazilian populations suggests that cpSSRs can be used as DNA barcodes for regional timber certification.     

Genetic and morphological analyses of tub gurnard Chelidonichthys lucerna populations in Turkish marine waters

Genetic and morphological structure of tub gurnard Chelidonichthys lucerna populations in Turkish marine waters were investigated with mtDNA sequencing of 16S rRNA and morphological characters. C. lucerna samples were collected from the Black Sea, Marmara, Aegean and northeastern Mediterranean coasts of Turkey. The lowest genetic diversity was found in the northeastern Mediterranean (Iskenderun Bay) population, while the highest was in the Marmara population with overall average value of genetic diversity within populations. A total of 14 haplotypes was found, and the highest haplotype diversity was in the Black Sea whereas the lowest was in the northeastern Mediterranean population (Iskenderun Bay). The Black Sea and Iskenderun Bay populations showed the least genetic divergence (0.001081), while the highest was between the Marmara Sea and northeastern Mediterranean (Antalya Bay) populations (0.002067). Pairwise comparisons of genetic distance revealed statistically significant differences (P < 0.05) between the Marmara and both the Aegean and northeastern Mediterranean (Antalya Bay) samples. Neighbour joining tree analyses clustered the northeastern Mediterranean populations (Antalya Bay and Iskenderun Bay) as genetically more interrelated populations, whereas the Aegean Sea population was clustered as most isolated one. Discriminant function analysis of morphological characters showed that only the Black Sea population is differentiated from the other populations.     

Genetic diversity of cultured populations of giant freshwater prawn (Macrobrachium rosenbergii) in China using mtDNA COI and 16S rDNA markers

Partial sequences of mitochondrial DNA 16S rDNA and COI genes (395 bp and 498 bp respectively) were sequenced from samples of ten cultured populations of Macrobrachium rosenbergii (Giant Freshwater Prawn – GFP) in Zhejiang, Guangdong and Guangxi Provinces in China and two wild populations of GFP from Mekong River and Dongnai River in Viet Nam. Five haplotypes of 16S rDNA were identified in the 360 samples. The wild populations displayed nucleotide diversity (π) of 0.0008 and 0.0003, and genetic diversity (h) of 0.3030 and 0.1310 in the Mekong River and Dongnai River respectively. The cultured populations displayed no significant genetic diversity. COI sequences identified 17 haplotypes based on 21 polymorphic sites. At this marker, the 12 populations showed a range of h from 0.1290 to 0.6940 and π from 0.0003 to 0.0073. The largest genetic distance (Da) among the 12 populations was 0.0065 (between ZJB and BT/DN populations) and the lowest Da was 0.0003 (between GDD and GDA populations). The wild populations had higher genetic diversity than the cultured populations, but three of cultured populations from Zhejiang (ZJA, ZJB and ZJC) had π higher than wild populations, because they originated from Thailand, Bangladesh and the Mekong River in Viet Nam.     

Mitochondrial sequence polymorphism in Chinese populations of <Emphasis Type="Italic">Hylobitelus xiaoi</Emphasis> (Coleoptera: Curculionidae)

The Chinese weevil, Hylobitelus xiaoi Zhang, is a major pest affecting the pine industry throughout southern China, but its dispersal is still poorly understood. We aimed to investigate its dispersal by assessing the genetic structure of seven geographical populations, using the distribution of two mitochondrial genes (cytochrome c oxidase subunit I and cytochrome b). A 1212-bp fragment of the two mitochondrial genes was sequenced for 70 individuals from seven geographical populations. The gene sequences included 861 conserved sites, 351 variable sites, 254 parsim-info sites, and 97 singleton sites. These polymorphic sites defined 24 haplotypes with a haplotype diversity of 0.825. Nucleotide diversity was low (0.04568). The differentiation parameter (0.619) was much greater than the coefficient of genetic differentiation (0.285). An analysis of molecular variance suggested that most of the variation was due to within-population differences (61.68%). This molecular data also demonstrated that there is significant divergence among the seven sampled populations of H. xiaoi. This is probably due to habitat fragmentation, preventing effective gene flow between sites. Our molecular evidence supports the view that the Chinese weevil is a native pest, and that outbreaks of H. xiaoi are caused by large plantations of the exotic slash pine.     

No genetic differentiation between geographically isolated populations of Clarias macrocephalus Günther in Malaysia revealed by sequences of mtDNA Cytochrome b and D‐loop gene regions

In the present study, we assessed the genetic variation of three Clarias macrocephalus Günther populations collected from Kedah, Perlis and Kelantan (Peninsular Malaysia) using sequences of partial mitochondrial cytochrome b (Cyt b) and D‐loop genes. A total of 57 individuals were sequenced and 1470 bp were obtained (1053 bp Cyt‐b; 417 bp D‐loop). The analysis revealed 21 haplotypes based on 81 polymorphic sites. Nucleotide diversity (π) was 0.003 in all populations while haplotype diversity ranged from 0.657 to 0.765. No significant genetic differentiation among the three populations was observed. Nevertheless, a number of private haplotypes was discovered, providing valuable information for selective breeding programs.     

Mitochondrial DNA variation in natural populations of endangered Indian Feather-Back Fish, Chitala chitala

Genetic variation at mitochondrial cytochrome b (cyt b) and D-loop region reveals the evidence of population sub-structuring in Indian populations of highly endangered primitive feather-back fish Chitala chitala. Samples collected through commercial catches from eight riverine populations from different geographical locations of India were analyzed for cyt b region (307 bp) and D-loop region (636–716 bp). The sequences of the both the mitochondrial regions revealed high haplotype diversity and low nucleotide diversity. The patterns of genetic diversity, haplotypes networks clearly indicated two distinct mitochondrial lineages and mismatch distribution strongly suggest a historical influence on the genetic structure of C. chitala populations. The baseline information on genetic variation and the evidence of population sub-structuring generated from this study would be useful for planning effective strategies for conservation and rehabilitation of this highly endangered species.     

Relationships among arctic and non-arctic haplotypes of the moss species Scorpidium cossonii and Scorpidium scorpioides (Calliergonaceae)

Relationships among Scorpidium cossonii and Scorpidium scorpioides haplotypes from most of the species’ distribution areas were analyzed based on ITS and rpl16. The haplotype networks were produced by TCS and were rooted by neighbor joining (ITS, recombination present) or maximum parsimony analysis. The haplotype closest to the root of the ITS network and some poorly represented haplotypes close to this have S. cossonii morphology and are arctic to subarctic, suggesting a northern origin of the complex. Additionally, two major lineages evolved from the root haplotype; one with S. cossonii morphology samples and one with all S. scorpioides morphology ones. The basal haplotypes in these major lineages include numerous temperate zone representatives, suggesting that adaptations for relatively warmer environments are present. Among haplotype groups that evolved later, one group in each major lineage is absent from America. rpl16 provides a scenario similar to that revealed by ITS for S. cossonii (no variation in S. scorpioides). Fossil and molecular evidence suggest that ancestral populations of S. cossonii evolved at least almost 4 Myr bp.     

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

为了解中国狼不同地理种群遗传多样性及系统发育情况,从中国境内狼的主要分布区青海、新疆、内蒙古和吉林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前。     

The Population Genetics of Cultivation: Domestication of a Traditional Chinese Medicine,Scrophularia ningpoensis Hemsl. (Scrophulariaceae)

Background

Domestic cultivation of medicinal plants is an important strategy for protecting these species from over harvesting. Some species of medicinal plants have been brought into cultivation for more than hundreds years. Concerns about severe loss of genetic diversity and sustainable cultivation can potentially limit future use of these valuable plants. Genetic studies with comprehensive sampling of multiple medicinal species by molecular markers will allow for assessment and management of these species. Here we examine the population genetic consequences of cultivation and domestication in Scrophularia ningpoensis Hemsl. We used chloroplast DNA and genomic AFLP markers to clarify not only the effects of domestication on genetic diversity, but also determine the geographic origins of cultivars and their genetic divergence from native populations. These results will allow both better management of cultivated populations, but also provide insights for crop improvement.

Results

Twenty-one cpDNA haplotypes of S. ningpoensis were identified. Wild populations contain all haplotypes, whereas only three haplotypes were found in cultivated populations with wild populations having twice the haplotype diversity of cultivated populations. Genetic differentiation between cultivated populations and wild populations was significant. Genomic AFLP markers revealed similar genetic diversity patterns. Furthermore, Structure analysis grouped all wild populations into two gene pools; two of which shared the same gene pool with cultivated S. ningpoensis. The result of Neighbor-Joining analysis was consistent with the structure analysis. In principal coordinate analysis, three cultivated populations from Zhejiang Province grouped together and were separated from other cultivated populations.

Conclusions

These results suggest that cultivated S. ningpoensis has experienced dramatic loss of genetic diversity under anthropogenic influence. We postulate that strong artificial selection for medicinal quality has resulted in genetic differentiation between cultivated and wild populations. Furthermore, it appears that wild populations in Jiangxi-Hunan area were involved in the origin of cultivated S. ningpoensis.     

Population genetic diversity of Chinese sea bass (Lateolabrax maculatus) from southeast coastal regions of China based on mitochondrial COI gene sequences

Chinese sea bass (Lateolabrax maculatus), an important commercial and recreational fishery species in China, whose natural resources have decreased dramatically for decades due to overfishing and environmental changes. To provide guidelines for the source conservation and management, genetic diversity and population structure of L. maculatus were analyzed based on the mitochondrial cytochrome oxidase subunit I (COI) gene. Totally 192 individuals were sampled from five locations: Qingdao (QD), Chongming (CM), Dongtou (DT), Lieyu (LY) and Fangcheng (FC). After sequencing of a 586 bp fragment of COI gene, 20 haplotypes were defined. H4 (haplotype 4) and H8 were the dominant haplotypes and existed in all populations. Haplotype diversity (Hd) and nucleotide diversity (π) of each population ranged from 0.572 to 0.721 and from 0.00129 to 0.00271, respectively. The highest Hd and π were found in DT and LY populations, while the lowest value of Hd and π were in the QD population. Genetic distance ranged from 0.0031 to 0.0027 within populations and from 0.0018 to 0.0035 between populations. Analysis of molecular variance (AMOVA) revealed that significant genetic divergence was found in QD and FC populations. In addition, neutrality tests and mismatch distribution analysis indicated that this species experienced potential population expansion events.     

Geographic differentiation of <Emphasis Type="Italic">Cobitis shikokuensis</Emphasis> inferred from mtDNA RFLP analysis

Mitochondrial DNA divergence among populations of the Japanese spinous loach Cobitis shikokuensis, endemic to Shikoku Island, was investigated by restricted fragment length polymorphism analysis. A total of 68 restriction sites on DNA fragments from the cytochrome b to D-loop regions and from the 12S rRNA to 16S rRNA regions, amplified by PCR, were analyzed. A total of 12 haplotypes (plus 6 in outgroups) were detected in 268 specimens collected from 19 localities in seven rivers (and 41 specimens from four localities in three rivers in outgroups). Three of the seven river populations of C. shikokuensis were shown to have unique haplotypes, and four of the seven river populations were monomorphic. The nested structure of the haplotype network for populations of C. shikokuensis exhibited two large clades corresponding to (1) populations from the Shimanto River and its neighbors and (2) two genetically divergent populations in the Shigenobu and Iwamatsu Rivers. The population from the Shimanto River, the largest river inhabited by C. shikokuensis, maintains great haplotype diversity as well as the allozyme diversity previously reported. On the other hand, populations from the Hiji River, the second largest river, which exhibited the highest allozyme diversity, were monomorphic in their mtDNA. The nested clade analysis (NCA) revealed that past fragmentation between the above two clades could occur in the initial distribution process of C. shikokuensis. The large genetic divergence of two river populations from the Shigenobu and Iwamatsu Rivers was inferred to be caused by a process of long distance colonization and fragmentation. MtDNA introgression into the Hiji River population from southern river populations was suggested. Taking genetic divergence into consideration, each river population of C. shikokuensis should be conserved separately as like a distinct species, and conservation programs for the small populations showing less genetic variability should be invoked as soon as possible.