Estimation of Evolutionary Rates of Mitochondrial DNA in Two Japanese Wood Mouse Species Based on Calibrations with Quaternary Environmental Changes

Region-Specific Genetic Diversity of Black Rats (Rattus rattus Complex) in Southeast and East Asia Shaped by Rapid Population Expansion Events

Among the six mitochondrial DNA lineages of the black rat (Rattus rattus Complex; RrC), lineages II and IV are widespread in Southeast and East Asia. This study explored their demographic history using 17 new sequences from the Miyako Islands in the Ryukyu archipelago, together with 178 publicly available cytochrome b sequences. We defined six and two haplotype groups showing rapid population expansion signals in Lineages II and IV, respectively. The six haplotype groups of Lineage II were represented by haplotypes from 1) Myanmar/Bangladesh/Northeast India, 2) Laos, 3) Thailand, 4) Indonesia/Philippines, 5) Vietnam/southern China, and 6) the Ryukyu archipelago. These expansion times were estimated using time-dependent evolutionary rates to be 115,300 years ago (ya), 128,500 ya, 9600 ya, 10,600 ya, 7200 ya, and 1400 ya, respectively, although all had large confidence intervals. The two groups of Lineage IV were recovered from the mainland and islands of Southeast Asia with predicted expansion times of 197,000 ya and 5800 ya, respectively. These results suggest that climatic fluctuations during the last 200,000 years of the Quaternary, affected the population dynamics in subtropical areas at different times. Furthermore, the results of the younger rapid expansion events of RrC suggest the possibility of agricultural advancement and dispersal of Neolithic farmers to different areas within the mainland and islands of Southeast Asia during the Holocene. A subset of rats from the Miyako Islands were found to have the same lineage IV haplotypes as those in Southeast Asia, suggesting a recent introduction of these new lineages.

Phylogenetics and Population Genetics of the Asian House Shrew, Suncus murinus-S. montanus Species Complex, Inferred From Whole-Genome and Mitochondrial DNA Sequences, with Special Reference to the Ryukyu Archipelago, Japan

The house shrew (Suncus murinus-S. montanus species complex) colonized regions across southern Asia and the Indian Ocean following human activity. The house shrew is distributed on islands of the Ryukyu Archipelago, the southernmost part of Japan, but the evolutionary history of the shrew on those islands and possible associations between these populations and humans remain unknown. In this study, we conducted phylogenetic and population genetic analyses based on both nuclear and mitochondrial genome sequences of house shrews. Phylogenetic analyses based on mitochondrial cytochrome b (cytb) sequences revealed that shrews from the Ryukyu Archipelago showed strong genetic affinity to Vietnamese and southern Chinese shrews. Demographic analyses of cytb sequences indicated a rapid population expansion event affecting the haplotype group in Vietnam, southern China, and the Ryukyu Archipelago 3300-7900 years ago. Furthermore, gene flow between Ryukyu (Yonaguni Island) and Taiwan and between Ryukyu and Vietnam inferred from f4 statistics of the nuclear genomes suggested repeated immigration to Ryukyu in recent years. The present study demonstrates that the Nagasaki population has a different origin from the Ryukyu population. These findings elucidate the complex pattern of genetic admixture in house shrews and provide insights into their evolutionary history.

Temporal dynamics of mildly deleterious nonsynonymous substitutions in mitochondrial gene sequences in rodents and moles

We have previously estimated the evolutionary rate (number of substitutions/site/million years) of mitochondrial cytochrome b gene (Cytb) sequences in rodents and moles to be about 0.11 at more recent divergence times of tens of thousands of years, and to decrease rapidly to about 0.03 at more distant divergence times. Because this time dependency is thought to be caused by the removal of mildly deleterious substitutions in later generations, we focused in this study on the abundance of nonsynonymous substitutions. We collected 23 haplogroups of Cytb with signals of late Quaternary population expansion events from rodents and moles and categorized them into three groups for comparison based on predicted expansion start time: 5,000-15,000 years ago (Group I), ca. 53,000 years ago (Group II) and 130,000-230,000 years ago (Group III). We counted the numbers of nonsynonymous and synonymous substitutions in all haplogroups. The rates of nonsynonymous substitutions were lowest in Groups II and III (0.08-0.22), whereas those in Group I varied markedly. We further classified Group I into two subgroups based on high (0.29-0.43) and low (0.09-0.20) nonsynonymous substitution rates, which were likely to be associated with the start of the expansion within 10,000 years and at around 15,000 years ago, respectively. The Group II and III networks had two- or three-step star-shaped structures and tended to exhibit frequent and less frequent nonsynonymous substitutions on exterior and interior branches, respectively. Based on temporal dynamics, nonsynonymous mitochondrial DNA (mtDNA) substitutions in small mammals accounted for at most 40% of all substitutions during the early evolutionary stage and then rapidly declined, dropping to approximately 15%. The results of this study provide a good explanation of the time-dependent trend in the mtDNA evolution rate predicted in previous work.

Histological analysis of implanted embryos in large Japanese field mouse (Apodemus speciosus) and estimation of developmental stage

The large Japanese field mouse (Apodemus speciosus) is a small rodent species endemic to Japan. The genetic characteristics of A. speciosus include different chromosome numbers within the same species. Furthermore, A. speciosus has been used in radiation and genetic research. In the present study, a pregnant A. speciosus was obtained, and histochemical analysis of the implanted embryos was performed and compared with the developmental stages of the mouse (Mus musculus). Although there were some differences, the structures of the implanted embryos, including the primitive streak and placenta of A. speciosus were similar to those of mouse. Our study will be important for the construction of a developmental atlas of A. speciosus.

Late Pleistocene Expansion of Small Murid Rodents across the Palearctic in Relation to the Past Environmental Changes

We investigated the evolutionary history of the striped field mouse to identify factors that initiated its past demographic changes and to shed light on the causes of its current genetic structure and trans-Eurasian distribution. We sequenced mitochondrial cyt b from 184 individuals, obtained from 35 sites in central Europe and eastern Mongolia. We compared genetic analyses with previously published historical distribution models and data on environmental and climatic changes. The past demographic changes displayed similar population trends in the case of recently expanded clades C1 and C3, with the glacial (MIS 3–4) expansion and postglacial bottleneck preceding the recent expansion initiated in the late Holocene and were related to environmental changes during the upper Pleistocene and Holocene. The past demographic trends of the eastern Asian clade C3 were correlated with changes in sea level and the formation of new land bridges formed by the exposed sea shelf during the glaciations. These data were supported by reconstructed historical distribution models. The results of our genetic analyses, supported by the reconstruction of the historical spatial distributions of the distinct clades, confirm that over time the local populations mixed as a consequence of environmental and climatic changes resulting from cyclical glaciation and the interglacial period during the Pleistocene.

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.

Genetic Diversity of Ligidium Isopods in Hokkaido and Niigata, Northern Japan, Based on Mitochondrial DNA Analysis

The genetic diversity of the genus Ligidium in Hokkaido and Niigata, northern Japan, was investigated by analyzing the cytochrome c oxidase subunit 1 (CO1) region in the mitochondrial DNA (mtDNA). The genetic diversity in Hokkaido was much lower than that in Niigata. Nine different operational taxonomic units (OTUs) were identified. Only a single OTU, most likely Ligidium japonicum, was found in Hokkaido, whereas all nine OTUs were found in Niigata. Using the mtDNA evolutionary rate determined for the marine invertebrate Haptosquilla pulchella (Miers, 1880), population expansion for OTU1 in Hokkaido was estimated to have occurred at 12,600 years BP, suggesting that Ligidium underwent a bottleneck due to glacial cooling, and the population then expanded after postglacial warming. Assuming that the expansion of the OTU1 population occurred at 9600 years BP, when the sea surface temperature rose offshore of Tokachi in the Northwestern Pacific, the evolutionary rate (µ) of the mtDNA CO1 region in Ligidium is calculated as: 0.087 (95% confidence intervals: min: 0.042-max: 0.12) (substitutions/site/million years). The presence of a haplotype common to Hokkaido and Niigata implies that the haplotype migrated across the Tsugaru Strait. Considering that geological evidence indicates that the Tsugaru Strait was continuously present even during the last glacial maximum when the sea level was at its lowest, accidental transport by human beings or animals might have been critical to the migration of Ligidium.

House mouse Mus musculus dispersal in East Eurasia inferred from 98 newly determined complete mitochondrial genome sequences

The Eurasian house mouse Mus musculus is useful for tracing prehistorical human movement related to the spread of farming. We determined whole mitochondrial DNA (mtDNA) sequences (ca. 16,000 bp) of 98 wild-derived individuals of two subspecies, M. m. musculus (MUS) and M. m. castaneus (CAS). We revealed directional dispersals reaching as far as the Japanese Archipelago from their homelands. Our phylogenetic analysis indicated that the eastward movement of MUS was characterised by five step-wise regional extension events: (1) broad spatial expansion into eastern Europe and the western part of western China, (2) dispersal to the eastern part of western China, (3) dispersal to northern China, (4) dispersal to the Korean Peninsula and (5) colonisation and expansion in the Japanese Archipelago. These events were estimated to have occurred during the last 2000-18,000 years. The dispersal of CAS was characterised by three events: initial divergences (ca. 7000-9000 years ago) of haplogroups in northernmost China and the eastern coast of India, followed by two population expansion events that likely originated from the Yangtze River basin to broad areas of South and Southeast Asia, including Sri Lanka, Bangladesh and Indonesia (ca. 4000-6000 years ago) and to Yunnan, southern China and the Japanese Archipelago (ca. 2000-3500). This study provides a solid framework for the spatiotemporal movement of the human-associated organisms in Holocene Eastern Eurasia using whole mtDNA sequences, reliable evolutionary rates and accurate branching patterns. The information obtained here contributes to the analysis of a variety of animals and plants associated with prehistoric human migration.

Late Quaternary Environmental and Human Impacts on the Mitochondrial DNA Diversity of Four Commensal Rodents in Myanmar

We addressed the spatiotemporal characteristics of four commensal rodent species occurring in Myanmar in comparison with other areas of the Indo-Malayan region. We examined sequence variations of the mitochondrial cytochrome b gene (Cytb) in the Pacific rat (Rattus exulans), roof rat (Rattus rattus complex, RrC), lesser bandicoot rat (Bandicota bengalensis), and house mouse (Mus musculus) using the recently developed time-dependent evolutionary rates of mtDNA. The Cytb sequences of RrC from Myanmar were shown to belong to RrC Lineage II, and their level of genetic diversity was relatively high compared to those of the other three species. RrC was found to have experienced bottleneck and rapid expansion events at least twice in the late Pleistocene period in Myanmar and a nearby region. Accordingly, paleoclimatic environmental fluctuations were shown to be an important factor affecting rodents in the subtropics of the Indo-Malayan region. Our results show that human activities during the last 10,000 years of the Holocene period affected the population dynamics of the rodent species examined, including introducing them to Myanmar from neighboring countries. Further study of these four commensal rodents in other geographic areas of the Indo-Malayan region would allow us to better understand the factors that drove their evolution and their ecological trends.

Molecular Characterization of Species of the Subgenus Mus from Myanmar

The species composition and genetic diversity of the subgenus Mus in Myanmar are not yet fully understood. In this study, mice were trapped in rural areas near the Ayeyarwady River basin, spanning five Myanmar cities from north to south: Mandalay, Bagan, Magway, Pyay, and Yangon. Mitochondrial cytochrome b (Cytb) and nuclear melanocortin 1 receptor (Mc1r) gene sequences were determined for mice sampled and revealed a widespread occurrence of Mus nitidulus in central Myanmar in addition to its previously known occurrence in the Yangon district of southern Myanmar. Analyses of Cytb revealed two haplogroups with a genetic distance of 2.0%, suggestive of divergence several hundred thousand years ago. Mus caroli and M. musculus were confined to Yangon and Mandalay/Bagan/Magway, respectively. Mice collected from a locality on the eastern side of the Ayeyarwady River in Pyay were identified using Cytb and Mc1r sequences as M. fragilicauda, which was previously identified only in Laos and Thailand. The species M. booduga and M. cervicolor previously predicted to be common in the study area were not found. These findings elucidate the species and genetic diversity of the subgenus Mus in the Indo-Malayan Region.

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.

Contrasting phylogeographic histories between the continent and islands of East Asia: Massive mitochondrial introgression and long-term isolation of hares (Lagomorpha: Lepus)

Hares of the genus Lepus are distributed worldwide, and introgressive hybridization is thought to be pervasive among species, leading to reticulate evolution and taxonomic confusion. Here, we performed phylogeographic analyses of the following species of hare across East Asia: L. timidus, L. mandshuricus, L. coreanus, and L. brachyurus collected from far-eastern Russia, South Korea, and Japan. Nucleotide sequences of one mitochondrial DNA and eight nuclear gene loci were examined, adding sequences of hares in China from databases. All nuclear DNA analyses supported the clear separation of three phylogroups: L. timidus, L. brachyurus, and the L. mandshuricus complex containing L. coreanus. On the other hand, massive mitochondrial introgression from two L. timidus lineages to the L. mandshuricus complex was suggested in continental East Asia. The northern population of the L. mandshuricus complex was mainly associated with introgression from the continental lineage of L. timidus, possibly since the last glacial period, whereas the southern population of the L. mandshuricus complex experienced introgression from another L. timidus lineage related to the Hokkaido population, possibly before the last glacial period. In contrast to continental hares, no evidence of introgression was found in L. brachyurus in the Japanese Archipelago, which showed the oldest divergence amongst East Asian hare lineages. Our findings suggest that glacial-interglacial climate changes in the circum-Japan Sea region promoted distribution shifts and introgressive hybridization among continental hare species, while the geographic structure of the region contributed to long-term isolation of hares on the islands, preventing inter-species gene flow.

Mitochondrial Genetic Diversity and Phylogeography of Mus musculus castaneus in Northern Punjab, Pakistan

Regions of Iran, Afghanistan, Pakistan and northwestern India have been proposed as the place of origin of Mus musculus castaneus. But despite the fact that Pakistan encompasses an important part of its range, M. m. castaneus populations in Pakistan have not been the subject of intensive genetic and biogeographic studies, except for a very small number of samples included in past studies. We studied genetic variation in M. m. castaneus (CAS) from northern Punjab Province, Pakistan, by using cytochrome b (Cytb) analysis in a sample of 98 individuals. Median-joining network revealed four well differentiated CAS sub-lineages coexisting within a small geographical region; these had previously been thought to have largely non-overlapping geographic distributions. Moreover, haplotypes from Pakistan occupied a central position in the network and all identified global haplotypes were also present in Pakistan. All identified CAS sub-lineages proved to be highly diverse on the basis of haplotype and nucleotide diversity indices. Tajima's D test and Fu's Fs tests of neutrality suggest recent population expansions in all sub-lineages. Expansion times were estimated as 21,760-134,930, 10,800-64,400 and 4950-30,665 ybp using substitution rates of 2.5%, 5% and 11%, respectively. Our results support the hypothesis that northern Punjab Province in Pakistan is the most likely source area for M. m. castaneus, and that the CAS sub-lineages in this region have undergone rapid population expansion events at different time periods, which appear to have benefitted from human-mediated transport, although one of them clearly predates the establishment of human settlements in this region.