Genetic recombination is directed away from functional genomic elements in mice

Genetic recombination occurs during meiosis, the key developmental program of gametogenesis. Recombination in mammals has been recently linked to the activity of a histone H3 methyl-transferase, PRDM91–6, the product of the only known speciation gene in mammals7. PRDM9 is thought to determine the preferred recombination sites – recombination hotspots – through sequence-specific binding of its highly polymorphic multi-Zn-finger domain8. Nevertheless, Prdm9 knockout mice are proficient at initiating recombination 9. Here we map and analyze the genome-wide distribution of recombination initiation sites in Prdm9 knockout mice and in two mouse strains with different Prdm9 alleles and their F1 hybrid. We show that PRDM9 determines the positions of practically all hotspots in the mouse genome, with the remarkable exception of the pseudoautosomal region – the only area of the genome that undergoes recombination in 100% of cells10. Surprisingly, hotspots are still observed in Prdm9 knockout mice and as in wild-type, these hotspots are found at H3K4 trimethylation marks. However, in the absence of PRDM9, the majority of recombination is initiated at promoters and at other sites of PRDM9-independent H3K4 trimethylation. Such sites are rarely targeted in wild-type mice indicating an unexpected role of the PRDM9 protein in sequestering the recombination machinery away from gene promoter regions and other functional genomic elements.

assess the extent to which hotspot location depends on PRDM9 we have employed this approach to generate high-resolution maps of recombination initiation sites in two nearly identical (congenic) mouse strains with different Prdm9 alleles (strains 9R and 13R, Supplementary Fig. 1). We found that although the number of hotspots in 9R and 13R is similar (14,869 and 15,481, respectively) there is practically no overlap in the locations of hotspots (Fig. 1a,b). We next mapped recombination initiation sites in the most commonly used laboratory mouse strain C57Bl/6 (18,313 hotspots), which has the same Prdm9 allele as 9R. We found that more than 98% of 9R hotspots were present in C57Bl/6 ( Fig. 1c). Taken together, these data constitute the first direct evidence that practically all recombination hotspots are PRDM9 dependent.
We next assessed the number and distribution of recombination hotspots in mice heterozygous for functional Prdm9 alleles, as such mice would have potentially twice as many proficient PRDM9 binding sites. We found that in F1 hybrids derived from a 9R × 13R cross the number of detected hotspots (15,677) is comparable to those in the parental strains (see Supplementary Methods). Strikingly, although the vast majority of hotspots in the F1 mice coincided with hotspots in the parental lines, 75% corresponded to 13R hotspots, while only 22% corresponded to those of 9R (Fig. 1d). Hotspots from the reciprocal cross showed a similar bias (see Supplementary Methods). The 13R-derived hotspots in F1 were also significantly stronger than those of 9R origin (median of 575 and 280 reads, respectively; p = 10 −182 , Wilcoxon rank sum test), but the relative strengths of hotspots within parental pools were retained in the F1 background. Importantly, F1 hotspots composed a stronger than average subset of hotspots in both respective parental strains (Fig.  1e). In aggregate, these data indicate that the number of recombination hotspots does not double in the presence of two different Prdm9 alleles, but is rather maintained at a constant level with mostly the strongest of the potential hotspots being utilized. The dominance of the 13R allele could result from higher affinity of the 13R PRDM9 protein variant to its preferred DNA sequence however other scenarios are also possible.
The complete discordance between the 9R and 13R hotspot maps suggested that different DNA sequences might be recognized by the 9R and 13R PRDM9 protein variants. Indeed, we identified a distinct motif co-centred with hotspots in each strain. These motifs are present in nearly all hotspots (9R: 94% ; 13R: 97%) and enriched at hotspot centres 4 and 8fold, respectively (Fig. 1f). The enrichment of better matches to motifs reaches 200-300-fold ( Supplementary Fig. 2), and motif quality correlates with hotspot strength ( Supplementary  Fig. 3). We next predicted the putative DNA binding sites for the Zn-finger arrays of the 9R and 13R PRDM9 protein variants and found that each motif aligns with the predicted binding site for its respective Prdm9 allele but not with that of the other (Fig. 1g). Both human 2 and mouse 12 hotspot motifs were previously found to align to the C-terminal part of the PRDM9 Zn finger array ( Supplementary Fig. 4). Our data suggest that DNA binding is not constrained to this region as the 9R motif aligns to N-terminal Zn fingers (Fig. 1g). Although these data are consistent with the sequence specificity of PRDM9 being a key determinant of hotspot positions in vivo, only a fraction of potential PRDM9 binding sites are utilized for hotspot formation. It is possible that sequence motifs do not capture some aspects of the binding site information such as di-or tri-nucleotide preferences of the PRDM9 protein. Since PRDM9 (or PRDM9-containing complex) must act within the constraints of the higher order chromatin organisation, it is also conceivable that the affinity of PRDM9 for its recognition sequence depends on the local chromatin environment.
In mammals, PRDM9 is the only histone H3 lysine 4 trimethyltransferase that has been reported to be specific to meiosis 9 . Since recombination hotspots are associated with meiosis-specific H3K4me3 marks 12,14,15 , the PRDM9 protein is likely responsible for their introduction 8 . If so, recombination hotspots in the strains with different Prdm9 alleles will correspond to H3K4me3 marks that are strain-specific. Indeed, we found that the vast majority of 9R and 13R hotspots overlap strain-specific H3K4me3 marks (9R: ≥83%; 13R: ≥89% - Supplementary Fig. 5a,b). Analysis of H3K4me3 in a Spo11 −/− mouse revealed that H3K4me3 marks are present at 89% of potential hotspot loci despite the lack of SPO11induced DSBs ( Supplementary Fig. 5c), and in agreement with previous findings 14, 15 . This indicates that hotspot-associated H3K4 trimethylation is not a consequence of DSB formation and/or repair, but a mark that is present prior to initiation of recombination.
Despite the clear critical role of PRDM9 in the initiation of genetic recombination meiotic DSBs are present in Prdm9 knockout mice 9 . Remarkably, we found that initiation of recombination in these mice is not random or uniform, but is still clustered in hotspots. The vast majority (99%) of these hotspots do not overlap the hotspots detected in any of the wild type strains, nevertheless, 94% of Prdm9 −/− hotspots still overlap H3K4me3 marks. Most of these marks (92%) are present in wild type testis and almost half are not specific to germline cells ( Supplementary Fig. 5d). As H3K4me3 marks outnumber DSB hotspots, they do not appear to be sufficient for hotspot formation ( Supplementary Fig. 6). H3K4me3 is a general mark of active transcription primarily associated with gene promoters 16 , enhancers 17,18 and possibly other functional genomic elements. We found that in the absence of PRDM9 and PRDM9-introduced H3K4me3 marks, recombination hotspots are re-routed to these alternative H3K4me3 sites (Fig. 2a). Almost half (44%) of recombination hotspots in the Prdm9 −/− mice localize to the promoters of annotated genes compared to just 3% in wild type. Unlike in wild type mice, promoter-associated hotspots in Prdm9 −/− are relatively strong ( Supplementary Fig. 7) and therefore, account for the majority of DSB-derived sequencing reads (62%; Fig. 2b). The H3K4me3 signal at promoters that overlap hotspots in Prdm9 −/− is almost 100-fold stronger than at those that do not (median: 350 tags and 5 tags, respectively, P < 10 −200 , Wilcoxon test). Furthermore, promoter-overlapping hotspots are slightly enriched at genes expressed in early meiotic prophase (P = 10 −138 , one-sided binomial test; Supplementary Fig. 8). Prdm9 −/− hotspots that are not associated with annotated promoters still likely represent DSBs formed at important genomic elements as they exhibit higher than expected sequence conservation ( Supplementary Fig. 9). Indeed, 40% of these hotspots overlap histone crotonylation sites, a recently discovered marker of promoters and enhancers 19 (see Supplementary Methods). The distribution of hotspots around transcription start sites (TSSs) in Prdm9 −/− mice is reminiscent of that in S. cerevisiae where most hotspots are found at promoters 20,21 . However, while DSBs at yeast promoters form in the upstream, nucleosome-depleted region 20,21 , hotspots in Prdm9 −/− mice form preferentially at the most frequently H3K4 trimethylated nucleosome at the +1 position, just downstream of the TSS (Fig. 2c).
Unique in the genome is a short region of homology at the end of X and Y chromosome called the pseudo-autosomal region (PAR) that undergoes recombination in every meiosis 10 . We have recently demonstrated that the only sequenced region of the mouse PAR, adjacent to the interstitial PAR boundary, contains a large cluster of overlapping DSB hotspots spanning at least 40 Kb 12 . We now find that this hotspot cluster is still present in Prdm9 −/− mice, and is shared between the strains with different Prdm9 alleles (Fig. 3). Akin to the hotspot cluster, the H3K4me3 signal in the PAR is far broader and stronger than any other H3K4me3 marks in the genome, but it is shifted slightly off the hotspot cluster, to within the PAR (Fig. 3). Recent studies have demonstrated that the structural organization of the PAR differs significantly from that of the rest of the chromosome 22 . It is conceivable that this region features constitutively accessible chromatin that allows initiation of recombination without participation of the PRDM9 protein. Initiation of homologous recombination at some of the hotspots in the PAR-flanking region is also defined by a PRDM9-independent mechanism. In this region, up to ~900 Kb upstream of PAR, there is an extraordinary concentration of hotspots shared between Prdm9 −/− mice and wild type strains with different Prdm9 alleles (C57Bl/6, 9R and 13R; Fig. 3). As expected, the H3K4me3 marks at these shared hotspots are also strain independent. The density of PRDM9-independent hotspots in the PAR-flanking region appears to increase with the proximity to the PAR (Fig.  3), consistent with a specific chromatin organization or particular epigenetic mark characteristic to the PAR propagating into the adjacent area. Yet, PRDM9-dependent hotspots are also present in the PAR-flanking area, therefore the DNA in this region is still accessible to the PRDM9 protein.
Why is the initiation of recombination in the PAR different from the rest of the genome? To ensure proper sex chromosome segregation during meiosis, at least one crossover must be formed inside the PAR 10 . Out of thousands of PRDM9-dependent hotspots only 200-400 are utilized in each meiosis and only 20-25 result in crossovers. Considering the short length of the PAR, a large fraction of cells would fail to introduce a crossover within this area simply by chance. It is possible that the mandatory crossover requirement is more efficiently fulfilled through an alternative mechanism of initiation, unique to the PAR. PRDM9independent initiation of homologous recombination in the PAR would also facilitate the rapid evolution of both PRDM9 23 and the PAR 24 by allowing them to evolve independently of each other. It is also conceivable that PRDM9-dependent and -independent recombination pathways utilize alternative isoforms of SPO11, as the number of meiotic DSBs in the PAR is diminished in mice expressing only the beta SPO11 isoform 22 .
Our study clearly demonstrates that the PRDM9 protein is not required for the initiation of genetic recombination per se, but rather, determines the positions (hotspots) where this initiation takes place. In the absence of PRDM9 the vast majority of DSB hotspots still coincide with H3K4me3 marks. Although a causal relationship is yet to be established, this suggests that H3K4me3 or/and mechanisms that induce H3K4me3 are necessary to provide a favourable chromatin environment for the initiation of meiotic recombination. Such an environment can be created with or without PRDM9, but only PRDM9-dependent H3K4me3 sites are associated with recombination hotspots in wild-type mice outside of the PAR. This indicates that PRDM9 is directly involved in recruitment of the recombination initiation machinery.
The sterility of Prdm9 −/− mice suggests that the initiation of recombination at H3K4me3marked functional elements may not be compatible with gametogenesis. Does recombination interfere with the transcription of essential meiotic genes? Since four chromatids are present before the first meiotic division it is unlikely that disruption of transcription on one of them will lead to a significant drop in overall transcription of any particular gene hit by a DSB. On the other hand, it is conceivable that the opposite is true and that transcription interferes with the repair of breaks introduced at promoters. The resulting aberrant recombination would account for the observation that in Prdm9 −/− mice persistent DSBs are present and only partial synapsis is seen between homologous chromosomes 9 ( Supplementary Fig. 10). The Prdm9 gene exhibits complex epistatic interactions culminating in sterility of F1 hybrid males carrying a particular combination of functional Prdm9 alleles 7 . Unlike in Prdm9 −/− mice, repair of meiotic DSBs in such males proceeds normally in a large fraction of cells with 40% of spermatocytes showing complete homologous synapsis 7 . While this suggests that distinct mechanisms are responsible for meiotic arrest in Prdm9 −/− mice and in F1 hybrid males, we cannot exclude the possibility that an alternative role of PRDM9 in meiotic progression accounts for both sterility phenotypes. Although a function of PRDM9 in transcription has not been demonstrated, PRDM9 contains conserved domains implicated in transcription 25,26 and could potentially be involved in executing the meiotic transcription program. It is also conceivable that PRDM9 is directly required for efficient progression of recombination. The potential importance of avoiding recombination at functional genomic elements is highlighted in the canid lineage, which has lost Prdm9 27,28 . Though dogs retain recombination hotspots 27 , preliminary analyses suggest that unlike in Prdm9 −/− mice, these hotspots do not appear to be enriched at promoters (Supplementary Information). A trivial explanation would be that dog Prdm9 resides in the unassembled parts of the genome, however, these data may also indicate that alternative mechanisms have evolved to target recombination to specialized genomic regions in species that lack Prdm9. Ultimately, we favour the hypothesis that PRDM9 prevents initiation of recombination in the areas of the genome where successful completion of recombination might be compromised by other nuclear processes (Fig. 4).
The re-routing of DSBs from important genomic elements may also play a protective role against potential mutagenic effects of recombination.

Methods Summary
SSDS 13 using antibodies against DMC1 was used to identify DSB hotspots in the following mice: B10.S-H2t4/(9R)/J, B10.F-H2pb1/(13R)J, B10.S-H2t4/(9R)/J×B10.F-H2pb1/(13R)J (F1), C57Bl/6J and B6;129P2-Prdm9 tm1Ymat /J. Chromatin immunoprecipitation using antibodies against H3K4me3 followed by Illumina sequencing was also performed on samples derived from each of these mice. DSB hotspots were defined by SISSRS 29 and post-processed using a bespoke pipeline. Peaks in H3K4me3 data were called using MACS v1.3.7 30 . Detailed methods are available in Supplementary Information.  PRDM9 redirects DSBs away from functional genomic elements. a. In Prdm9 −/− mice DSBs are formed at functional genomic elements marked by H3K4me3 (such as promoters, red boxes). These sites are refractory to DSB formation in wild type mice, where breaks are redirected to PRDM9-dependent H3K4me3 marks (*). b. Most DSBs occur in regions around TSSs in Prdm9 −/− mice but not in wild type mice. c. The centre of Prdm9 −/− hotspots at TSS coincides precisely with the H3K4me3 mark on the +1 nucleosome (red circle: −1 nucleosome; white circle: nucleosome free region; green circle: +1 nucleosome). Mean H3K4me3 coverage is grey, mean forward strand ssDNA coverage is blue and reverse strand coverage is red. The intersection point of the forward and reverse strand coverage is the mean hotspot centre ( Supplementary Fig. 11).