Estimating the molecular evolutionary rates of mitochondrial genes referring to Quaternary ice age events with inferred population expansions and dispersals in Japanese Apodemus.

BACKGROUND: Determining reliable evolutionary rates of molecular markers is essential in illustrating historical episodes with phylogenetic inferences. Although emerging evidence has suggested a high evolutionary rate for intraspecific genetic variation, it is unclear how long such high evolutionary rates persist because a recent calibration point is rarely available. Other than using fossil evidence, it is possible to estimate evolutionary rates by relying on the well-established temporal framework of the Quaternary glacial cycles that would likely have promoted both rapid expansion events and interisland dispersal events. RESULTS: We examined mitochondrial cytochrome b (Cytb) and control region (CR) gene sequences in two Japanese wood mouse species, Apodemus argenteus and A. speciosus, of temperate origin and found signs of rapid expansion in the population from Hokkaido, the northern island of Japan. Assuming that global warming after the last glacial period 7-10 thousand years before present (kyr BP) was associated with the expansion, the evolutionary rates (sites per million years, myr) of Cytb and CR were estimated as 11-16 % and 22-32 %, respectively, for A. argenteus, and 12-17 % and 17-24 %, respectively, for A. speciosus. Additionally, the significant signature of rapid expansion detected in the mtDNA sequences of A. speciosus from the remaining southern main islands, Honshu, Shikoku, and Kyushu, provided an estimated Cytb evolutionary rate of 3.1 %/site/myr under the assumption of a postglacial population expansion event long ago, most probably at 130 kyr BP. Bayesian analyses using the higher evolutionary rate of 11-17 %/site/myr for Cytb supported the recent demographic or divergence events associated with the Last Glacial Maximum. However, the slower evolutionary rate of 3.1 %/site/myr would be reasonable for several divergence events that were associated with glacial periods older than 130 kyr BP. CONCLUSIONS: The faster and slower evolutionary rates of Cytb can account for divergences associated with the last and earlier glacial maxima, respectively, in the phylogenetic inference of murine rodents. The elevated evolutionary rate seemed to decline within 100,000 years.

Spatial and temporal aspects of occurrence of Mogera species in the Japanese islands inferred from mitochondrial and nuclear gene sequences.

We assessed dispersal and vicariant events in four species of Japanese moles in the genera Mogera and Euroscaptor to better understand the factors shaping intra- and interspecific differentiation in Japanese moles. We used the combined viewpoints of molecular phylogeny and historical geology using nucleotide sequences of mitochondrial (cytochrome b; Cytb) and nuclear (A2ab, Bmp4, Tcf25, vWf) genes. The divergence times estimated from the molecular data were verified with available geological data on the chronology of fluctuations in sea level in the Korea Strait, assuming sequential migration and speciation events. This produced possible migration times of 5.6, 3.5, 2.4, and 1.3 million years ago for four species of Japanese moles, Euroscaptor mizura, Mogera tokudae, M. imaizumii, and M. wogura, respectively. For the western Japanese mole M. wogura, Cytb sequences revealed four major phylogroups with strong geographic affinities in southwestern Central Honshu (I), western Honshu/Shikoku (II), Kyushu/westernmost Honshu (III), and Korea/Russian Primorye (IV). The nuclear gene sequences supported the distinctiveness of phylogroups I and IV, indicating long, independent evolutionary histories. In contrast, phylogroups II and III were merged into a single geographic group based on the nuclear gene data. Intraspecific divergences in M. imaizumii and M. tokudae were rather apparent in Cytb but not in nuclear gene sequences. The results suggest that repeated dispersal events have occurred between the Asian continent and the Japanese Islands, and intensive vicariant events associated with abiotic and biotic factors have created higher levels of species and genetic diversities in moles occurring on the Japanese Islands.

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.

Influence of Quaternary environmental changes on mole populations inferred from mitochondrial sequences and evolutionary rate estimation.

Quaternary environmental changes fundamentally influenced the genetic diversity of temperate-zone terrestrial animals, including those in the Japanese Archipelago. The genetic diversity of present-day populations is taxon- and region-specific, but its determinants are poorly understood. Here, we analyzed cytochrome b gene (Cytb) sequences (1140 bp) of mitochondrial DNA (mtDNA) to elucidate the factors determining the genetic variation in three species of large moles: Mogera imaizumii and Mogera wogura, which occur in central and southern mainland Japan (Honshu, Shikoku, and Kyushu), and Mogera robusta, which occurs on the nearby Asian continent. Network construction with the Cytb sequences revealed 10 star-shaped clusters with apparent geographic affinity. Mismatch distribution analysis showed that modes of pairwise nucleotide differences (τ values) were grouped into five classes in terms of the level, implying the occurrence of five stages for rapid expansion. It is conceivable that severe cold periods and subsequent warm periods during the late Quaternary were responsible for the population expansion events. The first and third oldest events included island-derived haplotypes, indicative of the involvement of land bridge formation between remote islands, hence suggesting an association of the ends of the penultimate (PGM, ca. 130,000 years ago) and last (LGM, ca. 15,000 years ago) glacial maxima, respectively. Since the third event was followed by the fourth, it is plausible that the termination of the Younger Dryas and subsequent abrupt warming ca. 11,500 years ago facilitated the fourth expansion event. The second event most likely corresponded to early marine isotope stage (MIS) 3 (ca. 53,000 years ago) when the glaciation and subsequent warming period were predicted to have influenced biodiversity. Utilization of the critical times of 130,000, 53,000, 15,000, and 11,500 years ago as calibration points yielded evolutionary rates of 0.03, 0.045, 0.10 and 0.10 substitutions/site/million years, respectively, showing a time-dependent manner whose pattern was similar to that seen in small rodents reported in our previous studies. The age of the fifth expansion event was calculated to be 5800 years ago with a rate of 0.10 substitutions/site/million years ago during the mid-Holocene, suggestive of the influence of humans or other unspecified reasons, such as the Jomon marine transgression.

Rodent-Borne Orthohantaviruses in Vietnam, Madagascar and Japan.

Hantaviruses are harbored by multiple small mammal species in Asia, Europe, Africa, and the Americas. To ascertain the geographic distribution and virus-host relationships of rodent-borne hantaviruses in Japan, Vietnam, Myanmar, and Madagascar, RNAlater™-preserved lung tissues of 981 rodents representing 40 species, collected in 2011-2017, were analyzed for hantavirus RNA by RT-PCR. Our data showed Hantaan orthohantavirus Da Bie Shan strain in the Chinese white-bellied rat (Niviventer confucianus) in Vietnam, Thailand; orthohantavirus Anjo strain in the black rat (Rattus rattus) in Madagascar; and Puumala orthohantavirus Hokkaido strain in the grey-sided vole (Myodes rufocanus) in Japan. The Hokkaido strain of Puumala virus was also detected in the large Japanese field mouse (Apodemus speciosus) and small Japanese field mouse (Apodemus argenteus), with evidence of host-switching as determined by co-phylogeny mapping.

Detection of recombinant haplotypes in wild mice (Mus musculus) provides new insights into the origin of Japanese mice.

Japanese house mice (Mus musculus molossinus) are thought to be a hybrid lineage derived from two prehistoric immigrants, the subspecies M. m. musculus of northern Eurasia and M. m. castaneus of South Asia. Mice of the western European subspecies M. m. domesticus have been detected in Japanese ports and airports only. We examined haplotype structuring of a 200 kb stretch on chromosome 8 for 59 mice from throughout Eurasia, determining short segments (approximately 370-600 bp) of eight nuclear genes (Fanca, Spire2, Tcf25, Mc1r, Tubb3, Def8, Afg3l1 and Dbndd1) which are intermittently arranged in this order. Where possible we identified the subspecies origin for individual gene alleles and then designated haplotypes for concatenated alleles. We recovered 11 haplotypes among 19 Japanese mice examined, identified either as 'intact' haplotypes derived from the subspecies musculus (57.9%), domesticus (7.9%), and castaneus (2.6%), or as 'recombinant' haplotypes (31.6%). We also detected recombinant haplotypes unique to Sakhalin. The complex nature of the recombinant haplotypes suggests ancient introduction of all three subspecies components into the peripheral part of Eurasia or complicated genomic admixture before the movement from source areas. 'Intact'domesticus and castaneus haplotypes in other Japanese wild mice imply ongoing stowaway introductions. The method has general utility for assessing the history of genetic admixture and for disclosing ongoing genetic contamination.

Genetic Diversity and Phylogeography of Thottapalayam thottimvirus (Hantaviridae) in Asian House Shrew (Suncus murinus) in Eurasia.

Murid and cricetid rodents were previously believed to be the principal reservoir hosts of hantaviruses. Recently, however, multiple newfound hantaviruses have been discovered in shrews, moles, and bats, suggesting a complex evolutionary history. Little is known about the genetic diversity and geographic distribution of the prototype shrew-borne hantavirus, Thottapalayam thottimvirus (TPMV), carried by the Asian house shrew (Suncus murinus), which is widespread in Asia, Africa, and the Middle East. Comparison of TPMV genomic sequences from two Asian house shrews captured in Myanmar and Pakistan with TPMV strains in GenBank revealed that the Myanmar TPMV strain (H2763) was closely related to the prototype TPMV strain (VRC66412) from India. In the L-segment tree, on the other hand, the Pakistan TPMV strain (PK3629) appeared to be the most divergent, followed by TPMV strains from Nepal, then the Indian-Myanmar strains, and finally TPMV strains from China. The Myanmar strain of TPMV showed sequence similarity of 79.3-96.1% at the nucleotide level, but the deduced amino acid sequences showed a high degree of conservation of more than 94% with TPMV strains from Nepal, India, Pakistan, and China. Cophylogenetic analysis of host cytochrome b and TPMV strains suggested that the Pakistan TPMV strain was mismatched. Phylogenetic trees, based on host cytochrome b and cytochrome c oxidase subunit I genes of mitochondrial DNA, and on host recombination activating gene 1 of nuclear DNA, suggested that the Asian house shrew and Asian highland shrew (Suncus montanus) comprised a species complex. Overall, the geographic-specific clustering of TPMV strains in Asian countries suggested local host-specific adaptation. Additional in-depth studies are warranted to ascertain if TPMV originated in Asian house shrews on the Indian subcontinent.

Molecular cloning and characterization of the repetitive DNA sequences that comprise the constitutive heterochromatin of the A and B chromosomes of the Korean field mouse (Apodemus peninsulae, Muridae, Rodentia).

Three novel families of repetitive DNA sequences were molecularly cloned from the Korean field mouse (Apodemus peninsulae) and characterized by chromosome in-situ hybridization and filter hybridization. They were all localized to the centromeric regions of all autosomes and categorized into major satellite DNA, type I minor, and type II minor repetitive sequences. The type II minor repetitive sequence also hybridized interspersedly in the non-centromeric regions. The major satellite DNA sequence, which consisted of 30 bp elements, was organized in tandem arrays and constituted the majority of centromeric heterochromatin. Three families of repetitive sequences hybridized with B chromosomes in different patterns, suggesting that the B chromosomes of A. peninsulae were derived from A chromosomes and that the three repetitive sequences were amplified independently on each B chromosome. The minor repetitive sequences are present in the genomes of the other seven Apodemus species. In contrast, the major satellite DNA sequences that had a low sequence homology are present only in a few species. These results suggest that the major satellite DNA was amplified with base substitution in A. peninsulae after the divergence of the genus Apodemus from the common ancestor and that the B chromosomes of A. peninsulae might have a species-specific origin.

Late Pleistocene climate change and population dynamics of Japanese Myodes voles inferred from mitochondrial cytochrome b sequences.

The Japanese archipelago is comprised of four main islands-Hokkaido, Honshu, Shikoku, and Kyushu-which contain high mountainous areas that likely allowed for lineage differentiation and population genetic structuring during the climatic changes of the late Pleistocene. Here, we assess the historical background of the evolutionary dynamics of herbivorous red-backed voles (Myodes) in Japan, examining the evolutionary trends of mitochondrial cytochrome b gene (Cytb) sequence variation. Four apparent signals from rapid expansion events were detected in three species, M. rufocanus and M. rutilus from Hokkaido and M. smithii from central Honshu. Taken together with results from previous studies on Japanese wood mice (Apodemus spp.), three of the expansion events were considered to be associated with predicted bottleneck events at the marine isotope stage (MIS) 4 period, in which glaciers are thought to have expanded extensively, especially at higher elevations. In the late Pleistocene, the possible candidates are transitions MIS 6/5, MIS 4/3, and MIS 2/1, which can be characterized by the cold periods of the penultimate glacial maximum, MIS 4, and the last glacial maximum, respectively. Our data further reveal the genetic footprints of repeated range expansion and contraction in the northern and southern lineages of the vole species currently found in central Honshu, namely M. andersoni and M. smithii, in response to climatic oscillation during the late Pleistocene. The time-dependent evolutionary rates of the mitochondrial Cytb presented here would provide a possible way for assessing population dynamics of cricetid rodents responding to the late Pleistocene environmental fluctuation.

Molecular Phylogeny of Mobatviruses (Hantaviridae) in Myanmar and Vietnam.

The discovery of highly divergent lineages of hantaviruses (family Hantaviridae) in shrews, moles, and bats of multiple species raises the possibility that non-rodent hosts may have played a significant role in their evolutionary history. To further investigate this prospect, total RNA was extracted from RNAlater®-preserved lung tissues of 277 bats (representing five families, 14 genera and 40 species), captured in Myanmar and Vietnam during 2013⁻2016. Hantavirus RNA was detected in two of 15 black-bearded tomb bats (Taphozous melanopogon) and two of 26 Pomona roundleaf bats (Hipposideros pomona) in Myanmar, and in three of six ashy leaf-nosed bats (Hipposideros cineraceus) in Vietnam. Pair-wise alignment and comparison of coding regions of the S, M, and L segments of hantaviruses from Taphozous and Hipposideros bats revealed high nucleotide and amino acid sequence similarities to prototype Láibin virus (LAIV) and Xuân SÆ¡n virus (XSV), respectively. Phylogenetic analyses, generated by maximum-likelihood and Bayesian methods, showed a geographic clustering of LAIV strains from China and Myanmar, but not of XSV strains from China and Vietnam. These findings confirm that the black-bearded tomb bat is the natural reservoir of LAIV, and that more than one species of Hipposideros bats can host XSV.

The evolution of female mole ovotestes evidences high plasticity of mammalian gonad development.

Previous studies of the reproductive biology and genetics of European moles (Talpa spp.) showed that all females of these species have ovotestes (gonads with testicular and ovarian tissue) instead of normal ovaries, a unique specialization among mammals. Females are fertile as their ovarian tissue is fully functional. Testicular tissue is abnormal and sterile, but produces high levels of testosterone. This phenomenon also characterizes other talpid species from Europe and North America. To study the origin of this singular reproductive specialization, we examined the gonads of several female specimens belonging to two critical taxa. Although large Japanese moles (Mogera wogura) posses ovotestes, greater Japanese shrew moles (Urotrichus talpoides) are characterized by normal ovaries. The results fit parsimoniously with a recent phylogenetic study that places Urotrichus relatively basal in the talpid tree and separate from the American shrew mole. Parsimony reconstruction on alternative phylogenetic hypotheses clearly indicates that reversal(s) must have occurred and suggests that a relatively simple genetic mechanism must be associated with the evolution of female hermaphroditism in moles.

Mapping of rRNA gene loci in the mice, Mus musculus molossinus (Japan) and Mus musculus musculus (Russia) by double color FISH.

Ribosomal RNA gene (rDNA) loci of Russian Mus musculus musculus and of Japanese Mus musuculus molossinus were mapped by double color FISH. The total number of rDNA loci was varied from 5 to 12, although the loci on chromosomes 12, 15, 16, 18, and 19 were common to all mice examined. Instead, polymorphisms of the rDNA loci were found on chromosomes 1, 5, 8, 9, 10, 11, 13 and 17. The novel rDNA loci of M. m. musculus were found in Nov/TUA strain on chromosomes 8 and 17. These observations, together with those of previous reports, suggest that the rDNA loci of Mus musculus species are in the evolutionary process of further translocation to other chromosomes.

Risk of accidental invasion and expansion of allochthonous mice in Tokyo metropolitan coastal areas in Japan.

House mouse (Mus musculus) is one of the perilous animal vectors for imported zoonosis such as a lymphocytic choriomeningitis (LCMV) infectious disease, and probably unknown emerging and/or re-emerging infectious diseases as well. It is necessary to prevent such diseases by regular surveys for behavioral trends of these allochthonous mice. However, such a trial has never been attempted in Japan. From 1998 to 2002, we analyzed partial sequences of the D-loop region in mtDNA, which provides powerful diagnostic SNPs for subspecies identification in the Mus musculus species, from 301 individuals of mice collected in 23 international bays or airports in Japan. We found that invasion of many allochthonous mice, which were identified as European subspecies, Mus musculus domesticus, occurred in Tokyo metropolitan coastal area. Based on the evidence, we warn that extensive invasion of allochthonous mice has occurred recently and, therefore, the risk of emerging and/or re-emerging infectious diseases invasion might be high in Tokyo metropolitan area.

Impairment of endothelium-dependent relaxation of aortas and pulmonary arteries from spontaneously hyperlipidemic mice (Apodemus sylvaticus).

We evaluated the endothelial function of thoracic aortas and pulmonary arteries in a population of European wood mice (Apodemus sylvaticus), which exhibit hypercholesterolemia. According to the plasma cholesterol level, mice were divided into two groups: hypercholesterolemic (AHL, total plasma cholesterol 200-300 mg/dl) and normocholesterolemic (ANL, total plasma cholesterol <200 mg/dl). Acetylcholine (ACh) caused endothelium-dependent relaxation of precontracted aortas and pulmonary arteries. Relaxation of the pulmonary artery is completely dependent on nitric oxide. This relaxation was inhibited in AHL pulmonary arteries. On the other hand, part of the ACh-induced relaxation of the thoracic aorta was resistant to N(omega)-nitro-L-arginine (L-NNA). L-NNA-sensitive and -resistant relaxation to ACh were also inhibited in AHL aortas. Inhibition of endothelium-dependent relaxation of the aortas was correlated with total plasma cholesterol level. Endothelium-independent relaxation to sodium nitroprusside (SNP) was similar in AHL and ANL pulmonary arteries, but in the thoracic aorta of AHL mice, the sensitivity to SNP was slightly decreased, without a change in maximal response to SNP. No morphological change was observed in the aortas and the pulmonary arteries from AHL and ANL mice. Thus, AHL mice are valuable as a new experimental model to study the relation of hyperlipidemia to vascular disease since the endothelial function is impaired in these mild hyperlipidemic animals.

Intraspecific differentiation in the lesser Japanese mole in eastern Honshu, Japan, indicated by nuclear and mitochondrial gene analyses.

The lesser Japanese mole, Mogera imaizumii, recognized by Motokawa and Abe (1996), occurs in eastern Honshu, western Honshu and Shikoku. Mitochondrial and nuclear DNA were analyzed for mole samples from eastern Honshu to elucidate intraspecific differentiation. Analyses of sequences of the mitochondrial cytochrome b gene (Cytb) and of a restriction fragment length polymorphism of the nuclear 28S ribosomal RNA gene spacer (rDNA-RFLP) revealed two genetic types, partially corresponding to Hutterer's (1993) taxa, M. wogura (= M. imaizumii) and M. minor. Most samples showed either of two combinations of mitochondrial / nuclear gene types. However, two specimens showed a different combination. This incongruent combination of mitochondrial and nuclear genes might have derived, in part, from an introgression event between genetically differentiated populations after secondary contact during the evolutionary history of the lesser Japanese mole in eastern Honshu.

Preliminary genetic characterization of two lineages of black rats (Rattus rattus sensu lato) in Japan, with evidence for introgression at several localities.

We conducted a pilot survey of genetic diversity among 37 karyotyped individuals of the black rat Rattus rattus (sensu lato) from six localities on the Japanese Islands, using complete gene sequences of mitochondrial cytochrome b (cyt b) and nuclear interphotoreceptor retinoid binding protein (IRBP). Our sampling included two previously documented karyotypic groups: 'Oceanian' with 2n = 38 and 'Asian' with 2n = 42. Cyt b sequences for most individuals clustered according to their karyotypic groups, with an average between-group divergence of 3.8%. One exception was that individuals from Kagoshima (Kyushu Island) showed 'Asian' karyotypes combined with a cyt b haplotype that differed by a single nucleotide substitution from the haplotype of the 'Oceanian' karyotypic group. Six IRBP haplotypes were identified. They belonged to three distinct IRBP lineages (I-III), with an average inter-lineage divergence of 1%. Among homozygous individuals, these lineages showed good association with the karyotypic groups: IRBP lineage I occurred only with 'Oceanian' karyotypes, while IRBP lineages II and III both occurred with 'Asian' karyotypes. Individuals from Kagoshima all possessed IRBP of 'Asian' lineages, despite the presence of an 'Oceanian' mitochondrial type. The Chichijima population (Ogasawara Islands) featured exclusively 'Asian' karyotypes and cyt b sequences, but various combinations of all three IRBP lineages. The Kagoshima and Chichijima populations thus provide strong evidence of viable hybridization and genetic introgression between the two karyotypic groups, but with variable genetic outcomes. Our results demonstrate the potential of combined analysis of karyotypes and mitochondrial and nuclear gene sequences to elucidate the complex dispersal and population history of the black rat.

Abnormal pairing of X and Y sex chromosomes during meiosis I in interspecific hybrids of Phodopus campbelli and P. sungorus.

Hybrid sterility plays an important role in the maintenance of species identity and promotion of speciation. Male interspecific hybrids from crosses between Campbell's dwarf hamster (Phodopus campbelli) and the Djungarian hamster (P. sungorus) exhibit sterility with abnormal spermatogenesis. However, the meiotic phenotype of these hybrids has not been well described. In the present work, we observed the accumulation of spermatocytes and apoptosis of spermatocyte-like cells in the testes of hybrids between P. campbelli females and P. sungorus males. In hybrid spermatocytes, a high frequency of asynapsis of X and Y chromosomes during the pachytene-like stage and dissociation of these chromosomes during metaphase I (MI) was observed. No autosomal univalency was observed during pachytene-like and MI stages in the hybrids; however, a low frequency of synapsis between autosomes and X or Y chromosomes, interlocking and partial synapsis between autosomal pairs, and γ-H2AFX staining in autosomal chromatin was observed during the pachytene-like stage. Degenerated MI-like nuclei were frequently observed in the hybrids. Most of the spermatozoa in hybrid epididymides exhibited head malformation. These results indicate that the pairing of X and Y chromosomes is more adversely affected than that of autosomes in Phodopus hybrids.

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.

Molecular phylogeny of Japanese Rhinolophidae based on variations in the complete sequence of the mitochondrial cytochrome b gene.

Microchiroptera have diversified into many species whose size and the shapes of the complicated ear and nose have been adapted to their echolocation abilities. Their speciation processes, and intra- and interspecies relationships are still under discussion. Here we report on the geographical variation of Japanese Rhinolophus ferrumequinum and R. cornutus using the complete sequence of the mitochondrial cytochrome b gene to clarify the phylogenetic positions of the 2 species as well as that of Rhinolophidae within the Microchiroptera. We have found that sequence divergence values within each of the 2 species are unexpectedly low (0.07%-0.94%). We have also found that there is no local specificity of their mtCytb alleles. On the other hand, the divergence values for Japanese Microchiroptera (12.7%-16.6%) are much higher than those for other mammalian genera. Similarly, the values among five genera of Vespertilionidae were 20.5%-27.3%. Phylogenetic analysis shows that the 2 species of family Rhinolophidae in the suborder Microchiroptera belong to the Megachiroptera cluster in the constructed maximum parsimony tree. These results suggest that the speciation of Rhinolophidae involved its divergence as an independent lineage from other Microchiroptera, and other microbats might be paraphyletic. In addition, the tree also shows that the order Chiroptera is monophylitic, and the closest group to Chiroptera is the ungulates.

Geographic variation and diversity of the cytochrome b gene in Japanese wild populations of medaka, Oryzias latipes.

We conducted a polymerase chain reaction--restriction fragment length polymorphism (PCR-RFLP) analysis of the mitochondrial cytochrome b gene to elucidate the detailed genetic population structure of Japanese wild populations of medaka, Oryzias latipes. The analysis of 1,225 specimens collected from 303 sites identified 67 mitotypes. Subsequently we determined the nucleotide sequences of the complete cytochrome b gene (1141-bp) to clarify the phylogenetic relationships among mitotypes. The phylogenetic tree based on nucleotide sequences indicated three major clades (A, B and C) that differed by 11.3-11.8%, corresponding to three clusters previously identified by RFLP analysis of entire mitochondrial DNAs. The geographic distribution of mitotypes in clades A and B was fully concordant with the Northern and Southern Populations defined by allozymes. Clade A could be subdivided into three subclades and clade B into eleven, with sequence divergences among subclades of 1.3-5.8%. Each distribution of mitotypes in subclades roughly corresponded to that of mtDNA haplotypes in subclusters previously identified. Mitotypes in clade C were found only in the Kanto district. The phylogenetic relationships and the estimated divergence times suggest that three Japanese clades originated from a common ancestor and were separated during the Pliocene, and that the regional differentiation of subclades was closely connected with the geological history of the Quaternary. This study has also demonstrated the possibility of artificial disturbance of natural distribution especially in the Kanto district and the superior efficacy of PCR-RFLP analysis as a simple method for detecting genetic variation and artificial gene flow of medaka.