Sex-specific variation in the genome-wide recombination rate

Karyotypic stasis and swarming influenced the evolution of viral tolerance in a species-rich bat radiation

The emergence of COVID-19 and severe acute respiratory syndrome (SARS) has prioritized understanding bats' viral tolerance. Myotis bats are exceptionally species rich and have evolved viral tolerance. They also exhibit swarming, a cryptic behavior where large, multi-species assemblages gather for mating, which has been hypothesized to promote interspecific hybridization. To resolve the coevolution of genome architecture and their unusual antiviral tolerance, we undertook a phylogenomic analysis of 60 Old World Myotis genomes. We demonstrate an extensive history of introgressive hybridization that has replaced the species phylogeny across 17%-93% of the genome except for pericentromeric regions of macrochromosomes. Introgression tracts were enriched on microchromosome regions containing key antiviral pathway genes overexpressed during viral challenge experiments. Together, these results suggest that the unusual Myotis karyotype may have evolved to selectively position immune-related genes in high recombining genomic regions prone to introgression of divergent alleles, including a diversity of interleukin loci responsible for the release of pro-inflammatory cytokines.

The DNA helicase FANCJ (BRIP1) functions in double strand break repair processing, but not crossover formation during prophase I of meiosis in male mice

Meiotic recombination between homologous chromosomes is initiated by the formation of hundreds of programmed double-strand breaks (DSBs). Approximately 10% of these DSBs result in crossovers (COs), sites of physical DNA exchange between homologs that are critical to correct chromosome segregation. Virtually all COs are formed by coordinated efforts of the MSH4/MSH5 and MLH1/MLH3 heterodimers, the latter representing the defining marks of CO sites. The regulation of CO number and position is poorly understood, but undoubtedly requires the coordinated action of multiple repair pathways. In a previous report, we found gene-trap disruption of the DNA helicase, FANCJ (BRIP1/BACH1), elicited elevated numbers of MLH1 foci and chiasmata. In somatic cells, FANCJ interacts with numerous DNA repair proteins including MLH1, and we hypothesized that FANCJ functions with MLH1 to regulate the major CO pathway. To further elucidate the meiotic function of FANCJ, we produced three new Fancj mutant mouse lines via CRISPR/Cas9 gene editing: a full-gene deletion, truncation of the N-terminal Helicase domain, and a C-terminal dual-tagged allele. We also generated an antibody against the C-terminus of the mouse FANCJ protein. Surprisingly, none of our Fancj mutants show any change in either MLH1 focus counts during pachynema or total CO number at diakinesis of prophase I. We find evidence that FANCJ and MLH1 do not interact in meiosis; further, FANCJ does not co-localize with MSH4, MLH1, or MLH3 in meiosis. Instead, FANCJ co-localizes with BRCA1 and TOPBP1, forming discrete foci along the chromosome cores beginning in early meiotic prophase I and densely localized to unsynapsed chromosome axes in late zygonema and to the XY chromosomes in early pachynema. Fancj mutants also exhibit a subtle persistence of DSBs in pachynema. Collectively, these data indicate a role for FANCJ in early DSB repair, but they rule out a role for FANCJ in MLH1-mediated CO events.

Modeling Recombination Rate as a Quantitative Trait Reveals New Insight into Selection in Humans

Meiotic recombination is both a fundamental biological process required for proper chromosomal segregation during meiosis and an important genomic parameter that shapes major features of the genomic landscape. However, despite the central importance of this phenotype, we lack a clear understanding of the selective pressures that shape its variation in natural populations, including humans. While there is strong evidence of fitness costs of low rates of recombination, the possible fitness costs of high rates of recombination are less defined. To determine whether a single lower fitness bound can explain the variation in recombination rates observed in human populations, we simulated the evolution of recombination rates as a sexually dimorphic quantitative trait. Under each scenario, we statistically compared the resulting trait distribution with the observed distribution of recombination rates from a published study of the Icelandic population. To capture the genetic architecture of recombination rates in humans, we modeled it as a moderately complex trait with modest heritability. For our fitness function, we implemented a hyperbolic tangent curve with several flexible parameters to capture a wide range of existing hypotheses. We found that costs of low rates of recombination alone are likely insufficient to explain the current variation in recombination rates in both males and females, supporting the existence of fitness costs of high rates of recombination in humans. With simulations using both upper and lower fitness boundaries, we describe a parameter space for the costs of high recombination rates that produces results consistent with empirical observations.

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.

Gene flow biases population genetic inference of recombination rate

Accurate estimates of the rate of recombination are key to understanding a host of evolutionary processes as well as the evolution of the recombination rate itself. Model-based population genetic methods that infer recombination rates from patterns of linkage disequilibrium in the genome have become a popular method to estimate rates of recombination. However, these linkage disequilibrium-based methods make a variety of simplifying assumptions about the populations of interest that are often not met in natural populations. One such assumption is the absence of gene flow from other populations. Here, we use forward-time population genetic simulations of isolation-with-migration scenarios to explore how gene flow affects the accuracy of linkage disequilibrium-based estimators of recombination rate. We find that moderate levels of gene flow can result in either the overestimation or underestimation of recombination rates by up to 20–50% depending on the timing of divergence. We also find that these biases can affect the detection of interpopulation differences in recombination rate, causing both false positives and false negatives depending on the scenario. We discuss future possibilities for mitigating these biases and recommend that investigators exercise caution and confirm that their study populations meet assumptions before deploying these methods.

Background selection under evolving recombination rates

Background selection (BGS), the effect that purifying selection exerts on sites linked to deleterious alleles, is expected to be ubiquitous across eukaryotic genomes. The effects of BGS reflect the interplay of the rates and fitness effects of deleterious mutations with recombination. A fundamental assumption of BGS models is that recombination rates are invariant over time. However, in some lineages, recombination rates evolve rapidly, violating this central assumption. Here, we investigate how recombination rate evolution affects genetic variation under BGS. We show that recombination rate evolution modifies the effects of BGS in a manner similar to a localized change in the effective population size, potentially leading to underestimation or overestimation of the genome-wide effects of selection. Furthermore, we find evidence that recombination rate evolution in the ancestors of modern house mice may have impacted inferences of the genome-wide effects of selection in that species.

Male and female recombination landscapes of diploid Arabidopsis arenosa

The number and placement of meiotic crossover events during meiosis have important implications for the fidelity of chromosome segregation as well as patterns of inheritance. Despite the functional importance of recombination, recombination landscapes vary widely among and within species, and this can have a strong impact on evolutionary processes. A good knowledge of recombination landscapes is important for model systems in evolutionary and ecological genetics, since it can improve interpretation of genomic patterns of differentiation and genome evolution, and provides an important starting point for understanding the causes and consequences of recombination rate variation. Arabidopsis arenosa is a powerful evolutionary genetic model for studying the molecular basis of adaptation and recombination rate evolution. Here, we generate genetic maps for 2 diploid A. arenosa individuals from distinct genetic lineages where we have prior knowledge that meiotic genes show evidence of selection. We complement the genetic maps with cytological approaches to map and quantify recombination rates, and test the idea that these populations might have distinct patterns of recombination. We explore how recombination differs at the level of populations, individuals, sexes and genomic regions. We show that the positioning of crossovers along a chromosome correlates with their number, presumably a consequence of crossover interference, and discuss how this effect can cause differences in recombination landscape among sexes or species. We identify several instances of female segregation distortion. We found that averaged genome-wide recombination rate is lower and sex differences subtler in A. arenosa than in Arabidopsis thaliana.

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.

Higher Intercellular Variation in Genome-Wide Recombination Rate in Female Mice

Meiotic recombination affects fertility, shuffles genomes, and modulates the effectiveness of natural selection. Despite conservation of the recombination pathway, the rate of recombination varies among individuals and along chromosomes. Recombination rate also differs among cells from the same organism, but this form of variation has received less attention. To identify patterns that characterize intercellular variation in the genome-wide recombination rate, we counted foci of the MLH1 recombination-associated protein in oocytes and spermatocytes from a panel of wild-derived inbred strains of house mice. Females show higher intercellular variation in MLH1 focus count than males from the same inbred strains. This pattern is consistent across strains from multiple subspecies, including 2 strains in which the average MLH1 focus count is higher in males. The sex difference in genome-wide recombination rate we report suggests that selection targeting recombination rate will be more efficient in males than in females.