Eurasian house mouse (Mus musculus L.) differentiation at microsatellite loci identifies the Iranian plateau as a phylogeographic hotspot
© Hardouin et al.; licensee BioMed Central. 2015
Received: 29 September 2014
Accepted: 12 February 2015
Published: 25 February 2015
The phylogeography of the house mouse (Mus musculus L.), an emblematic species for genetic and biomedical studies, is only partly understood, essentially because of a sampling bias towards its most peripheral populations in Europe, Asia and the Americas. Moreover, the present-day phylogeographic hypotheses stem mostly from the study of mitochondrial lineages. In this article, we complement the mtDNA studies with a comprehensive survey of nuclear markers (19 microsatellite loci) typed in 963 individuals from 47 population samples, with an emphasis on the putative Middle-Eastern centre of dispersal of the species.
Based on correspondence analysis, distance and allele-sharing trees, we find a good coherence between geographical origin and genetic make-up of the populations. We thus confirm the clear distinction of the three best described peripheral subspecies, M. m. musculus, M. m. domesticus and M. m. castaneus. A large diversity was found in the Iranian populations, which have had an unclear taxonomic status to date. In addition to samples with clear affiliation to M. m. musculus and M. m. domesticus, we find two genetic groups in Central and South East Iran, which are as distinct from each other as they are from the south-east Asian M. m. castaneus. These groups were previously also found to harbor distinct mitochondrial haplotypes.
We propose that the Iranian plateau is home to two more taxonomic units displaying complex primary and secondary relationships with their long recognized neighbours. This central region emerges as the area with the highest known diversity of mouse lineages within a restricted geographical area, designating it as the focal place to study the mechanisms of speciation and diversification of this species.
KeywordsHouse mouse Microsatellites Phylogeographic hotspot Iranian plateau
The house mouse (Mus musculus L.) has long been viewed has an excellent model for the study of evolution and its genome has been one of the first to be nearly completely sequenced [1,2]. Moreover, its dispersal capacity through commensalism has ranked it as one of the “100 world worst most invasive alien species” (ISSG), therefore offering various possibilities to study adaptation to various environments . At the same time, this species is one of the most studied vertebrates due to its use as a prominent laboratory model, but its phylogeography and population genetics is so far only partly understood . The current knowledge has slowly accumulated over the last 30 years in a non-optimal fashion, since its most peripheral populations in Europe, Asia and the Americas have been studied before insights were gained for those from the Middle-Eastern centre of its distribution .
It is now widely recognised that Mus musculus L. constitutes a complex assembly of more or less well separated populations and subspecies. The term “subspecies” in itself is taken here in its broad sense of “genetically recognisable entities” but this does not imply on our part any deeper statement about the actual level of isolation among these entities. The last 45 years of literature on systematics of the house mouse revealed that nomenclatorial issues have been quite controversial, with the use of many terms ranging from biochemical groups, subspecies, semi-species to full species to designate the same entities. Here, we follow the generally held view that the more widely distributed populations are grouped into three different subspecies: Mus musculus musculus in Eastern Europe, Central and North East Asia, Mus musculus domesticus in Northern Africa and Western Europe, and Mus musculus castaneus in South East Asia. These last two subspecies have further expanded in modern times to the Americas, Australia and Oceania [6-9]. In addition, Mus musculus molossinus, a hybrid between M. m. musculus and M. m. castaneus found in Japan  is often considered as a subspecies on its own. Closer to the centre of the distribution, Mus musculus gentilulus has been identified in the eastern part of the Arabic peninsula on the basis of its mitochondrial DNA lineage  while from the same type of data [12,13] it has been shown that certain populations considered as M. m. castaneus in Iran, Pakistan and Afghanistan should probably be considered as belonging to further sub-specific groups. Moreover, another completely independent lineage has recently been identified on this basis in Nepal . Hence, the taxonomic situation close to the Middle-Eastern centre is far from being fully clarified. Since taxonomy reflects history, this clarification is a prerequisite if we want to further study the evolutionary mechanisms accounting for the species’ differentiation.
Our main findings support the global picture of a species having initially differentiated in several genetic entities in the Middle-East forming the extant subspecies. Several of these subsequently expanded outwards because of their propensity to engage in commensalism and finally colonise the entire world. Of particular interest is the situation of the populations inhabiting the Iranian plateau where, despite an important level of secondary admixture in this region, four main genetic groups could be identified, in partial congruence with the mitochondrial analysis.
Geographic origin of samples and their genetic diversity parameters at 19 microsatellite loci
Average number of allele
“South East Iran”
“South East Iran”
“South East Iran”
Iran North East
Individual allele-sharing distance tree
Summary of STRUCTURE results as displayed with STRUCTURE Harvester
As was to be expected, the population samples from the best described subspecies differentiate clearly from each other, and can be unambiguously assigned to M. m. musculus, M. m. domesticus or M. m. castaneus in keeping with previous studies (see  for a review). These long recognised entities are separated by the three first canonical axes of the CA in Figure 2, and constitute the three longest branches of the distance tree in Figure 3 and occupy different sectors of the circle in the allele-sharing tree of Figure 4. For the sake of comparison with the literature, an AMOVA performed among these three entities only indicated a FCT value of 0.70 for an average within group differentiation of FSC = 0.12.
When calculated inside each of our sampling of the peripheral subspecies, this becomes 0.10, 0.13 and 0.09 for M. m. castaneus, M. m. domesticus and M. m. musculus respectively, thus indicating a consequent amount of internal variation inside each of them. However, when the whole collection was submitted to the Bayesian assignation of STRUCTURE, only the partition between M. m. domesticus and all the rest appeared stable (Figure 5A).
One noticeable point is that there is a quite good geographical coherence in each group. If one considers M. m. domesticus for instance, the fan-shaped sub-tree of Figure 3 separates well the Near-East samples from the European ones. The Senegalese sample appears clearly as a recent offshoot of European origin, while this is not the case for Northern-African ones (Morocco and Tunisia). For M. m. musculus, the feather-like sub-tree indicate remote branching for the Chinese and Kazak samples, a fact to be put in relation with the probable recent eastwards expansion of this subspecies from the Caucasus region [26,30]. As to the M. m. castaneus cluster, it is the one where unexpected assignations occur. As reported above, the Malagasy sample clearly belongs to this group despite it possessing M. m. gentilulus mtDNA. This is another example of the possibility for maternal lineage capture during an expansion process. Madagascar is thought to have been populated by mice quite recently (ca 1,000 years, ) with the development of trade between India and Indonesia (were the nuclear genome most likely comes from), the Arabic peninsula (for the M. m. gentilulus maternal lineage ) and Africa. Of interest also is the position of the Kenyan mice, both from the coast (Mombasa) and the interior (Nairobi) which appear also as a recent M. m. castaneus offshoot, despite harbouring a mixture of M. m. castaneus and M. m. domesticus matrilines  as well as supplementary matrilines (see Additional file 2). Another noticeable point is, as expected, that samples known to be introgressed as the above-mentioned Abkhazian one, tend to “move” toward the centre of the tree. This is also the case for the Armenian sample from Megri, which is M. m. musculus for 2/3 of its nuclear genome  and was already shown to be under multiple genetic influences since it also contains for instance at the ABP locus three different alleles, each one ascribed to a given subspecies .
Status of Iranian plateau populations
The most interesting results of this study concern the samples of the Iranian plateau, which have sometimes been coined “central populations” because of many uncertainties as to their taxonomic affiliation. Putting aside the samples unambiguously referable to M. m. musculus (Khakh-Qaene and Iran North-East) and M. m. domesticus (Ahvaz and Hamedan), we are left with two clearly separated entities as described in the Results section. Interestingly, this fits rather well with the mitochondrial description recently provided  where each of the entities recognised in Iran is associated predominantly with one haplogroup (HG) essentially not found elsewhere. Namely, when recompiling the aforementioned data  for the populations’ samples of this study, the Central Iran group is associated with mitochondrial haplogroup Hg1B (78%) and South-East Iranian is associated with haplogroup Hg3 at 88%, while these two HGs are practically non-existent outside Iran where M. m. castaneus is almost completely associated with Hg2. The mitochondrial tree of all studied populations is provided in Additional file 1: Figure S3 to illustrate this point. As discussed by Rajabi-Maham et al. , the separate coalescence of these phylogenetically independent haplogroups have most likely arisen in allopatry during past periods of geographical isolation, and our nuclear genome results fit quite well with this view as they also point towards the existence of two independent groups in Iran, one in the center, and one in the South-East. These two groups are loosely related to M. m. castaneus, as are their matrilines (Additional file 1: Figure S3). Furthermore, the separation of the different entities in Iran seems to fit the topography of the country. M. m. domesticus samples are separated from the Central Iran individuals by the Zagros Mountains, Central Iran from M. m. musculus through the Kavir desert and from south-east Iran through the Lut desert (Figure 1B).
The central and south-east Iranian phylogroups are as distinct from each other as they are from the south-east Asian M. m. castaneus, even if many signs of admixture and secondary contacts are to be found, be it on the nuclear genome or the mitochondrial data. This is likely the principal reason why the STRUCTURE clustering did not converge to stable partitions for values of K larger than 2; probably because certain loci may have been more prone to be exchanged than other after secondary exchange. Hence, the present genomic make-up of Iranian populations is likely to be a mosaic of locally evolved haplotypes (like the mitochondrial matrilines) and segments imported from their neighbours in proportions that remain to be estimated. Nevertheless, they necessarily have evolved as geographical isolates during a certain amount of time and cannot be considered as resulting from simple admixture of surrounding populations. Hence, we may consider them as new independent entities, even if monophyly will rarely be achieved because of the retention of ancestral polymorphisms and secondary gene flow.
Attempts at setting the time frame for these exchanges with ABC methods have been recently performed , but these did not formally consider the existence of independent Iranian entities. The ancient origin of most Iranian phylogroups has however been recently reinforced by independent data showing they may have retained ancestral morphological features . Only more sophisticated model-based analyses relying on genome-wide multilocus data will be able to tell what kind of genetic exchanges they are still able to maintain. Such genome-wide data exist currently solely for the comparison between the peripheral M. m. domesticus populations in Germany and France, versus M. m. musculus in Czech Republic and Kazakhstan  and between wild-derived lines representing each of the three peripheral subspecies . These studies have revealed complex patterns of mutual introgression but cannot enlighten us as to the history and status of the Iranian phylogroups.
From the taxonomical standpoint, this situation deserves further scrutiny, since the distinctiveness of these three phylogroups should prevent the use of a single Latin trinomial like M. m. castaneus. However, solving the problem would require having access to hypothetical voucher specimens and a better delineation of present day limits of these phylogroups. This is not an easy task, because the species as a whole has undergone a phase of post-glacial expansion triggered by its association with Neolithic humans which may have induced a phase of complex hybridization, introgression and de-differentiation. This de-differentiation reaches some clear limits when distantly related subspecies like M. m. domesticus and M. m. musculus are concerned since their contact amounts to hybrid zones with restricted gene flow (see for instance  and the literature cited therein), but is not so clear what happens between the four non-domesticus entities considered here.
Our analysis of a Eurasian collection of wild mice encompassing both peripheral populations and a large number of central ones is the first comprehensive approach aiming at providing a global view of the relationships among the taxa constitutive of this complex subspecific assemblage. We propose here that the Iranian plateau is home to two more units displaying complex primary and secondary relationships with their long recognized neighbours M. m. domesticus, M. m. musculus and M. m. castaneus. If one adds to this the M. m. gentilulus from the Arabic peninsula and the new clade described from Nepal , the taxonomic diversity is clearly higher than previously thought. This is not a surprise if one considers the highly tormented landscape of the central regions where the complex species originated from, where there was clearly space for more than three geographical isolates before its worldwide expansion. The fact that this happens in and around the Iranian plateau is consistent with the key geographical position of this region which constitute an obligate passage between the lowlands of central Asia to the North, the Fertile Crescent and the Near-East to the West, and the Indian subcontinent to the East. Although the present-day distribution of distinct phylogroups does not necessarily indicate where the ancestral species was originally located, it becomes clear from the above results that the Iranian plateau should be included in the candidate regions together with the Indo-Pak sub-continent and neighboring Afghanistan. Indeed, in this large Middle-Eastern region, no less than six distinct phylogroups have differentiated and presently interact with each other. This should now be taken into account when building evolutionary scenarios.
Most of the individuals included in this study come from the DNA collection established at ISE-M between 1985 and 2005 with the exception of the samples from Massif Central, Cologne-Bonn, Czech Republic and Kazakhstan that come from the collection held at the MPI-Plön . All these samples have been described and used in previous publications, the references of which are given in Table 1, together with their geographical locations.
19 microsatellites were chosen from previous studies: PP8E11 (Chr2), PP6E09 (Chr3), D1EnsmusG22992 (Chr3),PP10E08 (Chr4), D5Mit149 (Chr5), D6Mit309 (Chr6), PP10A02 (Chr8), D9Mit54 (Chr9), D9Mit330 (Chr9), PP3A02 (Chr10), PP4A02 (Chr11), D13Mit61 (Chr13), D14Mit203 (Chr14), E1EnsmusG46849 (Chr15), D15Mit98 (Chr15), PP7B08 (Chr16), PP8A05 (Chr17), D19Mit39 (Chr19) and, PP2A01 (ChrX) [37,38]. Forward primers were labeled with FAM or HEX dye on the 5’end and the reaction conditions were as follow: denaturation, 95°C for 15 min, 28 cycles with 0.30s at 95°C, 1.30 min at 60°C and 1.30 min at 72°C, final elongation 10 min at 72°C. Primers (Additional file 1: Table S1) were pooled in 5 different reactions with a final reaction volume of 5 μl using Multiplex PCR kit (Qiagen). PCR products were then diluted 1/20 in Millipore water. 1 μL of this dilution was added to a mix of 0.1 μL ROX Standard (Applied Biosystems) and 10 μL HiDi Formamid. Heating for the denaturation step was 2 min at 90°C followed by 5 min at 20°C. Microsatellites were scored using GeneMapper (Applied Bioscience).
Basic diversity parameters were obtained using Genetix  and Arlequin . Several ways of dissecting the genetic variability present at the 19 microsatellites loci were used concurrently. We first performed a Correspondence Analysis (CA) using the AFC-3D procedure of Genetix. This is an unsupervised method which allows representation of the differentiation of samples along independent factorial axes. We did this by taking the centroids of samples as active elements (and individuals as supplementary elements) considering alternatively the whole collection (47 population samples, 963 individuals) with those samples previously described unambiguously as belonging to the peripheral subspecies, either M. m. domesticus (22 populations, 367 individuals), M. m. musculus (11 populations, 172 individuals) or M. m. castaneus (7 populations, 133 individuals) considered individually, or the samples from Iran (11 populations, 291 individuals) including individuals from all the different sub-species.
The pairwise Reynold’s distances were computed using the Gendist program of the Phylip package  and 1,000 bootstrapped datasets were generated with Seqboot, followed by the Neighbor and Consense procedures from the same package to obtain a Neighbour-Joining consensus population tree. Additionally, an Allele Sharing distance tree considering all individuals was calculated using the software Populations 1.2.32  and visualized with MEGA 6 .
Finally, we used the STRUCTURE program  to assign the collection of individuals to several putatively reproductively independent groups. The run parameters were as follow: a burn-in period of 1,000,000 simulations followed by a run length of 1,000,000 MCMC simulations and ten iterations for K (number of clusters) equals 1 to 5. The runs were performed using the admixture model with the Loc Prior option. Results were summarized using Structure Harvester . K was chosen using the criterion of Evanno et al. . To draw the structure diagrams the softwares CLUMPP  and Distruct  were used.
The mitochondrial D-loop sequences tree from Additional file 1: Figure S3 was produced with MEGA6 . The inference was performed according to the maximum likelihood method based on the Tamura-Nei model . The boostrap values (150 replicates) are shown. Initial tree(s) for the heuristic search were obtained by applying the Neighbor-Joining method to a matrix of pairwise distances estimated using the Maximum Composite Likelihood (MCL) approach. A discrete Gamma distribution was used to model evolutionary rate differences among sites (5 categories (+G, parameter = 0.1440)). The rate variation model allowed for some sites to be evolutionarily invariable ([+I], 62.0668% sites). The tree is drawn to scale, with branch lengths measured in the number of substitutions per site. The analysis involved 867 nucleotide sequences. All positions with less than 95% site coverage were eliminated. That is, fewer than 5% alignment gaps, missing data, and ambiguous bases were allowed at any position. There were a total of 849 positions in the final dataset.
Availability of supporting data
The microsatellites dataset supporting the results of this article on Dryad: doi:10.5061/dryad.ck276.
We are grateful to the numerous colleagues and students that have helped us to collect these samples over more than 30 years, with a special mention of Hassan Rajabi-Maham and Roohollah Siahsarvie for their help with Iranian mice. Special thanks goes to Bettina Harr for having designed and run the first microsatellite typing, to Raphael Leblois and Mathieu Gautier for having spent a considerable amount of time to convince us that STRUCTURE and GENELAND did not converge on our data set, and to Heinke Buhtz and Conny Burghardt for excellent technical assistance. John Stewart helped improve the English language. We are also grateful to the two anonymous reviewers and the editor for helping to improve the manuscript. This is publication ISEM 2015-021.
- Mouse Genome Sequencing Consortium. Initial sequencing and comparative analysis of the mouse genome. Nature. 2002;420:520–62.View ArticleGoogle Scholar
- Church DM, Goodstadt L, Hillier LW, Zody MC, Goldstein S, She X, et al. Lineage-specific biology revealed by a finished genome assembly of the mouse. PLoS Biol. 2009;7(5):e1000112. doi:10.1371/journal.pbio.1000112.View ArticlePubMed CentralPubMedGoogle Scholar
- Teschke M, Buentge A, Tautz D. Tracing recent adaptations in natural populations of the house mouse. In: Macholán M, Baird SJE, Munclinger P, Piálek J, editors. Evolution of the House mouse (Cambridge Studies in Morphology and Molecules: New Paradigms in Evolutionary Biology). UK: Cambridge University Press; 2012. p. 315–33.View ArticleGoogle Scholar
- Guénet J-L, Bonhomme F. Wild mice: an ever-increasing contribution to a popular mammalian model. Trends Genet. 2003;19:24–31.View ArticlePubMedGoogle Scholar
- Bonhomme F, Searle JB. House mouse phylogeography. In: Macholán M, Baird SJE, Munclinger P, Piálek J, editors. Evolution of the House mouse (Cambridge Studies in Morphology and Molecules: New Paradigms in Evolutionary Biology). UK: Cambridge University Press; 2012. p. 278–96.View ArticleGoogle Scholar
- Searle JB, Jamieson PM, Gündüz I, Stevens MI, Jones EP, Gemmill CEC, et al. The diverse origins of New Zealand house mice. Proc R Soc B. 2009;276:209–17. doi:10.1098/rspb.2008.0959.View ArticlePubMed CentralPubMedGoogle Scholar
- Hardouin EA, Chapuis J-L, Stevens MI, Van Vuuren JB, Quillfeldt P, Scavetta RJ, et al. House mouse colonization patterns on the sub-Antarctic Kerguelen Archipelago suggest singular primary invasions and resilience against re-invasion. BMC Evol Biol. 2010;10:325.View ArticlePubMed CentralPubMedGoogle Scholar
- Gabriel SI, Stevens MI, Da Luz MM, Searle JB. Of mice and ‘convicts’: origin of the Australian house mouse. Mus musculus. PLoS One. 2011;6(12):e28622. doi:10.1371/journal.pone.0028622.View ArticlePubMed CentralPubMedGoogle Scholar
- Jones EP, Jóhannesdóttir F, Gündüz I, Richards MB, Searle JB. The expansion of the house mouse into northern-western Europe. J Zool. 2011;283:257–68.View ArticleGoogle Scholar
- Suzuki H, Aplin KP. Phylogeny and biogeography of the genus Mus in Eurasia. In: Macholán M, Baird SJE, Munclinger P, Piálek J, editors. Evolution of the House mouse (Cambridge Studies in Morphology and Molecules: New Paradigms in Evolutionary Biology). UK: Cambridge University Press; 2012. p. 35–64.View ArticleGoogle Scholar
- Prager EM, Orrego C, Sage RD. Genetic variation and phylogeography of central Asian and other house mice, including a major new mitochondrial lineage in Yemen. Genetics. 1998;150:835–61.PubMed CentralPubMedGoogle Scholar
- Rajabi-Maham H, Orth A, Siahsarvie R, Boursot P, Darvish J, Bonhomme F. The south-eastern house mouse Mus musculus castaneus (Rodentia: Muridae) is a polytypic subspecies. Biol J Linn Soc. 2012;107:295–306.View ArticleGoogle Scholar
- Suzuki H, Nunome M, Kinoshita G, Aplin KP, Vogel P, Kryukov AP, et al. Evolutionary and dispersal history of Eurasian house mice Mus musculus clarified by more extensive geographic sampling of mitochondrial DNA. Heredity. 2013;111(5):375–90.View ArticlePubMed CentralPubMedGoogle Scholar
- Cucchi T, Eszter Kovács Z, Berthon R, Orth A, Bonhomme F, Evin A, et al. On the trail of Neolithic mice and men towards Transcaucasia: zooarchaeological clues from Nakhchivan (Azerbaijan). Biol J Linn Soc. 2013;108:917–28.View ArticleGoogle Scholar
- Duplantier J-M, Orth A, Catalan J, Bonhomme F. Evidence for a mitochondrial lineage originating from the Arabian peninsula in the Madagascar House Mouse (Mus musculus). Heredity. 2002;89:154–8.View ArticlePubMedGoogle Scholar
- Orth A, Belkhir K, Britton-Davidian J, Boursot P, Benazzou T, Bonhomme F. Hybridation naturelle entre deux espèces sympatriques de souris M. musculus domesticus L. et M. spretus Lataste. C R Biol. 2002;325:89–97.View ArticlePubMedGoogle Scholar
- Bonhomme F, Orth A, Cucchi T, Rajabi-Maham H, Catalan J, Boursot P, et al. Genetic differentiation of the house mouse around the Mediterranean basin: matrilineal footprints of early and late colonization. Proc R Soc B. 2011;278:1034–43.View ArticlePubMed CentralPubMedGoogle Scholar
- Orth A, Lyapunova E, Kandaurov A, Boissinot S, Boursot P, Vorontsov N, et al. L’espèce polytypique Mus musculus en Transcaucasie. C R Biol. 1996;319:435–41.Google Scholar
- Ihle S, Ravaoarimanana I, Thomas M, Tautz D. An analysis of signatures of selective sweeps in natural populations of the house mouse. Mol Biol Evol. 2006;23:790–7.View ArticlePubMedGoogle Scholar
- Din W, Anand R, Boursot P, Darviche D, Dod B, Jouvin-Marche E, et al. Origin and radiation of the house mouse: Clues from nuclear genes. J Evol Biol. 1996;9:519–39.View ArticleGoogle Scholar
- Rajabi-Maham H, Orth A, Bonhomme F. Phylogeography and postglacial expansion of Mus musculus domesticus inferred from mitochondrial DNA coalescent, from Iran to Europe. Mol Ecol. 2008;17:627–41.View ArticlePubMedGoogle Scholar
- Auffray J-C, Orth A, Catalan J, Gonzalez J, Desmarais E, Bonhomme F. Phylogenetic position and description of a new species of subgenus Mus (Rodentia, Mammalia) from Thailand. Zool Scr. 2003;32:119–27.View ArticleGoogle Scholar
- Cucchi T, Orth A, Auffray J-C, Renaud S, Fabre L, Catalan J, et al. A new endemic species of the subgenus Mus (Rodentia, Mammalia) on the Island of Cyprus. Zootaxa. 2006;1241:1–36.Google Scholar
- Teschke M, Mukabayire O, Wiehe T, Tautz D. Identification of selective sweeps in closely related populations of the house mouse based on microsatellite scans. Genetics. 2008;180:1537–45.View ArticlePubMed CentralPubMedGoogle Scholar
- Ihle S. Detecting genes involved in selective sweeps within populations of the house mouse species complex-a multi-locus candidate gene approach. PhD thesis. University of Cologne;2004. http://kups.ub.uni-koeln.de/view/creators/Ihle=3ASonja=3A=3A.html.
- Darvish J, Orth A, Bonhomme F. Genetic transition in the house mouse, Mus musculus of Eastern Iranian Plateau. Folia Zool. 2006;55:349–57.Google Scholar
- Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155:945–59.PubMed CentralPubMedGoogle Scholar
- Evanno G, Regnaut S, Goudet J. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol. 2005;14:2611–20.View ArticlePubMedGoogle Scholar
- Gao H, Bryc B, Bustamente CD. On Identifying the optimal number of population clusters via the deviance information criterion. PLoS One. 2011;6:e21014. doi:10.1371/journal.pone.0021014.View ArticlePubMed CentralPubMedGoogle Scholar
- Nunome M, Ishimori C, Aplin KP, Tsuchiya K, Yonekawa H, Moriwaki K, et al. Detection of recombinant haplotypes in wild mice (Mus musculus) provides new insights into the origin of Japanese mice. Mol Ecol. 2010;19:2474–89.PubMedGoogle Scholar
- Karn RC, Orth A, Bonhomme F, Boursot P. The complex history of a gene proposed to participate in a sexual isolation mechanism in house mice. Mol Biol Evol. 2002;19:462–71.View ArticlePubMedGoogle Scholar
- Duvaux L, Belkhir K, Boulesteix M, Boursot P. Isolation and gene flow: inferring the speciation history of European house mouse. Mol Ecol. 2011;20:5248–64.View ArticlePubMedGoogle Scholar
- Siahsarvie R, Auffray J-C, Darvish J, Rajabi-Maham H, Yu HT, Agret S, et al. Patterns of morphological evolution in the mandible of the house mouse Mus musculus (Rodentia: Muridae). Biol J Linn Soc. 2012;105:635–47.View ArticleGoogle Scholar
- Staubach F, Lorenc A, Messer PW, Tang K, Petrov DA, Tautz D. Genome patterns of selection and introgression of haplotypes in natural populations of the house mouse (Mus musculus). PLoS Genet. 2012;8:e1002891. doi:10.1371/journal.pgen.1002891.View ArticlePubMed CentralPubMedGoogle Scholar
- Phifer-Rixey M, Bomhoff M, Nachman M. Genome-wide patterns of differentiation among house mouse subspecies. Genetics. 2014;198:283–97.View ArticlePubMedGoogle Scholar
- Baird SJE, Macholan M. What can the Mus musculus musculus/M. m. domesticus hybrid zone tell us about speciation? In: Macholán M, Baird SJE, Munclinger P, Piálek J, editors. Evolution of the House mouse (Cambridge Studies in Morphology and Molecules: New Paradigms in Evolutionary Biology). UK: Cambridge University Press; 2012. p. 334–72.View ArticleGoogle Scholar
- Thomas M, Moeller F, Wiehe T, Tautz D. A pooling approach to detect signatures of selective sweeps in genome scans using microsatellites. Mol Ecol Notes. 2007;7:400–3.View ArticleGoogle Scholar
- Bayer T. Duplicated genes and their relevance in the process of speciation. PhD thesis. University of Cologne;2009. http://kups.ub.uni-koeln.de/2751/.
- Belkhir K, Borsa P, Chikhi L, Raufaste N, Bonhomme F. GENETIX 4.05, logiciel sous Windows TM pour la génétique des populations. Laboratoire Génome, Populations, Interactions, CNRS UMR 5171. Montpellier (France): Université de Montpellier II; 1996.Google Scholar
- Excoffier L, Lischer HEL. Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour. 2010;10:564–7.View ArticlePubMedGoogle Scholar
- Felsenstein J. PHYLIP (Phylogeny Inference Package) version 3.6. Distributed by the author. Seattle: Department of Genome Sciences, University of Washington; 2005.Google Scholar
- Langella O. Populations, 1.2.30 Copyright (C). CNRS UPR9034; 1999. Available at http://bioinformatics.org/~tryphon/populations/.Google Scholar
- Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6. Molecular Evolutionary Genetics Analysis Version 6.0; 2013 http://www.megasoftware.net/.
- Earl DA, von Holdt BM. STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv Genet Resour. 2012;4:359–61.View ArticleGoogle Scholar
- Jakobsson M, Rosenberg NA. CLUMPP: a cluster matching and permutation program for dealing with label switching and multimodality in analysis of population structure. Bioinformatics. 2004;23:1801–6.View ArticleGoogle Scholar
- Rosenberg NA. Distruct: a program for the graphical display of population structure. Mol Ecol Notes. 2004;4:137–8.View ArticleGoogle Scholar
- Tamura K, Nei M. Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol Biol Evol. 1993;10:512–26.PubMedGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.