Mitochondrial phylogeography of the Japanese pond turtle,Mauremys japonica (Testudines,Geoemydidae)

The phylogeography of the Japanese pond turtle, Mauremys japonica (Temminck and Schlegel, 1835), which is an endemic species in Japan, was studied by analyzing the variation in two mitochondrial DNA sequences, the cytochrome b gene and the control region. These analyses suggest that M. japonica comprises two major groups. The first one was found in the more eastern region, eastern Honshu Island and Shikoku Island, while the second was found in a western region, Kyushu Island and the Chugoku District (the westernmost part of Honshu Island). The boundary between the two groups is located in the Chugoku District. The nucleotide and haplotype diversities were very low, and these low diversities seem to have been caused by a bottleneck in the last glacial age. These results suggest that this species survived the last glacial period in two refugia, one in the central part of Honshu Island and the other one in Kyushu Island. Subsequently, population expansion took place in the postglacial period, and the groups from the two refugia extended their distribution ranges to the present boundary in the Chugoku District which represents a secondary contact zone.     

Geographic allozyme variation in the Japanese clouded salamander, Hynobius nebulosus (Amphibia: Urodela)

Using starch gel electrophoresis, geographical enzyme variation was surveyed in 46 populations of the Japanese salamander Hynobius nebulosus , which occurs widely in western Japan. This species exhibits substantial local genetic differentiation and is more diverse in inland regions than on small islands. In several different analyses, two groups of populations, one from Kyushu Island to the westernmost part of Chugoku district, Honshu (western group), and another from Shikoku through Kinki to Chubu district, Honshu (eastern group), were consistently recognized. Most of the remaining populations were placed in two less clearly-defined groups: the montane group from the Chugoku Mountains and the Chugoku group from the coastal regions of Chugoku district. The divergence patterns of H. nebulosus are thought to be related to the geological history of its range. The populations near the boundary of the western group were found to share several alleles that were predominant on both sides of the boundary, indicating past secondary contact. Differentiation in the montane group did not follow the isolation by distance model, probably because several populations included were genetically mixed with lowland populations of the Chugoku group.  © 2006 The Linnean Society of London, Biological Journal of the Linnean Society , 2006, 89 , 311–330.     

Morphological discrimination of two genetic groups of a Japanese salamander, Hynobius naevius (Amphibia, Caudata)

Hynobius naevius, distributed on western Honshu, Shikoku, and Kyushu Islands of Japan, includes two genetically distinct groups (Groups A and B) that have never been delimited morphologically. Using specimens from the entire species range, we investigated the possibility of distinguishing these groups morphologically. Multivariate analyses of morphometric characters resulted in recognition of two groups that corresponded well to the two genetic groups. One (Group A) was characterized by larger body, compressed tail, shallower vomerine tooth series, bluish- or reddish-purple ground color, and pale-white lateral markings. In contrast, another (Group B) was characterized by smaller body, cylindrical tail, longer vomerine tooth series, reddish-brown ground color, and white lateral markings. Group A was composed of populations from the Chugoku District of Honshu and northern Kyushu, and could not be divided into subgroups, while Group B encompassed populations from the Chubu and Kinki Districts of Honshu, Shikoku and Kyushu, and was subdivided into three local subgroups that are geographically separated by marine straits. Morphometric differentiation in Group A is presumed to have been less affected by genetic factors than by other factors, such as ecological relationships with other, coexisting species. Differentiation in Group B is assumed to have been enhanced not only by genetic but also by climatological factors.     

Chloroplast DNA phylogeography of Photinia glabra (Rosaceae) in Japan

Climate changes during glacial periods have had significant effects on the present geographic distribution of plant species. To elucidate the evolutionary history of a plant species with a disjunctive distribution, we investigated the geographic distribution patterns of cpDNA haplotypes in Photinia glabra (Rosaceae) growing in broadleaved evergreen forests in Japan. We examined cpDNA in 42 populations of P. glabra sampled over a geographic range that included Kinki and its surrounding areas and the disjunctive regions in the Amakusa Islands. Both areas had unique cpDNA haplotypes. Moreover, the AMOVA revealed that a large proportion of the total variance (51%, P < 0.001) could be explained by differences among these regions. These results suggest a past fragmentation of this plant species into two separate refugia: southwestern Kyushu and Kinki, including the surrounding area, during the Quaternary glacial periods. A particularly interesting result was that in the southern disjunct distribution in the Amakusa Islands, the genetic subdivision (Φ(CT) = 1.00, P < 0.001) appears to lie between the populations from nearly contiguous islands located across a fairway only approximately 80 to 150 m in width.     

Parsimony analysis of the distribution pattern of Japanese primary freshwater fishes, and its application to the distribution of the bagrid catfishes

The geographic distribution of Japanese primary and some secondary freshwater fishes was analyzed using parsimony analysis of endemicity (PAE). Analysis of 73 taxa (species and intraspecific forms) on the four main islands of Japan (divided into 25 geographic areas), produced 44 most parsimonious cladograms. In all of the latter, a total of eight single and compound areas were recognized as endemic areas in nested relationships. The area cladograms showed Japan as comprising a middle-eastern Hokkaido area plus southern areas, the latter containing mainly a northeastern-Honshu endemic area and more heterogeneous southwestern areas, including four endemic areas (western Kyushu, southeastern Chugoku, middle Kinki and Tokai around Ise Bay) and several peripheral areas. Some patterns, e.g., the distinction in fauna across Fossa Magna, noted by previous studies, were supported by these results. Even though the analysis had some problems (e.g., not all geographic ranges of taxa could be included), it provided evidence for the detection of general distribution patterns, because the relationships or similarities among areas were clearly defined by shared taxa. To demonstrate the historical implications of the analysis, the allopatric distribution of, four bagrid catfishes was reconsidered in the area cladogram. The general pattern implied secondary extinction ofPseudobagrus nudiceps around lse Bay, which was in keeping with the fossil record.     

Cytogeography and meiotic chromosome configurations of six intraspecific aneuploids ofCarex conica Boott (Cyperaceae) in Japan

Chromosome numbers were determined for 340 plants ofCarex conica from 83 populations in Japan. Six aneuploids, 2n=32, 33, 34, 36, 37 and 38, were found. Plants with even diploid chromosome numbers 2n=32, 34, and 36 were the most common and had different geographical distributions. Individuals with 2n=32 were from islands in the Seto Inland Sea and nearby coastal areas of the Chugoku District of Honshu; those with 2n=34 were from the Kanto, Chubu and Kinki Districts of Honshu; those with 2n=36 were from the mountainous areas of Chugoku, Shikoku and Kyushu Districts. Canonical discriminant analysis of 17 morphological characters demonstrated that the plants with 2n=32 were clearly distinct from those with 2n=34 or 36. All four aneuploids with even chromosome numbers showed normal bivalent pairing at meiotic metaphase I and probably represent cytogenetically stable cytodemes. Plants with 2n=33 had one heteromorphic trivalent and 15 bivalents, indicating a structural mutation. At mitotic metaphase I, one chromosome was markedly larger than the others, suggesting that the 2n=33 plants arose from 2n=34 plants by fusion of two chromsomes. The plant with 2n=37 was intermediate in morphology betweenCarex conica (2n=36) andC. morrowii (2n=38) and probably originated as an interspecific hybrid between these species.     

Phylogeography of Japanese population of Phelotrupes auratus (Coleoptera, Geotrupidae) inferred from mitochondrial DNA sequences

A partial sequence of the mitochondrial cytochrome c oxidase subunit I (COI) gene (745 bp) was determined for 57 specimens of a geotrupid beetle (Phelotrupes auratus) from throughout the Japanese archipelago. Of the 57 beetles examined, 42 haplotypes were identified. Phylogenetic trees inferred using maximum parsimony, neighbor joining, and Bayesian inference methods were highly congruent. Reconstructed phylogenetic relationships indicated that P. auratus from the Japanese archipelago was separated into two distinct lineages: Group A, which consisted of 35 haplotypes from Honshu, Shikoku, and Hokkaido Islands, and Group B, which consisted of seven haplotypes from Kyushu and Yakushima Islands. In addition, two sublineages were also recognized within Group A: Subgroup A-1, which consisted of 11 haplotypes from eastern Honshu and Hokkaido, and Subgroup A-2, which consisted of 10 haplotypes from western Honshu and Shikoku. Average genetic distances within Group A were positively correlated to geographic distance between sampling localities. Phylogenetic relationships among haplotypes did not correspond to subspecies classification.     

Refugia of Potentilla matsumurae (Rosaceae) located at high mountains in the Japanese archipelago

Molecular phylogeographic studies have revealed the genetic patterns and glacial-interglacial history of many plant and animal species. To infer the Quaternary history of alpine plants in the Japanese archipelago, which is poorly known, we investigated 203 individuals of Potentilla matsumurae and its varieties from 22 populations. We found 11 haplotypes based on approximately 1400 bp of two intergenetic spacers in chloroplast DNA (trnT-L and rpl20-rps20). The distribution of these haplotypes was geographically structured, which was supported by haplotype composition, principal component analysis, and unweighted pair group method with arithmetic mean (UPGMA), and N(ST) (0.71) was significantly greater than G(ST) (0.68). In addition to the positive correlation between genetic and geographic distance (Mantel test, r = 0.497, P < 0.001), an abrupt genetic change was detected between mountains in central Honshu and the Tohoku region. This genetic boundary was further supported by analysis of molecular variance (AMOVA), and high variation (54.0%) was explained by differences on either side of this boundary. Moreover, haplotypes in central Honshu were thought to have diverged, based on an outgroup comparison. These results suggest that mountains in central Honshu served as refugia during the Quaternary climatic oscillation, although the results could not reveal the history of most mountains in the Tohoku region and Hokkaido. Nevertheless, following floristic studies, our results indicate that alpine plants in Japan experienced a history different from that in Europe; i.e. they retreated into refugia during warm periods to avoid forest development, rather than glaciers.     

Spatial framework of nine distinct local populations of the Japanese dormouse Glirulus japonicus based on matrilineal cytochrome b and patrilineal SRY gene sequences

We previously revealed the presence of six genetically distinct matrilineal populations of the Japanese dormouse Glirulus japonicus in the distribution range of Honshu, Shikoku, and Kyushu islands. In this study, we extended this analysis using mitochondrial cytochrome b gene sequences (n = 96) and Y-chromosome-specific SRY gene sequences (n = 22) from individuals collected from Honshu, Shikoku, Kyushu, and Oki Dogo I. The cytochrome b sequence data allowed us to define precise geographic ranges of the six previously known and three newly found distinct matrilineal lineages: northeastern Honshu (I), east-central Honshu (II), west-central Honshu and the Kii Peninsula (III), the western part of Honshu (IV), Shikoku (V), westernmost Honshu and Kyushu (VI), the northern part of central Honshu (VII), the southern part of central Honshu (VIII), and Oki Dogo I. (IX). Our inference of geographic borders suggests that regions of lower and higher altitudes in the mountain systems played important roles in driving the hosting and separation of lineages, respectively. Six matrilineal lineages (I, II, V, VI, VIII, and XI) were shown to possess their own SRY haplotypes, while lineages III and IV shared one haplotype. These data together with our previous observation of nuclear ribosomal RNA gene variation indicate advanced populational subdivision in this species. It is thus evident that each of the populations, including those living at high latitudes and in limited geographic spaces, have survived for several million years. A specific ability to tolerate cold may have permitted G. japonicus to preserve anciently diverged lineages in each locality.     

Phylogeography of Arcterica nana (Ericaceae) suggests another range expansion history of Japanese alpine plants

We conducted a phylogeographic study on the alpine plant Arcterica nana based on haplotypes of chloroplast DNA. Using a sequence of approximately 1,071 bp of intergenic spacers of chloroplast DNA (trnT-L, psbB–psbF), we detected 13 haplotypes among 193 individuals sampled from 22 populations. Two dominant haplotypes were distributed over the entire range of this species in Japan, and we found several local private haplotypes. An analysis of molecular variance (AMOVA) indicated no geographic structure within the haplotype distribution. In addition, the genetic distance was not related to its corresponding geographic distance (Mantel test: r=−0.049, P=0.66), indicating a homogeneous geographic structure throughout the entire distribution range in the Japanese archipelago. The most parsimonious explanation for this geographic structure is that A. nana spread across its extant distribution range in the Japanese archipelago through a recent range expansion event. However, this pattern is inconsistent with the previous phylogeography of Japanese alpine plants, which reveals that haplotypes in central Honshu are differentiated from those in more northern regions. Arcterica nana may have experienced a different history from other alpine plants, suggesting that the history of Japanese alpine flora may include at least two different geographic radiation patterns.     

Multiple origins of Locusta migratoria (Orthoptera: Acrididae) in the Japanese Archipelago and the presence of two major clades in the world: evidence from a molecular approach

Phylogenetic relationships among migratory locust (Locusta migratoria) populations in different climatic regions were analysed by sequencing four mitochondrial DNA regions, with special reference to the origin of Japanese populations. The populations are clearly separated into two clades: one consists of individuals from temperate and cold‐temperate areas of Japan and the Chinese continent, and the other comprises those from subtropical islands of Japan, Hainan Island in China, Timol Leste, Australia, Ethiopia, France, and some individuals from Tsushima Island and Honshu of Japan. The divergence time between the two clades is estimated to be 0.86–1.89 Mya. The phylogenetic analysis revealed that Japanese L. migratoria populations were composed of individuals of six different origins: (1) Hokkaido populations possibly from the Russian continent; (2) Honshu–Kyushu populations from the Chinese continent; (3) Southwest Island populations from Hainan Island or adjacent areas; (4) Ogasawara populations that might have originated from Micronesia; (5) part of the Tsushima population that originated from somewhere in the Asian tropics; and (6) a possible relict population of ancient southern haplotypes that exists in western areas of northern Honshu. The Tsugaru Straits and Tokara Straits have acted as effective geographical barriers, as in other organisms, isolating locust populations for a few thousand years. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 99 , 570–581.     

Phylogenetic relationships and biogeography of the Japanese clawed salamander, Onychodactylus japonicus (Amphibia: Caudata: Hynobiidae), and its congener inferred from the mitochondrial cytochrome b gene

Using the mitochondrial cytochrome b gene, we investigated phylogenetic relationships between and within the Japanese clawed salamander, Onychodactylus japonicus, and its close continental relative O. fischeri. Monophyly of O. japonicus was well supported, and O. japonicus was clearly distinguished from O. fischeri. However, O. fischeri comprises genetically distinct populations from Russia, NE China, and Korea that do not form a monophyletic group. Within O. japonicus, four major clades were clearly recognized: Clade I from northern Tohoku district, Clade II from southern Tohoku district and the Tsukuba Mountains, Clade III from southwestern Honshu, and Clade IV from Kinki and Chugoku districts in Honshu and from Shikoku. Although genetic distances among these clades were large (5.5-9.6%), relationships among the clades were unresolved. All clades except Clade I contained two or three distinct subclades. In several localities in Kinki and Chugoku, Clades III and IV were sympatric. The estimated divergence times and available geohistorical data suggest that O. japonicus began to differentiate in the Upper Late Miocene and that the pattern of genetic differentiation of this species has been affected strongly by climate changes and geohistorical events such as volcanic activity and mountain formation. Our results suggest that both O. fischeri and O. japonicus comprise multiple cryptic species.     

Phylogeography of the leaf beetle Chrysolina virgata in wetlands of Japan inferred from the distribution of mitochondrial haplotypes

The genetic differentiation among populations of the leaf beetle Chrysolina virgata living in wetlands of Japan was studied based on the sequence data of the mitochondrial cytochrome oxidase subunit I gene region (750 bp). Two distinct lineages of mitochondrial haplotypes were found: one (clade A) consisted of 26 haplotypes distributed over the distribution range of C. virgata between north‐east Honshu and Kyushu, whereas the other (clade B) was monotypic and confined to a small region in north‐east Honshu where it coexisted with clade A. Nested clade analysis for these haplotypes suggested that range expansion and following differentiation due to isolation by distance might have resulted in the present distribution pattern of the haplotypes in clade A. We discuss the evolutionary process leading to the occurrence of two distinct haplotype clades in Japan in terms of repeated colonization from the continent and range expansion and contraction during climatic changes.     

Onset of cryptic vicariance in the Japanese dormouse Glirulus japonicus (Mammalia, Rodentia) in the Late Tertiary, inferred from mitochondrial and nuclear DNA analysis

The sequences of the mitochondrial cytochrome b gene and restriction site variation in the spacer region of the nuclear ribosomal RNA gene [rDNA-restriction fragment length polymorphism (RFLP)] were analysed to determine the phylogeographic structure of the Japanese dormouse ( Glirulus japonicus ), which is threatened by deforestation and has been designated an endangered species in Japan. The phylogenetic tree of cytochrome b grouped G. japonicus into six geographical populations: north-eastern Honshu (I), central Honshu (II), west-central Honshu/Kii Peninsula (III), western Honshu (IV), Shikoku (V), and westernmost Honshu/Kyushu (VI); the genetic distances among these groups suggest divergence in the Late Tertiary. The lineage of group VI was located at the basal position in the phylogenetic tree, followed by the radiation of the other lineages. An rDNA-RFLP analysis of 15 restriction sites roughly supported such genetic isolation; groups I, II, III, IV, V and VI have five, two, one, one, one and four unique restriction sites, respectively, revealing four geographic groups as cryptic species: I, II, III + IV + V and VI. Our results reveal the ancient divergences of the local population, which has a complicated evolutionary history, and should be useful in developing a framework for the conservation of this species.     

Genetic variation of the white-spotted longicorn beetle <Emphasis Type="Italic">Anoplophora</Emphasis> spp. (Coleoptera: Cerambycidae) in Japan detected by mitochondrial DNA sequence

Genetic variation of the white-spotted longicorn beetle Anoplophora spp., distributed on the Japanese mainland and the Ryukyu Islands, was examined using nucleotide sequences of mitochondrial DNA. Two fragments, a 1.2-kb-long fragment containing portions of the cytochrome oxidase subunit I and II genes and a 1.4-kb-long fragment containing portions of 16S and 12S rDNAs, were sequenced. In phylogenetic analyses based on the sequences, 294 individuals collected at 75 localities were divided into two major groups (groups A and B) and then split into seven subgroups (A1–A4 and B1–B3). In group A, the closely related subgroups A1, A2, and A3 were mainly distributed on the Japanese mainland and were roughly separated among geographic areas, although the range of A3 spread significantly from Kyushu to the Central Ryukyu Islands. The sequence of A4, detected from one individual collected in eastern Honshu, was almost the same as that reported for A. chinensis (Forster) on the Chinese continent. In group B, subgroups B1 and B2 were restricted to the Central and Southern Ryukyu Islands, respectively, while B3 was distributed widely in both regions. Based on the results, we discuss the geographic origin of the haplotypes and movement of the insects between geographic areas.     

Comments on cephalometrical configurations of the Japanese

Cephalometrical configurations of Japanese people were investigated on two hundred and sixteen groups living on the four main islands of Japan: Hokkaido, Honshu, Shikoku, and Kyushu. The Japanese in Honshu were divided geographically into a narrow-featured eastern major group and a broad-featured western major group. A physical characteristic of Japanese living in the western mid-Honshu was the lengthened feature of their cephalic measurement, an influence of a long-headed group of Korean people; Japanese residents in eastern mid-Honshu, however, showed the result of mixing with a short-headed group of Korean people. Japanese people living in Hokkaido, the northernmost island, and in the Ryukyu Islands in the South of Japan, presented no important finding which suggests their influence on the physical characteristics of the Japanese. Ongoing studies indicating physical relationships between Kyushu or Shikoku residents and Honshu residents, are as yet insufficient.     

Differential geographic patterns of mitochondrial DNA variation in two sympatric species of Japanese wood mice, Apodemus speciosus and A. argenteus

We examined the gene sequences of mitochondrial cytochrome b (cyt b) in two Japanese wood mouse species, Apodemus speciosus (n = 89) and A. argenteus (n = 46), which are distributed on the four main islands of Japan (Hokkaido, Honshu, Shikoku, and Kyushu) and on the small islands surrounding them. Apodemus speciosus, the larger of the two species, showed substantial genetic variation, with a maximum of 3% sequence divergence, and remarkable phylogenetic subdivision with two major clades. One clade represents haplotypes from a central region, including Honshu, Shikoku, Kyushu, and their adjacent islands; the other clade includes haplotypes from Hokkaido and the peripheral islands, forming four subclades: a) Hokkaido, b) Sado Island, c) Satsunan Islands, and d) the Izu Islands. Sequence divergence among the four subclades was 1.0 to 1.5%, implying that A. speciosus colonized these geographic regions 0.2 to 0.3 million years ago, assuming a substitution rate of 2.4% per million years. The population on the Izu Islands has preserved haplotypes that are distinct from those in any other region, providing good evidence for the natural colonization of the volcanic islands of the Izu Islands. The cyt b sequence variation had no relation to the karyotypic dimorphism for the eastern (2n = 48) and western (2n = 46) geographic groups, between which a strict border exists at central Honshu. On the other hand, Apodemus argenteus, the smaller of the two species, showed a similar level of sequence divergence (maximum of 3%) but no substantial geographic differentiation: populations in Hokkaido, Sado, and Yakushima shared similar haplotypes with each of the central populations, suggesting that genetic exchanges occurred among the localities in the last 0.15 million years. The apparent genetic structure of the mitochondrial DNA found in the A. speciosus population might be caused solely by long-term existence in insular regions, presumably due to ecological superiority relative to A. argenteus.     

Phylogeography of the endangered tideland snail Batillaria zonalis in the Japanese and Ryukyu Islands

The phylogeography of the endangered tideland snail Batillaria zonalis in the Japanese and Ryukyu Islands was analyzed on the basis of nucleotide sequences of a mitochondrial gene for cytochrome oxidase c subunit I (COI). Extremely low genetic diversity was found in populations at both the northern and southern boundaries of the geographic distribution of this species in Japan, i.e., Sendai Bay and Iriomotejima Island, respectively, which might be attributed to the population bottleneck due to historical environmental variations and/or the recent foundation of populations in the marginal part of the inhabitable range. Most populations contained unique rare haplotypes, and significant genetic differentiation on the whole was shown, while no clear geographic genetic structure was detected between the Japanese and Ryukyu Islands or over the distribution area of B. zonalis in Japan, with the exception of significant genetic divergences in Ago Bay in the central part of Honshu and the southern part of Okinawajima Island.     

Tracing the genetic roots of the sika deer Cervus nippon naturalized in Germany and Austria

To determine the geographical origin of the sika deer (Cervus nippon) naturalized in Germany and Austria, we sequenced the mitochondrial control region for 214 individuals. Adding these sequences to previously published data from native sika deer across its natural geographic range, the total comes to 245, extending what is already known about the geographical variation in this sequence in Cervus nippon. From these sequences, a neighbour-joining tree was constructed. This tree showed that the 49 different mitochondrial (mt)DNA types are grouped into three distinct phylogenetic clusters, which correspond to different geographic areas. Similarities between sequences of the naturalized sika deer and those described from native sika deer from both southern Honshu, Kyushu with associated islands, and northern Honshu suggest that the ancestors of the sika deer populations in Germany and Austria originated from the Japanese archipelago. In contrast, there is no evidence that female sika deer of Chinese, Taiwanese or north Vietnamese origin were involved in the ancestry of the present sika population in Germany and Austria.     

Genetic diversification of the Kanehira bitterling Acheilognathus rhombeus inferred from mitochondrial DNA,with comments on the phylogenetic relationship with its sister species Acheilognathus barbatulus

The Kanehira bitterling, Acheilognathus rhombeus, is a freshwater fish, discontinuously distributed in western Japan and the Korean Peninsula. Unusually among bitterling it is an autumn-spawning species and shows developmental diapause. Consequently, the characterization of its evolutionary history is significant not only in the context of the fish assemblage of East Asia, but also for understanding life-history evolution. This study aimed to investigate the phylogeography of A. rhombeus and its sister species Acheilognathus barbatulus, distributed in China, using a mitochondrial analysis of the ND1 gene from 311 samples collected from 50 localities in Japan and continental Asia. Phylogenetic analysis revealed that A. barbatulus is included in A. rhombeus and genetically closer to Japanese A. rhombeus than to Korean A. rhombeus. Divergence of Korean A. rhombeus and A. barbatulus from Japanese A. rhombeus was estimated to be from the late Pliocene (3.44 Mya) and the early Pleistocene (1.98 Mya), respectively. Each event closely coincided with the time of the Japan Sea opening. Japanese A. rhombeus comprised seven lineages: three in Honshu and four in Kyushu. One lineage in central Kyushu was genetically closer to the Honshu lineages than to other lineages in northern Kyushu. Divergence of Japanese lineages was estimated to be from the early to middle Pleistocene (0.55–0.93 Mya), during a period of geological and paleoclimatic change, including volcanic activity. Population expansion in the late Pleistocene (<0.10 Ma) was suggested in many of the lineages, which accords with other freshwater fishes. Biogeographically the ancestral A. rhombeus/A. barbatulus was likely to have repeatedly colonized Japan from the continent through land bridges in the late Pliocene and the early Pleistocene. However, the close genetic relationship between Japanese A. rhombeus and A. barbatulus suggests another possibility, with the second colonization occurring in reverse, from Japan to China. The small genetic distance between them indicates that the colonization occurred later than colonization events of other freshwater fishes, including other bitterling species.