# The PHD finger protein Spp1 has distinct functions in the Set1 and the meiotic DSB formation complexes

Céline Adam, Raphaël Guérois, Anna Citarella, Laura Verardi, Florine Adolphe, Claire Béneut, Vérane Sommermeyer, Claire Ramus, Jérôme Govin, Yohann Couté, et al.

► To cite this version:

Céline Adam, Raphaël Guérois, Anna Citarella, Laura Verardi, Florine Adolphe, et al. The PHD finger protein Spp1 has distinct functions in the Set1 and the meiotic DSB formation complexes. PLoS Genetics, 2018, 14 (2), pp.e1007223. 10.1371/journal.pgen.1007223. hal-01738379

HAL Id: hal-01738379

https://hal.sorbonne-universite.fr/hal-01738379

Submitted on 20 Mar 2018

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

# The PHD finger protein Spp1 has distinct functions in the Set1 and the meiotic DSB formation complexes

Céline Adam<sup>1,2</sup>, Raphaël Guérois<sup>3</sup>, Anna Citarella<sup>1,2</sup>, Laura Verardi<sup>1,2</sup>, Florine Adolphe<sup>1,2</sup>, Claire Béneut<sup>1,2</sup>, Vérane Sommermeyer<sup>1,2</sup>, Claire Ramus<sup>4</sup>, Jérôme Govin<sup>4</sup>, Yohann Coute<sup>4</sup>, Valérie Borde<sup>1,2*</sup>

1 Institut Curie, PSL Research University, CNRS, UMR3244, Paris, France, 2 Université Pierre et Marie Curie (UPMC), Paris, France, 3 Institute for Integrative Biology of the Cell (I2BC), CEA, CNRS, Université Paris-Sud, Université Paris-Saclay, Gif-sur-Yvette, France, 4 Univ. Grenoble Alpes, CEA, INSERM, BIG-BGE, Grenoble, France

* valerie.borde@curie.fr

<figure><img src="image_2.png" /></figure>

## Abstract

Histone H3K4 methylation is a feature of meiotic recombination hotspots shared by many organisms including plants and mammals. Meiotic recombination is initiated by programmed double-strand break (DSB) formation that in budding yeast takes place in gene promoters and is promoted by histone H3K4 di/trimethylation. This histone modification is recognized by Spp1, a PHD finger containing protein that belongs to the conserved histone H3K4 methyltransferase Set1 complex. During meiosis, Spp1 binds H3K4me3 and interacts with a DSB protein, Mer2, to promote DSB formation close to gene promoters. How Set1 complex- and Mer2- related functions of Spp1 are connected is not clear. Here, combining genome-wide localization analyses, biochemical approaches and the use of separation of function mutants, we show that Spp1 is present within two distinct complexes in meiotic cells, the Set1 and the Mer2 complexes. Disrupting the Spp1-Set1 interaction mildly decreases H3K4me3 levels and does not affect meiotic recombination initiation. Conversely, the Spp1-Mer2 interaction is required for normal meiotic recombination initiation, but dispensable for Set1 complex-mediated histone H3K4 methylation. Finally, we provide evidence that Spp1 preserves normal H3K4me3 levels independently of the Set1 complex. We propose a model where Spp1 works in three ways to promote recombination initiation: first by depositing histone H3K4 methylation (Set1 complex), next by "reading" and protecting histone H3K4 methylation, and finally by making the link with the chromosome axis (Mer2-Spp1 complex). This work deciphers the precise roles of Spp1 in meiotic recombination and opens perspectives to study its functions in other organisms where H3K4me3 is also present at recombination hotspots.

<figure><img src="image_3.png" /></figure>

### OPEN ACCESS

**Citation:** Adam C, Guérois R, Citarella A, Verardi L, Adolphe F, Béneut C, et al. (2018) The PHD finger protein Spp1 has distinct functions in the Set1 and the meiotic DSB formation complexes. PLoS Genet 14(2): e1007223. https://doi.org/10.1371/journal.pgen.1007223

**Editor:** Michael Lichten, National Cancer Institute, UNITED STATES

**Received:** August 23, 2017

**Accepted:** January 25, 2018

**Published:** February 14, 2018

**Copyright:** © 2018 Adam et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** The ChIP-chip data generated in this study have been deposited in the Gene Expression Omnibus database, accession number GSE102790. Processed data for all chromosomes are provided in S3 Table. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository via the url: http://www.ebi.ac.uk/pride/archive/projects/PXD007590.

**Funding:** This work was supported by Institut Curie, CNRS, Labex DEEP (ANR-11-LBX-0044)
part of the IDEX Idex PSL (ANR-10-IDEX-0001-02 PSL), by ANR-12-BSV6-0009, and by Fondation pour la Recherche Médicale (DEP20131128517). Proteomic experiments were partly supported by the Proteomics French Infrastructure (ANR-10-INBS-08-01 grant). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing interests:** The authors have declared that no competing interests exist.

## Author summary

Meiotic recombination is a conserved pathway of sexual reproduction that is required to faithfully segregate homologous chromosomes and produce viable gametes. Recombination events between homologous chromosomes are triggered by the programmed formation of DNA breaks, which occur preferentially at places called hotspots. In many organisms, these hotspots are located close to a particular chromatin modification, the methylation of lysine 4 of histone H3 (H3K4me3). It was previously shown in the budding yeast model that one protein, Spp1, plays an important function in this process. We further explored the functional link between Spp1 and its interacting partners, and show that Spp1 shows genetically separable functions, by depositing the H3K4me3 mark on the chromatin, "reading" and protecting it, and linking it to the recombination proteins. We provide evidence that Spp1 is in distinct complexes to perform these functions. This work opens perspectives for understanding the process in other eukaryotes such as mammals, where most of the proteins involved are conserved.

## Introduction

In sexually reproducing organisms, recombination between homologous chromosomes at meiosis plays an important role to form gametes with a normal chromosome content. Meiotic recombination is initiated by the programmed formation of DNA double-strand breaks (DSBs) catalyzed by the conserved Spo11 protein together with largely conserved accessory DSB proteins [1-3] and is required to promote genetic diversity and accurate homolog segregation [4]. Histone modifications are key players of the chromatin organization. Among these, histone H3K4 trimethylation (H3K4me3) is able to recruit downstream effectors such as chromatin remodelers [5]. In meiosis, Set1, the subunit of the Set1 complex that catalyzes histone H3K4 methylation, is required for normal levels and distribution of meiotic DSBs [6, 7]. In budding yeast, the vast majority of meiotic DSBs occur in the ~ 140 bp nucleosome depleted regions (NDR) at gene promoters, close to nucleosomes bearing the histone H3K4me3 modification [6, 8]. The tendency of meiotic recombination to occur in gene promoters and H3K4me3 regions is conserved among many organisms such as dogs, plants and birds [9-11]. In mice and humans, meiotic recombination also occurs at H3K4me3 sites, but only those created by the meiotic-specific histone H3K4 methyltransferase, PRDM9 [12-14]. PRDM9 recognizes specific DNA sequences with its array of zinc finger motifs, where it methylates H3K4 and promotes recombination, away from gene promoters [15]. Remarkably, the mutation of *Prdm9* in mice redirects meiotic recombination events towards gene promoters and H3K4me3 [16] as if PRDM9 had a dominant role over the default promoter/histone H3K4me3 conserved pathway.

The link between histone H3K4 methylation and meiotic DSB formation has recently been explained in budding yeast by the role of the PHD finger protein, Spp1, in spatially linking DSB sites to the recombination initiation machinery [17, 18]. During meiotic prophase, chromosomes adopt a specific three-dimensional structure formed of chromatin loops anchored at their basis to a chromosome axis [19]. DSBs are formed into loop DNA sequences, whereas the DSB proteins are located on the chromosome axis, implying a spatial contact between these two physically distant chromosomal regions during DSB formation [3, 20-23]. In meiosis, Spp1, a member of the Set1 complex, is, like Set1, required for normal DSB levels [17, 18]. Spp1 is specifically important for H3K4 trimethylation, and in its absence, H3K4me3 levels are reduced to about 20% of wild type [18]. This has been attributed to Spp1 being important for opening the catalytic site of Set1 and allowing trimethylation [24]. Spp1 also interacts with
Mer2, one of the axis-associated DSB proteins required for DSB formation, and is preferen-
tially located on the chromosome axis [1, 17, 18, 23]. The PHD finger of Spp1 interacts also
with H3K4me2/me3 at +1 nucleosomes and is required for normal DSB formation, and thus
Spp1 makes the physical link between gene promoters close to H3K4me2/3 sites and the DSB
formation machinery [17, 18]. Thus, Spp1 may facilitate or stabilize the interaction between
these distant regions, which triggers DSB formation by Spo11, the protein that bears the cata-
lytic DSB forming activity [17, 18].

In vegetative cells, Spp1 belongs to the Set1 complex and its distribution mirrors that of
RNA pol II [18]. By contrast, in meiosis, the chromosomal distribution of Spp1 shows no spa-
tial correlation with that of RNA pol II [18], raising the question whether Spp1 is still part of
the Set1 complex in meiosis, and if so, how it distributes between the Set1 complex and the
DSB proteins. In addition, given that Spp1 is required for normal levels of H3K4me3, known
to recruit downstream chromatin remodelers, it is not clear as well if the functional role of
Spp1 within the Set1 complex for H3K4 methylation can be separated from its implication in
DSB formation through its interaction with Mer2.

In this paper, we show that Spp1 interacts both with the Set1 complex and Mer2 in meiotic cells. However, Set1 complex does not associate with chromosome axes in meiosis, and its sub-units do not interact with Mer2, revealing that Spp1 is present in two distinct complexes. Next, we show that surprisingly, the presence of Spp1 in the Set1 complex is not important for maintaining H3K4 trimethylation levels and that Spp1 acts independently of the Set1 complex to promote meiotic DSB formation. Finally, we show that a mutant of *MER2* that no longer interacts with Spp1 but binds normally to chromosome axes is impaired for DSB formation. This demonstrates that solely affecting Spp1 interaction with Mer2 is sufficient to impair DSB formation, independently of any H3K4 methylation-related change in chromatin. Finally this work is relevant for understanding meiotic DSB formation in mammals and other organisms, for which a mechanism linking H3K4me3 and the DSB machinery likely exists but has not yet been elucidated.

**Results**

Spp1 is associated with both the DSB protein Mer2 and the Set1 complex in meiosis

During vegetative growth, Spp1 is a member of the eight-unit Set1 complex, and is predomi-
nantly located at highly transcribed genes, consistent with the Set1 complex being associated
with elongating RNA pol II [18, 25]. However, in meiosis, Spp1 interacts with the Mer2 protein
and is predominantly located on chromosomes axes [18]. We thus investigated how Spp1 dis-
tributes between the Set1 complex and its interaction with the axis-associated Mer2 DSB pro-
tein in meiotic cells. For this, we affinity-purified Spp1-TAP from cells at 3.5 hr in meiotic
prophase, the expected time of DSB formation [26]. As shown before for Myc-tagged Spp1, the
Spp1-TAP fusion protein associated with chromosome axis sites during meiosis (Fig 1A) [18].
Spp1-TAP interacting proteins were purified and identified by label-free mass-spectrometry-
based quantitative proteomics. Analyses revealed that Spp1-TAP interacts with the whole Set1
complex (composed of Bre2, Sdc1, Spp1, Set1, Swd1, Swd2, Swd3, Shg1) and with Mer2 in
meiotic cells (Fig 1B and 1C and S1 Table). No other proteins significantly purified with Spp1.

Distinct requirements govern Spp1 association with highly transcribed
genes and with chromosome axis in meiotic cells

Our finding that Spp1 is associated with the Set1 complex in meiosis is in apparent contradic-
tion with the fact that Spp1 appears mainly associated with regions of axis attachment, largely
<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

**Fig 1.** Spp1 interacts with the Set1 complex and the Mer2 DSB protein in meiotic cells. (A) ChIP-qPCR of Spp1-TAP during meiosis, showing its association with the chromosome axis at 3–4 hr, the expected time for DSB formation. Strain: VBD1266. Sites used for negative control (NFT1 gene) and axis (chr 3, nt 232942 to 233010) are the same in all figures. (B) Silver-stained gel of TAP eluates performed at 3.5 hr in meiosis in untagged (ORD7339) or Spp1-TAP (VBD1266). (C) Mass spectrometry analysis of proteins pulled down by Spp1-TAP in meiosis (t = 3.5 hr). SPP1-TAP: VBD1266 –no tag: ORD7339. The volcano plot indicates in red (Mer2) and in blue (Set1 complex subunits) the proteins that significantly co-purify with Spp1 (log2(Fold change)>3 and log10(p-value)<4). The experiment was done in three biological replicates, and average Fold Change values over the untagged control are presented with their corresponding p-value (see Methods).

https://doi.org/10.1371/journal.pgen.1007223.g001

co-localizing with Mer2, as assessed by genome wide mapping [18]. The contribution of Mer2 to the localization of Spp1 was assessed by mapping genome-wide Spp1 binding sites in *mer2Δ* meiotic cells. Whereas in wild type, Spp1 binding was not correlated with that of RNA pol II, in the absence of *MER2*, Spp1 location became positively correlated with that of RNA pol II (Pcorr = 0.57) (Fig 2A). Indeed, in the absence of Mer2, the strongest Spp1 peaks were now at highly transcribed genes, similar to the distribution of Spp1 in vegetative cells, likely reflecting its association with the Set1 complex (Fig 2B).
<figure><img src="image_7.png" /><figcaption>**Fig 2.** Spp1 binds distinct Set1- and Mer2-dependent sites in meiotic cells. (A) Comparison of Spp1 binding profile with that of RNA polymerase II, in vegetative cells, or in meiosis (t = 3 hr), in wild-type or *mer2Δ* mutant. Spp1-Myc and RNA pol II binding data are from [18]. Spp1-Myc *mer2Δ*: VBD1220. ChIPchip (Chromatin immunoprecipitation on chip) profiles are shown for chromosome 10. Centromeres are indicated by green dots, and ratios are plotted after smoothing with a 1 kb window. (B) Mean mRNA levels at the 100 strongest Spp1 peaks or in the whole genome. Spp1 peaks were determined from the experiments shown in A (see Methods). mRNA levels are from published data in SK1 diploid vegetative cells [57] (upper panel) or SK1 meiotic cells at t = 4 hr [58] (lower panel). Error bars: S.E.M. (C) Spp1 and Rec8 binding profile in *SET1* or *set1Δ* cells in meiosis (t = 3 hr). Spp1-Myc *set1Δ*: VBD1209. *SET1* Spp1-Myc and Rec8 binding data are from [18]. ChIPchip profiles are shown for chromosome 10. Centromere is indicated as a green dot, and ratios are plotted after smoothing with a 1 kb window. (D) qPCR analysis of Spp1 binding to an axis site and two highly transcribed genes in meiosis, *ACS1* and *CIT2* (t = 3 hr). No tag: ORD7339; WT: VBD1187; *mer2Δ*: VBD1220; *set1Δ*: VBD1209. Each dot represents a biological replicate and the bar indicates the mean.</figcaption></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

https://doi.org/10.1371/journal.pgen.1007223.g002

The localization of Spp1 was confirmed by ChIP-qPCR in WT (wild-type) and *mer2Δ* cells (Fig 2D). In WT cells, Spp1 was present at two highly transcribed genes, *ACS1* and *CIT2*, like in *mer2Δ*, in addition to its association with chromosome axes (Fig 2D). By contrast, Spp1-Myc fully associated with meiotic chromosomes axes in the absence of *SET1*, as shown by its positive correlation with that of the Rec8 protein (Pcorr = 0.66), similar to the correlation between Spp1 and Rec8 in *SET1* cells (Pcorr = 0.69 [18]) (Fig 2C and S1 Fig). However, the signal at the two highly transcribed genes *ACS1* and *CIT2* was lost (Fig 2D).

We propose that in wild-type cells, the association of Spp1 with Mer2 occurs at the highly localized axis association sites, resulting in strong peaks that mask its association with chromatin through the Set1 complex, which may be more diffuse along chromatin, and less concentrated in strongly defined peaks. Our findings are consistent with the existence of two classes
of Spp1 binding sites, ones with axis sites, dependent on *MER2*, and ones with transcription sites, dependent on *SET1*.

## The Set1 complex does not bind chromosome axes during meiosis and does not associate with Mer2

We next asked if when Spp1 is associated with Mer2 on chromosomes axes, it is within the Set1 complex. To answer this, we mapped the sites of chromatin association of Set1 and Swd1, another subunit of the Set1 complex, in meiosis. Set1 and Swd1 bound the two highly transcribed genes, *ACS1* and *CIT2*, during meiosis, but only weakly a chromosome axis-associated site (Fig 3A and S2 Fig). This was confirmed genome-wide, where the location of Swd1 in

<figure><img src="image_10.png" /><figcaption>D Swd1-HA ChIP</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>E Swd1-HA IP</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>B Chromosome 10</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>F Mer2-Flag IP</figcaption></figure>

<figure><img src="image_14.png" /></figure>

<figure><img src="image_15.png" /></figure>

<figure><img src="image_16.png" /></figure>

**Fig 3.** Set1 complex Set1 and Swd1 subunits associate in meiosis with highly transcribed genes, but not with chromosome axis sites. (A) ChIP-qPCR of Set1, Swd1 and Spp1 during meiosis (t = 3 hr) comparing their association with axis and two highly transcribed genes. HA-Set1: VBD1378; Swd1-HA: VBD1399; Spp1-Myc: VBD1187. See also S2 Fig. Each dot represents a biological replicate and the bar indicates the mean. (B) Chromosome profiles of Swd1 and RNA Pol II binding in meiosis (t = 3 hr). RNA Pol II: data are from [18]. Swd1: strain VBD1399 at t = 3 hr. ChIPchip profiles are shown for chromosome 10. Centromere is indicated as a green dot, and ratios are plotted after smoothing with a 1 kb window. (C) ChIPchip signal at the indicated features. The mean signal at the 200 strongest Red1 (axis) and DSB peaks (DSB) is represented, as defined in Methods and in ref [56]. Red1 and DSB data are from previously published studies [23, 59]. Ratios after smoothing with a 2 kb window are plotted. Boxplots indicate median (line), 25th-75th percentile (box) 61.5 times the interquartile range (whiskers). Non-overlapping notches of two boxes are indicative that medians are statistically different. (D) mean mRNA levels at the 100 strongest Swd1 peaks or in the whole genome. mRNA levels are from SK1 diploid meiotic cells at t = 4 hr [58]. Error bars: S.E.M. (E) Co-immunoprecipitation by the core Set1 complex Swd1 protein from cells at 3.5 hr in meiosis analyzed by western blot using anti-HA, anti-Myc or anti-Flag antibody as indicated. Swd1-HA Spp1-Myc Mer2-Flag: VBD1401; Spp1-Myc Mer2-Flag: VBD1395. IP: immunoprecipitation. The asterisk indicates non-specific cross-hybridizing band. (F) Co-immunoprecipitation by the DSB protein Mer2 from cells at 3.5 hr in meiosis. Same antibodies as in (E). Swd1-HA Spp1-Myc Mer2-Flag: VBD1401; Spp1-Myc Swd1-HA: VBD1400. The asterisks indicate non-specific cross-hybridizing bands unrelated to Swd1-HA.

https://doi.org/10.1371/journal.pgen.1007223.g003
meiotic cells was positively correlated with that of RNA pol II and Swd1 bound only weakly to chromosome axis sites (Fig 3B and 3C, respectively). Consistent with these findings, the strongest Swd1 peaks in the genome were at highly transcribed regions (Fig 3D).

In agreement with our genome-wide localization results, Swd1 immunoprecipitated Spp1 from meiotic cells, but not Mer2, and reciprocally, Mer2 pulled down Spp1, but not Swd1 (Fig 3E and 3F, respectively). The latter result was confirmed by mass spectrometry-based analysis of Mer2-TAP interactome (S1 Table). Together, our experiments point toward Spp1 being located in two physically distinct complexes.

## A mutant that disrupts the Set1-Spp1 interaction reveals that Spp1 works independently of the Set1 complex to promote DSB formation

We next set out to determine if the presence of Spp1 in the Set1 complex is required for its DSB formation function. For this, we designed a mutation to disrupt the Set1-Spp1 interaction without affecting Spp1's PHD finger or the interaction of Spp1 with Mer2.

The domain of Spp1 that interacts with Set1 has not been determined, and we failed to identify a Spp1 mutant that would disrupt its interaction with Set1. On the Set1 side, a short region (amino acids 762-794) close to the regulatory nSet domain of Set1 is sufficient for the interaction with Spp1 in a two-hybrid test [27]. We searched for conserved amino acids in this region that could potentially be involved in protein-protein interactions [28]. We identified a conserved negatively charged acidic motif (AIKDEEDM) that we mutated into a neutral motif (ASKSSSSM) (*set1 sid* mutant for Spp1 interaction-deficient) (Fig 4A and S3 Fig). Remarkably, the mutated Set1 protein lost all interaction with Spp1 in a two-hybrid assay, while keeping its interaction with its two other known direct binding partners, Shg1 and Swd2 (Fig 4B) [29]. Furthermore, the *set1 sid* mutation induced the loss of detectable interaction of Spp1 with the Set1 complex *in vivo*, as assessed by the loss of interaction between Spp1 and the Swd1 subunit, both in vegetative and in meiotic cells (Fig 4C) and by mass spectrometry analysis of Spp1-TAP co-purified proteins from meiotic *set1 sid* cells (S4 Fig and S1 Table). Thus, no other subunit of the Set1 complex seems able to retain Spp1 on its own. In addition, in vegetative or in meiotic *set1 sid* cells, Spp1-Myc was no longer detected at the *ADH1*, *ACS1* and *CIT2* highly transcribed genes, whereas it still associated with chromosome axis in meiotic cells (Fig 4D). This *set1 sid* mutant allowed us to assess if DSB formation occurs normally even if Spp1 is not associated with the Set1 complex. Meiotic DSB frequency at two Spo11 DSB hotspots, CYS3 and *DEP1* (8), in the *set1 sid* mutant was indistinguishable from wild-type (Fig 4E, upper panel). In situations where H3K4 methylation or Spp1 is absent, such as in *set1Δ* or *spp1Δ* mutants, respectively, DSB frequencies are generally reduced, but DSB formation is increased at a few sites, including the *PES4* gene promoter [6, 18]. Importantly, in the *set1 sid* mutant, no DSBs were induced at the *PES4* promoter (Fig 4E, lower panel). These results clearly indicate that the presence of Spp1 in the Set1 complex is dispensable for DSB formation. In agreement with these results, the transient interaction of Spp1 with two DSB hotspots during meiosis was still detected in the *set1 sid* mutant (S5 Fig). This indicates a different behavior of Spp1 than that at the *ACS1* and *CIT2* genes, which are not DSB hotspots [8], where Spp1 interaction is lost in the *set1 sid* mutant (Fig 4D).

## Spp1 maintains H3K4me3 levels without interacting with the Set1 complex

Since Spp1 is believed to facilitate the Set1 catalysis of H3K4 trimethylation, we would expect that in the *set1 sid* mutant where Spp1 is no longer associated with the Set1 complex, H3K4me3 levels should decrease similarly to *spp1Δ* mutant, to about 20% of WT levels ([18] and Fig 5A and 5B and S4 Table). Surprisingly, the *set1 sid* mutant still showed high levels of
<figure><img src="image_17.png" /><figcaption>**Fig 4.** A mutant in an acidic patch of Set1 abolishes its interaction with Spp1 but does not affect meiotic DSB formation. (A) Scheme of the Set1 protein domains, with the Spp1 interacting domain (762-794) [27], RRM1 and RRM2 RNA recognizing motifs [63], the SET catalytic and nSET regulatory domains [64] and the mutations creating the *set1_sid* mutant. (B) 2-hybrid assays for the interaction between Set1 and Spp1. Growth on the -His medium indicates an interaction between the two tested proteins. (C) Interaction of Spp1 with the Set1 complex subunit Swd1 in vegetative or meiotic cells examined by Western blot as indicated. Proteins pulled down by Spp1-TAP were released in the eluate after Tev cleavage of the TAP tag. Anti-HA and anti-TAP antibodies were used. The TAP antibody still recognizes the part of the tag left after Tev cleavage. *SET1 SPP1-TAP SWD1-HA*: VBD1745; *SET1 SWD1-HA*: VBD1742; *set1_sid SPP1-TAP SWD1-HA*: VBD1836. The asterisk indicates non-specific cross-hybridizing band. (D) Binding of Spp1 detected by ChIP-qPCR in vegetative cells or in meiotic cells at t = 3 hr. No tag: ORD7339; WT: VBD1187; *set1Δ*: VBD1209; *set1_sid*: VBD1868. Each dot represents a biological replicate and the bar indicates the mean. (E) Meiotic DSB formation monitored in *dmc1Δ* cells by Southern blot at CYS3 and *DEP1* DSB (upper panel), or at the *spp1Δ*-specific *PES4* DSB (lower panel). DSB sites are indicated by an arrow. WT: ORD7354; *set1_sid*: VBD1854; *spp1Δ*: VBD1748; *set1Δ*: ORD9624. Graph shows the DSB quantification relative to the level in WT (for CYS3, *DEP1*) or *spp1Δ* cells (for *PES4* site). DSB were quantified at the 5 hr time point, with the additional 7 hr time point for *set1Δ*.</figcaption></figure>

https://doi.org/10.1371/journal.pgen.1007223.g004

H3K4me3, as assessed by Western blot and chromatin immunoprecipitation at the *ADH1* locus (Fig 5A and 5B and S4 Table). Furthermore, H3K4me3 levels were significantly reduced in the double *set1_sid spp1Δ* compared to the single *set1_sid* mutant, indicating that Spp1 still promotes H3K4me3 levels outside of the Set1 complex (Fig 5A and 5B and S4 Table). One possibility is that its binding to H3K4me3 stabilizes this mark and/or protects it from active demethylation. We asked if this function could be through its PHD finger, known to recognize H3K4 methylation in vitro [30]. A deletion of the PHD finger or the point mutation *spp1W45A* of the PHD finger mildly affects global levels of H3K4me3 ([17] and Fig 5A) and H3K4me3 levels at the *ADH1* locus (Fig 5B). However, a *set1_sid spp1W45A* double mutant showed a reduction of H3K4me3 similar to that of the *set1_sid spp1Δ* mutant (Fig 5A and 5B). These results uncover a previously unknown role of the PHD finger of Spp1 in the
<figure><img src="image_18.png" /><figcaption>**B** H3K4me3 ChIP vegetative cells</figcaption></figure>

<figure><img src="image_19.png" /></figure>

**Fig 5.** In the *set1 sid* mutant, Spp1 is still important to maintain H3K4me3 levels. (A) Histone H3K4 methylation levels in vegetatively growing cells detected by Western blot. Anti-Spp1, anti H3K4me2, anti-H3K4me3 or anti-Pgk1 antibodies were used, as indicated. A representative experiment is shown. WT: ORT4601; *set1Δ*: ORT4784; *spp1Δ*: VBH152; *set1 sid*: VBH1881; *set1 sid spp1Δ*: VBH1972; *spp1W45A*: VBH1419; *set1 sid spp1W45A*: VBH2021. The bar graph on the right indicates histone modification levels normalized to Pgk1 levels and relative to the WT strain. Values are the mean ± S.E.M. of the normalized relative levels from 3 to 6 replicates for each strain, except for *spp1W45A*, where only 2 replicates are available and the error bars indicate the range. See also S4 Table. (B) Histone H3K4me3 levels in vegetatively growing cells detected by ChIP at the highly transcribed *ADH1* gene. Same strains as in A. Values are expressed as % of input DNA, and are the mean ± S.E.M. of six independent experiments.

https://doi.org/10.1371/journal.pgen.1007223.g005

maintenance of H3K4me3 levels that takes place even without stable interaction with the Set1 complex.

## A *mer2* mutant that no longer interacts with Spp1 mimics *spp1Δ* for DSB formation

Likewise, we asked if the *spp1Δ* phenotype for DSB formation is solely due to absence of Spp1 within the Mer2 complex, and not to indirect effects on Set1 complex function. For this, we designed a mutant that would disrupt the Mer2-Spp1 interaction but preserves the Spp1-Set1 interaction. The domain of Spp1 interacting with Mer2 lies in the last 131 amino acids of Spp1, and the deletion of four amino acids C263 to C266 in this domain was proposed to abolish Mer2-Spp1 interaction *in vivo* and be important for new DSB targeting by a Gal4BD-Spp1 fusion [17]. However, this Spp1 mutant retained two-hybrid interaction with Mer2, and when inserted behind its endogenous promoter, was largely proficient in DSB formation at the tested hotspots (S6 Fig). We thus looked for altering the region of Mer2 that interacts with Spp1, with as little as possible alteration of Mer2's other functions. This region has been mapped to amino acids 165 to 232 [17], which corresponds to one of the two major coiled coils of Mer2 predicted structure (Fig 6A). Combining analyses of structure prediction and conservation of amino acids in this region, we identified several conserved amino acids predicted to be surface
exposed (Fig 6A and S7 Fig). The mutation of one of these, V195 to D, totally abolished the two-hybrid interaction of Mer2 with Spp1 (*mer2 sided* mutant, Fig 6B). When tested *in vivo*, the mutant Mer2_síd protein was still recruited to the chromosome axis like the non-mutated protein (Fig 6C). However, the Mer2_síd mutant protein was totally deficient for interaction with Spp1 in meiotic cells, in agreement with the two-hybrid results (Fig 6D). Furthermore, like in *mer2Δ* cells, Spp1 interaction with the chromosome axis was lost in the *mer2 sided* mutant (Fig 6E). We examined the meiotic phenotypes of this Spp1 interaction-defective *mer2 sided* mutant.

<figure><img src="image_20.png" /><figcaption>**Fig 6.** A *mer2* point mutant affected for its interaction with Spp1 mimics the meiotic phenotype of *spp1Δ*. (A) Scheme of the Mer2 protein, with its predicted coil-coil structures (Hhpred) and the Spp1 interacting domain (165–232) [17]. Below is the Mer2 protein sequence in the predicted coiled-coiled heptad structure. Buried positions required for coiled-coil pairing are underlined. Amino acids conserved and predicted to be surface-exposed are colored. The *mer2 sided* mutation (V195D) is indicated by an arrow. (B) 2-hybrid assays for the interaction between Mer2 and Spp1. Growth on the-His medium indicates an interaction between the two tested proteins. (C) Association of Mer2 and Mer2_síd mutant proteins with chromosome axis during meiosis (t = 3 hr) monitored by ChIP. No tag: ORD7339; Mer2-Flag: VBD1251; Mer2_síd Flag: VBD1843. Each dot represents a biological replicate and the bar indicates the mean. (D) Interaction of Spp1 with Mer2 and the Set1 complex subunit Swd1 in meiotic cells examined by Western blot. Proteins pulled down by Spp1-TAP were released in the eluate after Tev cleavage of the Tap tag. Anti-HA, anti-Flag and anti-TAP antibodies were used as indicated. MER2-FLAG SPP1-TAP SWD1-HA: VBD1745; *mer2 sided*-FLAG SPP1-TAP SWD1-HA: VBD1852; SPP1-TAP SWD1-HA: VBD1742. (E) Association of Spp1 during meiosis (t = 3 hr) monitored by ChIP. No tag: ORD7339; WT: VBD1187; *mer2Δ*: VBD1220; *mer2 sided*: VBD1924. Each dot represents a biological replicate and the bar indicates the mean. (F) Meiotic DSB formation monitored in *dmc1Δ* cells by Southern blot at CYS3 and DEP1 DSB (upper panel), or at the *spp1Δ*-specific PES4 DSB. DSB sites are indicated by arrows. WT: ORD7354; *mer2 sided*: VBD1879; *spp1Δ*: VBD1748. Graph shows the DSB quantification relative to the level in WT (for CYS3, DEP1) or *spp1Δ* cells (for PES4 site). DSB were quantified at the 5 hr time point. Each dot represents a biological replicate and the bar indicates the mean. (G) Meiotic progression as assessed by DAPI staining of strains with the indicated genotype. WT: ORD7339; *mer2 sided*: VBD1880; *spp1Δ*: VBD1769; *spo11Y135F*: VBD1291; *spo11Y135F spp1Δ*: VBD1233. Representative experiments are shown. See also S8 Fig.</figcaption></figure>
The *mer2吸入 mutant had reduced DSB formation at the two Spo11 hotspots, CYS3 and DEP1,* like *spp1Δ* (Fig 6F, upper panel). In addition, DSB formation at the *PES4* site was detected at an increased level in the *mer2吸入 mutant, similar to spp1Δ (Fig 6F, lower panel).* We conclude that the *mer2吸入 mutation recapitulates all the meiotic DSB phenotypes of spp1Δ.*

Next, we checked meiotic progression, since *spp1Δ* was reported to have a meiotic delay [18]. The *mer2吸入 mutant had wild-type meiotic progression (Fig 6G and S8 Fig).* The meiotic delay of *spp1Δ* may thus be unrelated to a recombination phenotype of *spp1Δ* but to defects related to the Set1 complex. To answer this question, we examined meiotic progression in a *spo11Y135F* DSB deficient strain. Indeed, *spp1Δ* still showed a meiotic delay in this context, which is thus DSB-independent (Fig 6G and S8 Fig).

Finally, *mer2吸入 showed a wild-type level of spore viability (97% among 100 tetrads), con- trary to mer2Δ, but similar to spp1Δ [18].* Since Mer2 is an essential DSB protein, this indicates that Mer2-sid keeps its core meiotic functions, apart from its interaction with Spp1.

Altogether, these results highlight that all meiotic DSB formation defects of *spp1Δ* can be attributed to its lack of interaction with Mer2, and are unrelated to any change in chromatin opening or accessibility that might be due to affecting the Set1 complex.

Discussion

We dissected the respective contributions of Spp1 to the Set1 and the DSB formation com-
plexes, and show that they act as independent complexes. Our data further strongly suggest
that the meiotic DSB phenotype of *spp1Δ* cells is solely due to the function of Spp1 for tether-
ing DSB sites to the chromosome axis, and not changes in chromatin structure.

Identification of a Spp1-Mer2 complex that does not comprise other stably
associated meiotic DSB proteins

Mer2 is one of the ten DSB proteins identified in S. cerevisiae, which are all connected by phys-
ical interactions, and was proposed to form a complex with two other DSB proteins, Mei4 and
Rec114 [31, 32]. However, we did not identify any other candidate than Mer2 and the Set1
complex in the meiotic Spp1 purifications. Likewise, we previously found that Spp1 does not
interact with Mei4 in co-immunoprecipitation experiments [18]. Furthermore, in a meiotic
Mer2-TAP pulldown, we retrieved Spp1 as the top candidate, but no peptide of Mei4 or
Rec114 (S1 Table). This indicates that the strength or the frequency of interaction between
Mer2 and Spp1 is higher than that of Mer2 with Mei4 and Rec114. Spp1 interacts strongly
with Mer2 in two-hybrid assays, and thus this interaction does not need meiosis-specific modi-
fication of Mer2 [17]. By contrast, the co-immunoprecipitation of Mei4 and Rec114 with Mer2
and the recruitment of Mei4 and Rec114 to the chromosome axis require Mer2 to be phos-
phorylated, upon DNA replication, by CDK and DDK [23, 33–35]. Thus, the interaction of
Mer2 with Mei4 and Rec114 is perhaps more transient than the Mer2-Spp1 interaction and
too weak to be detected in our analysis. Finally, our mer2吸入 mutant, despite losing interac-
tion with Spp1, still likely keeps interaction with Mei4 and Rec114, which are essential DSB
proteins [1]. We thus propose two separable functions of Mer2: one in recruiting the essential
DSB proteins Mei4 and Rec114, and the other in increasing the tethering of DSB sites to the
axis, through Spp1, in order to favor cleavage by Spo11. The first function is regulated by repli-
cation and Mer2 phosphorylation, whereas the latter is a “constitutive” function of Mer2.

Spp1 plays a role outside the Set1 complex for maintaining H3K4me3 levels

Histone H3K4me3 levels are slightly decreased in the *set1吸入 mutant, consistent with the*
hypothesis that the presence of Spp1 in the Set1 complex is important for normal H3K4me3
<figure><img src="image_21.png" /><figcaption>**Fig 7.** Illustration of the different functions of Spp1 for H3K4me3 and meiotic DSB formation. 1) in the Set1 complex, Spp1 has a role to allow catalysis of H3K4 trimethylation by Set1, but this function is not essential, since the *set1-sid* mutant still maintains high levels of H3K4me3. 2) In addition, Spp1 maintains H3K4me3 levels, not by stimulating Set1 catalytic activity, but likely by binding H3K4me3 with its PHD finger. This can take place without interaction with the Set1 complex. We propose this may protect H3K4me3 from active demethylation, by the Jhd2 enzyme. Other possible explanations are described in the text. 3) Finally, the simultaneous binding of Spp1 to H3K4me3 and to the axis-associated Mer2 protein is essential to promote efficient DSB formation by Spo11. It has to be noted that a PHD finger mutant of Spp1 (W45A) is still able to bind Mer2 [18], so recognition of H3K4me3 by Spp1 PHD finger is not a prerequisite for its subsequent binding to Mer2. NDR: nucleosome-depleted region; Black circle: first nucleosome of genes; H3R2: arginine 2 of histone H3, in its non-asymmetrically methylated form. The blue square represents H3K4me3.</figcaption></figure>

https://doi.org/10.1371/journal.pgen.1007223.g007

levels to stimulate Set1 catalytic activity for H3K4me3 deposition as proposed before [24]. However, it is surprising that histone H3K4me3 levels were much less affected in the *set1-sid* mutant than in the *spp1Δ* mutant, revealing an additional function of Spp1 for histone H3K4me3 levels in the absence of detectable interaction with the Set1 complex. The *in vivo* function of Spp1 recognition of H3K4me3 by its PHD finger has been poorly studied. Within the Set1 complex, it may restrict Set1 activity to the +1 nucleosomes of genes, which harbor the combination of H3K4me2/3 and H3R2 not asymmetrically methylated that is specifically recognized by Spp1 PHD finger [36]. However, this explanation does not hold for the effect we saw in the *set1_sid* mutant. We thus propose that both in wild type *SET1* cells and in the *set1_sid* mutant, the PHD finger module of Spp1 binds and protects H3K4me3 from demethylation by the Jhd2/JARID1 enzyme, which specifically demethylates H3K4me3 *in vivo* [37, 38] (Fig 7). This binding may involve a small proportion of the total Spp1 population or occur only transiently during the cell cycle, for instance in S phase, when Jhd2 was proposed to act [39]. This would explain why we failed to detect it at the tested loci in the *set1_sid* mutant. An alternative explanation might be that Spp1 directly interacts with Jhd2 and inhibits its histone H3K4 demethylation activity, in a way that requires Spp1 PHD finger. Further studies will be required to test these hypotheses and the crosstalk between Spp1 and Jhd2 in regulating H3K4me3.

The closest homolog of Spp1 in Mammals, CXXC1 (CFP1) has a N-terminal PHD finger that is able to interact with histone H3K4 methylated peptide *in vitro*, and seems important to allow Set1 complex binding to chromatin *in vitro* [40]. However, the *in vivo* function of this
PHD finger has not been determined. In addition, CXXC1 contains a motif that binds unmethylated CpG islands, allowing CXXC1 to restrict H3K4 methylation by the Set1 complex to these sites at promoters [41, 42]. It is not known if CXXC1 may associate with these regions without the associated Set1 complex. Finally, CXXC1 was reported to interact with the DNMT1 methyltransferase, through a region distinct from its Set1 interacting region [43]. However, it was not investigated if *in vivo* this interaction occurs when CXXC1 is in the Set1 complex. To conclude, so far, no Set1 complex independent function in H3K4me3 levels have been described for the homolog of Spp1 in mammals.

## Why use the same protein to regulate deposition and read H3K4me3, and not another PHD finger protein?

Although Spp1 is required both for deposition, reading of H3K4me3 and tethering to DSB proteins, we show here that surprisingly, DSB formation is as proficient when Spp1 is physically separated from the Set1 complex as in the wild-type conditions. One could thus wonder why Spp1 divided its tasks in meiosis to both participate to H3K4me3 deposition and work as a reader and tether of this mark. There are 18 other PHD finger proteins in *S. cerevisiae*, among which 8 are able to bind H3K4me2/3 *in vitro* [30]. Our findings indicate that any PHD-finger motif that could interact with Mer2 would be able to promote DSB formation as Spp1. Maybe it is just by chance that evolution selected the same protein for fulfilling its complementary functions. However, the fact that this configuration may also be conserved in mammals makes it unlikely (see below). Since we revealed here that Spp1 is able to bind H3K4me3 independently of the Set1 complex, it is possible that in wild-type cells, Spp1 “persists” on H3K4me3 after the passage of the Set1 complex, and this constitutes a first step in the tethering of DSB sites to the chromosome axis (Fig 7). There may thus be an evolutionary advantage to the same protein having three separate but sequential activities, without all occurring in the same physical complex.

## H3K4me3 function for meiotic recombination is not related to a change in chromatin accessibility but solely to anchor meiotic DSB sites to the chromosome axis

It could be argued that H3K4me3 could only be a sign of open chromatin regions, and if it had any role, it could just be to make chromatin more accessible for DSB formation. Our present data argue against this. Indeed, in our *mer2_sid* mutant, the Set1 complex is intact, and thus chromatin modifications resulting from its activity are expected to be unaltered compared to wild-type. Nevertheless, the DSB phenotype recapitulates the one seen in the absence of H3K4 methylation or Spp1, suggesting that impaired tethering of DSB sites to the chromosome axis is the only cause of the DSB reduction and redistribution and that open chromatin structure is not sufficient for normal DSB formation.

## Conservation of Spp1 functions in other organisms and unifying role of H3K4me3 or CpG islands at meiotic recombination sites

Outside of budding yeast, histone H3K4me3 is emerging as being widely used for meiotic recombination. Two main pathways for choosing sites of meiotic recombination have been identified [44]. The first one, which occurs in mice, humans, cattle and a certain number of other vertebrates involves PRDM9, which promotes both H3K4me3 and meiotic DSB formation at its binding sites [12–15]. The other one uses open chromatin at promoters, enriched in H3K4me3 and/or CpG islands, to direct meiotic recombination [6, 9–11, 45]. Although the
role of H3K4me3 in meiotic recombination has not been proven apart from in budding yeast, interesting findings showed recently that CXXC1, the homolog of Spp1 in mammals, is able to interact in two-hybrid tests both with PRDM9 and IHO1, the proposed functional homolog of Mer2 [3, 46, 47]. Like Spp1, CXXC1 is a member of the Set1 complex, and has a PHD finger that may recognize H3K4me3. In addition, it also binds unmethylated CpG at promoters. A model was proposed in which PRDM9, by interacting with CXXC1, would direct it away from promoters, and would tether the PRDM9-bound sites to the axis for Spo11 cleavage. In *Prdm9/-* mice or in organisms that do not have PRDM9, CXXC1 would bind H3K4me3 and/or CpG islands and tether these sites to chromosome axis-bound IHO1, exactly as in budding yeast. Since we have shown that in budding yeast, Spp1 is able to bind H3K4me3 outside of the Set1 complex, CXXC1 may similarly recognize H3K4me3 deposited by PRDM9 outside of the Set1 complex with which it is normally associated. This likely conservation is fascinating, and indicates that H3K4me3 and/or CpG islands could be an evolutionary conserved tether for meiotic recombination sites.

# Methods

## Yeast manipulations

All yeast strains are derivatives of the SK1 background and are listed in S2 Table. For synchronous meiosis, cells were grown in SPS presporulation medium and transferred to 1% potassium acetate with vigorous shaking at 30°C as described [48]. For strain constructions and spore viability measurements, sporulation was performed on solid sporulation medium for two days at 30°C.

## Yeast strain constructions

Yeast strains were obtained by direct transformation or crossing to get the desired genotype. All transformants were confirmed using PCR (Polymerase Chain Reaction) discriminating between correct and incorrect integrations and sequencing for epitope tag insertion or mutagenesis.

Spp1 and Mer2 were fused at their C-terminus with a TAP-tag at their endogenous loci using the pBS1539 plasmid [49]. Set1 was tagged with 6 copies of HA at its N-terminus at its endogenous locus by using plasmid pOM10 [50] and Cre-Lox excision of the marker between the *SET1* promoter and the tag. Swd1 was tagged at its C-terminus at its endogenous locus by 3 copies of HA [51]. Site directed mutagenesis was performed using PCR. The mutagenic PCR was performed on the region of interest where the gene was flanked by a selectable marker and transformed into yeast. For making the *spp1W45A* and *spp1Δ263-266* mutants without tag, we first introduced an hphMX drug resistant cassette behind the 3'UTR of *SPP1*. This construct was fully functional for meiotic DSB formation (S6 Fig). Next, we used this construct to introduce the desired mutation by PCR and transformation of the fragment containing the mutated gene and its 3'UTR hphMX cassette. For *mer2-sided-Flag* we used genomic DNA from a strain containing the *MER2-FLAG-kanMX* allele for PCR mutagenesis. For introducing the *set1_ssd* and the *mer2_ssd* mutations without an associated tag or marker, we first deleted the *SET1* or *MER2* gene with kanMX cassette by yeast transformation. We next used CRISPR-Cas9 mediated cleavage, using a plasmid encoding Cas9 and expressing a guide RNA targeted to the kanMX cassette (plasmid generously provided by G. Zhao and B. Futcher), co-transformed together with a healing *set1* or *mer2* fragment containing the desired mutation.

## Two hybrid assays

Yeast two-hybrid assays were performed exactly as described [52]. *SET1*, *SHG1*, *SPP1* and *SWD2* ORFs were PCR-amplified from SK1 genomic DNA. *MER2* cDNA sequence was
amplified from pCA5-MER2 containing *MER2* cDNA, given by Scott Keeney [31]. Briefly, PCR products were cloned in plasmids derived from the 2 hybrid vectors pGADT7 (GAL4-activating domain) and pGBKT7 (GAL4-binding domain) creating N terminal fusions and transformed in yeast haploid strains Y187 and AH109 (Clontech), respectively. Interactions were scored, after mating and diploid selection on dropout medium without leucine and tryptophan, as growth on dropout medium without leucine, tryptophan and histidine.

## Southern blot to monitor DSB formation

Cells bearing the *dmc1Δ* mutation to accumulate DSBs were harvested from meiotic time courses at each time point. Genomic DNA was prepared in low melting temperature agarose plugs and digested with the *AflII* restriction enzyme as described [48]. Southern blotting and signal quantification was performed as described [18]. Probes used were from nt 123046 to 124295 chr1 for *CYS3* and *DEP1* DSBs, and from nt 194846 to 196285, chr6 for *PES4* DSB (R64-2-1 version, Nov. 2014 *S. cerevisiae* genome coordinates).

## Tap tag purification from meiotic cells

Spp1-TAP or Mer2-TAP purification (strain VBD1266 for Spp1-TAP, VBD1877 for Spp1-TAP *set1 sid*, VBD1402 for Mer2-TAP or ORD7339 for the untagged control) was performed from 1L (2.10¹⁰ cells) of a synchronous meiotic culture at t = 3.5 hr in meiosis. Each purification was performed in parallel with a control untagged strain, ORD7339. The protocol was essentially as described in [49] with the following modifications: PMSF (phenylmethylsulfonyl fluoride) to a final concentration of 1 mM was added to the culture prior harvesting the cells. Cells were washed with TAP lysis buffer (50 mM Tris/HCl pH 7.5; 1 mM EDTA; 0.5% NP-40; 10% glycerol; 300 mM NaCl; 1 mM PMSF; 10 mM NEM; 1X Complete Mini EDTA-Free (Roche); 1X PhosSTOP (Roche)), resuspended in about 2 ml of the same buffer and frozen as noodles in liquid nitrogen. For lysis, cells were ground in a mortar in liquid nitrogen, and lysis was performed in TAP Lysis buffer plus 1 mM PMSF and 1X Complete EDTA-free protease inhibitor cocktail (Roche).

## Mass spectrometry analysis

Protein preparation for mass spectrometry-based proteomic analyses was as described [53]. Briefly, extracted proteins were stacked in the top of a SDS-PAGE gel (NuPAGE 4–12%, Invitrogen) before in-gel digestion using trypsin (Promega, sequencing grade). Resulting peptides were analysed by online nanoLC-MS/MS (UltiMate 3000 coupled to LTQ-Orbitrap Velos Pro and Ultimate 3000 RSLCnano coupled to Q-Exactive Plus, Thermo Scientific, for Spp1 and Mer2 interactomes analyses, respectively) using a 120-min gradient. Peptides and proteins were identified and quantified using MaxQuant (version 1.5.8.3 [54]) and SwissProt database (June 2017 version, *Saccharomyces cerevisiae* S288c taxonomy). Only proteins identified with a minimum of two unique + razor peptides were taken into account for further analyses. Statistical analyses were performed using ProStaR [55]. Proteins identified in the reverse and contaminant databases, proteins only identified by site and proteins exhibiting less than 3 intensity values in one condition were discarded from the list. After log₂ transformation, intensity values were normalized by median centering before missing value imputation (replacing missing values by the 2.5 percentile value of each column); statistical testing was conducted using *limma* t-test. Differentially recovered proteins were sorted out using a log₂(fold change) cut-off of 7 and a FDR threshold on remaining p-values of 1% using the Benjamini-Hochberg method.
## Co-immunoprecipitation and western blot

6 x 10⁸ cells were washed with PBS (Phosphate Buffered Saline), and lysed in 1.5 ml lysis buffer (20 mM HEPES/KOH pH7.5; 150 mM NaCl; 0.5% Triton X-100; 10% Glycerol; 1 mM MgCl2; 2 mM EDTA; 1 mM PMSF; 1X Complete Mini EDTA-Free (Roche); 1X PhosSTOP (Roche)) and 125 U/mL benzonase (Sigma) and glass beads three times for 30 s in a Fastprep instrument (MP Biomedicals). The lysate was incubated 1 hr at 4°C with 125 U/ml benzonase, and cleared by centrifugation at 13,000 g for 5 minutes. For Spp1-TAP purifications, magnetic PanMouse IgG beads (Life Technologies) were added and incubated overnight at 4°C. Beads were washed 4 times with lysis buffer and precipitated proteins were eluted by Tev cleavage (2 µg in 20 mM Tris pH8; 150 mM NaCl; 0.1% NP-40; 5% glycerol; 1 mM MgCl2; 0.5 mM EDTA; 1 mM DTT for 1 hr at room temperature). One volume of 2 x SDS protein sample buffer was added and sample were denatured 10 min at 95°C before electrophoresis. For Swd1-HA IP, 25 µl of Protein G magnetic beads (New England Biolabs) and 5 µg of mouse HA monoclonal antibody 16B12 (Covance)) were added. After overnight incubation at 4°C, beads were washed 4 times with lysis buffer and resuspended in 30 µl of 2 x SDS protein sample buffer. The beads were heated at 95 °C for 10 min. For Mer2-Flag IP, Sigma Anti-FLAG M2 magnetic beads were added to the lysate and incubated overnight at 4°C. Beads were washed twice with lysis buffer, and eluted for 2 hrs at 4°C with 12.5 µg of Flag peptide in elution buffer (20 mM Tris pH8; 150 mM NaCl; 0.1% Tween; 10% Glycerol; 5mM MgCl2; 0.5 mM EDTA). One volume of 2 x SDS protein sample buffer was added and sample were denatured 10 min at 95°C before electrophoresis. Protein eluates were loaded onto a 4–12% SDS-polyacrylamide gel and blotted to PVDF (Polyvinylidene difluoride) membrane. Antibodies used were as follows: anti-H3K4me3 (MC315, Millipore, 1/5000); anti-H3K4me2 (07–030, Millipore, 1/5000); anti-PGK1 (Invitrogen, 1/20000); anti-HA (Roche, 12CA5, 1/750); anti-Myc (Santa Cruz, 9E10, 1/500); anti-Flag (Sigma, 1/1000); anti-TAP (Invitrogen, 1/2000); anti-Spp1 (rabbit polyclonal, 1/2000). The Spp1 polyclonal antibody was made in rabbit against the full-length Spp1 fused at its N-terminus with a 6His-MBP tag. Signal was detected using the SuperSignal West Pico Chemiluminescent Substrate (ThermoFisher) and a Chemidoc touch system (Biorad). Signal was quantified using Image J software.

## ChIP-qPCR and ChIP-chip

For each meiotic time point, 2 x 10⁸ cells were processed as described [52]. For HA-Set1, we used 1 µg of monoclonal 16B12 anti-HA antibody (Covance) and 50 µL PanMouse IgG magnetic beads (Life Technologies). For H3K4me3, we used 2 µl of MC315 H3K4me3 monoclonal antibody (Millipore) and 30 µL Protein G magnetic beads (New England Biolabs). For Spp1-Myc, in ChIPchip experiment we used 0.8 µg of c-Myc monoclonal antibody (9E10, Santa Cruz) and 30 µL Protein G magnetic beads, and for ChIP qPCR, we used 1.6 µg of c-Myc monoclonal antibody (9E10, Santa Cruz) and 50 µL PanMouse IgG magnetic beads. Quantitative PCR was performed from the immunoprecipitated DNA or the whole-cell extract using a 7900HT Fast Real-Time PCR System (Applied Biosystems) and SYBR Green PCR master mix (Applied Biosystems) as described [6]. Primers for *ADH1*, *BUD23*, *GAT1*, axis (chr3, nt 232942 to 233010) and control site (*NFT1*, chr11 nt 654468 to 654549) have been described [18, 56]. Primers for *ACS1* and *CIT2* genes amplified fragments with the coordinates: chr1, nt 44778–44843 and chr3, nt 122043–122105, respectively. Results were expressed as % of DNA in the total input present in the immunoprecipitated sample and normalized to the negative control site in the middle of *NFT1*, a 3.5 kb long gene (indicated “control site” in the figures). For microarray hybridizations, whole-cell extract or immunoprecipitated DNA was amplified and labeled with either Cy3 (whole-cell extract) or Cy5 (immunoprecipitated sample) and
hybridized on an Agilent 44K yeast whole-genome oligonucleotide array as described [6]. Microarray images were read using an Axon 4000B scanner and analyzed using GenePix Pro 6.0 software (Axon Instruments, Inc.). Files were converted to text files and analyzed using the R software. The signal was normalized, denoised, smoothed and peak calling was done exactly as described before [18]. The R code has been deposited to the github depository at https://github.com/ValerieBorde/ChIPchip-analysis. In all figures, ChIP ratio is the decile-normalized ratios after denoising and smoothing using a 1 kb (Figs 2A, 2C and 3B) or 2 kb (Fig 3C) window. To each probe of the array, we attributed the mRNA level of the corresponding gene determined from vegetative cells [57] or cells at 4 hr in meiosis [58]. For probes lying in a promoter, the mRNA level of the downstream gene was attributed, and for probes in divergent promoters, the mean value of the two divergent genes was attributed. For generating the graphs in Figs 2B and 3D, the average signal of the probes at the summit of the top 100 peaks (100 values) was calculated, and compared to the signal for the entire microarray (41479 values). For generating the boxplot graphs in Fig 3C, we considered the ChIP signal at the probes corresponding to the summit of the top 200 Red1 (for axis) or DSB peaks. The definition of Red1 and DSB peaks from previously published datasets [23, 59] has been described previously [56].

## Accession numbers

The ChIP-chip data generated in this study have been deposited in the Gene Expression Omnibus database, accession number GSE102790. Processed data for all chromosomes are provided in S3 Table.

The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository [60] with the dataset identifier PXD007590.

## Supporting information

**S1 Fig. Spp1 binds chromosome axis independently of Set1.** ChIP-qPCR of Spp1 during meiosis in *SET1* (VBD1187) or *set1Δ* (VBD1209) diploids. Values represent mean ± range of two independent experiments.
(TIFF)

**S2 Fig. Set1 binds highly transcribed genes and centromere during meiosis.** ChIP-qPCR of HA-Set1 during meiosis (VBD1378 strain). Two biological replicates are shown.
(TIFF)

**S3 Fig. Multiple alignments between fungal Set1 proteins in the Spp1 binding region (residues 761 to 795 in S. cerevisiae Set1).** Above is indicated the predicted structure (H: or C:) with its confidence level.
(TIFF)

**S4 Fig. TAP purification of Spp1 meiotic complexes in wild type Set1 or in Set1 sid cells.** Silver-stained gel of TAP eluates performed at 3.5 hrs in meiosis in an untagged (ORD7339), SPP1-TAP SET1 (VBD1266) or SPP1-TAP *set1 sid* (VBD1877) strain. Untagged and Spp1-TAP samples are from the same gel as in Fig 1B.
(TIFF)

**S5 Fig. Spp1 binds two DSB hotspots independently of its association with the Set1 complex.** ChIP-qPCR of Spp1 during meiosis in *spo11Y135F* (VBD1248) or *spo11Y135F set1 sid* (VBD1944) diploids. Values represent mean ± range of two independent experiments.
(TIFF)
**S6 Fig. The spp1Δ263-266 mutant is only partially affected for interaction with Mer2 and meiotic DSB formation.** (A) Schematic structure of the Spp1 protein sequence with the position of Δ263-266 mutation.

(B) Meiotic DSB formation in *dmc1Δ* cells by Southern blot at *CYS3* and *DEP1* DSB (upper panel), or at the *spp1Δ*-specific *PES4* DSB (lower panel). WT: VBD1689; *spp1Δ263-266*: VBD1737; *spp1Δ*: VBD1748. Graph shows the DSB quantification relative to the level in WT (for *CYS3*, *DEP1*) or *spp1Δ* cells (for *PES4* site). DSB were quantified at the 5 hr time point.

(C) Two-hybrid interaction between Spp1Δ263-266 and Set1 or Mer2 proteins. Growth on-His indicates an interaction between the two tested proteins.

(TIFF)

**S7 Fig. Sequence analysis of Mer2 sequence from *Saccharomyces cerevisiae*.** (A) Coil-coiled prediction probability vs Mer2 sequence using HHblits to build the multiple sequence alignment of Mer2 and Pcoil to run coiled-coil predictions [61]. The dashed red square indicates the delimitations of the minimal Spp1 interacting domain (Spp1 ID) of Mer2 as reported [17].

(B) Multiple sequence alignment of Mer2 homologs focused on the region comprising the Spp1 interacting domain (segment 165-232 of Mer2).

(C) Sequence of the Spp1 interacting domain of Mer2, with the top sequence color code of reporting the conservation grade from gray (variable) to red (highly conserved) as calculated by Consurf server [62] and the bottom sequence reporting the predicted location of residues in the coiled-coil heptad, with positions "a" and "d" matching the most buried residues indicated by squared shaded boxes.

(TIFF)

**S8 Fig. Meiotic progression of *spp1Δ* and *mer2 sid* mutants.** Meiotic progression as assessed by DAPI staining of strains with the indicated genotype. WT: ORD7339; *mer2 sid*: VBD1880; *spp1Δ*: VBD1769; *spo11Y135F*: VBD1291; *spo11Y135F spp1Δ*: VBD1233. Error bars indicate range of two independent experiments.

(TIFF)

**S1 Table. Table of TAP tag purification results.**

(XLSX)

**S2 Table. List of strains used in this study.**

(DOCX)

**S3 Table. Processed microarray data.** Average decile-normalized ratios for each condition are included, as well as denoised ratios after smoothing with a 2 kb window. RefNumber column corresponds to reference provided by the manufacturer (Agilent) to allow easy alignment with feature characteristics or other data using the same microarray platform. The raw and processed data have been deposited to Genome Omnibus Database (GEO GSE102790).

(XLSX)

**S4 Table. Data in figures and supplemental figures.**

(XLSX)

## Acknowledgments

We thank Scott Keeney and Bruce Futcher for reagents. We thank Corinne Grey and Arnaud De Muyt for critical reading of the manuscript, and Daniel Holoch and Raphaël Margueron for discussions. We thank the recombinant proteins platform of Institut Curie for the Spp1 antibody, and the support of the discovery platform and informatics group at EDyP.
## Author Contributions

**Conceptualization:** Raphaël Guérois, Valérie Borde.

**Data curation:** Yohann Couté, Valérie Borde.

**Formal analysis:** Raphaël Guérois, Yohann Couté, Valérie Borde.

**Funding acquisition:** Jérôme Govin, Valérie Borde.

**Investigation:** Céline Adam, Raphaël Guérois, Anna Citarella, Laura Verardi, Florine Adolphe, Claire Béneut, Vérane Sommermeyer, Claire Ramus, Jérôme Govin, Yohann Couté, Valérie Borde.

**Methodology:** Raphaël Guérois, Jérôme Govin, Yohann Couté, Valérie Borde.

**Project administration:** Valérie Borde.

**Supervision:** Jérôme Govin, Yohann Couté, Valérie Borde.

**Visualization:** Céline Adam, Raphaël Guérois, Anna Citarella, Yohann Couté, Valérie Borde.

**Writing – original draft:** Valérie Borde.

**Writing – review & editing:** Céline Adam, Raphaël Guérois, Jérôme Govin, Yohann Couté, Valérie Borde.

## References

1. Cole F, Keeney S, Jasin M. Evolutionary conservation of meiotic DSB proteins: more than just Spo11. *Genes Dev.* 2010; 24(12):1201–7. https://doi.org/10.1101/gad.1944710 PMID: 20551169

2. Robert T, Vrielynck N, Mezard C, de Massy B, Grelon M. A new light on the meiotic DSB catalytic complex. *Semin Cell Dev Biol.* 2016; 54:165–76. https://doi.org/10.1016/j.semcdb.2016.02.025 PMID: 26995551

3. Stanzione M, Baumann M, Papanikos F, Dereli I, Lange J, Ramlal A, et al. Meiotic DNA break formation requires the unsynapsed chromosome axis-binding protein IHO1 (CCDC36) in mice. *Nat Cell Biol.* 2016; 18(11):1208–20. https://doi.org/10.1038/ncb3417 PMID: 27723721

4. Hunter N. Meiotic recombination: The essence of heredity. *Cold Spring Harb Perspect Biol.* 2015; 7(12).

5. Ruthenburg AJ, Allis CD, Wysocka J. Methylation of lysine 4 on histone H3: intricacy of writing and reading a single epigenetic mark. *Mol Cell.* 2007; 25(1):15–30. https://doi.org/10.1016/j.molcel.2006.12.014 PMID: 17218268

6. Borde V, Robine N, Lin W, Bonfils S, Geli V, Nicolas A. Histone H3 lysine 4 trimethylation marks meiotic recombination initiation sites. *EMBO J.* 2009; 28(2):99–111. https://doi.org/10.1038/emboj.2008.257 PMID: 19078966

7. Sollier J, Lin W, Soustelle C, Suhre K, Nicolas A, Geli V, et al. Set1 is required for meiotic S-phase onset, double-strand break formation and middle gene expression. *EMBO J.* 2004; 23(9):1957–67. https://doi.org/10.1038/sj.emboj.7600204 PMID: 15071505

8. Pan J, Sasaki M, Kniwel R, Murakami H, Blitzblau HG, Tischfield SE, et al. A hierarchical combination of factors shapes the genome-wide topography of yeast meiotic recombination initiation. *Cell.* 2011; 144(5):719–31. https://doi.org/10.1016/j.cell.2011.02.009 PMID: 21376234

9. Auton A, Rui Li Y, Kidd J, Oliveira K, Nadel J, Holloway JK, et al. Genetic recombination is targeted towards gene promoter regions in dogs. *PLoS Genet.* 2013; 9(12):e1003984. https://doi.org/10.1371/journal.pgen.1003984 PMID: 24348265

10. Choi K, Zhao X, Kelly KA, Venn O, Higgins JD, Yelina NE, et al. *Arabidopsis* meiotic crossover hot spots overlap with H2A.Z nucleosomes at gene promoters. *Nat Genet.* 2013; 45(11):1327–36. https://doi.org/10.1038/ng.2766 PMID: 24056716

11. Singhal S, Leffler EM, Sannareddy K, Turner I, Venn O, Hooper DM, et al. Stable recombination hot-spots in birds. *Science.* 2015; 350(6263):928–32. https://doi.org/10.1126/science.aad0843 PMID: 26586757

12. Baudat F, Buard J, Grey C, Fledel-Alon A, Ober C, Przeworski M, et al. PRDM9 is a major determinant of meiotic recombination hotspots in humans and mice. *Science.* 2010; 327(5967):836–40. https://doi.org/10.1126/science.1183439 PMID: 20044539
13. Myers S, Bowden R, Tumian A, Bontrop RE, Freeman C, MacFie TS, et al. Drive against hotspot motifs in primates implicates the PRDM9 gene in meiotic recombination. Science. 2010; 327(5967):876-9. https://doi.org/10.1126/science.1182363 PMID: 20044541

14. Parvanov ED, Petkov PM, Paigen K. Prdm9 controls activation of mammalian recombination hotspots. Science. 2010; 327(5967):835. https://doi.org/10.1126/science.1181495 PMID: 20044538

15. Baudat F, Imai Y, de Massy B. Meiotic recombination in mammals: localization and regulation. Nat Rev Genet. 2013; 14(11):794-806. https://doi.org/10.1038/nrg3573 PMID: 24136506

16. Brick K, Smagulova F, Khil P, Camerini-Otero RD, Petukhova GV. Genetic recombination is directed away from functional genomic elements in mice. Nature. 2012; 485(7400):642-5. https://doi.org/10.1038/nature11089 PMID: 22660327

17. Acquaviva L, Szekvolgyi L, Dichtl B, Dichtl BS, de La Roche Saint Andre C, Nicolas A, et al. The COMPASS subunit Spp1 links histone methylation to initiation of meiotic recombination. Science. 2013; 339(6116):215-8. https://doi.org/10.1126/science.1225739 PMID: 23160953

18. Sommermeyer V, Beneut C, Chaplais E, Serrentino ME, Borde V. Spp1, a member of the Set1 Complex, promotes meiotic DSB formation in promoters by tethering histone H3K4 methylation sites to chromosome axes. Mol Cell. 2013; 49(1):43-54. https://doi.org/10.1016/j.molcel.2012.11.008 PMID: 23246437

19. Zickler D, Kleckner N. Meiotic chromosomes: integrating structure and function. Annu Rev Genet. 1999; 33:603-754. PMID: 10690419

20. Blat Y, Protacio RU, Hunter N, Kleckner N. Physical and functional interactions among basic chromosomal organizational features govern early steps of meiotic chiasma formation. Cell. 2002; 111(6):791-802. PMID: 12526806

21. Borde V, de Massy B. Programmed induction of DNA double strand breaks during meiosis: setting up communication between DNA and the chromosome structure. Curr Opin Genet Dev. 2013; 23(2):147-55. https://doi.org/10.1016/j.gde.2012.12.002 PMID: 23313097

22. Kumar R, Bourbon HM, de Massy B. Functional conservation of Mei4 for meiotic DNA double-strand break formation from yeasts to mice. Genes Dev. 2010; 24(12):1266-80. https://doi.org/10.1101/gad.571710 PMID: 20551173

23. Panizza S, Mendoza MA, Berlinger M, Huang L, Nicolas A, Shirahige K, et al. Spo11-accessory proteins link double-strand break sites to the chromosome axis in early meiotic recombination. Cell. 2011; 146(3):372-83. https://doi.org/10.1016/j.cell.2011.07.003 PMID: 21816273

24. Takahashi YH, Lee JS, Swanson SK, Saraf A, Florens L, Washburn MP, et al. Regulation of H3K4 trimethylation via Cps40 (Spp1) of COMPASS is monoubiquitination independent: implication for a Phe/Tyr switch by the catalytic domain of Set1. Mol Cell Biol. 2009; 29(13):3478-86. https://doi.org/10.1128/MCB.00013-09 PMID: 19398585

25. Ng HH, Robert F, Young RA, Struhl K. Targeted recruitment of Set1 histone methylase by elongating Pol II provides a localized mark and memory of recent transcriptional activity. Mol Cell. 2003; 11(3):709-19. PMID: 12667453

26. Borde V, Goldman AS, Lichten M. Direct coupling between meiotic DNA replication and recombination initiation. Science. 2000; 290(5492):806-9. PMID: 11052944

27. Halbach A, Zhang H, Wengi A, Jablonska Z, Gruber IM, Halbeisen RE, et al. Cotranslational assembly of the yeast SET1C histone methyltransferase complex. EMBO J. 2009; 28(19):2959-70. https://doi.org/10.1038/emboj.2009.240 PMID: 19713935

28. Andreani J, Faure G, Guerois R. Versatility and invariance in the evolution of homologous heteromeric interfaces. PLoS Comput Biol. 2012; 8(8):e1002677. https://doi.org/10.1371/journal.pcbi.1002677 PMID: 22952442

29. Kim J, Kim JA, McGinty RK, Nguyen UT, Muir TW, Allis CD, et al. The n-SET domain of Set1 regulates H2B ubiquitylation-dependent H3K4 methylation. Mol Cell. 2013; 49(6):1121-33. https://doi.org/10.1016/j.molcel.2013.01.034 PMID: 23453808

30. Shi X, Kachirskai A, Walter KL, Kuo JH, Lake A, Davrazou F, et al. Proteome-wide analysis in Saccharomyces cerevisiae identifies several PHD fingers as novel direct and selective binding modules of histone H3 methylated at either lysine 4 or lysine 36. J Biol Chem. 2007; 282(4):2450-5. https://doi.org/10.1074/jbc.C600286200 PMID: 17142463

31. Arora C, Kee K, Maleki S, Keeney S. Antiviral protein Ski8 is a direct partner of Spo11 in meiotic DNA break formation, independent of its cytoplasmic role in RNA metabolism. Mol Cell. 2004; 13(4):549-59. PMID: 14992724

32. Li J, Hooker GW, Roeder GS. Saccharomyces cerevisiae Mer2, Mei4 and Rec114 form a complex required for meiotic double-strand break formation. Genetics. 2006; 173(4):1969-81. https://doi.org/10.1534/genetics.106.058768 PMID: 16783010
33. Henderson KA, Kee K, Maleki S, Santini PA, Keeney S. Cyclin-dependent kinase directly regulates initiation of meiotic recombination. Cell. 2006; 125(7):1321-32. https://doi.org/10.1016/j.cell.2006.04.039 PMID: 16814718

34. Sasanuma H, Hirota K, Fukuda T, Kakusho N, Kugou K, Kawasaki Y, et al. Cdc7-dependent phosphorylation of Mer2 facilitates initiation of yeast meiotic recombination. Genes Dev. 2008; 22(3):398-410. https://doi.org/10.1101/gad.1626608 PMID: 18245451

35. Wan L, Niu H, Futcher B, Zhang C, Shokat KM, Boulton SJ, et al. Cdc28-Clb5 (CDK-S) and Cdc7-Dbf4 (DDK) collaborate to initiate meiotic recombination in yeast. Genes Dev. 2008; 22(3):386-97. https://doi.org/10.1101/gad.1626408 PMID: 18245450

36. Kirmizis A, Santos-Rosa H, Penkett CJ, Singer MA, Vermeulen M, Mann M, et al. Arginine methylation at histone H3R2 controls deposition of H3K4 trimethylation. Nature. 2007; 449(7164):928-32. https://doi.org/10.1038/nature06160 PMID: 17898715

37. Liang G, Klose RJ, Gardner KE, Zhang Y. Yeast Jhd2p is a histone H3 Lys4 trimethyl demethylase. Nat Struct Mol Biol. 2007; 14(3):243-5. https://doi.org/10.1038/nsmb1204 PMID: 17310254

38. Seward DJ, Cubberley G, Kim S, Schonewald M, Zhang L, Tripet B, et al. Demethylation of trimethylated histone H3 Lys4 *in vivo* by JARID1 JmJC proteins. Nat Struct Mol Biol. 2007; 14(3):240-2. https://doi.org/10.1038/nsmb1200 PMID: 17310255

39. Radman-Livaja M, Liu CL, Friedman N, Schreiber SL, Rando OJ. Replication and active demethylation represent partially overlapping mechanisms for erasure of H3K4me3 in budding yeast. Plos Genet. 2010; 6(2):e1000837. https://doi.org/10.1371/journal.pgen.1000837 PMID: 20140185

40. Eberl HC, Spruijt CG, Kelstrup CD, Vermeulen M, Mann M. A map of general and specialized chromatin readers in mouse tissues generated by label-free interaction proteomics. Mol Cell. 2013; 49(2):368-78. https://doi.org/10.1016/j.molcel.2012.10.026 PMID: 23201125

41. Clouaire T, Webb S, Skene P, Illingworth R, Kerr A, Andrews R, et al. Cfp1 integrates both CpG content and gene activity for accurate H3K4me3 deposition in embryonic stem cells. Genes Dev. 2012; 26 (15):1714-28. https://doi.org/10.1101/gad.194209.112 PMID: 22855832

42. Thomson JP, Skene PJ, Selfridge J, Clouaire T, Guy J, Webb S, et al. CpG islands influence chromatin structure via the CpG-binding protein Cfp1. Nature. 2010; 464(7291):1082-6. https://doi.org/10.1038/nature08924 PMID: 20393567

43. Butler JS, Lee JH, Skalnik DG. CFP1 interacts with DNMT1 independently of association with the Setd1 Histone H3K4 methyltransferase complexes. DNA Cell Biol. 2008; 27(10):533-43. https://doi.org/10.1089/dna.2007.0714 PMID: 18680430

44. Lichten M. Molecular biology. Putting the breaks on meiosis. Science. 2015; 350(6263):913. https://doi.org/10.1126/science.aad5404 PMID: 26586748

45. Lam I, Keeney S. Nonparadoxical evolutionary stability of the recombination initiation landscape in yeast. Science. 2015; 350(6263):932-7. https://doi.org/10.1126/science.aad0814 PMID: 26586758

46. Imai Y, Baudat F, Taillepierre M, Stanzione M, Toth A, de Massy B. The PRDM9 KRAB domain is required for meiosis and involved in protein interactions. Chromosoma. 2017; 126(6):681-695. https://doi.org/10.1007/s00412-017-0631-z PMID: 28527011

47. Parvanov ED, Tian H, Billings T, Saxl RL, Spruce C, Aithal R, et al. PRDM9 interactions with other proteins provide a link between recombination hotspots and the chromosomal axis in meiosis. Mol Biol Cell. 2017; 28(3):488-99. https://doi.org/10.1091/mbc.E16-09-0686 PMID: 27932493

48. Murakami H, Borde V, Nicolas A, Keeney S. Gel electrophoresis assays for analyzing DNA double-strand breaks in *Saccharomyces cerevisiae* at various spatial resolutions. Methods Mol Biol. 2009; 557:117-42. https://doi.org/10.1007/978-1-59745-527-5_9 PMID: 19799180

49. Puig O, Caspary F, Rigaut G, Rutz B, Bouveret E, Bragado-Nilsson E, et al. The tandem affinity purification (TAP) method: a general procedure of protein complex purification. Methods. 2001; 24(3):218-29. https://doi.org/10.1006/meth.2001.1183 PMID: 11403571

50. Gauss R, Trautwein M, Sommer T, Spang A. New modules for the repeated internal and N-terminal epitope tagging of genes in *Saccharomyces cerevisiae*. Yeast. 2005; 22(1):1-12. https://doi.org/10.1002/yea.1187 PMID: 15565729

51. Longtine MS, McKenzie A 3rd, Demarini DJ, Shah NG, Wach A, Brachat A, et al. Additional modules for versatile and economical PCR-based gene deletion and modification in *Saccharomyces cerevisiae*. Yeast. 1998; 14(10):953-61. PMID: 9717241

52. Duroc Y, Kumar R, Ranjha L, Adam C, Guerois R, Md Muntaz K, et al. Concerted action of the MutLbeta heterodimer and Mer3 helicase regulates the global extent of meiotic gene conversion. Elife. 2017; 6, e21900. https://doi.org/10.7554/eLife.21900 PMID: 28051769
53. Milbradt J, Kraut A, Hutterer C, Sonntag E, Schmeiser C, Ferro M, et al. Proteomic analysis of the multi-meric nuclear egress complex of human cytomegalovirus. Mol Cell Proteomics. 2014; 13(8):2132–46. https://doi.org/10.1074/mcp.M113.035782 PMID: 24969177

54. Cox J, Mann M. MaxQuant enables high peptide identification rates, individualized p.p.b.-range mass accuracies and proteome-wide protein quantification. Nat Biotechnol. 2008; 26(12):1367–72. https://doi.org/10.1038/nbt.1511 PMID: 19029910

55. Wieczorek S, Combes F, Lazar C, Giai Gianetto Q, Gatto L, Dorffer A, et al. DAPAR & ProStaR: software to perform statistical analyses in quantitative discovery proteomics. Bioinformatics. 2017; 33(1):135–6. https://doi.org/10.1093/bioinformatics/btw580 PMID: 27605098

56. Serrentino ME, Chaplais E, Sommermeyer V, Borde V. Differential association of the conserved SUMO ligase Zip3 with meiotic double-strand break sites reveals regional variations in the outcome of meiotic recombination. PLoS Genet. 2013; 9(4):e1003416. https://doi.org/10.1371/journal.pgen.1003416 PMID: 23593021

57. Williams RM, Primig M, Washburn BK, Winzeler EA, Bellis M, Sarrauste de Menthiere C, et al. The Ume6 regulon coordinates metabolic and meiotic gene expression in yeast. Proc Natl Acad Sci U S A. 2002; 99(21):13431–6. https://doi.org/10.1073/pnas.202495299 PMID: 12370439

58. Primig M, Williams RM, Winzeler EA, Tevzadze GG, Conway AR, Hwang SY, et al. The core meiotic transcriptome in budding yeasts. Nat Genet. 2000; 26(4):415–23. https://doi.org/10.1038/82539 PMID: 11101837

59. Buhler C, Borde V, Lichten M. Mapping meiotic single-strand DNA reveals a new landscape of DNA double-strand breaks in *Saccharomyces cerevisiae*. PLoS Biol. 2007; 5(12):e324. https://doi.org/10.1371/journal.pbio.0060104 PMID: 18076285

60. Vizcaino JA, Csordas A, Del-Toro N, Dianes JA, Griss J, Lavidas I, et al. 2016 update of the PRIDE database and its related tools. Nucleic Acids Res. 2016; 44(22):11033. https://doi.org/10.1093/nar/gkw880 PMID: 27683222

61. Alva V, Nam SZ, Soding J, Lupas AN. The MPI bioinformatics Toolkit as an integrative platform for advanced protein sequence and structure analysis. Nucleic Acids Res. 2016; 44(W1):W410–5. https://doi.org/10.1093/nar/gkw348 PMID: 27131380

62. Ashkenazy H, Abadi S, Martz E, Chay O, Mayrose I, Pupko T, et al. ConSurf 2016: an improved methodology to estimate and visualize evolutionary conservation in macromolecules. Nucleic Acids Res. 2016; 44(W1):W344–50 https://doi.org/10.1093/nar/gkw408 PMID: 27166375

63. Tresaugues L, Dehe PM, Guerois R, Rodriguez-Gil A, Varlet I, Salah P, et al. Structural characterization of Set1 RNA recognition motifs and their role in histone H3 lysine 4 methylation. J Mol Biol. 2006; 359(5):1170–81. https://doi.org/10.1016/j.jmb.2006.04.050 PMID: 16787775

64. Roguev A, Schaft D, Shevchenko A, Pijnappel WW, Wilm M, Aasland R, et al. The *Saccharomyces cerevisiae* Set1 complex includes an Ash2 homologue and methylates histone 3 lysine 4. EMBO J. 2001; 20(24):7137–48. https://doi.org/10.1093/emboj/20.24.7137 PMID: 11742990