Genetics and evolution of hybrid male sterility in house mice

Inference of selective forces on house mouse genomes during secondary contact in East Asia

The house mouse (Mus musculus), which is commensal to humans, has spread globally via human activities, leading to secondary contact between genetically divergent subspecies. This pattern of genetic admixture can provide insights into the selective forces at play in this well-studied model organism. Our analysis of 163 house mouse genomes, with a particular focus on East Asia, revealed substantial admixture between the subspecies castaneus and musculus, particularly in Japan and southern China. We revealed, despite the different level of autosomal admixture among regions, that all Y Chromosomes in the East Asian samples belonged to the musculus-type haplogroup, potentially explained by genomic conflict under sex-ratio distortion owing to varying copy numbers of ampliconic genes on sex chromosomes, Slx and Sly Our computer simulations, designed to replicate the observed scenario, show that the preferential fixation of musculus-type Y Chromosomes can be achieved with a slight increase in the male-to-female birth ratio. We also investigated the influence of selection on the posthybridization of the subspecies castaneus and musculus in Japan. Even though the genetic background of most Japanese samples closely resembles the subspecies musculus, certain genomic regions overrepresented the castaneus-like genetic components, particularly in immune-related genes. Furthermore, a large genomic block (∼2 Mbp) containing a vomeronasal/olfactory receptor gene cluster predominantly harbored castaneus-type haplotypes in the Japanese samples, highlighting the crucial role of olfaction-based recognition in shaping hybrid genomes.

Into the Wild: A novel wild-derived inbred strain resource expands the genomic and phenotypic diversity of laboratory mouse models

The laboratory mouse has served as the premier animal model system for both basic and preclinical investigations for over a century. However, laboratory mice capture only a subset of the genetic variation found in wild mouse populations, ultimately limiting the potential of classical inbred strains to uncover phenotype-associated variants and pathways. Wild mouse populations are reservoirs of genetic diversity that could facilitate the discovery of new functional and disease-associated alleles, but the scarcity of commercially available, well-characterized wild mouse strains limits their broader adoption in biomedical research. To overcome this barrier, we have recently developed, sequenced, and phenotyped a set of 11 inbred strains derived from wild-caught Mus musculus domesticus. Each of these "Nachman strains" immortalizes a unique wild haplotype sampled from one of five environmentally distinct locations across North and South America. Whole genome sequence analysis reveals that each strain carries between 4.73-6.54 million single nucleotide differences relative to the GRCm39 mouse reference, with 42.5% of variants in the Nachman strain genomes absent from current classical inbred mouse strain panels. We phenotyped the Nachman strains on a customized pipeline to assess the scope of disease-relevant neurobehavioral, biochemical, physiological, metabolic, and morphological trait variation. The Nachman strains exhibit significant inter-strain variation in >90% of 1119 surveyed traits and expand the range of phenotypic diversity captured in classical inbred strain panels. These novel wild-derived inbred mouse strain resources are set to empower new discoveries in both basic and preclinical research.

Different complex regulatory phenotypes underlie hybrid male sterility in divergent rodent crosses

Hybrid incompatibilities are a critical component of species barriers and may arise due to negative interactions between divergent regulatory elements in parental species. We used a comparative approach to identify common themes in the regulatory phenotypes associated with hybrid male sterility in two divergent rodent crosses, dwarf hamsters and house mice. We investigated three potential characteristic regulatory phenotypes in hybrids including the propensity towards over or underexpression relative to parental species, the influence of developmental stage on the extent of misexpression, and the role of the sex chromosomes on misexpression phenotypes. In contrast to near pervasive overexpression in hybrid house mice, we found that misexpression in hybrid dwarf hamsters was dependent on developmental stage. In both house mouse and dwarf hamster hybrids, however, misexpression increased with the progression of spermatogenesis, although to varying extents and with potentially different consequences. In both systems, we detected sex-chromosome specific overexpression in stages of spermatogenesis where inactivated X chromosome expression was expected, but the hybrid overexpression phenotypes were fundamentally different. Importantly, misexpression phenotypes support the presence of multiple histological blocks to spermatogenesis in dwarf hamster hybrids, including a potential role of meiotic stalling early in spermatogenesis. Collectively, we demonstrate that while there are some similarities in hybrid regulatory phenotypes of house mice and dwarf hamsters, there are also clear differences that point towards unique mechanisms underlying hybrid male sterility in each system. Our results highlight the potential of comparative approaches in helping to understand the importance of disrupted gene regulation in speciation.

Into the Wild: A novel wild-derived inbred strain resource expands the genomic and phenotypic diversity of laboratory mouse models

The laboratory mouse has served as the premier animal model system for both basic and preclinical investigations for a century. However, laboratory mice capture a narrow subset of the genetic variation found in wild mouse populations. This consideration inherently restricts the scope of potential discovery in laboratory models and narrows the pool of potentially identified phenotype-associated variants and pathways. Wild mouse populations are reservoirs of predicted functional and disease-associated alleles, but the sparsity of commercially available, well-characterized wild mouse strains limits their broader adoption in biomedical research. To overcome this barrier, we have recently imported, sequenced, and phenotyped a set of 11 wild-derived inbred strains developed from wild-caught Mus musculus domesticus. Each of these "Nachman strains" immortalizes a unique wild haplotype sampled from five environmentally diverse locations across North and South America: Saratoga Springs, New York, USA; Gainesville, Florida, USA; Manaus, Brazil; Tucson, Arizona, USA; and Edmonton, Alberta, Canada. Whole genome sequence analysis reveals that each strain carries between 4.73-6.54 million single nucleotide differences relative to the mouse reference assembly, with 42.5% of variants in the Nachman strain genomes absent from classical inbred mouse strains. We phenotyped the Nachman strains on a customized pipeline to assess the scope of disease-relevant neurobehavioral, biochemical, physiological, metabolic, and morphological trait variation. The Nachman strains exhibit significant inter-strain variation in >90% of 1119 surveyed traits and expand the range of phenotypic diversity captured in classical inbred strain panels alone. Taken together, our work introduces a novel wild-derived inbred mouse strain resource that will enable new discoveries in basic and preclinical research. These strains are currently available through The Jackson Laboratory Repository under laboratory code NachJ.

Hybrid breakdown in male reproduction between recently diverged Drosophila melanogaster populations has a complex and variable genetic architecture

Secondary contact between formerly isolated populations may result in hybrid breakdown, in which untested allelic combinations in hybrids are maladaptive and limit genetic exchange. Studying early-stage reproductive isolation may yield key insights into the genetic architectures and evolutionary forces underlying the first steps toward speciation. Here, we leverage the recent worldwide expansion of Drosophila melanogaster to test for hybrid breakdown between populations that diverged within the last 13,000 years. We found clear evidence for hybrid breakdown in male reproduction, but not female reproduction or viability, supporting the prediction that hybrid breakdown affects the heterogametic sex first. The frequency of non-reproducing F2 males varied among different crosses involving the same southern African and European populations, as did the qualitative effect of cross direction, implying a genetically variable basis of hybrid breakdown and a role for uniparentally inherited factors. The levels of breakdown observed in F2 males were not recapitulated in backcrossed individuals, consistent with the existence of incompatibilities with at least three partners. Thus, some of the very first steps toward reproductive isolation could involve incompatibilities with complex and variable genetic architectures. Collectively, our findings emphasize this system's potential for future studies on the genetic and organismal basis of early-stage reproductive isolation.

Genic and Chromosomal Components of Prdm9-Driven Hybrid Male Sterility in Mice (Mus musculus)

Hybrid sterility contributes to speciation by preventing gene flow between related taxa. Prdm9, the first and only hybrid male sterility (HMS) gene known in vertebrates, predetermines the sites of recombination between homologous chromosomes and their synapsis in early meiotic prophase. The asymmetric binding of PRDM9 to heterosubspecific homologs of Mus m. musculus x Mus m. domesticus F1 hybrids and increase of PRDM9-independent DNA double-strand break (DSB) hotspots results in difficult to repair DSBs, incomplete synapsis of homologous chromosomes and meiotic arrest at the first meiotic prophase. Here we show that Prdm9 behaves as a major HMS gene in mice outside the Mus m. musculus x Mus m. domesticus F1 hybrids, in the genomes composed of Mus m. castaneus and Mus m. musculus chromosomes segregating on the Mus m. domesticus background. The Prdm9cst/dom2 (castaneus/domesticus) allelic combination secures meiotic synapsis, testes weight and sperm count within physiological limits, while the Prdm9msc1/dom2 (musculus/domesticus) males show a range of fertility impairment. Out of five quantitative trait loci contributing to the Prdm9msc1/dom2-related infertility, four control either meiotic synapsis or fertility phenotypes and one controls both, synapsis and fertility. Whole-genome genotyping of individual chromosomes showed preferential involvement of nonrecombinant musculus chromosomes in asynapsis in accordance with the chromosomal character of HMS. Moreover, we show that the overall asynapsis rate can be estimated solely from the genotype of individual males by scoring the effect of nonrecombinant musculus chromosomes. Prdm9-controlled HMS represents an example of genetic architecture of HMS consisting of genic and chromosomal components.

Rapid Divergence of Key Spermatogenesis Genes in nasuta-Subgroup of Drosophila

The crosses between closely related Drosophila species usually produce sterile hybrid males with spermatogenesis disrupted at post-meiotic phase, especially in sperm individualization stage than the pre-meiotic stage. This is possibly due to the rapid interspecies divergence of male sex and reproduction-related genes. Here we annotated 11 key spermatogenesis genes in 35 strains of species belonging to nasuta-subgroup of Drosophila, where many interspecies crosses produce sterile males. We characterized the divergence and polymorphism in the protein coding regions by employing gene-wide, codon-wide, and lineage-specific selection analysis to test the mode and strength of selection acting on these genes. Our analysis showed signature of positive selection at bag of marbles (bam) and benign gonial cell neoplasma (bgcn) despite the selection constrains and the absence of endosymbiont infection which could potentially drive rapid divergence due to an arms race while roughex (rux) showed lineage-specific rapid divergence in frontal sheen complex of nasuta-subgroup. cookie monster (comr) showed rapid divergence consistent with the possibility of meiotic arrest observed in sterile hybrids of Drosophila species. Rapid divergence observed at don juan (dj) and Mst98Ca-like was consistent with fused sperm-tail abnormality observed in the hybrids of Drosophila nasuta and Drosophila albomicans. These findings highlight the potential role of rapid nucleotide divergence in bringing about hybrid incompatibility in the form of male sterility; however, additional genetic manipulation studies can widen our understanding of hybrid incompatibilities. Furthermore, our study emphasizes the importance of young species belonging to nasuta-subgroup of Drosophila in studying post-zygotic reproductive isolation mechanisms.

Y chromosomal noncoding RNAs regulate autosomal gene expression via piRNAs in mouse testis

Background

Deciphering the functions of Y chromosome in mammals has been slow owing to the presence of repeats. Some of these repeats transcribe coding RNAs, the roles of which have been studied. Functions of the noncoding transcripts from Y chromosomal repeats however, remain unclear. While a majority of the genes expressed during spermatogenesis are autosomal, mice with different deletions of the long arm of the Y chromosome (Yq) were previously also shown to be characterized by subfertility, sterility and sperm abnormalities, suggesting the presence of effectors of spermatogenesis at this location. Here we report a set of novel noncoding RNAs from mouse Yq and explore their connection to some of the autosomal genes expressed in testis.

Results

We describe a set of novel mouse male-specific Y long arm (MSYq)-derived long noncoding (lnc) transcripts, named Pirmy and Pirmy-like RNAs. Pirmy shows a large number of splice variants in testis. We also identified Pirmy-like RNAs present in multiple copies at different loci on mouse Y chromosome. Further, we identified eight differentially expressed autosome-encoded sperm proteins in a mutant mouse strain, XY^RIIIqdel (2/3 Yq-deleted). Pirmy and Pirmy-like RNAs have homology to 5′/3′UTRs of these deregulated autosomal genes. Several lines of experiments show that these short homologous stretches correspond to piRNAs. Thus, Pirmy and Pirmy-like RNAs act as templates for several piRNAs. In vitro functional assays reveal putative roles for these piRNAs in regulating autosomal genes.

Conclusions

Our study elucidates a set of autosomal genes that are potentially regulated by MSYq-derived piRNAs in mouse testis. Sperm phenotypes from the Yq-deleted mice seem to be similar to that reported in inter-specific male-sterile hybrids. Taken together, this study provides novel insights into possible role of MSYq-derived ncRNAs in male sterility and speciation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01125-x.

Altering the binding properties of PRDM9 partially restores fertility across the species boundary

Sterility or subfertility of male hybrid offspring is commonly observed. This phenomenon contributes to reproductive barriers between the parental populations, an early step in the process of speciation. One frequent cause of such infertility is a failure of proper chromosome pairing during male meiosis. In subspecies of the house mouse, the likelihood of successful chromosome synapsis is improved by the binding of the histone methyltransferase PRDM9 to both chromosome homologues at matching positions. Using genetic manipulation, we altered PRDM9 binding to occur more often at matched sites, and find that chromosome pairing defects can be rescued, not only in an inter-subspecific cross, but also between distinct species. Using different engineered variants, we demonstrate a quantitative link between the degree of matched homologue binding, chromosome synapsis and rescue of fertility in hybrids between Mus musculus and Mus spretus. The resulting partial restoration of fertility reveals additional mechanisms at play that act to lock-in the reproductive isolation between these two species.

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.

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.

Characterization of black spot resistance in diploid roses with QTL detection, meta-analysis and candidate-gene identification

Two environmentally stable QTLs linked to black spot disease resistance in the Rosa wichurana genetic background were detected, in different connected populations, on linkage groups 3 and 5. Co-localization between R-genes and defense response genes was revealed via meta-analysis. The widespread rose black spot disease (BSD) caused by the hemibiotrophic fungus Diplocarpon rosae Wolf. is efficiently controlled with fungicides. However, in the actual context of reducing agrochemical use, the demand for rose bushes with higher levels of resistance has increased. Qualitative resistance conferred by major genes (Rdr genes) has been widely studied but quantitative resistance to BSD requires further investigation. In this study, segregating populations connected through the BSD resistant Rosa wichurana male parent were phenotyped for disease resistance over several years and locations. A pseudo-testcross approach was used, resulting in six parental maps across three populations. A total of 45 individual QTLs with significant effect on BSD resistance were mapped on the male maps (on linkage groups (LG) B3, B4, B5 and B6), and 12 on the female maps (on LG A1, A2, A3, A4 and A5). Two major regions linked to BSD resistance were identified on LG B3 and B5 of the male maps and were integrated into a consensus map built from all three of the male maps. A meta-analysis was used to narrow down the confidence intervals of individual QTLs from three populations by generating meta-QTLs. Two 'hot spots' or meta-QTLs were found per LG, enabling reduction of the confidence interval to 10.42 cM for B3 and 11.47 cM for B5. An expert annotation of NBS-LRR encoding genes of the genome assembly of Hibrand et al. was performed and used to explore potential co-localization with R-genes. Co-localization with defense response genes was also investigated.

Sexual dimorphism in the meiotic requirement for PRDM9: A mammalian evolutionary safeguard

In many mammals, genomic sites for recombination are determined by the histone methyltransferase PRMD9. Some mouse strains lacking PRDM9 are infertile, but instances of fertility or semifertility in the absence of PRDM9 have been reported in mice, canines, and a human female. Such findings raise the question of how the loss of PRDM9 is circumvented to maintain fertility. We show that genetic background and sex-specific modifiers can obviate the requirement for PRDM9 in mice. Specifically, the meiotic DNA damage checkpoint protein CHK2 acts as a modifier allowing female-specific fertility in the absence of PRDM9. We also report that, in the absence of PRDM9, a PRDM9-independent recombination system is compatible with female meiosis and fertility, suggesting sex-specific regulation of meiotic recombination, a finding with implications for speciation.

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.

A high-throughput method for unbiased quantitation and categorization of nuclear morphology†

The physical arrangement of chromatin in the nucleus is cell type and species-specific, a fact particularly evident in sperm, in which most of the cytoplasm has been lost. Analysis of the characteristic falciform (“hook shaped”) sperm in mice is important in studies of sperm development, hybrid sterility, infertility, and toxicology. However, quantification of sperm shape differences typically relies on subjective manual assessment, rendering comparisons within and between samples difficult.

We have developed an analysis program for morphometric analysis of asymmetric nuclei and characterized the sperm of mice from a range of inbred, outbred, and wild-derived mouse strains. We find that laboratory strains have elevated sperm shape variability both within and between samples in comparison to wild-derived inbred strains, and that sperm shape in F1 offspring from a cross between CBA and C57Bl6J strains is subtly affected by the direction of the cross. We further show that hierarchical clustering can discriminate distinct sperm shapes with greater efficiency and reproducibility than even experienced manual assessors, and is useful both to distinguish between samples and also to identify different morphological classes within a single sample.

Our approach allows for the analysis of nuclear shape with unprecedented precision and scale and will be widely applicable to different species and different areas of biology.

Instability of the Pseudoautosomal Boundary in House Mice

Faithful segregation of homologous chromosomes at meiosis requires pairing and recombination. In taxa with dimorphic sex chromosomes, pairing between them in the heterogametic sex is limited to a narrow interval of residual sequence homology known as the pseudoautosomal region (PAR). Failure to form the obligate crossover in the PAR is associated with male infertility in house mice (Mus musculus) and humans. Yet despite this apparent functional constraint, the boundary and organization of the PAR is highly variable in mammals, and even between subspecies of mice. Here, we estimate the genetic map in a previously documented expansion of the PAR in the M. musculus castaneus subspecies and show that the local recombination rate is 100-fold higher than the autosomal background. We identify an independent shift in the PAR boundary in the M. musculus musculus subspecies and show that it involves a complex rearrangement, but still recombines in heterozygous males. Finally, we demonstrate pervasive copy-number variation at the PAR boundary in wild populations of M. m. domesticus, M. m. musculus, and M. m. castaneus Our results suggest that the intensity of recombination activity in the PAR, coupled with relatively weak constraints on its sequence, permit the generation and maintenance of unusual levels of polymorphism in the population of unknown functional significance.

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 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.

Relationship Between Sequence Homology, Genome Architecture, and Meiotic Behavior of the Sex Chromosomes in North American Voles

In most mammals, the X and Y chromosomes synapse and recombine along a conserved region of homology known as the pseudoautosomal region (PAR). These homology-driven interactions are required for meiotic progression and are essential for male fertility. Although the PAR fulfills key meiotic functions in most mammals, several exceptional species lack PAR-mediated sex chromosome associations at meiosis. Here, we leveraged the natural variation in meiotic sex chromosome programs present in North American voles (Microtus) to investigate the relationship between meiotic sex chromosome dynamics and X/Y sequence homology. To this end, we developed a novel, reference-blind computational method to analyze sparse sequencing data from flow-sorted X and Y chromosomes isolated from vole species with sex chromosomes that always (Microtus montanus), never (Microtus mogollonensis), and occasionally synapse (Microtus ochrogaster) at meiosis. Unexpectedly, we find more shared X/Y homology in the two vole species with no and sporadic X/Y synapsis compared to the species with obligate synapsis. Sex chromosome homology in the asynaptic and occasionally synaptic species is interspersed along chromosomes and largely restricted to low-complexity sequences, including a striking enrichment for the telomeric repeat sequence, TTAGGG. In contrast, homology is concentrated in high complexity, and presumably euchromatic, sequence on the X and Y chromosomes of the synaptic vole species, M. montanus Taken together, our findings suggest key conditions required to sustain the standard program of X/Y synapsis at meiosis and reveal an intriguing connection between heterochromatic repeat architecture and noncanonical, asynaptic mechanisms of sex chromosome segregation in voles.

Spermatogenesis and the Evolution of Mammalian Sex Chromosomes

Developmental constraint and sexual conflict shape the evolution of heteromorphic sex chromosomes. These contrasting forces are perhaps strongest during spermatogenesis in species with XY males. In this review, we consider how the unique regulatory environment and selective pressures of spermatogenesis interact to impact sex chromosome evolution in mammals. We explore how each developmental phase of spermatogenesis influences sex chromosome gene content, structure, and rate of molecular evolution, and how these attributes may contribute to speciation. We argue that a developmental context is fundamental to understanding sex chromosome evolution and that an evolutionary perspective can shed new light on our understanding of sperm development.

Genome-wide association for testis weight in the diversity outbred mouse population

Testis weight is a genetically mediated trait associated with reproductive efficiency across numerous species. We sought to evaluate the genetically diverse, highly recombinant Diversity Outbred (DO) mouse population as a tool to identify and map quantitative trait loci (QTLs) associated with testis weight. Testis weights were recorded for 502 male DO mice and the mice were genotyped on the GIGAMuga array at ~ 143,000 SNPs. We performed a genome-wide association analysis and identified one significant and two suggestive QTLs associated with testis weight. Using bioinformatic approaches, we developed a list of candidate genes and identified those with known roles in testicular size and development. Candidates of particular interest include the RNA demethylase gene Alkbh5, the cyclin-dependent kinase inhibitor gene Cdkn2c, the dynein axonemal heavy chain gene Dnah11, the phospholipase D gene Pld6, the trans-acting transcription factor gene Sp4, and the spermatogenesis-associated gene Spata6, each of which has a human ortholog. Our results demonstrate the utility of DO mice in high-resolution genetic mapping of complex traits, enabling us to identify developmentally important genes in adult mice. Understanding how genetic variation in these genes influence testis weight could aid in the understanding of mechanisms of mammalian reproductive function.

Genetic mapping of species differences via in vitro crosses in mouse embryonic stem cells

Discovering the genetic changes underlying species differences is a central goal in evolutionary genetics. However, hybrid crosses between species in mammals often suffer from hybrid sterility, greatly complicating genetic mapping of trait variation across species. Here, we describe a simple, robust, and transgene-free technique to generate "in vitro crosses" in hybrid mouse embryonic stem (ES) cells by inducing random mitotic cross-overs with the drug ML216, which inhibits the DNA helicase Bloom syndrome (BLM). Starting with an interspecific F1 hybrid ES cell line between the Mus musculus laboratory mouse and Mus spretus (∼1.5 million years of divergence), we mapped the genetic basis of drug resistance to the antimetabolite tioguanine to a single region containing hypoxanthine-guanine phosphoribosyltransferase (Hprt) in as few as 21 d through "flow mapping" by coupling in vitro crosses with fluorescence-activated cell sorting (FACS). We also show how our platform can enable direct study of developmental variation by rederiving embryos with contribution from the recombinant ES cell lines. We demonstrate how in vitro crosses can overcome major bottlenecks in mouse complex trait genetics and address fundamental questions in evolutionary biology that are otherwise intractable through traditional breeding due to high cost, small litter sizes, and/or hybrid sterility. In doing so, we describe an experimental platform toward studying evolutionary systems biology in mouse and potentially in human and other mammals, including cross-species hybrids.

Craniomandibular form and body size variation of first generation mouse hybrids: A model for hominin hybridization

Hybridization occurs in a number of mammalian lineages, including among primate taxa. Analyses of ancient genomes have shown that hybridization between our lineage and other archaic hominins in Eurasia occurred numerous times in the past. However, we still have limited empirical data on what a hybrid skeleton looks like, or how to spot patterns of hybridization among fossils for which there are no genetic data. Here we use experimental mouse models to supplement previous studies of primates. We characterize size and shape variation in the cranium and mandible of three wild-derived inbred mouse strains and their first generation (F1) hybrids. The three parent taxa in our analysis represent lineages that diverged over approximately the same period as the human/Neanderthal/Denisovan lineages and their hybrids are variably successful in the wild. Comparisons of body size, as quantified by long bone measurements, are also presented to determine whether the identified phenotypic effects of hybridization are localized to the cranium or represent overall body size changes. The results indicate that hybrid cranial and mandibular sizes, as well as limb length, exceed that of the parent taxa in all cases. All three F1 hybrid crosses display similar patterns of size and form variation. These results are generally consistent with earlier studies on primates and other mammals, suggesting that the effects of hybridization may be similar across very different scenarios of hybridization, including different levels of hybrid fitness. This paper serves to supplement previous studies aimed at identifying F1 hybrids in the fossil record and to introduce further research that will explore hybrid morphologies using mice as a proxy for better understanding hybridization in the hominin fossil record.

The genomic ancestry, landscape genetics and invasion history of introduced mice in New Zealand

The house mouse (Mus musculus) provides a fascinating system for studying both the genomic basis of reproductive isolation, and the patterns of human-mediated dispersal. New Zealand has a complex history of mouse invasions, and the living descendants of these invaders have genetic ancestry from all three subspecies, although most are primarily descended from M. m. domesticus. We used the GigaMUGA genotyping array (approximately 135 000 loci) to describe the genomic ancestry of 161 mice, sampled from 34 locations from across New Zealand (and one Australian city-Sydney). Of these, two populations, one in the south of the South Island, and one on Chatham Island, showed complete mitochondrial lineage capture, featuring two different lineages of M. m. castaneus mitochondrial DNA but with only M. m. domesticus nuclear ancestry detectable. Mice in the northern and southern parts of the North Island had small traces (approx. 2-3%) of M. m. castaneus nuclear ancestry, and mice in the upper South Island had approximately 7-8% M. m. musculus nuclear ancestry including some Y-chromosomal ancestry-though no detectable M. m. musculus mitochondrial ancestry. This is the most thorough genomic study of introduced populations of house mice yet conducted, and will have relevance to studies of the isolation mechanisms separating subspecies of mice.

Three representative inter and intra-subspecific crosses reveal the genetic architecture of reproductive isolation in rice

Systematic characterization of genetic and molecular mechanisms in the formation of hybrid sterility is of fundamental importance in understanding reproductive isolation and speciation. Using ultra-high-density genetic maps, 43 single-locus quantitative trait loci (QTLs) and 223 digenic interactions for embryo-sac, pollen, and spikelet fertility are depicted from three crosses between representative varieties of japonica and two varietal groups of indica, which provide an extensive archive for investigating the genetic basis of reproductive isolation in rice. Ten newly detected single-locus QTLs for inter- and intra-subspecific fertility are identified. Three loci for embryo-sac fertility are detected in both Nip × ZS97 and Nip × MH63 crosses, whereas QTLs for pollen fertility are not in common between the two crosses thus leading to fertility variation. Five loci responsible for fertility and segregation distortion are observed in the ZS97 × MH63 cross. The importance of two-locus interactions on fertility are quantified in the whole genome, which identify that three types of interaction contribute to fertility reduction in the hybrid. These results construct the genetic architecture with respect to various forms of reproductive barriers in rice, which have significant implications in utilization of inter-subspecific heterosis along with improvement in the fertility of indica-indica hybrids at single- and multi-locus level.

Bounds to parapatric speciation: A Dobzhansky-Muller incompatibility model involving autosomes, X chromosomes, and mitochondria

We investigate the conditions for the origin and maintenance of postzygotic isolation barriers, so called (Bateson-)Dobzhansky-Muller incompatibilities or DMIs, among populations that are connected by gene flow. Specifically, we compare the relative stability of pairwise DMIs among autosomes, X chromosomes, and mitochondrial genes. In an analytical approach based on a continent-island framework, we determine how the maximum permissible migration rates depend on the genomic architecture of the DMI, on sex bias in migration rates, and on sex-dependence of allelic and epistatic effects, such as dosage compensation. Our results show that X-linkage of DMIs can enlarge the migration bounds relative to autosomal DMIs or autosome-mitochondrial DMIs, in particular in the presence of dosage compensation. The effect is further strengthened with male-biased migration. This mechanism might contribute to a higher density of DMIs on the X chromosome (large X-effect) that has been observed in several species clades. Furthermore, our results agree with empirical findings of higher introgression rates of autosomal compared to X-linked loci.

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.

Meiotic Consequences of Genetic Divergence Across the Murine Pseudoautosomal Region

The production of haploid gametes during meiosis is dependent on the homology-driven processes of pairing, synapsis, and recombination. On the mammalian heterogametic sex chromosomes, these key meiotic activities are confined to the pseudoautosomal region (PAR), a short region of near-perfect sequence homology between the X and Y chromosomes. Despite its established importance for meiosis, the PAR is rapidly evolving, raising the question of how proper X/Y segregation is buffered against the accumulation of homology-disrupting mutations. Here, I investigate the interplay of PAR evolution and function in two interfertile house mouse subspecies characterized by structurally divergent PARs, Mus musculus domesticus and M. m. castaneus. Using cytogenetic methods to visualize the sex chromosomes at meiosis, I show that intersubspecific F₁ hybrids harbor an increased frequency of pachytene spermatocytes with unsynapsed sex chromosomes. This high rate of asynapsis is due, in part, to the premature release of synaptic associations prior to completion of prophase I. Further, I show that when sex chromosomes do synapse in intersubspecific hybrids, recombination is reduced across the paired region. Together, these meiotic defects afflict ∼50% of spermatocytes from F₁ hybrids and lead to increased apoptosis in meiotically dividing cells. Despite flagrant disruption of the meiotic program, a subset of spermatocytes complete meiosis and intersubspecific F₁ males remain fertile. These findings cast light on the meiotic constraints that shape sex chromosome evolution and offer initial clues to resolve the paradox raised by the rapid evolution of this functionally significant locus.

Hybrid Sterility Locus on Chromosome X Controls Meiotic Recombination Rate in Mouse

Meiotic recombination safeguards proper segregation of homologous chromosomes into gametes, affects genetic variation within species, and contributes to meiotic chromosome recognition, pairing and synapsis. The Prdm9 gene has a dual role, it controls meiotic recombination by determining the genomic position of crossover hotspots and, in infertile hybrids of house mouse subspecies Mus m. musculus (Mmm) and Mus m. domesticus (Mmd), it further functions as the major hybrid sterility gene. In the latter role Prdm9 interacts with the hybrid sterility X 2 (Hstx2) genomic locus on Chromosome X (Chr X) by a still unknown mechanism. Here we investigated the meiotic recombination rate at the genome-wide level and its possible relation to hybrid sterility. Using immunofluorescence microscopy we quantified the foci of MLH1 DNA mismatch repair protein, the cytological counterparts of reciprocal crossovers, in a panel of inter-subspecific chromosome substitution strains. Two autosomes, Chr 7 and Chr 11, significantly modified the meiotic recombination rate, yet the strongest modifier, designated meiotic recombination 1, Meir1, emerged in the 4.7 Mb Hstx2 genomic locus on Chr X. The male-limited transgressive effect of Meir1 on recombination rate parallels the male-limited transgressive role of Hstx2 in hybrid male sterility. Thus, both genetic factors, the Prdm9 gene and the Hstx2/Meir1 genomic locus, indicate a link between meiotic recombination and hybrid sterility. A strong female-specific modifier of meiotic recombination rate with the effect opposite to Meir1 was localized on Chr X, distally to Meir1. Mapping Meir1 to a narrow candidate interval on Chr X is an important first step towards positional cloning of the respective gene(s) responsible for variation in the global recombination rate between closely related mouse subspecies.

During differentiation of germ cells into gametes, a maternal and a paternal copy of each chromosome have to find each other, pair, and synapse in order to ensure proper chromosome segregation into the gametes. Because of the unique ability to identify homologous DNA sequences between homologous chromosomes, meiotic recombination is an essential step in proper chromosome pairing and synapsis in the majority of species. However, when the paternal and maternal sets of chromosomes come from different (sub)species, the recognition of homologs can be disturbed and result in sterility of male hybrids. In this study we investigated the genetic control of variation in the global recombination rate between two closely related mouse subspecies with regard to the known infertility of their F1 hybrids. We show that the variation in the global recombination rate between both subspecies is under the control of three genomic loci. The strongest one appeared within the hybrid sterility X2 genomic locus on Chromosome X. Our findings will allow positional cloning of the gene and will shed new light on the role of meiotic recombination in reproductive isolation between closely related species.

The Founder Strains of the Collaborative Cross Express a Complex Combination of Advantageous and Deleterious Traits for Male Reproduction

Surveys of inbred strains of mice are standard approaches to determine the heritability and range of phenotypic variation for biomedical traits. In addition, they may lead to the identification of novel phenotypes and models of human disease. Surprisingly, male reproductive phenotypes are among the least-represented traits in the Mouse Phenome Database. Here we report the results of a broad survey of the eight founder inbred strains of both the Collaborative Cross (CC) and the Diversity Outbred populations, two new mouse resources that are being used as platforms for systems genetics and sources of mouse models of human diseases. Our survey includes representatives of the three main subspecies of the house mice and a mix of classical and wild-derived inbred strains. In addition to standard staples of male reproductive phenotyping such as reproductive organ weights, sperm counts, and sperm morphology, our survey includes sperm motility and the first detailed survey of testis histology. As expected for such a broad survey, heritability varies widely among traits. We conclude that although all eight inbred strains are fertile, most display a mix of advantageous and deleterious male reproductive traits. The CAST/EiJ strain is an outlier, with an unusual combination of deleterious male reproductive traits including low sperm counts, high levels of morphologically abnormal sperm, and poor motility. In contrast, sperm from the PWK/PhJ and WSB/EiJ strains had the greatest percentages of normal morphology and vigorous motility. Finally, we report an abnormal testis phenotype that is highly heritable and restricted to the WSB/EiJ strain. This phenotype is characterized by the presence of a large, but variable, number of vacuoles in at least 10% of the seminiferous tubules. The onset of the phenotype between 2 and 3 wk of age is temporally correlated with the formation of the blood-testis barrier. We speculate that this phenotype may play a role in high rates of extinction in the CC project and in the phenotypes associated with speciation in genetic crosses that use the WSB/EiJ strain as representative of the Mus muculus domesticus subspecies.

Mechanisms Underlying Mammalian Hybrid Sterility in Two Feline Interspecies Models

The phenomenon of male sterility in interspecies hybrids has been observed for over a century, however, few genes influencing this recurrent phenotype have been identified. Genetic investigations have been primarily limited to a small number of model organisms, thus limiting our understanding of the underlying molecular basis of this well-documented "rule of speciation." We utilized two interspecies hybrid cat breeds in a genome-wide association study employing the Illumina 63 K single-nucleotide polymorphism array. Collectively, we identified eight autosomal genes/gene regions underlying associations with hybrid male sterility (HMS) involved in the function of the blood-testis barrier, gamete structural development, and transcriptional regulation. We also identified several candidate hybrid sterility regions on the X chromosome, with most residing in close proximity to complex duplicated regions. Differential gene expression analyses revealed significant chromosome-wide upregulation of X chromosome transcripts in testes of sterile hybrids, which were enriched for genes involved in chromatin regulation of gene expression. Our expression results parallel those reported in Mus hybrids, supporting the "Large X-Effect" in mammalian HMS and the potential epigenetic basis for this phenomenon. These results support the value of the interspecies feline model as a powerful tool for comparison to rodent models of HMS, demonstrating unique aspects and potential commonalities that underpin mammalian reproductive isolation.

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.

Direct Gamete Sequencing Reveals No Evidence for Segregation Distortion in House Mouse Hybrids

Understanding the molecular basis of species formation is an important goal in evolutionary genetics, and Dobzhansky-Muller incompatibilities are thought to be a common source of postzygotic reproductive isolation between closely related lineages. However, the evolutionary forces that lead to the accumulation of such incompatibilities between diverging taxa are poorly understood. Segregation distorters are believed to be an important source of Dobzhansky-Muller incompatibilities between hybridizing species of Drosophila as well as hybridizing crop plants, but it remains unclear if these selfish genetic elements contribute to reproductive isolation in other taxa. Here, we collected viable sperm from first-generation hybrid male progeny of Mus musculus castaneus and M. m. domesticus, two subspecies of rodent in the earliest stages of speciation. We then genotyped millions of single nucleotide polymorphisms in these gamete pools and tested for a skew in the frequency of parental alleles across the genome. We show that segregation distorters are not measurable contributors to observed infertility in these hybrid males, despite sufficient statistical power to detect even weak segregation distortion with our novel method. Thus, reduced hybrid male fertility in crosses between these nascent species is attributable to other evolutionary forces.

Insights into mammalian biology from the wild house mouse Mus musculus

The house mouse, Mus musculus, was established in the early 1900s as one of the first genetic model organisms owing to its short generation time, comparatively large litters, ease of husbandry, and visible phenotypic variants. For these reasons and because they are mammals, house mice are well suited to serve as models for human phenotypes and disease. House mice in the wild consist of at least three distinct subspecies and harbor extensive genetic and phenotypic variation both within and between these subspecies. Wild mice have been used to study a wide range of biological processes, including immunity, cancer, male sterility, adaptive evolution, and non-Mendelian inheritance. Despite the extensive variation that exists among wild mice, classical laboratory strains are derived from a limited set of founders and thus contain only a small subset of this variation. Continued efforts to study wild house mice and to create new inbred strains from wild populations have the potential to strengthen house mice as a model system.

DOI:http://dx.doi.org/10.7554/eLife.05959.001

Sex ratio meiotic drive as a plausible evolutionary mechanism for hybrid male sterility

Biological diversity on Earth depends on the multiplication of species or speciation, which is the evolution of reproductive isolation such as hybrid sterility between two new species. An unsolved puzzle is the exact mechanism(s) that causes two genomes to diverge from their common ancestor so that some divergent genes no longer function properly in the hybrids. Here we report genetic analyses of divergent genes controlling male fertility and sex ratio in two very young fruitfly species, Drosophila albomicans and D. nasuta. A majority of the genetic divergence for both traits is mapped to the same regions by quantitative trait loci mappings. With introgressions, six major loci are found to contribute to both traits. This genetic colocalization implicates that genes for hybrid male sterility have evolved primarily for controlling sex ratio. We propose that genetic conflicts over sex ratio may operate as a perpetual dynamo for genome divergence. This particular evolutionary mechanism may largely contribute to the rapid evolution of hybrid male sterility and the disproportionate enrichment of its underlying genes on the X chromosome – two patterns widely observed across animals.

Millions of species live on Earth, thanks to an evolutionary process that splits one species to two or more new species. The formation of new species is benchmarked by the evolution of reproductive isolation (RI) such as hybrid sterility between new species. The fundamental question of how RI evolves, however, remains largely unknown. In a pair of very young fruitfly species, we localized six loci expressing dual functions of hybrid male sterility (HMS) and sex ratio distortion, implicating an evolutionary causal link between these two traits. The rapid evolution of HMS widely observed across animal taxa can be attributed to the rapid evolution of genes controlling sex chromosome segregation. All genes in a genome are not equal. This study suggests that conflicts among various parts of a genome might confer strong evolutionary pressure—a mechanism that has hitherto been regarded as rare and could actually be more ubiquitous than currently appreciated.

Impact of the number of Robertsonian chromosomes on germ cell death in wild male house mice

Previous studies in the house mouse have shown that the presence of Robertsonian (Rb) metacentric chromosomes in heterozygous condition affects the process of spermatogenesis. This detrimental effect mainly depends on the number of metacentrics involved and the complexity of the resulting meiotic figures. In this study, we aimed at elucidating the relationship between the chromosomal composition and spermatogenesis impairment in mice present in an area of chromosomal polymorphism (the so-called Barcelona system BRbS) in which Rb mice are surrounded by all acrocentric animals, no established metacentric races are present and the level of structural heterozygosity is relatively low. Using the terminal deoxinucleotidyl transferase-mediated dUTP nick end labelling (TUNEL) assay, we report higher frequency of apoptotic spermatogenetic cells in mice carrying six pairs of metacentrics at the homozygous state than in those carrying two or three fusions at the heterozygous state. Specifically, we detected a higher frequency of TUNEL-positive (T+) tubules and of T+ cells per tubule cross section and also a lower spermatid/spermatocyte ratio. These results indicate that the number of metacentrics at the homozygous state is more influential in determining apoptotic germ cell death than that of moderate chromosome heterozygosity. The percentage of germ cell death lower than 50 % found in our samples and the geographic distribution of the set of metacentrics within the BRbS indicate that although the spermatogenic alterations detected in this area could act as a partial barrier to gene flow, they are not sufficient to prevent Rb chromosomes from spreading in nature.

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.

Interspecific tests of allelism reveal the evolutionary timing and pattern of accumulation of reproductive isolation mutations

Despite extensive theory, little is known about the empirical accumulation and evolutionary timing of mutations that contribute to speciation. Here we combined QTL (Quantitative Trait Loci) analyses of reproductive isolation, with information on species evolutionary relationships, to reconstruct the order and timing of mutations contributing to reproductive isolation between three plant (Solanum) species. To evaluate whether reproductive isolation QTL that appear to coincide in more than one species pair are homologous, we used cross-specific tests of allelism and found evidence for both homologous and lineage-specific (non-homologous) alleles at these co-localized loci. These data, along with isolation QTL unique to single species pairs, indicate that >85% of isolation-causing mutations arose later in the history of divergence between species. Phylogenetically explicit analyses of these data support non-linear models of accumulation of hybrid incompatibility, although the specific best-fit model differs between seed (pairwise interactions) and pollen (multi-locus interactions) sterility traits. Our findings corroborate theory that predicts an acceleration ('snowballing') in the accumulation of isolation loci as lineages progressively diverge, and suggest different underlying genetic bases for pollen versus seed sterility. Pollen sterility in particular appears to be due to complex genetic interactions, and we show this is consistent with a snowball model where later arising mutations are more likely to be involved in pairwise or multi-locus interactions that specifically involve ancestral alleles, compared to earlier arising mutations.

Genome-wide patterns of differentiation among house mouse subspecies

One approach to understanding the genetic basis of speciation is to scan the genomes of recently diverged taxa to identify highly differentiated regions. The house mouse, Mus musculus, provides a useful system for the study of speciation. Three subspecies (M. m. castaneus, M. m. domesticus, and M. m. musculus) diverged ∼350 KYA, are distributed parapatrically, show varying degrees of reproductive isolation in laboratory crosses, and hybridize in nature. We sequenced the testes transcriptomes of multiple wild-derived inbred lines from each subspecies to identify highly differentiated regions of the genome, to identify genes showing high expression divergence, and to compare patterns of differentiation among subspecies that have different demographic histories and exhibit different levels of reproductive isolation. Using a sliding-window approach, we found many genomic regions with high levels of sequence differentiation in each of the pairwise comparisons among subspecies. In all comparisons, the X chromosome was more highly differentiated than the autosomes. Sequence differentiation and expression divergence were greater in the M. m. domesticus-M. m. musculus comparison than in either pairwise comparison with M. m. castaneus, which is consistent with laboratory crosses that show the greatest reproductive isolation between M. m. domesticus and M. m. musculus. Coalescent simulations suggest that differences in estimates of effective population size can account for many of the observed patterns. However, there was an excess of highly differentiated regions relative to simulated distributions under a wide range of demographic scenarios. Overlap of some highly differentiated regions with previous results from QTL mapping and hybrid zone studies points to promising candidate regions for reproductive isolation.

Genetic dissection of a genomic region with pleiotropic effects on domestication traits in maize reveals multiple linked QTL

The domesticated crop maize and its wild progenitor, teosinte, have been used in numerous experiments to investigate the nature of divergent morphologies. This study examines a poorly understood region on the fifth chromosome of maize associated with a number of traits under selection during domestication, using a quantitative trait locus (QTL) mapping population specific to the fifth chromosome. In contrast with other major domestication loci in maize where large-effect, highly pleiotropic, single genes are responsible for phenotypic effects, our study found the region on chromosome five fractionates into multiple-QTL regions, none with singularly large effects. The smallest 1.5-LOD support interval for a QTL contained 54 genes, one of which was a MADS MIKC(C) transcription factor, a family of proteins implicated in many developmental programs. We also used simulated trait data sets to investigate the power of our mapping population to identify QTL for which there is a single underlying causal gene. This analysis showed that while QTL for traits controlled by single genes can be accurately mapped, our population design can detect no more than ∼4.5 QTL per trait even when there are 100 causal genes. Thus when a trait is controlled by ≥5 genes in the simulated data, the number of detected QTL can represent a simplification of the underlying causative factors. Our results show how a QTL region with effects on several domestication traits may be due to multiple linked QTL of small effect as opposed to a single gene with large and pleiotropic effects.

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.

Anolis sex chromosomes are derived from a single ancestral pair

To explain the frequency and distribution of heteromorphic sex chromosomes in the lizard genus Anolis, we compared the relative roles of sex chromosome conservation versus turnover of sex-determining mechanisms. We used model-based comparative methods to reconstruct karyotype evolution and the presence of heteromorphic sex chromosomes onto a newly generated Anolis phylogeny. We found that heteromorphic sex chromosomes evolved multiple times in the genus. Fluorescent in situ hybridization (FISH) of repetitive DNA showed variable rates of Y chromosome degeneration among Anolis species and identified previously undetected, homomorphic sex chromosomes in two species. We confirmed homology of sex chromosomes in the genus by performing FISH of an X-linked bacterial artificial chromosome (BAC) and quantitative PCR of X-linked genes in multiple Anolis species sampled across the phylogeny. Taken together, these results are consistent with long-term conservation of sex chromosomes in the group. Our results pave the way to address additional questions related to Anolis sex chromosome evolution and describe a conceptual framework that can be used to evaluate the origins and evolution of heteromorphic sex chromosomes in other clades.

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.

Genome-wide dissection of hybrid sterility in Drosophila confirms a polygenic threshold architecture

To date, different studies about the genetic basis of hybrid male sterility (HMS), a postzygotic reproductive barrier thoroughly investigated using Drosophila species, have demonstrated that no single major gene can produce hybrid sterility without the cooperation of several genetic factors. Early work using hybrids between Drosophila koepferae (Dk) and Drosophila buzzatii (Db) was consistent with the idea that HMS requires the cooperation of several genetic factors, supporting a polygenic threshold (PT) model. Here we present a genome-wide mapping strategy to test the PT model, analyzing serially backcrossed fertile and sterile males in which the Dk genome was introgressed into the Db background. We identified 32 Dk-specific markers significantly associated with hybrid sterility. Our results demonstrate 1) a strong correlation between the number of segregated sterility markers and males' degree of sterility, 2) the exchangeability among markers, 3) their tendency to cluster into low-recombining chromosomal regions, and 4) the requirement for a minimum number (threshold) of markers to elicit sterility. Although our findings do not contradict a role for occasional major hybrid-sterility genes, they conform more to the view that HMS primarily evolves by the cumulative action of many interacting genes of minor effect in a complex PT architecture.

Diversity of Prdm9 zinc finger array in wild mice unravels new facets of the evolutionary turnover of this coding minisatellite

In humans and mice, meiotic recombination events cluster into narrow hotspots whose genomic positions are defined by the PRDM9 protein via its DNA binding domain constituted of an array of zinc fingers (ZnFs). High polymorphism and rapid divergence of the Prdm9 gene ZnF domain appear to involve positive selection at DNA-recognition amino-acid positions, but the nature of the underlying evolutionary pressures remains a puzzle. Here we explore the variability of the Prdm9 ZnF array in wild mice, and uncovered a high allelic diversity of both ZnF copy number and identity with the caracterization of 113 alleles. We analyze features of the diversity of ZnF identity which is mostly due to non-synonymous changes at codons -1, 3 and 6 of each ZnF, corresponding to amino-acids involved in DNA binding. Using methods adapted to the minisatellite structure of the ZnF array, we infer a phylogenetic tree of these alleles. We find the sister species Mus spicilegus and M. macedonicus as well as the three house mouse (Mus musculus) subspecies to be polyphyletic. However some sublineages have expanded independently in Mus musculus musculus and M. m. domesticus, the latter further showing phylogeographic substructure. Compared to random genomic regions and non-coding minisatellites, none of these patterns appears exceptional. In silico prediction of DNA binding sites for each allele, overlap of their alignments to the genome and relative coverage of the different families of interspersed repeated elements suggest a large diversity between PRDM9 variants with a potential for highly divergent distributions of recombination events in the genome with little correlation to evolutionary distance. By compiling PRDM9 ZnF protein sequences in Primates, Muridae and Equids, we find different diversity patterns among the three amino-acids most critical for the DNA-recognition function, suggesting different diversification timescales.

Failed sperm development as a reproductive isolating barrier between species

Hybrid male sterility is a common reproductive isolating barrier between species. Yet, little is known about the actual developmental causes of this phenomenon, especially in naturally hybridizing species. We sought to evaluate the developmental causes of hybrid male sterility, using spadefoot toads as our study system. Plains spadefoot toads (Spea bombifrons) and Mexican spadefoot toads (S. multiplicata) hybridize where they co-occur in the southwestern USA. Hybrids are viable, but hybrid males suffer reduced fertility. We compared testes size and developmental stages of sperm cell maturation between hybrid males and males of each species. We found that testes of hybrid males did not differ in mean size from pure-species males. However, hybrids showed a greater range of within-individual variation in testes size than pure-species males. Moreover, although hybrids produced similar numbers of early stage sperm cells, hybrids produced significantly fewer mature spermatozoids than pure-species males. Interestingly, an introgressed individual produced numbers of live sperm comparable to pure-species males, but the majority of these sperm cells were abnormally shaped and non-motile. These results indicate that hybrid incompatibilities in late sperm development serve as a reproductive isolating barrier between species. The nature of this breakdown highlights the possibilities that hybrid males may vary in fertility and that fertility could possibly be recovered in introgressed males.

A novel approach identifying hybrid sterility QTL on the autosomes of Drosophila simulans and D. mauritiana

When species interbreed, the hybrid offspring that are produced are often sterile. If only one hybrid sex is sterile, it is almost always the heterogametic (XY or ZW) sex. Taking this trend into account, the predominant model used to explain the genetic basis of F₁ sterility involves a deleterious interaction between recessive sex-linked loci from one species and dominant autosomal loci from the other species. This model is difficult to evaluate, however, as only a handful of loci influencing interspecies hybrid sterility have been identified, and their autosomal genetic interactors have remained elusive. One hindrance to their identification has been the overwhelming effect of the sex chromosome in mapping studies, which could ‘mask’ the ability to accurately map autosomal factors. Here, we use a novel approach employing attached-X chromosomes to create reciprocal backcross interspecies hybrid males that have a non-recombinant sex chromosome and recombinant autosomes. The heritable variation in phenotype is thus solely caused by differences in the autosomes, thereby allowing us to accurately identify the number and location of autosomal sterility loci. In one direction of backcross, all males were sterile, indicating that sterility could be entirely induced by the sex chromosome complement in these males. In the other direction, we identified nine quantitative trait loci that account for a surprisingly large amount (56%) of the autosome-induced phenotypic variance in sterility, with a large contribution of autosome-autosome epistatic interactions. These loci are capable of acting dominantly, and thus could contribute to F₁ hybrid sterility.