Mitochondrial DNA in the hybrid zone between Mus musculus musculus and Mus musculus domesticus: a comparison of two transects

We studied mtDNA introgression across the contact zone between Mus musculus musculus and M. m. domesticus in two independent transects in the Czech Republic and Bavaria, Germany. A total of 1270 mice from 98 localities in the Czech transect and 456 mice from 41 localities in the Bavarian transect were examined for presence or absence of a Bam HI restriction site in the mt-Nd1 gene. Using this simple mtDNA marker, variants that belonged to the M. m. domesticus lineage (presence of restriction site) could be unequivocally distinguished from those belonging to the M. m. musculus lineage (absence of restriction site). The extent of introgression of mtDNA, three autosomal allozymes and the X chromosome was compared. The introgression of X markers was more limited than was that of the allozymes and mtDNA. In the Czech transect, the centre for the mtDNA cline was shifted about 3.6 km to the west relative to the X chromosome cline, with asymmetric introgression from M. m. musculus to M. m. domesticus . Interestingly, in the Bavarian transect, the centre of the mtDNA cline was shifted about 10.9 km to the east relative to the X chromosome cline, with asymmetric introgression from M. m. domesticus to M. m. musculus, opposite in direction to that observed in the Czech transect.  © 2005 The Linnean Society of London, Biological Journal of the Linnean Society , 2005, 84 , 363–378.     

Higher differentiation among subspecies of the house mouse (Mus musculus) in genomic regions with low recombination

In the early stages of reproductive isolation, genomic regions of reduced recombination are expected to show greater levels of differentiation, either because gene flow between species is reduced in these regions or because the effects of selection at linked sites within species are enhanced in these regions. Here, we study the patterns of DNA sequence variation at 27 autosomal loci among populations of Mus musculus musculus, M. m. domesticus, and M. m. castaneus, three subspecies of house mice with collinear genomes. We found that some loci exhibit considerable shared variation among subspecies, while others exhibit fixed differences. We used an isolation-with-gene-flow model to estimate divergence times and effective population sizes (N(e) ) and to disentangle ancestral variation from gene flow. Estimates of divergence time indicate that all three subspecies diverged from one another within a very short period of time approximately 350,000 years ago. Overall, N(e) for each subspecies was associated with the degree of genetic differentiation: M. m. musculus had the smallest N(e) and the greatest proportion of monophyletic gene genealogies, while M. m. castaneus had the largest N(e) and the smallest proportion of monophyletic gene genealogies. M. m. domesticus and M. m. musculus were more differentiated from each other than either were from M. m. castaneus, consistent with greater reproductive isolation between M. m. domesticus and M. m. musculus. F(ST) was significantly greater at loci experiencing low recombination rates compared to loci experiencing high recombination rates in comparisons between M. m. castaneus and M. m. musculus or M. m. domesticus. These results provide evidence that genomic regions with less recombination show greater differentiation, even in the absence of chromosomal rearrangements.     

Mosaic genealogy of the Mus musculus genome revealed by 21 nuclear genes from its three subspecies

Patterns of genetic variation provide insight into the evolutionary history of a species. Mouse (Mus musculus) is a good model for this purpose. Here we present the analysis of genealogies of the 21 nuclear loci and one mitochondrial DNA region in M. musculus based on our nucleotide sequences of nine inbred strains from three M. musculus subspecies (musculus, domesticus, and castaneus) and one M. spicilegus strain as an outgroup. The mitochondrial DNA gene genealogy of those strains confirmed the introgression pattern of one musculus strain. When all the nuclear DNA data were concatenated to produce a phylogenetic tree of nine strains, musculus and domesticus strains formed monophyletic clusters with each other, while the two castaneus strains were paraphyletic. When each DNA region was treated independently, the phylogenetic networks revealed an unnegligibly high level of subspecies admixture and the mosaic nature of their genome. Estimation of ancestral and derived population sizes and migration rates suggests the effects of ancestral polymorphism and gene flow on the pattern of genetic variation of the current subspecies. Gene genealogies of Fut4 and Dfy loci also suggested existence of the gene flow between M. musculus and M. spicilegus or other distant species.     

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. Taxon-specific 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. wagnery 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. abbottii). The estimates of genetic similarity and the taxonomic 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.     

The mitochondrial genome sequence of Mus terricolor: comparison with Mus musculus domesticus and implications for xenomitochondrial mouse modeling

Knowledge of the mitochondrial DNA (mtDNA) sequence of divergent murine species is critical from both a phylogenetic perspective and in understanding nuclear-mitochondrial interactions, particularly as the latter influences our xenocybrid models of mitochondrial disease. To this end, the sequence of the mitochondrial genome of the murine species Mus terricolor (formerly Mus dunni) is reported and compared with the published sequence for the common laboratory mouse Mus musculus domesticus strain C57BL/6J. These species are of interest because xenomitochondrial cybrid mice were created that harbor M. terricolor mtDNA in a M. m. domesticus nuclear background. Although the total of 1763 nucleotide substitutions represents striking heterogeneity, the majority of these are silent, leading to highly conserved protein sequences with only 159 amino acid differences. Moreover, 58% of these amino acid differences represented conservative substitutions. All of the tRNA genes and rRNA genes have homology of 91% or greater. The control region shows the greatest heterogeneity, as expected, with 85% homology overall. Regions of 100% homology were found for Conserved Sequence Block I, Conserved Sequence Block III and the L-strand origin of replication. Complex I genes showed the greatest degree of difference among protein-coding genes with amino acid homology of 91-97% among the seven mitochondrial genes. Complexes III and IV genes show high homology ranging from 98-100%. From these data, complex I differences appear most critical for the viability of M. m. domesticus: M. terricolor cybrids. Moreover, the sequence information reported here should be useful in identifying critical regions for mitochondrial transfer between species, for furthering the understanding of mitochondrial dynamics and pathology in transmitochondrial organisms, and for the study of Mus genus origins.     

Genome Patterns of Selection and Introgression of Haplotypes in Natural Populations of the House Mouse (Mus musculus)

General parameters of selection, such as the frequency and strength of positive selection in natural populations or the role of introgression, are still insufficiently understood. The house mouse (Mus musculus) is a particularly well-suited model system to approach such questions, since it has a defined history of splits into subspecies and populations and since extensive genome information is available. We have used high-density single-nucleotide polymorphism (SNP) typing arrays to assess genomic patterns of positive selection and introgression of alleles in two natural populations of each of the subspecies M. m. domesticus and M. m. musculus. Applying different statistical procedures, we find a large number of regions subject to apparent selective sweeps, indicating frequent positive selection on rare alleles or novel mutations. Genes in the regions include well-studied imprinted loci (e.g. Plagl1/Zac1), homologues of human genes involved in adaptations (e.g. alpha-amylase genes) or in genetic diseases (e.g. Huntingtin and Parkin). Haplotype matching between the two subspecies reveals a large number of haplotypes that show patterns of introgression from specific populations of the respective other subspecies, with at least 10% of the genome being affected by partial or full introgression. Using neutral simulations for comparison, we find that the size and the fraction of introgressed haplotypes are not compatible with a pure migration or incomplete lineage sorting model. Hence, it appears that introgressed haplotypes can rise in frequency due to positive selection and thus can contribute to the adaptive genomic landscape of natural populations. Our data support the notion that natural genomes are subject to complex adaptive processes, including the introgression of haplotypes from other differentiated populations or species at a larger scale than previously assumed for animals. This implies that some of the admixture found in inbred strains of mice may also have a natural origin.     

Hybrid origin of Japanese mice "Mus musculus molossinus": evidence from restriction analysis of mitochondrial DNA

The Japanese mouse, Mus musculus molossinus, has long been considered an independent subspecies of the house mouse. A survey of restriction- site haplotypes of mitochondrial DNA (mtDNA) showed that Japanese mice have two main maternal lineages. The most common haplotype is closely related to the mtDNA of the European subspecies M. m. musculus. The other common haplotype and two minor ones are closely related to each other and to the mtDNA of an Asiatic subspecies, M. m. castaneus. Two other rare variants are probably the result of recent contamination by European M. m. domesticus. The musculus type of mtDNA is found in the southern two-thirds of Japan, whereas the common castaneus type is found in the northern third and the minor variants are found sporadically throughout Japan. The castaneus mtDNA lineage had a few minor variants, whereas the musculus lineage was completely monomorphic. By contrast, the native population of M. m. castaneus and the Chinese and Korean musculus populations were highly polymorphic. These results suggest that M. m. molossinus is a hybrid between ancestral colonies, possibly very small, of M. m. musculus and M. m. castaneus, rather than an independent subspecies.      

Natural hybridization between 2 sympatric species of mice, Mus musculus domesticus L. and Mus spretus Lataste

Using protein loci and DNA markers, we show by a multilocus genetic analysis that certain populations of the two sympatric mouse species Mus musculus domesticus and Mus spretus show clear signs of partial introgression. Given the sterility of F1 males and the known partial genetic incompatibilities between the genomes of the two species, our finding does not invalidate the biological species complex, but allows to think that very limited genetic exchanges remain possible even long after the divergence of taxa. This may have some consequences on the dynamics of certain kinds of invasive or advantageous DNAs like transposable elements or pathogen resistance genes.     

Reduced nucleotide variability at an androgen-binding protein locus (Abpa) in house mice: evidence for positive natural selection.

Previous work has shown that the gene for the alpha subunit of androgen-binding protein, Abpa, may be involved in premating isolation between different subspecies of the house mouse, Mus musculus. We investigated patterns of DNA sequence variation at Abpa within and between species of mice to test several predictions of a model of neutral molecular evolution. Intraspecific variation among 10 Mus musculus domesticus alleles was compared with divergence between M. m. domesticus and M. caroli for Abpa and two X-linked genes, Glra2 and Amg. No variation was observed at Abpa within M. m. domesticus. The ratio of polymorphism to divergence was significantly lower at Abpa than at Glra2 and Amg, despite the fact that all three genes experience similar rates of recombination. Interspecific comparisons among M. m. domesticus, Mus musculus musculus, Mus musculus castaneus, Mus spretus, Mus spicilegus, and Mus caroli revealed that the ratio of nonsynonymous substitutions to synonymous substitutions on a per-site basis (Ka/Ks) was generally greater than one. The combined observations of no variation at Abpa within M. m: domesticus and uniformly high Ka/Ks values between species suggest that positive directional selection has acted recently at this locus.     

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.     

The role of salivary androgen-binding protein in reproductive isolation between two subspecies of house mouse: Mus musculus musculus and Mus musculus domesticus

Previous behavioural studies using inbred lines have suggested that the gene ( Abpa ) for the alpha subunit of salivary androgen-binding protein (ABP) plays a role in prezygotic isolation between house mouse Mus musculus subspecies. We tested this hypothesis in animals from wild allopatric (121 individuals from four samples) and parapatric (320 animals from 15 samples) populations sampled on the Czech–Bavarian transect across the hybrid zone between M. m. domesticus and M. m. musculus . The study did not reveal a consistent statistically significant trend of homosubspecific preferences in individual allopatric and parapatric populations. Nonetheless, the whole pattern of preference was skewed toward homosubspecific preference mostly on the M. m. musculus side of the hybrid zone. The pattern of homosubspecific preferences was stronger for the time spent sniffing than it was for the first choice of the signal (the ratio of homosubspecific vs. heterosubspecific preferences for both sexes was 6 : 2 in allopatric and 21 : 9 in parapatric populations, while the same rates were 4 : 4 and 16 : 14 for the first choice). To the extent that Y-maze tests reflect preference under wild conditions, we suggest that this slight preference may not in itself be sufficient to impede gene flow between the two subspecies and thus act as a reproductive barrier. ABP most probably participates in a complex system of subspecies-specific recognition in the hybrid zone, but the picture is far too complex at this time to allow a conclusive evaluation of the importance of this role.  © 2005 The Linnean Society of London, Biological Journal of the Linnean Society , 2005, 84 , 349–361.     

The genomics of speciation: investigating the molecular correlates of X chromosome introgression across the hybrid zone between Mus domesticus and Mus musculus

Understanding the genetic details of reproductive isolation is a key goal in the study of speciation. Hybrid zones, geographical regions where two species meet and exchange genes, can provide insight into the genetic basis of reproductive isolation. This is especially true in species with mapped molecular markers because patterns of gene flow can be compared among different genomic regions. Even greater insight can be obtained in species with complete genome sequences because gene identity, gene number and other features of interest can be assessed for genomic regions with different patterns of introgression. Here, we review recent studies on the well-characterized hybrid zone between Mus domesticus and M. musculus , including a detailed survey of patterns of introgression for 13 markers on the X chromosome. We then compare levels of introgression for these 13 regions to a number of genomic attributes inferred from the complete sequence of the X chromosome, with two purposes. First, we identify candidate genes for reproductive isolation by finding genes that map to an X-linked region of reduced introgression and that are only expressed in the male germ line or that show high rates of protein evolution in comparison with rat. Second, we ask whether patterns of gene flow are correlated with recombination rate, gene density, base composition, CpG island density, mutation rate and the rate of protein evolution, as might be expected if many genes contribute to reproductive isolation. We identify seven candidate genes for reproductive isolation between M. domesticus and M. musculus , and our analyses reveal no general correlations between levels of introgression and other measured sequence characteristics. These results underline the utility of the house mouse as a model system for the study of speciation.  © 2005 The Linnean Society of London, Biological Journal of the Linnean Society , 2005, 84 , 523–534.     

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.     

Two types of mouse (Mus musculus domesticus) Y chromosomes in Quebec.

A Y chromosomal repetitive sequence identified two types of Y chromosomes in mice (Mus musculus domesticus) caught near Ste. Anne de Bellevue, Quebec. One type is apparently identical to the Y chromosome found in Maryland, Delaware, and California, whereas the other type is similar, but not identical, to the Y chromosome present in M.m. poschiavinus, an Alpine race of M.m. domesticus. These findings suggest that the domesticus Y chromosome is highly polymorphic and thus useful for elucidating the relationships among American and European house mouse populations.     

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.     

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.     

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.     

Frequency and distribution of t-haplotypes in the Southeast Asian house mouse (Mus musculus castaneus) in Taiwan

t-haplotypes are a meiotic drive system found on the 17th chromosome of the house mouse (Mus musculus). They can be found in wild populations of all four genetically differentiated subspecies. The drive phenomenon is male-specific, such that heterozygous males (+/t) show non-Mendelian transmission and transmit the t-chromosome to > 90% of their offspring. So far the most comprehensive studies on the frequencies of t-haplotypes in natural populations have been on just one of the subspecies (M. musculus domesticus). We applied molecular methods to accurately screen t-haplotypes in a large number of populations of a second subspecies (M. musculus castaneus) distributed in Taiwan. We found that the overall t-haplotype frequency is low in M. m. castaneus (0.108), and the chromosomes are patchily distributed among its populations, closely resembling the situation found in M. m. domesticus. Further, we found the frequencies of t-haplotypes in our sample did not differ in relation to the sex or age of mice. This resemblance in the frequency and distribution among populations of the two distinct subspecies suggests that similar general mechanisms might be responsible for the low frequencies in both subspecies.     

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.