Parental origin of chromosomes influences crossover activity within the Kcnq1 transcriptionally imprinted domain of Mus musculus

Background Among the three functions of DNA, mammalian replication and transcription can be subject to epigenetic imprinting specified by the parental origin of chromosomes, and although there is suggestive indication that this is also true for meiotic recombination, no definitive evidence has yet been reported. Results We have now obtained such evidence on mouse chromosome 7 by assaying meiotic recombination as it occurs in reciprocal F1 mice. A 166 kb region near the Kcnq1 transcriptionally imprinted domain showed significantly higher recombination activity in the CAST×B6 parental direction (p < 0.03). Characterizing hotspots within this domain revealed a cluster of three hotspots lying within a 100 kb span, among these hotspots, Slc22a18 showed a definitive parent of origin effect on recombination frequency (p < 0.02). Comparing recombination activity in the mouse Kcnq1 and neighboring H19-Igf2 imprinted domains with their human counterparts, we found that elevated recombination activity in these domains is a consequence of their chromosomal position relative to the telomere and not an intrinsic characteristic of transcriptionally imprinted domains as has been previously suggested. Conclusion Similar to replication and transcription, we demonstrate that meiotic recombination can be subjected to epigenetic imprinting and hotspot activity can be influenced by the parental origin of chromosomes. Furthermore, transcriptionally imprinted regions exhibiting elevated recombination activity are likely a consequence of their chromosomal location rather than their transcriptional characteristic.


Background
DNA serves three major functions; it is replicated, providing the material for hereditary transmission from one generation to the next; it is transcribed, expressing its stored information as a variety of RNA products; and it undergoes meiotic recombination, generating population varia-tion and substrates for evolution. Two of these processes, replication [1,2] and transcription [3], can be subject to parentally determined epigenetic modification, generally known as genomic imprinting, and recombination has also been postulated to be affected by imprinting [4][5][6].
Transcriptional imprinting is characterized by the silencing of one parental-specific allele in the offspring [7][8][9], and asynchronous DNA replication due to the parental origin of chromosomes has been observed in some transcriptionally imprinted domains [1,2,10]. While the mechanisms of genomic imprinting and its epigenetic regulation of gene expression and replication have been extensively studied (for a review see [3]), a direct demonstration of epigenetic controls in meiotic recombination has not been available.
Among the three functions of DNA, meiotic recombination is the only process that involves physical interactions between the two parental chromosomes. Meiotic recombination begins in meiosis I with a Spo11 mediated double strand break (DSB) [11]. This break is subsequently repaired via either the double strand break repair pathway, where DNA sequences flanking the DSB are exchanged between non-sister chromatids resulting in a crossover, or the synthesis-dependent strand annealing pathway in which the DSB is repaired using a non-sister chromatid as the copying template generating a short region of gene conversion without the exchange of flanking sequences (a non-crossover). For a review of the recombination pathways see [11,12].
In mammals, where recombination initiates prior to synapsis, meiotic crossovers occur at preferred 1-2 kb regions, known as recombination hotspots. Currently, no comprehensive model exists that accounts for the location and activity of recombination hotspots. Regional differences in recombination frequencies have been associated with both transcriptionally imprinted regions and regions near telomeres. Telomeric regions in both humans and mice generally have higher male than female recombination activity [13][14][15] and imprinted regions located near telomeres, such as the human H19-Igf2 and Kcnq1 domains, show higher male recombination activity [16].
The mouse H19-Igf2 and Kcnq1 imprinted domains, located on the distal region of chromosome 7, provide a model for examining if meiotic recombination hotspot can be subjected to epigenetic imprinting where recombination activity is affected by the parental origin of the recombining chromatids. Mapping recombination rates in this region, we found one region near the Kcnq1 transcriptionally imprinted domain subject to meiotic recombination imprinting as evidenced by recombination being much elevated in F1 animals arising from one parental direction of the cross versus the opposite. Fine mapping of recombination activity within this region to hotspot level resolution showed a cluster of three hotspots whose meiotic recombination activities are imprinted.
Additionally, while the imprinting mechanisms are identical in the H19-Igf2 and Kcnq1 imprinted domains in humans and mice [17,18], the genomic arrangement of these two domains relative to the telomere is inverted in the mouse genome. Our analysis of this region showed that in mice recombination activity is elevated only in the Kcnq1 domain, suggesting that recombination activities of hotspots within transcriptionally imprinted domains are influenced by their chromosomal position and not all transcriptionally imprinted domains exhibit elevated recombination activity.

High recombination was observed within the Kcnq1 imprinted domain but not within H19-Igf2 imprinted domain
We mapped crossover events within the H19-Igf2 and Kcnq1 transcriptionally imprinted domains on mouse chromosome 7 (142.366 -143.400 Mb, NCBI Build 36) in 5,914 meioses of a cross between the C57/BL6J and CAST/EiJ strains ( Figure 1). Only three crossovers were observed within the H19-Igf2 domain (142.366 -142.727 Mb), yielding a sex-averaged recombination rate of 0.14 cM/Mb, approximately 4 times lower than the genome average of 0.5 to 0.6 cM/Mb [6,15,19] (Figure  2A). In contrast, the Kcnq1 imprinted domain (142.727-143.400 Mb), encompassing genes between Th and Nap1l4, contained 119 recombinants yielding a crossover activity of 2.98 cM/Mb, five times the mouse genomewide average.
We further mapped the recombination activity within this latter region to an average resolution of about 172 kb and compared the recombination activity between reciprocal parental crosses (B6Ǩ × CASTǩ vs. CASTǨ × B6ǩ) to detect if the chromosomal origin of the parental chromosome affects meiotic recombination. A meiotic recombination imprinting effect was found in a 166 kb region near distal side of the Kcnq1 transcriptionally imprinted domain ( Figure 2B). This region contained a total of 31 recombinants with 7 recombinants arising in B6 × CAST F1 mice and 24 recombinants arising in CAST × B6 F1 mice (P FET = 0.027 after Bonferroni correction for multiple testing). None of the other intervals within the H19-Igf2 and Kcnq1 transcriptionally imprinted domains showed any significant differences in recombination activity between the two reciprocal parental crosses.

Fine mapping revealed five novel recombination hotspots within the Kcnq1 transcriptionally imprinted domain
To better understand whether meiotic imprinting influences recombination at a regional or hotspot level, we mapped the crossovers occurring in the Kcnq1 domain to hotspot resolution, revealing five highly active hotspots. These hotspots were named according to their closest neighboring gene ( Figure 3, Table 1). Hotspot Th was mapped to an 11 kb interval (the limit of resolution possible with available SNPs), while hotspots Kcnq1, Cdkn1c, Slc22a18 and Nap1l4 were mapped to a resolution of 2.3 kb or less. We considered the 11 kb Th interval as a single hotspot as it is unlikely that multiple hotspots will be found within such a short interval [6]. Of the five hotspots, three (Th, Cdkn1c and Nap114) were located in intergenic regions and two hotspots were found within introns; hotspot Kcnq1 within a large intron of 65 kb while hotspot Slc22a18 is located within a short 4 kb intron.
The interval displaying meiotic recombination imprinting encompassed hotspots Cdkn1c, Slc22a18 and Nap1l4. These hotspots clustered within a 100 kb region (NCBI build 36, 143.27 -143.37 Mb), and, as a cluster, all three hotspots showed elevated recombination activity in CASTǨ × B6ǩ F1 mice regardless of whether the recombination occurs in the male or female germline (P FET = 0.005; the number of recombinants is slightly less than those detected at the regional level as some crossover events are not associated with any particular hotspot). This differential recombination activity between reciprocal crosses was statistically significant for Slc22a18 (P FET = 0.02 without Bonferroni correction and P FET = 0.10 with Bonferroni correction), however, the number of crossovers at each of flanking hotspots, although consistent with a parent of origin effect, was not sufficient to provide low p values (Table 2). Individually, hotspot Cdkn1c and Nap1l4 showed similar rate of recombination when the recombination event occurs in the female F1 animals,

Meiotic recombination imprinting and possible chromatid of origin effects
One possible site of action for meiotic recombination imprinting is the choice of the chromatid on which the initiating DSB occurs. When DSBs are initiated on one chromatid, the subsequent repair processes use the homologous chromatid as a template, resulting in conversion of short sequence stretch at the center of the hotspot to the donor sequence [20]. This makes it possible to determine the chromatid on which the initiating double strand break occurs when informative SNPs across the hotspot are available. For hotspot Kcnq1, which is not recombinationally imprinted, the initiating chromatid showed no evidence of parental imprinting. However, at hotspot Cdkn1c, which is recombinationally imprinted, the limited number of crossovers suggests a possible initiation bias on the paternal chromosome.
At hotspot Kcnq1, the molecular details of crossing over were typed among 49 offspring using three additional SNP markers. Among progeny of both reciprocal parents, the central marker was three times more likely to be of CAST genotype in both parental directions suggesting that DSB initiation occurs preferentially on the B6 chromosome regardless of its parental origin; thus, no meiotic recombination imprinting was detectable at this hotspot in our study (Table 3 and 4). Among the hotspots where meiotic recombination imprinting was detected, only Cdkn1c contained informative SNPs spanning across the region of crossover exchange, and the seven recombinants at this hotspot were typed for these internal SNPs (Table 3 and 5). Limited by the small number of crossovers occurring at this hotspot, we found only suggestive confirmation that direction of the parental cross influenced the choice of the initiating chromosome with DSBs preferentially initiating on the paternal chromosome (P-value < 0.13, Table 3).

Parental origin of chromosomes can influence the recombination activity at specific hotspots
When recombination was mapped along the entirety of mouse chromosome 1, a significant excess of recombina-tion intervals was observed in which preferential recombination occurred in one parental direction of the cross or the other, but no single interval showed a statistically significant preference by itself [6].
Mapping the recombination activity on chromosome 7 at a regional scale, we found an interval showing definitive evidence of meiotic recombination imprinting. Fine mapping revealed meiotic recombination imprinting influences at the hotspot level involving hotspots Cdkn1c, Slc22a18 and Nap1l4. Their close proximity, clustering within a 100 kb region, suggests that the parental origin of chromosomes exerts a regional influence on meiotic recombination activity (Figure 3). While these hotspots with imprinted meiotic recombination activity are located in a transcriptionally imprinted region, it will require additional examples of recombination imprinting to clarify the extent to which meiotic recombination imprinting overlaps with and/or is related to transcriptional imprinting. It is quite possible that some hotspots with imprinted meiotic recombination activity are not located within a transcriptionally imprinted domain. Another possibility, though remote, is that an imprinted trans-acting factor controls the recombination activities of hotspots showing parent of origin effects.
While DNA methylation is a well studied epigenetic marker controlling transcription of imprinted genes, the methylation pattern is reportedly erased and reestablished prior to the beginning of meiosis in male [21-23] but erased and not yet established prior to meiotic recombination in females [24]. The observed recombination activity imprinted in the CASTǨ × B6ǩ parental direction is unlikely to be a consequence of DNA methylation, and another imprinting epigenetic marker may be responsible for differentiating the two parental chromosomes [25,26]. In considering possible mechanisms of meiotic recombination imprinting, it is also worth noting that it is only during recombination, not replication or transcription, Hotspots Slc22a18 and its flanking hotspots, Cdkn1c and Nap1l4, clustered within 100 kb and showed a regionally significant increase in crossover activity in F1 animals from CASTǨ × B6ǩ parents. Both sexes showed higher number of recombinant in the CASTǨ × B6ǩ parental cross within the 100 kb cluster with a C×B:B×C ratio of 3.8 for males and 2.5 for females. that homologous chromosomes come into physical contact with each other. For transcriptional imprinting and replication imprinting there is no indication of a sex bias where the imprinted pattern is influenced by the sex of an individual, and similarly we do did not find sex to influence meiotic recombination imprinting.

DSB initiation may be influenced by epigenetic imprinting
If DSB initiation occurs on both parental chromosomes at equal frequencies, the number of DSBs will be constant regardless of the parental origin of chromosomes, and no influence on the overall recombination activity would be observed. For example, if the DSB initiating recombination occurs preferentially on the B6 chromatid but epigenetic imprinting does not play a role, one would expect equal rates of recombination in B6Ǩ × CASTǩ and CASTǨ × B6ǩ F1 mice. Conversely, if the DSB is equally likely to occur on a B6 or CAST chromatid and epigenetic imprinting does occur, suppressing DSB formation on the maternal chromosome, the decrease in DSBs would affect both the maternal B6 chromatid in B6Ǩ × CASTǩ mice and the maternal CAST chromatid in CASTǨ × B6ǩ mice equally and no change in recombination frequency would be observed. Thus, detection of meiotic recombination imprinting requires both a strong preference for recombination to initiate on one specific parental chromatid and that this region of the genome is subject to epigenetic imprinting. In this latter case, if DSBs initiate preferentially on the B6 chromatid (indeed, such preference has been previously observed in human and mouse [6,27]) and epigenetic imprinting suppresses DSB initiation on the maternal chromatid, one would observed a decrease in recombination activity in the B6Ǩ × CASTǩ mice as the number of DSB initiation events on the B6 maternal chromatid (and the subsequent recombination activity) is reduced. In contrast, the number of recombination events will remain unchanged in the CASTǨ × B6ǩ mice where Examining 49 recombinant offspring at the Kcnq1 hotspot showed DSBs preferentially occurs on the B6 chromosome, regardless of the chromosome's parental origin. SNP position indicated in parentheses (base pairs -Build 36).  the maternal CAST chromatid is suppressed but DSB initiation on the paternal B6 chromatid in unaffected.
Although our data is limited by the number of recombinants observed at hotspot Cdkn1c, even with 5,914 meioses, they do lean in the direction of suggesting that meiotic recombination imprinting regulates recombination activity by directing DSBs preferentially towards one chromatid. It is possible that DSB initiation is favored towards the paternal chromosome at all of the imprinted hotspots (Cdkn1c, Slc22a18 and Nap1l4) clustered within 100 kb; unfortunately the lack of suitable SNPs precluded testing the remaining hotspots within this region. This is the first detailed examination of meiotic recombination imprinting at a specific hotspot.
We should also point out two alternative theoretical possibilities that, although unlikely, could explain our results. There could exist a trans-acting gene that is both subject to imprinting and whose B6 and CAST alleles differ in their ability to control recombination in the Slc22a18 region. This would require imprinting of gene function during meiosis, shifting the imprinted region from Slc22a18 to another site. It is also possible that a trans-acting gene controlling recombination in the Slc22a18 region is located on the X chromosome as reciprocal F1 males are X B6 Y CAST v. X CAST Y B6 ; again, the B6 and CAST alleles of this putative gene would have to differ in their activating ability. In this context, we should point out that trans control of some hotspot activities, albeit by an autosomal locus, has been reported [28,29].
We have expanded our analysis to over 10 Mb near the distal end of chromosome 7, but no other region showed a meiotic recombination imprinting effect bias similar to those observed near the Kcnq1 domain (data not shown). It is likely that meiotic imprinted hotspots are not a common occurrence, and the challenge is now to find additional hotspots subject to meiotic recombination imprinting and characterize their molecular mechanism.

Overall recombination activity in the H19-Igf2 and Kcnq1 imprinted domains is dependent upon genomic positioning and not transcriptional imprinting
Human transcriptionally imprinted regions have been associated with elevated recombination frequencies [4,5,16,30]. In humans, both the H19-Igf2 and Kcnq1 transcriptionally imprinted domains showed elevated male recombination activities [16]. The low crossover activity at the mouse H19-Igf2 region is in marked contrast to the human H19-Igf2 recombination pattern. In addition, the difference in recombination activity between the two transcriptionally imprinted domains in mice indicates that elevated recombination is not necessarily correlated with transcriptional imprinting.
It is likely that the location of hotspots relative to the telomeres has a greater influence on recombination activity. It appears that sex-biased recombination activity (or any recombination activity) is not a general characteristic of mammalian transcriptionally imprinted domains, and recombination within these imprinted regions is likely influenced by their relative position to the telomeres.

Conclusion
We mapped five new recombination hotspots in the Kcnq1 imprinted domain, and comparing recombination activity of F1 animals from reciprocal parents obtained evidence that meiotic recombination imprinting can influence crossover activity at a cluster of three closely spaced hotspots. These results suggest that, like replication and transcription, meiotic recombination can be subjected to epigenetic imprinting. However, any epigenetic marker for recombination is likely independent of DNA methylation as the methylation pattern of germ cells are reported to be erased and re-established prior to meiosis. Recombination activity in the H19-Igf2 and Kcnq1 chromosomal regions showed that transcriptionally imprinted mammalian regions do not necessarily have elevated recombination activity. The higher recombination activity observed at Kcnq1 is likely a consequence of their chromosomal location rather than their imprinted characteristic.

Mouse Strain and Spleen DNA extraction
Mice were purchased from The Jackson Laboratory. The DNA samples used in this study were previously described elsewhere [6].