Does geography matter in hybrid sterility in house mice?

The contribution of sex chromosome conflict to disrupted spermatogenesis in hybrid house mice

Incompatibilities on the sex chromosomes are important in the evolution of hybrid male sterility, but the evolutionary forces underlying this phenomenon are unclear. House mice (Mus musculus) lineages have provided powerful models for understanding the genetic basis of hybrid male sterility. X chromosome-autosome interactions cause strong incompatibilities in M. musculus F1 hybrids, but variation in sterility phenotypes suggests a more complex genetic basis. In addition, XY chromosome conflict has resulted in rapid expansions of ampliconic genes with dosage-dependent expression that is essential to spermatogenesis. Here, we evaluated the contribution of XY lineage mismatch to male fertility and stage-specific gene expression in hybrid mice. We performed backcrosses between two house mouse subspecies to generate reciprocal Y-introgression strains and used these strains to test the effects of XY mismatch in hybrids. Our transcriptome analyses of sorted spermatid cells revealed widespread overexpression of the X chromosome in sterile F1 hybrids independent of Y chromosome subspecies origin. Thus, postmeiotic overexpression of the X chromosome in sterile F1 mouse hybrids is likely a downstream consequence of disrupted meiotic X-inactivation rather than XY gene copy number imbalance. Y chromosome introgression did result in subfertility phenotypes and disrupted expression of several autosomal genes in mice with an otherwise nonhybrid genomic background, suggesting that Y-linked incompatibilities contribute to reproductive barriers, but likely not as a direct consequence of XY conflict. Collectively, these findings suggest that rapid sex chromosome gene family evolution driven by genomic conflict has not resulted in strong male reproductive barriers between these subspecies of house mice.

Stage-specific disruption of X chromosome expression during spermatogenesis in sterile house mouse hybrids

Hybrid sterility is a complex phenotype that can result from the breakdown of spermatogenesis at multiple developmental stages. Here, we disentangle two proposed hybrid male sterility mechanisms in the house mice, Mus musculus domesticus and M. m. musculus, by comparing patterns of gene expression in sterile F1 hybrids from a reciprocal cross. We found that hybrid males from both cross directions showed disrupted X chromosome expression during prophase of meiosis I consistent with a loss of meiotic sex chromosome inactivation (MSCI) and Prdm9-associated sterility, but that the degree of disruption was greater in mice with an M. m. musculus X chromosome consistent with previous studies. During postmeiotic development, gene expression on the X chromosome was only disrupted in one cross direction, suggesting that misexpression at this later stage was genotype-specific and not a simple downstream consequence of MSCI disruption which was observed in both reciprocal crosses. Instead, disrupted postmeiotic expression may depend on the magnitude of earlier disrupted MSCI, or the disruption of particular X-linked genes or gene networks. Alternatively, only hybrids with a potential deficit of Sly copies, a Y-linked ampliconic gene family, showed overexpression in postmeiotic cells, consistent with a previously proposed model of antagonistic coevolution between the X- and Y-linked ampliconic genes contributing to disrupted expression late in spermatogenesis. The relative contributions of these two regulatory mechanisms and their impact on sterility phenotypes await further study. Our results further support the hypothesis that X-linked hybrid sterility in house mice has a variable genetic basis, and that genotype-specific disruption of gene regulation contributes to overexpression of the X chromosome at different stages of development. Overall, these findings underscore the critical role of epigenetic regulation of the X chromosome during spermatogenesis and suggest that these processes are prone to disruption in hybrids.

Hybrid sterility genes in mice (Mus musculus): a peculiar case of PRDM9 incompatibility

Hybrid sterility is a critical step in the evolution of reproductive barriers between diverging taxa during the process of speciation. Recent studies of young subspecies of the house mouse revealed a multigenic nature and frequent polymorphism of hybrid sterility genes as well as the recurrent engagement of the meiosis-specific gene PR domain-containing 9 (Prdm9) and X-linked loci. Prdm9-controlled hybrid sterility is essentially chromosomal in nature, conditioned by the sequence divergence between subspecies. Depending on the Prdm9 interallelic interactions and the X-linked Hstx2 locus, the same homologs either regularly recombine and synapse, or show impaired DNA DSB repair, asynapsis, and early meiotic arrest. Thus, Prdm9-dependent hybrid sterility points to incompatibilities affecting meiotic recombination as a possible mechanism of reproductive isolation between (sub)species.

Disrupted Gene Networks in Subfertile Hybrid House Mice

The Dobzhansky–Muller (DM) model provides a widely accepted mechanism for the evolution of reproductive isolation: incompatible substitutions disrupt interactions between genes. To date, few candidate incompatibility genes have been identified, leaving the genes driving speciation mostly uncharacterized. The importance of interactions in the DM model suggests that gene coexpression networks provide a powerful framework to understand disrupted pathways associated with postzygotic isolation. Here, we perform weighted gene coexpression network analysis to infer gene interactions in hybrids of two recently diverged European house mouse subspecies, Mus mus domesticus and M. m. musculus, which commonly show hybrid male sterility or subfertility. We use genome-wide testis expression data from 467 hybrid mice from two mapping populations: F₂s from a laboratory cross between wild-derived pure subspecies strains and offspring of natural hybrids captured in the Central Europe hybrid zone. This large data set enabled us to build a robust consensus network using hybrid males with fertile phenotypes. We identify several expression modules, or groups of coexpressed genes, that are disrupted in subfertile hybrids, including modules functionally enriched for spermatogenesis, cilium and sperm flagellum organization, chromosome organization, and DNA repair, and including genes expressed in spermatogonia, spermatocytes, and spermatids. Our network-based approach enabled us to hone in on specific hub genes likely to be influencing module-wide gene expression and hence potentially driving large-effect DM incompatibilities. A disproportionate number of hub genes lie within sterility loci identified previously in the hybrid zone mapping population and represent promising candidate barrier genes and targets for future functional analysis.

Prdm9 Intersubspecific Interactions in Hybrid Male Sterility of House Mouse

The classical definition posits hybrid sterility as a phenomenon when two parental taxa each of which is fertile produce a hybrid that is sterile. The first hybrid sterility gene in vertebrates, Prdm9, coding for a histone methyltransferase, was identified in crosses between two laboratory mouse strains derived from Mus mus musculus and M. m. domesticus subspecies. The unique function of PRDM9 protein in the initiation of meiotic recombination led to the discovery of the basic molecular mechanism of hybrid sterility in laboratory crosses. However, the role of this protein as a component of reproductive barrier outside the laboratory model remained unclear. Here, we show that the Prdm9 allelic incompatibilities represent the primary cause of reduced fertility in intersubspecific hybrids between M. m. musculus and M. m. domesticus including 16 musculus and domesticus wild-derived strains. Disruption of fertility phenotypes correlated with the rate of failure of synapsis between homologous chromosomes in meiosis I and with early meiotic arrest. All phenotypes were restored to normal when the domesticus Prdm9dom2 allele was substituted with the Prdm9dom2H humanized variant. To conclude, our data show for the first time the male infertility of wild-derived musculus and domesticus subspecies F1 hybrids controlled by Prdm9 as the major hybrid sterility gene. The impairment of fertility surrogates, testes weight and sperm count, correlated with increasing difficulties of meiotic synapsis of homologous chromosomes and with meiotic arrest, which we suppose reflect the increasing asymmetry of PRDM9-dependent DNA double-strand breaks.

Sperm quality, aggressiveness and generation turnover may facilitate unidirectional Y chromosome introgression across the European house mouse hybrid zone

The widespread and locally massive introgression of Y chromosomes of the eastern house mouse (Mus musculus musculus) into the range of the western subspecies (M. m. domesticus) in Central Europe calls for an explanation of its underlying mechanisms. Given the paternal inheritance pattern, obvious candidates for traits mediating the introgression are characters associated with sperm quantity and quality. We can also expect traits such as size, aggression or the length of generation cycles to facilitate the spread. We have created two consomic strains carrying the non-recombining region of the Y chromosome of the opposite subspecies, allowing us to study introgression in both directions, something impossible in nature due to the unidirectionality of introgression. We analyzed several traits potentially related to male fitness. Transmission of the domesticus Y onto the musculus background had negative effects on all studied traits. Likewise, domesticus males possessing the musculus Y had, on average, smaller body and testes and lower sperm count than the parental strain. However, the same consomic males tended to produce less- dissociated sperm heads, to win more dyadic encounters, and to have shorter generation cycles than pure domesticus males. These data suggest that the domesticus Y is disadvantageous on the musculus background, while introgression in the opposite direction can confer a recognizable, though not always significant, selective advantage. Our results are thus congruent with the unidirectional musculus → domesticus Y chromosome introgression in Central Europe. In addition to some previous studies, they show this to be a multifaceted phenomenon demanding a multidisciplinary approach.

Phenotypic effects of the Y chromosome are variable and structured in hybrids among house mouse recombinant lines

Hybrid zones between divergent populations sieve genomes into blocks that introgress across the zone, and blocks that do not, depending on selection between interacting genes. Consistent with Haldane's rule, the Y chromosome has been considered counterselected and hence not to introgress across the European house mouse hybrid zone. However, recent studies detected massive invasion of M. m. musculus Y chromosomes into M. m. domesticus territory. To understand mechanisms facilitating Y spread, we created 31 recombinant lines from eight wild-derived strains representing four localities within the two mouse subspecies. These lines were reciprocally crossed and resulting F1 hybrid males scored for five phenotypic traits associated with male fitness. Molecular analyses of 51 Y-linked SNPs attributed ~50% of genetic variation to differences between the subspecies and 8% to differentiation within both taxa. A striking proportion, 21% (frequencies of sperm head abnormalities) and 42% (frequencies of sperm tail dissociations), of phenotypic variation was explained by geographic Y chromosome variants. Our crossing design allowed this explanatory power to be examined across a hierarchical scale from subspecific to local intrastrain effects. We found that divergence and variation were expressed diversely in different phenotypic traits and varied across the whole hierarchical scale. This finding adds another dimension of complexity to studies of Y introgression not only across the house mouse hybrid zone but potentially also in other contact zones.

Genetic vs environment influences on house mouse hybrid zone in Iran

Genetic divergence and environment influence on speciation process are the great deal studies over recent decades. One of the best ways for exploring the interaction of geography and genetics is the evaluation of hybrids in a contact zone. To understand if there is one or more hybrid zone between house mouse subspecies in Iran and what are the differences comparing these zones with European well-known hybrid zone, we performed this approach. Samples were live-trapped from Ilam city in west for sensu lato M. m. domesticus subspecies, and Neishabur city in north-east of Iran for sensu lato M. m. musculus subspecies. In five experimental groups, male and female mice of the two subspecies were crossed reciprocally to generate F1 hybrids, and then F1 offspring males and females were crossed also reciprocally between siblings to generate F2 hybrids. In the same manner as seen in European hybrid zone, hybridization between female M. m. musculus and male M. m. domesticus of all five groups showed male sterility in F1 generation, but intact female offspring. These sterile males comparing with a parent or healthy males showed low count and more abnormal sperm percentage in morphological and testis histological section studies. Comparing the results from this study with numerous studies carried out during several years on the European hybrid zone showed an equal condition of contact between two subspecies. Genetical elements have kept their same influence on postzygotic reproductive isolation more than environmental effects far from the Europe, here in Iran.

Genetic Dissection of Hybrid Male Sterility Across Stages of Spermatogenesis

Hybrid sterility is a common form of reproductive isolation between nascent species. Although hybrid sterility is routinely documented and genetically dissected in speciation studies, its developmental basis is rarely examined, especially in generations beyond the F1 generation. To identify phenotypic and genetic determinants of hybrid male sterility from a developmental perspective, we characterized testis histology in 312 F2 hybrids generated by intercrossing inbred strains of Mus musculus domesticus and M. m. musculus, two subspecies of house mice. Hybrids display a range of histologic abnormalities that indicate defective spermatogenesis. Among these abnormalities, we quantified decreased testis size, reductions in spermatocyte and spermatid number, increased apoptosis of meiosis I spermatocytes, and more multinucleated syncytia. Collectively, our phenotypic data point to defects in meiosis I as a primary barrier to reproduction. We identified seven quantitative trait loci (QTL) controlling five histologic traits. A region of chromosome 17 that contains Prdm9, a gene known to confer F1 hybrid male sterility, affects multinucleated syncytia and round spermatids, potentially extending the phenotypic outcomes of this incompatibility. The X chromosome also plays a key role, with loci affecting multinucleated syncytia, apoptosis of round spermatids, and round spermatid numbers. We detected an epistatic interaction between QTL on chromosomes 17 and X for multinucleated syncytia. Our results refine the developmental basis of a key reproductive barrier in a classic model system for speciation genetics.

The Evolution of Polymorphic Hybrid Incompatibilities in House Mice

Resolving the mechanistic and genetic bases of reproductive barriers between species is essential to understanding the evolutionary forces that shape speciation. Intrinsic hybrid incompatibilities are often treated as fixed between species, yet there can be considerable variation in the strength of reproductive isolation between populations. The extent and causes of this variation remain poorly understood in most systems. We investigated the genetic basis of variable hybrid male sterility (HMS) between two recently diverged subspecies of house mice, Mus musculus domesticus and Mus musculus musculus We found that polymorphic HMS has a surprisingly complex genetic basis, with contributions from at least five autosomal loci segregating between two closely related wild-derived strains of M. m. musculus One of the HMS-linked regions on chromosome 4 also showed extensive introgression among inbred laboratory strains and transmission ratio distortion (TRD) in hybrid crosses. Using additional crosses and whole genome sequencing of sperm pools, we showed that TRD was limited to hybrid crosses and was not due to differences in sperm motility between M. m. musculus strains. Based on these results, we argue that TRD likely reflects additional incompatibilities that reduce hybrid embryonic viability. In some common inbred strains of mice, selection against deleterious interactions appears to have unexpectedly driven introgression at loci involved in epistatic hybrid incompatibilities. The highly variable genetic basis to F1 hybrid incompatibilities between closely related mouse lineages argues that a thorough dissection of reproductive isolation will require much more extensive sampling of natural variation than has been commonly utilized in mice and other model systems.

Inconsistent reproductive isolation revealed by interactions between Catostomus fish species

Interactions between species are central to evolution and ecology, but we do not know enough about how outcomes of interactions between species vary across geographic locations, in heterogeneous environments, or over time. Ecological dimensions of interactions between species are known to vary, but evolutionary interactions such as the establishment and maintenance of reproductive isolation are often assumed to be consistent across instances of an interaction between species. Hybridization among Catostomus fish species occurs over a large and heterogeneous geographic area and across taxa with distinct evolutionary histories, which allows us to assess consistency in species interactions. We analyzed hybridization among six Catostomus species across the Upper Colorado River basin (US mountain west) and found extreme variation in hybridization across locations. Different hybrid crosses were present in different locations, despite similar species assemblages. Within hybrid crosses, hybridization varied from only first generation hybrids to extensive hybridization with backcrossing. Variation in hybridization outcomes might result from uneven fitness of hybrids across locations, polymorphism in genetic incompatibilities, chance, unidentified historical contingencies, or some combination thereof. Our results suggest caution in assuming that one or a few instances of hybridization represent all interactions between the focal species, as species interactions vary substantially across locations.

Molecular heterogeneity in major urinary proteins of Mus musculus subspecies: potential candidates involved in speciation

When hybridisation carries a cost, natural selection is predicted to favour evolution of traits that allow assortative mating (reinforcement). Incipient speciation between the two European house mouse subspecies, Mus musculus domesticus and M.m.musculus, sharing a hybrid zone, provides an opportunity to understand evolution of assortative mating at a molecular level. Mouse urine odours allow subspecific mate discrimination, with assortative preferences evident in the hybrid zone but not in allopatry. Here we assess the potential of MUPs (major urinary proteins) as candidates for signal divergence by comparing MUP expression in urine samples from the Danish hybrid zone border (contact) and from allopatric populations. Mass spectrometric characterisation identified novel MUPs in both subspecies involving mostly new combinations of amino acid changes previously observed in M.m.domesticus. The subspecies expressed distinct MUP signatures, with most MUPs expressed by only one subspecies. Expression of at least eight MUPs showed significant subspecies divergence both in allopatry and contact zone. Another seven MUPs showed divergence in expression between the subspecies only in the contact zone, consistent with divergence by reinforcement. These proteins are candidates for the semiochemical barrier to hybridisation, providing an opportunity to characterise the nature and evolution of a putative species recognition signal.

The Pace of Hybrid Incompatibility Evolution in House Mice

Hybrids between species are often sterile or inviable. This form of reproductive isolation is thought to evolve via the accumulation of mutations that interact to reduce fitness when combined in hybrids. Mathematical formulations of this "Dobzhansky-Muller model" predict an accelerating buildup of hybrid incompatibilities with divergence time (the "snowball effect"). Although the Dobzhansky-Muller model is widely accepted, the snowball effect has only been tested in two species groups. We evaluated evidence for the snowball effect in the evolution of hybrid male sterility among subspecies of house mice, a recently diverged group that shows partial reproductive isolation. We compared the history of subspecies divergence with patterns of quantitative trait loci (QTL) detected in F2 intercrosses between two pairs of subspecies (Mus musculus domesticus with M. m. musculus and M. m. domesticus with M. m. castaneus). We used a recently developed phylogenetic comparative method to statistically measure the fit of these data to the snowball prediction. To apply this method, QTL were partitioned as either shared or unshared in the two crosses. A heuristic partitioning based on the overlap of QTL confidence intervals produced unambiguous support for the snowball effect. An alternative approach combining data among crosses favored the snowball effect for the autosomes, but a linear accumulation of incompatibilities for the X chromosome. Reasoning that the X chromosome analyses are complicated by low mapping resolution, we conclude that hybrid male sterility loci have snowballed in house mice. Our study illustrates the power of comparative genetic mapping for understanding mechanisms of speciation.

Highly variable reproductive isolation among pairs of Catostomus species

Hybridization between diverged taxa tests the strength of reproductive isolation and can therefore reveal mechanisms of reproductive isolation. However, it remains unclear how consistent reproductive isolation is across species' ranges and to what extent reproductive isolation might remain polymorphic as species diverge. To address these questions, we compared outcomes of hybridization across species pairs of Catostomus fishes in three rivers in the Upper Colorado River basin, where an introduced species, C. commersoni, hybridizes with at least two native species, C. discobolus and C. latipinnis. We observed substantial heterogeneity in outcomes of hybridization, both between species pairs and across geographically separate rivers within each species pair. We also observed hybridization of additional related species with our focal species, suggesting that reproductive isolation in this group involves interactions of multiple evolutionary and ecological factors. These findings suggest that a better understanding of the determinants of variation in reproductive isolation is needed and that studies of reproductive isolation in hybrids should consider how the dynamics and mechanisms of reproductive isolation vary over ecological space and over evolutionary time. Our results also have implications for the conservation and management of native catostomids in the Colorado River basin. Heterogeneity in outcomes of hybridization suggests that the threat posed by hybridization and genetic introgression to the persistence of native species probably varies with extent of reproductive isolation, both across rivers and across species pairs.

Maternal-foetal genomic conflict and speciation: no evidence for hybrid placental dysplasia in crosses between two house mouse subspecies

Interspecific hybridization between closely related mammalian species, including various species of the genus Mus, is commonly associated with abnormal growth of the placenta and hybrid foetuses, a phenomenon known as hybrid placental dysplasia (HPD). The role of HPD in speciation is anticipated but still poorly understood. Here, we studied placental and foetal growth in F1 crosses between four inbred mouse strains derived from two house mouse subspecies, Mus musculus musculus and Mus musculus domesticus. These subspecies are in the early stage of speciation and still hybridize in nature. In accordance with the maternal-foetal genomic conflict hypothesis, we found different parental influences on placental and foetal development, with placental weight most affected by the father's body weight and foetal weight by the mother's body weight. After removing the effects of parents' body weight, we did not find any significant differences in foetal or placental weights between intra-subspecific and inter-subspecific F1 crosses. Nevertheless, we found that the variability in placental weight in inter-subspecific crosses is linked to the X chromosome, similarly as for HPD in interspecific mouse crosses. Our results suggest that maternal-foetal genomic conflict occurs in the house mouse system, but has not yet diverged sufficiently to cause abnormalities in placental and foetal growth in inter-subspecific crosses. HPD is thus unlikely to contribute to speciation in the house mouse system. However, we cannot rule out that it might have contributed to other speciation events in the genus Mus, where differences in the levels of polyandry exist between the species.

Functional organization of the genome may shape the species boundary in the house mouse

Genomic features such as rate of recombination and differentiation have been suggested to play a role in species divergence. However, the relationship of these phenomena to functional organization of the genome in the context of reproductive isolation remains unexplored. Here, we examine genomic characteristics of the species boundaries between two house mouse subspecies (Mus musculus musculus/M. m. domesticus). These taxa form a narrow semipermeable zone of secondary contact across Central Europe. Due to the incomplete nature of reproductive isolation, gene flow in the zone varies across the genome. We present an analysis of genomic differentiation, rate of recombination, and functional composition of genes relative to varying amounts of introgression. We assessed introgression using 1,316 autosomal single nucleotide polymorphism markers, previously genotyped in hybrid populations from three transects. We found a significant relationship between amounts of introgression and both genomic differentiation and rate of recombination with genomic regions of reduced introgression associated with higher genomic differentiation and lower rates of recombination, and the opposite for genomic regions of extensive introgression. We also found a striking functional polarization of genes based on where they are expressed in the cell. Regions of elevated introgression exhibit a disproportionate number of genes involved in signal transduction functioning at the cell periphery, among which olfactory receptor genes were found to be the most prominent group. Conversely, genes expressed intracellularly and involved in DNA binding were the most prevalent in regions of reduced introgression. We hypothesize that functional organization of the genome is an important driver of species divergence.

Prdm9 incompatibility controls oligospermia and delayed fertility but no selfish transmission in mouse intersubspecific hybrids

PR-domain 9 (Prdm9) is the first hybrid sterility gene identified in mammals. The incompatibility between Prdm9 from Mus musculus domesticus (Mmd; the B6 strain) and the Hstx2 region of chromosome (Chr) X from M. m. musculus (Mmm; the PWD strain) participates in the complete meiotic arrest of mouse intersubspecific (PWD×B6)F1 hybrid males. Other studies suggest that also semisterile intersubspecific hybrids are relevant for mouse speciation, but the genes responsible remain unknown. To investigate the causes of this semisterility, we analyzed the role of Prdm9 and Chr X in hybrids resulting from the crosses of PWK, another Mmm-derived inbred strain. We demonstrate that Prdm9 and Chr X control the partial meiotic arrest and reduced sperm count in (PWK×B6)F1 males. Asynapsis of heterosubspecific chromosomes and semisterility were partially suppressed by removal of the B6 allele of Prdm9. Polymorphisms between PWK and PWD on Chr X but not in the Prdm9 region were responsible for the modification of the outcome of Prdm9-Chr X F1 hybrid incompatibility. Furthermore, (PWK×B6)F1 hybrid males displayed delayed fertility dependent on the Prdm9 incompatibility. While the Drosophila hybrid sterility gene Overdrive causes both delayed fertility and increased transmission of its own chromosome to the offspring, the segregation of Chr X and the Prdm9 region from the mouse (PWK×B6)F1 males was normal. Our results indicate extended functional consequences of Prdm9-Chr X intersubspecific incompatibility on the fertility of hybrids and should influence the design of fertility analyses in hybrid zones and of laboratory crosses between Mmm and Mmd strains.

Genomic networks of hybrid sterility

Hybrid dysfunction, a common feature of reproductive barriers between species, is often caused by negative epistasis between loci ("Dobzhansky-Muller incompatibilities"). The nature and complexity of hybrid incompatibilities remain poorly understood because identifying interacting loci that affect complex phenotypes is difficult. With subspecies in the early stages of speciation, an array of genetic tools, and detailed knowledge of reproductive biology, house mice (Mus musculus) provide a model system for dissecting hybrid incompatibilities. Male hybrids between M. musculus subspecies often show reduced fertility. Previous studies identified loci and several X chromosome-autosome interactions that contribute to sterility. To characterize the genetic basis of hybrid sterility in detail, we used a systems genetics approach, integrating mapping of gene expression traits with sterility phenotypes and QTL. We measured genome-wide testis expression in 305 male F2s from a cross between wild-derived inbred strains of M. musculus musculus and M. m. domesticus. We identified several thousand cis- and trans-acting QTL contributing to expression variation (eQTL). Many trans eQTL cluster into eleven 'hotspots,' seven of which co-localize with QTL for sterility phenotypes identified in the cross. The number and clustering of trans eQTL-but not cis eQTL-were substantially lower when mapping was restricted to a 'fertile' subset of mice, providing evidence that trans eQTL hotspots are related to sterility. Functional annotation of transcripts with eQTL provides insights into the biological processes disrupted by sterility loci and guides prioritization of candidate genes. Using a conditional mapping approach, we identified eQTL dependent on interactions between loci, revealing a complex system of epistasis. Our results illuminate established patterns, including the role of the X chromosome in hybrid sterility. The integrated mapping approach we employed is applicable in a broad range of organisms and we advocate for widespread adoption of a network-centered approach in speciation genetics.

Sperm-related phenotypes implicated in both maintenance and breakdown of a natural species barrier in the house mouse

The house mouse hybrid zone (HMHZ) is a species barrier thought to be maintained by a balance between dispersal and natural selection against hybrids. While the HMHZ is characterized by frequency discontinuities for some sex chromosome markers, there is an unexpected large-scale regional introgression of a Y chromosome across the barrier, in defiance of Haldane's rule. Recent work suggests that a major force maintaining the species barrier acts through sperm traits. Here, we test whether the Y chromosome penetration of the species barrier acts through sperm traits by assessing sperm characteristics of wild-caught males directly in a field laboratory set up in a Y introgression region of the HMHZ, later calculating the hybrid index of each male using 1401 diagnostic single nucleotide polymorphisms (SNPs). We found that both sperm count (SC) and sperm velocity were significantly reduced across the natural spectrum of hybrids. However, SC was more than rescued in the presence of the invading Y. Our results imply an asymmetric advantage for Y chromosome introgression consistent with the observed large-scale introgression. We suggest that selection on sperm-related traits probably explains a large component of patterns observed in the natural hybrid zone, including the Y chromosome penetration.

Interallelic and intergenic incompatibilities of the Prdm9 (Hst1) gene in mouse hybrid sterility

The Dobzhansky-Muller model of incompatibilities explains reproductive isolation between species by incorrect epistatic interactions. Although the mechanisms of speciation are of great interest, no incompatibility has been characterized at the gene level in mammals. The Hybrid sterility 1 gene (Hst1) participates in the arrest of meiosis in F₁ males of certain strains from two Mus musculus subspecies, e.g., PWD from M. m. musculus and C57BL/6J (henceforth B6) from M. m. domesticus. Hst1 has been identified as a meiotic PR-domain gene (Prdm9) encoding histone 3 methyltransferase in the male offspring of PWD females and B6 males, (PWD×B6)F₁. To characterize the incompatibilities underlying hybrid sterility, we phenotyped reproductive and meiotic markers in males with altered copy numbers of Prdm9. A partial rescue of fertility was observed upon removal of the B6 allele of Prdm9 from the azoospermic (PWD×B6)F₁ hybrids, whereas removing one of the two Prdm9 copies in PWD or B6 background had no effect on male reproduction. Incompatibility(ies) not involving Prdm9^B6 also acts in the (PWD×B6)F₁ hybrids, since the correction of hybrid sterility by Prdm9^B6 deletion was not complete. Additions and subtractions of Prdm9 copies, as well as allelic replacements, improved meiotic progression and fecundity also in the progeny-producing reciprocal (B6×PWD)F₁ males. Moreover, an increased dosage of Prdm9 and reciprocal cross enhanced fertility of other sperm-carrying male hybrids, (PWD×B6-C3H.Prdm9)F₁, harboring another Prdm9 allele of M. m. domesticus origin. The levels of Prdm9 mRNA isoforms were similar in the prepubertal testes of all types of F₁ hybrids of PWD with B6 and B6-C3H.Prdm9 despite their different prospective fertility, but decreased to 53% after removal of Prdm9^B6. Therefore, the Prdm9^B6 allele probably takes part in posttranscriptional dominant-negative hybrid interaction(s) absent in the parental strains.

Disturbed gametogenesis in the progeny of two fertile parental forms is called hybrid sterility; it is an important part of reproductive barriers between species. The Dobzhansky-Muller model of incompatibilities explains reproductive isolation between species by incorrect interactions between genes. Hybrid sterility 1 (Hst1) is one of the genes causing meiotic arrest in F₁ male hybrids between certain Mus musculus musculus (e.g., the PWD strain) and M. m. domesticus (C57BL/6J etc.) mice. Hst1, the first mammalian candidate for a speciation gene, was identified as a meiotic PR/SET-domain gene, Prdm9, but the mechanism causing sterility has remained unknown. While the F₁ male offspring of C57BL/6J males and PWD females produce no sperm, the males from the reciprocal cross using PWD males and C57BL/6J females yield progeny. Here we show that the meiotic progress and fertility of hybrid males from both F₁ crosses improved by removal as well as overexpression of the C57BL/6J allele of Prdm9, suggesting that Prdm9 interactions not present in the parental species (incompatibilities) play a role in hybrid sterility. Furthermore, the Prdm9 dosage also controlled fecundity in other F₁ hybrids, indicating that this gene is an important regulator of mouse hybrid fertility.

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 F1 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 F1 autosomal background, and we measured a suite of male reproductive traits. Reproductive deficits were observed in most F1 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 F1 male sterility appears to be caused mainly by X-autosome incompatibilities, X-Y incompatibilities contribute to some aspects of sterility.

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.

Genome-wide architecture of reproductive isolation in a naturally occurring hybrid zone between Mus musculus musculus and M. m. domesticus

Studies of a hybrid zone between two house mouse subspecies (Mus musculus musculus and M. m. domesticus) along with studies using laboratory crosses reveal a large role for the X chromosome and multiple autosomal regions in reproductive isolation as a consequence of disrupted epistasis in hybrids. One limitation of previous work has been that most of the identified genomic regions have been large. The goal here is to detect and characterize precise genomic regions underlying reproductive isolation. We surveyed 1401 markers evenly spaced across the genome in 679 mice collected from two different transects. Comparisons between transects provide a means for identifying common patterns that likely reflect intrinsic incompatibilities. We used a genomic cline approach to identify patterns that correspond to epistasis. From both transects, we identified contiguous regions on the X chromosome in which markers were inferred to be involved in epistatic interactions. We then searched for autosomal regions showing the same patterns and found they constitute about 5% of autosomal markers. We discovered substantial overlap between these candidate regions underlying reproductive isolation and QTL for hybrid sterility identified in laboratory crosses. Analysis of gene content in these regions suggests a key role for several mechanisms, including the regulation of transcription, sexual conflict and sexual selection operating at both the postmating prezygotic and postzygotic stages of reproductive isolation. Taken together, these results indicate that speciation in two recently diverged (c. 0.5 Ma) house mouse subspecies is complex, involving many genes dispersed throughout the genome and associated with distinct functions.

Dissecting the genetic architecture of F1 hybrid sterility in house mice

Hybrid sterility as a postzygotic reproductive isolation mechanism has been studied for over 80 years, yet the first identifications of hybrid sterility genes in Drosophila and mouse are quite recent. To study the genetic architecture of F(1) hybrid sterility between young subspecies of house mouse Mus m. domesticus and M. m. musculus, we conducted QTL analysis of a backcross between inbred strains representing these two subspecies and probed the role of individual chromosomes in hybrid sterility using the intersubspecific chromosome substitution strains. We provide direct evidence that the asymmetry in male infertility between reciprocal crosses is conferred by the middle region of M. m. musculus Chr X, thus excluding other potential candidates such as Y, imprinted genes, and mitochondrial DNA. QTL analysis identified strong hybrid sterility loci on Chr 17 and Chr X and predicted a set of interchangeable autosomal loci, a subset of which is sufficient to activate the Dobzhansky-Muller incompatibility of the strong loci. Overall, our results indicate the oligogenic nature of F(1) hybrid sterility, which should be amenable to reconstruction by proper combination of chromosome substitution strains. Such a prefabricated model system should help to uncover the gene networks and molecular mechanisms underlying hybrid sterility.

Incompatibility between X chromosome factor and pericentric heterochromatic region causes lethality in hybrids between Drosophila melanogaster and its sibling species

The Dobzhansky-Muller model posits that postzygotic reproductive isolation results from the evolution of incompatible epistatic interactions between species: alleles that function in the genetic background of one species can cause sterility or lethality in the genetic background of another species. Progress in identifying and characterizing factors involved in postzygotic isolation in Drosophila has remained slow, mainly because Drosophila melanogaster, with all of its genetic tools, forms dead or sterile hybrids when crossed to its sister species, D. simulans, D. sechellia, and D. mauritiana. To circumvent this problem, we used chromosome deletions and duplications from D. melanogaster to map two hybrid incompatibility loci in F(1) hybrids with its sister species. We mapped a recessive factor to the pericentromeric heterochromatin of the X chromosome in D. simulans and D. mauritiana, which we call heterochromatin hybrid lethal (hhl), which causes lethality in F(1) hybrid females with D. melanogaster. As F(1) hybrid males hemizygous for a D. mauritiana (or D. simulans) X chromosome are viable, the lethality of deficiency hybrid females implies that a dominant incompatible partner locus exists on the D. melanogaster X. Using small segments of the D. melanogaster X chromosome duplicated onto the Y chromosome, we mapped a dominant factor that causes hybrid lethality to a small 24-gene region of the D. melanogaster X. We provide evidence suggesting that it interacts with hhl(mau). The location of hhl is consistent with the emerging theme that hybrid incompatibilities in Drosophila involve heterochromatic regions and factors that interact with the heterochromatin.

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.

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.

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.

The microevolution of V1r vomeronasal receptor genes in mice

Vomeronasal sensitivity is important for detecting intraspecific pheromonal cues as well as environmental odorants and is involved in mating, social interaction, and other daily activities of many vertebrates. Two large families of seven-transmembrane G-protein-coupled receptors, V1rs and V2rs, bind to various ligands to initiate vomeronasal signal transduction. Although the macroevolution of V1r and V2r genes has been well characterized throughout vertebrates, especially mammals, little is known about their microevolutionary patterns, which hampers a clear understanding of the evolutionary forces behind the rapid evolutionary turnover of V1r and V2r genes and the great diversity in receptor repertoire across species. Furthermore, the role of divergent vomeronasal perception in enhancing premating isolation and maintaining species identity has not been evaluated. Here we sequenced 44 V1r genes and 25 presumably neutral noncoding regions in 14 wild-caught mice belonging to Mus musculus and M. domesticus, two closely related species with strong yet incomplete reproductive isolation. We found that nucleotide changes in V1rs are generally under weak purifying selection and that only ∼5% of V1rs may have been subject to positive selection that promotes nonsynonymous substitutions. Consistent with the low functional constraints on V1rs, 18 of the 44 V1rs have null alleles segregating in one or both species. Together, our results demonstrate that, despite occasional actions of positive selection, the evolution of V1rs is in a large part shaped by purifying selection and random drift. These findings have broad implications for understanding the driving forces of rapid gene turnovers that are often observed in the evolution of large gene families.

Reinforcement selection acting on the European house mouse hybrid zone

Behavioural isolation may lead to complete speciation when partial postzygotic isolation acts in the presence of divergent-specific mate-recognition systems. These conditions exist where Mus musculus musculus and M. m. domesticus come into contact and hybridize. We studied two mate-recognition signal systems, based on urinary and salivary proteins, across a Central European portion of the mouse hybrid zone. Introgression of the genomic regions responsible for these signals: the major urinary proteins (MUPs) and androgen binding proteins (ABPs), respectively, was compared to introgression at loci assumed to be nearly neutral and those under selection against hybridization. The preference of individuals taken from across the zone regarding these signals was measured in Y mazes, and we develop a model for the analysis of the transition of such traits under reinforcement selection. The strongest assortative preferences were found in males for urine and females for ABP. Clinal analyses confirm nearly neutral introgression of an Abp locus and two loci closely linked to the Abp gene cluster, whereas two markers flanking the Mup gene region reveal unexpected introgression. Geographic change in the preference traits matches our reinforcement selection model significantly better than standard cline models. Our study confirms that behavioural barriers are important components of reproductive isolation between the house mouse subspecies.

A candidate subspecies discrimination system involving a vomeronasal receptor gene with different alleles fixed in M. m. domesticus and M. m. musculus

Assortative mating, a potentially efficient prezygotic reproductive barrier, may prevent loss of genetic potential by avoiding the production of unfit hybrids (i.e., because of hybrid infertility or hybrid breakdown) that occur at regions of secondary contact between incipient species. In the case of the mouse hybrid zone, where two subspecies of Mus musculus (M. m. domesticus and M. m. musculus) meet and exchange genes to a limited extent, assortative mating requires a means of subspecies recognition. We based the work reported here on the hypothesis that, if there is a pheromone sufficiently diverged between M. m. domesticus and M. m. musculus to mediate subspecies recognition, then that process must also require a specific receptor(s), also sufficiently diverged between the subspecies, to receive the signal and elicit an assortative mating response. We studied the mouse V1R genes, which encode a large family of receptors in the vomeronasal organ (VNO), by screening Perlegen SNP data and identified one, Vmn1r67, with 24 fixed SNP differences most of which (15/24) are nonsynonymous nucleotide substitutions between M. m. domesticus and M. m. musculus. We observed substantial linkage disequilibrium (LD) between Vmn1r67 and Abpa27, a mouse salivary androgen-binding protein gene that encodes a proteinaceous pheromone (ABP) capable of mediating assortative mating, perhaps in conjunction with its bound small lipophilic ligand. The LD we observed is likely a case of association rather than residual physical linkage from a very recent selective sweep, because an intervening gene, Vmn1r71, shows significant intra(sub)specific polymorphism but no inter(sub)specific divergence in its nucleotide sequence. We discuss alternative explanations of these observations, for example that Abpa27 and Vmn1r67 are coevolving as signal and receptor to reinforce subspecies hybridization barriers or that the unusually divergent Vmn1r67 allele was not a product of fast positive selection, but was derived from an introgressed allele, possibly from Mus spretus.

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.

Fine haplotype structure of a chromosome 17 region in the laboratory and wild mouse

Extensive linkage disequilibrium among classical laboratory strains represents an obstacle in the high-resolution haplotype mapping of mouse quantitative trait loci (QTL). To determine the potential of wild-derived mouse strains for fine QTL mapping, we constructed a haplotype map of a 250-kb region of the t-complex on chromosome 17 containing the Hybrid sterility 1 (Hst1) gene. We resequenced 33 loci from up to 80 chromosomes of five mouse (sub)species. Trans-species single-nucleotide polymorphisms (SNPs) were rare between Mus m. musculus (Mmmu) and Mus m. domesticus (Mmd). The haplotypes in Mmmu and Mmd differed and therefore strains from these subspecies should not be combined for haplotype-associated mapping. The haplotypes of t-chromosomes differed from all non-t Mmmu and Mmd haplotypes. Half of the SNPs and SN indels but only one of seven longer rearrangements found in classical laboratory strains were useful for haplotype mapping in the wild-derived M. m. domesticus. The largest Mmd haplotype block contained three genes of a highly conserved synteny. The lengths of the haplotype blocks deduced from 36 domesticus chromosomes were in tens of kilobases, suggesting that the wild-derived Mmd strains are suitable for fine interval-specific mapping.

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.

Genome-wide patterns of gene flow across a house mouse hybrid zone

Hybrid zones between closely related species or subspecies provide useful settings for studying the genetic architecture of speciation. Using markers distributed throughout the mouse genome, we use a hybrid zone between two recently diverged species of house mice (Mus musculus and Mus domesticus) as a natural mapping experiment to identify genomic regions that may be involved in reproductive isolation. Using cline analysis we document a nearly 50-fold variation in level of introgression among markers. Some markers have extremely narrow cline widths; these genomic regions may contribute to reproductive isolation. Biological processes associated with these narrow clines include physiological and immune responses to the environment as well as physiological and behavioral aspects of reproduction. Other autosomal markers exhibit asymmetrically broad clines, usually with high frequencies of M. domesticus alleles on the M. musculus side of the hybrid zone. These markers identify genome regions likely housing genes with alleles that are spreading from one species to the other. Biological processes associated with these wide clines include cell signaling, olfaction, and pheromone response. These processes play important roles in survival and reproduction, and associated genes are likely targets of selection. Patterns of linkage disequilibrium in the center of the hybrid zone suggest that isolation may be caused by multiple epistatic interactions between sets of genes. These data highlight the complex genetic architecture underlying speciation even at early stages of divergence and point to some of the biological processes that may govern this architecture.

Extensive additivity of gene expression differentiates subspecies of the house mouse

We have studied different subspecies of the house mouse and their reciprocal F(1) hybrids to estimate the within-locus mode of inheritance for subspecies differences in gene expression in three tissues (brain, liver, and testis) of male mice. This study investigates the mode of inheritance in crosses at a larger taxonomic distance than has been previously systematically investigated. We found the vast majority of transcripts to be additively expressed with only a few transcripts showing dominance or overdominance in expression, except for one direction of one cross, which showed large mis-expression in the testis. We suggest that, as time passes, more genes come to influence expression, and if there is no directional dominance, additivity becomes increasingly more likely, up to a threshold beyond which there is F(1) hybrid breakdown. Some previous studies on different organisms have found a large degree of dominance, commonly at shorter taxonomic differences. We surveyed these findings and show that the most consistent association exists between the amount of additivity detected in a study and the expression analysis method (in particular microarray platform), suggesting that at least some of the differences among studies might be methodological. Most studies agree with ours in that within-locus additivity seems to be general mode of inheritance for transcript expression. Differentially expressed transcripts identified in our screen among subspecies of house mice are candidate genes that could be involved in reproductive isolation between these subspecies.

Disruption of genetic interaction between two autosomal regions and the X chromosome causes reproductive isolation between mouse strains derived from different subspecies

Reproductive isolation that initiates speciation is likely caused by incompatibility among multiple loci in organisms belonging to genetically diverging populations. Laboratory C57BL/6J mice, which predominantly originated from Mus musculus domesticus, and a MSM/Ms strain derived from Japanese wild mice (M. m. molossinus, genetically close to M. m. musculus) are reproductively isolated. Their F1 hybrids are fertile, but successive intercrosses result in sterility. A consomic strain, C57BL/6J-ChrX(MSM), which carries the X chromosome of MSM/Ms in the C57BL/6J background, shows male sterility, suggesting a genetic incompatibility of the MSM/Ms X chromosome and other C57BL/6J chromosome(s). In this study, we conducted genomewide linkage analysis and subsequent QTL analysis using the sperm shape anomaly that is the major cause of the sterility of the C57BL/6J-ChrX(MSM) males. These analyses successfully detected significant QTL on chromosomes 1 and 11 that interact with the X chromosome. The introduction of MSM/Ms chromosomes 1 and 11 into the C57BL/6J-ChrX(MSM) background failed to restore the sperm-head shape, but did partially restore fertility. This result suggests that this genetic interaction may play a crucial role in the reproductive isolation between the two strains. A detailed analysis of the male sterility by intracytoplasmic sperm injection and zona-free in vitro fertilization demonstrated that the C57BL/6J-ChrX(MSM) spermatozoa have a defect in penetration through the zona pellucida of eggs.

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

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