Phylogeography of Chinese house mice (Mus musculus musculus/castaneus): distribution, routes of colonization and geographic regions of hybridization

The effectiveness of selection in a species affects the direction of amino acid frequency evolution

Nearly neutral theory predicts that species with higher effective population size (N e ) are better able to purge slightly deleterious mutations. We compare evolution in high-N e vs. low-N e vertebrates to reveal which amino acid frequencies are subject to subtle selective preferences. We take three complementary approaches, two measuring flux and one measuring outcomes. First, we fit non-stationary substitution models of amino acid flux using maximum likelihood, comparing the high-N e clade of rodents and lagomorphs to its low-N e sister clade of primates and colugos. Second, we compare evolutionary outcomes across a wider range of vertebrates, via correlations between amino acid frequencies and N e . Third, we dissect the details of flux in human, chimpanzee, mouse, and rat, as scored by parsimony - this also enables comparison to a historical paper. All three methods agree on which amino acids are preferred under more effective selection. Preferred amino acids tend to be smaller, less costly to synthesize, and to promote intrinsic structural disorder. Parsimony-induced bias in the historical study produces an apparent reduction in structural disorder, perhaps driven by slightly deleterious substitutions. Within highly exchangeable pairs of amino acids, arginine is strongly preferred over lysine, and valine over isoleucine, consistent with more effective selection preferring a marginally larger free energy of folding. These two preferences match differences between thermophiles and mesophilic relatives. These results reveal the biophysical consequences of mutation-selection-drift balance, and demonstrate the utility of nearly neutral theory for understanding protein evolution.

Genealogical characterization of regional populations and dorsal coat color variation in the house mouse Mus musculus from Asia based on haplotype structure analysis of a gene-rich region harboring Mc1r

We analyzed 196 haplotype sequences from a gene-rich region (250 kb) that includes Mc1r, a gene involved in coat color regulation, to gain insight into the evolution of coat color variation in subspecies of the house mouse Mus musculus. Phylogenetic networks revealed haplotype groups from the major subspecies of M. m. castaneus (CAS), M. m. domesticus (DOM), and M. m. musculus (MUS). Using haplotype sequences assigned to each of CAS and MUS through phylogenetic analysis, we proposed migration routes associated with prehistoric humans from west to east across Eurasia. Comparing nucleotide diversity among subspecies-specific haplotypes in different geographic areas showed a marked reduction during migration, particularly in MUS-derived haplotypes from Korea and Japan, suggesting intensive population bottlenecks during migration. We found that a C>T polymorphism at site 302 (c.302C>T) in the Mc1r coding region correlated with a darkening of dorsal fur color in both CAS and MUS. However, C/C homozygous mice in MUS showed marked variation in lightness, indicating the possibility of another genetic determinant that affects the lightness of dorsal fur color. Detailed sequence comparisons of haplotypes revealed that short fragments assigned to DOM were embedded in CAS-assigned fragments, indicating ancient introgression between subspecies. The estimated age of c.302C>T also supports the hypothesis that genetic interaction between subspecies occurred in ancient times. This suggests that the genome of M. musculus evolved through gene flow between subspecies over an extended period before the movement of the species in conjunction with prehistoric humans.

Geographic isolation and climatic heterogeneity drive population differentiation of Rosa chinensis var. spontanea complex

Global biodiversity is contracting rapidly due to potent anthropogenic activities and severe climate change. Wild populations of Rosa chinensis var. spontanea and Rosa lucidissima are rare species endemic to China, as well as important germplasm resources for rose breeding. However, these populations are at acute risk of extinction and require urgent action to ensure their preservation. We harnessed 16 microsatellite loci to 44 populations of these species and analysed population structure and differentiation, demographic history, gene flow and barrier effect. In addition, a niche overlap test and potential distribution modelling in different time periods were also carried out. The data indicate that: (1) R. lucidissima cannot be regarded as a separate species from R. chinensis var. spontanea; (2) the Yangtze River and the Wujiang River function as barriers in population structure and differentiation, and precipitation in the coldest quarter may be the key factor for niche divergence of R. chinensis var. spontanea complex; (3) historical gene flow showed a converse tendency to current gene flow, indicating that alternate migration events of R. chinensis var. spontanea complex between south and north were a response to climate oscillations; and (4) extreme climate change will decrease the distribution range of R. chinensis var. spontanea complex, whereas the opposite will occur under a moderate scenario for the future. Our results resolve the relationship between R. chinensis var. spontanea and R. lucidissima, highlight the pivotal roles of geographic isolation and climate heterogeneity in their population differentiation, and provide an important reference for comparable conservation studies on other endangered species.

Comparative phylogeography of two commensal rat species (Rattus tanezumi and Rattus norvegicus) in China: Insights from mitochondrial DNA, microsatellite, and 2b-RAD data

Rattus norvegicus and Rattus tanezumi are dominant species of Chinese house rats, but the colonization and demographic history of two species in China have not been thoroughly explored. Phylogenetic analyses with mitochondrial DNA including 486 individuals from 31 localities revealed that R. norvegicus is widely distributed in China, R. tanezumi is mainly distributed in southern China with currently invading northward; northeast China was the natal region of R. norvegicus, while the spread of R. tanezumi in China most likely started from the southeast coast. A total of 123 individuals from 18 localities were subjected to 2b-RAD analyses. In neighbor-joining tree, individuals of R. tanezumi grouped into geographic-specific branches, and populations from southeast coast were ancestral groups, which confirmed the colonization route from southeast coast to central and western China. However, individuals of R. norvegicus were generally grouped into two clusters instead of geographic-specific branches. One cluster comprised inland populations, and another cluster included both southeast coast and inland populations, which indicated that spread history of R. norvegicus in China was complex; in addition to on-land colonization, shipping transportation also have played great roles. ADMIXTURE and principal component analyses provided further supports for the colonization history. Demographic analyses revealed that climate changes at ~40,000 to 18,000 years ago and ~4000 years ago had led to population declines of both species; the R. norvegicus declined rapidly while the population of R. tanezumi continuously expanded since ~1500 years ago, indicating the importance of interspecies' competition in their population size changes. Our study provided a valuable framework for further investigation on phylogeography of two species in China.

Genetic Diversity and Population Structure of Spirobolus bungii as Revealed by Mitochondrial DNA Sequences

Soil macrofauna, such as Spirobolus bungii, are an important component of ecosystems. However, systematic studies of the genetic diversity, population genetic structure, and the potential factors affecting the genetic differentiation of S. bungii are lacking. We performed a population genetic study of 166 individuals from the mountains to the south of the Yangtze River, north of the Yangtze River in Nanjing city, and near Tianjin city, in order to investigate the correlations between geographical distance and genetic diversity. A total of 1182 bp of COX2 and Cytb gene sequences of mitochondrial DNA, and 700 bp of the 18S rRNA gene sequence were analyzed. There were two haplotypes and one variable site in the 18S rRNA gene, and 28 haplotypes and 78 variable sites in the COX2 and Cytb genes. In this study, the 18S rRNA gene was used for species identification, and mtDNA (concatenated sequences with Cytb and COX2) was used for population genetic analysis. Structure cluster analysis indicated that the genetic structures of the different populations of S. bungii tended to be consistent at small geographical scales. Phylogenetic trees revealed that the haplotypes were clearly divided into three branches: the area south of the Yangtze River, the area to the north of the Yangtze River in Nanjing, and the area in Tianjin. Large geographical barriers and long geographical distance significantly blocked gene flow between populations of S. bungii. Our results provide a basic theoretical basis for subsequent studies of millipede taxonomy and population genetic evolution.

Genic and Chromosomal Components of Prdm9-Driven Hybrid Male Sterility in Mice (Mus musculus)

Hybrid sterility contributes to speciation by preventing gene flow between related taxa. Prdm9, the first and only hybrid male sterility (HMS) gene known in vertebrates, predetermines the sites of recombination between homologous chromosomes and their synapsis in early meiotic prophase. The asymmetric binding of PRDM9 to heterosubspecific homologs of Mus m. musculus x Mus m. domesticus F1 hybrids and increase of PRDM9-independent DNA double-strand break (DSB) hotspots results in difficult to repair DSBs, incomplete synapsis of homologous chromosomes and meiotic arrest at the first meiotic prophase. Here we show that Prdm9 behaves as a major HMS gene in mice outside the Mus m. musculus x Mus m. domesticus F1 hybrids, in the genomes composed of Mus m. castaneus and Mus m. musculus chromosomes segregating on the Mus m. domesticus background. The Prdm9cst/dom2 (castaneus/domesticus) allelic combination secures meiotic synapsis, testes weight and sperm count within physiological limits, while the Prdm9msc1/dom2 (musculus/domesticus) males show a range of fertility impairment. Out of five quantitative trait loci contributing to the Prdm9msc1/dom2-related infertility, four control either meiotic synapsis or fertility phenotypes and one controls both, synapsis and fertility. Whole-genome genotyping of individual chromosomes showed preferential involvement of nonrecombinant musculus chromosomes in asynapsis in accordance with the chromosomal character of HMS. Moreover, we show that the overall asynapsis rate can be estimated solely from the genotype of individual males by scoring the effect of nonrecombinant musculus chromosomes. Prdm9-controlled HMS represents an example of genetic architecture of HMS consisting of genic and chromosomal components.

Insights into Mus musculus population structure across Eurasia revealed by whole-genome analysis

For more than 100 years, house mice (Mus musculus) have been used as a key animal model in biomedical research. House mice are genetically diverse, yet their genetic background at the global level has not been fully understood. Previous studies have suggested that they originated in South Asia and diverged into three major subspecies, almost simultaneously, approximately 110,000–500,000 years ago; however, they have spread across the world with the migration of modern humans in prehistoric and historic times (∼10,000 years ago to the present day) and have undergone secondary contact, which has complicated the genetic landscape of wild house mice. In this study, we sequenced the whole-genome sequences of 98 wild house mice collected from Eurasia, particularly East Asia, Southeast Asia, and South Asia. Although wild house mice were found to consist of three major genetic groups corresponding to the three major subspecies, individuals representing admixtures between subspecies were more prevalent in East Asia than has been previously recognized. Furthermore, several samples exhibited an incongruent pattern of genealogies between mitochondrial and autosomal genomes. Using samples that likely retained the original genetic components of subspecies with the least admixture, we estimated the pattern and timing of divergence among the subspecies. The estimated divergence time of the three subspecies was 187,000–226,000 years ago. These results will help us to understand the genetic diversity of wild mice on a global scale, and the findings will be particularly useful in future biomedical and evolutionary studies involving laboratory mice established from such wild mice.

New Perspective on the Geographic Distribution and Evolution of Lymphocytic Choriomeningitis Virus, Central Europe

Lymphocytic choriomeningitis virus (LCMV) is an Old World mammarenavirus found worldwide because of its association with the house mouse. When LCMV spills over to immunocompetent humans, the virus can cause aseptic meningitis; in immunocompromised persons, systemic infection and death can occur. Central Europe is a strategic location for the study of LCMV evolutionary history and host specificity because of the presence of a hybrid zone (genetic barrier) between 2 house mouse subspecies, Mus musculus musculus and M. musculus domesticus. We report LCMV prevalence in natural mouse populations from a Czech Republic–Germany transect and genomic characterization of 2 new LCMV variants from the Czech Republic. We demonstrate that the main division in the LCMV phylogenetic tree corresponds to mouse host subspecies and, when the virus is found in human hosts, the mouse subspecies found at the spillover location. Therefore, LCMV strains infecting humans can be predicted by the genetic structure of house mice.

Demographic History and Genomic Response to Environmental Changes in a Rapid Radiation of Wild Rats

For organisms to survive and prosper in a harsh environment, particularly under rapid climate change, poses tremendous challenges. Recent studies have highlighted the continued loss of megafauna in terrestrial ecosystems and the subsequent surge of small mammals, such as rodents, bats, lagomorphs, and insectivores. However, the ecological partitioning of these animals will likely lead to large variation in their responses to environmental change. In the present study, we investigated the evolutionary history and genetic adaptations of white-bellied rats (Niviventer Marshall, 1976), which are widespread in the natural terrestrial ecosystems in Asia but also known as important zoonotic pathogen vectors and transmitters. The southeastern Qinghai-Tibet Plateau was inferred as the origin center of this genus, with parallel diversification in temperate and tropical niches. Demographic history analyses from mitochondrial and nuclear sequences of Niviventer demonstrated population size increases and range expansion for species in Southeast Asia, and habitat generalists elsewhere. Unexpectedly, population increases were seen in N. eha, which inhabits the highest elevation among Niviventer species. Genome scans of nuclear exons revealed that among the congeneric species, N. eha has the largest number of positively selected genes. Protein functions of these genes are mainly related to olfaction, taste, and tumor suppression. Extensive genetic modification presents a major strategy in response to global changes in these alpine species.

Distribution and Host Selection of Tropical Rat Mite, Ornithonyssus bacoti, in Yunnan Province of Southwest China

Simple Summary

The tropical rat mite (Ornithonyssus bacoti) is a transmission vector of rickettsia pox and a potential vector of hemorrhagic fever with renal syndrome (HFRS). This article reports the distribution and host selection of O. bacoti in Yunnan Province of Southwest China. The original data came from the investigations in 39 counties of Yunnan. The prevalence (P_M), mean abundance (MA) and mean intensity (MI) were calculated to reflect the infestations of the dominant rat hosts with O. bacoti mites. The patchiness index and Taylor’s power law were used to measure the spatial distribution of the mites. A total of 4121 O. bacoti mites were identified from 15 species of small mammal hosts in 27 of the 39 investigated counties, and 99.20% of them (4088/4121) were found on rodents. The majority of total O. bacoti mites was found in the flatland landscape (91.28%) and indoor habitat (73.48%). Moreover, 51.78% and 40.09% of O. bacoti mites were identified from Rattus tanezumi and R. norvegicus, the two synanthropic rat species. The mites had some host-specificity, with a preference to two dominant hosts (R. tanezumi and R. norvegicus), and they were of aggregated distribution on R. tanezumi.

Abstract

(1) Background: As a species of gamasid mite, the tropical rat mite (Ornithonyssus bacoti) is a common ectoparasite on rodents and some other small mammals. Besides stinging humans to cause dermatitis, O. bacoti can be a vector of rickettsia pox and a potential vector of hemorrhagic fever with renal syndrome (HFRS). (2) Objective: The present study was conducted to understand the host selection of O. bacoti on different animal hosts and the distribution in different environmental gradients in Yunnan Province of Southwest China. (3) Methods: The original data came from the investigations in 39 counties of Yunnan, between 1990 and 2015. The animal hosts, rodents and some other small mammals were mainly trapped with mouse traps. The O. bacoti mites on the body surface of animal hosts were collected and identified in a conventional way. The constituent ratio (Cr), prevalence (P_M), mean abundance (MA) and mean intensity (MI) were used to reflect infestations of animal hosts with O. bacoti mites. The patchiness index and Taylor’s power law were used to measure the spatial distribution pattern of O. bacoti mites on their hosts. (4) Results: A total of 4121 tropical rat mites (O. bacoti) were identified from 15 species and 14,739 individuals of hosts, and 99.20% of them were found on rodents. More than half of O. bacoti mites (51.78%) were identified from the Asian house rat (Rattus tanezumi), and 40.09% of the mites from the Norway rat (R. norvegicus) (p < 0.05). The infestations of R. tanezumi (P_M = 7.61%, MA = 0.40 and MI = 5.31) and R. norvegicus (P_M = 10.98, MA = 1.14 and MI = 10.39) with O. bacoti mites were significantly higher than those of other host species (p < 0.05). The infestations of two dominant rat hosts (R. tanezumi and R. norvegicus) with O. bacoti mites varied in different environmental gradients (latitudes, longitudes, altitudes, landscapes and habitats) and on different sexes and ages of the hosts. The prevalence of juvenile R. norvegicus rats with O. bacoti mites (P_M = 12.90%) was significantly higher than that of adult rats (P_M = 9.62%) (p < 0.05). The prevalence (P_M = 38.46%) and mean abundance (MA = 2.28 mites/host) of R. tanezumi rats with O. bacoti mites in the high latitude were higher than those in the low latitudes (p < 0.05). The majority of the total collected 4121 O. bacoti mites was found in the flatland landscape (91.28%) and indoor habitat (73.48%) (p < 0.05). The P_M (10.66%) and MA (0.49 mites/host) of R. tanezumi rats with O. bacoti mites were significantly higher in the indoor habitat than in the outdoor habitat (p < 0.05). The tropical rat mites showed an aggregated distribution pattern on their first dominant host, R. tanezumi. Conclusion: The tropical rat mite (O. bacoti) is a widely distributed species of gamasid mite in Yunnan Province, Southwest China, and its dominant hosts are two synanthropic species of rats, R. tanezumi and R. norvegicus. It is mainly distributed in the flatland landscape and indoor habitat. It has some host-specificity, with a preference to rodents, especially R. tanezumi and R. norvegicus. The O. bacoti mites are of aggregated distribution on R. tanezumi rats.

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.

Mus musculus populations in Western Australia lack VKORC1 mutations conferring resistance to first generation anticoagulant rodenticides: Implications for conservation and biosecurity

Background

Humans routinely attempt to manage pest rodent populations with anticoagulant rodenticides (ARs). We require information on resistance to ARs within rodent populations to have effective eradication programs that minimise exposure in non-target species. Mutations to the VKORC1 gene have been shown to confer resistance in rodents with high proportions of resistance in mice found in all European populations tested. We screened mutations in Mus musculus within Western Australia, by sampling populations from the capital city (Perth) and a remote island (Browse Island). These are the first Australian mouse populations screened for resistance using this method. Additionally, the mitochondrial D-loop of house mice was sequenced to explore population genetic structure, identify the origin of Western Australian mice, and to elucidate whether resistance was linked to certain haplotypes.

Results

No resistance-related VKORC1 mutations were detected in either house mouse population. A genetic introgression in the intronic sequence of the VKORC1 gene of Browse Island house mouse was detected which is thought to have originated through hybridisation with the Algerian mouse (Mus spretus). Analysis of the mitochondrial D-loop reported two haplotypes in the house mouse population of Perth, and two haplotypes in the population of Browse Island.

Conclusions

Both house mouse populations exhibited no genetic resistance to ARs, in spite of free use of ARs in Western Australia. Therefore weaker anticoagulant rodenticides can be employed in pest control and eradication attempts, which will result in reduced negative impacts on non-target species. Biosecurity measures must be in place to avoid introduction of resistant house mice, and new house mouse subspecies to Western Australia.

Further resolution of the house mouse (Mus musculus) phylogeny by integration over isolation-with-migration histories

BACKGROUND

The three main subspecies of house mice, Mus musculus castaneus, Mus musculus domesticus, and Mus musculus musculus, are estimated to have diverged ~ 350-500KYA. Resolution of the details of their evolutionary history is complicated by their relatively recent divergence, ongoing gene flow among the subspecies, and complex demographic histories. Previous studies have been limited to some extent by the number of loci surveyed and/or by the scope of the method used. Here, we apply a method (IMa3) that provides an estimate of a population phylogeny while allowing for complex histories of gene exchange.

RESULTS

Results strongly support a topology with M. m. domesticus as sister to M. m. castaneus and M. m. musculus. In addition, we find evidence of gene flow between all pairs of subspecies, but that gene flow is most restricted from M. m. musculus into M. m. domesticus. Estimates of other key parameters are dependent on assumptions regarding generation time and mutation rate in house mice. Nevertheless, our results support previous findings that the effective population size, Ne, of M. m. castaneus is larger than that of the other two subspecies, that the three subspecies began diverging ~ 130 - 420KYA, and that the time between divergence events was short.

CONCLUSIONS

Joint demographic and phylogenetic analyses of genomic data provide a clearer picture of the history of divergence in house mice.

Colonization and Differentiation Traits of the Japanese House Mouse, Mus musculus (Rodentia, Muridae), Inferred from Mitochondrial Haplotypes and External Body Characteristics

To evaluate the colonization histories of the Japanese house mice (Mus musculus), phenotypic and genotypic admixtures of the subspecific traits were studied by evaluation of external body characteristics and mitochondrial gene elements. We analyzed mitochondrial Cytb gene and coat colorations and body dimensions as subspecific characteristics in mice from four areas of the Japanese Islands, the Sorachi, Ishikari and Iburi areas of Hokkaido, the Hidaka area of Hokkaido, and northeastern and central Honshu. Three occurrence patterns of the subspecific haplotypes of Cytb-the castaneus type only, the musculus type only, and the castaneus, musculus, and domesticus types together-were observed in the study areas. In central Honshu, the properties of haplotypes were in accord with the external characteristics as reported in previous findings. In contrast, complicated external characteristics were observed in the Hidaka area, where mice showed multiple haplotype properties. In addition, in northeastern Honshu, coat colorations were not in accord with haplotype properties and such discordance was also observed in most mice in the Sorachi, Ishikari and Iburi areas of Hokkaido. These complexities and discordances suggest that the genetic and phenotypic properties have been caused by different processes, not only through founder effects by migrations and subsequent subspecific hybridizations but also through differentiation in each study area.

Disjunct distribution and distinct intraspecific diversification of Eothenomys melanogaster in South China

BACKGROUND

South China encompasses complex and diverse landforms, giving rise to high biological diversity and endemism from the Hengduan Mountains to Taiwan Island. Many species are widely distributed across South China with similar disjunct distribution patterns. To explore the causes of these disjunct distribution patterns and their genetic consequences, we investigated the endemic species Père David's Chinese Vole (Eothenomys melanogaster) by integrating geological and ecological factors. We analysed the genetic structure and divergence time of E. melanogaster based on fast-evolving mitochondrial and nuclear markers using Bayesian trees and coalescent species tree approaches. Historical scenarios of distribution range and demography were reconstructed based on spatial interpolations of genetic diversity and distance, extended Bayesian skyline plots, phylogeographic diffusion analysis, and ecological niche modelling (ENM) during different periods. We also assessed the relationships between geographical distance/ecological vicariance and genetic distance (isolation by distance, IBD; isolation by environment, IBE).

RESULTS

The genetic analysis revealed three deeply divergent clades-Southeast, Southwest and Central clades, centred on the Wuyi Mountains, the Yunnan-Guizhou Plateau (YGP) and the mountains around the Sichuan Basin, respectively-that have mostly developed since the Pleistocene. IBD played an important role in early divergence, and geological events (sedimentation of plains and linking of palaeo-rivers) and IBE further reinforced genetic differentiation. ENM shows the importance of suitable habitats and elevations.

CONCLUSIONS

Our results suggest that the primary cause of the disjunct distribution in E. melanogaster is the high dependence on middle-high-altitude habitat in the current period. Mountains in the occurence range have served as "sky islands" for E. melanogaster and hindered gene flow. Pleistocene climatic cycles facilitated genetic admixture in cold periods and genetic diversification in warm periods for inland clades. During cold periods, these cycles led to multiple colonization events between the mainland and Taiwan and erased genetic differentiation.

Phylogeography of the Rice Spittle Bug (Callitettix versicolor) Implies Two Long-Term Mountain Barriers in South China

South China is a region of remarkable topographic complexity. However, the impact of climate fluctuations in the Pleistocene on the local fauna and especially insects has not been extensively studied. We integrated mitochondrial DNA (mtDNA) and microsatellite data of the rice spittle bug, Callitettix versicolor, to determine the genetic structure, potential biogeographic barriers, and historical demography of this species. The mtDNA data revealed two distinct lineages (Western and Eastern) congruent with the geographically separated western and eastern sub-regions of the Hengduan Mountains. The Eastern lineage was subdivided into two sub-lineages, E1 and E2, congruent with the geographically separated northern and southern sub-regions of the Dabie Mountains. E2 was further subdivided into two sub-groups, E2-1 and E2-2, with a hybrid zone (Guizhou and Hunan Provinces) in which their areas were contiguous. The genetic structures constructed using mtDNA were corroborated by four clusters (G1-G4) of microsatellite data. The populations of each cluster were nearly consistent with a sub-lineage of the mtDNA gene tree (G1-G4 corresponded to the Western, E1, E2-2 and E2-1 lineages, respectively). The divergence time estimated between the Western and Eastern lineages was 1.17 (0.50-2.37) to 0.89 (0.39-1.78) Mya, indicating that the lineages diversified on both geographic and temporal scales. The historical demography of the Eastern lineage showed continuous population growth after the Last Interglacial (LIG) and a stable population during the Last Glacial Maximum (LGM) period. However, the Western lineage remained largely unchanged during the LIG and LGM periods. This suggests that the historical demography of C. versicolor is probably related not only to the paleoclimate of South China, but also to the geological restriction and specific habitat preferences of species.

Geographical distribution of Toxoplasma gondii genotypes in Asia: A link with neighboring continents

Defining the pattern of genetic diversity of Toxoplasma gondii is important to understand its worldwide distribution. During the last decades, a large number of studies have been published on Toxoplasma genotypes circulating in Europe, in North and South America. Two continents are still largely unexplored, Africa and, to a less extent, Asia. In this last continent, an increasing number of publications reported genotypes circulating in diverse provinces of China, but very few data are available for other Asian countries. After a systematic database search, 47 papers related to T. gondii genotypes in Asia were analyzed. Genetic characterization of DNA was performed by microsatellite markers, or more usually by a multiplex PCR using 11 PCR-RFLP markers, allowing data comparison to draw a first global picture of the population structure of this parasite throughout Asia. Overall, 390 isolates or DNA extracts were completely typed by PCR-RFLP and/or microsatellite marker methods, revealing 36 different PCR-RFLP or equivalent microsatellite genotypes: 15 genotypes identified by a ToxoDB number and 21 atypical or unique genotypes. The most common genotype found in Asia is the genotype ToxoDB#9 (Chinese 1). The clonal types I, II and II variant, and III were also commonly found in Asia. The geographical distribution of these genotypes across Asia may reflect either a continuum with Europe for the western part of Asia (presence of Type II), or the circulation of strains through animal migration or human activities between Africa and the Southwestern part of Asia (Africa 1 genotype in Turkey or ToxoDB#20 both I Sri-Lanka and in Ethiopia or Egypt). Although there are some indications of a genetic population structure in Southeast Asian countries different from the rest of Asia, more studies in this tropical part of Asia will be necessary for a region which represent as well as Africa one of the missing links of the T. gondii genetic diversity.

First molecular evidence of Mus musculus bactrianus in Nepal inferred from the mitochondrial DNA cytochrome B gene sequences

To identify the house mice collected in Pokhara and Lumbini of Nepal at the subspecies level, morphological and molecular analyses were carried out. Morphologically, two populations collected in Pokhara and Lumbini were distinguished by fur colour, but there was no significant difference in external measurements (p > .05). The phylogenetic analysis results revealed that the haplotypes sequences of mitochondrial DNA (mtDNA) Cytochrome B (CytB) gene distinguished into two distinct clades on a phylogenetic tree representing two subspecies, Mus musculus bactrianus and M. m. castaneus in Pokhara and Lumbini, respectively. In Nepal, the subspecies M. m. bactrianus was not reported before this study. These findings concluded that at least two subspecies, M. m. bactrianus and M. m. castaneus currently exist in Nepal. We estimated that these two subspecies could have introduced together with human migration, while further study is required to understand their evolutionary history and current distribution.

Yangtze River, an insignificant genetic boundary in tufted deer (Elaphodus cephalophus): the evidence from a first population genetics study

Great rivers were generally looked at as the geographical barrier to gene flow for many taxonomic groups. The Yangtze River is the third largest river in the world, and flows across South China and into the East China Sea. Up until now, few studies have been carried out to evaluate its effect as a geographical barrier. In this study, we attempted to determine the barrier effect of the Yangtze River on the tufted deer (Elaphodus cephalophus) using the molecular ecology approach. Using mitochondrial DNA control region (CR) sequences and 13 nuclear microsatellite loci, we explored the genetic structure and gene flow in two adjacent tufted deer populations (Dabashan and Wulingshan populations), which are separated by the Yangtze River. Results indicated that there are high genetic diversity levels in the two populations, but no distinguishable haplotype group or potential genetic cluster was detected which corresponded to specific geographical population. At the same time, high gene flow was observed between Wulingshan and Dabashan populations. The tufted deer populations experienced population decrease from 0.3 to 0.09 Ma BP, then followed by a distinct population increase. A strong signal of recent population decline (T = 4,396 years) was detected in the Wulingshan population by a Markov-Switching Vector Autoregressions(MSVAR) process population demography analysis. The results indicated that the Yangtze River may not act as an effective barrier to gene flow in the tufted deer. Finally, we surmised that the population demography of the tufted deer was likely affected by Pleistocene climate fluctuations and ancient human activities.

Genomic resources for wild populations of the house mouse, Mus musculus and its close relative Mus spretus

Wild populations of the house mouse (Mus musculus) represent the raw genetic material for the classical inbred strains in biomedical research and are a major model system for evolutionary biology. We provide whole genome sequencing data of individuals representing natural populations of M. m. domesticus (24 individuals from 3 populations), M. m. helgolandicus (3 individuals), M. m. musculus (22 individuals from 3 populations) and M. spretus (8 individuals from one population). We use a single pipeline to map and call variants for these individuals and also include 10 additional individuals of M. m. castaneus for which genomic data are publically available. In addition, RNAseq data were obtained from 10 tissues of up to eight adult individuals from each of the three M. m. domesticus populations for which genomic data were collected. Data and analyses are presented via tracks viewable in the UCSC or IGV genome browsers. We also provide information on available outbred stocks and instructions on how to keep them in the laboratory.

The Pace of Hybrid Incompatibility Evolution in House Mice

Hybrids between species are often sterile or inviable. This form of reproductive isolation is thought to evolve via the accumulation of mutations that interact to reduce fitness when combined in hybrids. Mathematical formulations of this "Dobzhansky-Muller model" predict an accelerating buildup of hybrid incompatibilities with divergence time (the "snowball effect"). Although the Dobzhansky-Muller model is widely accepted, the snowball effect has only been tested in two species groups. We evaluated evidence for the snowball effect in the evolution of hybrid male sterility among subspecies of house mice, a recently diverged group that shows partial reproductive isolation. We compared the history of subspecies divergence with patterns of quantitative trait loci (QTL) detected in F2 intercrosses between two pairs of subspecies (Mus musculus domesticus with M. m. musculus and M. m. domesticus with M. m. castaneus). We used a recently developed phylogenetic comparative method to statistically measure the fit of these data to the snowball prediction. To apply this method, QTL were partitioned as either shared or unshared in the two crosses. A heuristic partitioning based on the overlap of QTL confidence intervals produced unambiguous support for the snowball effect. An alternative approach combining data among crosses favored the snowball effect for the autosomes, but a linear accumulation of incompatibilities for the X chromosome. Reasoning that the X chromosome analyses are complicated by low mapping resolution, we conclude that hybrid male sterility loci have snowballed in house mice. Our study illustrates the power of comparative genetic mapping for understanding mechanisms of speciation.