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

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

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

House mice (Mus musculus) are human commensals and have served as a primary model in biomedical, ecological and evolutionary research. Although there is detailed knowledge of the biogeography of house mice in Europe, little is known of the history of house mice in China, despite the fact that China encompasses an enormous portion of their range. In the present study, 535 house mice caught from 29 localities in China were studied by sequencing the mitochondrial D‐loop and genotyping 10 nuclear microsatellite markers distributed on 10 chromosomes. Phylogenetic analyses revealed two evolutionary lineages corresponding to Mus musculus castaneus and Mus musculus musculus in the south and north, respectively, with the Yangtze River approximately representing the boundary. More detailed analyses combining published sequence data from mice sampled in neighbouring countries revealed the migration routes of the two subspecies into China: M. m. castaneus appeared to have migrated through a southern route (Yunnan and Guangxi), whereas M. m. musculus entered China from Kazakhstan through the north‐west border (Xinjiang). Bayesian analysis of mitochondrial sequences indicated rapid population expansions in both subspecies, approximately 4650–9300 and 7150–14 300 years ago for M. m. castaneus and M. m. musculus, respectively. Interestingly, the migration routes of Chinese house mice coincide with the colonization routes of modern humans into China, and the expansion times of house mice are consistent with the development of agriculture in southern and northern China, respectively. Finally, our study confirmed the existence of a hybrid zone between M. m. castaneus and M. m. musculus in China. Further study of this hybrid zone will provide a useful counterpart to the well‐studied hybrid zone between M. m. musculus and Mus musculus domesticus in central Europe.     

A molecular characterization of the charismatic Faroe house mouse

Faroe house mice are a ‘classic’ system of rapid and dramatic morphological divergence highlighted by J. S. Huxley during the development of the Modern Synthesis. In the present study, we characterize these charismatic mice using modern molecular techniques, examining specimens from all Faroe islands occupied by mice. The aims were to classify the mice within the modern house mouse taxonomy (i.e. as either Mus musculus domesticus or Mus musculus musculus) using four molecular markers and a morphological feature, and to examine the genetic diversity and possible routes of colonization using mitochondrial (mt) control region DNA sequences and microsatellite data (15 loci). Mice on the most remote islands were characterized as M. m. domesticus and exhibited exceptionally low genetic diversity, whereas those on better connected islands were more genetically diverse and had both M. m. musculus and M. m. domesticus genetic elements, including one population which was morphologically M. m. musculus‐like. The mtDNA data indicate that the majority of the mice had their origins in south‐western Norway (or possibly southern Denmark/northern Germany), and probably arrived with the Vikings, earlier than suggested by Huxley. The M. m. musculus genetic component appears to derive from recent mouse immigration from Denmark. © 2011 The Linnean Society of London, Biological Journal of the Linnean Society, 2011, 102 , 471–482.     

Polymorphic HSRs in chromosome 1 of the two semispecies Mus musculus musculus and M. m. domesticus have a common origin in an ancestral population

HSRs (homogeneously staining regions) are the cytological correlates of DNA amplification. In the house mouse, Mus musculus, many populations are polymorphic for the presence or absence of HSRs on chromosome 1. In the semispecies M. m. domesticus the amplified DNA is present within one HSR, whereas in M. m. musculus chromosomes 1 with two HSRs are found. Hybridization of HSR-specific probes to Southern blots of HSR-carrying genomic DNAs from different localities and semispecies revealed similar complex band patterns. the remaining variation is restricted to sequences with a low degree of amplification. Variation is higher between semispecies than within one semispecies. It is assumed that HSRs are derived from one original amplification event and that unequal recombination is the mechanism underlying the length variation of HSRs present today in both semispecies. Evidence from G-banding and in situ hybridization shows that the two HSRs of M. m. musculus originated from a single HSR by means of a paracentric inversion, where one break-point was located within the single HSR and the second outside the HSR. As a consequence of the paracentric inversion the two HSRs of M. m. musculus are permanently linked together. Since exchange of genes between the two semispecies is restricted to a narrow hybrid zone the amplification that gave rise to the HSR most probably occurred prior to the divergence into the semispecies M. m. domesticus and M. m. musculus about 1 million years ago.by D. Schweizer     

Location of a contact zone between Mus musculus domesticus and M. m. domesticus with M. m. castaneus mtDNA in southern New Zealand

The house mouse, Mus musculus, was first introduced into New Zealand in significant numbers in the early to mid nineteenth century, with genomic components from different sources of the three subspecies M. m. domesticus, M. m. musculus and M. m. castaneus. M. m. domesticus is now widely distributed in New Zealand, with genomic and morphological evidence of M. m. musculus in a few scattered locations. M. m. domesticus/M. m. castaneus hybrids are dominant in the southern third of the South Island. We anticipated that there should be a definable southern contact zone between pure M. m. domesticus and M. m. domesticus/M. m. castaneus hybrids. We tested this hypothesis by screening 170 DNA samples from mice collected in the southern South Island, using a PCR technique which rapidly distinguishes the mitochondrial genomes of the three subspecies.All mice sampled from in or north of Lincoln (43.63° S) had only M. m. domesticus mtDNA, whereas all those from or further south than Hook (44.68° S) had M. m. castaneus mtDNA. Between the two sites, mice carrying mtDNA of both subspecies were found, sometimes in the same building. On present data, this contact zone extends approximately 50 km north to south and some 30 km inland. Classical tests with three nuclear DNA markers confirmed earlier work showing that the nuclear genomes of all mice appeared to be predominantly domesticus-like.We conclude that if purebred M. m. castaneus mice did originally reach New Zealand, extensive backcrossing with M. m. domesticus has made the castaneus nuclear genome virtually undetectable with the tests that we employ.     

Polymerase chain reaction multiplexing of microsatellites and single nucleotide polymorphism markers for quantitative trait loci mapping of wild house mice

We tested 96 microsatellites and 10 single nucleotide polymorphisms for their allelic distribution in two subspecies of the house mouse, Mus musculus musculus and M. m. domesticus. Sixty‐two microsatellites discriminated strain‐specific differences among nine wild‐derived ‘musculus’ and ‘domesticus’ and three ‘classical’ laboratory strains. For efficient genotyping, we optimized multiplex conditions using five microsatellites per polymerase chain reaction. All 10 single nucleotide polymorphisms were also optimized for simultaneous analysis in one reaction using SNaPshot multiplex. The uniform distribution of markers on autosomes and on the X chromosome makes these panels potentially useful tools for quantitative trait loci mapping of wild house mice.     

Genetic differentiation of subspecies of the house mouse Mus musculus and their taxonomic relationships inferred from RAPD-PCR data

Genetic differentiation of six subspecies of the house mouse Mus musculus (Mus musculus musculus, M. m. domesticus, M. m. castaneus, M. m. gansuensis, M. m. wagneri, and M. m. ssp. (bactrianus?) was examined using RAPD-PCR analysis. In all, 373 loci of total length of about 530 kb were identified. Taxonspecific molecular markers were detected and the levels of genetic differences among the subspecies were estimated. Different degree of subspecific genetic differentiation was shown. The most similar subspecies pairs were M. m. castaneus-M. m. domesticus and M. m. musculus-M. m. gansuensis. In our phylogenetic reconstruction, M. m. wagneri proved to be most different from all the other subspecies. Genetic distances between it and other subspecies were two-to threefold higher than those between the “good”species of the subgenus Mus (e.g., between M. m. musculus and M. spicilegus, M. musculus and M. abbotti). The estimates of genetic similarity and the phylogenetic relationships between six house mouse subspecies inferred from RAPD partially conformed to the results based on cytogenetic and allozyme data. However, they were considerably different from phylogenetic reconstructions based on sequencing of the control mtDNA region, which reflects mutual inconsistency of different systems of inheritance.     

Gastrointestinal microbiota of wild and inbred individuals of two house mouse subspecies assessed using high‐throughput parallel pyrosequencing

The effects of gastrointestinal tract microbiota (GTM) on host physiology and health have been the subject of considerable interest in recent years. While a variety of captive bred species have been used in experiments, the extent to which GTM of captive and/or inbred individuals resembles natural composition and variation in wild populations is poorly understood. Using 454 pyrosequencing, we performed 16S rDNA GTM barcoding for 30 wild house mice (Mus musculus) and wild‐derived inbred strain mice belonging to two subspecies (M. m. musculus and M. m. domesticus). Sequenced individuals were selected according to a 2 × 2 experimental design: wild (14) vs. inbred origin (16) and M. m. musculus (15) vs. M. m. domesticus (15). We compared alpha diversity (i.e. number of operational taxonomic units – OTUs), beta diversity (i.e. interindividual variability) and microbiota composition across the four groups. We found no difference between M. m. musculus and M. m. domesticus subspecies, suggesting low effect of genetic differentiation between these two subspecies on GTM structure. Both inbred and wild populations showed the same level of microbial alpha and beta diversity; however, we found strong differentiation in microbiota composition between wild and inbred populations. Relative abundance of ~ 16% of OTUs differed significantly between wild and inbred individuals. As laboratory mice represent the most abundant model for studying the effects of gut microbiota on host metabolism, immunity and neurology, we suggest that the distinctness of laboratory‐kept mouse microbiota, which differs from wild mouse microbiota, needs to be considered in future biomedical research.     

Using the Mus musculus hybrid zone to assess covariation and genetic architecture of limb bone lengths

Two subspecies of the house mouse, Mus musculus domesticus and Mus musculus musculus, meet in a narrow contact zone across Europe. Mice in the hybrid zone are highly admixed, representing the full range of mixed ancestry from the two subspecies. Given the distinct morphologies of these subspecies, these natural hybrids can be used for genomewide association mapping at sufficiently high resolution to directly infer candidate genes. We focus here on limb bone length differences, which is of special interest for understanding the evolution of developmentally correlated traits. We used 172 first‐generation descendants of wild‐caught mice from the hybrid zone to measure the length of stylopod (humerus/femur), zeugopod (ulna/tibia) and autopod (metacarpal/metatarsal) elements in skeletal CT scans. We find phenotypic covariation between limb elements in the hybrids similar to patterns previously described in Mus musculus domesticus inbred strains, suggesting that the hybrid genotypes do not influence the covariation pattern in a major way. Mapping was performed using 143,592 SNPs and identified several genomic regions associated with length differences in each bone. Bone length was found to be highly polygenic. None of the candidate regions include the canonical genes known to control embryonic limb development. Instead, we are able to identify candidate genes with known roles in osteoblast differentiation and bone structure determination, as well as recently evolved genes of, as yet, unknown function.     

The uncharacterized gene 1700093K21Rik and flanking regions are correlated with reproductive isolation in the house mouse,Mus musculus

Reproductive barriers exist between the house mouse subspecies, Mus musculus musculus and M. m. domesticus, members of the Mus musculus species complex, primarily as a result of hybrid male infertility, and a hybrid zone exists where their ranges intersect in Europe. Using single nucleotide polymorphisms (SNPs) diagnostic for the two taxa, the extent of introgression across the genome was previously compared in these hybrid populations. Sixty-nine of 1316 autosomal SNPs exhibited reduced introgression in two hybrid zone transects suggesting maladaptive interactions among certain loci. One of these markers is within a region on chromosome 11 that, in other studies, has been associated with hybrid male sterility of these subspecies. We assessed sequence variation in a 20 Mb region on chromosome 11 flanking this marker, and observed its inclusion within a roughly 150 kb stretch of DNA showing elevated sequence differentiation between the two subspecies. Four genes are associated with this genomic subregion, with two entirely encompassed. One of the two genes, the uncharacterized 1700093K21Rik gene, displays distinguishing features consistent with a potential role in reproductive isolation between these subspecies. Along with its expression specifically within spermatogenic cells, we present various sequence analyses that demonstrate a high rate of molecular evolution of this gene, as well as identify a subspecies amino acid variant resulting in a structural difference. Taken together, the data suggest a role for this gene in reproductive isolation.     

Patterns of morphological evolution in the mandible of the house mouse Mus musculus (Rodentia: Muridae)

The worldwide distributed house mouse, Mus musculus, is subdivided into at least three lineages, Mus musculus musculus, Mus musculus domesticus, and Mus musculus castaneus. The subspecies occur parapatrically in a region considered to be the cradle of the species in Southern Asia (‘central region’), as well as in the rest of the world (‘peripheral region’). The morphological evolution of this species in a phylogeographical context is studied using a landmark‐based approach on mandible morphology of different populations of the three lineages. The morphological variation increases from central to peripheral regions at the population and subspecific levels, confirming a centrifugal sub‐speciation within this species. Furthermore, the outgroup comparison with sister species suggests that M. musculus musculus and populations of all subspecies inhabiting the Iranian plateau have retained a more ancestral mandible morphology, suggesting that this region may represent one of the relevant places of the origin of the species. Mus musculus castaneus, both from central and peripheral regions, is morphologically the most variable and divergent subspecies. Finally, the results obtained in the present study suggest that the independent evolution to commensalism in the three lineages is not accompanied by a convergence detectable on jaw morphology. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 105 , 635–647.     

Counterselection on sex chromosomes in the Mus musculus European hybrid zone

The extent to which alleles can disperse across a hybrid zone depends on the selection they are subjected to in the hybrid genetic background or, for those that are selectively neutral, on their ability to escape from the unfavourable environment by recombination. Three markers spanning a 45 cM segment in the center of the X chromosome were used to investigate the degree to which selection against X chromosome linked genes helps to maintain the barrier to gene flow in the hybrid zone between Mus musculus domesticus and M. m. musculus in Denmark. The introgression of all the sex chromosome specific markers was more limited than that of the autosomal enzymes (Idh1, Amy, Gpd1, Pgm1, Es1, Es2, Mpi, Np1, Es10, Sod1) and the mitochondrial DNA. The cline for DXPas2, which is in the center of the X chromosome, is extremely steep and shows that certain genes located in this region are strongly selected against in the hybrid background. The clines of the other two X-linked markers, Hprt and DXPas1, and of the Y chromosome are not as abrupt and all three have similar asymmetric introgression patterns. Although the musculus variants appear to behave in much the same way as those of the autosomal genes, the domesticus variants do not introgress. The results show that X-linked and to a lesser extent Y-linked genes are more strongly selected against in the hybrid genome than the mitochondrial genome or the different autosomal loci. This suggests that co-adapted gene systems involving the sex chromosomes may play an important role in the hybrid breakdown between the two subspecies.     

Origin and radiation of the house mouse: mitochondrial DNA phylogeny

On the basis of patterns of allele frequency variation in nuclear genes (Din et al., in press) it has been proposed that the house mouse M. musculus originated in the northern Indian subcontinent, from where it radiated in several directions to form the well-described peripheral subspecies (M. m. domesticus, M. m. musculus and M. m. castaneus). Here we use a mitochondrial DNA (mtDNA) phylogeny to test this hypothesis and to analyse the historical and demographic events that have accompanied this differentiation. This marker also provides a powerful means to check for genetic continuity between the central and peripheral populations. We studied restriction site polymorphism of samples from India and the Middle East as well as samples from the rest of Eurasia and northern Africa. M. m. domesticus and M. m. musculus are both monophyletic for mtDNA and belong to the subspecies-specific mtDNA lineages that have been described previously. Average nucleotide diversity is low in M. m. musculus (0.2–5%). It is not only higher in M. m. domesticus (0.7–0.9%) but the distribution of pairwise divergence is wider, and the rate of evolution in this branch appears to be higher than in M. m. musculus. The nucleotide diversity found in M. m. castaneus (0.4%) is due to the existence of two rather divergent linages with little intralineage variation. These two lineages are part of a diversified bush of the phylogenetic tree that also comprises several previously undescribed branches and includes all samples from the northern Indian subcontinent and Iran. The degree of diversity found in each of the samples from this region is high (1.2–2.4%) although they come from small geographic areas. This agrees well with the idea that the origin of the radiation was in the northern Indian subcontinent. However, as neither haplotypes on the M. m. domesticus nor on the M. m. musculus branches were found in this region, there appear to be important phylogeographic discontinuities between this central region and these peripherial subspecies. On the basis of the present result and the nuclear data (Din et al., in press), we propose that M. musculus originated in the north of the indian subcontinent. Our calibration of the evolutionary rate of mtDNA in mice suggests that the mouse settlement in this region could be as old as 900 000 years. Possibly from there, a first radiation could have reach the Middle East and the Caspian Sea, where the M. m. domesticus and M. m. musculus lineages, respectively, would have started to differentiate a few hundred thousand years ago, and from where they could have colonised the peripheral part of their ranges only recently.M. m. castaneus appears from its mtDNA to be recent offshoot of the northern Indian population. This multiple and gradual radiation ultimately led to recent peripheral secondary contacts, such as the well-known European hybrid zone.     

Female House Mice do not Differ in Their Exploratory Behaviour from Males

Dispersal is an important ecological trait affecting genetic variation and dynamics of populations. Hence, the exploratory behaviour prior to actual dispersal may be crucial for potentially dispersing individuals. In mammals, females are traditionally seen as the more philopatric sex and dispersal as male‐biased behaviour, and so behavioural strategies related to the exploration of novel resources should be differentially expressed in males and females. In addition, due to sexual selection exploratory strategies may be expected to vary according to females’ reproductive phase. We employed a standard open‐field test as an approximation of the first phase of dispersal, using adult house mice representing two subspecies, M. m. musculus and M. m. domesticus. We tested the prediction that exploration of neutral area varies in females during different phases of the oestrus cycle and is different between both sexes and subspecies. We expected to find higher exploration in males, as the more dispersing sex and less pronounced subspecies‐specific differences in females than in males. We found no significant effect of the oestrous phase on any of the parameters of the exploratory behaviour measured. Sexual dimorphism was found only in latency to enter the arena in M. m. domesticus where females hesitated longer to enter a new area than males. Significant subspecies‐specific differences were found in three of four tested exploration parameters, so we conclude that females of both subspecies follow similar strategies to those displayed by males. Musculus mice show shorter latency to enter a new area, but once inside, domesticus mice explore the arena significantly longer, with less frequent retreats to a shelter. Our results thus highlight that the role of female dispersal in interdemic gene flow should not be neglected.     

Variation in mandible shape and body size of house mice Mus musculus in five separate New Zealand forest habitats

This study investigates variation in house mouse Mus musculus body size and mandible shape across New Zealand, using geometric morphometrics and biomechanical advantage analyses. The Mus phylogroups currently known in New Zealand include Mus musculus domesticus, M. m. musculus and M. m. castaneus. We examined samples of house mice inhabiting five different podocarp and beech forest environments across the North and South Islands (Pureora Forest, Zealandia Wildlife Sanctuary, Craigieburn Forest Park, Eglinton Valley and Hollyford Valley). Significant variation in mandible shape and body size was found between all five forest populations. South Island mice had larger bodies and greater mechanical advantage in the temporalis muscle compared with their North Island counterparts. Zealandia Sanctuary mouse mandibles were broader and shorter than South Island mouse mandibles, and had greater masseter muscle advantage. Centroid size and body weight, but not head-body length, varied significantly with two distinct genetic haplotypes. Finally, annual rainfall was the most significant covariate with mandible shape. Higher rainfall locations were generally associated with soft-food related mandible shapes, while lower rainfall correlated with hard-food mandible shapes. This preliminary investigation provides the framework for further research into mandible shape and body size variation in New Zealand house mice.     

Transgressive segregation in a behavioural trait? Explorative strategies in two house mouse subspecies and their hybrids

Hybrid zones between genetically diverged populations are widespread among animals and plants. Their dynamics usually depend on selection against admixture and dispersal of parental forms in the zone. Although indirect estimates of selection have been the target of many studies, dispersal has been neglected. In this study we carried out open field experiments to test whether males of two house mouse subspecies, Mus musculus musculus and Mus musculus domesticus, differ in their propensity to disperse and in their character of exploration. We tested wild‐caught males and males of two wild‐derived inbred strains. In addition, we examined reciprocal F₁ crosses to test the prediction that these hybrids display intermediate behaviours. We revealed that M. m. musculus males were less hesitant to enter the experimental arena than were M. m. domesticus males, but once inside the arena their movements were more timid. F₁ males differed from both parental strains, with longer latencies to enter the arena, but explored the arena in a similar fashion as the M. m. domesticus males, thus displaying transgressive behavioural phenotypes. These results contribute to our knowledge of behavioural divergence between the mouse subspecies, and add a new facet to the study of speciation. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, ●●, ●●–●●.     

DISSECTING THE GENETIC ARCHITECTURE OF F1 HYBRID STERILITY IN HOUSE MICE

Hybrid sterility as a postzygotic reproductive isolation mechanism has been studied for over 80 years, yet the first identifications of hybrid sterility genes in Drosophila and mouse are quite recent. To study the genetic architecture of F₁ hybrid sterility between young subspecies of house mouse Mus m. domesticus and M. m. musculus, we conducted QTL analysis of a backcross between inbred strains representing these two subspecies and probed the role of individual chromosomes in hybrid sterility using the intersubspecific chromosome substitution strains. We provide direct evidence that the asymmetry in male infertility between reciprocal crosses is conferred by the middle region of M. m. musculus Chr X, thus excluding other potential candidates such as Y, imprinted genes, and mitochondrial DNA. QTL analysis identified strong hybrid sterility loci on Chr 17 and Chr X and predicted a set of interchangeable autosomal loci, a subset of which is sufficient to activate the Dobzhansky–Muller incompatibility of the strong loci. Overall, our results indicate the oligogenic nature of F₁ hybrid sterility, which should be amenable to reconstruction by proper combination of chromosome substitution strains. Such a prefabricated model system should help to uncover the gene networks and molecular mechanisms underlying hybrid sterility.     

Of mice and the ‘Age of Discovery’: the complex history of colonization of the Azorean archipelago by the house mouse (Mus musculus) as revealed by mitochondrial DNA variation

Humans have introduced many species onto remote oceanic islands. The house mouse (Mus musculus) is a human commensal and has consequently been transported to oceanic islands around the globe as an accidental stowaway. The history of these introductions can tell us not only about the mice themselves but also about the people that transported them. Following a phylogeographic approach, we used mitochondrial D‐loop sequence variation (within an 849‐ to 864‐bp fragment) to study house mouse colonization of the Azores. A total of 239 sequences were obtained from all nine islands, and interpretation was helped by previously published Iberian sequences and 66 newly generated Spanish sequences. A Bayesian analysis revealed presence in the Azores of most of the D‐loop clades previously described in the domesticus subspecies of the house mouse, suggesting a complex colonization history of the archipelago as a whole from multiple geographical origins, but much less heterogeneity (often single colonization?) within islands. The expected historical link with mainland Portugal was reflected in the pattern of D‐loop variation of some of the islands but not all. A more unexpected association with a distant North European source area was also detected in three islands, possibly reflecting human contact with the Azores prior to the 15th century discovery by Portuguese mariners. Widening the scope to colonization of the Macaronesian islands as a whole, human linkages between the Azores, Madeira, the Canaries, Portugal and Spain were revealed through the sharing of mouse sequences between these areas. From these and other data, we suggest mouse studies may help resolve historical uncertainties relating to the ‘Age of Discovery’.