Variation in histone configurations correlates with gene expression across nine inbred strains of mice.

The Diversity Outbred (DO) mice and their inbred founders are widely used models of human disease. However, although the genetic diversity of these mice has been well documented, their epigenetic diversity has not. Epigenetic modifications, such as histone modifications and DNA methylation, are important regulators of gene expression and, as such, are a critical mechanistic link between genotype and phenotype. Therefore, creating a map of epigenetic modifications in the DO mice and their founders is an important step toward understanding mechanisms of gene regulation and the link to disease in this widely used resource. To this end, we performed a strain survey of epigenetic modifications in hepatocytes of the DO founders. We surveyed four histone modifications (H3K4me1, H3K4me3, H3K27me3, and H3K27ac), as well as DNA methylation. We used ChromHMM to identify 14 chromatin states, each of which represents a distinct combination of the four histone modifications. We found that the epigenetic landscape is highly variable across the DO founders and is associated with variation in gene expression across strains. We found that epigenetic state imputed into a population of DO mice recapitulated the association with gene expression seen in the founders, suggesting that both histone modifications and DNA methylation are highly heritable mechanisms of gene expression regulation. We illustrate how DO gene expression can be aligned with inbred epigenetic states to identify putative cis-regulatory regions. Finally, we provide a data resource that documents strain-specific variation in the chromatin state and DNA methylation in hepatocytes across nine widely used strains of laboratory mice.

Differential effects of two catalytic mutations on full-length PRDM9 and its isolated PR/SET domain reveal a case of pseudomodularity.

PRDM9 is a DNA-binding histone methyltransferase that designates and activates recombination hotspots in mammals by locally trimethylating lysines 4 and 36 of histone H3. In mice, we recently reported two independently produced point mutations at the same residue, Glu360Pro (Prdm9EP) and Glu360Lys (Prdm9EK), which severely reduce its H3K4 and H3K36 methyltransferase activities in vivo. Prdm9EP is slightly less hypomorphic than Prdm9EK, but both mutations reduce both the number and amplitude of PRDM9-dependent H3K4me3 and H3K36me3 peaks in spermatocytes. While both mutations cause infertility with complete meiotic arrest in males, Prdm9EP, but not Prdm9EK, is compatible with some female fertility. When we tested the effects of these mutations in vitro, both Prdm9EP and Prdm9EK abolished H3K4 and H3K36 methyltransferase activity in full-length PRDM9. However, in the isolated PRDM9 PR/SET domain, these mutations selectively compromised H3K36 methyltransferase activity, while leaving H3K4 methyltransferase activity intact. The difference in these effects on the PR/SET domain vs the full-length protein shows that PRDM9 is not an intrinsically modular enzyme; its catalytic domain is influenced by its tertiary structure and possibly by its interactions with DNA and other proteins in vivo. These two informative mutations illuminate the enzymatic chemistry of PRDM9, and potentially of PR/SET domains in general, reveal the minimal threshold of PRDM9-dependent catalytic activity for female fertility, and potentially have some practical utility for genetic mapping and genomics.

Sexual dimorphism in the meiotic requirement for PRDM9: A mammalian evolutionary safeguard.

In many mammals, genomic sites for recombination are determined by the histone methyltransferase PRMD9. Some mouse strains lacking PRDM9 are infertile, but instances of fertility or semifertility in the absence of PRDM9 have been reported in mice, canines, and a human female. Such findings raise the question of how the loss of PRDM9 is circumvented to maintain fertility. We show that genetic background and sex-specific modifiers can obviate the requirement for PRDM9 in mice. Specifically, the meiotic DNA damage checkpoint protein CHK2 acts as a modifier allowing female-specific fertility in the absence of PRDM9. We also report that, in the absence of PRDM9, a PRDM9-independent recombination system is compatible with female meiosis and fertility, suggesting sex-specific regulation of meiotic recombination, a finding with implications for speciation.

HELLS and PRDM9 form a pioneer complex to open chromatin at meiotic recombination hot spots.

Chromatin barriers prevent spurious interactions between regulatory elements and DNA-binding proteins. One such barrier, whose mechanism for overcoming is poorly understood, is access to recombination hot spots during meiosis. Here we show that the chromatin remodeler HELLS and DNA-binding protein PRDM9 function together to open chromatin at hot spots and provide access for the DNA double-strand break (DSB) machinery. Recombination hot spots are decorated by a unique combination of histone modifications not found at other regulatory elements. HELLS is recruited to hot spots by PRDM9 and is necessary for both histone modifications and DNA accessibility at hot spots. In male mice lacking HELLS, DSBs are retargeted to other sites of open chromatin, leading to germ cell death and sterility. Together, these data provide a model for hot spot activation in which HELLS and PRDM9 form a pioneer complex to create a unique epigenomic environment of open chromatin, permitting correct placement and repair of DSBs.

Histone methyltransferase PRDM9 is not essential for meiosis in male mice.

A hallmark of meiosis is the rearrangement of parental alleles to ensure genetic diversity in the gametes. These chromosome rearrangements are mediated by the repair of programmed DNA double-strand breaks (DSBs) as genetic crossovers between parental homologs. In mice, humans, and many other mammals, meiotic DSBs occur primarily at hotspots, determined by sequence-specific binding of the PRDM9 protein. Without PRDM9, meiotic DSBs occur near gene promoters and other functional sites. Studies in a limited number of mouse strains showed that functional PRDM9 is required to complete meiosis, but despite its apparent importance, Prdm9 has been repeatedly lost across many animal lineages. Both the reason for mouse sterility in the absence of PRDM9 and the mechanism by which Prdm9 can be lost remain unclear. Here, we explore whether mice can tolerate the loss of Prdm9 By generating Prdm9 functional knockouts in an array of genetic backgrounds, we observe a wide range of fertility phenotypes and ultimately demonstrate that PRDM9 is not required for completion of male meiosis. Although DSBs still form at a common subset of functional sites in all mice lacking PRDM9, meiotic outcomes differ substantially. We speculate that DSBs at functional sites are difficult to repair as a crossover and that by increasing the efficiency of crossover formation at these sites, genetic modifiers of recombination rates can allow for meiotic progression. This model implies that species with a sufficiently high recombination rate may lose Prdm9 yet remain fertile.

Tissue-Specific Trans Regulation of the Mouse Epigenome.

The epigenetic landscape varies greatly among cell types. Although a variety of writers, readers, and erasers of epigenetic features are known, we have little information about the underlying regulatory systems controlling the establishment and maintenance of these features. Here, we have explored how natural genetic variation affects the epigenome in mice. Studying levels of H3K4me3, a histone modification at sites such as promoters, enhancers, and recombination hotspots, we found tissue-specific trans-regulation of H3K4me3 levels in four highly diverse cell types: male germ cells, embryonic stem cells, hepatocytes, and cardiomyocytes. To identify the genetic loci involved, we measured H3K4me3 levels in male germ cells in a mapping population of 59 BXD recombinant inbred lines. We found extensive trans-regulation of H3K4me3 peaks, including six major histone quantitative trait loci (QTL). These chromatin regulatory loci act dominantly to suppress H3K4me3, which at hotspots reduces the likelihood of subsequent DNA double-strand breaks. QTL locations do not correspond with genes encoding enzymes known to metabolize chromatin features. Instead their locations match clusters of zinc finger genes, making these possible candidates that explain the dominant suppression of H3K4me3. Collectively, these data describe an extensive, set of chromatin regulatory loci that control the epigenetic landscape.

PRDM9 and Its Role in Genetic Recombination.

PRDM9 is a zinc finger protein that binds DNA at specific locations in the genome where it trimethylates histone H3 at lysines 4 and 36 at surrounding nucleosomes. During meiosis in many species, including humans and mice where PRDM9 has been most intensely studied, these actions determine the location of recombination hotspots, where genetic recombination occurs. In addition, PRDM9 facilitates the association of hotspots with the chromosome axis, the site of the programmed DNA double-strand breaks (DSBs) that give rise to genetic exchange between chromosomes. In the absence of PRDM9 DSBs are not properly repaired. Collectively, these actions determine patterns of genetic linkage and the possibilities for chromosome reorganization over successive generations.

PRDM9 interactions with other proteins provide a link between recombination hotspots and the chromosomal axis in meiosis.

In mammals, meiotic recombination occurs at 1- to 2-kb genomic regions termed hotspots, whose positions and activities are determined by PRDM9, a DNA-binding histone methyltransferase. We show that the KRAB domain of PRDM9 forms complexes with additional proteins to allow hotspots to proceed into the next phase of recombination. By a combination of yeast-two hybrid assay, in vitro binding, and coimmunoprecipitation from mouse spermatocytes, we identified four proteins that directly interact with PRDM9's KRAB domain, namely CXXC1, EWSR1, EHMT2, and CDYL. These proteins are coexpressed in spermatocytes at the early stages of meiotic prophase I, the limited period when PRDM9 is expressed. We also detected association of PRDM9-bound complexes with the meiotic cohesin REC8 and the synaptonemal complex proteins SYCP3 and SYCP1. Our results suggest a model in which PRDM9-bound hotspot DNA is brought to the chromosomal axis by the action of these proteins, ensuring the proper chromatin and spatial environment for subsequent recombination events.

The Meiotic Recombination Activator PRDM9 Trimethylates Both H3K36 and H3K4 at Recombination Hotspots In Vivo.

In many mammals, including humans and mice, the zinc finger histone methyltransferase PRDM9 performs the first step in meiotic recombination by specifying the locations of hotspots, the sites of genetic recombination. PRDM9 binds to DNA at hotspots through its zinc finger domain and activates recombination by trimethylating histone H3K4 on adjacent nucleosomes through its PR/SET domain. Recently, the isolated PR/SET domain of PRDM9 was shown capable of also trimethylating H3K36 in vitro, raising the question of whether this reaction occurs in vivo during meiosis, and if so, what its function might be. Here, we show that full-length PRDM9 does trimethylate H3K36 in vivo in mouse spermatocytes. Levels of H3K4me3 and H3K36me3 are highly correlated at hotspots, but mutually exclusive elsewhere. In vitro, we find that although PRDM9 trimethylates H3K36 much more slowly than it does H3K4, PRDM9 is capable of placing both marks on the same histone molecules. In accord with these results, we also show that PRDM9 can trimethylate both K4 and K36 on the same nucleosomes in vivo, but the ratio of K4me3/K36me3 is much higher for the pair of nucleosomes adjacent to the PRDM9 binding site compared to the next pair further away. Importantly, H3K4me3/H3K36me3-double-positive nucleosomes occur only in regions of recombination: hotspots and the pseudoautosomal (PAR) region of the sex chromosomes. These double-positive nucleosomes are dramatically reduced when PRDM9 is absent, showing that this signature is PRDM9-dependent at hotspots; the residual double-positive nucleosomes most likely come from the PRDM9-independent PAR. These results, together with the fact that PRDM9 is the only known mammalian histone methyltransferase with both H3K4 and H3K36 trimethylation activity, suggest that trimethylation of H3K36 plays an important role in the recombination process. Given the known requirement of H3K36me3 for double strand break repair by homologous recombination in somatic cells, we suggest that it may play the same role in meiosis.

Health and population effects of rare gene knockouts in adult humans with related parents.

Examining complete gene knockouts within a viable organism can inform on gene function. We sequenced the exomes of 3222 British adults of Pakistani heritage with high parental relatedness, discovering 1111 rare-variant homozygous genotypes with predicted loss of function (knockouts) in 781 genes. We observed 13.7% fewer homozygous knockout genotypes than we expected, implying an average load of 1.6 recessive-lethal-equivalent loss-of-function (LOF) variants per adult. When genetic data were linked to the individuals' lifelong health records, we observed no significant relationship between gene knockouts and clinical consultation or prescription rate. In this data set, we identified a healthy PRDM9-knockout mother and performed phased genome sequencing on her, her child, and control individuals. Our results show that meiotic recombination sites are localized away from PRDM9-dependent hotspots. Thus, natural LOF variants inform on essential genetic loci and demonstrate PRDM9 redundancy in humans.

Multimer Formation Explains Allelic Suppression of PRDM9 Recombination Hotspots.

Genetic recombination during meiosis functions to increase genetic diversity, promotes elimination of deleterious alleles, and helps assure proper segregation of chromatids. Mammalian recombination events are concentrated at specialized sites, termed hotspots, whose locations are determined by PRDM9, a zinc finger DNA-binding histone methyltransferase. Prdm9 is highly polymorphic with most alleles activating their own set of hotspots. In populations exhibiting high frequencies of heterozygosity, questions remain about the influences different alleles have in heterozygous individuals where the two variant forms of PRDM9 typically do not activate equivalent populations of hotspots. We now find that, in addition to activating its own hotspots, the presence of one Prdm9 allele can modify the activity of hotspots activated by the other allele. PRDM9 function is also dosage sensitive; Prdm9+/- heterozygous null mice have reduced numbers and less active hotspots and increased numbers of aberrant germ cells. In mice carrying two Prdm9 alleles, there is allelic competition; the stronger Prdm9 allele can partially or entirely suppress chromatin modification and recombination at hotspots of the weaker allele. In cell cultures, PRDM9 protein variants form functional heteromeric complexes which can bind hotspots sequences. When a heteromeric complex binds at a hotspot of one PRDM9 variant, the other PRDM9 variant, which would otherwise not bind, can still methylate hotspot nucleosomes. We propose that in heterozygous individuals the underlying molecular mechanism of allelic suppression results from formation of PRDM9 heteromers, where the DNA binding activity of one protein variant dominantly directs recombination initiation towards its own hotspots, effectively titrating down recombination by the other protein variant. In natural populations with many heterozygous individuals, allelic competition will influence the recombination landscape.

Affinity-seq detects genome-wide PRDM9 binding sites and reveals the impact of prior chromatin modifications on mammalian recombination hotspot usage.

BACKGROUND: Genetic recombination plays an important role in evolution, facilitating the creation of new, favorable combinations of alleles and the removal of deleterious mutations by unlinking them from surrounding sequences. In most mammals, the placement of genetic crossovers is determined by the binding of PRDM9, a highly polymorphic protein with a long zinc finger array, to its cognate binding sites. It is one of over 800 genes encoding proteins with zinc finger domains in the human genome. RESULTS: We report a novel technique, Affinity-seq, that for the first time identifies both the genome-wide binding sites of DNA-binding proteins and quantitates their relative affinities. We have applied this in vitro technique to PRDM9, the zinc-finger protein that activates genetic recombination, obtaining new information on the regulation of hotspots, whose locations and activities determine the recombination landscape. We identified 31,770 binding sites in the mouse genome for the PRDM9(Dom2) variant. Comparing these results with hotspot usage in vivo, we find that less than half of potential PRDM9 binding sites are utilized in vivo. We show that hotspot usage is increased in actively transcribed genes and decreased in genomic regions containing H3K9me2/3 histone marks or bound to the nuclear lamina. CONCLUSIONS: These results show that a major factor determining whether a binding site will become an active hotspot and what its activity will be are constraints imposed by prior chromatin modifications on the ability of PRDM9 to bind to DNA in vivo. These constraints lead to the presence of long genomic regions depleted of recombination.

Nuclear localization of PRDM9 and its role in meiotic chromatin modifications and homologous synapsis.

Developmental progress of germ cells through meiotic phases is closely tied to ongoing meiotic recombination. In mammals, recombination preferentially occurs in genomic regions known as hotspots; the protein that activates these hotspots is PRDM9, containing a genetically variable zinc finger (ZNF) domain and a PR-SET domain with histone H3K4 trimethyltransferase activity. PRDM9 is required for fertility in mice, but little is known about its localization and developmental dynamics. Application of spermatogenic stage-specific markers demonstrates that PRDM9 accumulates in male germ cell nuclei at pre-leptonema to early leptonema but is no longer detectable in nuclei by late zygonema. By the pachytene stage, PRDM9-dependent histone H3K4 trimethyl marks on hotspots also disappear. PRDM9 localizes to nuclei concurrently with the deposition of meiotic cohesin complexes, but is not required for incorporation of cohesin complex proteins into chromosomal axial elements, or accumulation of normal numbers of RAD51 foci on meiotic chromatin by late zygonema. Germ cells lacking PRDM9 exhibit inefficient homology recognition and synapsis, with aberrant repair of meiotic DNA double-strand breaks and transcriptional abnormalities characteristic of meiotic silencing of unsynapsed chromatin. Together, these results on the developmental time course for nuclear localization of PRDM9 establish its direct window of function and demonstrate the independence of chromosome axial element formation from the concurrent PRDM9-mediated activation of recombination hotspots.

A multi-megabase copy number gain causes maternal transmission ratio distortion on mouse chromosome 2.

Significant departures from expected Mendelian inheritance ratios (transmission ratio distortion, TRD) are frequently observed in both experimental crosses and natural populations. TRD on mouse Chromosome (Chr) 2 has been reported in multiple experimental crosses, including the Collaborative Cross (CC). Among the eight CC founder inbred strains, we found that Chr 2 TRD was exclusive to females that were heterozygous for the WSB/EiJ allele within a 9.3 Mb region (Chr 2 76.9 - 86.2 Mb). A copy number gain of a 127 kb-long DNA segment (designated as responder to drive, R2d) emerged as the strongest candidate for the causative allele. We mapped R2d sequences to two loci within the candidate interval. R2d1 is located near the proximal boundary, and contains a single copy of R2d in all strains tested. R2d2 maps to a 900 kb interval, and the number of R2d copies varies from zero in classical strains (including the mouse reference genome) to more than 30 in wild-derived strains. Using real-time PCR assays for the copy number, we identified a mutation (R2d2WSBdel1) that eliminates the majority of the R2d2WSB copies without apparent alterations of the surrounding WSB/EiJ haplotype. In a three-generation pedigree segregating for R2d2WSBdel1, the mutation is transmitted to the progeny and Mendelian segregation is restored in females heterozygous for R2d2WSBdel1, thus providing direct evidence that the copy number gain is causal for maternal TRD. We found that transmission ratios in R2d2WSB heterozygous females vary between Mendelian segregation and complete distortion depending on the genetic background, and that TRD is under genetic control of unlinked distorter loci. Although the R2d2WSB transmission ratio was inversely correlated with average litter size, several independent lines of evidence support the contention that female meiotic drive is the cause of the distortion. We discuss the implications and potential applications of this novel meiotic drive system.

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 Prdm9Cst allele in two Mus musculus subspecies, M.m. castaneus, in which Prdm9Cst arose, and M.m. domesticus, into which Prdm9Cst 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 PRDM9Cst 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.

PRDM9 binding organizes hotspot nucleosomes and limits Holliday junction migration.

In mammals, genetic recombination during meiosis is limited to a set of 1- to 2-kb regions termed hotspots. Their locations are predominantly determined by the zinc finger protein PRDM9, which binds to DNA in hotspots and subsequently uses its SET domain to locally trimethylate histone H3 at lysine 4 (H3K4me3). This sets the stage for double-strand break (DSB) formation and reciprocal exchange of DNA between chromatids, forming Holliday junctions. Here we report genome-wide analyses of PRDM9-dependent histone modifications using two inbred mouse strains differing only in their PRDM9 zinc finger domain. We show that PRDM9 binding actively reorganizes nucleosomes into a symmetrical pattern, creating an extended nucleosome-depleted region. These regions are centered by a consensus PRDM9 binding motif, whose location and identity was confirmed in vitro. We also show that DSBs are centered over the PRDM9 binding motif within the nucleosome-depleted region. Combining these results with data from genetic crosses, we find that crossing-over is restricted to the region marked by H3K4me3. We suggest that PRDM9-modified nucleosomes create a permissible environment that first directs the location of DSBs and then defines the boundaries of Holliday junction branch migration.

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 Prdm9Cst and one activated by Prdm9Dom2, 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 Prdm9Cst activated hotspots. The binding specificity of each position in the Hlx1 binding site, activated by Prdm9Cst, 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.

Meiotic DSBs and the control of mammalian recombination.

The laboratories of Galina Petukhova and R Daniel Camerini-Otero have achieved significant technical advances in determining the genome-wide sites of DNA double-strand breaks (DSBs) where the process of genetic exchange between chromatids during meiosis begins. Applying the new approaches to male meiosis in mice, their experimental results considerably increase our insights into the nature and regulation of these processes.

Clusters of ancestrally related genes that show paralogy in whole or in part are a major feature of the genomes of humans and other species.

Arrangements of genes along chromosomes are a product of evolutionary processes, and we can expect that preferable arrangements will prevail over the span of evolutionary time, often being reflected in the non-random clustering of structurally and/or functionally related genes. Such non-random arrangements can arise by two distinct evolutionary processes: duplications of DNA sequences that give rise to clusters of genes sharing both sequence similarity and common sequence features and the migration together of genes related by function, but not by common descent. To provide a background for distinguishing between the two, which is important for future efforts to unravel the evolutionary processes involved, we here provide a description of the extent to which ancestrally related genes are found in proximity.Towards this purpose, we combined information from five genomic datasets, InterPro, SCOP, PANTHER, Ensembl protein families, and Ensembl gene paralogs. The results are provided in publicly available datasets (http://cgd.jax.org/datasets/clustering/paraclustering.shtml) describing the extent to which ancestrally related genes are in proximity beyond what is expected by chance (i.e. form paraclusters) in the human and nine other vertebrate genomes, as well as the D. melanogaster, C. elegans, A. thaliana, and S. cerevisiae genomes. With the exception of Saccharomyces, paraclusters are a common feature of the genomes we examined. In the human genome they are estimated to include at least 22% of all protein coding genes. Paraclusters are far more prevalent among some gene families than others, are highly species or clade specific and can evolve rapidly, sometimes in response to environmental cues. Altogether, they account for a large portion of the functional clustering previously reported in several genomes.

Mammalian recombination hot spots: properties, control and evolution.

Recombination, together with mutation, generates the raw material of evolution, is essential for reproduction and lies at the heart of all genetic analysis. Recent advances in our ability to construct genome-scale, high-resolution recombination maps and new molecular techniques for analysing recombination products have substantially furthered our understanding of this important biological phenomenon in humans and mice: from describing the properties of recombination hot spots in male and female meiosis to the recombination landscape along chromosomes. This progress has been accompanied by the identification of trans-acting systems that regulate the location and relative activity of individual hot spots.