Development of unique house mouse resources suitable for evolutionary studies of speciation

Two house mouse subspecies, Mus musculus domesticus and Mus musculus musculus, form a hybrid zone in Europe and represent a suitable model for inferring the genes contributing to isolation barriers between parental taxa. Despite long-term intensive studies of this hybrid zone, we still know relatively little about the causes and mechanisms maintaining the 2 taxa as separate subspecies; therefore, to gain insight into this process, we developed 8 wild-derived inbred house mouse strains. In order to produce strains as pure domesticus or musculus genomes as possible, the individuals used to establish the breeding colony for the 3 domesticus and 2 of the musculus strains were captured in the Czech Republic from wild populations at extreme western and eastern edges of the subspecific contact zone, respectively. The remaining 3 musculus strains were bred from mice captured about 250 km east of the hybrid zone. Genetic analysis based on 361 microsatellite loci showed that 82% of these markers are diagnostic for either the musculus or the domesticus strains. In order to demonstrate the potential utility of this genetic differentiation in such strains, phenotypic variation was scored for 2 strains from opposite edges of the hybrid zone and significant differences in morphology, reproductive performance, in vitro immune responses, mate choice based on urinary signals, and aggressiveness were found. In addition, the 3 strains derived from musculus populations far from the hybrid zone display significant differences in polymorphism in hybrid male sterility when crossed with the laboratory strains C57BL/6 or C57BL/10, which have a predominantly domesticus genome. Although further studies will be necessary to demonstrate intersubspecific differences, all analyses presented here indicate that these newly developed house mouse strains represent a powerful tool for elucidating the genetic basis of isolation barriers in hybrid zones and for studying speciation in general.     

Strong premating divergence in a unimodal hybrid zone between two subspecies of the house mouse

Although selection against hybridization is expected to generate prezygotic divergence in unimodal hybrid zones, such a pattern has been seldom described. This study aims to better understand how prezygotic mechanisms may evolve in such zones. We investigated prezygotic divergence between populations of two subspecies of mice (Mus musculus musculus and M. m. domesticus) located at the edges of their unimodal hybrid zone in Denmark, and we developed an original multiple-population choice-test design, which allows assessment of within and between subspecies variation. Our study demonstrates that a strong assortative preference characterises one of the two subspecies (musculus) and that urinary signals are involved in this subspecies recognition. Taking into account the specific genetic and geographical characteristics of the Danish hybrid zone, we discuss the influence of the above pattern on its fate and the mechanisms that could have favoured this prezygotic divergence, among which the role of recombined populations constituting the core of the zone.     

Chromosomal introgression in house mice from the hybrid zone between M. m. domesticus and M. m. musculus in Denmark

The recent discovery of Robertsonian (Rb) translocations in Danish mice from the hybrid zone between Mus musculus musculus and M. m. domesticus stimulated the chromosomal analysis of populations along a north-south transect through this zone. G-Banding identified the Rb fusions as Rb(3.8), Rb(2.5) and Rb(6.9). The cytogenetic results show that there is a gradual decrease in the number of fusions as one proceeds north, the translocations abruptly ending in populations from the centre of the hybrid zone determined by seven diagnostic allozymic markers. These results indicate that Rb fusions are present only in domesticus or predominantly domesticus-genotype mice and that they do not introgress into M. m. musculus . To test if genie incompatibilities between the musculus genetic background and Rb fusions were involved in the systematic elimination of the latter, predominantly musculus mice from the hybrid zone were crossed with Rb domesticus mice carrying Rb(3.8). The karyotypic analysis of the progeny showed no distortion of the transmission ratio of this fusion.
The chromosomal and allozymic analysis of these mice further indicates that (i) recombination is not suppressed between metacentrics and their acrocentric homologues and (ii) specific domesticus chromosomal segments are tolerated in the musculus genomes whereas the Rb centromeres are not.     

Empirical study of hybrid zone movement

Hybrid zones are 'natural laboratories' for studying the origin, maintenance and demise of species. Theory predicts that hybrid zones can move in space and time, with significant consequences for both evolutionary and conservation biology, though such movement is often perceived as rare. Here, a review of empirical studies of moving hybrid zones in animals and plants shows 23 examples with observational evidence for movement, and a further 16 where patterns of introgression in molecular markers could be interpreted as signatures of movement. The strengths and weaknesses of methods used for detecting hybrid zone movement are discussed, including long-term replicated sampling, historical surveys, museum/herbarium collections, patterns of relictual populations and introgression of genetic markers into an advancing taxon. Factors governing hybrid zone movement are assessed in the light of the empirical studies, including environmental selection, competition, asymmetric hybridization, dominance drive, hybrid fitness, human activity and climate change. Hybrid zone movement means that untested assumptions of stability in evolutionary studies on hybrid zone can lead to mistaken conclusions. Movement also means that conservation effort aimed at protecting against introgression could unwittingly favour an invading taxon. Moving hybrid zones are of wide interest as examples of evolution in action and possible indicators of environmental change. More long-term experimental studies are needed that incorporate reciprocal transplants, hybridization experiments and surveys of molecular markers and population densities on a range of scales.     

Historical hybrid zone movement: More pervasive than appreciated

Hybrid zones are established where two divergent populations meet and interbreed, but experience some reproductive isolation. If one population expands its range at the expense of the other, their hybrid zone moves. While hybrid zone movement is generally considered to be uncommon and insignificant, recent studies challenge this idea. The commonality of contemporary hybrid zone movement—with shifts in hybrid zones tracked over years to decennia—cannot be disputed, given the many examples available. Cases of historical hybrid zone movement—covering centuries or millennia of mobility—are accumulating, with movement having been inferred from five lines of evidence: (1) range shifts documented in the fossil/pollen record; (2) distribution dynamics derived from species distribution modelling; (3) enclaves of a displaced population persisting inside the range of an expanding one; (4) a peak of linkage disequilibrium at the leading edge of a moving hybrid zone; and (5) genome‐wide genetic traces of a displaced population, left behind in an expanding one. While most of these lines of evidence are not straightforward to interpret and/or broadly applicable, the latter—a genomic footprint of hybrid zone movement—promises to be particularly suitable to determine whether a hybrid zone has been on the move since its inception. I argue that historical hybrid zone movement is likely to be prevalent and deserves wider acknowledgement in historical biogeography.     

Differential patterns of introgression across the X chromosome in a hybrid zone between two species of house mice

A complete understanding of the speciation process requires the identification of genomic regions and genes that confer reproductive barriers between species. Empirical and theoretical research has revealed two important patterns in the evolution of reproductive isolation in animals: isolation typically arises as a result of disrupted epistatic interactions between multiple loci and these disruptions map disproportionately to the X chromosome. These patterns suggest that a targeted examination of natural gene flow between closely related species at X-linked markers with known positions would provide insight into the genetic basis of speciation. We take advantage of the existence of genomic data and a well-documented European zone of hybridization between two species of house mice, Mus domesticus and M. musculus, to conduct such a survey. We evaluate patterns of introgression across the hybrid zone for 13 diagnostic X-linked loci with known chromosomal positions using a maximum likelihood model. Interlocus comparisons clearly identify one locus with reduced introgression across the center of the hybrid zone, pinpointing a candidate region for reproductive isolation. Results also reveal one locus with high frequencies of M. domesticus alleles in populations on the M. musculus side of the zone, suggesting the possibility that positive selection may act to drive the spread of alleles from one species on to the genomic background of the other species. Finally, cline width and cline center are strongly positively correlated across the X chromosome, indicating that gene flow of the X chromosome may be asymmetrical. This study highlights the utility of natural populations of hybrids for mapping speciation genes and suggests that the middle of the X chromosome may be important for reproductive isolation between species of house mice.     

A genetically distinct hybrid zone occurs for two globally invasive mosquito fish species with striking phenotypic resemblance

Hybrid zones allow for the investigation of incipient speciation and related evolutionary processes of selection, gene flow, and migration. Interspecific dynamics, like competition, can impact the size, shape, and directional movement of species in hybrid zones. Hybrid zones contribute to a paradox for the biological species concept because interbreeding between species occurs while parental forms remain distinct. A long‐standing zone of intergradation or introgression exists for eastern and western mosquito fish (Gambusia holbrooki and G. affinis) around Mobile Bay, AL. The region has been studied episodically, over decades, making it perfect for addressing temporal dynamics and for providing a deeper understanding of the genetics of these periodically reclassified fishes (as species or subspecies). We used six microsatellite markers to assess the current population structure and gene flow patterns across 19 populations of mosquito fish and then compared our results with historical data. Genetic evidence demonstrates that the current hybrid zone is located in a similar geographic region as the historical one, even after three decades. Hybrid fish, however, demonstrate relatively low heterozygosity and are genetically distinct from western and eastern mosquito fish populations. Fin ray counts, sometimes used to distinguish the two species from one another, demonstrate more eastern (G. holbrooki) phenotype fish within the molecular genetic hybrid zone today. Mosquito fish are globally invasive, often found on the leading edge of flooded waters that they colonize, so the impact of hurricanes in the wake of climate change was also evaluated. An increase in the frequency and intensity of hurricanes in the hybrid region has occurred, and this point warrants further attention since hurricanes are known to move these aggressive, invasive species into novel territory. This work contributes to our classical understanding of hybrid zone temporal dynamics, refines our understanding of mosquito fish genetics in their native range, evaluates important genotype–phenotype relationships, and identifies a potential new impact of climate change.     

Signatures of reproductive isolation in patterns of single nucleotide diversity across inbred strains of mice

Reproductive isolation is often caused by the disruption of genic interactions that evolve in geographically separate populations. Identifying the genomic regions and genes involved in these interactions, known as "Dobzhansky-Muller incompatibilities," can be challenging but is facilitated by the wealth of genetic markers now available in model systems. In recent years, the complete genome sequence and thousands of single nucleotide polymorphisms (SNPs) from laboratory mice, which are largely genetic hybrids between Mus musculus and M. domesticus, have become available. Here, we use these resources to locate genomic regions that may underlie reproductive isolation between these two species. Using genotypes from 332 SNPs that differ between wild-derived strains of M. musculus and M. domesticus, we identified several physically unlinked SNP pairs that show exceptional gametic disequilibrium across the lab strains. Conspecific alleles were associated in a disproportionate number of these cases, consistent with the action of natural selection against hybrid gene combinations. As predicted by the Dobzhansky-Muller model, this bias was differentially attributable to locus pairs for which one hybrid genotype was missing. We assembled a list of potential Dobzhansky-Muller incompatibilities from locus pairs that showed extreme associations (only three gametic types) among conspecific alleles. Two SNPs in this list map near known hybrid sterility loci on chromosome 17 and the X chromosome, allowing us to nominate partners for disrupted interactions involving these genomic regions for the first time. Together, these results indicate that patterns produced by speciation between M. musculus and M. domesticus are visible in the genomes of lab strains of mice, underscoring the potential of these genetic model organisms for addressing general questions in evolutionary biology.     

Testing for selection on the androgen-binding protein in the Danish mouse hybrid zone

The three peripheral subspecies of the house mouse Mus musculus have fixed specific variants of the androgen-binding protein (ABP) that have been proposed to be part of a recognition mechanism that could participate in sexual isolation between the subspecies. We tested for selection on Abpa by characterizing the pattern of Abpa introgression across a transect of the hybrid zone between M. m. musculus and M. m. domesticus in Jutland. On the musculus side, the cline for Abpa resembled that of a nearly neutral allozyme more than that of strongly selected X and Y chromosome markers. However, the high central step which displaces the tail of introgression of Abpa to higher frequencies was best accounted for by linkage to a locus under strong selection against hybrids. Still, we cannot exclude that this pattern results from weak selection against Abpa in the tail of introgression, which would be compatible with an assortative choice mechanism. On the domesticus side there was little introgression close to the hybrid zone, presumably due to a geographical barrier to migration. However, substantial frequencies of musculus alleles occurred further away, suggesting mixed colonization patterns as well as flow across the hybrid zone.  © 2005 The Linnean Society of London, Biological Journal of the Linnean Society , 2005, 84 , 447–459.     

Searching the genomes of inbred mouse strains for incompatibilities that reproductively isolate their wild relatives

Identification of the genes that underlie reproductive isolation provides important insights into the process of speciation. According to the Dobzhansky-Muller model, these genes suffer disrupted interactions in hybrids due to independent divergence in separate populations. In hybrid populations, natural selection acts to remove the deleterious heterospecific combinations that cause these functional disruptions. When selection is strong, this process can maintain multilocus associations, primarily between conspecific alleles, providing a signature that can be used to locate incompatibilities. We applied this logic to populations of house mice that were formed by hybridization involving two species that show partial reproductive isolation, Mus domesticus and Mus musculus. Using molecular markers likely to be informative about species ancestry, we scanned the genomes of 1) classical inbred strains and 2) recombinant inbred lines for pairs of loci that showed extreme linkage disequilibria. By using the same set of markers, we identified a list of locus pairs that displayed similar patterns in both scans. These genomic regions may contain genes that contribute to reproductive isolation between M. domesticus and M. musculus. This hypothesis can now be tested using laboratory crosses and surveys of introgression in the wild.     

Genetic analysis of autosomal and X-linked markers across a mouse hybrid zone

In this paper, we present results of the first comprehensive study of the introgression of both autosomal and sex-chromosome markers across the central European portion of the hybrid zone between two house mouse subspecies, Mus musculus musculus and M. m. domesticus. More than 1800 individuals sampled from 105 sites were analyzed with a set of allozyme loci (hopefully representing neutral or nearly neutral markers) and X-linked loci (which are assumed to be under selection). The zone center is best modeled as a single straight line independent of fine-scale local geographic or climatic conditions, being maintained by a balance between dispersal and selection against hybrids. The width (w) of the multilocus autosomal cline was estimated as 9.6 km whereas the estimate for the compound X-chromosome cline was about 4.6 km only. As the former estimate is comparable to that of the Danish portion of the zone (assumed to be much younger than the central European one), zone width does not appear to be related to its age. The strength (B) of the central barrier was estimated as about 20 km; with dispersal (sigma) of about 1 km/gen(1/2), this means effective selection (s*) is approximately 0.06-0.09 for autosomal loci and about 0.25 for X-linked loci. The number of loci under selection was estimated as N= 56-99 for autosomes and about 380 for X-linked loci. Finally, we highlight some potential pitfalls in hybrid zone analyses and in comparisons of different transects. We suggest that conclusions about parts of the mouse genome involved in reproductive isolation and speciation should be drawn with caution and that analytical approaches always providing some estimates should not be used without due care regarding the support or confidence of such estimates, especially if conclusions are based on the difference between these estimates. Finally, we recommend that analysis in two-dimensional space, dense sampling, and rigorous treatment of data, including inspection of likelihood profiles, are essential for hybrid zone studies.     

A comparison of chromosomal and allozymal variation across a narrow hybrid zone in the grasshopper Caledia captiva

A hybrid zone between the Moreton and Torresian taxa of the grasshopper Caledia captiva in S.E. Queensland has been characterised in terms of allozyme and chromosome variation within the same individuals. — On chromosomal criteria (pericentric rearrangements), the zone is asymmetrical with evidence of high levels of introgression of Torresian chromosomes into the Moreton taxon. This is apparent from the analysis of two independent transects across the hybrid zone. Major changes in chromosomal frequency occur over distances of less than 0.5 km. and the level of introgression differs between the two transects, with much higher levels in the northern Moreton populations, characterised by an acrocentric X-chromosome, when compared with the southern metacentric-X Moreton populations. Chromosome analysis of samples taken from the same transect over two years has revealed no major changes in the structure of the zone. Moreover, a Moreton population located only 0.5 km. from the null point was found to be stable over 6 generations with evidence for a new balanced genome having originated following the differential incorportation of Torresian chromosomes. — Contrary to the chromosomal situation, the same hybrid zone was found to be symmetrical with respect to allozyme variation with evidence of movement of diagnostic alleles in both directions across the zone. The alleles are independent and not tightly linked to any of the pericentric rearrangements. Thus these 5 alleles are acting as markers of the background genome and reveal the relatively free movement of genes which are located outside the pericentric rearrangements. — It is proposed that the hybrid zone in Caledia captiva is unstable and is moving slowly in a westerly direction into the Torresian territory. This is due to the ability of the Moreton taxon to incorporate more readily into its genome those Torresian chromosomes or chromosome segments which increase the fitness of the Moreton taxon. On chromosomal criteria, the Torresian taxon does not share the same capacity. — It is suggested that, so long as the two taxa retain their ability to hybridise with subsequent asymmetrical introgression, the zone will continue to move westwards and eventually lead to the selective incorporation of the Torresian genome into the Moreton taxon. This will result in a polymorphic situation with clinal variation in chromosomal frequencies. The structure of the zone is dependent upon a fine balance between genomic reorganisation in recombinant genotypes and the relative dispersal capacities of the two hybridising taxa.     

Inferences of selection and migration in the Danish house mouse hybrid zone

We analysed the patterns of allele frequency change for ten diagnostic autosomal allozyme loci in the hybrid zone between the house mouse subspecies Mus musculus domesticus and M. m. musculus in central Jutland. After determining the general orientation of the clines of allele frequencies, we analysed the cline shapes along the direction of maximum gradient. Eight of the ten clines are best described by steep central steps with coincident positions and an average width of 8.9 km (support limits 7.6–12.4) flanked by tails of introgression, indicating the existence of a barrier to gene flow and only weak selection on the loci studied. We derived estimates of migration from linkage disequilibrium in the centre of the zone, and by applying isolation by distance methods to microsatellite data from some of these populations. These give concordant estimates of σ =  0.5–0.8 km generation     . The barrier to gene flow is of the order of 20 km (support limits 14–28), and could be explained by selection of a few per cent at 43–120 underdominant loci that reduces the mean fitness in the central populations to 0.45. Some of the clines appear symmetrical, whereas others are strongly asymmetrical, and two loci appear to have escaped the central barrier to gene flow, reflecting the differential action of selection on different parts of the genome. Asymmetry is always in the direction of more introgression into musculus , indicating either a general progression of domesticus into the musculus territory, possibly mediated by differential behaviour, or past movement of the hybrid zone in the opposite direction, impeded by potential geographical barriers to migration in domesticus territory.  © 2005 The Linnean Society of London, Biological Journal of the Linnean Society , 2005, 84 , 593–616.     

Cross‐decades stability of an avian hybrid zone

Hybrid zones are particularly valuable for understanding the evolution of partial reproductive isolation between differentiated populations. An increasing number of hybrid zones have been inferred to move over time, but in most such cases zone movement has not been tested with long‐term genomic data. The hybrid zone between Townsend's Warblers (Setophaga townsendi) and Hermit Warblers (S. occidentalis) in the Washington Cascades was previously inferred to be moving from northern S. townsendi southwards towards S. occidentalis, based on plumage and behavioural patterns as well as a 2000‐km genetic wake of hermit mitochondrial DNA (mtDNA) in coastal Townsend's Warblers. We directly tested whether hybrid zone position has changed over 2–3 decades by tracking plumage, mtDNA and nuclear genomic variation across the hybrid zone over two sampling periods (1987–94 and 2015–16). Surprisingly, there was no significant movement in genomic or plumage cline centres between the two time periods. Plumage cline widths were narrower than expected by neutral diffusion, consistent with a ‘tension zone’ model, in which selection against hybrids is balanced by movement of parental forms into the zone. Our results indicate that this hybrid zone is either stable in its location or moving at a rate that is not detectable over 2–3 decades. Despite considerable gene flow, the stable clines in multiple phenotypic and genotypic characters over decades suggest evolutionary stability of this young pair of sister species, allowing divergence to continue. We propose a novel biogeographic scenario to explain these patterns: rather than the hybrid zone having moved thousands of kilometres to its current position, inland Townsend's met coastal Hermit Warbler populations along a broad front of the British Columbia and Alaska coast and hybridization led to replacement of the Hermit Warbler plumage with Townsend's Warbler plumage patterns along this coastline. Hence, hybrid zones along British Columbia and Alaska moved only a short distance from the inland to the coast, whereas the Hermit Warbler phenotype appears stable in Washington and further south. This case provides an example of the complex biogeographic processes that have led to the distribution of current phenotypes within and among closely related species.     

Diverse wild mouse origins of xenotropic, mink cell focus-forming, and two types of ecotropic proviral genes

We analyzed wild mouse DNAs for the number and type of proviral genes related to the env sequences of various murine leukemia viruses (MuLVs). Only Mus species closely related to laboratory mice carried these retroviral sequences, and the different subclasses of viral env genes tended to be restricted to specific taxonomic groups. Only Mus musculus molossinus carried proviral genes which cross-reacted with the inbred mouse ecotropic MuLV env gene. The ecotropic viral env sequence associated with the Fv-4 resistance gene was found in the Asian mice M. musculus molossinus and Mus musculus castaneus and in California mice from Lake Casitas (LC). Both M. musculus castaneus and LC mice carried many additional Fv-4 env-related proviruses, two of which are common to both mouse populations, which suggests that these mice share a recent common ancestry. Xenotropic and mink cell focus-forming (MCF) virus env sequences were more widely dispersed in wild mice than the ecotropic viral env genes, which suggests that nonecotropic MuLVs were integrated into the Mus germ line at an earlier date. Xenotropic MuLVs represented the major component of MuLV env-reactive genes in Asian and eastern European mice classified as M. musculus molossinus, M. musculus castaneus, and Mus musculus musculus, whereas Mus musculus domesticus from western Europe, the Mediterranean, and North America contained almost exclusively MCF virus env copies. M. musculus musculus mice from central Europe trapped near the M. musculus domesticus/M. musculus musculus hybrid zone carried multiple copies of both types of env genes. LC mice also carried both xenotropic and MCF viral env genes, which is consistent with the above conclusion that they represent natural hybrids of M. musculus domesticus and M. musculus castaneus.     

Relative warp analysis of skull shape across the hybrid zone of the house mouse (Mus musculus) in Denmark.

The changes in skull shape that occur in the hybrid zone between the two European subspecies of the house mouse ( Mus musculus musculus and M. m. domesticus ) were studied by relative warp analyses. Landmarks observed on the ventral view of the skulls of 269 mice sampled in 18 localities in Denmark were analysed. Each population was also characterized by a hybrid index estimated from allozymic data. Although no clear pattern of within-population variability of shape could be established across the hybrid zone, the shape changes among the 18 populations were found to be correlated with the allozymic introgression but not with geographical location. Finally, the cline obtained for skull shape seems to be steeper than that corresponding to the average allozyme hybrid index which suggests that skull development could be slightly perturbed in hybrids.     

Quantifying the spatiotemporal dynamics in a chorus frog (Pseudacris) hybrid zone over 30 years

Although theory suggests that hybrid zones can move or change structure over time, studies supported by direct empirical evidence for these changes are relatively limited. We present a spatiotemporal genetic study of a hybrid zone between Pseudacris nigrita and P. fouquettei across the Pearl River between Louisiana and Mississippi. This hybrid zone was initially characterized in 1980 as a narrow and steep “tension zone,” in which hybrid populations were inferior to parentals and were maintained through a balance between selection and dispersal. We reanalyzed historical tissue samples and compared them to samples of recently collected individuals using microsatellites. Clinal analyses indicate that the cline has not shifted in roughly 30 years but has widened significantly. Anthropogenic and natural changes may have affected selective pressure or dispersal, and our results suggest that the zone may no longer best be described as a tension zone. To the best of our knowledge, this study provides the first evidence of significant widening of a hybrid cline but stasis of its center. Continued empirical study of dynamic hybrid zones will provide insight into the forces shaping their structure and the evolutionary potential they possess for the elimination or generation of species.     

Admixture in European Populus hybrid zones makes feasible the mapping of loci that contribute to reproductive isolation and trait differences

The use of admixed human populations to scan the genome for chromosomal segments affecting complex phenotypic traits has proved a powerful analytical tool. However, its potential in other organisms has not yet been evaluated. Here, we use DNA microsatellites to assess the feasibility of this approach in hybrid zones between two members of the 'model tree' genus Populus: Populus alba (white poplar) and Populus tremula (European aspen). We analyzed samples of both species and a Central European hybrid zone (N=544 chromosomes) for a genome-wide set of 19 polymorphic DNA microsatellites. Our results indicate that allele frequency differentials between the two species are substantial (mean delta=0.619+/-0.067). Background linkage disequilibrium (LD) in samples of the parental gene pools is moderate and should respond to sampling schemes that minimize drift and account for rare alleles. LD in hybrids decays with increasing number of backcross generations as expected from theory and approaches background levels of the parental gene pools in advanced generation backcrosses. Introgression from P. tremula into P. alba varies strongly across marker loci. For several markers, alleles from P. tremula are slightly over-represented relative to neutral expectations, whereas a single locus exhibits evidence of selection against P. tremula genotypes. We interpret our results in terms of the potential for admixture mapping in these two ecologically divergent Populus species, and we validate a modified approach of studying genotypic clines in 'mosaic' hybrid zones.     

The diverse origins of New Zealand house mice

Molecular markers and morphological characters can help infer the colonization history of organisms. A combination of mitochondrial (mt) D-loop DNA sequences, nuclear DNA data, external measurements and skull characteristics shows that house mice (Mus musculus) in New Zealand and its outlying islands are descended from very diverse sources. The predominant genome is Mus musculus domesticus (from western Europe), but Mus musculus musculus (from central Europe) and Mus musculus castaneus (from southern Asia) are also represented genetically. These subspecies have hybridized to produce combinations of musculus and domesticus nuclear DNA coupled with domesticus mtDNA, and castaneus or musculus mtDNA with domesticus nuclear DNA. The majority of the mice with domesticus mtDNA that we sampled had D-loop sequences identical to two haplotypes common in Britain. This is consistent with long-term British-New Zealand cultural linkages. The origins of the castaneus mtDNA sequences widespread in New Zealand are less easy to identify.