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.     

Asymmetrical reproductive character displacement in the house mouse

Our study addressed reproductive character displacement between two subspecies of the house mouse, Mus musculus musculus and Mus musculus domesticus that hybridize in Europe along a zone where selection against hybridization is known to occur. Based on a multi-population approach, we investigated spatial patterns of divergence of mate preference in the two taxa. Mate preference was significantly higher in the contact zone than in allopatry in both subspecies, suggesting that reproductive character displacement occurs. Moreover, patterns of preference were stronger in M. m. musculus than in M. m. domesticus, indicating an asymmetrical divergence between the two. In the context of selection against hybridization, our results may provide empirical support for the hypothesis of reinforcement in a parapatric hybrid zone. We discuss factors that could explain the asymmetrical pattern of divergence and the possible impact of a unimodal structure on the maintenance of premating divergence between the two subspecies.     

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.     

Amplified ribosomal spacer sequence: structure and evolutionary origin

A novel class of repeated sequences consisting of tandem arrays of ribosomal spacer sequence has been discovered in a mouse genome. Comparison to normal ribosomal DNA reveals that one repeat unit consists of two separate parts of spacer sequence. This amplified spacer sequence has a pseudogene-like structure but is distinct from the previously reported pseudogenes and orphons in regions lacking coding sequences. So far the amplified spacer sequence has been found only in the BALB/c mouse genome but not in ten other laboratory strains and several wild-type mouse stocks. Surprisingly, a part of the amplified spacer sequence unit had a higher homology to the corresponding part of the ribosomal DNA sequence of Mus musculus molossinus, a Japanese wild-type mouse, than to the corresponding part of the rDNA of the BALB/c mouse. These findings suggest that the amplified spacer sequence of the BALB/c mouse might have partly originated in M. m. molossinus or in a related subspecies.     

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.     

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.     

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.     

Chromosomal dynamics of nucleolar organizer regions (NORs) in the house mouse: micro-evolutionary insights

Variation in the number and chromosomal location of nucleolar organizer regions (NORs) was studied in the house mouse, Mus musculus (2n=40). From an origin in Western Asia, this species colonized the Middle East, Europe and Asia. This expansion was accompanied by diversification into five subspecies. NOR diversity was revealed by fluorescence in situ hybridization using 18S and 28S probes on specimens spanning Asia to Western Europe. The results showed that the house mouse genome possessed a large number of NOR-bearing autosomes and a surprisingly high rate of polymorphism for the presence/absence of rRNA genes on all these chromosomes. All NOR sites were adjacent to the centromere except for two that were telomeric. Subspecific differentiation established from the NOR frequency data was concordant with the overall pattern of radiation proposed from molecular studies, but highlighted several discrepancies that need to be further addressed. NOR diversity in M. musculus consisted of a large number of polymorphic NORs that were common to at least two subspecies, and a smaller number of NORs that were unique to one subspecies. The most parsimonious scenario argues in favor of a subspecific differentiation by lineage sorting of ancestral NOR polymorphisms; only the unique NORs would have appeared by inter-chromosomal transposition, except for the two telomeric ones that may have originated by hybridization with another species. Such a scenario provides an alternative view from the one prevailing in most systematic and phylogenetic analyses that NORs have a high transposition rate due to concerted evolution of rRNA genes.     

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.     

B1 insertions as easy markers for mouse population studies

Few simple, easy-to-score PCR markers are available for studying genetic variation in wild mice populations belonging to Mus musculus at the population and subspecific levels. In this study, we show the abundant B1 family of short interspersed DNA elements (SINEs) is a very promising source of such markers. Thirteen B1 sequences from different regions of the genome were retrieved on the basis of their high degree of homology to a mouse consensus sequence, and the presence of these elements was screened for in wild derived mice representing M. spretus, macedonicus and spicilegus and the different subspecies of M. musculus. At five of these loci, varying degrees of insertion polymorphism were found in M. m. domesticus mice. These insertions were almost totally absent in the mice representing the other subspecies and species. Six other B1 elements were fixed in all the Mus species tested. At these loci, polymorphism associated with three restriction sites in the B1 consensus sequence was found in M. musculus. Most of these polymorphisms appear to be ancestral as they are shared by at least one of the other Mus species tested. Both insertion and restriction polymorphism revealed differences between five inbred laboratory strains considered to be of mainly domesticus origin, and at the six restriction loci a surprising number of these strains carried restriction variants that were either not found or very infrequent in domesticus. This suggests that in this particular group of loci, alleles of far Eastern origin are more frequent than expected.     

A TEST OF THE NEUTRAL MODEL OF EXPRESSION CHANGE IN NATURAL POPULATIONS OF HOUSE MOUSE SUBSPECIES

Changes in expression of genes are thought to contribute significantly to evolutionary divergence. To study the relative role of selection and neutrality in shaping expression changes, we analyzed 24 genes in three different tissues of the house mouse ( Mus musculus ). Samples from two natural populations of the subspecies M. m. domesticus and M. m. musculus were investigated using quantitative PCR assays and sequencing of the upstream region. We have developed an approach to quantify expression polymorphism within such populations and to disentangle technical from biological variation in the data. We found a correlation between expression polymorphism within populations and divergence between populations. Furthermore, we found a correlation between expression polymorphism and sequence polymorphism of the respective genes. These data are most easily interpreted within a framework of a predominantly neutral model of gene expression change, where only a fraction of the changes may have been driven by positive selection. Although most genes investigated were expressed in all three tissues analyzed, significant changes of expression levels occurred predominantly in a single tissue only. This adds to the notion that enhancer-specific effects or transregulatory effects can modulate the evolution of gene expression in a tissue-specific way.     

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.     

Characterization of the 23S and 5S rRNA genes and 23S-5S intergenic spacer region (ITS-2) of Photobacterium damselae

The 23S ribosomal RNA (rRNA) gene has been sequenced in strains of the fish pathogens Photobacterium damselae subsp. damselae (ATCC 33539) and subsp. piscicida (ATCC 29690), showing that 3 nucleotide positions are clearly different between subspecies. In addition, the 5S rRNA gene plus the intergenic spacer region between the 23S and 5S rRNA genes (ITS-2) were amplified, cloned and sequenced for the 2 reference strains as well as the field isolates RG91 (subsp. damselae) and DI21 (subsp. piscicida). A 100% similarity was found for the consensus 5S rRNA gene sequence in the 2 subspecies, although some microheterogeneity was detected as inter-cistronic variability within the same chromosome. Sequence analysis of the spacer region between the 23S and 5S rRNA genes revealed 2 conserved and 3 variable nucleotide sequence blocks, and 4 different modular organizations were found. The ITS-2 spacer region exhibited both inter-subspecies and intercistronic polymorphism, with a mosaic-like structure. The EMBL accession numbers for the 23S, 5S and ITS-2 sequences are: P. damselae subsp. piscicida 5S gene (AJ274379), P. damselae subsp. damselae 23S gene (Y18520), subsp. piscicida 23S gene (Y17901), P. damselae subsp. piscicida ITS-2 (AJ250695, AJ250696), P. damselae subsp. damselae ITS-2 (AJ250697, AJ250698).     

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.     

Evolutionary implication of heterogeneity of the nontranscribed spacer region of ribosomal DNA repeating units in various subspecies of Mus musculus

Genetic variability of the nontranscribed spacer (NTS) region within ribosomal DNA repeating units in the various subspecies of Mus musculus was determined. Mice belonging to several laboratory mouse strains were examined by means of Southern blot hybridization with a mouse ribosomal DNA probe. This probe encompasses the 3' end of the 28S ribosomal RNA (rRNA) gene and the following spacer. Restriction enzyme digestions of the liver DNAs from various wild mice revealed that each of the subspecies has a unique pattern in the spacer encompassing a distance approximately 10 kb downstream from the ribosomal gene. These restriction patterns permit the classification of mouse subspecies and also provide insights into the origin of the laboratory mouse strains.     

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.     

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.     

Organization of ribosomal RNA gene repeats of the mouse.

The organization of the ribosomal RNA (rRNA) genes of the mouse was determined by Southern blot hybridization using cloned rDNA fragments as probes, which could encompass the entire spacer region between two rRNA gene regions. The rRNA genes are organized into tandem repeats of nearly uniform length of about 44 kb. The heterogeneity detected in the nontranscribed spacer appears to be caused by its sequence rather than its length difference. At least three kinds of repetitive sequences are present in the non-transcribed spacer region; two of them are located 13 kb upstream from the 5'-end of 18S RNA gene and the other located 1 to 4 kb downstream from the 3'-end of 28S RNA gene.