PRDM9 Drives Evolutionary Erosion of Hotspots in Mus musculus through Haplotype-Specific Initiation of Meiotic Recombination

Meiotic recombination generates new genetic variation and assures the proper segregation of chromosomes in gametes. PRDM9, a zinc finger protein with histone methyltransferase activity, initiates meiotic recombination by binding DNA at recombination hotspots and directing the position of DNA double-strand breaks (DSB). The DSB repair mechanism suggests that hotspots should eventually self-destruct, yet genome-wide recombination levels remain constant, a conundrum known as the hotspot paradox. To test if PRDM9 drives this evolutionary erosion, we measured activity of the Prdm9 ^Cst allele in two Mus musculus subspecies, M.m. castaneus, in which Prdm9^Cst arose, and M.m. domesticus, into which Prdm9^Cst was introduced experimentally. Comparing these two strains, we find that haplotype differences at hotspots lead to qualitative and quantitative changes in PRDM9 binding and activity. Using Mus spretus as an outlier, we found most variants affecting PRDM9^Cst binding arose and were fixed in M.m. castaneus, suppressing hotspot activity. Furthermore, M.m. castaneus×M.m. domesticus F1 hybrids exhibit novel hotspots, with large haplotype biases in both PRDM9 binding and chromatin modification. These novel hotspots represent sites of historic evolutionary erosion that become activated in hybrids due to crosstalk between one parent''s Prdm9 allele and the opposite parent''s chromosome. Together these data support a model where haplotype-specific PRDM9 binding directs biased gene conversion at hotspots, ultimately leading to hotspot erosion.     

DNA binding specificities of the long zinc-finger recombination protein PRDM9

Background

Meiotic recombination ensures proper segregation of homologous chromosomes and creates genetic variation. In many organisms, recombination occurs at limited sites, termed ''hotspots'', whose positions in mammals are determined by PR domain member 9 (PRDM9), a long-array zinc-finger and chromatin-modifier protein. Determining the rules governing the DNA binding of PRDM9 is a major issue in understanding how it functions.

Results

Mouse PRDM9 protein variants bind to hotspot DNA sequences in a manner that is specific for both PRDM9 and DNA haplotypes, and that in vitro binding parallels its in vivo biological activity. Examining four hotspots, three activated by Prdm9^Cst and one activated by Prdm9^Dom2, we found that all binding sites required the full array of 11 or 12 contiguous fingers, depending on the allele, and that there was little sequence similarity between the binding sites of the three Prdm9^Cst activated hotspots. The binding specificity of each position in the Hlx1 binding site, activated by Prdm9^Cst, was tested by mutating each nucleotide to its three alternatives. The 31 positions along the binding site varied considerably in the ability of alternative bases to support binding, which also implicates a role for additional binding to the DNA phosphate backbone.

Conclusions

These results, which provide the first detailed mapping of PRDM9 binding to DNA and, to our knowledge, the most detailed analysis yet of DNA binding by a long zinc-finger array, make clear that the binding specificities of PRDM9, and possibly other long-array zinc-finger proteins, are unusually complex.     

Mouse PRDM9 DNA-binding specificity determines sites of histone H3 lysine 4 trimethylation for initiation of meiotic recombination

Meiotic recombination generates reciprocal exchanges between homologous chromosomes (also called crossovers, COs) that are essential for proper chromosome segregation during meiosis and are a major source of genome diversity by generating new allele combinations. COs have two striking properties: they occur at specific sites, called hotspots, and these sites evolve rapidly. In mammals, the Prdm9 gene, which encodes a meiosis-specific histone H3 methyltransferase, has recently been identified as a determinant of CO hotspots. Here, using transgenic mice, we show that the sole modification of PRDM9 zinc fingers leads to changes in hotspot activity, histone H3 lysine 4 trimethylation (H3K4me3) levels, and chromosome-wide distribution of COs. We further demonstrate by an in vitro assay that the PRDM9 variant associated with hotspot activity binds specifically to DNA sequences located at the center of the three hotspots tested. Remarkably, we show that mutations in cis located at hotspot centers and associated with a decrease of hotspot activity affect PRDM9 binding. Taken together, these results provide the direct demonstration that Prdm9 is a master regulator of hotspot localization through the DNA binding specificity of its zinc finger array and that binding of PRDM9 at hotspots promotes local H3K4me3 enrichment.     

Prdm9 Polymorphism Unveils Mouse Evolutionary Tracks

PR/SET domain containing 9 (Prdm9) mediates histone modifications such as H3K4me3 and marks hotspots of meiotic recombination. In many mammalian species, the Prdm9 gene is highly polymorphic. Prdm9 polymorphism is assumed to play two critical roles in evolution: to diversify the spectrum of meiotic recombination hotspots and to cause male hybrid sterility, leading to reproductive isolation and speciation. Nevertheless, information about Prdm9 sequences in natural populations is very limited. In this study, we conducted a comprehensive population survey on Prdm9 polymorphism in the house mouse, Mus musculus. Overall M. musculus Prdm9 displays an extraordinarily high level of polymorphism, particularly in regions encoding zinc finger repeats, which recognize recombination hotspots. Prdm9 alleles specific to various M. musculus subspecies dominate in subspecies territories. Moreover, introgression into other subspecies territories was found for highly divergent Prdm9 alleles associated with t-haplotype. The results of our phylogeographical analysis suggest that the requirement for hotspot diversity depends on geographical range and time span in mouse evolution, and that Prdm9 polymorphism has not been maintained by a simple balanced selection in the population of each subspecies.     

Diversity of Prdm9 Zinc Finger Array in Wild Mice Unravels New Facets of the Evolutionary Turnover of this Coding Minisatellite

In humans and mice, meiotic recombination events cluster into narrow hotspots whose genomic positions are defined by the PRDM9 protein via its DNA binding domain constituted of an array of zinc fingers (ZnFs). High polymorphism and rapid divergence of the Prdm9 gene ZnF domain appear to involve positive selection at DNA-recognition amino-acid positions, but the nature of the underlying evolutionary pressures remains a puzzle. Here we explore the variability of the Prdm9 ZnF array in wild mice, and uncovered a high allelic diversity of both ZnF copy number and identity with the caracterization of 113 alleles. We analyze features of the diversity of ZnF identity which is mostly due to non-synonymous changes at codons −1, 3 and 6 of each ZnF, corresponding to amino-acids involved in DNA binding. Using methods adapted to the minisatellite structure of the ZnF array, we infer a phylogenetic tree of these alleles. We find the sister species Mus spicilegus and M. macedonicus as well as the three house mouse (Mus musculus) subspecies to be polyphyletic. However some sublineages have expanded independently in Mus musculus musculus and M. m. domesticus, the latter further showing phylogeographic substructure. Compared to random genomic regions and non-coding minisatellites, none of these patterns appears exceptional. In silico prediction of DNA binding sites for each allele, overlap of their alignments to the genome and relative coverage of the different families of interspersed repeated elements suggest a large diversity between PRDM9 variants with a potential for highly divergent distributions of recombination events in the genome with little correlation to evolutionary distance. By compiling PRDM9 ZnF protein sequences in Primates, Muridae and Equids, we find different diversity patterns among the three amino-acids most critical for the DNA-recognition function, suggesting different diversification timescales.     

Prdm9 Intersubspecific Interactions in Hybrid Male Sterility of House Mouse

The classical definition posits hybrid sterility as a phenomenon when two parental taxa each of which is fertile produce a hybrid that is sterile. The first hybrid sterility gene in vertebrates, Prdm9, coding for a histone methyltransferase, was identified in crosses between two laboratory mouse strains derived from Mus mus musculus and M. m. domesticus subspecies. The unique function of PRDM9 protein in the initiation of meiotic recombination led to the discovery of the basic molecular mechanism of hybrid sterility in laboratory crosses. However, the role of this protein as a component of reproductive barrier outside the laboratory model remained unclear. Here, we show that the Prdm9 allelic incompatibilities represent the primary cause of reduced fertility in intersubspecific hybrids between M. m. musculus and M. m. domesticus including 16 musculus and domesticus wild-derived strains. Disruption of fertility phenotypes correlated with the rate of failure of synapsis between homologous chromosomes in meiosis I and with early meiotic arrest. All phenotypes were restored to normal when the domesticus Prdm9^dom2 allele was substituted with the Prdm9^dom2H humanized variant. To conclude, our data show for the first time the male infertility of wild-derived musculus and domesticus subspecies F1 hybrids controlled by Prdm9 as the major hybrid sterility gene. The impairment of fertility surrogates, testes weight and sperm count, correlated with increasing difficulties of meiotic synapsis of homologous chromosomes and with meiotic arrest, which we suppose reflect the increasing asymmetry of PRDM9-dependent DNA double-strand breaks.     

What are the genomic drivers of the rapid evolution of PRDM9?

Mammalian Prdm9 has been proposed to be a key determinant of the positioning of chromosome double-strand breaks during meiosis, a contributor to speciation processes, and the most rapidly evolving gene in human, and other animal, genomes. Prdm9 genes often exhibit substantial variation in their numbers of encoded zinc fingers (ZFs), not only between closely related species but also among individuals of a species. The near-identity of these ZF sequences appears to render them very unstable in copy number. The rare sequence differences, however, cluster within ZF sites that determine the DNA-binding specificity of PRDM9, and these substitutions are frequently positively selected. Here, possible drivers of the rapid evolution of Prdm9 are discussed, including selection for efficient pairing of homologous chromosomes or for recombination of deleterious linked alleles, and selection against depletion of recombination hotspots or against disease-associated genome rearrangement.     

The Red Queen Model of Recombination Hotspots Evolution in the Light of Archaic and Modern Human Genomes

Recombination is an essential process in eukaryotes, which increases diversity by disrupting genetic linkage between loci and ensures the proper segregation of chromosomes during meiosis. In the human genome, recombination events are clustered in hotspots, whose location is determined by the PRDM9 protein. There is evidence that the location of hotspots evolves rapidly, as a consequence of changes in PRDM9 DNA-binding domain. However, the reasons for these changes and the rate at which they occur are not known. In this study, we investigated the evolution of human hotspot loci and of PRDM9 target motifs, both in modern and archaic human lineages (Denisovan) to quantify the dynamic of hotspot turnover during the recent period of human evolution. We show that present-day human hotspots are young: they have been active only during the last 10% of the time since the divergence from chimpanzee, starting to be operating shortly before the split between Denisovans and modern humans. Surprisingly, however, our analyses indicate that Denisovan recombination hotspots did not overlap with modern human ones, despite sharing similar PRDM9 target motifs. We further show that high-affinity PRDM9 target motifs are subject to a strong self-destructive drive, known as biased gene conversion (BGC), which should lead to the loss of the majority of them in the next 3 MYR. This depletion of PRDM9 genomic targets is expected to decrease fitness, and thereby to favor new PRDM9 alleles binding different motifs. Our refined estimates of the age and life expectancy of human hotspots provide empirical evidence in support of the Red Queen hypothesis of recombination hotspots evolution.     

High diversity at PRDM9 in chimpanzees and bonobos

Background

The PRDM9 locus in mammals has increasingly attracted research attention due to its role in mediating chromosomal recombination and possible involvement in hybrid sterility and hence speciation processes. The aim of this study was to characterize sequence variation at the PRDM9 locus in a sample of our closest living relatives, the chimpanzees and bonobos.

Methodology/Principal Findings

PRDM9 contains a highly variable and repetitive zinc finger array. We amplified this domain using long-range PCR and determined the DNA sequences using conventional Sanger sequencing. From 17 chimpanzees representing three subspecies and five bonobos we obtained a total of 12 alleles differing at the nucleotide level. Based on a data set consisting of our data and recently published Pan PRDM9 sequences, we found that at the subspecies level, diversity levels did not differ among chimpanzee subspecies or between chimpanzee subspecies and bonobos. In contrast, the sample of chimpanzees harbors significantly more diversity at PRDM9 than samples of humans. Pan PRDM9 shows signs of rapid evolution including no alleles or ZnFs in common with humans as well as signals of positive selection in the residues responsible for DNA binding.

Conclusions and Significance

The high number of alleles specific to the genus Pan, signs of positive selection in the DNA binding residues, and reported lack of conservation of recombination hotspots between chimpanzees and humans suggest that PRDM9 could be active in hotspot recruitment in the genus Pan. Chimpanzees and bonobos are considered separate species and do not have overlapping ranges in the wild, making the presence of shared alleles at the amino acid level between the chimpanzee and bonobo species interesting in view of the hypothesis that PRDM9 plays a universal role in interspecific hybrid sterility.     

A Highlights from MBoC Selection: EWSR1 affects PRDM9-dependent histone 3 methylation and provides a link between recombination hotspots and the chromosome axis protein REC8

Meiotic recombination in most mammals requires recombination hotspot activation through the action of the histone 3 Lys-4 and Lys-36 methyltransferase PRDM9 to ensure successful double-strand-break initiation and repair. Here we show that EWSR1, a protein whose role in meiosis was not previously clarified in detail, binds to both PRDM9 and pREC8, a phosphorylated meiosis-specific cohesin, in male meiotic cells. We created a Ewsr1 conditional knockout mouse model to deplete EWSR1 before the onset of meiosis and found that absence of EWSR1 causes meiotic arrest with decreased histone trimethylation at meiotic hotspots, impaired DNA double-strand-break repair, and reduced crossover number. Our results demonstrate that EWSR1 is essential for promoting PRDM9-dependent histone methylation and normal meiotic progress, possibly by facilitating the linking between PRDM9-bound hotspots and the nascent chromosome axis through its component cohesin pREC8.     

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

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

Trans-Regulation of Mouse Meiotic Recombination Hotspots by Rcr1

Meiotic recombination is required for the orderly segregation of chromosomes during meiosis and for providing genetic diversity among offspring. Among mammals, as well as yeast and higher plants, recombination preferentially occurs at highly delimited chromosomal sites 1–2 kb long known as hotspots. Although considerable progress has been made in understanding the roles various proteins play in carrying out the molecular events of the recombination process, relatively little is understood about the factors controlling the location and relative activity of mammalian recombination hotspots. To search for trans-acting factors controlling the positioning of recombination events, we compared the locations of crossovers arising in an 8-Mb segment of a 100-Mb region of mouse Chromosome 1 (Chr 1) when the longer region was heterozygous C57BL/6J (B6) × CAST/EiJ (CAST) and the remainder of the genome was either similarly heterozygous or entirely homozygous B6. The lack of CAST alleles in the remainder of the genome resulted in profound changes in hotspot activity in both females and males. Recombination activity was lost at several hotspots; new, previously undetected hotspots appeared; and still other hotspots remained unaffected, indicating the presence of distant trans-acting gene(s) whose CAST allele(s) activate or suppress the activity of specific hotspots. Testing the activity of three activated hotspots in sperm samples from individual male progeny of two genetic crosses, we identified a single trans-acting regulator of hotspot activity, designated Rcr1, that is located in a 5.30-Mb interval (11.74–17.04 Mb) on Chr 17. Using an Escherichia coli cloning assay to characterize the molecular products of recombination at two of these hotspots, we found that Rcr1 controls the appearance of both crossover and noncrossover gene conversion events, indicating that it likely controls the sites of the double-strand DNA breaks that initiate the recombination process.     

Transmission Distortion Affecting Human Noncrossover but Not Crossover Recombination: A Hidden Source of Meiotic Drive

Meiotic recombination ensures the correct segregation of homologous chromosomes during gamete formation and contributes to DNA diversity through both large-scale reciprocal crossovers and very localised gene conversion events, also known as noncrossovers. Considerable progress has been made in understanding factors such as PRDM9 and SNP variants that influence the initiation of recombination at human hotspots but very little is known about factors acting downstream. To address this, we simultaneously analysed both types of recombinant molecule in sperm DNA at six highly active hotspots, and looked for disparity in the transmission of allelic variants indicative of any cis-acting influences. At two of the hotspots we identified a novel form of biased transmission that was exclusive to the noncrossover class of recombinant, and which presumably arises through differences between crossovers and noncrossovers in heteroduplex formation and biased mismatch repair. This form of biased gene conversion is not predicted to influence hotspot activity as previously noted for SNPs that affect recombination initiation, but does constitute a powerful and previously undetected source of recombination-driven meiotic drive that by extrapolation may affect thousands of recombination hotspots throughout the human genome. Intriguingly, at both of the hotspots described here, this drive favours strong (G/C) over weak (A/T) base pairs as might be predicted from the well-established correlations between high GC content and recombination activity in mammalian genomes.     

Prdm9, a major determinant of meiotic recombination hotspots, is not functional in dogs and their wild relatives, wolves and coyotes

Meiotic recombination is a fundamental process needed for the correct segregation of chromosomes during meiosis in sexually reproducing organisms. In humans, 80% of crossovers are estimated to occur at specific areas of the genome called recombination hotspots. Recently, a protein called PRDM9 was identified as a major player in determining the location of genome-wide meiotic recombination hotspots in humans and mice. The origin of this protein seems to be ancient in evolutionary time, as reflected by its fairly conserved structure in lineages that diverged over 700 million years ago. Despite its important role, there are many animal groups in which Prdm9 is absent (e.g. birds, reptiles, amphibians, diptera) and it has been suggested to have disruptive mutations and thus to be a pseudogene in dogs. Because of the dog's history through domestication and artificial selection, we wanted to confirm the presence of a disrupted Prdm9 gene in dogs and determine whether this was exclusive of this species or whether it also occurred in its wild ancestor, the wolf, and in a close relative, the coyote. We sequenced the region in the dog genome that aligned to the last exon of the human Prdm9, containing the entire zinc finger domain, in 4 dogs, 17 wolves and 2 coyotes. Our results show that the three canid species possess mutations that likely make this gene non functional. Because these mutations are shared across the three species, they must have appeared prior to the split of the wolf and the coyote, millions of years ago, and are not related to domestication. In addition, our results suggest that in these three canid species recombination does not occur at hotspots or hotspot location is controlled through a mechanism yet to be determined.     

Characterization of Prdm9 in Equids and Sterility in Mules

Prdm9 (Meisetz) is the first speciation gene discovered in vertebrates conferring reproductive isolation. This locus encodes a meiosis-specific histone H3 methyltransferase that specifies meiotic recombination hotspots during gametogenesis. Allelic differences in Prdm9, characterized for a variable number of zinc finger (ZF) domains, have been associated with hybrid sterility in male house mice via spermatogenic failure at the pachytene stage. The mule, a classic example of hybrid sterility in mammals also exhibits a similar spermatogenesis breakdown, making Prdm9 an interesting candidate to evaluate in equine hybrids. In this study, we characterized the Prdm9 gene in all species of equids by analyzing sequence variation of the ZF domains and estimating positive selection. We also evaluated the role of Prdm9 in hybrid sterility by assessing allelic differences of ZF domains in equine hybrids. We found remarkable variation in the sequence and number of ZF domains among equid species, ranging from five domains in the Tibetan kiang and Asiatic wild ass, to 14 in the Grevy’s zebra. Positive selection was detected in all species at amino acid sites known to be associated with DNA-binding specificity of ZF domains in mice and humans. Equine hybrids, in particular a quartet pedigree composed of a fertile mule showed a mosaic of sequences and number of ZF domains suggesting that Prdm9 variation does not seem by itself to contribute to equine hybrid sterility.     

Zinc Finger Binding Motifs Do Not Explain Recombination Rate Variation within or between Species of Drosophila

In humans and mice, the Cys₂His₂ zinc finger protein PRDM9 binds to a DNA sequence motif enriched in hotspots of recombination, possibly modifying nucleosomes, and recruiting recombination machinery to initiate Double Strand Breaks (DSBs). However, since its discovery, some researchers have suggested that the recombinational effect of PRDM9 is lineage or species specific. To test for a conserved role of PRDM9-like proteins across taxa, we use the Drosophila pseudoobscura species group in an attempt to identify recombination associated zinc finger proteins and motifs. We leveraged the conserved amino acid motifs in Cys₂His₂ zinc fingers to predict nucleotide binding motifs for all Cys₂His₂ zinc finger proteins in Drosophila pseudoobscura and identified associations with empirical measures of recombination rate. Additionally, we utilized recombination maps from D. pseudoobscura and D. miranda to explore whether changes in the binding motifs between species can account for changes in the recombination landscape, analogous to the effect observed in PRDM9 among human populations. We identified a handful of potential recombination-associated sequence motifs, but the associations are generally tenuous and their biological relevance remains uncertain. Furthermore, we found no evidence that changes in zinc finger DNA binding explains variation in recombination rate between species. We therefore conclude that there is no protein with a DNA sequence specific human-PRDM9-like function in Drosophila. We suggest these findings could be explained by the existence of a different recombination initiation system in Drosophila.     

Evidence of positive selection and concerted evolution in the rapidly evolving PRDM9 zinc finger domain in goats and sheep

Meiotic recombination contributes to augmentation of genetic diversity, exclusion of deleterious alleles and proper segregation of chromatids. PRDM9 has been identified as the gene responsible for specifying the location of recombination hotspots during meiosis and is also the only known vertebrate gene associated with reproductive isolation between species. PRDM9 encodes a protein with a highly variable zinc finger (ZF) domain that varies between as well as within species. In the present study, the ZF domain of PRDM9 on chromosome 1 was characterized for the first time in 15 goat breeds and 25 sheep breeds of India. A remarkable variation in the number and sequence of ZF domains was observed. The number of ZF repeats in the ZF array varied from eight to 12 yielding five homozygous and 10 heterozygous genotypes. The number of different ZF domains was 84 and 52 producing 36 and 26 unique alleles in goats and sheep respectively. The posterior mean of dN/dS or omega values were calculated using the codeml tool of pamlx to identify amino acids that are evolving positively in goats and sheep, as positions ?1, +3 and +6 in the ZF domain have been reported to experience strong positive selection across different lineages. Our study identified sites ?5, ?1, +3, +4 and +6 to be experiencing positive selection. Small ruminant zinc fingers were also found to be evolving under concerted evolution. Our results demonstrate the existence of a vast diversity of PRDM9 in goats and sheep, which is in concert with reports in many metazoans.     

Genetic Recombination Is Targeted towards Gene Promoter Regions in Dogs

The identification of the H3K4 trimethylase, PRDM9, as the gene responsible for recombination hotspot localization has provided considerable insight into the mechanisms by which recombination is initiated in mammals. However, uniquely amongst mammals, canids appear to lack a functional version of PRDM9 and may therefore provide a model for understanding recombination that occurs in the absence of PRDM9, and thus how PRDM9 functions to shape the recombination landscape. We have constructed a fine-scale genetic map from patterns of linkage disequilibrium assessed using high-throughput sequence data from 51 free-ranging dogs, Canis lupus familiaris. While broad-scale properties of recombination appear similar to other mammalian species, our fine-scale estimates indicate that canine highly elevated recombination rates are observed in the vicinity of CpG rich regions including gene promoter regions, but show little association with H3K4 trimethylation marks identified in spermatocytes. By comparison to genomic data from the Andean fox, Lycalopex culpaeus, we show that biased gene conversion is a plausible mechanism by which the high CpG content of the dog genome could have occurred.