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.     

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.     

A pronounced evolutionary shift of the pseudoautosomal region boundary in house mice

The pseudoautosomal region (PAR) is essential for the accurate pairing and segregation of the X and Y chromosomes during meiosis. Despite its functional significance, the PAR shows substantial evolutionary divergence in structure and sequence between mammalian species. An instructive example of PAR evolution is the house mouse Mus musculus domesticus (represented by the C57BL/6J strain), which has the smallest PAR among those that have been mapped. In C57BL/6J, the PAR boundary is located just ~700 kb from the distal end of the X chromosome, whereas the boundary is found at a more proximal position in Mus spretus, a species that diverged from house mice 2-4 million years ago. In this study we used a combination of genetic and physical mapping to document a pronounced shift in the PAR boundary in a second house mouse subspecies, Mus musculus castaneus (represented by the CAST/EiJ strain), ~430 kb proximal of the M. m. domesticus boundary. We demonstrate molecular evolutionary consequences of this shift, including a marked lineage-specific increase in sequence divergence within Mid1, a gene that resides entirely within the M. m. castaneus PAR but straddles the boundary in other subspecies. Our results extend observations of structural divergence in the PAR to closely related subspecies, pointing to major evolutionary changes in this functionally important genomic region over a short time period.     

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.     

Adaptive evolution and effective population size in wild house mice

Estimates of the proportion of amino acid substitutions that have been fixed by selection (α) vary widely among taxa, ranging from zero in humans to over 50% in Drosophila. This wide range may reflect differences in the efficacy of selection due to differences in the effective population size (N(e)). However, most comparisons have been made among distantly related organisms that differ not only in N(e) but also in many other aspects of their biology. Here, we estimate α in three closely related lineages of house mice that have a similar ecology but differ widely in N(e): Mus musculus musculus (N(e) ～ 25,000-120,000), M. m. domesticus (N(e) ～ 58,000-200,000), and M. m. castaneus (N(e) ～ 200,000-733,000). Mice were genotyped using a high-density single nucleotide polymorphism array, and the proportions of replacement and silent mutations within subspecies were compared with those fixed between each subspecies and an outgroup, Mus spretus. There was significant evidence of positive selection in M. m. castaneus, the lineage with the largest N(e), with α estimated to be approximately 40%. In contrast, estimates of α for M. m. domesticus (α = 13%) and for M. m. musculus (α = 12 %) were much smaller. Interestingly, the higher estimate of α for M. m. castaneus appears to reflect not only more adaptive fixations but also more effective purifying selection. These results support the hypothesis that differences in N(e) contribute to differences among species in the efficacy of selection.     

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.      

Nucleotide variation in wild and inbred mice

The house mouse is a well-established model organism, particularly for studying the genetics of complex traits. However, most studies of mice use classical inbred strains, whose genomes derive from multiple species. Relatively little is known about the distribution of genetic variation among these species or how variation among strains relates to variation in the wild. We sequenced intronic regions of five X-linked loci in large samples of wild Mus domesticus and M. musculus, and we found low levels of nucleotide diversity in both species. We compared these data to published data from short portions of six X-linked and 18 autosomal loci in wild mice. We estimate that M. domesticus and M. musculus diverged <500,000 years ago. Consistent with this recent divergence, some gene genealogies were reciprocally monophyletic between these species, while others were paraphyletic or polyphyletic. In general, the X chromosome was more differentiated than the autosomes. We resequenced classical inbred strains for all 29 loci and found that inbred strains contain only a small amount of the genetic variation seen in wild mice. Notably, the X chromosome contains proportionately less variation among inbred strains than do the autosomes. Moreover, variation among inbred strains derives from differences between species as well as from differences within species, and these proportions differ in different genomic regions. Wild mice thus provide a reservoir of additional genetic variation that may be useful for mapping studies. Together these results suggest that wild mice will be a valuable complement to laboratory strains for studying the genetics of complex traits.     

Genetic studies of house mice from the hybrid area of the Primorskii Region

Cytogenetic and allozyme differentiation of house mice Mus musculus L. was studied in 21 samples of Primorskii krai (region) from 1984 to 1998. The use of Q-H- and C-banding made it possible to reveal a high variation in the content and distribution of pericentromeric heterochromatin and high frequencies of marker autosomes 5, 7, 10, 12, 15, and 17-19 in the karyotypes of mice from some samples. The presence of marker alleles that are typical of the subspecies groups castaneus and musculus was shown for protein loci Idh, Sod-1, Aat-1, and Hbb. The data obtained indicate that the fauna of house mice of Primorye is characterized by subspecific mosaicism and hybrid origin. Apparently, subspecies M. m. castaneus, M. musculus ssp., M. M. homourus, and M. m. gansuensis contributed, to different degrees, to the formation of this fauna.     

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.     

Genetic diversity of the house mouse Mus musculus and geographic distribution of its subspecies-specific RAPD markers on the territory of Russia

Genetic diversity and geographic distribution of taxon-specific RAPD markers was examined in ten local populations of the house mouse Mus musculus (n = 42). The house mice were generally characterized by moderate genetic variation: polymorphism P99 = 60%, P95 = 32.57%; heterozygosity H = 0.12; the observed allele number n(a) = 1.6; the effective allele number n(e) = 1.18; the within-population differentiation Theta = 0.388; and Shannon index I = 0.19. The degree of genetic isolation of individual local populations was greatly variable. The genetic subdivision index G(st) varied from 0.162 to 0.770 at the gene flow of Nm = 2.58-0.149, while the among-population distances D(N) varied from 0.026 to 0.178. of the largest part of the genetic diversity was found among the populations (H(T) = 0.125), while the within-population diversity was twice lower (H(S) = 0.06). The samples examined were well discriminated relative to the sets of RAPD markers. The character distribution pattern provided conditional subdivision of the mice into the "western" and the "eastern" groups with the putative boarder along the Baikal Lake. The first group was characterized by the prevalence of the markers typical of M. m. musculus and M. m. domesticus. The second group was characterized by the prevalence of the markers typical of M. m. musculus, M. m. gansuensis, M. m. castaneus, M. m. domesticus, and m. m. wagneri. The genotype of the nominative subspecies M. m. musculus was background for all populations. In the populations examined some of earlier described subspecies-specific molecular markers were found at different frequencies, pointing to the involvement of several subspecies of M. musculus in the process of hybridization.     

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.     

Restriction fragment length variations and chromosome mapping of two mouse metallothionein genes,Mt-1 andMt-2

Restriction endonuclease fragment length variations (RFLVs) were found through the use of cDNA probes for metallothionein genes 1 (Mt-1) and 2 (Mt-2) in the mouse. RFLVs were detected in restriction patterns generated by BglII and XbaI in the Mt-1 gene and by PvuII in the Mt-2 gene. All laboratory strains carry the Mt-1a and Mt-2a alleles. Among strains of wild origin, some Western European subspecies (Mus mus domesticus and M. m. brevirostris) also carry the Mt-1a and Mt-2a alleles. In contrast, a European subspecies (M. m. musculus) and the great majority of subspecies from East Asian countries (M. m. molossinus, Chinese mice of wild origin, and M. m. yamashinai) carry the Mt-1b and Mt-2b alleles. A domesticus strain from Bulgaria and two castaneus strains from Thailand and Philippines carry the intermediate combination of Mt-1b and Mt-2a alleles. Using the RFLVs, we mapped the Mt-1 and Mt-2 genes on chromosome 8, and they appear to be very closely linked since no recombination was observed between them in any of the mice examined. Data from three-point cross tests showed that the recombination frequencies are 4.31% between Os and Mt, 15.52% between Mt and Prt-2, and 19.83% between Os and Prt-2. The gene order of Os-Mt-1,Mt-2-Prt-2 has been confirmed.     

Genetic analysis of the hybridization zone between two subspecies Mus musculus domesticus and Mus musculus musculus in Bulgaria

The hybrid zone between the two subspecies of mice Mus musculus domesticus and Mus musculus musculus, which has been studied extensively in Denmark, crosses Europe to the Black Sea through the Alps and the Balkans. Two hundred and seventy-nine animals were captured in 22 localities along a transect across the Balkans. The animals were characterized for seven diagnostic nuclear loci by protein electrophoresis and by restriction pattern analysis of their mitochondrial DNA. The nuclear data show a sharp transition between the two subspecies, most of the variations in allele frequencies (from 0.9 to 0.1) occurring within a 36-km section of the transect. The introgression varies from one locus to the other and is more pronounced, in terms of distance, in M. m. musculus territory. Mitochondrial DNA introgression is important but occurs in one direction only, i.e. from M. m. musculus to M. m. domesticus, while a cytoplasmic transfer from M. m. domesticus to M. m. musculus has been reported. A previous study showed that no Y chromosome introgression occurs. The different behaviour of these three types of markers could be due to the interaction between selection against hybrid genomes and meiotic recombination. Objectively, it would appear that the genes that can introgress are neutral or nearly so and have been separated from deleterious genes they were linked to by recombination. This could explain the differential introgression between autosomal loci. The mitochondrial and Y chromosomes undergo no or very little recombination and each is transmitted as a whole. Their degree of introgression is thus indicative of the intensity of selection resulting from the amount of functional differentiation between the two taxa, which seems to be strong for the Y chromosome and weak for mitochondrial DNA. We propose that the asymmetry of nuclear introgression is due to different population structures. As M. m. musculus is relatively less structured, the rapid spreading of introgressed genes would be favoured. Such a scheme, however, can hardly account for the unidirectionality of the mitochondrial flow, which could be due to sex-dependent behaviour.     

Isolation and gene flow: inferring the speciation history of European house mice

Inferring the history of isolation and gene flow during species differentiation can inform us on the processes underlying their formation. Following their recent expansion in Europe, two subspecies of the house mouse (Mus musculus domesticus and Mus musculus musculus) have formed a hybrid zone maintained by hybrid incompatibilities and possibly behavioural reinforcement, offering a good model of incipient speciation. We reconstruct the history of their divergence using an approximate Bayesian computation framework and sequence variation at 57 autosomal loci. We find support for a long isolation period preceding the advent of gene flow around 200,000 generations ago, much before the formation of the European hybrid zone a few thousand years ago. The duration of the allopatric episode appears long enough (74% of divergence time) to explain the accumulation of many post-zygotic incompatibilities expressed in the present hybrid zone. The ancient contact inferred could have played a role in mating behaviour divergence and laid the ground for further reinforcement. We suggest that both subspecies originally colonized the Middle East from the northern Indian subcontinent, domesticus settling on the shores of the Persian Gulf and musculus on those of the Caspian Sea. Range expansions during interglacials would have induced secondary contacts, presumably in Iran, where they must have also interacted with Mus musculus castaneus. Future studies should incorporate this possibility, and we point to Iran and its surroundings as a hot spot for house mouse diversity and speciation studies.     

Mitochondrial DNA Sequence Variation in the Eastern House Mouse, Mus Musculus: Comparison with Other House Mice and Report of a 75-Bp Tandem Repeat

The control region and flanking tRNAs were sequenced from 139 Mus musculus mitochondrial DNAs (mtDNAs) from mice collected at 44 localities extending from Germany to Japan. Among the 36 types of M. musculus mtDNA resolved, five have an added 75-bp direct repeat; the two copies within an individual differ by two to four base substitutions. Among 90 M. domesticus mtDNAs sequenced, 12 new types were found; 96 M. domesticus types have now been identified by sequencing this segment. Representative mtDNAs from M. castaneus, M. macedonicus, M. spicilegus and M. spretus were also sequenced. A parsimony tree for the M. musculus mtDNAs is about half as deep as the tree for the M. domesticus mtDNAs, which is consistent with the idea that M. musculus is genetically less diverse and younger than M. domesticus. The patterns of variation as a function of position are similar but not identical in M. musculus and M. domesticus mtDNAs. M. castaneus and M. musculus mtDNAs are allied, at a tree depth about three times as great as the start of intra-M. musculus divergence. The coalescence of the M. musculus and M. castaneus mtDNAs is about half as deep as their coalescence with the M. domesticus mtDNA lineages. The mtDNAs of the aboriginal M. macedonicus and M. spicilegus are each other's closest relatives, at a tree depth greater than the deepest intracommensal node. The mtDNA results support the view that the aboriginal M. spretus is the sister group of the other five species.     

Genetics and evolution of hybrid male sterility in house mice

Comparative genetic mapping provides insights into the evolution of the reproductive barriers that separate closely related species. This approach has been used to document the accumulation of reproductive incompatibilities over time, but has only been applied to a few taxa. House mice offer a powerful system to reconstruct the evolution of reproductive isolation between multiple subspecies pairs. However, studies of the primary reproductive barrier in house mice-hybrid male sterility-have been restricted to a single subspecies pair: Mus musculus musculus and Mus musculus domesticus. To provide a more complete characterization of reproductive isolation in house mice, we conducted an F(2) intercross between wild-derived inbred strains from Mus musculus castaneus and M. m. domesticus. We identified autosomal and X-linked QTL associated with a range of hybrid male sterility phenotypes, including testis weight, sperm density, and sperm morphology. The pseudoautosomal region (PAR) was strongly associated with hybrid sterility phenotypes when heterozygous. We compared QTL found in this cross with QTL identified in a previous F(2) intercross between M. m. musculus and M. m. domesticus and found three shared autosomal QTL. Most QTL were not shared, demonstrating that the genetic basis of hybrid male sterility largely differs between these closely related subspecies pairs. These results lay the groundwork for identifying genes responsible for the early stages of speciation in house mice.     

Mouse subspecies differentiation and H-2 polymorphism

Four major groups of Mus musculus subspecies are analysed from the genetical view point. A number of genetic traits exhibited polymorphism. Most were common among domesticus, bactrianus and castaneus subspecies, while others were musculus subspecies-specific. This suggests that the earliest differentiation was between domesticus and musculus . Mouse class I MHC was compared in these two lines. Immunological and molecular analyses of H-2K molecules indicated that these two groups of mouse have several common haplotypes of class I molecules. Phenotype frequencies of these haplotypes are considerably different in each subspecies. "Overdominance selection" in class I H-2 is one of the possible explanations for this.     

Salivary androgen-binding protein variation in Mus and other rodents.

We have searched for genetic variation in the expression of salivary androgen-binding protein (ABP) in a wide variety of mice and other rodents. ABP was present in the salivas of mice of all species and subspecies studied. Genetic studies have identified three common variants of the ABP Alpha subunit (Abpaa, Abpab, and Abpac) in Mus musculus populations with distributions that correspond roughly to those of the subspecies studied (domesticus, musculus, and castaneus, respectively). It appears that the ABP a and b polymorphisms conform to the hybrid zone between the domesticus and musculus subspecies characterized by others. Our studies suggest that the presence of Abpab in inbred strains may be due to a M. m. musculus contribution, perhaps via oriental fancy mice bred to European mice in the early lines leading to the common inbred strains. The relatively common occurrence of the ABP a type in other Mus species leads us to conclude that it is the ancestral type in mice. Further, the observation of what amounts to unique alleles in the three different subspecies indicates that microevolution of the protein has occurred. In a broader survey, ABP was also found in the salivas of Murid and Cricetid rodents generally. These findings suggest that ABP has an important functional role in rodent salivas.     

Reduced X-linked diversity in derived populations of house mice

Contrasting patterns of X-linked vs. autosomal diversity may be indicative of the mode of selection operating in natural populations. A number of observations have shown reduced X-linked (or Z-linked) diversity relative to autosomal diversity in various organisms, suggesting a large impact of genetic hitchhiking. However, the relative contribution of other forces such as population bottlenecks, variation in reproductive success of the two sexes, and differential introgression remains unclear. Here, we survey 13 loci, 6 X-linked and 7 autosomal, in natural populations of the house mouse (Mus musculus) subspecies complex. We studied seven populations of three different subspecies, the eastern house mouse M. musculus castaneus, the central house mouse M. m. musculus, and the western house mouse M. m. domesticus, including putatively ancestral and derived populations for each. All populations display lower diversity on the X chromosomes relative to autosomes, and this effect is most pronounced in derived populations. To assess the role of demography, we fit the demographic parameters that gave the highest likelihood of the data using coalescent simulations. We find that the reduction in X-linked diversity is too large to be explained by a simple demographic model in at least two of four derived populations. These observations are also not likely to be explained by differences in reproductive success between males and females. They are consistent with a greater impact of positive selection on the X chromosome, and this is supported by the observation of an elevated K(A) and elevated K(A)/K(S) ratios on the rodent X chromosome. A second contribution may be that the X chromosome less readily introgresses across subspecies boundaries.