Genetic analysis of genome-scale recombination rate evolution in house mice

Genetic background affects the strength of crossover interference in house mice

Meiotic recombination is required for faithful chromosome segregation in most sexually reproducing organisms and shapes the distribution of genetic variation in populations. Both the overall rate and the spatial distribution of crossovers vary within and between species. Adjacent crossovers on the same chromosome tend to be spaced more evenly than expected at random, a phenomenon known as crossover interference. Although interference has been observed in many taxa, the factors that influence the strength of interference are not well understood. We used house mice ( Mus musculus ), a well-established model system for understanding recombination, to study the effects of genetics and age on recombination rate and interference in the male germline. We analyzed crossover positions in 503 progeny from reciprocal F1 hybrids between inbred strains representing the three major subspecies of house mice. Consistent with previous studies, autosomal alleles from M. m. musculus tend to increase recombination rate, while inheriting a M. m. musculus X chromosome decreases recombination rate. Old males transmit an average of 0.6 more crossovers per meiosis (5.0%) than young males, though the effect varies across genetic backgrounds. We show that the strength of crossover interference depends on genotype, providing a rare demonstration that interference evolves over short timescales. Differences between reciprocal F1s suggest that X-linked factors modulate the strength of interference. Our findings motivate additional comparisons of interference among recently diverged species and further examination of the role of paternal age in determining the number and positioning of crossovers.

Understanding the Genetic Basis of Variation in Meiotic Recombination: Past, Present, and Future

Meiotic recombination is a fundamental feature of sexually reproducing species. It is often required for proper chromosome segregation and plays important role in adaptation and the maintenance of genetic diversity. The molecular mechanisms of recombination are remarkably conserved across eukaryotes, yet meiotic genes and proteins show substantial variation in their sequence and function, even between closely related species. Furthermore, the rate and distribution of recombination shows a huge diversity within and between chromosomes, individuals, sexes, populations, and species. This variation has implications for many molecular and evolutionary processes, yet how and why this diversity has evolved is not well understood. A key step in understanding trait evolution is to determine its genetic basis-that is, the number, effect sizes, and distribution of loci underpinning variation. In this perspective, I discuss past and current knowledge on the genetic basis of variation in recombination rate and distribution, explore its evolutionary implications, and present open questions for future research.

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.

Investigating the evolutionary dynamics of diploid number variation in Ctenomys (Ctenomyidae, Rodentia)

Contrary to predictions from classical hybrid sterility models of chromosomal speciation, some organisms display high rates of karyotype variation. Ctenomys are the current mammals with the greatest interspecific and intraspecific chromosomal variation. A large number of species have been studied cytogenetically. The diploid numbers of chromosomes range from 2n = 10 to 2n = 70. Here, we analyzed karyotype evolution in Ctenomys using comparative phylogenetic methods. We found a strong phylogenetic signal with chromosome number. This refutes the chromosomal megaevolution model, which proposes the independent accumulation of multiple chromosomal rearrangements in each closely related species. We found that Brownian motion (BM) described the observed characteristic changes more thoroughly than the Ornstein-Uhlenbeck and Early-Burst models. This suggests that the evolution of chromosome numbers occurs by a random walk along phylogenetic clades. However, our data indicate that the BM model alone does not fully characterize the chromosomal evolution of Ctenomys.

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.

Rapid Evolution of the Fine-scale Recombination Landscape in Wild House Mouse (Mus musculus) Populations

Meiotic recombination is an important evolutionary force and an essential meiotic process. In many species, recombination events concentrate into hotspots defined by the site-specific binding of PRMD9. Rapid evolution of Prdm9's zinc finger DNA-binding array leads to remarkably abrupt shifts in the genomic distribution of hotspots between species, but the question of how Prdm9 allelic variation shapes the landscape of recombination between populations remains less well understood. Wild house mice (Mus musculus) harbor exceptional Prdm9 diversity, with >150 alleles identified to date, and pose a particularly powerful system for addressing this open question. We employed a coalescent-based approach to construct broad- and fine-scale sex-averaged recombination maps from contemporary patterns of linkage disequilibrium in nine geographically isolated wild house mouse populations, including multiple populations from each of three subspecies. Comparing maps between wild mouse populations and subspecies reveals several themes. First, we report weak fine- and broad-scale recombination map conservation across subspecies and populations, with genetic divergence offering no clear prediction for recombination map divergence. Second, most hotspots are unique to one population, an outcome consistent with minimal sharing of Prdm9 alleles between surveyed populations. Finally, by contrasting aggregate hotspot activity on the X versus autosomes, we uncover evidence for population-specific differences in the degree and direction of sex dimorphism for recombination. Overall, our findings illuminate the variability of both the broad- and fine-scale recombination landscape in M. musculus and underscore the functional impact of Prdm9 allelic variation in wild mouse populations.

Meiotic pairing and double-strand break formation along the heteromorphic threespine stickleback sex chromosomes

Double-strand break repair during meiosis is normally achieved using the homologous chromosome as a repair template. Heteromorphic sex chromosomes share little sequence homology, presenting unique challenges to the repair of double-strand breaks. Our understanding of how heteromorphic sex chromosomes behave during meiosis has been focused on ancient sex chromosomes, where the X and Y differ markedly in overall structure and gene content. It remains unclear how more recently evolved sex chromosomes that share considerably more sequence homology with one another pair and form double-strand breaks. One possibility is barriers to pairing evolve rapidly. Alternatively, recently evolved sex chromosomes may exhibit pairing and double-strand break repair that more closely resembles that of their autosomal ancestors. Here, we use the recently evolved X and Y chromosomes of the threespine stickleback fish (Gasterosteus aculeatus) to study patterns of pairing and double-stranded break formation using molecular cytogenetics. We found that the sex chromosomes of threespine stickleback fish did not pair exclusively in the pseudoautosomal region. Instead, the chromosomes fully paired in a non-homologous fashion. To achieve this, the X chromosome underwent synaptic adjustment during pachytene to match the axis length of the Y chromosome. Double-strand break formation and repair rate also matched that of the autosomes. Our results highlight that recently evolved sex chromosomes exhibit meiotic behavior that is reminiscent of autosomes and argues for further work to identify the homologous templates that are used to repair double-strand breaks on the X and Y chromosomes.

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.

High male specific contribution of the X-chromosome to individual global recombination rate in dairy cattle

Background

Meiotic recombination plays an important role in reproduction and evolution. The individual global recombination rate (GRR), measured as the number of crossovers (CO) per gametes, is a complex trait that has been shown to be heritable. The sex chromosomes play an important role in reproduction and fertility related traits. Therefore, variants present on the X-chromosome might have a high contribution to the genetic variation of GRR that is related to meiosis and to reproduction.

Results

We herein used genotyping data from 58,474 New Zealand dairy cattle to estimate the contribution of the X-chromosome to male and female GRR levels. Based on the pedigree-based relationships, we first estimated that the X-chromosome accounted for 30% of the total additive genetic variance for male GRR. This percentage was equal to 19.9% when the estimation relied on a SNP-BLUP approach assuming each SNP has a small contribution. We then carried out a haplotype-based association study to map X-linked QTL, and subsequently fine-mapped the identified QTL with imputed sequence variants. With this approach we identified three QTL with large effect accounting for 7.7% of the additive genetic variance of male GRR. The associated effects were equal to + 0.79, − 1.16 and + 1.18 CO for the alternate alleles. In females, the estimated contribution of the X-chromosome to GRR was null and no significant association with X-linked loci was found. Interestingly, two of the male GRR QTL were associated with candidate genes preferentially expressed in testis, in agreement with a male-specific effect. Finally, the most significant QTL was associated with PPP4R3C, further supporting the important role of protein phosphatase in double-strand break repair by homologous recombination.

Conclusions

Our study illustrates the important role the X-chromosome can have on traits such as individual recombination rate, associated with testis in males. We also show that contribution of the X-chromosome to such a trait might be sex dependent.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08328-8.

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

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

Sexual dimorphism and the effect of wild introgressions on recombination in cassava (Manihot esculenta Crantz) breeding germplasm

Recombination has essential functions in meiosis, evolution, and breeding. The frequency and distribution of crossovers dictate the generation of new allele combinations and can vary across species and between sexes. Here, we examine recombination landscapes across the 18 chromosomes of cassava (Manihot esculenta Crantz) with respect to male and female meioses and known introgressions from the wild relative Manihot glaziovii. We used SHAPEIT2 and duoHMM to infer crossovers from genotyping-by-sequencing data and a validated multigenerational pedigree from the International Institute of Tropical Agriculture cassava breeding germplasm consisting of 7020 informative meioses. We then constructed new genetic maps and compared them to an existing map previously constructed by the International Cassava Genetic Map Consortium. We observed higher recombination rates in females compared to males, and lower recombination rates in M. glaziovii introgression segments on chromosomes 1 and 4, with suppressed recombination along the entire length of the chromosome in the case of the chromosome 4 introgression. Finally, we discuss hypothesized mechanisms underlying our observations of heterochiasmy and crossover suppression and discuss the broader implications for plant breeding.

Diet effects on mouse meiotic recombination: a warning for recombination studies

Meiotic recombination is a critical process for sexually reproducing organisms. This exchange of genetic information between homologous chromosomes during meiosis is important not only because it generates genetic diversity, but also because it is often required for proper chromosome segregation. Consequently, the frequency and distribution of crossovers are tightly controlled to ensure fertility and offspring viability. However, in many systems, it has been shown that environmental factors can alter the frequency of crossover events. Two studies in flies and yeast point to nutritional status affecting the frequency of crossing over. However, this question remains unexplored in mammals. Here, we test how crossover frequency varies in response to diet in Mus musculus males. We use immunohistochemistry to estimate crossover frequency in multiple genotypes under two diet treatments. Our results indicate that while crossover frequency was unaffected by diet in some strains, other strains were sensitive even to small composition changes between two common laboratory chows. Therefore, recombination is both resistant and sensitive to certain dietary changes in a strain-dependent manner and, hence, this response is genetically determined. Our study is the first to report a nutrition effect on genome-wide levels of recombination. Moreover, our work highlights the importance of controlling diet in recombination studies and may point to diet as a potential source of variability among studies, which is relevant for reproducibility.

Evolution and Plasticity of Genome-Wide Meiotic Recombination Rates

Sex, as well as meiotic recombination between homologous chromosomes, is nearly ubiquitous among eukaryotes. In those species that use it, recombination is important for chromosome segregation during gamete production, and thus for fertility. Strikingly, although in most species only one crossover event per chromosome is required to ensure proper segregation, recombination rates vary considerably above this minimum and show variation within and among species. However, whether this variation in recombination is adaptive or neutral and what might shape it remain unclear. Empirical studies and theory support the idea that recombination is generally beneficial but can also have costs. Here, we review variation in genome-wide recombination rates, explore what might cause this, and discuss what is known about its mechanistic basis. We end by discussing the environmental sensitivity of meiosis and recombination rates, how these features may relate to adaptation, and their implications for a broader understanding of recombination rate evolution. Expected final online publication date for the Annual Review of Genetics, Volume 55 is November 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

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.

Sex Differences in the Recombination Landscape

Sex differences in overall recombination rates are well known, but little theoretical or empirical attention has been given to how and why sexes differ in their recombination landscapes: the patterns of recombination along chromosomes. In the first scientific review of this phenomenon, we find that recombination is biased toward telomeres in males and more uniformly distributed in females in most vertebrates and many other eukaryotes. Notable exceptions to this pattern exist, however. Fine-scale recombination patterns also frequently differ between males and females. The molecular mechanisms responsible for sex differences remain unclear, but chromatin landscapes play a role. Why these sex differences evolve also is unclear. Hypotheses suggest that they may result from sexually antagonistic selection acting on coding genes and their regulatory elements, meiotic drive in females, selection during the haploid phase of the life cycle, selection against aneuploidy, or mechanistic constraints. No single hypothesis, however, can adequately explain the evolution of sex differences in all cases. Sex-specific recombination landscapes have important consequences for population differentiation and sex chromosome evolution.

Molecular evolution of the meiotic recombination pathway in mammals

Meiotic recombination shapes evolution and helps to ensure proper chromosome segregation in most species that reproduce sexually. Recombination itself evolves, with species showing considerable divergence in the rate of crossing-over. However, the genetic basis of this divergence is poorly understood. Recombination events are produced via a complicated, but increasingly well-described, cellular pathway. We apply a phylogenetic comparative approach to a carefully selected panel of genes involved in the processes leading to crossovers-spanning double-strand break formation, strand invasion, the crossover/non-crossover decision, and resolution-to reconstruct the evolution of the recombination pathway in eutherian mammals and identify components of the pathway likely to contribute to divergence between species. Eleven recombination genes, predominantly involved in the stabilization of homologous pairing and the crossover/non-crossover decision, show evidence of rapid evolution and positive selection across mammals. We highlight TEX11 and associated genes involved in the synaptonemal complex and the early stages of the crossover/non-crossover decision as candidates for the evolution of recombination rate. Evolutionary comparisons to MLH1 count, a surrogate for the number of crossovers, reveal a positive correlation between genome-wide recombination rate and the rate of evolution at TEX11 across the mammalian phylogeny. Our results illustrate the power of viewing the evolution of recombination from a pathway perspective.

Parallel Molecular Evolution in Pathways, Genes, and Sites in High-Elevation Hummingbirds Revealed by Comparative Transcriptomics

High-elevation organisms experience shared environmental challenges that include low oxygen availability, cold temperatures, and intense ultraviolet radiation. Consequently, repeated evolution of the same genetic mechanisms may occur across high-elevation taxa. To test this prediction, we investigated the extent to which the same biochemical pathways, genes, or sites were subject to parallel molecular evolution for 12 Andean hummingbird species (family: Trochilidae) representing several independent transitions to high elevation across the phylogeny. Across high-elevation species, we discovered parallel evolution for several pathways and genes with evidence of positive selection. In particular, positively selected genes were frequently part of cellular respiration, metabolism, or cell death pathways. To further examine the role of elevation in our analyses, we compared results for low- and high-elevation species and tested different thresholds for defining elevation categories. In analyses with different elevation thresholds, positively selected genes reflected similar functions and pathways, even though there were almost no specific genes in common. For example, EPAS1 (HIF2α), which has been implicated in high-elevation adaptation in other vertebrates, shows a signature of positive selection when high-elevation is defined broadly (>1,500 m), but not when defined narrowly (>2,500 m). Although a few biochemical pathways and genes change predictably as part of hummingbird adaptation to high-elevation conditions, independent lineages have rarely adapted via the same substitutions.

Inference of recombination maps from a single pair of genomes and its application to ancient samples

Understanding the causes and consequences of recombination landscape evolution is a fundamental goal in genetics that requires recombination maps from across the tree of life. Such maps can be obtained from population genomic datasets, but require large sample sizes. Alternative methods are therefore necessary to research organisms where such datasets cannot be generated easily, such as non-model or ancient species. Here we extend the sequentially Markovian coalescent model to jointly infer demography and the spatial variation in recombination rate. Using extensive simulations and sequence data from humans, fruit-flies and a fungal pathogen, we demonstrate that iSMC accurately infers recombination maps under a wide range of scenarios-remarkably, even from a single pair of unphased genomes. We exploit this possibility and reconstruct the recombination maps of ancient hominins. We report that the ancient and modern maps are correlated in a manner that reflects the established phylogeny of Neanderthals, Denisovans, and modern human populations.

The Past and Future of Experimental Speciation

Speciation is the result of evolutionary processes that generate barriers to gene flow between populations, facilitating reproductive isolation. Speciation is typically studied via theoretical models and snapshot tests in natural populations. Experimental speciation enables real-time direct tests of speciation theory and has been long touted as a critical complement to other approaches. We argue that, despite its promise to elucidate the evolution of reproductive isolation, experimental speciation has been underutilised and lags behind other contributions to speciation research. We review recent experiments and outline a framework for how experimental speciation can be implemented to address current outstanding questions that are otherwise challenging to answer. Greater uptake of this approach is necessary to rapidly advance understanding of speciation.

A first genetic portrait of synaptonemal complex variation

The synaptonemal complex (SC) is a proteinaceous scaffold required for synapsis and recombination between homologous chromosomes during meiosis. Although the SC has been linked to differences in genome-wide crossover rates, the genetic basis of standing variation in SC structure remains unknown. To investigate the possibility that recombination evolves through changes to the SC, we characterized the genetic architecture of SC divergence on two evolutionary timescales. Applying a novel digital image analysis technique to spermatocyte spreads, we measured total SC length in 9,532 spermatocytes from recombinant offspring of wild-derived mouse strains with differences in this fundamental meiotic trait. Using this large dataset, we identified the first known genomic regions involved in the evolution of SC length. Distinct loci affect total SC length divergence between and within subspecies, with the X chromosome contributing to both. Joint genetic analysis of MLH1 foci—immunofluorescent markers of crossovers—from the same spermatocytes revealed that two of the identified loci also confer differences in the genome-wide recombination rate. Causal mediation analysis suggested that one pleiotropic locus acts early in meiosis to designate crossovers prior to SC assembly, whereas a second locus primarily shapes crossover number through its effect on SC length. One genomic interval shapes the relationship between SC length and recombination rate, likely modulating the strength of crossover interference. Our findings pinpoint SC formation as a key step in the evolution of recombination and demonstrate the power of genetic mapping on standing variation in the context of the recombination pathway.

During the first stages of meiosis, the chromosome axes are organized along a protein scaffold in preparation for recombination and their subsequent segregation. This scaffold, known as the synaptonemal complex (SC), is critical for the regular progression of recombination. A complex relationship exists between the organization of the SC, the frequency of recombination, and the likelihood of improper chromosome segregation. In this study, we investigate the genetics of synaptonemal complex variation in the house mouse and connect it with variation in the rate of recombination. We found five loci and several compelling candidate genes responsible for the evolution of synaptonemal complex length within and between mouse subspecies. Several of these loci also affect recombination rate, and our joint analyses of the phenotypes suggest an order by which their effects manifest within the recombination pathway. Our results show that evolution of SC length is crucial to recombination rate divergence. Our work here also demonstrates that genetic analysis of additional meiotic phenotypes can help explain the evolution of recombination, a fundamental evolutionary force.

Conservation of the genome-wide recombination rate in white-footed mice

Despite being linked to the fundamental processes of chromosome segregation and offspring diversification, meiotic recombination rates vary within and between species. Recent years have seen progress in quantifying recombination rate evolution across multiple temporal and genomic scales. Nevertheless, the level of variation in recombination rate within wild populations-a key determinant of evolution in this trait-remains poorly documented on the genomic scale. To address this notable gap, we used immunofluorescent cytology to quantify genome-wide recombination rates in males from a wild population of the white-footed mouse, Peromyscus leucopus. For comparison, we measured recombination rates in a second population of male P. leucopus raised in the laboratory and in male deer mice from the subspecies Peromyscus maniculatus bairdii. Although we found differences between individuals in the genome-wide recombination rate, levels of variation were low-within populations, between populations, and between species. Quantification of synaptonemal complex length and crossover positions along chromosome 1 using a novel automated approach also revealed conservation in broad-scale crossover patterning, including strong crossover interference. We propose stabilizing selection targeting recombination or correlated processes as the explanation for these patterns.

Interbreed variation in meiotic recombination rate and distribution in the domestic chicken Gallus gallus

The efficiency of natural and artificial selection is critically dependent on the recombination rate. However, interbreed and individual variation in recombination rate in poultry remains unknown. Conventional methods of analysis of recombination such as genetic linkage analysis, sperm genotyping and chiasma count at lampbrush chromosomes are expensive and time-consuming. In this study, we analyzed the number and distribution of recombination nodules in spermatocytes of the roosters of six chicken breeds using immunolocalization of key proteins involved in chromosome pairing and recombination. We revealed significant effects of breed (R 2 = 0.17 ; p < 0.001 ) and individual (R 2 = 0.28 ; p < 0.001 ) on variation in the number of recombination nodules. Both interbreed and individual variations in recombination rate were almost entirely determined by variation in recombination density on macrochromosomes, because almost all microchromosomes in each breed had one recombination nodule. Despite interbreed differences in the density of recombination nodules, the patterns of their distribution along homologous chromosomes were similar. The breeds examined in this study showed a correspondence between the age of the breed and its recombination rate. Those with high recombination rates (Pervomai, Russian White and Brahma) are relatively young breeds created by crossing several local breeds. The breeds displaying low recombination rate are ancient local breeds: Cochin (Indo-China), Brown Leghorn (Tuscany, Italy) and Russian Crested (the European part of Russia).

Divergent Fine-Scale Recombination Landscapes between a Freshwater and Marine Population of Threespine Stickleback Fish

Meiotic recombination is a highly conserved process that has profound effects on genome evolution. At a fine-scale, recombination rates can vary drastically across genomes, often localized into small recombination "hotspots" with highly elevated rates, surrounded by regions with little recombination. In most species studied, the location of hotspots within genomes is highly conserved across broad evolutionary timescales. The main exception to this pattern is in mammals, where hotspot location can evolve rapidly among closely related species and even among populations within a species. Hotspot position in mammals is controlled by the gene, Prdm9, whereas in species with conserved hotspots, a functional Prdm9 is typically absent. Due to a limited number of species where recombination rates have been estimated at a fine-scale, it remains unclear whether hotspot conservation is always associated with the absence of a functional Prdm9. Threespine stickleback fish (Gasterosteus aculeatus) are an excellent model to examine the evolution of recombination over short evolutionary timescales. Using a linkage disequilibrium-based approach, we found recombination rates indeed varied at a fine-scale across the genome, with many regions organized into narrow hotspots. Hotspots had highly divergent landscapes between stickleback populations, where only ∼15% of these hotspots were shared. Our results indicate that fine-scale recombination rates may be diverging between closely related populations of threespine stickleback fish. Interestingly, we found only a weak association of a PRDM9 binding motif within hotspots, which suggests that threespine stickleback fish may possess a novel mechanism for targeting recombination hotspots at a fine-scale.

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.

Genetic variation in thyroid folliculogenesis influences susceptibility to hypothyroidism-induced hearing impairment

Maternal and fetal sources of thyroid hormone are important for the development of many organ systems. Thyroid hormone deficiency causes variable intellectual disability and hearing impairment in mouse and man, but the basis for this variation is not clear. To explore this variation, we studied two thyroid hormone-deficient mouse mutants with mutations in pituitary-specific transcription factors, POU1F1 and PROP1, that render them unable to produce thyroid stimulating hormone. DW/J-Pou1f1^dw/dw mice have profound deafness and both neurosensory and conductive hearing impairment, while DF/B-Prop1^df/df mice have modest elevations in hearing thresholds consistent with developmental delay, eventually achieving normal hearing ability. The thyroid glands of Pou1f1 mutants are more severely affected than those of Prop1^df/df mice, and they produce less thyroglobulin during the neonatal period critical for establishing hearing. We previously crossed DW/J-Pou1f1^dw/+ and Cast/Ei mice and mapped a major locus on Chromosome 2 that protects against hypothyroidism-induced hearing impairment in Pou1f1^dw/dw mice: modifier of dw hearing (Mdwh). Here we refine the location of Mdwh by genotyping 196 animals with 876 informative SNPs, and we conduct novel mapping with a DW/J-Pou1f1^dw/+ and 129/P2 cross that reveals 129/P2 mice also have a protective Mdwh locus. Using DNA sequencing of DW/J and DF/B strains, we determined that the genes important for thyroid gland function within Mdwh vary in amino acid sequence between strains that are susceptible or resistant to hypothyroidism-induced hearing impairment. These results suggest that the variable effects of congenital hypothyroidism on the development of hearing ability are attributable to genetic variation in postnatal thyroid gland folliculogenesis and function.

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.

High-resolution crossover mapping reveals similarities and differences of male and female recombination in maize

Meiotic crossovers (COs) are not uniformly distributed across the genome. Factors affecting this phenomenon are not well understood. Although many species exhibit large differences in CO numbers between sexes, sex-specific aspects of CO landscape are particularly poorly elucidated. Here, we conduct high-resolution CO mapping in maize. Our results show that CO numbers as well as their overall distribution are similar in male and female meioses. There are, nevertheless, dissimilarities at local scale. Male and female COs differ in their locations relative to transcription start sites in gene promoters and chromatin marks, including nucleosome occupancy and tri-methylation of lysine 4 of histone H3 (H3K4me3). Our data suggest that sex-specific factors not only affect male–female CO number disparities but also cause fine differences in CO positions. Differences between male and female CO landscapes indicate that recombination has distinct implications for population structure and gene evolution in male and in female meioses.

Sex-specific meiotic crossover (CO) landscapes have been identified in multiple species. Here, the authors show that male and female meioses in maize have similar CO landscapes, and differences between COs in the two sexes only exists in their location relative to transcription start sites and some chromatin marks.

Connecting theory and data to understand recombination rate evolution

Meiotic recombination is necessary for successful gametogenesis in most sexually reproducing organisms and is a fundamental genomic parameter, influencing the efficacy of selection and the fate of new mutations. The molecular and evolutionary functions of recombination should impose strong selective constraints on the range of recombination rates. Yet, variation in recombination rate is observed on a variety of genomic and evolutionary scales. In the past decade, empirical studies have described variation in recombination rate within genomes, between individuals, between sexes, between populations and between species. At the same time, theoretical work has provided an increasingly detailed picture of the evolutionary advantages to recombination. Perhaps surprisingly, the causes of natural variation in recombination rate remain poorly understood. We argue that empirical and theoretical approaches to understand the evolution of recombination have proceeded largely independently of each other. Most models that address the evolution of recombination rate were created to explain the evolutionary advantage of recombination rather than quantitative differences in rate among individuals. Conversely, most empirical studies aim to describe variation in recombination rate, rather than to test evolutionary hypotheses. In this Perspective, we argue that efforts to integrate the rich bodies of empirical and theoretical work on recombination rate are crucial to moving this field forward. We provide new directions for the development of theory and the production of data that will jointly close this gap.This article is part of the themed issue 'Evolutionary causes and consequences of recombination rate variation in sexual organisms'.

Mammalian Comparative Genomics Reveals Genetic and Epigenetic Features Associated with Genome Reshuffling in Rodentia

Understanding how mammalian genomes have been reshuffled through structural changes is fundamental to the dynamics of its composition, evolutionary relationships between species and, in the long run, speciation. In this work, we reveal the evolutionary genomic landscape in Rodentia, the most diverse and speciose mammalian order, by whole-genome comparisons of six rodent species and six representative outgroup mammalian species. The reconstruction of the evolutionary breakpoint regions across rodent phylogeny shows an increased rate of genome reshuffling that is approximately two orders of magnitude greater than in other mammalian species here considered. We identified novel lineage and clade-specific breakpoint regions within Rodentia and analyzed their gene content, recombination rates and their relationship with constitutive lamina genomic associated domains, DNase I hypersensitivity sites and chromatin modifications. We detected an accumulation of protein-coding genes in evolutionary breakpoint regions, especially genes implicated in reproduction and pheromone detection and mating. Moreover, we found an association of the evolutionary breakpoint regions with active chromatin state landscapes, most probably related to gene enrichment. Our results have two important implications for understanding the mechanisms that govern and constrain mammalian genome evolution. The first is that the presence of genes related to species-specific phenotypes in evolutionary breakpoint regions reinforces the adaptive value of genome reshuffling. Second, that chromatin conformation, an aspect that has been often overlooked in comparative genomic studies, might play a role in modeling the genomic distribution of evolutionary breakpoints.

The Recombination Landscape in Wild House Mice Inferred Using Population Genomic Data

Characterizing variation in the rate of recombination across the genome is important for understanding several evolutionary processes. Previous analysis of the recombination landscape in laboratory mice has revealed that the different subspecies have different suites of recombination hotspots. It is unknown, however, whether hotspots identified in laboratory strains reflect the hotspot diversity of natural populations or whether broad-scale variation in the rate of recombination is conserved between subspecies. In this study, we constructed fine-scale recombination rate maps for a natural population of the Eastern house mouse, Mus musculus castaneus We performed simulations to assess the accuracy of recombination rate inference in the presence of phase errors, and we used a novel approach to quantify phase error. The spatial distribution of recombination events is strongly positively correlated between our castaneus map, and a map constructed using inbred lines derived predominantly from M. m. domesticus Recombination hotspots in wild castaneus show little overlap, however, with the locations of double-strand breaks in wild-derived house mouse strains. Finally, we also find that genetic diversity in M. m. castaneus is positively correlated with the rate of recombination, consistent with pervasive natural selection operating in the genome. Our study suggests that recombination rate variation is conserved at broad scales between house mouse subspecies, but it is not strongly conserved at fine scales.

Evolutionary mechanisms of runaway chromosome number change in Agrodiaetus butterflies

Despite predictions of the classic, hybrid-sterility model of chromosomal speciation, some organisms demonstrate high rate of karyotype evolution. This rate is especially impressive in Agrodiaetus butterflies that rapidly evolved the greatest chromosome number diversity known in animal kingdom within a single subgenus. Here we analyzed karyotype evolution in Agrodiaetus using phylogenetic comparative methods. We found that chromosome numbers possess a strong phylogenetic signal. This disproves the chromosome megaevolution model that proposes multiple chromosome rearrangements to accumulate independently in each of closely related species. We found that Brownian motion gives a more adequate description of observed trait changes than Ornstein-Uhlenbeck model. This indicates that chromosome numbers evolve via random walk along branches of the phylogeny. We discovered a correlation between karyotype changes and phylogeny branch lengths. This gradual pattern is inconsistent with the hybrid-sterility model which, due to association of major chromosome changes with cladogenetic events, predicts a high degree of punctualism in karyotype evolution. Thus, low underdominace of chromosomal rearrangements and/or prevalence of the recombination-suppression model over the hybrid-sterility model of chromosome speciation are the most common engines of the runaway chromosome number change observed.

Genetics of Genome-Wide Recombination Rate Evolution in Mice from an Isolated Island

Recombination rate is a heritable quantitative trait that evolves despite the fundamentally conserved role that recombination plays in meiosis. Differences in recombination rate can alter the landscape of the genome and the genetic diversity of populations. Yet our understanding of the genetic basis of recombination rate evolution in nature remains limited. We used wild house mice (Mus musculus domesticus) from Gough Island (GI), which diverged recently from their mainland counterparts, to characterize the genetics of recombination rate evolution. We quantified genome-wide autosomal recombination rates by immunofluorescence cytology in spermatocytes from 240 F2 males generated from intercrosses between GI-derived mice and the wild-derived inbred strain WSB/EiJ. We identified four quantitative trait loci (QTL) responsible for inter-F2 variation in this trait, the strongest of which had effects that opposed the direction of the parental trait differences. Candidate genes and mutations for these QTL were identified by overlapping the detected intervals with whole-genome sequencing data and publicly available transcriptomic profiles from spermatocytes. Combined with existing studies, our findings suggest that genome-wide recombination rate divergence is not directional and its evolution within and between subspecies proceeds from distinct genetic loci.

Structural Variation Shapes the Landscape of Recombination in Mouse

Meiotic recombination is an essential feature of sexual reproduction that ensures faithful segregation of chromosomes and redistributes genetic variants in populations. Multiparent populations such as the Diversity Outbred (DO) mouse stock accumulate large numbers of crossover (CO) events between founder haplotypes, and thus present a unique opportunity to study the role of genetic variation in shaping the recombination landscape. We obtained high-density genotype data from [Formula: see text] DO mice, and localized 2.2 million CO events to intervals with a median size of 28 kb. The resulting sex-averaged genetic map of the DO population is highly concordant with large-scale (order 10 Mb) features of previously reported genetic maps for mouse. To examine fine-scale (order 10 kb) patterns of recombination in the DO, we overlaid putative recombination hotspots onto our CO intervals. We found that CO intervals are enriched in hotspots compared to the genomic background. However, as many as [Formula: see text] of CO intervals do not overlap any putative hotspots, suggesting that our understanding of hotspots is incomplete. We also identified coldspots encompassing 329 Mb, or [Formula: see text] of observable genome, in which there is little or no recombination. In contrast to hotspots, which are a few kilobases in size, and widely scattered throughout the genome, coldspots have a median size of 2.1 Mb and are spatially clustered. Coldspots are strongly associated with copy-number variant (CNV) regions, especially multi-allelic clusters, identified from whole-genome sequencing of 228 DO mice. Genes in these regions have reduced expression, and epigenetic features of closed chromatin in male germ cells, which suggests that CNVs may repress recombination by altering chromatin structure in meiosis. Our findings demonstrate how multiparent populations, by bridging the gap between large-scale and fine-scale genetic mapping, can reveal new features of the recombination landscape.

Germline and reproductive tract effects intensify in male mice with successive generations of estrogenic exposure

The hypothesis that developmental estrogenic exposure induces a constellation of male reproductive tract abnormalities is supported by experimental and human evidence. Experimental data also suggest that some induced effects persist in descendants of exposed males. These multi- and transgenerational effects are assumed to result from epigenetic changes to the germline, but few studies have directly analyzed germ cells. Typically, studies of transgenerational effects have involved exposing one generation and monitoring effects in subsequent unexposed generations. This approach, however, has limited human relevance, since both the number and volume of estrogenic contaminants has increased steadily over time, intensifying rather than reducing or eliminating exposure. Using an outbred CD-1 mouse model, and a sensitive and quantitative marker of germline development, meiotic recombination, we tested the effect of successive generations of exposure on the testis. We targeted the germline during a narrow, perinatal window using oral exposure to the synthetic estrogen, ethinyl estradiol. A complex three generation exposure protocol allowed us to compare the effects of individual, paternal, and grandpaternal (ancestral) exposure. Our data indicate that multiple generations of exposure not only exacerbate germ cell exposure effects, but also increase the incidence and severity of reproductive tract abnormalities. Taken together, our data suggest that male sensitivity to environmental estrogens is increased by successive generations of exposure.

X-Chromosome Control of Genome-Scale Recombination Rates in House Mice

Sex differences in recombination are widespread in mammals, but the causes of this pattern are poorly understood. Previously, males from two interfertile subspecies of house mice, Mus musculus musculus and M. m. castaneus, were shown to exhibit a ∼30% difference in their global crossover frequencies. Much of this crossover rate divergence is explained by six autosomal loci and a large-effect locus on the X chromosome. Intriguingly, the allelic effects at this X-linked locus are transgressive, with the allele conferring increased crossover rate being transmitted by the low crossover rate M. m. castaneus parent. Despite the pronounced divergence between males, females from these subspecies exhibit similar crossover rates, raising the question of how recombination is genetically controlled in this sex. Here, I analyze publicly available genotype data from early generations of the Collaborative Cross, an eight-way panel of recombinant inbred strains, to estimate crossover frequencies in female mice with sex-chromosome genotypes of diverse subspecific origins. Consistent with the transgressive influence of the X chromosome in males, I show that females inheriting an M. m. castaneus X possess higher average crossover rates than females lacking the M. m. castaneus X chromosome. The differential inheritance of the X chromosome in males and females provides a simple genetic explanation for sex-limited evolution of this trait. Further, the presence of X-linked and autosomal crossover rate modifiers with antagonistic effects hints at an underlying genetic conflict fueled by selection for distinct crossover rate optima in males and females.

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.

Recombination rate variation in mice from an isolated island

Recombination rate is a heritable trait that varies among individuals. Despite the major impact of recombination rate on patterns of genetic diversity and the efficacy of selection, natural variation in this phenotype remains poorly characterized. We present a comparison of genetic maps, sampling 1212 meioses, from a unique population of wild house mice (Mus musculus domesticus) that recently colonized remote Gough Island. Crosses to a mainland reference strain (WSB/EiJ) reveal pervasive variation in recombination rate among Gough Island mice, including subchromosomal intervals spanning up to 28% of the genome. In spite of this high level of polymorphism, the genomewide recombination rate does not significantly vary. In general, we find that recombination rate varies more when measured in smaller genomic intervals. Using the current standard genetic map of the laboratory mouse to polarize intervals with divergent recombination rates, we infer that the majority of evolutionary change occurred in one of the two tested lines of Gough Island mice. Our results confirm that natural populations harbour a high level of recombination rate polymorphism and highlight the disparities in recombination rate evolution across genomic scales.

Variation and Evolution of the Meiotic Requirement for Crossing Over in Mammals

The segregation of homologous chromosomes at the first meiotic division is dependent on the presence of at least one well-positioned crossover per chromosome. In some mammalian species, however, the genomic distribution of crossovers is consistent with a more stringent baseline requirement of one crossover per chromosome arm. Given that the meiotic requirement for crossing over defines the minimum frequency of recombination necessary for the production of viable gametes, determining the chromosomal scale of this constraint is essential for defining crossover profiles predisposed to aneuploidy and understanding the parameters that shape patterns of recombination rate evolution across species. Here, I use cytogenetic methods for in situ imaging of crossovers in karyotypically diverse house mice (Mus musculus domesticus) and voles (genus Microtus) to test how chromosome number and configuration constrain the distribution of crossovers in a genome. I show that the global distribution of crossovers in house mice is thresholded by a minimum of one crossover per chromosome arm, whereas the crossover landscape in voles is defined by a more relaxed requirement of one crossover per chromosome. I extend these findings in an evolutionary metaanalysis of published recombination and karyotype data for 112 mammalian species and demonstrate that the physical scale of the genomic crossover distribution has undergone multiple independent shifts from one crossover per chromosome arm to one per chromosome during mammalian evolution. Together, these results indicate that the chromosomal scale constraint on crossover rates is itself a trait that evolves among species, a finding that casts light on an important source of crossover rate variation in mammals.

Genetic Background, Maternal Age, and Interaction Effects Mediate Rates of Crossing Over in Drosophila melanogaster Females

Meiotic recombination is a genetic process that is critical for proper chromosome segregation in many organisms. Despite being fundamental for organismal fitness, rates of crossing over vary greatly between taxa. Both genetic and environmental factors contribute to phenotypic variation in crossover frequency, as do genotype-environment interactions. Here, we test the hypothesis that maternal age influences rates of crossing over in a genotypic-specific manner. Using classical genetic techniques, we estimated rates of crossing over for individual Drosophila melanogaster females from five strains over their lifetime from a single mating event. We find that both age and genetic background significantly contribute to observed variation in recombination frequency, as do genotype-age interactions. We further find differences in the effect of age on recombination frequency in the two genomic regions surveyed. Our results highlight the complexity of recombination rate variation and reveal a new role of genotype by maternal age interactions in mediating recombination rate.

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.

Homoploid hybrid speciation and genome evolution via chromosome sorting

Genomes of numerous diploid plant and animal species possess traces of interspecific crosses, and many researches consider them as support for homoploid hybrid speciation (HHS), a process by which a new reproductively isolated species arises through hybridization and combination of parts of the parental genomes, but without an increase in ploidy. However, convincing evidence for a creative role of hybridization in the origin of reproductive isolation between hybrid and parental forms is extremely limited. Here, through studying Agrodiaetus butterflies, we provide proof of a previously unknown mode of HHS based on the formation of post-zygotic reproductive isolation via hybridization of chromosomally divergent parental species and subsequent fixation of a novel combination of chromosome fusions/fissions in hybrid descendants. We show that meiotic segregation, operating in the hybrid lineage, resulted in the formation of a new diploid genome, drastically rearranged in terms of chromosome number. We also demonstrate that during the heterozygous stage of the hybrid species formation, recombination was limited between rearranged chromosomes of different parental origin, representing evidence that the reproductive isolation was a direct consequence of hybridization.

Correlations between Synaptic Initiation and Meiotic Recombination: A Study of Humans and Mice

Meiotic recombination is initiated by programmed double strand breaks (DSBs), only a small subset of which are resolved into crossovers (COs). The mechanism determining the location of these COs is not well understood. Studies in plants, fungi, and insects indicate that the same genomic regions are involved in synaptic initiation and COs, suggesting that early homolog alignment is correlated with the eventual resolution of DSBs as COs. It is generally assumed that this relationship extends to mammals, but little effort has been made to test this idea. Accordingly, we conducted an analysis of synaptic initiation sites (SISs) and COs in human and mouse spermatocytes and oocytes. In contrast to our expectation, we observed remarkable sex- and species-specific differences, including pronounced differences between human males and females in both the number and chromosomal location of SISs. Further, the combined data from our studies in mice and humans suggest that the relationship between SISs and COs in mammals is a complex one that is not dictated by the sites of synaptic initiation as reported in other organisms, although it is clearly influenced by them.

No Evidence that Infection Alters Global Recombination Rate in House Mice

Recombination rate is a complex trait, with genetic and environmental factors shaping observed patterns of variation. Although recent studies have begun to unravel the genetic basis of recombination rate differences between organisms, less attention has focused on the environmental determinants of crossover rates. Here, we test the effect of one ubiquitous environmental pressure-bacterial infection-on global recombination frequency in mammals. We applied MLH1 mapping to assay global crossover rates in male mice infected with the pathogenic bacterium Borrelia burgdorferi, the causative agent of Lyme Disease, and uninfected control animals. Despite ample statistical power to identify biologically relevant differences between infected and uninfected animals, we find no evidence for a global recombination rate response to bacterial infection. Moreover, broad-scale patterns of crossover distribution, including the number of achiasmate bivalents, are not affected by infection status. Although pathogen exposure can plastically increase recombination in some species, our findings suggest that recombination rates in house mice may be resilient to at least some forms of infection stress. This negative result motivates future experiments with alternative house mouse pathogens to evaluate the generality of this conclusion.

Cattle Sex-Specific Recombination and Genetic Control from a Large Pedigree Analysis

Meiotic recombination is an essential biological process that generates genetic diversity and ensures proper segregation of chromosomes during meiosis. From a large USDA dairy cattle pedigree with over half a million genotyped animals, we extracted 186,927 three-generation families, identified over 8.5 million maternal and paternal recombination events, and constructed sex-specific recombination maps for 59,309 autosomal SNPs. The recombination map spans for 25.5 Morgans in males and 23.2 Morgans in females, for a total studied region of 2,516 Mb (986 kb/cM in males and 1,085 kb/cM in females). The male map is 10% longer than the female map and the sex difference is most pronounced in the subtelomeric regions. We identified 1,792 male and 1,885 female putative recombination hotspots, with 720 hotspots shared between sexes. These hotspots encompass 3% of the genome but account for 25% of the genome-wide recombination events in both sexes. During the past forty years, males showed a decreasing trend in recombination rate that coincided with the artificial selection for milk production. Sex-specific GWAS analyses identified PRDM9 and CPLX1 to have significant effects on genome-wide recombination rate in both sexes. Two novel loci, NEK9 and REC114, were associated with recombination rate in both sexes, whereas three loci, MSH4, SMC3 and CEP55, affected recombination rate in females only. Among the multiple PRDM9 paralogues on the bovine genome, our GWAS of recombination hotspot usage together with linkage analysis identified the PRDM9 paralogue on chromosome 1 to be associated in the U.S. Holstein data. Given the largest sample size ever reported for such studies, our results reveal new insights into the understanding of cattle and mammalian recombination.

The blue butterfly Polyommatus (Plebicula) atlanticus (Lepidoptera, Lycaenidae) holds the record of the highest number of chromosomes in the non-polyploid eukaryotic organisms

The blue butterfly species Polyommatus (Plebicula) atlanticus (Elwes, 1906) (Lepidoptera, Lycaenidae) is known to have a very high haploid number of chromosomes (n= circa 223). However, this approximate count made by Hugo de Lesse 45 years ago was based on analysis of a single meiotic I metaphase plate, not confirmed by study of diploid chromosome set and not documented by microphotographs. Here I demonstrate that (1) Polyommatus atlanticus is a diploid (non-polyploid) species, (2) its meiotic I chromosome complement includes at least 224-226 countable chromosome bodies, and (3) all (or nearly all) chromosome elements in meiotic I karyotype are represented by bivalents. I also provide the first data on the diploid karyotype and estimate the diploid chromosome number as 2n=ca448-452. Thus, Polyommatus atlanticus is confirmed to possess the highest chromosome number among all the non-polyploid eukaryotic organisms.

Effects of crossovers between homeologs on inheritance and population genomics in polyploid-derived salmonid fishes

A whole genome duplication occurred in the ancestor of all salmonid fishes some 50-100 million years ago. Early inheritance studies with allozymes indicated that loci in the salmonid genome are inherited disomically in females. However, some pairs of duplicated loci showed patterns of inheritance in males indicating pairing and recombination between homeologous chromosomes. Nearly 20% of loci in the salmonid genome are duplicated and share the same alleles (isoloci), apparently due to homeologous recombination. Half-tetrad analysis revealed that isoloci tend to be telomeric. These results suggested that residual tetrasomic inheritance of isoloci results from homeologous recombination near chromosome ends and that continued disomic inheritance resulted from homologous pairing of centromeric regions. Many current genetic maps of salmonids are based on single nucleotide polymorphisms and microsatellites that are no longer duplicated. Therefore, long sections of chromosomes on these maps are poorly represented, especially telomeric regions. In addition, preferential multivalent pairing of homeologs from the same species in F1 hybrids results in an excess of nonparental gametes (so-called pseudolinkage). We consider how not including duplicated loci has affected our understanding of population and evolutionary genetics of salmonids, and we discuss how incorporating these loci will benefit our understanding of population genomics.

Sex-specific recombination maps for individual macrochromosomes in the Japanese quail (Coturnix japonica)

Meiotic recombination in the Japanese quail was directly studied by immunolocalization of mutL homolog 1 (MLH1), a mismatch repair protein of mature recombination nodules. In total, 15,862 crossovers were scored along the autosomal synaptonemal complexes in 308 meiotic nuclei from males and females. Crossover frequencies calculated from MLH1 foci show wide similitude between males and females with slightly higher number of foci in females. From this analysis, we predict that the sex-averaged map length of the Japanese quail is 2580 cM, with a genome-wide recombination rate of 1.9 cM/Mb. MLH1 focus mapping along the six largest bivalents showed few intersex differences in the distribution of crossovers along with variant patterns in metacentric and acrocentric macrobivalents. These results provide valuable information to complement linkage map analysis in the species while providing insight into our understanding of the mechanisms of crossover distribution along chromosome arms.

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.