Meiotic recombination involves the exchange of DNA between two homologous chromosomes, forming cross-overs and gene conversion events. The cross-over process is important for the proper segregation of chromosomes during meiosis, and drives genetic diversity. Human hotspots are enriched for a 13-bp motif, CCNCCNTNNCCNC; a close match to this motif occurs in about 40% of our cross-over hotspots. A DNA binding protein called PRDM9, having histone trimethyltransferase (H3K4me3) activity, binds the motif and is becoming established as a major determinant of recombination hotspots (narrow regions with high cross-over activity). This research aimed to understand the mechanisms involved in promoting PRDM9 binding to its target sites, and subsequently, initiating cross-over hotspot activity. We first explored the relationship between PRDM9 binding and DNA sequence, to directly confirm whether PRDM9 binds to the 13-bp hotspot motif using in-vitro gel-shift assays, and found that it does bind sequence specifically to the canonical 13-mer motif. PRDM9 is able to bind the motif in a highly selective manner, with certain single base pair changes abolishing binding. However, we observe that it is also able to tolerate degeneracy in its binding sites, as demonstrated by strong in-vitro binding to degenerate versions of the 13-bp motif. Hence, these results confirmed that PRDM9 is able to directly bind to the 13-bp hotspot motifs, and given that it can also tolerate degeneracy, this raised the question of why PRDM9 is able to bind only a subset of all such potential binding sites in the genome. To address this, a ChIP-seq analysis was performed to identify genome wide binding sites for PRDM9. This information also helped us to characterise binding sites and investigate if factors such as the local chromatin environment play a role in specifying PRDM9 binding tar- gets and hotspot formation. We were able to identify over 170,000 PRDM9 binding sites in the genome. Surprisingly, these binding sites were also enriched in promoter regions, however, bound sites in these regulatory regions showed low recombination activity. We found that PRDM9 is able to confer the H3K4me3 mark on all bound sites, even those without a pre-existing H3K4me2 mark. We also investigated the role of other chromatin related marks on PRDM9 binding and found that binding occurs in chromatin accessible, but nucleosome rich regions, whereas heterochromatin regions tend to inhibit binding. Further, for hotspot formation, it was seen that less chromatin accessible, nucleosome dense regions away from transcribed sites, are preferred. Hotspots tend to avoid regions marked by transcription activating histone modifications, however, these regions do not appear to inhibit PRDM9 binding itself. These results show how PRDM9 binding in the genome is dependent on both primary DNA sequence and the surrounding epigenetic factors. Together these factors promote binding and, with additional downstream factors, positioning of hotspot locations in the human genome.

Previous evidence has shown that the simple sequences microsatellites and poly-purine/poly-pyrimidine tracts (PPTs) could be both a cause, and an effect, of meiotic recombination. The causal link between simple sequences and recombination has not been much explored, however, probably because other evidence has cast doubt on its generality, though this evidence has never been conclusive. Several questions have remained unanswered in the literature, and I have addressed aspects of three of them in my thesis. First, what is the scale and magnitude of the association between simple sequences and recombination? I found that microsatellites and PPTs are strongly associated with meiotic double-strand break (DSB) hotspots in yeast, and that PPTs are generally more common in human recombination hotspots, particularly in close proximity to hotspot central regions, in which recombination events are markedly more frequent. I also showed that these associations can't be explained by coincidental mutual associations between simple sequences, recombination and other factors previously shown to correlate with both. A second question not conclusively answered in the literature is whether simple sequences, or their high levels of polymorphism, are an effect of recombination. I used three methods to address this question. Firstly, I investigated the distributions of two-copy tandem repeats and short PPTs in relation to yeast DSB hotspots in order to look for evidence of an involvement of recombination in simple sequence formation. I found no significant associations. Secondly, I compared the fraction of simple sequences containing polymorphic sites between human recombination hotspots and coldspots. The third method I used was generalized linear model analysis, with which I investigated the correlation between simple sequence variation and recombination rate, and the influence on the correlation of additional factors with potential relevance including GC-content and gene density. Both the direct comparison and correlation methods showed a very weak and inconsistent effect of recombination on simple sequence polymorphism in the human genome.Whether simple sequences are an important cause of recombination events is a third question that has received relatively little previous attention, and I have explored one aspect of it. Simple sequences of the types I studied have previously been shown to form non-B-DNA structures, which can be recombinagenic in model systems. Using a previously described sodium bisulphite modification assay, I tested for the presence of these structures in sequences amplified from the central regions of hotspots and cloned into supercoiled plasmids. I found significantly higher sensitivity to sodium bisulphite in humans in than in chimpanzees in three out of six genomic regions in which there is a hotspot in humans but none in chimpanzees. In the DNA2 hotspot, this correlated with a clear difference in numbers of molecules showing long contiguous strings of converted cytosines, which are present in previously described intramolecular quadruplex and triplex structures. Two out of the five other hotspots tested show evidence for secondary structure comparable to a known intramolecular triplex, though with similar patterns in humans and chimpanzees. In conclusion, my results clearly motivate further investigation of a functional link between simple sequences and meiotic recombination, including the putative role of non-B-DNA structures.

Previous evidence has shown that the simple sequences microsatellites and poly-purine/poly-pyrimidine tracts (PPTs) could be both a cause, and an effect, of meiotic recombination. The causal link between simple sequences and recombination has not been much explored, however, probably because other evidence has cast doubt on its generality, though this evidence has never been conclusive. Several questions have remained unanswered in the literature, and I have addressed aspects of three of them in my thesis. First, what is the scale and magnitude of the association between simple sequences and recombination? I found that microsatellites and PPTs are strongly associated with meiotic double-strand break (DSB) hotspots in yeast, and that PPTs are generally more common in human recombination hotspots, particularly in close proximity to hotspot central regions, in which recombination events are markedly more frequent. I also showed that these associations can't be explained by coincidental mutual associations between simple sequences, recombination and other factors previously shown to correlate with both. A second question not conclusively answered in the literature is whether simple sequences, or their high levels of polymorphism, are an effect of recombination. I used three methods to address this question. Firstly, I investigated the distributions of two-copy tandem repeats and short PPTs in relation to yeast DSB hotspots in order to look for evidence of an involvement of recombination in simple sequence formation. I found no significant associations. Secondly, I compared the fraction of simple sequences containing polymorphic sites between human recombination hotspots and coldspots. The third method I used was generalized linear model analysis, with which I investigated the correlation between simple sequence variation and recombination rate, and the influence on the correlation of additional factors with potential relevance including GC-content and gene density. Both the direct comparison and correlation methods showed a very weak and inconsistent effect of recombination on simple sequence polymorphism in the human genome.Whether simple sequences are an important cause of recombination events is a third question that has received relatively little previous attention, and I have explored one aspect of it. Simple sequences of the types I studied have previously been shown to form non-B-DNA structures, which can be recombinagenic in model systems. Using a previously described sodium bisulphite modification assay, I tested for the presence of these structures in sequences amplified from the central regions of hotspots and cloned into supercoiled plasmids. I found significantly higher sensitivity to sodium bisulphite in humans in than in chimpanzees in three out of six genomic regions in which there is a hotspot in humans but none in chimpanzees. In the DNA2 hotspot, this correlated with a clear difference in numbers of molecules showing long contiguous strings of converted cytosines, which are present in previously described intramolecular quadruplex and triplex structures. Two out of the five other hotspots tested show evidence for secondary structure comparable to a known intramolecular triplex, though with similar patterns in humans and chimpanzees. In conclusion, my results clearly motivate further investigation of a functional link between simple sequences and meiotic recombination, including the putative role of non-B-DNA structures.

Meiotic recombination is a fundamental biological process in sexually reproducing organisms, enabling offspring to inherit novel combinations of mutations, and ensuring even segregation of chromosomes into gametes. Recombination is initiated by programmed Double Strand Breaks (DSBs), the genomic locations of which are determined in most mammals by PRDM9, a rapidly evolving DNA-binding protein. In crosses between different mouse subspecies, certain Prdm9 alleles cause infertility in hybrid males, implying a critical role in fertility and speciation. Upon binding to DNA, PRDM9 deposits a histone modification (H3K4me3) typically found in the promoters of expressed genes, suggesting that binding might alter the expression of nearby genes. Many other questions have remained about how PRDM9 initiates recombination, how it causes speciation, and why it evolves so rapidly. This body of work investigates these questions using complementary experimental and analytical methodologies. By generating a map of human PRDM9 binding sites and applying novel sequence analysis methods, I uncovered new DNA-binding modalities of PRDM9 and identified sequence-independent factors that predict binding and recombination outcomes. I also confirmed that PRDM9 can affect gene expression by binding to promoters, identifying candidate regulatory targets in meiosis. Furthermore, I showed that PRDM9’s DNA-binding domain also mediates strong protein-protein interactions that produce PRDM9 multimers, which may play an important functional role. Finally, by generating high-resolution maps of PRDM9 binding in hybrid mice, I provide evidence for a mechanism to explain PRDM9-mediated speciation as a consequence of the joint evolution of PRDM9 and its binding targets. This work reveals that PRDM9 binding on one chromosome strongly impacts DSB formation and/or repair on the homologue, suggesting a novel role for PRDM9 in promoting efficient homology search and DSB repair, both critical for meiotic progression and fertility. One consequence is that PRDM9 may play a wider role in mammalian speciation.

Les positions des évènements de recombinaison s’agrègent ensemble, formant des hotspots déterminés en partie par la protéine à évolution rapide PRDM9. En particulier, ces positions de hotspots sont déterminées par le domaine de doigts de zinc (ZnF) de PRDM9 qui reconnait certains motifs d’ADN. Les allèles de PRDM9 contenant le ZnF de type k ont été préalablement associés avec une cohorte de patients affectés par la leucémie aigüe lymphoblastique. Les allèles de PRDM9 sont difficiles à identifier à partir de données de séquençage de nouvelle génération (NGS), en raison de leur nature répétitive. Dans ce projet, nous proposons une méthode permettant la caractérisation d’allèles de PRDM9 à partir de données de NGS, qui identifie le nombre d’allèles contenant un type spécifique de ZnF. Cette méthode est basée sur la corrélation entre les profils représentant le nombre de séquences nucléotidiques uniques à chaque ZnF retrouvés chez les lectures de NGS simulées sans erreur d’une paire d’allèles et chez les lectures d’un échantillon. La validité des prédictions obtenues par notre méthode est confirmée grâce à analyse basée sur les simulations. Nous confirmons également que la méthode peut correctement identifier le génotype d’allèles de PRDM9 qui n’ont pas encore été identifiés. Nous conduisons une analyse préliminaire identifiant le génotype des allèles de PRDM9 contenant un certain type de ZnF dans une cohorte de patients atteints de glioblastomes multiforme pédiatrique, un cancer du cerveau caractérisé par les mutations récurrentes dans le gène codant pour l’histone H3, la cible de l’activité épigénétique de PRDM9. Cette méthode ouvre la possibilité d’identifier des associations entre certains allèles de PRDM9 et d’autres types de cancers pédiatriques, via l’utilisation de bases de données de NGS de cellules tumorales. / The positions of recombination events cluster tightly together in recombination hotspots, which are determined in part by the rapidly evolving protein PRDM9 via its tri- methyltransferase activity. The locations of hotspots are determined by the repetitive ZnF array of PRDM9, which binds to DNA. Alleles of PRDM9 containing the k-ZnF have previously been associated with patients affected with childhood acute lymphoblastic leukaemia. PRDM9 alleles are notoriously difficult to type due to the repetitive nature of the ZnF arrays. Here, we propose a method to characterize the alleles of PRDM9 from next- generation sequencing samples, by identifying the number of alleles containing a specific ZnF type. Our method is based on the correlation between profiles from the sample, representing the counts of nucleotide sequences unique to each ZnF, and from ideal sets of short reads representing an allele pair. We conduct a simulation analysis to examine the validity of the predictions obtained by our method with all pairs of known alleles. We confirm that the method can accurately genotype previously unobserved PRDM9 alleles. We also conducted a preliminary analysis to identify the PRDM9 k-ZnF genotype in a cohort of paediatric glioblastoma (pGBM), a childhood cancer characterized by the recurrent mutations in the coding sequence of the histone H3, the target of the enzymatic activity of PRDM9. Although no associations of k-ZnF containing PRDM9 alleles is found in our pGBM cohort, this method opens the possibility of identifying associations between certain PRDM9 alleles with other types of early onset childhood cancers, through a data-mining effort in public cancer databases.