Advertisement

High‐resolution sperm typing of meiotic recombination in the mouse MHC Eβ gene

C.L. Yauk, P.R.J Bois, A.J. Jeffreys

Author Affiliations

  1. C.L. Yauk*,2,
  2. P.R.J Bois3 and
  3. A.J. Jeffreys1
  1. 1 Department of Genetics, University of Leicester, Leicester, LE1 7RH, UK
  2. 2 Present address: Environmental Health Centre, Health Canada, Tunney's Pasture, Ottawa, Ontario, Canada, K1A 0L2
  3. 3 Present address: St Jude Children's Research Hospital, Department of Genetics, 332 N. Lauderdale, Memphis, TN, 38105, USA
  1. *Corresponding author. E-mail: Carole_Yauk{at}hc-sc.gc.ca
View Full Text

Abstract

Meiotic crossovers detected by pedigree analysis in the mouse MHC cluster into hotspots. To explore the properties of hotspots, we subjected the class II Eβ gene to high‐resolution sperm crossover analysis. We confirm the presence of a highly localized hotspot 1.0–1.6 kb wide in the second intron of Eβ and show that it is flanked by DNA which is almost completely recombinationally inert. Mice heterozygous for haplotype s and another MHC haplotype show major haplotype‐dependant variation in crossover rate but always the same hotspot, even in crosses including the highly diverged p haplotype. Crossovers in reciprocal orientations occur at similar rates but show different distributions across the hotspot, with the position of centre points in the two orientations shifted on average by 400 bp. This asymmetry results in crossover products showing biased gene conversion in favour of hotspot markers from the non‐initiating haplotype, and supports the double‐strand break repair model of recombination, with haplotype s as the most efficient crossover initiator. The detailed behaviour of the Eβ hotspot, including evidence for highly localized recombination initiation, is strikingly similar to human hotspots.

Introduction

There is growing evidence that meiotic crossovers are not distributed randomly along mammalian chromosomes. Rather, mammalian genomes appear, at least to some extent, to be a mosaic of recombinationally suppressed regions separated by intervals that undergo very elevated and targeted crossover activity (Pittman and Schimenti, 1998; Petes, 2001; Yu et al., 2001). Recently, high‐resolution sperm typing has shown that crossovers in the human MHC tend to cluster into very narrow (1–2 kb) recombination hotspots separated by long stretches of DNA exhibiting strong linkage disequilibrium (Jeffreys et al., 2000, 2001). Similarly, meiotic recombination events in the mouse MHC also show strong clustering. Eight mouse hotspots have been identified so far by breakpoint mapping in recombinant inbred lines and in pedigrees. These hotspots vary in the degree to which they have been characterized, from the EαC4 region at 300 kb in length (Lafuse et al., 1989), down to the Eα hotspot ∼0.5 kb wide (Khambata et al., 1996).

Mouse studies have also shown that hotspot location and/or intensity can be influenced by MHC haplotype, as well as genetic background (strain of mouse) and sex (Steinmetz et al., 1987; Shiroishi et al., 1991, 1993, 1995). For example, the very active hotspot in the class II Eβ gene, first characterized by Steinmetz et al. (1982), is ∼4 kb wide. Its activity varies between hybrid strains, and appears to be completely suppressed in heterozygotes carrying the p haplotype, in which crossovers occur instead preferentially at the Eα gene (Lafuse and David, 1986; Zimmerer and Passmore, 1991). Heine et al. (1994) found that these Eα haplotype p recombinants only result when p is from the intra‐H2 recombinant strain B10F(13R), suggesting that distal modifiers of crossover can influence hotspot activity. Similarly, crossover rates at the Psmb9 (formerly known as Lmp2) hotspot are governed not only by haplotype and genetic background, but also by the sex of the animal, with high female rates occurring only in some strains (Shiroishi et al., 1993, 1995).

While crossover breakpoint mapping in mice has provided valuable information about meiotic recombination, the use of recombinant inbred lines and back‐crossed offspring to identify MHC‐recombinant progeny is extremely laborious. The limited numbers of crossovers that can be detected prevent the detailed definition of hotspot morphology and of the recombination processes operating within hotspots. Similarly, little is known about the recombinational activity of DNA outside these mouse hotspots. These limitations can be overcome by recovering crossover products directly from sperm DNA (Jeffreys et al., 1998, 2001), as shown by Guillon and De Massy (2002), who mapped 69 sperm recombinants at the Psmb9 locus to show that crossover density was maximal in an interval of 210 bp and decreased roughly symmetrically on both sides, in a manner similar to human recombination hotspots. We now extend this approach to analyse recombination at the mouse Eβ locus and to investigate the influence of haplotype on male recombination rates and on the distribution and reciprocity of crossover events within the hotspot. We show that the Eβ hotspot shares remarkably similar properties with its counterparts in the human MHC, implying similar recombination processes in these two very different species.

Results

Detection of meiotic crossovers in sperm DNA

Re‐sequencing across the Eβ gene did not reveal any previously unidentified polymorphisms across the crossover hotspot region in the second intron (see Materials and methods). Mice homozygous for MHC haplotypes b, d, k or p were crossed with haplotype s homozygotes. Sperm DNA from a mature F1 hybrid from each cross was prepared and Eβ‐recombinant DNA molecules recovered from multiple pools of sperm DNA by nested repulsion‐phase allele‐specific PCR directed to heterozygous single nucleotide polymorphism (SNP) sites, as shown in Figure 1. Nested PCR gave sufficient allele specificity to allow the selective amplification of crossover molecules even from large pools of sperm DNA (up to 30 000 amplifiable molecules of each progenitor haplotype per PCR). Parallel analysis of brain DNA yielded much lower levels of PCR products, establishing that the sperm crossovers recovered were genuine products of meiotic recombination (Jeffreys et al., 2000) (Figure 1). Crossovers in reciprocal orientation were analysed either by using haplotype s 5′ allele‐specific primers plus non‐s 3′ primers (crossovers in orientation A) or 5′ non‐s primers plus 3′ s haplotype primers (orientation B). We analysed a total of 6.5 × 106 amplifiable molecules from the four crosses and isolated 516 crossover molecules.

Figure 1.

Recovery of sperm crossovers. (A) Batches of sperm DNA from a hybrid mouse with heterozygosities (black and white circles) across the test region are amplified using two rounds of allele‐specific PCR in repulsion phase (primers shown as black and white arrows) to selectively amplify recombinant molecules. Grey circles represent alleles of unknown status. (B) SNP typing is used to localize the breakpoint interval for each recombinant molecule amplified. (C) Examples of crossover molecules recovered at the Eβ locus from s × k sperm DNA at 280 molecules per reaction, electrophoresed next to a 1 kb marker (Gibco) (M), and alongside 840 molecules per reaction of brain DNA.

Sperm crossover rates in the Eβ region varied significantly between different crosses (Fisher's exact test, P <0.001 for each pair‐wise comparison; Table I), but not between individuals from the same cross (data not shown). The s × k heterozygotes showed the highest crossover rate (1.6 × 10−3 per sperm) and s × p showing the lowest (0.021 × 10−3). These rate differences were consistent both for orientation A crossovers and for orientation B exchanges (Table I), and are therefore reproducible. This establishes that recombination rates decrease in the order k > d > b >> p in hybrids carrying haplotype s.

View this table:
Table 1. Summary of sperm crossover data at the mouse MHC Eβ gene for separate orientation A and B crossovers as well as for combined A + B exchanges

Allele‐specific oligonucleotide mapping of crossover breakpoints

The selected polymorphisms (mainly SNPs) used to map crossover breakpoints by hybridizing amplified recombinant molecules with allele‐specific oligonucleotides (ASOs) (Wood et al., 1985) are shown in Table II. All 516 sperm crossovers were simple (Figure 1), with breakpoints mapping to a single inter‐marker interval. There were no instances of more complex exchanges involving switching back and forth between haplotypes near the site of exchange. Crossover breakpoints were non‐randomly distributed across the Eβ region in all hybrids (Figure 2). Nearly all crossovers (99%) occurred within the second intron of the Eβ gene and were highly concentrated in the 3′ half, demonstrating the presence of a highly localized hotspot. Crossover frequencies in each interval were used to determine the recombination activity per unit length in cM/Mb (Figure 2). Each hybrid showed an approximately normal distribution of combined orientation A + B crossover breakpoints across the hotspot (Figure 3A). All four heterozygotes showed the centre of the hotspot at a very similar location (Table I), within an MT element (Figure 2). The width of the hotspot was also similar in all crosses (1.0–1.6 kb), although the peak activity varied dramatically from 3.5 cM/Mb in the s × p hybrid to 170 cM/Mb in s × k. The same hotspot therefore functions in male recombination in all heterozygotes, despite the substantial rate variation. Only one crossover, located upstream of position 6114 bp, was found outside the hotspot. The mean recombination activity over all crosses in DNA flanking the hotspot is therefore extremely low, at 0.005 cM/Mb (95% CI 0.0002–0.02 cM/Mb) and ∼100‐fold below the mouse genome average rate of 0.5 cM/Mb.

Figure 2.

Mouse Eβ sperm crossover distributions for orientation A (s/non‐s) and orientation B (non‐s/s) recombinants. The haplotypes in each mouse analysed are given at top. The numbers of recombinants recovered in each interval, shown above the histograms, were used to calculate the crossover frequency per unit length of DNA for each interval in cM/Mb. Note the considerable differences in recombination scaling due to substantial variation in crossover activity across the hybrids tested. The Eβ gene is shown beneath the histograms. The locations of the allele‐specific primers used to recover crossovers are shown as vertical lines and the locations of the hotspot and MT element are indicated. Locations are relative to clone BAC 135G15 (DDBJ/EMBL/GenBank accession No. AF050157).

Figure 3.

Distribution of sperm crossovers across the Eβ hotspot. (A) Cumulative frequency distributions of sperm crossovers in the four hybrids analysed, for the combined orientation A + B recombinants. The least‐squares best‐fit cumulative frequency distribution assumes that crossovers are normally distributed across the hotspot (centre at 7260 bp, standard deviation 350 bp). (B) Cumulative frequency distributions of separate orientation A (filled squares) and B (open squares) crossovers. The best‐fit line for the combined A + B distributions is indicated. The B distributions are shifted 3′ relative to A distributions in all crosses.

View this table:
Table 2. List of ASO sequences

Reciprocal crossover asymmetry and transmission distortion

Orientation A and B crossovers were detected in each cross at similar frequencies, consistent with these sperm recombinants being the products of reciprocal crossover. Minor differences in A versus B rates within a cross (at most only 2.2‐fold for s × d; Table I) show no consistent bias towards A or B exchanges, and could easily have arisen from variation in PCR efficiency between different primer sets. While sperm crossovers appear to be reciprocal in rate, the separate A and B crossovers showed different distributions, with orientation A (s→non‐s) crossovers consistently displaced 5′‐ward relative to orientation B (non‐ss) exchanges in all four hybrids (Table I; Figures 2 and 3B). This displacement varied from 220 bp for s × k to 540 bp for s × p, and averaged ∼400 bp over all four heterozygotes.

This reciprocal crossover asymmetry results in transmission ratio distortion of alleles into crossover products. Without asymmetry, s and non‐s alleles should show a 50:50 Mendelian ratio in combined A and B crossovers. However, crossover asymmetry results in highly significant over‐transmission of non‐s alleles within the hotspot (Figure 4). The most extreme distortion was seen for marker 3.2 at 7247 bp, very close to the centre of the hotspot, which in s × p animals show a transmission ratio of 81:19 for non‐s:s alleles in A + B crossover products. This distortion also affects neighbouring markers, with the strength of distortion decreasing with distance from the centre of the hotspot. Superimposition of data from all four hybrids (Figure 4, bottom) showed that the transmission distortion profiles were similar in all animals. However, the maximum degree of distortion varied between crosses in the order p > (b, d) > k; this order is the inverse of the order of crossover rates in these heterozygotes (Table I).

Figure 4.

Transmission distortion in Eβ crossover products in each heterozygote. The proportions of non‐s alleles in equal numbers of orientation A and B sperm crossovers are shown for each marker, along with 95% CIs. The bottom graph shows all hybrids superimposed. The positions of the markers with respect to the Eβ gene are shown at bottom, and the centre of the hotspot is indicated with arrows.

Discussion

Male recombination at the Eβ hotspot

This work describes the application of high‐resolution sperm typing to analyse meiotic recombination in the class II Eβ gene in four mouse hybrid crosses. Our data confirm the presence of a recombination hotspot in the second intron of the gene, centred on an MT element and only 1.0–1.6 kb wide, and within the interval determined by prior mapping of recombinant inbred lines (Steinmetz et al., 1982, 1987; Bryda et al., 1992). This hotspot is flanked by DNA that is almost totally suppressed for recombination. Male crossovers are symmetrically distributed across the hotspot and occur at a rate strongly influenced by MHC haplotype. Sperm recombination rates for the s × b, s × d and s × k crosses (0.2–1.6 × 10−3 per sperm) are comparable to crossover rates determined from pedigrees (1.0 × 10−3 across the Eβ region; Passmore et al., 1987; Shiroishi et al., 1995), establishing that sperm analysis yields reliable quantitative data on male recombination.

Sequence determinants of recombination rates

It has previously been suggested that recombination rates across mouse MHC hotspots, as in other species, can be reduced by sequence divergence between the inter acting haplotypes (Sant‘Angelo et al., 1992; Yoshino et al., 1995). Datta et al. (1997) found that even a single mismatch over a 350 bp sequence between other wise identical sequences in yeast can inhibit mitotic recombination up to 4‐fold; this anti‐recombination was exerted through the action of mismatch repair proteins. Our data are consistent with the hypothesis that sequence divergence results in decreased recombination efficiency. The degree of divergence between the s haplotype and the other haplotypes over the central 200–400 bp of the hotspot is least for haplotype k (0.25%) and by far the greatest for haplotype p (2.5%) (Table III). The order of divergence from haplotype s over the centre of the hotspot (k < d < b << p) inversely correlates with the relative rates of recombination in the heterozygotes (k > d > b >> p), consistent with the hypothesis that recombination resulting in crossover is favoured when nucleotide identity is at its greatest at the Eβ hotspot (Bryda et al., 1992; Sant'Angelo et al., 1992). However, it is possible that more distal flanking polymorphisms, as well as genetic background, could influence hotspot activity. For example, cis‐acting elements have been shown to be involved in controlling recombination rates and sex‐specificity at the mouse Psmb9 hotspot (Shiroishi et al., 1991; Heine et al., 1994).

View this table:
Table 3. Number of sequence differences between haplotype s and each of the non‐s haplotypes, around the centre of the Eβ hotspot or over the whole second intron

It has also been suggested that recombination may be influenced by an imperfect four‐base repeat sequence (AGGC)n located just upstream of the third exon (8130 bp) and outside the hotspot (Kobori et al., 1986). This sequence shows similarity to the lambda phage chi sequence and 80% sequence identity to the human minisatellite core sequence (Jeffreys et al., 1985), and has been postulated to provide a recombination signal with the number of repeats governing the activity of the hotspot (Zimmerer and Passmore, 1991). Haplotype s has 11 copies of this repeat sequence, compared with 18 in haplotype d, 10 in b and k, and six in haplotype p. There is no obvious correlation between relative recombination rates (Table I) and either the absolute number of repeats or the difference in repeat copy number compared with haplotype s. There is therefore no clear evidence that this tetramer plays a role in regulating recombination activity.

Reciprocal crossover asymmetry and meiotic drive

Previous evidence of reciprocal crossover asymmetry between the k and b haplotypes at the Eβ hotspot was based on the mapping of only seven recombinants and was therefore very speculative (Zimmerer and Passmore, 1991). The large numbers of sperm crossovers analysed in the present study confirm that asymmetry in the location of reciprocal A and B crossover breakpoints does occur in all hybrids carrying haplotype s. It remains to be seen whether other heterozygous crosses will show a similar phenomenon. The asymmetry results in non‐s alleles being preferentially over‐transmitted to crossover products in all heterozygotes, with transmission distortion being strongest for markers closest to the centre of the hotspot (Figure 4). This phenomenon is remarkably similar to the asymmetry seen in one of the human MHC crossover hotspots (Jeffreys and Neumann, 2002). The degree of transmission distortion in human sperm crossover products (up to 87:13) is very similar to that seen in the Eβ hotspot (up to 89:11), as is the degree of displacement of orientation A and B crossovers (on average 400 bp in Eβ, compared with 430 bp in the human hotspot).

Reciprocal crossover asymmetry and transmission distortion can be readily explained by the double‐strand break (DSB) repair model of recombination (Szostak et al., 1983; Pâques and Haber, 1999; Jeffreys and Neumann, 2002). This model is supported by previous data at human hotspots (Jeffreys et al., 2000, 2001; Jeffreys and Neumann, 2002) and at the Psmb9 hotspot in mice (Guillon and De Massy, 2002) which show that crossovers are approximately symmetrical and normally distributed across a hotspot, compatible with the presence of initiation site(s) located at or close to the centre of the hotspot. Initiation creates a DSB, which then undergoes resection to create 3′ overhanging segments. The 3′ overhangs invade the undamaged chromosome where mismatch removal from the invading strand may occur (Alani et al., 1994), and repair is carried out using information from the undamaged (donor) chromosome (Jeffreys and Neumann, 2002). This process results in sites of crossover resolution occurring at distances from the initiation site determined, in part at least, by the extent of resection, and contributes to the observed width of hotspots.

If crossovers are initiated on both chromosomes with equal frequency, then breakpoints will be identically distributed in orientation A and B crossovers. However, if initiation rates vary between haplotypes, then asymmetry and transmission distortion will arise, increasing in strength with the degree of initiation suppression on one of the haplotypes (Jeffreys and Neumann, 2002). If a hybrid carries a high and low efficiency haplotype, then the majority of crossovers will arise from initiations on the active haplotype. DSB repair will therefore result in the replacement of markers from the initiating chromosome with those from the non‐initiating donor chromosome, resulting in biased gene conversion in favour of the donor chromosome. The consistent over‐transmission of non‐s alleles in all heterozygous mice therefore suggests that haplotype s is most active in initiation. The relative strengths of transmission distortion in different crosses (Table I; Figure 4) further suggest that the initiation efficiency on different haplotypes decreases in the order s > k > (b, d) > p. The extreme over‐transmission of the non‐s allele (up to 89:11 in the s × p cross), combined with the inverse correlation between recombination rate and the strength of over‐transmission, provides strong evidence in favour of this initiation model. Alternative models involving biased mismatch repair of heteroduplex DNA following strand invasion and subsequent branch migration are less plausible (for a fuller discussion see Jeffreys and Neumann, 2002).

There are, however, differences in the properties of reciprocal crossover asymmetry in the Eβ hotspot and in the human MHC hotspot. In the latter, the trigger for asymmetry has been traced to a single SNP at the centre of the hotspot that, when heterozygous, appears to be sufficient to cause asymmetry, presumably by directly influencing the efficiency of DSB formation (Jeffreys and Neumann, 2002). One SNP heterozygosity inside the Eβ hotspot is shared by all four mouse hybrids (marker 3.1, located 100 bp 3′ of the average hotspot centre, with allele C shared by haplotypes b, d, k and p, and the alternative allele T on haplotype s), and the C allele is over‐transmitted in every cross. However, without recombinant haplotype information, it is impossible to determine whether this SNP is triggering asymmetry. Given the variation in transmission distortion between different crosses, it is unlikely that this one SNP is exclusively responsible for asymmetry. Other hotspot SNPs might influence initiation in different heterozygotes, although this seems unlikely since first, they would always have to favour the s haplotype, and secondly, most SNPs within human MHC hotspots do not cause asymmetry (Jeffreys and Neumann, 2002). Instead, it seems likely that more distal elements influence the efficiency of recombination initiation in the Eβ hotspot.

The second difference between the mouse and human hotspots is revealed by recombination rates. In the human hotspot, heterozygosity for the suppressing variant causes asymmetry but has little influence on crossover rate (in theory, the recombination rate in a heterozygote with one totally suppressed chromosome should be 50% of that in a person homozygous for the initiation‐proficient variant) (Jeffreys and Neumann, 2002). In mice, recombination rates can be strongly suppressed in hybrids, particularly in s × p heterozygotes. This implies that in such hybrids, initiating events may be suppressed in trans on the s haplotype or, if they occur, fail to form a productive recombination complex. Perhaps sequence mismatches between haplotypes result in single strands from the s haplotype, produced by gap resection, that either fail to invade the other chromosome or are rejected on invasion. Subsequent DSB recombination repair may then occur preferentially from the sister chromatid, resulting in a sister chromatid exchange. Alternatively, mismatched invasions might be preferentially resolved to yield gene conversions rather than crossovers. The decrease of sperm crossover rates with increasing divergence from the s haplotype is consistent with these suggestions.

Exchanges involving haplotype p

Recombination has never before been detected between p and the other haplotypes at the Eβ hotspot (Lafuse and David, 1986; Shiroishi et al., 1993). The present data indicate that this hotspot does function in s × p heterozygotes, but at a very low efficiency, apparently via recombination initiation occurring preferentially, if not exclusively, on the s haplotype. Yoshino et al. (1995) suggested that insertions and deletions (in their case at a tandem repeat locus at the Psmb9 hotspot) may play a significant role in influencing recombination. In this respect, haplotype p is unique in containing a 597 bp ancestral DNA segment in the 5′ half of the second intron, outside the hotspot, that is replaced by a 1634 bp retroposon element in all other haplotypes (Zimmerer and Passmore, 1991). It is not clear how this insertion/deletion might influence Eβ recombination. It has been suggested that this retroposon element creates a recombination initiation activator on non‐p haplotypes (Zimmerer and Passmore, 1991). This would be very difficult to test, given that it would require analysis of p/non‐p recombinant haplotypes either carrying or devoid of the retroposon.

Similarities between mammalian recombination hotspots

Several mammalian hotspots have been characterized using pedigree approaches, and more recently by high resolution sperm typing. The Eβ hotspot shares remarkably similar features with human hotspots (Jeffreys et al., 1998, 2000, 2001), including the same width, the extreme suppression of recombination in the immediate flanking DNA, the simple structure of the vast majority of sperm crossovers and the evidence from crossover asymmetry for highly localized sites of initiation at the centres of hotspots. These shared features suggest that very similar recombination mechanisms operate in humans and in mice. Over‐transmission of markers from recombination‐suppressed haplotypes into crossover progeny also seems to be a common phenomenon in hotspots, having now been seen in one of six human MHC hotspots (Jeffreys and Neumann, 2002), in a human minisatellite‐associated hotspot (Jeffreys et al., 1998) and in the mouse Eβ hotspot. The resulting meiotic drive of markers from recombinationally suppressed chromosomes, which can include markers that actually suppress recombination, can be sufficiently strong to influence the population frequencies of alleles and to favour the eventual extinction of hotspot activity (Boulton et al., 1997; Jeffreys and Neumann, 2002). The prevalence and evolutionary significance of this recombination‐based meiotic drive will require further high‐resolution studies of meiotic crossover in a range of species.

Materials and methods

Mouse strains and breeding

Mice from strains A.SW, B10.S/SgmedJ (Eβ haplotype s), C57BL/10J (Eβ haplotype b), B10.D2 (Eβ haplotype d), B10.A/SgSnJ (Eβ haplotype k) and P/J (Eβ haplotype p) were purchased from the Jackson Laboratory (Bar Harbor, ME) (haplotype and strain information available on the Jackson Laboratory website: http://www.jax.org). All relevant sequence information has been published previously (Saito et al., 1983; Widera and Flavell, 1984; Kobori et al., 1986; Padgett et al., 1991; Zimmerer and Passmore, 1991; Bryda et al., 1992). Mice with MHC haplotypes b, d, k and p were crossed with haplotype s mice. Caudal epidydimus and brains were collected from 8‐ to 10‐week‐old mature offspring. All crosses were generated at the University of Leicester under guidance issued by the MRC in ‘Responsibility in the use of animals for medical research’ and Home Office project license No. PPL 80/1353.

Genotyping and SNP verification

We used PCR primers designed from the sequenced clone BAC 135G15 (DDBJ/EMBL/GenBank accession No. AF050157) to identify and verify SNPs in the Eβ gene of all strains in the study. We re‐sequenced using BigDye Terminators (ABI) on an ABI 377 Automated Sequencer and identified SNPs using ABI AutoAssembler software. Selected SNPs were used to map recombination breakpoints in subsequent sperm analyses. The ASO probes used for genotyping are shown in Table II.

DNA extraction and quantification

Sperm DNA was prepared using procedures similar to those described by Jeffreys et al. (1994). To minimize the risk of contamination, all manipulations were carried out in a laminar flow hood. Caudal epidydimus from one mouse was chopped in 1 ml of phosphate‐buffered saline (138 mM NaCl, 2.7 mM KCl, 5.4 mM Na2HPO4, 0.2 mM KH2PO4), filtered and centrifuged at 10 000 g. The cell pellet was resuspended in 1× SSC and somatic cells were lysed at room temperature by the addition of SDS to 0.15%. The lysate was centrifuged and the resulting sperm pellet resuspended in 1 ml 0.2× SSC, 0.1% SDS, 1 M 2‐mercaptoethanol and digested with 200 μg/ml proteinase K (Sigma) for 1 h at 37°C. DNA was recovered by phenol/chloroform extraction and ethanol precipitation. DNA from brain was prepared by tissue homogenization in 1× SSC, followed by SDS lysis, phenol/chloroform extraction and ethanol precipitation. All DNAs were dissolved in 5 mM Tris–HCl (pH 7.5). Approximately 20 μg of DNA was digested with enzymes (HpaI, BssHII or XhoI) that cleave outside the test interval analysed, to render genomic DNA fully soluble prior to dilution. DNA concentrations were determined by UV spectrometry and by agarose gel electrophoresis against λ phage DNA.

Sperm crossover analyses

We designed allele‐specific primers 17–20 nt long with at least 50% GC content for SNP sites 5′ and 3′ of the putative Eβ hotspot to amplify across a maximum of 7.7 kb. 5′ primers were located between 3227 and 4301 bp on the reference clone sequence, and 3′ primers between 9855 and 10952 bp. The 5′ allele‐specific primers, named on the basis of distance from the first ASO typed (assigned as 0.0), were:

−1.0p, 5′‐GGT GGG GAG GAT CTG AAA G‐3′, used for haplotypes d and p;

−0.6g, 5′‐TAA CTG TCC AGC CTG GGG‐3′, used for haplotypes b, d and k;

−0.56A, 5′‐TGC CCT GTT AGT TGT GGA‐3′, used for haplotypes b and k;

−0.56T, 5′‐TGC CCC GTT AGT TGT GGT‐3′, used for haplotype s;

−0.57T, 5′‐TTG GAG CTG GGG CAT GT‐3′, used for haplotype s;

−0.5C, 5′‐GGG TCT GGT CGG TCA TTC‐3′, used for haplotype s;

−0.5T, 5′‐CTT GGG TCT GGT CAT TT‐3′, used for haplotype p;

−0.2C, 5′‐GGA GCA AAA GCG GGC CGC‐3′, used for haplotype s;

0.2F, 5′‐CCT ACC TAC ACG GTG TGC GG‐3′, used for haplotype s.

The 3′ primers were:

5.8T, 5′‐CTG ACT GCT TCA GGC CTC T‐3′, used for haplotype s;

5.8g, 5′‐CTG ACT GCT TCA GGC CTC G‐3′, used for haplotypes p and d;

6.75C, 5′‐GGC CCA CAG CAA CAG CTT C‐3′, used for haplotype s;

6.75T, 5′‐GGC CCA CAG CAA CAT CTT T‐3′, used for haplotypes b and k;

6.85C, 5′‐CTC TTC CCA TCC TGC TGC‐3′, used for haplotypes b, d, k and p;

6.85T, 5′‐CTC TTC CCA TCC TGC TGT‐3′, used for haplotype s.

The optimal annealing temperature for each primer was determined by amplification of inbred parental DNA samples from each strain, using the allele‐specific primer in combination with an internal universal primer that can amplify from all haplotypes. Primers were selected that were both very efficient and highly specific for the allele of interest. Primers specific to the 5′ end of haplotype s were used in combination with primers specific to the 3′ end of the non‐s haplotype to amplify orientation A crossovers, and vice versa for orientation B.

We amplified recombinant molecules using procedures similar to those described previously (Jeffreys et al., 2000, 2001). We amplified multiple batches of restriction‐digested and diluted DNA samples, from sperm and brain of a single mouse per cross, each containing on average between 0.4 and 2 amplifiable crossover molecules per reaction. We tested two heterozygous mice for the b and p crosses, and one mouse each for the other two crosses. Since data for mice from the same cross were not significantly different (data not shown), such data were pooled. Brain DNA was used as a negative control for PCR artefacts and/or bleed‐through amplification of one or other progenitor haplotype. DNA inputs depended on crossover rates and were determined in pilot experiments using 300–30 000 amplifiable molecules of each haplotype per reaction; we assume a single molecule PCR efficiency of 85% (Jeffreys et al., 2000). Poisson analysis of the number of positive and negative reactions (Jeffreys et al., 2000, 2001) was used to estimate the number of amplifiable crossover molecules. In the full crossover analyses, inputs varied from 216 ng sperm DNA per reaction (30 000 molecules) for s × p hybrids down to 2 ng DNA per reaction (280 molecules per reaction) for s × k crosses. DNA was amplified in two rounds of nested repulsion‐phase PCR to selectively amplify recombinant molecules (Jeffreys et al., 1998, 2000) (Figure 1). Touch‐down PCRs were used [typically 96°C for 1 min, followed by 26 rounds of 96°C for 20 s, 60–56°C for 30 s (step‐down) and 63°C for 7–8 min (depending on length of amplicon)] with two allele‐specific primers in repulsion phase. Primary PCR products were digested with S1 nuclease to remove single‐stranded DNA and PCR artefacts (Jeffreys et al., 1998) and re‐amplified with a second round of repulsion‐phase allele‐specific PCR as above. Secondary PCR products were electrophoresed on an agarose gel and visualized with ethidium bromide and UV light. Positive reactions were re‐amplified using internal universal primers. PCR products were dot blotted and SNPs typed by ASO hybridization (Wood et al., 1985) as modified by Jeffreys et al. (2000) using the probes listed in Table II. Poisson estimations of the number of amplifiable recombinant molecules, to correct for instances of multiple recombinants recovered in a single positive reaction, were performed as described previously (Jeffreys et al., 1998).

Acknowledgements

We thank Dr John Stead, Dr Celia May, Liisa Kauppi and four anonymous reviewers for helpful comments, and Rita Neumann for help with ASO typing. This work was supported by grants to A.J.J. from the Royal Society and the Wellcome Trust (ref. 058084/Z/99/Z).

References

View Abstract