Genetic dissection of a key reproductive barrier between nascent species of house mice

Reproductive isolation between species is often caused by deleterious interactions among loci in hybrids. Finding the genes involved in these incompatibilities provides insight into the mechanisms of speciation. With recently diverged subspecies, house mice provide a powerful system for understanding the genetics of reproductive isolation early in the speciation process. Although previous studies have yielded important clues about the genetics of hybrid male sterility in house mice, they have been restricted to F1 sterility or incompatibilities involving the X chromosome. To provide a more complete characterization of this key reproductive barrier, we conducted an F2 intercross between wild-derived inbred strains from two subspecies of house mice, Mus musculus musculus and Mus musculus domesticus. We identified a suite of autosomal and X-linked QTL that underlie measures of hybrid male sterility, including testis weight, sperm density, and sperm morphology. In many cases, the autosomal loci were unique to a specific sterility trait and exhibited an effect only when homozygous, underscoring the importance of examining reproductive barriers beyond the F1 generation. We also found novel two-locus incompatibilities between the M. m. musculus X chromosome and M. m. domesticus autosomal alleles. Our results reveal a complex genetic architecture for hybrid male sterility and suggest a prominent role for reproductive barriers in advanced generations in maintaining subspecies integrity in house mice.     

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.     

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.     

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.     

Asymmetry and polymorphism of hybrid male sterility during the early stages of speciation in house mice

House mice offer a powerful system for dissecting the genetic basis of phenotypes that isolate species in the early stages of speciation. We used a series of reciprocal crosses between wild-derived strains of Mus musculus and M. domesticus to examine F(1) hybrid male sterility, one of the primary phenotypes thought to isolate these species. We report four main results. First, we found significantly smaller testes and fewer sperm in hybrid male progeny of most crosses. Second, in some crosses hybrid male sterility was asymmetric and depended on the species origin of the X chromosome. These observations confirm and extend previous findings, underscoring the central role that the M. musculus X chromosome plays in reproductive isolation. Third, comparisons among reciprocal crosses revealed polymorphism at one or more hybrid incompatibilities within M. musculus. Fourth, the spermatogenic phenotype of this polymorphic interaction appears distinct from previously described hybrid incompatibilities between these species. These data build on previous studies of speciation in house mice and show that the genetic basis of hybrid male sterility is fairly complex, even at this early stage of divergence.     

A complex genetic basis to X-linked hybrid male sterility between two species of house mice

The X chromosome plays a central role in the evolution of reproductive isolation, but few studies have examined the genetic basis of X-linked incompatibilities during the early stages of speciation. We report the results of a large experiment focused on the reciprocal introgression of the X chromosome between two species of house mice, Mus musculus and M. domesticus. Introgression of the M. musculus X chromosome into a wild-derived M. domesticus genetic background produced male-limited sterility, qualitatively consistent with previous experiments using classic inbred strains to represent M. domesticus. The genetic basis of sterility involved a minimum of four X-linked factors. The phenotypic effects of major sterility QTL were largely additive and resulted in complete sterility when combined. No sterility factors were uncovered on the M. domesticus X chromosome. Overall, these results revealed a complex and asymmetric genetic basis to X-linked hybrid male sterility during the early stages of speciation in mice. Combined with data from previous studies, we identify one relatively narrow interval on the M. musculus X chromosome involved in hybrid male sterility. Only a handful of spermatogenic genes are within this region, including one of the most rapidly evolving genes on the mouse X chromosome.     

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.     

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.     

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.     

Reduced male fertility is common but highly variable in form and severity in a natural house mouse hybrid zone

Barriers to gene flow between naturally hybridizing taxa reveal the initial stages of speciation. Reduced hybrid fertility is a common feature of reproductive barriers separating recently diverged species. In house mice (Mus musculus), hybrid male sterility has been studied extensively using experimental crosses between subspecies. Here, we present the first detailed picture of hybrid male fertility in the European M. m. domesticus-M. m. musculus hybrid zone. Complete sterility appears rare or absent in natural hybrids but a large proportion of males (~30%) have sperm count or relative testis weight below the range in pure subspecies, and likely suffer reduced fertility. Comparison of a suite of traits related to fertility among subfertile males indicates reduced hybrid fertility in the contact zone is highly variable among individuals and ancestry groups in the type, number, and severity of spermatogenesis defects present. Taken together, these results suggest multiple underlying genetic incompatibilities are segregating in the hybrid zone, which likely contribute to reproductive isolation between subspecies.     

Do changes in gene expression contribute to sexual isolation and reinforcement in the house mouse?

Expression divergence, rather than sequence divergence, has been shown to be important in speciation, particularly in the early stages of divergence of traits involved in reproductive isolation. In the two European subspecies of house mice, Mus musculus musculus and Mus musculus domesticus, earlier studies have demonstrated olfactory‐based assortative mate preference in populations close to their hybrid zone. It has been suggested that this behaviour evolved following the recent secondary contact between the two taxa (~3,000 years ago) in response to selection against hybridization. To test for a role of changes in gene expression in the observed behavioural shift, we conducted a RNA sequencing experiment on mouse vomeronasal organs. Key candidate genes for pheromone‐based subspecies recognition, the vomeronasal receptors, are expressed in these organs. Overall patterns of gene expression varied significantly between samples from the two subspecies, with a large number of differentially expressed genes between the two taxa. In contrast, only ~200 genes were found repeatedly differentially expressed between populations within M. m. musculus that did or did not display assortative mate preferences (close to or more distant from the hybrid zone, respectively), with an overrepresentation of genes belonging to vomeronasal receptor family 2. These receptors are known to play a key role in recognition of chemical cues that handle information about genetic identity. Interestingly, four of five of these differentially expressed receptors belong to the same phylogenetic cluster, suggesting specialization of a group of closely related receptors in the recognition of odorant signals that may allow subspecies recognition and assortative mating.     

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.     

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.     

Experimental trypanosomiasis of natural hybrids between house mouse subspecies.

This study characterises the extent of the susceptibility to parasites (first demonstrated with helminths) of hybrids between Mus musculus domesticus and Mus musculus musculus. Experimental infections with Trypanosoma musculi of M. m. domesticus, M. m. musculus and their natural hybrids have been performed to compare their level of resistance/susceptibility. It appears that contrary to the results with helminths, hybrid mice present the same level of resistance/susceptibility to the trypanosome as M. m. musculus and M. m. domesticus individuals. This result is interpreted in the light of the modalities of host parasite interactions and leads us to hypothesise on the role of parasitism in the evolution of the house mouse hybrid zone.     

The complex history of a gene proposed to participate in a sexual isolation mechanism in house mice

Previous behavioral experiments showed that mouse salivary androgen-binding protein (ABP) was involved in interindividual recognition and might play a role in sexual isolation between house mouse (Mus musculus) subspecies. The pattern of evolution of Abpa, the gene for the alpha subunit of ABP, was found to be consistent with this hypothesis. Abpa apparently diverged rapidly between species and subspecies with a large excess of nonsynonymous substitutions, a lack of exon polymorphism within each of the three subspecies, and a lack of intron polymorphism in the one subspecies studied (M. musculus domesticus). Here we characterized the intron and exon sequence variations of this gene in house mouse populations from central Eurasia, a region yet unsampled and thought to be close to the cradle of the radiation of the subspecies. We also determined the intron and exon sequences in seven other species of the genus Mus. We confirmed the general pattern of rapid evolution by essentially nonsynonymous substitutions, both inter- and intraspecifically, supporting the idea that Darwinian selection has driven the evolution of this gene. We also observed a uniform intron sequence in five samples of M. musculus musculus, suggesting that a selective sweep might have occurred for that allele. In contrast to previous results, however, we found extensive intron and exon polymorphism in some house mouse populations from central Eurasia. We also found evidence for secondary admixture of the subspecies-specific alleles in regions of transition between the subspecies in central Eurasia. Furthermore, an abnormal intron phylogeny suggested that interspecific exchanges had occurred between the house mouse subspecies and three other Palearctic species. These observations appear to be at variance with the simple hypothesis that Abpa is involved in reproductive isolation. Although we do not rule out a role in recognition, the situation appears to be more complex than previously thought. Thus the selective mechanism behind the evolution of Abpa remains to be resolved, and we suggest that it may have changed during the recent colonization history of the house mouse.     

The contribution of the y chromosome to hybrid male sterility in house mice

Hybrid sterility in the heterogametic sex is a common feature of speciation in animals. In house mice, the contribution of the Mus musculus musculus X chromosome to hybrid male sterility is large. It is not known, however, whether F(1) male sterility is caused by X-Y or X-autosome incompatibilities or a combination of both. We investigated the contribution of the M. musculus domesticus Y chromosome to hybrid male sterility in a cross between wild-derived strains in which males with a M. m. musculus X chromosome and M. m. domesticus Y chromosome are partially sterile, while males from the reciprocal cross are reproductively normal. We used eight X introgression lines to combine different X chromosome genotypes with different Y chromosomes on an F(1) autosomal background, and we measured a suite of male reproductive traits. Reproductive deficits were observed in most F(1) males, regardless of Y chromosome genotype. Nonetheless, we found evidence for a negative interaction between the M. m. domesticus Y and an interval on the M. m. musculus X that resulted in abnormal sperm morphology. Therefore, although F(1) male sterility appears to be caused mainly by X-autosome incompatibilities, X-Y incompatibilities contribute to some aspects of sterility.     

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.     

Measures of linkage disequilibrium among neighbouring SNPs indicate asymmetries across the house mouse hybrid zone

Theory predicts that naturally occurring hybrid zones between genetically distinct taxa can move over space and time as a result of selection and/or demographic processes, with certain types of hybrid zones being more or less likely to move. Determining whether a hybrid zone is stationary or moving has important implications for understanding evolutionary processes affecting interactions in hybrid populations. However, direct observations of hybrid zone movement are difficult to make unless the zone is moving rapidly. Here, evidence for movement in the house mouse Mus musculus domesticus × Mus musculus musculus hybrid zone is provided using measures of LD and haplotype structure among neighbouring SNP markers from across the genome. Local populations of mice across two transects in Germany and the Czech Republic were sampled, and a total of 1301 mice were genotyped at 1401 markers from the nuclear genome. Empirical measures of LD provide evidence for extinction and (re)colonization in single populations and, together with simulations, suggest hybrid zone movement because of either geography-dependent asymmetrical dispersal or selection favouring one subspecies over the other.     

There is no heterotic effect upon developmental stability in the ventral side of the skull within the house mouse hybrid zone

Fluctuating asymmetry (FA) of a body is expected to be related to an organism's developmental instability. We studied patterns of FA in the ventral side of the skull along a transect across the central-European portion of the hybrid zone between two house mouse subspecies, Mus musculus musculus and Mus musculus domesticus. We found that FA was not significantly different between introgression classes relative to differences between localities within those classes and that the within-class differences were much larger than differences between individual-signed asymmetries within localities. However, if year was added to the same analysis as another factor, FA was not significantly different among localities within the same introgression class. When individual asymmetries were plotted against individual hybrid indices, hybrids appeared more asymmetric than individuals from outside of the zone. Thus contrary to previous studies, we did not find lower FA indicating heterotic effect in hybrids for the traits studied. It is suggested that the impact of hybridization on FA in the ventral side of the mouse skull is negligible or overwhelmed by other factors.