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

The role of gene presence-absence variations on genetic incompatibility in Asian rice

Genetic incompatibilities are widespread between species. However, it remains unclear whether they all originated after population divergence as suggested by the Bateson-Dobzhansky-Muller model, and if not, what is their prevalence and distribution within populations. The gene presence-absence variations (PAVs) provide an opportunity for investigating gene-gene incompatibility. Here, we searched for the repulsion of coexistence between gene PAVs to identify the negative interaction of gene functions separately in two Oryza sativa subspecies. Many PAVs are involved in subspecies-specific negative epistasis and segregate at low-to-intermediate frequencies in focal subspecies but at low or high frequencies in the other subspecies. Incompatible PAVs are enriched in two functional groups, defense response and protein phosphorylation, which are associated with plant immunity and consistent with autoimmunity being a known mechanism of hybrid incompatibility in plants. Genes in the two enriched functional groups are older and seldom directly interact with each other. Instead, they interact with other younger gene PAVs with diverse functions. Our results illustrate the landscape of genetic incompatibility at gene PAVs in rice, where many incompatible pairs have already segregated as polymorphisms within subspecies, and many are novel negative interactions between older defense-related genes and younger genes with diverse functions.

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.

Variable hybridization between two Lake Tanganyikan cichlid species in recent secondary contact

Closely related taxa frequently exist in sympatry before the evolution of robust reproductive barriers, which can lead to substantial gene flow. Post-divergence gene flow can promote several disparate trajectories of divergence ranging from the erosion of distinctiveness and eventual collapse of the taxa to the strengthening of reproductive isolation. Among many relevant factors, understanding the demographic history of divergence (e.g. divergence time, extent of historical gene flow) can be particularly informative when examining contemporary gene flow between closely related taxa because this history can influence gene flow's prevalence and consequences. Here, we used genotyping-by-sequencing data to investigate speciation and contemporary hybridization in two closely related and sympatrically distributed Lake Tanganyikan cichlid species in the genus Petrochromis. Demographic modeling supported a speciation scenario involving divergence in isolation followed by secondary contact with bidirectional gene flow. Further investigation of this recent gene flow found evidence of ongoing hybridization between the species that varied in extent between different co-occurring populations. Relationships between abundance and the degree of admixture across populations suggest that the availability of conspecific mates may influence patterns of hybridization. These results, together with the observation that sets of recently diverged cichlid taxa are generally geographically separated in the lake, suggest that ongoing speciation in Lake Tanganyikan cichlids relies on initial spatial isolation. Additionally, the spatially heterogeneous patterns of admixture between the Petrochromis species illustrates the complexities of hybridization when species are in recent secondary contact.

Unraveling patterns of disrupted gene expression across a complex tissue

Whole tissue RNASeq is the standard approach for studying gene expression divergence in evolutionary biology and provides a snapshot of the comprehensive transcriptome for a given tissue. However, whole tissues consist of diverse cell types differing in expression profiles, and the cellular composition of these tissues can evolve across species. Here, we investigate the effects of different cellular composition on whole tissue expression profiles. We compared gene expression from whole testes and enriched spermatogenesis populations in two species of house mice, Mus musculus musculus and M. m. domesticus, and their sterile and fertile F1 hybrids, which differ in both cellular composition and regulatory dynamics. We found that cellular composition differences skewed expression profiles and differential gene expression in whole testes samples. Importantly, both approaches were able to detect large-scale patterns such as disrupted X chromosome expression, although whole testes sampling resulted in decreased power to detect differentially expressed genes. We encourage researchers to account for histology in RNASeq and consider methods that reduce sample complexity whenever feasible. Ultimately, we show that differences in cellular composition between tissues can modify expression profiles, potentially altering inferred gene ontological processes, insights into gene network evolution, and processes governing gene expression evolution.

The microbiome impacts host hybridization and speciation

Microbial symbiosis and speciation profoundly shape the composition of life's biodiversity. Despite the enormous contributions of these two fields to the foundations of modern biology, there is a vast and exciting frontier ahead for research, literature, and conferences to address the neglected prospects of merging their study. Here, we survey and synthesize exemplar cases of how endosymbionts and microbial communities affect animal hybridization and vice versa. We conclude that though the number of case studies remain nascent, the wide-ranging types of animals, microbes, and isolation barriers impacted by hybridization will likely prove general and a major new phase of study that includes the microbiome as part of the functional whole contributing to reproductive isolation. Though microorganisms were proposed to impact animal speciation a century ago, the weight of the evidence supporting this view has now reached a tipping point.

X-Linked Signature of Reproductive Isolation in Humans is Mirrored in a Howler Monkey Hybrid Zone

Reproductive isolation is a fundamental step in speciation. While sex chromosomes have been linked to reproductive isolation in many model systems, including hominids, genetic studies of the contribution of sex chromosome loci to speciation for natural populations are relatively sparse. Natural hybrid zones can help identify genomic regions contributing to reproductive isolation, like hybrid incompatibility loci, since these regions exhibit reduced introgression between parental species. Here, we use a primate hybrid zone (Alouatta palliata × Alouatta pigra) to test for reduced introgression of X-linked SNPs compared to autosomal SNPs. To identify X-linked sequence in A. palliata, we used a sex-biased mapping approach with whole-genome re-sequencing data. We then used genomic cline analysis with reduced-representation sequence data for parental A. palliata and A. pigra individuals and hybrids (n = 88) to identify regions with non-neutral introgression. We identified ~26 Mb of non-repetitive, putatively X-linked genomic sequence in A. palliata, most of which mapped collinearly to the marmoset and human X chromosomes. We found that X-linked SNPs had reduced introgression and an excess of ancestry from A. palliata as compared to autosomal SNPs. One outlier region with reduced introgression overlaps a previously described "desert" of archaic hominin ancestry on the human X chromosome. These results are consistent with a large role for the X chromosome in speciation across animal taxa and further, suggest shared features in the genomic basis of the evolution of reproductive isolation in primates.

The genomic consequences of hybridization

In the past decade, advances in genome sequencing have allowed researchers to uncover the history of hybridization in diverse groups of species, including our own. Although the field has made impressive progress in documenting the extent of natural hybridization, both historical and recent, there are still many unanswered questions about its genetic and evolutionary consequences. Recent work has suggested that the outcomes of hybridization in the genome may be in part predictable, but many open questions about the nature of selection on hybrids and the biological variables that shape such selection have hampered progress in this area. We synthesize what is known about the mechanisms that drive changes in ancestry in the genome after hybridization, highlight major unresolved questions, and discuss their implications for the predictability of genome evolution after hybridization.

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.

Age and Genetic Background Modify Hybrid Male Sterility in House Mice

Hybrid male sterility (HMS) contributes to reproductive isolation commonly observed among house mouse (Mus musculus) subspecies, both in the wild and in laboratory crosses. Incompatibilities involving specific Prdm9 alleles and certain Chromosome (Chr) X genotypes are known determinants of fertility and HMS, and previous work in the field has demonstrated that genetic background modifies these two major loci. We constructed hybrids that have identical genotypes at Prdm9 and identical X chromosomes, but differ widely across the rest of the genome. In each case, we crossed female PWK/PhJ mice representative of the M. m. musculus subspecies to males from a classical inbred strain representative of M. m. domesticus: 129S1/SvImJ, A/J, C57BL/6J, or DBA/2J. We detected three distinct trajectories of fertility among the hybrids using breeding experiments. The PWK129S1 males were always infertile. PWKDBA2 males were fertile, despite their genotypes at the major HMS loci. We also observed age-dependent changes in fertility parameters across multiple genetic backgrounds. The PWKB6 and PWKAJ males were always infertile before 12 weeks and after 35 weeks. However, some PWKB6 and PWKAJ males were transiently fertile between 12 and 35 weeks. This observation could resolve previous contradictory reports about the fertility of PWKB6. Taken together, these results point to multiple segregating HMS modifier alleles, some of which have age-related modes of action. The ultimate identification of these alleles and their age-related mechanisms will advance understanding both of the genetic architecture of HMS and of how reproductive barriers are maintained between house mouse subspecies.

Further resolution of the house mouse (Mus musculus) phylogeny by integration over isolation-with-migration histories

BACKGROUND

The three main subspecies of house mice, Mus musculus castaneus, Mus musculus domesticus, and Mus musculus musculus, are estimated to have diverged ~ 350-500KYA. Resolution of the details of their evolutionary history is complicated by their relatively recent divergence, ongoing gene flow among the subspecies, and complex demographic histories. Previous studies have been limited to some extent by the number of loci surveyed and/or by the scope of the method used. Here, we apply a method (IMa3) that provides an estimate of a population phylogeny while allowing for complex histories of gene exchange.

RESULTS

Results strongly support a topology with M. m. domesticus as sister to M. m. castaneus and M. m. musculus. In addition, we find evidence of gene flow between all pairs of subspecies, but that gene flow is most restricted from M. m. musculus into M. m. domesticus. Estimates of other key parameters are dependent on assumptions regarding generation time and mutation rate in house mice. Nevertheless, our results support previous findings that the effective population size, Ne, of M. m. castaneus is larger than that of the other two subspecies, that the three subspecies began diverging ~ 130 - 420KYA, and that the time between divergence events was short.

CONCLUSIONS

Joint demographic and phylogenetic analyses of genomic data provide a clearer picture of the history of divergence in house mice.

The importance of intrinsic postzygotic barriers throughout the speciation process

Intrinsic postzygotic barriers can play an important and multifaceted role in speciation, but their contribution is often thought to be reserved to the final stages of the speciation process. Here, we review how intrinsic postzygotic barriers can contribute to speciation, and how this role may change through time. We outline three major contributions of intrinsic postzygotic barriers to speciation. (i) reduction of gene flow: intrinsic postzygotic barriers can effectively reduce gene exchange between sympatric species pairs. We discuss the factors that influence how effective incompatibilities are in limiting gene flow. (ii) early onset of species boundaries via rapid evolution: intrinsic postzygotic barriers can evolve between recently diverged populations or incipient species, thereby influencing speciation relatively early in the process. We discuss why the early origination of incompatibilities is expected under some biological models, and detail how other (and often less obvious) incompatibilities may also serve as important barriers early on in speciation. (iii) reinforcement: intrinsic postzygotic barriers can promote the evolution of subsequent reproductive isolation through processes such as reinforcement, even between relatively recently diverged species pairs. We incorporate classic and recent empirical and theoretical work to explore these three facets of intrinsic postzygotic barriers, and provide our thoughts on recent challenges and areas in the field in which progress can be made. This article is part of the theme issue 'Towards the completion of speciation: the evolution of reproductive isolation beyond the first barriers'.

Reverse genetics reveals single gene of every candidate on Hybrid sterility, X Chromosome QTL 2 (Hstx2) are dispensable for spermatogenesis

F1 hybrid progenies between related subspecies often show hybrid sterility (HS) or inviability. HS is caused by failure of meiotic chromosome synapsis and sex body formation in house mouse. Previous studies identified two HS critical genomic regions named Hstx2 on Chr X and Hst1 on Chr 17 by murine forward genetic approaches. HS gene on Hst1 was reported to be Prdm9. Intersubspecific polymorphisms of Prdm9 induce HS in hybrids, and Prdm9 null mutation leads to sterility in the inbred strain. However, HS gene on Hstx2 remains unknown. Here, using knock-out studies, we showed that HS candidate genes on Hstx2 are not individually essential for spermatogenesis in B6 strain. We examined 12 genes on Hstx2: Ctag2, 4930447F04Rik, Mir743, Mir465d, Mir465c-2, Mir465b-1, Mir465c-1, Mir465, Gm1140, Gm14692, 4933436I01Rik, and Gm6812. These genes were expressed in adult testes, and showed intersubspecific polymorphisms on expressed regions. This first reverse genetic approach to identify HS gene on Hstx2 suggested that the loss of function of any one HS candidate gene does not cause complete sterility, unlike Prdm9. Thus, the mechanism(s) of HS by the HS gene on Hstx2 might be different from that of Prdm9.

Variable hybridization outcomes in trout are predicted by historical fish stocking and environmental context

Hybridization can profoundly affect the genomic composition and phenotypes of closely related species, and provides an opportunity to identify mechanisms that maintain reproductive isolation between species. Recent evidence suggests that hybridization outcomes within a species pair can vary across locations. However, we still do not know how variable outcomes of hybridization are across geographic replicates, and what mechanisms drive that variation. In this study, we described hybridization outcomes across 27 locations in the North Fork Shoshone River basin (Wyoming, USA) where native Yellowstone cutthroat trout and introduced rainbow trout co-occur. We used genomic data and hierarchical Bayesian models to precisely identify ancestry of hybrid individuals. Hybridization outcomes varied across locations. In some locations, only rainbow trout and advanced backcrossed hybrids towards rainbow trout were present, while trout in other locations had a broader range of ancestry, including both parental species and first-generation hybrids. Later-generation intermediate hybrids were rare relative to backcrossed hybrids and rainbow trout individuals. Using an individual-based simulation, we found that outcomes of hybridization in the North Fork Shoshone River basin deviate substantially from what we would expect under null expectations of random mating and no selection against hybrids. Since this deviation implies that some mechanisms of reproductive isolation function to maintain parental taxa and a diversity of hybrid types, we then modelled hybridization outcomes as a function of environmental variables and stocking history that are likely to affect prezygotic barriers to hybridization. Variables associated with history of fish stocking were the strongest predictors of hybridization outcomes, followed by environmental variables that might affect overlap in spawning time and location.

Premeiotic and meiotic failures lead to hybrid male sterility in the Anopheles gambiae complex

Hybrid male sterility (HMS) contributes to speciation by restricting gene flow between related taxa. Detailed cytological characterization of reproductive organs in hybrid males is important for identifying phenotypes that can help guide searches of speciation genes. To investigate possible cellular causes of HMS, we performed crosses between closely related species of the Anopheles gambiae complex: An. merus with An. gambiae or An. coluzzii. We demonstrate that HMS in African malaria mosquitoes involves two defects in the reciprocal crosses: a premeiotic arrest of germline stem cells in degenerate testes and a failure of the reductional meiotic division of primary spermatocytes in normal-like testes. The premeiotic arrest in degenerate testes of hybrids is accompanied by a strong suppression of meiotic and postmeiotic genes. Unlike pure species, sex chromosomes in normal-like testes of F1 hybrids are largely unpaired during meiotic prophase I and all chromosomes show various degrees of insufficient condensation. Instead of entering reductional division in meiosis I, primary spermatocytes prematurely undergo an equational mitotic division producing non-motile diploid sperm. Thus, our study identified cytogenetic errors in interspecies hybrids that arise during the early stages of postzygotic isolation.

Parallel pattern of differentiation at a genomic island shared between clinal and mosaic hybrid zones in a complex of cryptic seahorse lineages

Diverging semi-isolated lineages either meet in narrow clinal hybrid zones, or have a mosaic distribution associated with environmental variation. Intrinsic reproductive isolation is often emphasized in the former and local adaptation in the latter, although both reduce gene flow between groups. Rarely are these two patterns of spatial distribution reported in the same study system. Here, we report that the long-snouted seahorse Hippocampus guttulatus is subdivided into discrete panmictic entities by both types of hybrid zones. Along the European Atlantic coasts, a northern and a southern lineage meet in the southwest of France where they coexist in sympatry-i.e., in the same geographical zone-with little hybridization. In the Mediterranean Sea, two lineages have a mosaic distribution, associated with lagoon-like and marine habitats. A fifth lineage was identified in the Black Sea. Genetic homogeneity over large spatial scales contrasts with isolation maintained in sympatry or close parapatry at a fine scale. A high variation in locus-specific introgression rates provides additional evidence that partial reproductive isolation must be maintaining the divergence. We find that fixed differences between lagoon and marine populations in the Mediterranean Sea belong to the most differentiated SNPs between the two Atlantic lineages, against the genome-wide pattern of structure that mostly follow geography. These parallel outlier SNPs cluster on a single chromosome-wide island of differentiation. Since Atlantic lineages do not map to lagoon-sea habitat variation, genetic parallelism at the genomic island suggests a shared genetic barrier contributes to reproductive isolation in contrasting contexts-i.e., spatial versus ecological. We discuss how a genomic hotspot of parallel differentiation could have evolved and become associated both with space and with a patchy environment in a single study system.

Hybrid Sterility with Meiotic Metaphase Arrest in Intersubspecific Mouse Crosses

Although organisms belonging to different species and subspecies sometimes produce fertile offspring, a hallmark of the speciation process is reproductive isolation, characterized by hybrid sterility (HS) due to failure in gametogenesis. In mammals, HS is usually exhibited by males, the heterogametic sex. The phenotypic manifestations of HS are complex. The most frequently observed are abnormalities in both autosomal and sex chromosome interactions that are linked to meiotic prophase arrest or postmeiotic spermiogenesis aberrations and lead to defective or absent gametes. The aim of this study was to determine the HS phenotypes in intersubspecific F1 mice produced by matings between Mus musculus molossinus-derived strains and diverse Mus musculus domesticus-inbred laboratory mouse strains. Most of these crosses produced fertile F1 offspring. However, when female BALB/cJ (domesticus) mice were mated to male JF1/MsJ (molossinus) mice, the (BALBdomxJF1mol)F1 males were sterile, whereas the (JF1molxBALBdom)F1 males produced by the reciprocal crossings were fertile; thus the sterility phenotype was asymmetric. The sterile (BALBdomxJF1mol) F1 males exhibited a high rate of meiotic metaphase arrest with misaligned chromosomes, probably related to a high frequency of XY dissociation. Intriguingly, in the sterile (BALBdomxJF1mol)F1 males we observed aberrant allele-specific expression of several meiotic genes, that play critical roles in important meiotic events including chromosome pairing. Together, these observations of an asymmetrical HS phenotype in intersubspecific F1 males, probably owing to meiotic defects in the meiotic behavior of the XY chromosomes pair and possibly also transcriptional misregulation of meiotic genes, provide new models and directions for understanding speciation mechanisms in mammals.

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.

Influence of Effective Population Size on Genes under Varying Levels of Selection Pressure

The ratio of diversities at amino acid changing (nonsynonymous) and neutral (synonymous) sites (ω = πN/πS) is routinely used to measure the intensity of selection pressure. It is well known that this ratio is influenced by the effective population size (Ne) and selection coefficient (s). Here, we examined the effects of effective population size on ω by comparing protein-coding genes from Mus musculus castaneus and Mus musculus musculus-two mouse subspecies with different Ne. Our results revealed a positive relationship between the magnitude of selection intensity and the ω estimated for genes. For genes under high selective constraints, the ω estimated for the subspecies with small Ne (M. m. musculus) was three times higher than that observed for that with large Ne (M. m. castaneus). However, this difference was only 18% for genes under relaxed selective constraints. We showed that the observed relationship is qualitatively similar to the theoretical predictions. We also showed that, for highly expressed genes, the ω of M. m. musculus was 2.1 times higher than that of M.m. castaneus and this difference was only 27% for genes with low expression levels. These results suggest that the effect of effective population size is more pronounced in genes under high purifying selection. Hence the choice of genes is important when ω is used to infer the effective size of a population.

A two-locus hybrid incompatibility is widespread, polymorphic, and active in natural populations of Mimulus

Reproductive isolation, which is essential for the maintenance of species in sympatry, is often incomplete between closely related species. In these taxa, reproductive barriers must evolve within species, without being degraded by ongoing gene flow. To better understand this dynamic, we investigated the frequency and geographic distribution of alleles underlying a two-locus, hybrid lethality system between naturally hybridizing species of monkeyflower (Mimulus guttatus and M. nasutus). We found that M. guttatus typically carries hybrid lethality alleles at one locus (hl13) and M. nasutus typically carries hybrid lethality alleles at the other locus (hl14). As a result, natural hybrids carry incompatible alleles at both loci, and express hybrid lethality in later generations. We also discovered considerable polymorphism at both hl13 and hl14 within both species. For M. guttatus, polymorphism at both loci occurs within populations, meaning that incompatible allele pairings likely arise through intraspecific gene flow. Genetic variation at markers linked to hl13 and hl14 suggest that introgression from M. nasutus is the primary driver of this polymorphism within M. guttatus. Additionally, patterns of introgression at the two hybrid lethality loci suggest that natural selection eliminates incompatible allele pairings, suggesting that even weak reproductive barriers might promote genomic divergence between species.

Genetic basis to hybrid inviability is more complex than hybrid male sterility in Caenorhabditis nematodes

Hybrid male sterility often evolves before female sterility or inviability of hybrids, implying that the accumulation of divergence between separated lineages should lead hybrid male sterility to have a more polygenic basis. However, experimental evidence is mixed. Here, we use the nematodes Caenorhabditis remanei and C. latens to characterize the underlying genetic basis of asymmetric hybrid male sterility and hybrid inviability. We demonstrate that hybrid male sterility is consistent with a simple genetic basis, involving a single X-autosome incompatibility. We also show that hybrid inviability involves more genomic compartments, involving diverse nuclear-nuclear incompatibilities, a mito-nuclear incompatibility, and maternal effects. These findings demonstrate that male sensitivity to genetic perturbation may be genetically simple compared to hybrid inviability in Caenorhabditis and motivates tests of generality for the genetic architecture of hybrid incompatibility across the breadth of phylogeny.

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.

Rapid neo-sex chromosome evolution and incipient speciation in a major forest pest

Genome evolution is predicted to be rapid following the establishment of new (neo) sex chromosomes, but it is not known if neo-sex chromosome evolution plays an important role in speciation. Here we combine extensive crossing experiments with population and functional genomic data to examine neo-XY chromosome evolution and incipient speciation in the mountain pine beetle. We find a broad continuum of intrinsic incompatibilities in hybrid males that increase in strength with geographic distance between reproductively isolated populations. This striking progression of reproductive isolation is coupled with extensive gene specialization, natural selection, and elevated genetic differentiation on both sex chromosomes. Closely related populations isolated by hybrid male sterility also show fixation of alternative neo-Y haplotypes that differ in structure and male-specific gene content. Our results suggest that neo-sex chromosome evolution can drive rapid functional divergence between closely related populations irrespective of ecological drivers of divergence.

The evolution of new sex chromosomes potentially generates reproductive isolation. Here, Bracewell et al. combine crossing experiments with population and functional genomics to characterize neo-sex chromosome evolution and incipient speciation in the mountain pine beetle, Dendroctonus ponderosae.

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.

The Composite Regulatory Basis of the Large X-Effect in Mouse Speciation

The disruption of meiotic sex chromosome inactivation (MSCI) has been proposed to be a major developmental mechanism underlying the rapid evolution of hybrid male sterility. We tested this idea by analyzing cell-specific gene expression across spermatogenesis in two lineages of house mice and their sterile and fertile reciprocal hybrids. We found pervasive disruption of sex chromosome gene expression in sterile hybrids at every stage of spermatogenesis. Failure of MSCI was developmentally preceded by increased silencing of autosomal genes, supporting the hypothesis that divergence at the hybrid incompatibility gene, Prdm9, results in increased rates of autosomal asynapsis which in turn triggers widespread silencing of unsynapsed chromatin. We also detected opposite patterns of postmeiotic overexpression or hyper-repression of the sex chromosomes in reciprocal hybrids, supporting the hypothesis that genomic conflict has driven functional divergence that leads to deleterious X-Y dosage imbalances in hybrids. Our developmental timeline also exposed more subtle patterns of mitotic misregulation on the X chromosome, a previously undocumented stage of spermatogenic disruption in this cross. These results indicate that multiple hybrid incompatibilities have converged on a common regulatory phenotype, the disrupted expression of the sex chromosomes during spermatogenesis. Collectively, these data reveal a composite regulatory basis to hybrid male sterility in mice that helps resolve the mechanistic underpinnings of the well-documented large X-effect in mice speciation. We propose that the inherent sensitivity of spermatogenesis to X-linked regulatory disruption has the potential to be a major driver of reproductive isolation in species with chromosomal sex determination.

Whole Genome Sequence of Two Wild-Derived Mus musculus domesticus Inbred Strains, LEWES/EiJ and ZALENDE/EiJ, with Different Diploid Numbers

Wild-derived mouse inbred strains are becoming increasingly popular for complex traits analysis, evolutionary studies, and systems genetics. Here, we report the whole-genome sequencing of two wild-derived mouse inbred strains, LEWES/EiJ and ZALENDE/EiJ, of Mus musculus domesticus origin. These two inbred strains were selected based on their geographic origin, karyotype, and use in ongoing research. We generated 14× and 18× coverage sequence, respectively, and discovered over 1.1 million novel variants, most of which are private to one of these strains. This report expands the number of wild-derived inbred genomes in the Mus genus from six to eight. The sequence variation can be accessed via an online query tool; variant calls (VCF format) and alignments (BAM format) are available for download from a dedicated ftp site. Finally, the sequencing data have also been stored in a lossless, compressed, and indexed format using the multi-string Burrows-Wheeler transform. All data can be used without restriction.

What, if anything, are hybrids: enduring truths and challenges associated with population structure and gene flow

Hybridization is a potent evolutionary process that can affect the origin, maintenance, and loss of biodiversity. Because of its ecological and evolutionary consequences, an understanding of hybridization is important for basic and applied sciences, including conservation biology and agriculture. Herein, we review and discuss ideas that are relevant to the recognition of hybrids and hybridization. We supplement this discussion with simulations. The ideas we present have a long history, particularly in botany, and clarifying them should have practical consequences for managing hybridization and gene flow in plants. One of our primary goals is to illustrate what we can and cannot infer about hybrids and hybridization from molecular data; in other words, we ask when genetic analyses commonly used to study hybridization might mislead us about the history or nature of gene flow and selection. We focus on patterns of variation when hybridization is recent and populations are polymorphic, which are particularly informative for applied issues, such as contemporary hybridization following recent ecological change. We show that hybridization is not a singular process, but instead a collection of related processes with variable outcomes and consequences. Thus, it will often be inappropriate to generalize about the threats or benefits of hybridization from individual studies, and at minimum, it will be important to avoid categorical thinking about what hybridization and hybrids are. We recommend potential sampling and analytical approaches that should help us confront these complexities of hybridization.

Gene regulation and speciation in house mice

One approach to understanding the process of speciation is to characterize the genetic architecture of post-zygotic isolation. As gene regulation requires interactions between loci, negative epistatic interactions between divergent regulatory elements might underlie hybrid incompatibilities and contribute to reproductive isolation. Here, we take advantage of a cross between house mouse subspecies, where hybrid dysfunction is largely unidirectional, to test several key predictions about regulatory divergence and reproductive isolation. Regulatory divergence between Mus musculus musculus and M. m. domesticus was characterized by studying allele-specific expression in fertile hybrid males using mRNA-sequencing of whole testes. We found extensive regulatory divergence between M. m. musculus and M. m. domesticus, largely attributable to cis-regulatory changes. When both cis and trans changes occurred, they were observed in opposition much more often than expected under a neutral model, providing strong evidence of widespread compensatory evolution. We also found evidence for lineage-specific positive selection on a subset of genes related to transcriptional regulation. Comparisons of fertile and sterile hybrid males identified a set of genes that were uniquely misexpressed in sterile individuals. Lastly, we discovered a nonrandom association between these genes and genes showing evidence of compensatory evolution, consistent with the idea that regulatory interactions might contribute to Dobzhansky-Muller incompatibilities and be important in speciation.

Speciation and reduced hybrid female fertility in house mice

In mammals, intrinsic postzygotic isolation has been well studied in males but has been less studied in females, despite the fact that female gametogenesis and pregnancy provide arenas for hybrid sterility or inviability that are absent in males. Here, we asked whether inviability or sterility is observed in female hybrids of Mus musculus domesticus and M. m. musculus, taxa which hybridize in nature and for which male sterility has been well characterized. We looked for parent-of-origin growth phenotypes by measuring adult body weights in F1 hybrids. We evaluated hybrid female fertility by crossing F1 females to a tester male and comparing multiple reproductive parameters between intrasubspecific controls and intersubspecific hybrids. Hybrid females showed no evidence of parent-of-origin overgrowth or undergrowth, providing no evidence for reduced viability. However, hybrid females had smaller litter sizes, reduced embryo survival, fewer ovulations, and fewer small follicles relative to controls. Significant variation in reproductive parameters was seen among different hybrid genotypes, suggesting that hybrid incompatibilities are polymorphic within subspecies. Differences in reproductive phenotypes in reciprocal genotypes were observed and are consistent with cyto-nuclear incompatibilities or incompatibilities involving genomic imprinting. These findings highlight the potential importance of reduced hybrid female fertility in the early stages of speciation.

Ultrasonic Vocalizations of Male Mice Differ among Species and Females Show Assortative Preferences for Male Calls

Male house mice (Mus musculus) emit ultrasonic vocalizations (USVs) during courtship, which attract females, and we aimed to test whether females use these vocalizations for species or subspecies recognition of potential mates. We recorded courtship USVs of males from different Mus species, Mus musculus subspecies, and populations (F1 offspring of wild-caught Mus musculus musculus, Mus musculus domesticus (and F1 hybrid crosses), and Mus spicilegus), and we conducted playback experiments to measure female preferences for male USVs. Male vocalizations contained at least seven distinct syllable types, whose frequency of occurrence varied among species, subspecies, and populations. Detailed analyses of multiple common syllable types indicated that Mus musculus and Mus spicilegus could be discriminated based on spectral and temporal characteristics of their vocalizations, and populations of Mus musculus were also distinctive regardless of the classification model used. Females were able to discriminate USVs from different species, and showed assortative preferences for conspecific males. We found no evidence that females discriminate USVs of males from a different subspecies or separate populations of the same species, even though our spectral analyses identified acoustic features that differ between species, subspecies, and populations of the same species. Our results provide the first comparison of USVs between Mus species or between Mus musculus subspecies, and the first evidence that male USVs potentially facilitate species recognition.

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.

Hybrid male sterility between Drosophila willistoni species is caused by male failure to transfer sperm during copulation

BACKGROUND

The biological concept of species stresses the importance of understanding what mechanisms maintain species reproductively isolated from each other. Often such mechanisms are divided into premating and postmating, with the latest being the result of either prezygotic or postzygotic isolation barriers. Drosophila willistoni quechua and Drosophila willistoni willistoni are two subspecies that experience reproductive isolation. When a D. w. quechua female is crossed with a D. w. willistoni male, the hybrid males (F1QW) are unable to father progeny; however, the reciprocal cross produces fertile hybrids. Thus, the mechanism of isolation is unidirectional hybrid male sterility. However, the sterile F1QW males contain large amounts of motile sperm. Here we explore whether pre-copulatory or post-copulatory pre-zygotic mechanisms serve as major deterrents in the ability of F1QW males to father progeny.

RESULTS

Comparisons of parental and hybrid males copulation durations showed no significant reduction in copulation duration of F1QW males. Interrupted copulations of the parental species confirmed that sperm transfer occurs before the minimum copulation duration registered for F1QW males. However, we found that when females mate with F1QW males, sperm is not present inside the female storage organs and that the lack of sperm in storage is due to failure to transfer sperm rather than spillage or active sperm dumping by females.

CONCLUSIONS

Sterility of F1QW hybrid males is primarily driven by their inability to transfer sperm during copulation.

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.

Analysis of intestinal microbiota in hybrid house mice reveals evolutionary divergence in a vertebrate hologenome

Recent evidence suggests that natural selection operating on hosts to maintain their microbiome contributes to the emergence of new species, that is, the ‘hologenomic basis of speciation’. Here we analyse the gut microbiota of two house mice subspecies, Mus musculus musculus and M. m. domesticus, across their Central European hybrid zone, in addition to hybrids generated in the lab. Hybrid mice display widespread transgressive phenotypes (that is, exceed or fall short of parental values) in a variety of measures of bacterial community structure, which reveals the importance of stabilizing selection operating on the intestinal microbiome within species. Further genetic and immunological analyses reveal genetic incompatibilities, aberrant immune gene expression and increased intestinal pathology associated with altered community structure among hybrids. These results provide unique insight into the consequences of evolutionary divergence in a vertebrate ‘hologenome’, which may be an unrecognized contributing factor to reproductive isolation in this taxonomic group.

Animal hosts and their associated microbes are largely the outcome of coevolution. Here, the authors show differences in the intestinal microbiome of hybrids compared with pure species of house mice, which suggests that host–microbiome interactions contribute to the evolution of host species.

Genomic evidence of geographically widespread effect of gene flow from polar bears into brown bears

Polar bears are an arctic, marine adapted species that is closely related to brown bears. Genome analyses have shown that polar bears are distinct and genetically homogeneous in comparison to brown bears. However, these analyses have also revealed a remarkable episode of polar bear gene flow into the population of brown bears that colonized the Admiralty, Baranof and Chichagof islands (ABC islands) of Alaska. Here, we present an analysis of data from a large panel of polar bear and brown bear genomes that includes brown bears from the ABC islands, the Alaskan mainland and Europe. Our results provide clear evidence that gene flow between the two species had a geographically wide impact, with polar bear DNA found within the genomes of brown bears living both on the ABC islands and in the Alaskan mainland. Intriguingly, while brown bear genomes contain up to 8.8% polar bear ancestry, polar bear genomes appear to be devoid of brown bear ancestry, suggesting the presence of a barrier to gene flow in that direction.

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.

Genome-wide mapping in a house mouse hybrid zone reveals hybrid sterility loci and Dobzhansky-Muller interactions

Mapping hybrid defects in contact zones between incipient species can identify genomic regions contributing to reproductive isolation and reveal genetic mechanisms of speciation. The house mouse features a rare combination of sophisticated genetic tools and natural hybrid zones between subspecies. Male hybrids often show reduced fertility, a common reproductive barrier between incipient species. Laboratory crosses have identified sterility loci, but each encompasses hundreds of genes. We map genetic determinants of testis weight and testis gene expression using offspring of mice captured in a hybrid zone between M. musculus musculus and M. m. domesticus. Many generations of admixture enables high-resolution mapping of loci contributing to these sterility-related phenotypes. We identify complex interactions among sterility loci, suggesting multiple, non-independent genetic incompatibilities contribute to barriers to gene flow in the hybrid zone.

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.

X-y interactions underlie sperm head abnormality in hybrid male house mice

The genetic basis of hybrid male sterility in house mice is complex, highly polygenic, and strongly X linked. Previous work suggested that there might be interactions between the Mus musculus musculus X and the M. m. domesticus Y with a large negative effect on sperm head morphology in hybrid males with an F1 autosomal background. To test this, we introgressed the M. m. domesticus Y onto a M. m. musculus background and measured the change in sperm morphology, testis weight, and sperm count across early backcross generations and in 11th generation backcross males in which the opportunity for X-autosome incompatibilities is effectively eliminated. We found that abnormality in sperm morphology persists in M. m. domesticus Y introgression males, and that this phenotype is rescued by M. m. domesticus introgressions on the X chromosome. In contrast, the severe reductions in testis weight and sperm count that characterize F1 males were eliminated after one generation of backcrossing. These results indicate that X-Y incompatibilities contribute specifically to sperm morphology. In contrast, X-autosome incompatibilities contribute to low testis weight, low sperm count, and sperm morphology. Restoration of normal testis weight and sperm count in first generation backcross males suggests that a small number of complex incompatibilities between loci on the M. m. musculus X and the M. m. domesticus autosomes underlie F1 male sterility. Together, these results provide insight into the genetic architecture of F1 male sterility and help to explain genome-wide patterns of introgression across the house mouse hybrid zone.

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.

X chromosome control of meiotic chromosome synapsis in mouse inter-subspecific hybrids

Hybrid sterility (HS) belongs to reproductive isolation barriers that safeguard the integrity of species in statu nascendi. Although hybrid sterility occurs almost universally among animal and plant species, most of our current knowledge comes from the classical genetic studies on Drosophila interspecific crosses or introgressions. With the house mouse subspecies Mus m. musculus and Mus m. domesticus as a model, new research tools have become available for studies of the molecular mechanisms and genetic networks underlying HS. Here we used QTL analysis and intersubspecific chromosome substitution strains to identify a 4.7 Mb critical region on Chromosome X (Chr X) harboring the Hstx2 HS locus, which causes asymmetrical spermatogenic arrest in reciprocal intersubspecific F1 hybrids. Subsequently, we mapped autosomal loci on Chrs 3, 9 and 13 that can abolish this asymmetry. Combination of immunofluorescent visualization of the proteins of synaptonemal complexes with whole-chromosome DNA FISH on pachytene spreads revealed that heterosubspecific, unlike consubspecific, homologous chromosomes are predisposed to asynapsis in F1 hybrid male and female meiosis. The asynapsis is under the trans- control of Hstx2 and Hst1/Prdm9 hybrid sterility genes in pachynemas of male but not female hybrids. The finding concurred with the fertility of intersubpecific F1 hybrid females homozygous for the Hstx2^Mmm allele and resolved the apparent conflict with the dominance theory of Haldane's rule. We propose that meiotic asynapsis in intersubspecific hybrids is a consequence of cis-acting mismatch between homologous chromosomes modulated by the trans-acting Hstx2 and Prdm9 hybrid male sterility genes.

Genomes of newly emerging species restrict their gene exchange with related taxa in order to secure integrity. Hybrid sterility is one of the reproductive isolation mechanisms restricting gene flow between closely related, sexually reproducing organisms. We showed that hybrid sterility between two closely related mouse subspecies is executed by a failure of meiotic synapsis of orthologous chromosomes in F1 hybrid males. The asynapsis of orthologous chromosomes occurred in meiosis of male and female hybrids, though only males were sterile due to trans-acting male-specific hybrid sterility genes. We located one of the two major hybrid sterility genes to a 4.7 Mb interval on Chromosome X, showed that it controls male sterility by modulating the extent of meiotic asynapsis and using the inter-subspecific chromosome substitution strains we refuted the simple interpretation of dominance theory of Haldane's rule. A new working hypothesis posits male sterility of mouse inter-subsubspecific F1 hybrids as a consequence of meiotic chromosome asynapsis caused by the cis-acting mismatch between orthologous chromosomes modulated by the trans-acting hybrid male sterility genes.

Sexual selection against natural hybrids may contribute to reinforcement in a house mouse hybrid zone

Sexual selection may hinder gene flow across contact zones when hybrid recognition signals are discriminated against. We tested this hypothesis in a unimodal hybrid zone between Mus musculus musculus and Mus musculus domesticus where a pattern of reinforcement was described and lower hybrid fitness documented. We presented mice from the border of the hybrid zone with a choice between opposite sex urine from the same subspecies versus hybrids sampled in different locations across the zone. While no preference was evidenced in domesticus mice, musculus males discriminated in favour of musculus signals and against hybrid signals. Remarkably, the pattern of hybrid unattractiveness did not vary across the hybrid zone. Moreover, allopatric populations tested in the same conditions did not discriminate against hybrid signals, indicating character displacement for signal perception or preference. Finally, habituation-discrimination tests assessing similarities between signals pointed out that hybrid signals differed from the parental ones. Overall, our results suggest that perception of hybrids as unattractive has evolved in border populations of musculus after the secondary contact with domesticus. We discuss the mechanisms involved in hybrid unattractiveness, and the potential impact of asymmetric sexual selection on the hybrid zone dynamics and gene flow between the two subspecies.

Wild mice as bountiful resources of novel genetic variants for quantitative traits

Most traits of biological importance, including traits for human complex diseases (e.g., obesity and diabetes), are continuously distributed. These complex or quantitative traits are controlled by multiple genetic loci called QTLs (quantitative trait loci), environments and their interactions. The laboratory mouse has long been used as a pilot animal model for understanding the genetic architecture of quantitative traits. Next-generation sequencing analyses and genome-wide SNP (single nucleotide polymorphism) analyses of mouse genomes have revealed that classical inbred strains commonly used throughout the world are derived from a few fancy mice with limited and non-randomly distributed genetic diversity that occurs in nature and also indicated that their genomes are predominantly Mus musculus domesticus in origin. Many QTLs for a huge variety of traits have so far been discovered from a very limited gene pool of classical inbred strains. However, wild M. musculus mice consisting of five subspecies widely inhabit areas all over the world, and hence a number of novel QTLs may still lie undiscovered in gene pools of the wild mice. Some of the QTLs are expected to improve our understanding of human complex diseases. Using wild M. musculus subspecies in Asia as examples, this review illustrates that wild mice are untapped natural resources for valuable QTL discovery.

Rapid male-specific regulatory divergence and down regulation of spermatogenesis genes in Drosophila species hybrids

In most crosses between closely related species of Drosophila, the male hybrids are sterile and show postmeiotic abnormalities. A series of gene expression studies using genomic approaches have found significant down regulation of postmeiotic spermatogenesis genes in sterile male hybrids. These results have led some to suggest a direct relationship between down regulation in gene expression and hybrid sterility. An alternative explanation to a cause-and-effect relationship between misregulation of gene expression and male sterility is rapid divergence of male sex regulatory elements leading to incompatible interactions in an interspecies hybrid genome. To test the effect of regulatory divergence in spermatogenesis gene expression, we isolated 35 fertile D. simulans strains with D. mauritiana introgressions in either the X, second or third chromosome. We analyzed gene expression in these fertile hybrid strains for a subset of spermatogenesis genes previously reported as significantly under expressed in sterile hybrids relative to D. simulans. We found that fertile autosomal introgressions can cause levels of gene down regulation similar to that of sterile hybrids. We also found that X chromosome heterospecific introgressions cause significantly less gene down regulation than autosomal introgressions. Our results provide evidence that rapid male sex gene regulatory divergence can explain misexpression of spermatogenesis genes in hybrids.

Low fertility of wild hybrid male flycatchers despite recent divergence

Postzygotic isolation may be important for maintaining species boundaries, particularly when premating barriers are incomplete. Little is known about the course of events leading from minor environmental mismatches affecting hybrid fitness to severe genetic incompatibilities causing sterility or inviability. We investigated whether reduced reproductive success of hybrid males was caused by suboptimal sperm traits or by more severe genetic incompatibilities in a hybrid zone of pied (Ficedula hypoleuca) and collared flycatchers (F. albicollis) on the island of Öland, Sweden. About 4 per cent hybridization is observed in this population and all female hybrids are sterile. We found no sperm in the ejaculates of most sampled hybrid males, and sperm with abnormal morphology in two hybrids. Furthermore, none of the hybrids sired any offspring because of high levels of hatching failure and extra-pair paternity in their nests. These results from a natural hybrid zone suggest that the spermatogenesis of hybrid males may become disrupted despite little genetic divergence between the parental species.

Fertilisation and early developmental barriers to hybridisation in field crickets

BACKGROUND

Post-mating interactions between the reproductive traits and gametes of mating individuals and among their genes within zygotes are invariably complex, providing multiple opportunities for reproduction to go awry. These interactions have the potential to act as barriers to gene flow between species, and may be important in the process of speciation. There are multiple post-mating barriers to interbreeding between the hybridising field crickets Gryllus bimaculatus and G. campestris. Female G. bimaculatus preferentially store sperm from conspecific males when mated to both conspecific and heterospecific partners. Additionally, conspecific males sire an even greater proportion of offspring than would be predicted from their sperm's representation in the spermatheca. The nature of these post-sperm-storage barriers to hybridisation are unknown. We use a fluorescent staining technique to determine whether barriers occur prior to, or during embryo development.

RESULTS

We show that eggs laid by G. bimaculatus females mated to G. campestris males are less likely to begin embryogenesis than eggs from conspecific mating pairs. Of the eggs that are successfully fertilised and start to develop, those from heterospecific mating pairs are more likely to arrest early, prior to blastoderm formation. We find evidence for bimodal variation among egg clutches in the number of developing embryos that subsequently arrest, indicating that there is genetic variation for incompatibility between mating individuals. In contrast to the pattern of early embryonic mortality, those hybrids reaching advanced stages of embryogenesis have survival rates equal to that of embryos from conspecific mating pairs.

CONCLUSIONS

Post-sperm-storage barriers to hybridisation show evidence of genetic polymorphism. They are sufficiently large, that if the species interbreed where they are sympatric, these barriers could play a role in the maintenance of reproductive isolation between them. The number of eggs that fail to develop represents a substantial cost of hybridization to G. bimaculatus females, and this cost could reinforce the evolution of barriers occurring earlier in the reproductive process.

Geographically multifarious phenotypic divergence during speciation

Speciation is an important evolutionary process that occurs when barriers to gene flow evolve between previously panmictic populations. Although individual barriers to gene flow have been studied extensively, we know relatively little regarding the number of barriers that isolate species or whether these barriers are polymorphic within species. Herein, we use a series of field and lab experiments to quantify phenotypic divergence and identify possible barriers to gene flow between the butterfly species Lycaeides idas and Lycaeides melissa. We found evidence that L. idas and L. melissa have diverged along multiple phenotypic axes. Specifically, we identified major phenotypic differences in female oviposition preference and diapause initiation, and more moderate divergence in mate preference. Multiple phenotypic differences might operate as barriers to gene flow, as shown by correlations between genetic distance and phenotypic divergence and patterns of phenotypic variation in admixed Lycaeides populations. Although some of these traits differed primarily between species (e.g., diapause initiation), several traits also varied among conspecific populations (e.g., male mate preference and oviposition preference).

Sex-specific clines support incipient speciation in a common European mammal

Hybrid zones provide excellent opportunities to study processes and mechanisms underlying reproductive isolation and speciation. Here we investigated sex-specific clines of molecular markers in hybrid zones of morphologically cryptic yet genetically highly-diverged evolutionary lineages of the European common vole (Microtus arvalis). We analyzed the position and width of four secondary contact zones along three independent transects in the region of the Alps using maternally (mitochondrial DNA) and paternally (Y-chromosome) inherited genetic markers. Given male-biased dispersal in the common vole, a selectively neutral secondary contact would show broader paternal marker clines than maternal ones. In a selective case, for example, involving a form of Haldane's rule, Y-chromosomal clines would not be expected to be broader than maternal markers because they are transmitted by the heterogametic sex and thus gene flow would be restricted. Consistent with the selective case, paternal clines were significantly narrower or at most equal in width to maternal clines in all contact zones. In addition, analyses using maximum likelihood cline-fitting detected a shift of paternal relative to maternal clines in three of four contact zones. These patterns suggest that processes at the contact zones in the common vole are not selectively neutral, and that partial reproductive isolation is already established between these evolutionary lineages. We conclude that hybrid zone movement, sexual selection and/or genetic incompatibilities are likely associated with an unusual unidirectional manifestation of Haldane's rule in this common European mammal.

Genetic disruption of the copulatory plug in mice leads to severely reduced fertility

Seminal fluid proteins affect fertility at multiple stages in reproduction. In many species, a male's ejaculate coagulates to form a copulatory plug. Although taxonomically widespread, the molecular details of plug formation remain poorly understood, limiting our ability to manipulate the structure and understand its role in reproduction. Here I show that male mice knockouts for transglutaminase IV (Tgm4) fail to form a copulatory plug, demonstrating that this gene is necessary for plug formation and lending a powerful new genetic tool to begin characterizing plug function. Tgm4 knockout males show normal sperm count, sperm motility, and reproductive morphology. However, very little of their ejaculate migrates into the female's reproductive tract, suggesting the plug prevents ejaculate leakage. Poor ejaculate migration leads to a reduction in the proportion of oocytes fertilized. However, Tgm4 knockout males fertilized between 3-11 oocytes, which should be adequate for a normal litter. Nevertheless, females mated to Tgm4 knockout males for approximately 14 days were significantly less likely to give birth to a litter compared to females mated to wild-type males. Therefore, it appears that the plug also affects post-fertilization events such as implantation and/or gestation. This study shows that a gene influencing the viscosity of seminal fluid has a major influence on male fertility.