Phylogeographic sympatry and isolation of the Eurasian badgers (Meles, Mustelidae, Carnivora): Implications for an alternative analysis using maternally as well as paternally inherited genes

In the present study, to further understand the phylogenetic relationships among the Eurasian badgers (Meles, Mustelidae, Carnivora), which are distributed widely in the Palearctic, partial sequences of the mitochondrial DNA (mtDNA) control region (539-545 base-pairs) as a maternal genetic marker, and the sex-determining region on the Y-chromosome gene (SRY: 1052-1058 base-pairs), as a paternal genetic marker, were examined. The present study revealed ten SRY haplotypes from 47 males of 112 individuals of the Eurasian Continent and Japan. In addition, 39 mtDNA haplotypes were identified from those animals. From the phylogeography of both the uniparentally inherited genes, four lineages were recognized as Japanese, eastern Eurasian, Caucasian, and western Eurasian. The distribution patterns of the mtDNA lineages showed the existence of a sympatric zone between the eastern and western Eurasian lineages around the Volga River in western Russia. Furthermore, the present study suggested that in the Japanese badgers, the larger genetic differentiation of the Shikoku population was attributable to geographic history in the Japanese islands.     

Phylogenetic Relationships Among Fragmented Asian Black Bear (Ursus Thibetanus) Populations in Western Japan

The number of Asian black bears (Ursus thibetanus) in Japan has been reduced and their habitats fragmented and isolated because of human activities. Our previous study examining microsatellite DNA loci revealed significant genetic differentiation among four local populations in the western part of Honshu. Here, an approximate 700-bp nucleotide sequence of mitochondrial DNA (mtDNA) control region was analysed in 119 bears to infer the evolutionary history of these populations. Thirteen variable sites and variation in the number of Ts at a T-repeat site were observed among the analysed sequences, which defined 20 mtDNA haplotypes with the average sequence divergence of 0.0051 (SD = 0.00001). The observed haplotype frequencies differed significantly among the four populations. Phylogeographic analysis of the haplotypes suggested that black bears in this region have gone through two different colonisation histories, since the observed haplotypes belonged to two major monophyletic lineages and the lineages were distributed with an apparent border. The spatial genetic structure revealed by using mtDNA was different from that observed using microsatellite DNA markers, probably due to female philopatry and male-biased dispersal. Since nuclear genetic diversity will be lost in the three western populations because of the small population size and genetic isolation, their habitats need to be preserved, and these four populations should be linked to each other by corridors to promote gene flow from the easternmost population with higher nuclear genetic diversity.     

High genetic diversity in the blue-listed British Columbia population of the purple martin maintained by multiple sources of immigrants

To assess genetic diversity in the blue-listed purple martin (Progne subis) population in British Columbia, we analysed mitochondrial control region sequences of 93 individuals from British Columbia and 121 individuals collected from seven localities of the western and eastern North American subspecies P. s. arboricola and P. s. subis, respectively. Of the 47 haplotypes we detected, 34 were found exclusively in western populations, and 12 were found only in eastern populations. The most common eastern haplotype (25) was also found in three nestlings in British Columbia and one in Washington. Another British Columbia nestling had a haplotype (35) that differed by a C to T transition from haplotype 25. Coalescent analysis indicated that these five nestlings are probably descendents of recent immigrants dispersing from east to the west, because populations were estimated to have diverged about 200,000–400,000 ybp, making ancestral polymorphism a less likely explanation. Maximum likelihood estimates of gene flow among all populations detected asymmetrical gene flow into British Columbia not only of rare migrants from the eastern subspecies in Alberta but also a substantial number of migrants from the adjacent Washington population, and progressively lower numbers from Oregon in an isolation-by distance pattern. The influx of migrants from different populations is consistent with the migrant-pool model of recolonization which has maintained high genetic diversity in the small recovering population in British Columbia. Thus, the risk to this population is not from genetic erosion or inbreeding following a severe population crash, but from demographic stochasticity and extinction in small populations.     

Phylogeography and introgressive hybridization of the ground beetle Carabus yamato in Japan based on mitochondrial gene sequences

To study the phylogeography of the ground beetle Carabus yamato in Japan, we compared 1,020-bp sequences of the mitochondrial NADH dehydrogenase subunit 5 (ND5) gene from 373 specimens from 37 localities with those of three parapatric species (C. albrechti, C. kimurai, and C. japonicus) that might share mitochondrial lineages with C. yamato through introgressive hybridization. We found 81 haplotypes from C. yamato. Of these, 17 haplotypes were considered to be of an introgressed lineage from C. albrechti, based on the phylogeny and geographic distribution. In addition, one haplotype of C. kimurai was likely an introgressant from C. yamato. Putative introgression events among the four species were restricted to these two directional cases. We analyzed the phylogeography of C. yamato using nested clade phylogeographical analysis and population genetic parameters. The mitochondrial lineages of C. yamato were estimated to have diverged no more than approximately 1.12 million years ago, implying that the estimated historical events occurred after the Early Pleistocene. Carabus yamato was inferred to have experienced a contraction of its distribution range, followed by recent range expansion. Populations in the western and eastern regions, segregated by Ise Bay and the Nobi and Okazaki Plains, diverged in the mitochondrial clades. The northern and most western populations possessed one clade only (except an introgressed lineage), whereas eastern and some southwestern populations possessed several diverged clades, which were considered to be ancestral; these populations may have been associated with refugia during glacial periods.     

Late Pleistocene divergence between eastern and western populations of wood ducks (Aix sponsa) inferred by the 'isolation with migration' coalescent method

During the Late Pleistocene, glaciers sundered many species into multiple glacial refugia where populations diverged in allopatry. Although deeply divergent mitochondrial DNA (mtDNA) lineages often reflect the number of refugia occupied, it is unlikely that populations that split during the recent Wisconsin glaciations will have reached reciprocal monophyly. We examined mtDNA control region sequences from eastern and western populations of wood ducks (Aix sponsa) to determine whether their current, disjunct distribution is consistent with the occupancy of two glacial refugia. We used the 'isolation with migration' coalescent method (im) to simultaneously estimate effective population sizes, maternal gene flow, and time since divergence. We found 24 unique haplotypes, none of which were shared between the eastern and western populations, but we did not find diagnostic monophyletic lineages suggestive of long-term isolation in multiple glacial refugia. However, a high Phi ST (0.31) indicates that eastern and western populations are well differentiated in mtDNA, and results from im suggest that these populations have been diverging, without extensive gene flow, for 10,000 to 124,000 years. Results from im further suggest that these populations most likely split about 34,000 years ago, and this time of divergence is consistent with the occupancy of multiple glacial refugia during the Late Wisconsin glaciation. Eastern wood ducks are characterized by high genetic diversity, a large effective population size, and a recent population expansion, while western wood ducks have much less genetic diversity, a smaller population size, and have not undergone a recent population expansion.     

Phylogeography of the endangered Otago skink, Oligosoma otagense: population structure, hybridisation and genetic diversity in captive populations

Climatic cooling and substantial tectonic activity since the late Miocene have had a pronounced influence on the evolutionary history of the fauna of New Zealand's South Island. However, many species have recently experienced dramatic range reductions due to habitat fragmentation and the introduction of mammalian predators and competitors. These anthropogenic impacts have been particularly severe in the tussock grasslands of the Otago region. The Otago skink (Oligosoma otagense), endemic to the region, is one of the most critically endangered vertebrates in New Zealand. We use mitochondrial DNA sequence data to investigate the evolutionary history of the Otago skink, examine its population genetic structure, and assess the level of genetic diversity in the individuals in the captive breeding program. Our data indicate that the Otago skink diverged from its closest relatives in the Miocene, consistent with the commencement of tectonic uplift of the Southern Alps. However, there is evidence for past introgression with the scree skink (O. waimatense) in the northern Otago-southern Canterbury region. The remnant populations in eastern Otago and western Otago are estimated to have diverged in the mid-Pliocene, with no haplotypes shared between these two regions. This divergence accounts for 95% of the genetic diversity in the species. Within both regions there is strong genetic structure among populations, although shared haplotypes are generally evident between adjacent localities. Although substantial genetic diversity is present in the captive population, all individuals originate from the eastern region and the majority had haplotypes that were not evident in the intensively managed populations at Macraes Flat. Our data indicate that eastern and western populations should continue to be regarded as separate management units. Knowledge of the genetic diversity of the breeding stock will act to inform the captive management of the Otago skink and contribute to a key recovery action for the species.     

Phylogeography of the Coastal Mosquito Aedes togoi across Climatic Zones: Testing an Anthropogenic Dispersal Hypothesis

The coastal mosquito Aedes togoi occurs more or less continuously from subarctic to subtropic zones along the coasts of the Japanese islands and the East Asian mainland. It occurs also in tropical Southeast Asia and the North American Pacific coast, and the populations there are thought to have been introduced from Japan by ship. To test this hypothesis, the genetic divergence among geographic populations of A. togoi was studied using one mitochondrial and three nuclear gene sequences. We detected 71 mitochondrial haplotypes forming four lineages, with high nucleotide diversity around temperate Japan and declining towards peripheral ranges. The major lineage (L1) comprised 57 haplotypes from temperate and subarctic zones in Japan and Southeast Asia including southern China and Taiwan. Two other lineages were found from subtropical islands (L3) and a subarctic area (L4) of Japan. The Canadian population showed one unique haplotype (L2) diverged from the other lineages. In the combined nuclear gene tree, individuals with mitochondrial L4 haplotypes diverged from those with the other mitochondrial haplotypes L1—L3; although individuals with L1—L3 haplotypes showed shallow divergences in the nuclear gene sequences, individuals from Southeast Asia and Canada each formed a monophyletic group. Overall, the genetic composition of the Southeast Asian populations was closely related to that of temperate Japanese populations, suggesting recent gene flow between these regions. The Canadian population might have originated from anthropogenic introduction from somewhere in Asia, but the possibility that it could have spread across the Beringian land bridge cannot be ruled out.     

Response of a southern temperate marsupial, the Tasmanian pademelon (Thylogale billardierii), to historical and contemporary forest fragmentation

Despite only limited Pleistocene glacial activity in the southern hemisphere, temperate forest species experienced complex distributional changes resulting from the combined effects of glaciation, sea level change and increased aridity. The effects of these historical processes on population genetic structure are now overlain by the effects of contemporary habitat modification. In this study, 10 microsatellites and 629 bp of the mitochondrial control region were used to assess the effects of historical forest fragmentation and recent anthropogenic habitat change on the broad-scale population genetic structuring of a southern temperate marsupial, the Tasmanian pademelon. A total of 200 individuals were sampled from seven sites across Tasmania and two islands in Bass Strait. High mitochondrial and nuclear genetic diversity indicated the maintenance of large historical population sizes. There was weak phylogeographical structuring of haplotypes, although all King Island haplotypes and three Tasmanian haplotypes formed a divergent clade implying the mid-Pleistocene isolation of a far northwestern population. Both the mitochondrial and nuclear data indicated a division of Tasmanian populations into eastern and western regions. This was consistent with a historical barrier resulting from increased aridity in the lowland 'midlands' region during glacial periods, and with a contemporary barrier resulting from recent habitat modification in that region. In Tasmania, gene flow appears to have been relatively unrestricted during glacial maxima in the west, while in the east there was evidence for historical expansion from at least one large glacial refuge and recolonization of Flinders Island.     

Evidence of a phylogeographic break in the Romanian brown bear (Ursus arctos) population from the Carpathians

We analysed 33 brown bears from the Romanian Carpathians and the Italian Apennines at sequences of the mitochondrial control region and nine polymorphic microsatellite loci with regard to genetic variability and haplotype distribution. The Italian brown bears were monomorphic for mtDNA sequences. The Romanian bears yielded the highest variability found so far in this species. Haplotypes of both previously identified mtDNA lineages (western and eastern) were found in Romania. In the eastern part of the Carpathians western and eastern haplotypes occurred sympatrically, the bears from the western part of the mountain range only exhibited western-type sequences. This pattern provides evidence of a mitochondrial phylogeographic break in the distribution of the eastern lineage within the Romanian Carpathians. Conservation implications of this finding are discussed.     

Postglacial population expansion of Japanese macaques (<Emphasis Type="Italic">Macaca fuscata</Emphasis>) inferred from mitochondrial DNA phylogeography

We investigated the diversity and phylogeography of mitochondrial DNA (mtDNA) in Japanese macaques (Macaca fuscata), an endemic species in Japan that has the northernmost distribution of any non-human primate species. DNA samples from 135 localities representing the entire range of this species were compared. A total of 53 unique haplotypes were observed for the 412-bp partial mtDNA control region sequence, with length variation distinguishing the two subspecies. Clustering analyses suggested two putative major haplogroups, of which one was geographically distributed in eastern Japan and the other in western Japan. The populations in the east showed lower mtDNA diversity than those in the west. Phylogeographical relationships of haplotypes depicted with minimum spanning network suggested differences in population structure. Population expansion was significant for the eastern but not the western population, suggesting establishment of the ancestral population was relatively long ago in the west and recent in the east. Based on fossil evidence and past climate and vegetation changes, we inferred that the postulated population expansion may have taken place after the last glacial period (after 15,000 years ago). Mitochondrial DNA showed contrasting results in both variability and phylogenetic status of local populations to those of previous studies using protein variations, particularly for populations in the periphery of the range, with special inference on habitat change during the glacial period in response to cold adaptation. Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.     

Genetic status of Asiatic black bear (Ursus thibetanus) reintroduced into South Korea based on mitochondrial DNA and microsatellite loci analysis

The Asiatic black bear is one of the most endangered mammals in South Korea owing to population declines resulting from human exploitation and habitat fragmentation. To restore the black bear population in South Korea, 27 bear cubs from North Korea and Russian Far East (Primorsky Krai) were imported and released into Jirisan National Park, a reservoir of the largest wild population in South Korea, in 2004. To monitor the success of this reintroduction, the genetic diversity and population structure of the reintroduced black bears were measured using both mitochondrial and nuclear DNA markers. Mitochondrial D-loop region DNA sequences (615 bp) of 43 Japanese black bears from previous study and 14 Southeast Asian black bears in this study were employed to obtain phylogenetic inference of the reintroduced black bears. The mitochondrial phylogeny indicated Asiatic black bear populations from Russian Far East and North Korea form a single evolutionary unit distinct from populations from Japan and Southeast Asia. Mean expected heterozygosity (H(E)) across 16 microsatellite loci was 0.648 for Russian and 0.676 for North Korean populations. There was a moderate but significant level of microsatellite differentiation (F(ST) = 0.063) between black bears from the 2 source areas. In addition, genetic evidences revealed that 2 populations are represented as diverging groups, with lingering genetic admixture among individuals of 2 source populations. Relatedness analysis based on genetic markers indicated several discrepancies with the pedigree records. Implication of the phylogenetic and genetic evidences on long-term management of Asiatic black bears in South Korea is discussed.     

Genetic variation and population structure of the carpet shell clam <Emphasis Type="Italic">Ruditapes decussatus</Emphasis> along the Tunisian coast inferred from mtDNA and ITS1 sequence analysis

Surveys of allozyme polymorphisms in the carpet shell clam Ruditapes decussatus have revealed sharp genetic differentiation of populations. Analysis of population structure in this species has now been extended to include nuclear and mitochondrial genes. A partial sequence of a mitochondrial COI gene and of the internal transcribed spacer region (ITS-1) were used to study haplotype distribution, the pattern of gene flow, and population genetic structure of R. decussatus. The samples were collected from twelve populations from the eastern and western Mediterranean coasts of Tunisia, one from Concarneau and one from Thau. A total of twenty and twenty-one haplotypes were detected in the examined COI and ITS1 regions respectively. The study revealed higher levels of genetic diversity for ITS1 compared to COI. The analysis of haplotype frequency distribution and molecular variation indicated that the majority of the genetic variation was distributed within populations (93% and 86% for COI and ITS1 respectively). No significant differentiation was found among eastern and western groups on either side of the Siculo-Tunisian strait. However, distinct and significant clinal changes in haplotypes frequencies between eastern and western samples were found at the most frequent COI haplotype and at three out of five major ITS1 haplotypes. These results suggest the relative importance of historical processes and contemporary hydrodynamic features on the observed patterns of genetic structure.     

Radiation and phylogeography in the Japanese macaque, Macaca fuscata

The Japanese macaque (Macaca fuscata) presumably differentiated from eastern rhesus macaque (Macaca mulatta) populations during the Pleistocene and the two species are closely related. In order to analyse speciation and subspeciation events in the Japanese macaque and to describe historical and current relationships among their populations, we sequenced and analysed a fragment of 392bp of mitochondrial DNA (mtDNA) control region in 50 individuals belonging to six populations of Japanese macaque and compared these sequences with 89 eastern rhesus macaque control region sequences from GenBank/EMBL database. There were high genetic similarities between both species and only two positions were fixed within each species, which supports the inclusion of the Japanese macaque in a single species with eastern populations of rhesus macaques. Japanese macaque ancestors colonised Japan after the separation of the two species, estimated at between 0.31 and 0.88 million years ago (Mya). The star-like phylogeny, multimodal mismatch distribution, and lack of correlation between geographic and genetic distances are in accordance with a rapid dispersion of macaques throughout the archipelago after the arrival into Japan. The species shows low genetic variation within populations and high levels of genetic differentiation among populations with no mtDNA haplotype shared across populations. Genetic distances between Yakushima macaques (Macaca fuscata yakui) and any other population of Macaca fuscata fuscata subspecies are comparable to the distances between populations of Honshu, Awajishima, and Kyushu, not supporting the classification of Yakushima macaques as a different subspecies.     

Population differentiation in the Pacific white-sided dolphin Lagenorhynchus obliquidens inferred from mitochondrial DNA and microsatellite analyses

We investigated genetic diversity and differentiation of the Pacific white-sided dolphin (Lagenorhynchus obliquidens) in Japanese coastal waters and offshore North Pacific by analyzing mitochondrial DNA and nuclear microsatellite variation. A total of 519 bp of the mitochondrial control region was sequenced and five microsatellite locus were genotyped for 59 individuals. A high level of haplotypic diversity (h=96.1%), moderate level of nucleotide diversity (pi=1.65%) and average expected heterozygosity (HE=0.66-0.76) were within an extent of those reported for other odontocetes. Consistent genetic difference between the samples from Japanese coastal Pacific-Sea of Japan and offshore North Pacific was indicated by analyses of molecular variance (AMOVAs) based on mtDNA and microsatellite variations, comparison of genetic variabilities, and geographical distributions of mtDNA haplotypes and microsatellite alleles. This result suggests that Pacific white-sided dolphins in each of the above two areas belong to different populations between which gene flow has been severely restricted. The low genetic diversity and mtDNA genealogy of the population in Japanese coastal waters suggest that it originated from a small population that colonized the Sea of Japan or that experienced population reduction when this Sea was isolated from the North Pacific during a glacial period in the Late Pleistocene.     

Phylogeography and genetic structure of northern populations of the yellow warbler (Dendroica petechia)

Phylogeographic patterns of intraspecific variation can provide insights into the population-level processes responsible for speciation and yield information useful for conservation purposes. To examine phylogeography and population structure in a migratory passerine bird at both continental and regional geographical scales, we analysed 344 bp of mitochondrial DNA (mtDNA) control region sequence from 155 yellow warblers (Dendroica petechia) collected from seven locations across Canada and from Alaska. There is a major subdivision between eastern (Manitoba to Newfoundland) and western (Alaska and British Columbia) populations which appears to have developed during the recent Pleistocene. Some localities within these two regions also differ significantly in their genetic composition, suggesting further subdivision on a regional geographical scale. Eastern and western birds form distinct phylogeographic entities and the clustering of all western haplotypes with two eastern haplotypes suggests that the western haplotypes may be derived from an eastern lineage. Analyses based on coalescent models support this explanation for the origin of western haplotypes. These results are consistent with important features of Mengel's model of warbler diversification. From a conservation perspective they also suggest that individual populations of migrant birds may form demographically isolated management units on a smaller scale than previously appreciated.     

Phylogeography of the white-tailed eagle, a generalist with large dispersal capacity

Aim Late Pleistocene glacial changes had a major impact on many boreal and temperate taxa, and this impact can still be detected in the present‐day phylogeographic structure of these taxa. However, only minor effects are expected in species with generalist habitat requirements and high dispersal capability. One such species is the white‐tailed eagle, Haliaeetus albicilla, and we therefore tested for the expected weak population structure at a continental level in this species. This also allowed us to describe phylogeographic patterns, and to deduce Ice Age refugia and patterns of postglacial recolonization of Eurasia. Location Breeding populations from the easternmost Nearctic (Greenland) and across the Palaearctic (Iceland, continental Europe, central and eastern Asia, and Japan). Methods Sequencing of a 500 base‐pair fragment of the mitochondrial DNA control region in 237 samples from throughout the distribution range. Results Our analysis revealed pronounced phylogeographic structure. Overall, low genetic variability was observed across the entire range. Haplotypes clustered in two distinct haplogroups with a predominantly eastern or western distribution, and extensive overlap in Europe. These two major lineages diverged during the late Pleistocene. The eastern haplogroup showed a pattern of rapid population expansion and colonization of Eurasia around the end of the Pleistocene. The western haplogroup had lower diversity and was absent from the populations in eastern Asia. These results suggest survival during the last glaciation in two refugia, probably located in central and western Eurasia, followed by postglacial population expansion and admixture. Relatively high genetic diversity was observed in northern regions that were ice‐covered during the last glacial maximum. This, and phylogenetic relationships between haplotypes encountered in the north, indicates substantial population expansion at high latitudes. Areas of glacial meltwater runoff and proglacial lakes could have provided suitable habitats for such population growth. Main conclusions This study shows that glacial climate fluctuations had a substantial impact on white‐tailed eagles, both in terms of distribution and demography. These results suggest that even species with large dispersal capabilities and relatively broad habitat requirements were strongly affected by the Pleistocene climatic shifts.     

Regionally high rates of hybridization and introgression in German wildcat populations (Felis silvestris, Carnivora, Felidae)

While the western populations of the wildcat ( Felis silvestris silvestris ) in Germany come into contact with wildcats in France and Switzerland, the eastern distribution area is geographically completely isolated and consists of scattered subpopulations. To investigate population structure, evolutionary relationships and degree of hybridization with domestic cats we analysed the mitochondrial control region of 86 cats in combination with 11 microsatellite loci of 149 cats. According to our microsatellite data, German wildcats are divided into two separate populations corresponding to the western and eastern distribution areas. We found no indication of a further subdivision of the eastern population. German wildcat populations are genetically distinct from domestic cats in the main, but we identified 18.4% of the whole wildcat sample as being of hybrid origin, corresponding to 4.2% of the eastern and 42.9% of the western wildcat population, and 2.7% of the domestic cat sample. The mitochondrial haplotypes form a network of three connected clusters and reveal a high level of genetic diversity, especially within the eastern population. Our findings are explained at best in terms of continuous introgression between domestic cats and wildcat populations and differing degrees of recent hybridization in the various populations. Future conservation efforts should focus on preserving the existing gene flow between the isolated distribution areas, but also on preventing the spread of hybrids and limiting the habitat alterations that lead to increased contact with domestic cats. In conclusion we discuss possible evolutionary reasons for the still traceable genetic integrity of the wildcat despite its long history of interbreeding.     

Effects of topography on dispersal of black-billed magpie Pica pica sericea revealed by population genetic analysis

This study was conducted using random amplified polymorphic DNA analysis to examine whether the effects of topographical structure on the dispersal of black-billed magpie Pica pica sericea are reflected in its population genetic structure. The black-billed magpie is typically seen in the lowlands and is remarkably sedentary in Korea and Japan. The unweighted pair-group method of clustering analysis showed two main clusters: five populations in the western region of the Bekdudegan mountains (WRBM) in one cluster and five populations in the eastern region of the Bekdudegan mountains (ERBM), including the Japanese population in the other cluster. The populations in WRBM appeared to be diverged from the populations in ERBM. Analysis of molecular variance revealed that the populations in ERBM were more genetically divergent from each other than were those in WRBM. The high-rising mountains, very rugged topographical features and the sea in ERBM may have resulted in a lower dispersal rate and larger genetic variation among populations in ERBM compared to those in WRBM. These results suggest that the topographical structure may have an influence on the dispersal and population genetic structure of the black-billed magpie. Electronic Publication     

Phylogeography of the Japanese greater horseshoe bat Rhinolophus nippon (Mammalia: Chiroptera) in Northeast Asia: New insight into the monophyly of the Japanese populations

The Japanese greater horseshoe bat (Rhinolophus nippon) is distributed widely in East Asia. Within the species, R. nippon in Northeast Asia is regarded as the lineage that diverged most recently. However, the monophyly of the Japanese populations is unclear due to insufficient data about phylogenetic relationship of the western Japanese populations. To test the monophyly of the Japanese populations of R. nippon, we sampled R. nippon from western Japan and performed a phylogeographic analysis based on mitochondrial DNA cytochrome b and the D‐loop. The Northeast Asian lineage consisted of three main clades in eastern Japan (clade I), western Japan (clade II), and the continent as well as the Kumamoto population in westernmost Japan (clade III). The results of this study do not support the monophyly of the Japanese population. The findings suggest the “reverse colonization” of R. nippon from the Japanese Archipelago to the Eurasian continent, and provide important insight into the role of the island system in creation and supply of diversity to the continent.     

Mitochondrial DNA Sequence-Based Population Genetic Structure of the Firefly, Pyrocoelia rufa (Coleoptera: Lampyridae)

The genetic divergence, population genetic structure, and possible speciation of the Korean firefly, Pyrocoelia rufa, were investigated on the midsouthern Korean mainland, coastal islets, a remote offshore island, Jedu-do, and Tsushima Island in Japan. Analysis of DNA sequences from the mitochondrial COI protein-coding gene revealed 20 mtDNA-sequence-based haplotypes with a maximum divergence of 5.5%. Phylogenetic analyses using PAUP, PHYLIP, and networks subdivided the P. rufa into two clades (termed clade A and B) and the minimum nucleotide divergence between them was 3.7%. Clade A occurred throughout the Korean mainland and the coastal islets and Tsushima Island in Japan, whereas clade B was exclusively found on Jeju-do Island. In the analysis of the population genetic structure, clade B formed an independent phylogeographic group, but clade A was further subdivided into three groups: two covering western and eastern parts of the Korean peninsula, respectively, and the other occupying one eastern coastal islet and Japanese Tsushima Island. Considering both phylogeny and population structure of P. rufa, the Jeju-do Island population is obviously differentiated from other P. rufa populations, but the Tsushima Island population was a subset of the Korean coastal islet, Geoje. We interpreted the isolation of the Jeju-do population and the grouping of Tsushima Island with Korean coastal islets in terms of Late Pleistocene-Holocene events. The eastern-western subdivision on the Korean mainland was interpreted partially by the presence of a large major mountain range, which bisects the midpart of the Korean peninsula into western and eastern parts.