Phylogeography and conservation genetics of Eld's deer (Cervus eldi)

Eld's deer (Cervus eldi) is a highly endangered cervid, distributed historically throughout much of South Asia and Indochina. We analysed variation in the mitochondrial DNA (mtDNA) control region for representatives of all three Eld's deer subspecies to gain a better understanding of the genetic population structure and evolutionary history of this species. A phylogeny of mtDNA haplotypes indicates that the critically endangered and ecologically divergent C. eldi eldi is related more closely to C. e. thamin than to C. e. siamensis, a result that is consistent with biogeographic considerations. The results also suggest a strong degree of phylogeographic structure both between subspecies and among populations within subspecies, suggesting that dispersal of individuals between populations has been very limited historically. Haplotype diversity was relatively high for two of the three subspecies (thamin and siamensis), indicating that recent population declines have not yet substantially eroded genetic diversity. In contrast, we found no haplotype variation within C. eldi eldi or the Hainan Island population of C. eldi siamensis, two populations which are known to have suffered severe population bottlenecks. We also compared levels of haplotype and nucleotide diversity in an unmanaged captive population, a managed captive population and a relatively healthy wild population. Diversity indices were higher in the latter two, suggesting the efficacy of well-designed breeding programmes for maintaining genetic diversity in captivity. Based on significant genetic differentiation among Eld's deer subspecies, we recommend the continued management of this species in three distinct evolutionarily significant units (ESUs). Where possible, it may be advisable to translocate individuals between isolated populations within a subspecies to maintain levels of genetic variation in remaining Eld's deer populations.     

植物保护遗传学

从遗传多样性、系统发育和地理系统发育角度对当前保护遗传学在植物保护中的作用进行了探讨。在植物保护过程中,要考虑多群遗传多样性的大小,在就地和迁地保护的过程要减少近交和远交衰退影响,并可利用遗传标记为提供关于种群大小、基因流动等方面的信息。系统发育和地理系统发育的研究在于了解物种进化的历史,以确定物种保护单元。同时结合两个具体生态学问题（生境片段化和外来种）对当前植物保护遗传学的研究进行了介绍,揭示了保护遗传学在植物保护上重要作用和不可取代性,为今后研究提供参考。     

Phylogeography of the Angolan black and white colobus monkey,Colobus angolesnsis palliatus,in Kenya and Tanzania

Little is known about genetic variation in the 6–8 subspecies of Colobus angolensis, currently distinguished by pelage differences. We present a comparative genetic analysis of one of these subspecies, C. a. palliatus, in Kenya and Tanzania that assesses evolutionary relationships and patterns of mitochondrial genetic diversity in 103 individuals across its geographic range. Fecal samples from approximately 156 individuals were collected in four localities: (1) Diani Forest, Kenya; (2) Shimoni, Kenya; (3) Udzungwa Mountains National Park, Eastern Arc Mountains, Tanzania; and (4) Mount Rungwe, Southern Highlands, Tanzania. These samples represent at least six groups, with 5–15 samples from each. Comparative sequence analysis of a 1,795 base pair mtDNA fragment revealed 19 unique haplotypes in four populations. Phylogenetic analyses suggest that sampled Kenyan haplotypes are paraphyletic, with one Kenyan haplotype basal to all other sampled haplotypes. Analysis of molecular variance (AMOVA) suggests high levels of genetic variation among populations (Φ_ST 0.72, P<0.001). Genetic data are concordant with a subspecies level differentiation between C. a. palliatus populations in Kenya and those in Central and southern Tanzania, as earlier suggested based on pelage differences. This study highlights the evolutionary distinctiveness of Kenyan populations of C. a. palliatus relative to Tanzanian populations. Although C. a. palliatus habitat in Tanzania is currently better protected than in Kenya, our results suggest Kenyan and Tanzanian populations should be considered distinct units, and the protection of C. a. palliatus habitat in Kenya, as well as habitat connectivity between Kenyan populations, should be prioritized for conservation and management. Am. J. Primatol. 72:715–724, 2010. © 2010 Wiley‐Liss, Inc.     

Relationships of Palearctic and Nearctic 'glacial relict' Myoxocephalus sculpins from mitochondrial DNA data

The relationships among Myoxocephalus quadricornis complex fish from Arctic coastal waters and from 'glacial relict' populations in Nearctic and Palearctic postglacial lakes were assessed using mtDNA sequence data (1978 bp). A principal phylogeographical split separated the North American continental deepwater sculpin (M. q. thompsonii) from a lineage of the Arctic marine and North European landlocked populations of the fourhorn sculpin (M. q. quadricornis). The North American continental invasion took place several glaciation cycles ago in the Early-to-Middle Pleistocene (0.9% sequence divergence); the divergence of the European and Arctic populations was somewhat later (0.5% divergence). The Nearctic-Palearctic freshwater vicariance in Myoxocephalus, however, appears clearly younger than in similarly distributed 'glacial relict' crustacean taxa; the phylogeographical structure is more similar to that in other northern Holarctic freshwater fish complexes.     

Haplotype Frequency Distribution in Northeastern European Saduria entomon (Crustacea: Isopoda) Populations. A Phylogeographic Approach

The distribution pattern of mtDNA haplotypes in distinct populations of the glacial relict crustacean Saduria entomon was examined to assess phylogeographic relationships among them. Populations from the Baltic, the White Sea and the Barents Sea were screened for mtDNA variation using PCR‐based RFLP analysis of a 1150 bp fragment containing part of the CO I and CO II genes. Five mtDNA haplotypes were recorded. An analysis of geographical heterogeneity in haplotype frequency distributions revealed significant differences among populations. The isolated populations of S. entomon have diverged since the retreat of the last glaciation. The geographical pattern of variation is most likely the result of stochastic (founder effect, genetic drift) mechanisms and suggests that the haplotype differentiation observed is probably older than the isolation of the Baltic and Arctic seas.     

Phylogeography, genetic structure and diversity in the endangered bearded vulture (Gypaetus barbatus, L) as revealed by mitochondrial DNA

Bearded vulture populations in the Western Palearctic have experienced a severe decline during the last two centuries that has led to the near extinction of the species in Europe. In this study we analyse the sequence variation at the mitochondrial control region throughout the species range to infer its recent evolutionary history and to evaluate the current genetic status of the species. This study became possible through the extensive use of museum specimens to study populations now extinct. Phylogenetic analysis revealed the existence of two divergent mitochondrial lineages, lineage A occurring mainly in Western European populations and lineage B in African, Eastern European and Central Asian populations. The relative frequencies of haplotypes belonging to each lineage in the different populations show a steep East-West clinal distribution with maximal mixture of the two lineages in the Alps and Greece populations. A genealogical signature for population growth was found for lineage B, but not for lineage A; futhermore the Clade B haplotypes in western populations and clade A haplo-types in eastern populations are recently derived, as revealed by their peripheral location in median-joining haplotype networks. This phylogeographical pattern suggests allopatric differentiation of the two lineages in separate Mediterranean and African or Asian glacial refugia, followed by range expansion from the latter leading to two secondary contact suture zones in Central Europe and North Africa. High levels of among-population differentiation were observed, although these were not correlated with geographical distance. Due to the marked genetic structure, extinction of Central European populations in the last century re-sulted in the loss of a major portion of the genetic diversity of the species. We also found direct evidence for the effect of drift altering the genetic composition of the remnant Pyrenean population after the demographic bottleneck of the last century. Our results argue for the management of the species as a single population, given the apparent ecological exchangeability of extant stocks, and support the ongoing reintroduction of mixed ancestry birds in the Alps and planned reintroductions in Southern Spain.     

Phylogeography of isolated freshwater three-spined stickleback Gasterosteus aculeatus populations in the Adriatic Sea basin

Analyses of mitochondrial (mt) DNA and microsatellite variation were carried out to examine the relationships between 10 freshwater populations of three-spined sticklebacks Gasterosteus aculeatus along the eastern coast of the Adriatic Sea. Partial sequences of the mtDNA control region and cytochrome b gene, in addition to 15 microsatellite loci, were used to analyse populations from four isolated river catchments. Results uncovered an Adriatic lineage that was clearly divergent from the European lineage, and confirmed that the most divergent and ancient populations are located within the Adriatic lineage as compared with other European populations. Two northern Adriatic populations formed independent clades within the European mitochondrial lineage, suggesting different colonization histories of the different Adriatic populations. Nuclear marker analyses also indicated deep divergence between Adriatic and European populations, albeit with some discordance between the mtDNA phylogeny of the northern Adriatic populations, further highlighting the strong differentiation among the Adriatic populations. The southern populations within the Adriatic lineage were further organized into distinct clades corresponding to respective river catchments and sub-clades corresponding to river tributaries, reflecting a high degree of population structuring within a small geographic region, concurrent with suggestions of existence of several microrefugia within the Balkan Peninsula. The highly divergent clades and haplotypes unique to the southern Adriatic populations further suggest, in accordance with an earlier, more limited survey, that southern Adriatic populations represent an important reservoir for ancient genetic diversity of G. aculeatus.     

报春花属的保育遗传学及亲缘地理学研究进展

现存物种的分布格局、遗传结构是当前因素和历史因素共同作用的结果。人类的过度开发、全球变暖等当前因素造成的生境片段化是目前许多报春花属植物濒危的一大原因。该文总结了近年来报春花属内的保育遗传学研究进展,期望以此为更好地保护报春花属植物提供一定的理论基础。应用亲缘地理学研究方法可以弥补古生物学难以研究报春花植物历史成因的不足,因此也总结该方法在报春花属内的研究进展,并初步整理不同地区间的报春花属植物的分化式样,同时期望这些研究成果能为为研究报春花属植物在应对全球气候变化的响应机制方面提供一些参考。     

Phylogeography and conservation genetics of the Columbia spotted frog (Rana luteiventris; Amphibia, Ranidae)

The Columbia spotted frog (Rana luteiventris) has a widespread distribution in western Canada and the western US, although the southern reach of its range is highly fragmented into several isolated populations. Threats from various factors have raised concerns regarding the long-term survival of many small, isolated populations. Here, we report a study designed to determine the phylogeographic and conservation genetic parameters of R. luteiventris in the western US. Mitochondrial DNA (mtDNA) sequences were examined for phylogeographic structuring using phylogenetic reconstruction methods, coupled with networking and nested clade analyses. These methods permitted a distinction to be made between historic and demographic forces acting to generate geographical patterning of genetic variation. Phylogenetic analysis revealed four geographically correlated monophyletic clades. Three of these clades correspond to well-defined, nonoverlapping geographical locations in the fragmented portion of the range. The other is comprised of all samples collected from the contiguous range and includes one isolate from northern Wyoming. Networking and nested clade analyses confirmed these results and revealed that historical processes, such as range expansion and vicariance, rather than recurrent gene flow are likely responsible for observed patterns of genetic variation. A measure of genetic variation (theta = 4N(e)mu) revealed that R. luteiventris populations in Utah have a relatively low amount of genetic variation compared with populations in the continuous portion of the range.     

Phylogeography of the endangered darkling beetle species of Pimelia endemic to Gran Canaria (Canary Islands)

Phylogenetic and geographical nested clade analysis (NCA) methods were applied to mitochondrial DNA sequences of Pimelia darkling beetles (Coleoptera, Tenebrionidae) endemic to Gran Canaria, an island in the Canary archipelago. The three species P. granulicollis, P. estevezi and P. sparsa occur on the island, the latter with three recognized subspecies. Another species, P. fernandezlopezi (endemic to the island of La Gomera) is a close relative of P. granulicollis based on partial Cytochrome Oxidase I mtDNA sequences obtained in a previous study. Some of these beetles are endangered, so phylogeographical structure within species and populations can help to define conservation priorities. A total of about 700 bp of Cytochrome Oxidase II were examined in 18 populations and up to 75 individuals excluding outgroups. Among them, 22 haplotypes were exclusive to P. granulicollis and P. estevezi and 31 were from P. sparsa. Phylogenetic analysis points to the paraphyly of Gran Canarian Pimelia, as the La Gomera P. fernandezlopezi haplotypes are included in them, and reciprocal monophyly of two species groups: one constituted by P. granulicollis, P. estevezi and P. fernandezlopezi (subgenus Aphanaspis), and the other by P. sparsa'sensu lato'. The two species groups show a remarkably high mtDNA divergence. Within P. sparsa, different analyses all reveal a common result, i.e. conflict between current subspecific taxonomic designations and evolutionary units, while P. estevezi and P. fernandezlopezi are very close to P. granulicollis measured at the mtDNA level. Geographical NCA identifies several cases of nonrandom associations between haplotypes and geography that may be caused by allopatric fragmentation of populations with some cases of restriction of gene flow or range expansion. Analyses of molecular variance and geographical NCA allow definition of evolutionary units for conservation purposes in both species-groups and suggest scenarios in which vicariance caused by geological history of the island may have shaped the pattern of the mitochondrial genetic diversity of these beetles.     

Phylogeography, population history and conservation genetics of jaguars (Panthera onca, Mammalia, Felidae)

The jaguar (Panthera onca), the largest felid in the American Continent, is currently threatened by habitat loss, fragmentation and human persecution. We have investigated the genetic diversity, population structure and demographic history of jaguars across their geographical range by analysing 715 base pairs of the mitochondrial DNA (mtDNA) control region and 29 microsatellite loci in approximately 40 individuals sampled from Mexico to southern Brazil. Jaguars display low to moderate levels of mtDNA diversity and medium to high levels of microsatellite size variation, and show evidence of a recent demographic expansion. We estimate that extant jaguar mtDNA lineages arose 280 000-510 000 years ago (95% CI 137 000-830 000 years ago), a younger date than suggested by available fossil data. No strong geographical structure was observed, in contrast to previously proposed subspecific partitions. However, major geographical barriers such as the Amazon river and the Darien straits between northern South America and Central America appear to have restricted historical gene flow in this species, producing measurable genetic differentiation. Jaguars could be divided into four incompletely isolated phylogeographic groups, and further sampling may reveal a finer pattern of subdivision or isolation by distance on a regional level. Operational conservation units for this species can be defined on a biome or ecosystem scale, but should take into account the historical barriers to dispersal identified here. Conservation strategies for jaguars should aim to maintain high levels of gene flow over broad geographical areas, possibly through active management of disconnected populations on a regional scale.     

Phylogeography of the pitviper clade Agkistrodon: historical ecology, species status, and conservation of cantils

We used mitochondrial DNA sequences from three gene regions and two tRNAs (ND4, tRNA-HIS-SER, 12S, and 16S rDNA) to investigate the historical ecology of the New World pitviper clade Agkistrodon, with emphasis on the disjunct subspecies of the cantil, A. bilineatus. We found strong evidence that the copperhead (A. contortrix) is basal to its congeners, and that the cottonmouth (A. piscivorus) is basal to cantils. Phylogeography and natural history of the living terminal taxa imply that Agkistrodon primitively occupied relatively temperate habitats, with subsequent evolution of tropicality in ancestral A. bilineatus. Our best supported phylogeny rejects three gulf arc scenarios for the biogeography of A. bilineatus. We find significant statistical support for an initial divergence between populations on the east and west coasts of México and subsequent occupancy of the Yucatán Peninsula, by way of subhumid corridors in northern Central America. Based on phylogenetic relationships, morphological and molecular divergence, and allopatry we elevate A. b. taylori of northeastern México to species status. Taylor's cantil is likely threatened by habitat destruction and small geographical range, and we offer recommendations for its conservation and management.     

Molecular phylogeography and conservation genetics of Sladenia celastrifolia inferred from chloroplast DNA sequence variation

The genetic variation and structure of Sladenia celastrifolia Kurz, a species of conservation concern, were investigated. Analyses of two chloroplast DNA loci (trnS-trnGand atpB-rbcL intergenic regions) were carried out for 24 populations of S. celastrifolia and five haplotypes were identified. High levels of genetic differentiation (G_ST = 1, F_ST = 1) were detected, which may be a result of limited gene flow caused by geographic isolation. Analysis of molecular variance suggests that the existence of marked phylogeographical structure within the haplotype distribution is probably due to geographic barriers among populations. The haplotype network and mismatch distribution analyses did not detect any signals for recent population expansions in S. celastrifolia. Thus, it can be inferred that the species likely persistedin situ during climatic oscillations. Considering its genetic diversity and uniqueness, conservation strategies are further discussed for this species.     

Conservation genetics and population history of the threatened European mink Mustela lutreola, with an emphasis on the west European population

In species of great conservation concern, special attention must be paid to their phylogeography, in particular the origin of animals for captive breeding and reintroduction. The endangered European mink lives now in at least three well-separated populations in northeast, southeast and west Europe. Our aim is to assess the genetic structure of these populations to identify 'distinct population segments' (DPS) and advise captive breeding programmes. First, the mtDNA control region was completely sequenced in 176 minks and 10 polecats. The analysis revealed that the western population is characterized by a single mtDNA haplotype that is closely related to those in eastern regions but nevertheless, not found there to date. The northeast European animals are much more variable (pi = 0.012, h = 0.939), with the southeast samples intermediate (pi = 0.0012, h = 0.469). Second, 155 European mink were genotyped using six microsatellites. The latter display the same trends of genetic diversity among regions as mtDNA [gene diversity and allelic richness highest in northeast Europe (H(E) = 0.539, R(S) = 3.76), lowest in west Europe (H(E) = 0.379, R(S) = 2.12)], and provide evidences that the southeast and possibly the west populations have undergone a recent bottleneck. Our results indicate that the western population derives from a few animals which recently colonized this region, possibly after a human introduction. Microsatellite data also reveal that isolation by distance occurs in the western population, causing some inbreeding because related individuals mate. As genetic data indicate that the three populations have not undergone independent evolutionary histories for long (no phylogeographical structure), they should not be considered as distinct DPS. In conclusion, the captive breeding programme should use animals from different parts of the species' present distribution area.     

Global population structure of the tope (Galeorhinus galeus) inferred by mitochondrial control region sequence data

In order to properly manage and conserve exploited shark species, detailed analyses of their population structure is needed. Global populations of Galeorhinus galeus are in decline due to the exploitation of the fishery over the past 80 years. Currently, the genetic structure of eastern Pacific populations of G. galeus is not known and recent observations in the northeastern Pacific suggest an increase in numbers. To evaluate gene flow among populations of G. galeus , 116 samples were collected and analysed from six geographically dispersed locations: Australia, North America, South Africa, South America (Argentina and Peru), and the UK. Analysis of 968 to 1006 bp of the 1068-bp mitochondrial control region revealed 38 unique haplotypes that were largely restricted to their collecting locality. Significant genetic structure was detected among populations (Φ_ST = 0.84; P  < 0.000001) and migration estimates were low ( Nm  = 0.05–0.97). Due to an apparent lack of migration, populations of G. galeus appear to be isolated from each other with little to no gene flow occurring among them. As a consequence of this isolation, increasing numbers of G. galeus in the northeastern Pacific can be best explained by local recruitment and not by input from geographically distant populations.     

Phylogeography of olive ridley turtles (Lepidochelys olivacea) on the east coast of India: implications for conservation theory

Orissa, on the east coast of India, is one of the three mass nesting sites in the world for olive ridley turtles (Lepidochelys olivacea). This population is currently under threat as a result of fishery-related mortality; more than 100 000 olive ridleys have been counted dead in the last 10 years in Orissa. In general, the globally distributed olive ridley turtle has received significantly less conservation attention than its congener, the Kemp's ridley turtle (L. kempi), because the latter is recognized as a distinct species consisting of a single endangered population. Our study of mitochondrial DNA haplotypes suggests that the ridley population on the east coast of India is panmictic, but distinct from all other populations including Sri Lanka. About 96% of the Indian population consisted of a distinct 'K' clade with haplotypes not found in any other population. Nested clade analysis and conventional analysis both supported range expansions and/or long-distance colonization from the Indian Ocean clades to other oceanic basins, which suggested that these are the ancestral source for contemporary global populations of olive ridley turtles. These data support the distinctiveness of the Indian Ocean ridleys, suggesting that conservation prioritization should be based on appropriate data and not solely on species designations.     

Phylogeography of the greater horseshoe bat, Rhinolophus ferrumequinum: contrasting results from mitochondrial and microsatellite data

Phylogeographical studies are typically based on haplotype data, occasionally on nuclear markers such as microsatellites, but rarely combine both. This is unfortunate because the use of markers with contrasting modes of inheritance and rates of evolution might provide a more accurate and comprehensive understanding of a species' history. Here we present a detailed study of the phylogeography of the greater horseshoe bat, Rhinolophus ferrumequinum , using 1098 bp of the mitochondrial ND2 gene from 45 localities from across its Palaearctic range to infer population history. In addition, we re-analysed a large microsatellite data set available for this species and compared the results of both markers to infer population relationships and the historical processes influencing them. We show that mtDNA, the most popular marker in phylogeography studies, yielded a misleading result, and would have led us to conclude erroneously that a single expansion had taken place in Europe. Only by combining the mitochondrial and microsatellite data sets are we able to reconstruct the species' history and show two colonization events in Europe, one before the Last Glacial Maximum (LGM) and one after it. Combining markers also revealed the importance of Asia Minor as an ancient refugium for this species and a source population for the expansion of the greater horseshoe bat into Europe before the LGM.     

Phylogeography and conservation genetics of Lake Qinghai scaleless carp Gymnocypris przewalskii

The objective of this study was to examine the spatial genetic relationships of the Lake Qinghai scaleless carp Gymnocypris przewalskii within the Lake Qinghai system, determining whether genetic evidence supports the current taxonomy of Gymnocypris przewalskii przewalskii and Gymnocypris przewalskii ganzihonensis and whether Gymnocypris przewalskii przewalskii are returning to their natal rivers to spawn. Comparison of mitochondrial (control region) variation (42 haplotypes in 203 fish) of G. przewalskii with the postulated ancestral species found in the Yellow River, Gymnocypris eckloni (10 haplotypes in 23 fish), indicated no haplotype sharing, but incomplete lineage sorting. Consistent with the sub-species status, an AMOVA indicated that the Ganzi River population was significantly different from all other river populations (F(ST) = 0·1671, P < 0·001). No genetic structure was found among the other rivers in the Lake Qinghai catchment. An AMOVA of amplified fragment length polymorphism (AFLP) loci, however, revealed significant genetic differences between most spawning populations (F(ST) = 0·0721, P < 0·001). Both mitochondrial and AFLP data found significant differences among G. p. przewalskii, G. p. ganzihonensis and G. eckloni (F(ST) values of 0·1959 and 0·1431, respectively, P < 0·001). Consistent with the incomplete lineage sorting, Structure analysis of AFLP loci showed evidence of five clusters. One cluster is shared among all sample locations, one is unique to G. p. ganzihonensis and G. eckloni, and the others are mostly found in G. p. przewalskii. Genetic evidence therefore supports the current taxonomy, including the sub-species status of G. p. ganzihonensis, and is consistent with natal homing of most Lake Qinghai populations. These findings have significant implications for the conservation and management of this unique and threatened species. The evidence suggests that G. p. przewalskii should be treated as a single population for conservation purposes. Exchangeability of the populations, however, should not be used to promote homogenization of fish spawning in the different rivers. As some degree of genetic divergence was detected in this study, it is recommended that the spawning groups be treated as separate management units.     

Phylogeography of the common goby, Pomatoschistus microps, with particular emphasis on the colonization of the Mediterranean and the North Sea

The phylogeographical patterns of a small marine fish, the common goby, Pomatoschistus microps, were assessed at 12 sites along the northeastern Atlantic coasts and the western Mediterranean Sea. A combination of two genetic markers was employed: cellulose acetate allozyme electrophoresis (CAGE) and sequence analysis of a 289 bp fragment of the mitochondrial locus cytochrome b. Both markers were congruent in revealing significant differences between samples (global FST = 0.247 for the allozymes and PhiST = 0.437 for the mitochondrial DNA data) and a pattern of isolation-by-distance. Phylogeographical analyses yielded a shallow branching structure with four groups. Three of those were confined to the Atlantic basin and showed a star-like pattern. The fourth group contained a central haplotype occurring at the edges of the species' distribution, accompanied by a few more rare variants, which were restricted to the Mediterranean Sea. A genetic break was observed around the British Isles, with distinct haplotypes dominating at either side of the English Channel. A significantly negative correlation between the degree of genetic diversity and latitude was recorded both for mitochondrial DNA (mtDNA) and allozymes in the Atlantic basin. Gene flow analysis suggested that recolonization of the North Sea and the coasts of western Scotland and Ireland may have taken place from a glacial refugium in the Southern Bight of the North Sea. These results are discussed in the perspective of possible postglacial migration routes of marine fish along the northeastern Atlantic coasts.     

Phylogeography of silver‐studded blue,Plebejus subsolanus (Lepidoptera: Lycaenidae), in the central mountainous regions of Japan

The silver‐studded blue, Plebejus subsolanus, is widely distributed in the Russian Altai mountains, northeastern China, the Korean Peninsula, and the Japanese archipelago. In Japan, the species is distributed across wide elevation ranges from the lowlands of Hokkaido to the subalpine zone of Honshu. Current subspecies classification in Japan is as follows: ssp. iburiensis, occurring in lowland grasslands in Hokkaido; ssp. yaginus in lower mountain grasslands in Honshu; and ssp. yarigadakeanus in higher mountain grasslands in Honshu. The habitat of this species has been markedly reduced due to recent habitat destruction and land‐use changes. Here, we undertook phylogeographic analyses of two subspecies, ssp. yaginus and yarigadakeanus in the central mountainous regions of Japan, based on two mitochondrial gene sequences, in order to collect information for establishing effective conservation strategies. From 57 samples from the four mountain ranges, we obtained a haplotype network comprised of 12 haplotypes. Because of the haplotype network topology, the geographic distribution of haplotypes and the correspondence of haplotype divergence to subspecies taxonomy, we provisionally divided the haplotypes into three haplogroups: YR1 and YR2, which comprised ssp. yarigadakeanus, and YG, which comprised ssp. yaginus. Mitochondrial DNA genetic differentiation generally agreed with morphological subspecies classification. The haplotype network suggested that ssp. yarigadakeanus populations had multiple origins, and the subspecies character of “bright blue of the male's wings” was assumed to have evolved independently in each subalpine meadow. We found that P. subsolanus was genetically differentiated depending upon the elevation at each mountain region, suggesting that each haplogroup should be a conservation unit.