Sex Differences in Polymorphism and Expression of AAT-1 in the Hiroshima Population of Buergeria buergeri (Anura, Rhacophoridae)

Buergeria buergeri is female heterozygous in sex determination; chromosome pair No. 7 in this species is a pair of sex chromosomes of the ZZ/ZW type. Genetic analysis of AAT-1 variants was carried out to elucidate the mode of inheritance of this locus by starch-gel electrophoresis using field-caught females and males and their offspring produced by artificial crossings. The results showed that the AAT-1 locus is sex-linked and that alleles are expressed on the Z chromosome, but not the W chromosome. It is evident that the AAT-1 gene of female offspring is hemizygous and that the allele present is on the Z chromosome, which is derived from the male parent.     

The origin and differentiation of the heteromorphic sex chromosomes Z, W, X, and Y in the frog Rana rugosa, inferred from the sequences of a sex-linked gene, ADP/ATP translocase

Sex chromosomes of the Japanese frog Rana rugosa are heteromorphic in the male (XX/XY) or in the female (ZZ/ZW) in two geographic forms, whereas they are still homomorphic in both sexes in two other forms (Hiroshima and Isehara types). To make clear the origin and differentiation mechanisms of the heteromorphic sex chromosomes, we isolated a sex-linked gene, ADP/ATP translocase, and constructed a phylogenetic tree of the genes derived from the sex chromosomes. The tree shows that the Hiroshima gene diverges first, and the rest form two clusters: one includes the Y and Z genes and the other includes the X, W, and Isehara genes. The Hiroshima gene shares more sequence similarity with the Y and Z genes than with the X, W, and Isehara genes. This suggests that the Y and Z sex chromosomes originate from the Hiroshima type, whereas the X and W chromosomes originate from the Isehara-type sex chromosome. Thus, we infer that hybridization between two ancestral forms, with the Hiroshima-type sex chromosome in one and the Isehara-type sex chromosome in the other, was the primary event causing differentiation of the heteromorphic sex chromosomes.      

Low genetic diversities in isolated populations of the Asian black bear (<Emphasis Type="Italic">Ursus thibetanus</Emphasis>) in Japan,in comparison with large stable populations

Populations of the Asian black bear (Ursus thibetanus) are relatively large and continuous in central Honshu, the main island of Japan, but they are isolated in western Honshu. To clarify the degree of genetic isolation of the populations in western Honshu, we compared the genetic diversities of four populations in western Honshu with that of one of the continuous populations of central Honshu. Three of the four western Honshu populations were isolated and the other was continuous with the central Honshu populations on a geographical distribution basis. The genotypes at 10 microsatellite loci of the sampled individuals were determined and the genetic structures of the populations examined. Genetic diversities were significantly lower in the isolated populations than in the continuous populations. The continuous population in central Honshu had high levels of genetic diversity, comparable to those in populations of the American black bear (Ursus americanus) and the brown bear (Ursus arctos). The genetic distances between the two continuous populations were smallest, even though their geographic distance was largest (>200 km) among all the pairs of neighboring populations examined. Low genetic diversity within the isolated populations suggested genetic drift due to the small population size; the genetic differentiation among the populations indicated low rates of gene flow among them.     

Y chromosomal DNA variation and the peopling of Japan.

Four loci mapping to the nonrecombining portion of the Y chromosome were genotyped in Japanese populations from Okinawa, the southernmost island of Japan; Shizuoka and Aomori on the main island of Honshu; and a small sample of Taiwanese. The Y Alu polymorphic (YAP) element is present in 42% of the Japanese and absent in the Taiwanese, confirming the irregular distribution of this polymorphism in Asia. Data from the four loci were used to determine genetic distances among populations, construct Y chromosome haplotypes, and estimate the degree of genetic diversity in each population and on different Y chromosome haplotypes. Evolutionary analysis of Y haplotypes suggests that polymorphisms at the YAP (DYS287) and DXYS5Y loci originated a single time, whereas restriction patterns at the DYS1 locus and microsatellite alleles at the DYS19 locus arose more than once. Genetic distance analysis indicated that the Okinawans are differentiated from Japanese living on Honshu. The data support the hypotheses that modern Japanese populations have resulted from distinctive genetic contributions involving the ancient Jomon people and Yayoi immigrants from Korea or mainland China, with Okinawans experiencing the least amount of admixture with the Yayoi. It is suggested that YAP+ chromosomes migrated to Japan with the Jomon people > 10,000 years ago and that a large infusion of YAP- chromosomes entered Japan with the Yayoi migration starting 2,300 years ago. Different degrees of genetic diversity carried by these two ancient chromosomal lineages may be explained by the different life-styles (hunter-gatherer versus agriculturalist). of the migrant groups, the size of the founding populations, and the antiquities of the founding events.     

Evolutionary and biogeographical implications of the karyological variations in the oviparous and viviparous forms of the lizard Lacerta (Zootoca) vivipara

The lizard Lacerta ( Zootoca ) vivipara has two modes of reproduction and is variable karyologically. We describe its karyological variation from literature data and from new data on two viviparous populations from France, on two oviparous populations from the Pyrenees in south-western France and on three oviparous populations recently discovered in Slovenia. Males have 36 chromosomes, whereas females have only 35 chromosomes in all viviparous populations and in the Pyrenean oviparous populations. This karyotype has been interpreted to result from a fusion of an ancestral sexual W chromosome with an autosome from the Zl or from the Z2 pair. The karyotype formula is 32 autosomes + ZIZ2W for the female and 32 autosomes + Z1Z1Z2Z2 for the male. The karyotype of the Slovenian oviparous populations, 34 autosomes + ZW in the male and 34 autosomes +ZW in the female, represents an evolutionary stage that preceded the chromosomal fusion. There is minor karyological variation, mainly concerning the W and Z2 chromosomes, within the Pyrenean oviparous populations. This parallels the geographic variation of the W-linked alleles of the MPI enzyme and suggests that allopatric differentiation of these oviparous populations might have occurred in the vicinity of the Pyrenees during the Pleistocene.
The viviparous populations from western Europe carry a metacentric W chromosome, whereas oviparous populations from south-western Europe and eastern viviparous populations both show an acrocentrie, or a subtelocentrie. W chromosome. This suggests that the acrocentric-subtelocentric W is a primitive character and that viviparity probably arose in an eastern lineage of the species.     

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.     

Reduced variation on the chicken Z chromosome

Understanding the population genetic factors that shape genome variability is pivotal to the design and interpretation of studies using large-scale polymorphism data. We analyzed patterns of polymorphism and divergence at Z-linked and autosomal loci in the domestic chicken (Gallus gallus) to study the influence of mutation, effective population size, selection, and demography on levels of genetic diversity. A total of 14 autosomal introns (8316 bp) and 13 Z-linked introns (6856 bp) were sequenced in 50 chicken chromosomes from 10 highly divergent breeds. Genetic variation was significantly lower at Z-linked than at autosomal loci, with one segregating site every 39 bp at autosomal loci (theta(W) = 5.8 +/- 0.8 x 10(-3)) and one every 156 bp on the Z chromosome (theta(W) = 1.4 +/- 0.4 x 10(-3)). This difference may in part be due to a low male effective population size arising from skewed reproductive success among males, evident both in the wild ancestor-the red jungle fowl-and in poultry breeding. However, this effect cannot entirely explain the observed three- to fourfold reduction in Z chromosome diversity. Selection, in particular selective sweeps, may therefore have had an impact on reducing variation on the Z chromosome, a hypothesis supported by the observation of heterogeneity in diversity levels among loci on the Z chromosome and the lower recombination rate on Z than on autosomes. Selection on sex-linked genes may be particularly important in organisms with female heterogamety since the heritability of sex-linked sexually antagonistic alleles advantageous to males is improved when fathers pass a Z chromosome to their sons.     

Sex and death: CHD1Z associated with high mortality in moorhens

Abstract.— Sex ratios in clutches of moorhens (Gallinula chloropus) in Britain were measured on 83 chicks using the sex-linked CHD1 gene (Chromo-helicase/ATPase-DNA binding protein 1). Among birds, the female is the hetero-gametic sex (Z and W chromosomes), and the male is homogametic (two copies of the Z chromosome). We report variation among the PCR-amplified fragments of the CHD1Z , and the death of nearly all heterozygous male chicks (92%). In contrast, survivorship among females and homozygote males was 54–60%. Mortality in male heterozygotes was significantly higher than that of male homozygotes (P < 0.001). Chick and egg biometrics were not significantly different between these males. The CHD1Z was unlikely to be directly responsible but may have been hitchhiked by the causal gene(s). The observations appear to follow a classic underdominance (heterozygote inferiority) pattern, but raise the paradoxical question of why one form of the Z chromosome has not been fixed, as is expected from evolutionary theory. We discuss possible explanations and include a survey of British populations based on skin specimens.     

Mitochondrial DNA variation among worldwide populations of gypsy moths, Lymantria dispar

Gypsy moth populations from Japan, mainland Asia, Europe, Tunisia, and North America were analyzed for variation in mitochondrial DNA (mtDNA) sequences from three gene regions. These samples resolve into four groups, representing gypsy moths from (1) Okinawa, Japan, (2) Hokkaido, Japan, (3) Honshu and Kyushu, Japan and mainland Asia, and (4) Europe, Tunisia, and North America. Some patterns of geographic variation observed for mtDNA (for example, the distinctiveness of gypsy moths from Hokkaido, Japan) coincide with those observed by Goldschmidt from analyses of morphology, life history, and intersexuality. Other patterns (relative sequence homogeneity across Asia, Honshu, and Kyushu and reduced levels of variation in mainland Japan) do not.     

Structural differences between XX and ZW sex lampbrush chromosomes inRana rugosa females (Anura: Ranidae)

In this study we investigated the morphology and pairing behavior of sex lampbrush chromosomes of XX and ZW females of Rana rugosa from five localities in Japan. Whereas lampbrush chromosomes of XX females from Hiroshima and Isehara had subterminally located centromeres and showed remarkable similarity, those of XX females from Hamakita had the centromeres in the middle. Analysis of landmark configurations revealed that chromosome Xq of Hamakita females closely resembled a part of Xq of Hiroshima and Isehara females, whereas Xp of Hamakita females was inverted compared with the other part of Xq of Hiroshima and Isehara females. Z chromosomes from Kanazawa and Niigata closely resembled the Hiroshima X, whereas the W closely resembled the Hamakita X. XX pairings from Hiroshima, Isehara, and Hamakita were found to be joined by one to four chiasmata at various points all along the axis in both the short and long arms, whereas chromosomal pairs from Kanazawa and Niigata showed only one chiasma between Zp and the distal region of Wq. From these findings we conclude that (1) both the W and the Hamakita X must have evolved from the more primitive Hiroshima and Isehara X chromosomes by a series of pericentric inversions; and (2) females distributed in Hamakita possess two X chromosomes similar to the W, suggesting that either sex-determining or sex-modifying genes on the Hamakita X are clearly different from those on the Kanazawa and Niigata W chromosome. Received: 27 February 1996; in revised form: 22 May 1996 / Accepted: 25 May 1996     

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.     

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.     

Geographical Differentiation in a Japanese Stream-Breeding Frog, Buergeria buergeri, Elucidated by Morphometric Analyses and Crossing Experiments

Hiroshima and Aomori populations of Buergeria buergeri (hereinafter abbreviated as HIROSHIMA and AOMORI, respectively) were morphologically differentiated in both sexually matured frogs and tadpoles. The mean snout-vent lengths of females were 67.4 mm in HIROSHIMA and 50.4 mm in AOMORI, and those of males were 42.9 mm in HIROSHIMA and 37.2 mm in AOMORI. The mean body weights of females of HIROSHIMA and AOMORI just after spawning were 18.0 g and 7.2 g, and those of males of HIROSHIMA and AOMORI were 4.9 g and 3.8 g, respectively. AOMORI tadpoles were rather stocky and their appearance seemed to be more adaptive to lentic water than HIROSHIMA tadpoles. The lower lip of the tadpoles at stage XIII consisted of 2 continuous and 1 broken teeth rows in HIROSHIMA, whereas 3 continuous and 1 broken ones in AOMORI. The growth rate of the embryos was higher in AOMORI than in HIROSHIMA. The embryos of AOMORI were more tolerant to high temperature and less tolerant to low temperature than those of HIROSHIMA. Hybrids between these two populations showed considerably reduced viability in either combination of reciprocal crosses.     

Geographic Variations in Reproductive Parameters of the Starspotted Dogfish, Mustelus manazo, from Five Localities in Japan and in Taiwan

Geographic variations in reproductive parameters of Mustelus manazo from five localities, four in Japan (Aomori, Tokyo Bay, Maizuru, and Shimonoseki) and one in Taiwan, were investigated from November 1994 to June 1996. Differences in age and length at sexual maturity from Aomori to Taiwan were approximately 3 years and 300mm in TL, respectively, for both sexes. The sharks appeared to mature at a later age and to a larger size in the most northern population, Aomori, and to be faster and smaller in the southern population, Taiwan. Mating, ovulation and fertilization periods were generally during May and August in the four localities of Tokyo Bay, Maizuru, Shimonoseki, and Taiwan. In Aomori, males showed no clear monthly variation in gonad index, and females in Aomori had a protracted reproductive period. Females became pregnant every year in all localities, except Aomori. The Aomori population apparently has a different reproductive cycle. The Tokyo Bay population is distinctive regarding growth and sex ratios of embryos. In all localities, the number of embryos per litter increased relative to total length of mother. There was no geographic variation in reproductive parameters between Maizuru and Shimonoseki. The coast between these last two sampling areas is a continuous steep slope with similar environmental conditions, and mixing of the two populations may occur.     

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.     

Management units of the endangered herb Primula sieboldii based on microsatellite variation among and within populations throughout Japan

To promote programs for the conservation and restoration of the endangered species Primula sieboldii, we examined genetic variation at eight microsatellite loci among and within 32 remnant wild populations throughout Japan. Total allelic diversity within a population was higher in larger populations, but not so after rarefaction adjustment. The positive relationship between population size and the inbreeding coefficient may suggest that more heterozygous genets tend to survive the habitat contraction possibly because of the higher fitness associated with heterozygosity. By principal coordinate analysis and Bayesian analysis, we detected four genetic groups (Hokkaido, northern Honshu, central Honshu, and western Japan), which could be recognized as management units of P. sieboldii. If supplementation with plants from other populations were planned, it should be conducted among populations which belong to the same management unit and which are likely to represent the same adaptive variation.     

Impact of geographical isolation on genetic differentiation in insular and mainland populations of Weigela coraeensis (Caprifoliaceae) on Honshu and the Izu Islands

Aim To provide insights into genetic differentiation between insular endemic Weigela coraeensis var. fragrans and its progenitor variety W. coraeensis var. coraeensis, the population genetic structure of both varieties was examined, and factors promoting genetic differentiation between the two taxa were explored. Location The natural range of W. coraeensis (sensu lato) throughout mainland Japan (Honshu) and the Izu Islands. Methods The analysis included 349 and 504 individuals across the mainland (Honshu) and the Izu Islands, respectively, using 10 allozyme and 10 microsatellite loci. The population genetic structure of W. coraeensis was assessed by analysing genetic diversity indices for each population, genetic differentiation among populations, model‐based Bayesian clustering or distance‐based clustering, and bottleneck tests. Results The level of genetic diversity in each of the populations on the Izu Islands was negatively correlated with geographical distance between each island and the mainland. The populations on the mainland and on the Izu Islands were genetically differentiated to a certain extent; however, the microsatellite analyses suggested that gene flow also occurred between the mainland and the islands, and among individual islands. These microsatellite analyses also suggested recent bottlenecks in several populations in both areas. Main conclusions The decrease in genetic diversity throughout the Izu Islands, which correlated with distance to the mainland, Honshu, may be the result of a repeated founder effect occurring at a series of inter‐island colonizations from north to south. The stepping stone‐like configuration of the islands may have played a role in the dispersal of the species. Geographical isolation by sea would effectively result in genetic differentiation of W. coraeensis between mainland Honshu and the Izu Islands, although some gene flow may still occur between Honshu and the northern Izu Islands. The differentiation process of the endemic plants on the Izu Islands is anagenetic but not completed, and the study of these plants will provide insightful knowledge concerning the evolution of insular endemics.     

Lineage admixture during postglacial range expansion is responsible for the increased gene diversity of Kalopanax septemlobus in a recently colonised territory

We aimed to reveal the effects of range expansion and subsequent lineage admixture from separated glacial refugia on genetic diversity of Kalopanax septemlobus in Japan, by combining nuclear microsatellite data and ecological niche modelling. Allelic richness and gene diversity were compared at the population and regional level. We also statistically examined these indices as a function of population accessibility to the last glacial maximum (LGM) palaeodistribution reconstructed by ecological niche modelling to test a simple range expansion scenario from glacial refugia. Genetic diversity was highest in the populations of southern Japan and gradually decreased towards the north. However, an additional centre of genetic diversity, when measured as gene diversity, was found in northern Honshu Island, where distinct lineages were shown to be in contact. Positive effects of population accessibility to the LGM range were detected in both diversity indices at different spatial scales. The combined data support independent postglacial range expansions towards the north from the edge populations on the exposed coastal shelf of Pacific and Sea of Japan in northern Honshu during the LGM, which subsequently resulted in markedly low genetic diversity in the northernmost extant range, Hokkaido. The regional increase in gene diversity in northern Honshu is likely to be the result of postglacial lineage admixture. Relative difference in the spatial scales best relating population genetic diversity with the LGM distribution can be explained by a higher rate of allelic richness diversity loss during range expansions and stronger effects of lineage admixture on gene diversity.     

Occurrence of the endangered species Nitellopsis obtusa (Charales,Charophyceae) in western Japan and the genetic differences within and among Japanese populations

Nitellopsis obtusa (Charales, Charophyceae) are widely distributed from Europe to Asia; however, this species has been recorded in only the five lakes in central Honshu in Japan. This species was thought to be extinct in Japan, but was rediscovered in limited areas of Lake Kawaguchi in central Honshu. More recently, we discovered more Japanese populations of N. obtusa in Lake Biwa in western Honshu, and it became clear that the species had a broader distribution in Japan than originally believed. In addition, although only male or female thalli have been collected at each lake, both male and female thalli were found from Lake Biwa. This is the first report of a potentially sexual population of N. obtusa in Japan. The DNA sequences of three chloroplast DNA markers, including both coding and non‐coding regions, were identical in all specimens from Lake Kawaguchi and Lake Nojiri (Central Honshu), and differed from those of Lake Biwa and German specimens. Although Japanese and German specimens were genetically similar, Japanese specimens displayed considerable genetic diversity according to locality.     

Phylogeography of the Japanese giant flying squirrel, Petaurista leucogenys (Rodentia: Sciuridae): implication of glacial refugia in an arboreal small mammal in the Japanese Islands

To test the association between temperate forest dynamics and glacial refugia for arboreal small mammals, we studied the phylogeography of the Japanese giant flying squirrel ( Petaurista leucogenys ) using complete mitochondrial cytochrome b gene sequences (1140 bp). This squirrel is endemic to three of Japan's main islands: Honshu, Shikoku, and Kyushu. We examined 58 specimens of P. leucogenys collected from 40 localities in Japan. Additionally, two individuals with unknown sampling localities were included in phylogenetic analyses. There were 54 haplotypes of P. leucogenys. We found five major phylogroups (Northern, Central, South-eastern, South-western, and Southern). These phylogroups may have originated from glacial refugia during the Late Pleistocene. After the last glaciation, the Northern phylogroup, widely distributed in eastern Japan, could have extensively expanded northward from its refugia. By contrast, in western Japan, population expansion was restricted to western Japan. All members of four phylogroups existed in western Japan during glaciations. The complicated phylogeographical pattern of P. leucogenys populations originating from western Japan may have resulted from the long history.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 98 , 47–60.