High resolution localization of recombination hot spots using sperm typing

Genome-wide recombination map construction from single sperm sequencing in cattle

BACKGROUND

Meiotic recombination is one of the important phenomena contributing to gamete genome diversity. However, except for human and a few model organisms, it is not well studied in livestock, including cattle.

RESULTS

To investigate their distributions in the cattle sperm genome, we sequenced 143 single sperms from two Holstein bulls. We mapped meiotic recombination events at high resolution based on phased heterozygous single nucleotide polymorphism (SNP). In the absence of evolutionary selection pressure in fertilization and survival, recombination events in sperm are enriched near distal chromosomal ends, revealing that such a pattern is intrinsic to the molecular mechanism of meiosis. Furthermore, we further validated these findings in single sperms with results derived from sequencing its family trio of diploid genomes and our previous studies of recombination in cattle.

CONCLUSIONS

To our knowledge, this is the first large-scale single sperm whole-genome sequencing effort in livestock, which provided useful information for future studies of recombination, genome instability, and male infertility.

Population-Specific Recombination Maps from Segments of Identity by Descent

Recombination rates vary significantly across the genome, and estimates of recombination rates are needed for downstream analyses such as haplotype phasing and genotype imputation. Existing methods for recombination rate estimation are limited by insufficient amounts of informative genetic data or by high computational cost. We present a method and software, called IBDrecomb, for using segments of identity by descent to infer recombination rates. IBDrecomb can be applied to sequenced population cohorts to obtain high-resolution, population-specific recombination maps. In simulated admixed data, IBDrecomb obtains higher accuracy than admixture-based estimation of recombination rates. When applied to 2,500 simulated individuals, IBDrecomb obtains similar accuracy to a linkage-disequilibrium (LD)-based method applied to 96 individuals (the largest number for which computation is tractable). Compared to LD-based maps, our IBD-based maps have the advantage of estimating recombination rates in the recent past rather than the distant past. We used IBDrecomb to generate new recombination maps for European Americans and for African Americans from TOPMed sequence data from the Framingham Heart Study (1,626 unrelated individuals) and the Jackson Heart Study (2,046 unrelated individuals), and we compare them to LD-based, admixture-based, and family-based maps.

Fast Estimation of Recombination Rates Using Topological Data Analysis

Accurate estimation of recombination rates is critical for studying the origins and maintenance of genetic diversity. Because the inference of recombination rates under a full evolutionary model is computationally expensive, we developed an alternative approach using topological data analysis (TDA) on genome sequences. We find that this method can analyze datasets larger than what can be handled by any existing recombination inference software, and has accuracy comparable to commonly used model-based methods with significantly less processing time. Previous TDA methods used information contained solely in the first Betti number ([Formula: see text]) of a set of genomes, which aims to capture the number of loops that can be detected within a genealogy. These explorations have proven difficult to connect to the theory of the underlying biological process of recombination, and, consequently, have unpredictable behavior under perturbations of the data. We introduce a new topological feature, which we call ψ, with a natural connection to coalescent models, and present novel arguments relating [Formula: see text] to population genetic models. Using simulations, we show that ψ and [Formula: see text] are differentially affected by missing data, and package our approach as TREE (Topological Recombination Estimator). TREE's efficiency and accuracy make it well suited as a first-pass estimator of recombination rate heterogeneity or hotspots throughout the genome. Our work empirically and theoretically justifies the use of topological statistics as summaries of genome sequences and describes a new, unintuitive relationship between topological features of the distribution of sequence data and the footprint of recombination on genomes.

Construction of PRDM9 allele-specific recombination maps in cattle using large-scale pedigree analysis and genome-wide single sperm genomics

PRDM9 contributes to hybrid sterility and species evolution. However, its role is to be confirmed in cattle, a major domesticated livestock species. We previously found an association near PRDM9 with cattle recombination features, but the causative variants are still unknown. Using millions of genotyped cattle with pedigree information, we characterized five PRDM9 alleles and generated allele-specific recombination maps. By examining allele-specific recombination patterns, we observed the impact of PRDM9 on global distribution of recombination, especially in the two ends of chromosomes. We also showed strong associations between recombination hotspot regions and functional mutations within PRDM9 zinc finger domain. More importantly, we found one allele of PRDM9 to be very different from others in both protein composition and recombination landscape, indicating the causative role of this allele on the association between PRDM9 and cattle recombination. When comparing recombination maps from sperm and pedigree data, we observed similar genome-wide recombination patterns, validating the quality of pedigree-based results. Collectively, these evidence supported PRDM9 alleles as causal variants for the reported association with cattle recombination. Our study comprehensively surveyed the bovine PRDM9 alleles, generated allele-specific recombination maps, and expanded our understanding of the role of PRDM9 on genome distribution of recombination.

Effects of Demographic History on the Detection of Recombination Hotspots from Linkage Disequilibrium

In some species, meiotic recombination is concentrated in small genomic regions. These "recombination hotspots" leave signatures in fine-scale patterns of linkage disequilibrium, raising the prospect that the genomic landscape of hotspots can be characterized from sequence variation. This approach has led to the inference that hotspots evolve rapidly in some species, but are conserved in others. Historic demographic events, such as population bottlenecks, are known to affect patterns of linkage disequilibrium across the genome, violating population genetic assumptions of this approach. Although such events are prevalent, demographic history is generally ignored when making inferences about the evolution of recombination hotspots. To determine the effect of demography on the detection of recombination hotspots, we use the coalescent to simulate haplotypes with a known recombination landscape. We measure the ability of popular linkage disequilibrium-based programs to detect hotspots across a range of demographic histories, including population bottlenecks, hidden population structure, population expansions, and population contractions. We find that demographic events have the potential to greatly reduce the power and increase the false positive rate of hotspot discovery. Neither the power nor the false positive rate of hotspot detection can be predicted without also knowing the demographic history of the sample. Our results suggest that ignoring demographic history likely overestimates the power to detect hotspots and therefore underestimates the degree of hotspot sharing between species. We suggest strategies for incorporating demographic history into population genetic inferences about recombination hotspots.

Recombination hotspots: Models and tools for detection

Recombination hotspots are the regions within the genome where the rate, and the frequency of recombination are optimum with a size varying from 1 to 2kb. The recombination event is mediated by the double-stranded break formation, guided by the combined enzymatic action of DNA topoisomerase and Spo 11 endonuclease. These regions are distributed non-uniformly throughout the human genome and cause distortions in the genetic map. Numerous lines of evidence suggest that the number of hotspots known in humans has increased manifold in recent years. A few facts about the hotspot evolutions were also put forward, indicating the differences in the hotspot position between chimpanzees and humans. In mice, recombination hot spots were found to be clustered within the major histocompatibility complex (MHC) region. Several models, that help explain meiotic recombination has been proposed. Moreover, scientists also developed some computational tools to locate the hotspot position and estimate their recombination rate in humans is of great interest to population and medical geneticists. Here we reviewed the molecular mechanisms, models and in silico prediction techniques of hot spot residues.

ARG-walker: inference of individual specific strengths of meiotic recombination hotspots by population genomics analysis

BACKGROUND

Meiotic recombination hotspots play important roles in various aspects of genomics, but the underlying mechanisms for regulating the locations and strengths of recombination hotspots are not yet fully revealed. Most existing algorithms for estimating recombination rates from sequence polymorphism data can only output average recombination rates of a population, although there is evidence for the heterogeneity in recombination rates among individuals. For genome-wide association studies (GWAS) of recombination hotspots, an efficient algorithm that estimates the individualized strengths of recombination hotspots is highly desirable.

RESULTS

In this work, we propose a novel graph mining algorithm named ARG-walker, based on random walks on ancestral recombination graphs (ARG), to estimate individual-specific recombination hotspot strengths. Extensive simulations demonstrate that ARG-walker is able to distinguish the hot allele of a recombination hotspot from the cold allele. Integrated with output of ARG-walker, we performed GWAS on the phased haplotype data of the 22 autosome chromosomes of the HapMap Asian population samples of Chinese and Japanese (JPT+CHB). Significant cis-regulatory signals have been detected, which is corroborated by the enrichment of the well-known 13-mer motif CCNCCNTNNCCNC of PRDM9 protein. Moreover, two new DNA motifs have been identified in the flanking regions of the significantly associated SNPs (single nucleotide polymorphisms), which are likely to be new cis-regulatory elements of meiotic recombination hotspots of the human genome.

CONCLUSIONS

Our results on both simulated and real data suggest that ARG-walker is a promising new method for estimating the individual recombination variations. In the future, it could be used to uncover the mechanisms of recombination regulation and human diseases related with recombination hotspots.

The 8p23 inversion polymorphism determines local recombination heterogeneity across human populations

For decades, chromosomal inversions have been regarded as fascinating evolutionary elements as they are expected to suppress recombination between chromosomes with opposite orientations, leading to the accumulation of genetic differences between the two configurations over time. Here, making use of publicly available population genotype data for the largest polymorphic inversion in the human genome (8p23-inv), we assessed whether this inhibitory effect of inversion rearrangements led to significant differences in the recombination landscape of two homologous DNA segments, with opposite orientation. Our analysis revealed that the accumulation of genetic differentiation is positively correlated with the variation in recombination profiles. The observed recombination dissimilarity between inversion types is consistent across all populations analyzed and surpasses the effects of geographic structure, suggesting that both structures (orientations) have been evolving independently over an extended period of time, despite being subjected to the very same demographic history. Aside this mainly independent evolution, we also identified a short segment (350 kb, <10% of the whole inversion) in the central region of the inversion where the genetic divergence between the two structural haplotypes is diminished. Although it is difficult to demonstrate it, this could be due to gene flow (possibly via double-crossing over events), which is consistent with the higher recombination rates surrounding this segment. This study demonstrates for the first time that chromosomal inversions influence the recombination landscape at a fine-scale and highlights the role of these rearrangements as drivers of genome evolution.

Reliable and Fast Estimation of Recombination Rates by Convergence Diagnosis and Parallel Markov Chain Monte Carlo

Genetic recombination is an essential event during the process of meiosis resulting in an exchange of segments between paired chromosomes. Estimating recombination rate is crucial for understanding the process of recombination. Experimental methods are normally difficult and limited to small scale estimations. Thus statistical methods using population genetics data are important for large-scale analysis. LDhat is an extensively used statistical method using rjMCMC algorithm to predict recombination rates. Due to the complexity of rjMCMC scheme, LDhat may take a long time for large SNP data sets. In addition, rjMCMC parameters should be manually defined in the original program which directly impact results. To address these issues, we designed an improved algorithm based on LDhat implementing MCMC convergence diagnostic algorithms to automatically predict values of parameters and monitor the mixing process. Then parallel computation methods were employed to further accelerate the new program. The new algorithms have been tested on ten samples from HapMap phase 2 data set. The results were compared with previous code and showed nearly identical output. However, our new methods achieved significant acceleration proving that they are more efficient and reliable for the estimation of recombination rates. The stand-alone package is freely available for download http://www.ntu.edu.sg/home/zhengjie/software/CPLDhat.

High-resolution linkage map for two honeybee chromosomes: the hotspot quest

Meiotic recombination is a fundamental process ensuring proper disjunction of homologous chromosomes and allele shuffling in successive generations. In many species, this cellular mechanism occurs heterogeneously along chromosomes and mostly concentrates in tiny fragments called recombination hotspots. Specific DNA motifs have been shown to initiate recombination in these hotspots in mammals, fission yeast and drosophila. The aim of this study was to check whether recombination also occurs in a heterogeneous fashion in the highly recombinogenic honeybee genome and whether this heterogeneity can be connected with specific DNA motifs. We completed a previous picture drawn from a routine genetic map built with an average resolution of 93 kb. We focused on the two smallest honeybee chromosomes to increase the resolution and even zoomed at very high resolution (3.6 kb) on a fragment of 300 kb. Recombination rates measured in these fragments were placed in relation with occurrence of 30 previously described motifs through a Poisson regression model. A selection procedure suitable for correlated variables was applied to keep significant motifs. These fine and ultra-fine mappings show that recombination rate is significantly heterogeneous although poorly contrasted between high and low recombination rate, contrarily to most model species. We show that recombination rate is probably associated with the DNA methylation state. Moreover, three motifs (CGCA, GCCGC and CCAAT) are good candidates of signals promoting recombination. Their influence is however moderate, doubling at most the recombination rate. This discovery extends the way to recombination dissection in insects.

Inferring genome-wide recombination landscapes from advanced intercross lines: application to yeast crosses

Accurate estimates of recombination rates are of great importance for understanding evolution. In an experimental genetic cross, recombination breaks apart and rejoins genetic material, such that the genomes of the resulting isolates are comprised of distinct blocks of differing parental origin. We here describe a method exploiting this fact to infer genome-wide recombination profiles from sequenced isolates from an advanced intercross line (AIL). We verified the accuracy of the method against simulated data. Next, we sequenced 192 isolates from a twelve-generation cross between West African and North American yeast Saccharomyces cerevisiae strains and inferred the underlying recombination landscape at a fine genomic resolution (mean segregating site distance 0.22 kb). Comparison was made with landscapes inferred for a similar cross between four yeast strains, and with a previous single-generation, intra-strain cross (Mancera et al., Nature 2008). Moderate congruence was identified between landscapes (correlation 0.58-0.77 at 5 kb resolution), albeit with variance between mean genome-wide recombination rates. The multiple generations of mating undergone in the AILs gave more precise inference of recombination rates than could be achieved from a single-generation cross, in particular in identifying recombination cold-spots. The recombination landscapes we describe have particular utility; both AILs are part of a resource to study complex yeast traits (see e.g. Parts et al., Genome Res 2011). Our results will enable future applications of this resource to take better account of local linkage structure heterogeneities. Our method has general applicability to other crossing experiments, including a variety of experimental designs.

Sequencing of isolated sperm cells for direct haplotyping of a human genome

There is increasing evidence that the phenotypic effects of genomic sequence variants are best understood in terms of variant haplotypes rather than as isolated polymorphisms. Haplotype analysis is also critically important for uncovering population histories and for the study of evolutionary genetics. Although the sequencing of individual human genomes to reveal personal collections of sequence variants is now well established, there has been slower progress in the phasing of these variants into pairs of haplotypes along each pair of chromosomes. Here, we have developed a distinct approach to haplotyping that can yield chromosome-length haplotypes, including the vast majority of heterozygous single-nucleotide polymorphisms (SNPs) in an individual human genome. This approach exploits the haploid nature of sperm cells and employs a combination of genotyping and low-coverage sequencing on a short-read platform. In addition to generating chromosome-length haplotypes, the approach can directly identify recombination events (averaging 1.1 per chromosome) with a median resolution of <100 kb.

Analysis of recombinants in female mouse meiosis

During meiosis, homologous chromosomes (homologs) undergo recombinational interactions, resulting in the formation of crossovers (COs) or noncrossovers (NCOs). Both COs and NCOs are initiated by the same event: programmed double-strand DNA breaks (DSBs), which occur preferentially at hotspots throughout the genome. COs contribute to the genetic diversity of gametes and are needed to promote proper meiotic chromosome segregation. Accordingly, their formation is tightly controlled. In the mouse, the sites of preferred CO formation differ between male and female chromosomes, both on a regional level and on the level of individual hotspots. Sperm typing using (half-sided) allele-specific PCR has proven a powerful technique to characterize COs and all detectable NCOs at hotspots on male human and mouse chromosomes. In contrast, very little is known about the properties of hotspots in female meiosis. This chapter describes an adaptation of sperm typing to analyze recombinants in a hotspot, using DNA isolated from an ovary cell suspension enriched for oocytes.

Isolation of meiotic recombinants from mouse sperm

Homologous recombination during meiosis is critical for the formation of gametes. Recombination is initiated by programmed DNA double-strand breaks which preferentially occur at hotspots dispersed throughout the genome. These double-strand breaks are repaired from the homolog, resulting in either a crossover or noncrossover product. Multiple noncrossover events are required for homolog pairing, and at least one crossover is critical for proper chromosome segregation at the first meiotic division. Consequently, homologous recombination in meiosis occurs at high frequencies. This chapter describes how to characterize crossovers and noncrossovers at a hotspot in mice using allele-specific PCR. Amplification of recombinant products directly from sperm DNA is a powerful approach to determine recombination frequencies and map recombination breakpoints, providing insight into homologous recombination mechanisms.

Detecting sequence polymorphisms associated with meiotic recombination hotspots in the human genome

Background

Meiotic recombination events tend to cluster into narrow spans of a few kilobases long, called recombination hotspots. Such hotspots are not conserved between human and chimpanzee and vary between different human ethnic groups. At the same time, recombination hotspots are heritable. Previous studies showed instances where differences in recombination rate could be associated with sequence polymorphisms.

Results

In this work we developed a novel computational approach, LDsplit, to perform a large-scale association study of recombination hotspots with genetic polymorphisms. LDsplit was able to correctly predict the association between the FG11 SNP and the DNA2 hotspot observed by sperm typing. Extensive simulation demonstrated the accuracy of LDsplit under various conditions. Applying LDsplit to human chromosome 6, we found that for a significant fraction of hotspots, there is an association between variations in intensity of historical recombination and sequence polymorphisms. From flanking regions of the SNPs output by LDsplit we identified a conserved 11-mer motif GGNGGNAGGGG, whose complement partially matches 13-mer CCNCCNTNNCCNC, a critical motif for the regulation of recombination hotspots.

Conclusions

Our result suggests that computational approaches based on historical recombination events are likely to be more powerful than previously anticipated. The putative associations we identified may be a promising step toward uncovering the mechanisms of recombination hotspots.

Contrasting methods of quantifying fine structure of human recombination

There has been considerable excitement over the ability to construct linkage maps based only on genome-wide genotype data for single nucleotide polymorphic sites (SNPs) in a population sample. These maps, which are derived from estimates of linkage disequilibrium (LD), rely on population genetics theory to relate the decay of LD to the local rate of recombination, but other population processes also come into play. Here we contrast these LD maps to the classically derived, pedigree-based human recombination maps. The LD maps have a level of resolution greatly exceeding that of the pedigree maps, and at this fine scale, sperm typing allows a means of validation. While at a gross level both the pedigree maps and the sperm typing methods generally agree with LD maps, there are significant local differences between them, and the fact that these maps measure different genetic features should be remembered when using them for other genetic inferences.

Genetic crossovers are predicted accurately by the computed human recombination map

Hotspots of meiotic recombination can change rapidly over time. This instability and the reported high level of inter-individual variation in meiotic recombination puts in question the accuracy of the calculated hotspot map, which is based on the summation of past genetic crossovers. To estimate the accuracy of the computed recombination rate map, we have mapped genetic crossovers to a median resolution of 70 Kb in 10 CEPH pedigrees. We then compared the positions of crossovers with the hotspots computed from HapMap data and performed extensive computer simulations to compare the observed distributions of crossovers with the distributions expected from the calculated recombination rate maps. Here we show that a population-averaged hotspot map computed from linkage disequilibrium data predicts well present-day genetic crossovers. We find that computed hotspot maps accurately estimate both the strength and the position of meiotic hotspots. An in-depth examination of not-predicted crossovers shows that they are preferentially located in regions where hotspots are found in other populations. In summary, we find that by combining several computed population-specific maps we can capture the variation in individual hotspots to generate a hotspot map that can predict almost all present-day genetic crossovers.

In eukaryotes genetic crossovers are responsible for generating genetic diversity and ensuring the proper segregation of chromosomes. Genetic crossovers are tightly clustered in hotspots. Although the existence of hotspots in humans is clearly proven, mechanisms of their formation and the regulation of meiotic recombination in general remain poorly understood. An additional complication in studies of meiotic recombination is the fact that the direct experimental mapping of human hotspots on a genome-wide scale is not feasible with current methods. The best available indirect methods compute the position of hotspots from patterns of historic associations between genetic markers in population samples. In this study we determined the positions of genetic crossovers in ten pedigrees of European origin and then compared the positions of crossovers with the hotspots computed from HapMap data. Importantly, we find that the population-averaged computed map is in close agreement with the observed distribution of genetic crossovers. We also find that cryptic hotspots that are not easily detected in the computed European map can be more effectively identified if other populations are included in the analysis. Our analysis shows that high-resolution recombination profiles are highly similar between distantly related populations and that by including computed hotspots from several populations we can predict nearly all crossovers.

Product length, dye choice, and detection chemistry in the bead-emulsion amplification of millions of single DNA molecules in parallel

The amplification of millions of single molecules in parallel can be performed on microscopic magnetic beads that are contained in aqueous compartments of an oil-buffer emulsion. These bead-emulsion amplification (BEA) reactions result in beads that are covered by almost-identical copies derived from a single template. The post-amplification analysis is performed using different fluorophore-labeled probes. We have identified BEA reaction conditions that efficiently produce longer amplicons of up to 450 base pairs. These conditions include the use of a Titanium Taq amplification system. Second, we explored alternate fluorophores coupled to probes for post-PCR DNA analysis. We demonstrate that four different Alexa fluorophores can be used simultaneously with extremely low crosstalk. Finally, we developed an allele-specific extension chemistry that is based on Alexa dyes to query individual nucleotides of the amplified material that is both highly efficient and specific.

Fraction of informative recombinations: a heuristic approach to analyze recombination rates

In this article we present a new heuristic approach (informative recombinations, InfRec) to analyze recombination density at the sequence level. InfRec is intuitive and easy and combines previously developed methods that (i) resolve genotypes into haplotypes, (ii) estimate the minimum number of recombinations, and (iii) evaluate the fraction of informative recombinations. We tested this approach in its sliding-window version on 117 genes from the SeattleSNPs program, resequenced in 24 African-Americans (AAs) and 23 European-Americans (EAs). We obtained population recombination rate estimates (rho(obs)) of 0.85 and 0.37 kb(-1) in AAs and EAs, respectively. Coalescence simulations indicated that these values account for both the recombinations and the gene conversions in the history of the sample. The intensity of rho(obs) varied considerably along the sequence, revealing the presence of recombination hotspots. Overall, we observed approximately 80% of recombinations in one-third and approximately 50% in only 10% of the sequence. InfRec performance, tested on published simulated and additional experimental data sets, was similar to that of other hotspot detection methods. Fast, intuitive, and visual, InfRec is not constrained by sample size limitations. It facilitates understanding data and provides a simple and flexible tool to analyze recombination intensity along the sequence.

Mammalian meiotic recombination hot spots

Our understanding of the details of mammalian meiotic recombination has recently advanced significantly. Sperm typing technologies, linkage studies, and computational inferences from population genetic data have together provided information in unprecedented detail about the location and activity of the sites of crossing-over in mice and humans. The results show that the vast majority of meiotic recombination events are localized to narrow DNA regions (hot spots) that constitute only a small fraction of the genome. The data also suggest that the molecular basis of hot spot activity is unlikely to be strictly determined by specific DNA sequence motifs in cis. Further molecular studies are needed to understand how hot spots originate, function and evolve.

Playing hide and seek with mammalian meiotic crossover hotspots

Crossovers (COs) are essential for meiosis and contribute to genome diversity by promoting the reassociation of alleles, and thus improve the efficiency of selection. COs are not randomly distributed but are found at specific regions, or CO hotspots. Recent results have revealed the historical recombination rates and the distribution of hotspots across the human genome. Surprisingly, CO hotspots are highly dynamic, as shown by differences in activity between individuals, populations and closely related species. We propose a role for DNA methylation in preventing the formation of COs, a regulation that might explain, in part, the correlation between recombination rates and GC content in mammals.

Combining sperm typing and linkage disequilibrium analyses reveals differences in selective pressures or recombination rates across human populations

A previous polymorphism survey of the type 2 diabetes gene CAPN10 identified a segment showing an excess of polymorphism levels in all population samples, coinciding with localized breakdown of linkage disequilibrium (LD) in a sample of Hausa from Cameroon, but not in non-African samples. This raised the possibility that a recombination hotspot is present in all populations and we had insufficient power to detect it in the non-African data. To test this possibility, we estimated the crossover rate by sperm typing in five non-African men; these estimates were consistent with the LD decay in the non-African, but not in the Hausa data. Moreover, resequencing the orthologous region in a sample of Western chimpanzees did not show either an excess of polymorphism level or rapid LD decay, suggesting that the processes underlying the patterns observed in humans operated only on the human lineage. These results suggest that a hotspot of recombination has recently arisen in humans and has reached higher frequency in the Hausa than in non-Africans, or that there is no elevation in crossover rate in any human population, and the observed variation results from long-standing balancing selection.

Molecular features of meiotic recombination hot spots

Meiotic recombination occurs preferentially at certain regions called hot spots and is important for generating genetic diversity and proper segregation of chromosomes during meiosis. Hot spots have been characterized most extensively in yeast, mice and humans. The development of methods based on sperm typing and population genetics has facilitated rapid and high-resolution mapping of hot spots in mice and humans in recent years. With increasing information becoming available on meiotic recombination in different species, it is now possible to compare several molecular features associated with hot-spot loci. Further, there have been advances in our knowledge of the factors influencing hot-spot activity and the role that they play in structuring the genome into haplotype blocks. We review the molecular features associated with hot spots in terms of their properties and mechanisms underlying their function and distribution. A large number of these features seem to be shared among hot spots from different species suggesting common mechanisms for their formation and function.

Positive selection can create false hotspots of recombination

Simulations of positive directional selection, under parameter values appropriate for approximating human genetic diversity and rates of recombination, reveal that the effects of strong selective sweeps on patterns of linkage disequilibrium (LD) mimic the pattern expected with recombinant hotspots.

Singleton births after routine preimplantation genetic diagnosis using exclusion testing (D4S43 and D4S126) for Huntington’s disease

OBJECTIVE

To develop exclusion testing protocols for Huntington's disease (HD) linkage markers suitable for use in a clinical preimplantation genetic diagnosis (PGD) setting for couples in whom a partner was at 50% risk of inheriting HD, but who choose not to undergo presymptomatic mutation testing.

DESIGN

Preimplantation genetic diagnosis using exclusion testing.

SETTING

In vitro fertilization clinic.

PATIENT(S)

Three couples with family histories of HD, two couples opposed to direct mutation testing.

INTERVENTION(S)

Development of single-cell polymerase chain reaction tests for PGD for the HD mutation and two HD gene-flanking markers (D4S43 and D4S126), allowing the identification of an individual embryo as being at either low or high risk for developing HD without being diagnostic of the presence of the mutation.

MAIN OUTCOME MEASURE(S)

D4S43, D4S126, and HD mutation.

RESULT(S)

After PGD for HD, couple 1 gave birth to a healthy girl after a frozen embryo transfer, and genetic status was confirmed by prenatal diagnosis to be very low risk for developing HD. Couple 2 gave birth to a healthy boy after their second cycle of PGD, and couple 3, after a third cycle, gave birth to a boy with congenital heart defects, which were successfully corrected with surgery at age 5 days. Both couples 2 and 3 declined prenatal testing, and therefore relinquished the opportunity to confirm the PGD.

CONCLUSION(S)

Preimplantation genetic diagnosis for HD using exclusion testing resulted in three live singleton births after six oocyte recovery procedures. The diagnostic protocol provided couples the opportunity to minimize the likelihood of disease transmission to their children, without the requirement for predictive testing.

HUMHOT: a database of human meiotic recombination hot spots

Meiotic recombination occurs preferentially at certain regions in the genome referred to as hot spots. The number of hot spots known in humans has increased manifold in recent years. The identification of these hot spots in humans is of great interest to population and medical geneticists since they influence the structure of Linkage Disequilibrium and Haplotype blocks in human populations, whose patterns have applications in mapping disease genes. HUMHOT is a web-based database of Human Meiotic Recombination Hot Spots. The database comprises DNA sequences corresponding to the hot spot regions from the literature that have been mapped to a high resolution (<4 kb) in humans. It also provides flanking sequence information for the hot spot region along with references describing the hot spot. The database can be queried based on hot spot identity, chromosome position or by homology to user-defined sequences. It is also updated with new hot spot sequences as they are discovered and provides hyperlinks to commonly used tools for estimating recombination rates, performing genetic analysis and new advances in our understanding of meiotic hot spots. Public access to the HUMHOT database is available at .

Allelic recombination and de novo deletions in sperm in the human beta-globin gene region

Meiotic recombination is of fundamental importance in creating haplotype diversity in the human genome and has the potential to cause genomic rearrangements by ectopic recombination between repeat sequences and through other changes triggered by recombination-initiating events. However, the relationship between allelic recombination and genome instability in the human germline remains unclear. We have therefore analysed recombination and DNA instability in the delta-, beta-globin gene region and its associated recombination hotspot. Sperm typing has for the first time accurately defined the hotspot and shown it to be the most active autosomal crossover hotspot yet described, although unusually inactive in non-exchange gene conversion. The hotspot just extends into a homology block shared by the delta- and beta-globin genes, within which ectopic exchanges can generate Hb Lepore deletions. We developed a physical selection method for recovering and validating extremely rare de novo deletions in human DNA and used it to characterize the dynamics of these Hb Lepore deletions in sperm as well as other deletions not arising from ectopic exchanges between homologous DNA sequences. Surprisingly, both classes of deletion showed breakpoints that avoided the beta-globin hotspot, establishing that it possesses remarkable fidelity and does not play a significant role in triggering these DNA rearrangements. This study also provides the first direct analysis of de novo deletion in the human germline and points to a possible deletion-controlling element in the beta-globin gene separate from the crossover hotspot.

Genome amplification of single sperm using multiple displacement amplification

Sperm typing is an effective way to study recombination rate on a fine scale in regions of interest. There are two strategies for the amplification of single meiotic recombinants: repulsion-phase allele-specific PCR and whole genome amplification (WGA). The former can selectively amplify single recombinant molecules from a batch of sperm but is not scalable for high-throughput operation. Currently, primer extension pre-amplification is the only method used in WGA of single sperm, whereas it has limited capacity to produce high-coverage products enough for the analysis of local recombination rate in multiple large regions. Here, we applied for the first time a recently developed WGA method, multiple displacement amplification (MDA), to amplify single sperm DNA, and demonstrated its great potential for producing high-yield and high-coverage products. In a 50 mul reaction, 76 or 93% of loci can be amplified at least 2500- or 250-fold, respectively, from single sperm DNA, and second-round MDA can further offer >200-fold amplification. The MDA products are usable for a variety of genetic applications, including sequencing and microsatellite marker and single nucleotide polymorphism (SNP) analysis. The use of MDA in single sperm amplification may open a new era for studies on local recombination rates.

Single-nucleotide-polymorphism genotyping for whole-genome-amplified samples using automated fluorescence correlation spectroscopy

Whole-genome amplification (WGA) methods were adopted for single-nucleotide-polymorphism (SNP) typing to minimize the amount of genomic DNA that has to be used in typing for thousands of different SNPs in large-scale studies; 5-10 ng of genomic DNA was amplified by a WGA method (improved primer-extension-preamplification-polymerase chain reaction (I-PEP-PCR), degenerated oligonucleotide primer-PCR (DOP-PCR), or multiple displacement amplification (MDA)). Using 1/100 to 1/500 amounts of the whole-genome-amplified products as templates, subsequent analyses were successfully performed. SNPs were genotyped by the sequence-specific primer (SSP)-PCR method followed by fluorescence correlation spectroscopy (FCS). The typing results were evaluated for four different SNPs on tumor necrosis factor receptor 1 and 2 genes (TNFR1 and TNFR2). The genotypes determined by the SSP-FCS method using the WGA products were 100% in concordance with those determined by nucleotide sequencing using genomic DNAs. We have already carried out typing of more than 300 different SNPs and are currently performing 7,500-10,000 typings per day using WGA samples from patients with several common diseases. WGA coupled with FCS allows specific and high-throughput genotyping of thousands of samples for thousands of different SNPs.

Absence of the TAP2 human recombination hotspot in chimpanzees

Recent experiments using sperm typing have demonstrated that, in several regions of the human genome, recombination does not occur uniformly but instead is concentrated in “hotspots” of 1–2 kb. Moreover, the crossover asymmetry observed in a subset of these has led to the suggestion that hotspots may be short-lived on an evolutionary time scale. To test this possibility, we focused on a region known to contain a recombination hotspot in humans, TAP2, and asked whether chimpanzees, the closest living evolutionary relatives of humans, harbor a hotspot in a similar location. Specifically, we used a new statistical approach to estimate recombination rate variation from patterns of linkage disequilibrium in a sample of 24 western chimpanzees (Pan troglodytes verus). This method has been shown to produce reliable results on simulated data and on human data from the TAP2 region. Strikingly, however, it finds very little support for recombination rate variation at TAP2 in the western chimpanzee data. Moreover, simulations suggest that there should be stronger support if there were a hotspot similar to the one characterized in humans. Thus, it appears that the human TAP2 recombination hotspot is not shared by western chimpanzees. These findings demonstrate that fine-scale recombination rates can change between very closely related species and raise the possibility that rates differ among human populations, with important implications for linkage-disequilibrium based association studies.

The human TAP2 recombination hotspot is absent from the homologous region in western chimpanzees, with important implications for association studies, the HapMap project and understanding fine-scale variation in recombination rates

Justified chauvinism: advances in defining meiotic recombination through sperm typing

Sperm typing offers an efficient means of studying the quantitative and qualitative aspects of meiotic recombination that are virtually unapproachable by pedigree analysis. Since the initial development of the technique >10 years ago, several salient findings based on empirically derived recombination data have been described. The precise rates and distributions of recombination have been reported for specific regions of the genome, serving as the prototype for high-resolution genome-wide recombination patterns. Identification and characterization of molecular genetic events, such as unequal crossing over, gene conversion and crossover asymmetry, are under close inspection for the first time as a result of this technology. The influence of these phenomena on the evolution of the genome is of primary interest from a scientific and medical perspective. In this article, we review the novel discoveries in mammalian meiotic recombination that have been revealed through sperm typing.

Comparative linkage-disequilibrium analysis of the beta-globin hotspot in primates

Recombination rates vary both across the genome and between different species, but little information is available about the temporal and physical scales over which such rates change. To shed light on these questions, we performed a high-resolution analysis of a genomic region within the beta-globin gene cluster that is known to experience elevated recombination rates in humans. For this purpose, we developed new linkage disequilibrium-based methods that thoroughly search for subsets of the data with unusually high or unusually low estimated values of the population-recombination parameter (4Nr, where N is the effective population size and r is the crossover rate between adjacent base pairs). By resequencing a 15-kb segment in a human population sample, we were able to narrow the recombinational hotspot to a segment <2 kb in length that coincides with the beta-globin replication origin. In addition, we analyzed the orthologous region in samples of rhesus macaques and common chimpanzees. Whereas the analysis of the chimpanzee data is complicated by the sample structure, the macaque data imply that this region may not be a hotspot in that species. These results suggest a time scale for the evolution of hotspots in primates. Furthermore, they allow us to propose diverged sequence elements that may contribute to the differences in the recombinational landscape in the two species.

Mapping candidate hotspots of meiotic recombination in segments of human DNA cloned in the yeast Saccharomyces cerevisiae

The hotspots of meiotic recombination in the human genome can be localized by genetic techniques. The resolution of these techniques is in the range of kilobases and depends on the density of the physical markers identifying allelic variants of the chromosomal loci. We thought it would be interesting to localize these sites with higher resolution. Assuming that some human chromosomal sites conserve their propensity for recombination when cloned in yeast, we localized the hotspots of recombination in several yeast artificial chromosomes (YACs) carrying human DNA. A number of potential recombination hotspots could be identified in the clones studied. Among them there are two classes of sites that are particularly recombination prone also in human meiotic cells: sites associated with CpG islands and sites located in the vicinity of long minisatellite sequences.

Haplotype blocks and linkage disequilibrium in the human genome

There is great interest in the patterns and extent of linkage disequilibrium (LD) in humans and other species. Characterizing LD is of central importance for gene-mapping studies and can provide insights into the biology of recombination and human demographic history. Here, we review recent developments in this field, including the recently proposed 'haplotype-block' model of LD. We describe some of the recent data in detail and compare the observed patterns to those seen in simulations.

Linkage disequilibrium analysis of the human adenosine deaminase (ada) gene provides evidence for a lack of correlation between hot spots of equal and unequal homologous recombination

The linkage disequilibrium (LD) pattern within the adenosine deaminase (ADA) gene was analyzed by studying 13 polymorphic loci in 137 families from two European and three African populations. Evidence for the presence of a 12-kb meiotic crossover hot spot, spanning part of the first and the second intron and flanked by regions of reduced recombination activity, was obtained. Moreover, segregation analysis of 113 informative meioses revealed two recombination events that are internal or overlap the 12-kb region, thus suggesting a recombination rate for the hot-spot region about 50-fold higher than the mean rate across the human genome. Within the hot spot, a 144-bp palindromic sequence was also identified and its possible involvement in the recombination process is discussed. The 12-kb region characterized by the low degree of LD does not include the 3.2-kb region that is deleted, as a result of recurrent unequal homologous recombination between two Alu elements, in patients affected by autosomal severe combined immunodeficiency. This observation provides the first evidence for an absence of correlation between hot spots of equal and unequal homologous recombination.

Spline methods for the comparison of physical and genetic maps

The first genetic maps were constructed by linkage analysis. Physical mapping techniques, such as radiation hybrids and complete sequencing, produce a different picture. For the purposes of population genetics, clinical genetics, and genetic epidemiology, it is important to harmonize and amalgamate existing genetic and physical maps. Among other things, comparisons of the two kinds of maps promotes better understanding of the wide variation in local recombination rates per unit physical length of DNA. The current paper presents methods for estimating recombination intensity as a function of physical distance along a chromosome. Genetic map distance is the integral of intensity. We derive fast reliable estimation algorithms based on a Poisson process model, penalized likelihoods, and cubic spline interpolation. Our methods provide a rigorous and statistically sound foundation for comparing physical and genetic maps. To illustrate the possibilities, we apply the methods to published recombination data on CEPH families and the complete sequences of chromosomes 21 and 22. Our results are in good agreement with previous studies and the biological data.

High-resolution patterns of meiotic recombination across the human major histocompatibility complex

Definitive characteristics of meiotic recombination events over large (i.e., >1 Mb) segments of the human genome remain obscure, yet they are essential for establishing the haplotypic structure of the genome and for efficient mapping of complex traits. We present a high-resolution map of recombination at the kilobase level across a 3.3-Mb interval encompassing the major histocompatibility complex (MHC). Genotyping of 20,031 single sperm from 12 individuals resulted in the identification and fine mapping of 325 recombinant chromosomes within genomic intervals as small as 7 kb. Several principal characteristics of recombination in this region were observed: (1) rates of recombination can differ significantly between individuals; (2) intense hot spots of recombination occur at least every 0.8 Mb but are not necessarily evenly spaced; (3) distribution in the location of recombination events can differ significantly among individuals; (4) between hot spots, low levels of recombination occur fairly evenly across 100-kb segments, suggesting the presence of warm spots of recombination; and (5) specific sequence motifs associate significantly with recombination distribution. These data provide a plausible model for recombination patterns of the human genome overall.

Meiotic recombination hot spots and cold spots

Meiotic recombination events are distributed unevenly throughout eukaryotic genomes. This inhomogeneity leads to distortions of genetic maps that can hinder the ability of geneticists to identify genes by map-based techniques. Various lines of evidence, particularly from studies of yeast, indicate that the distribution of recombination events might reflect, at least in part, global features of chromosome structure, such as the distribution of modified nucleosomes.

Bipolar disorder and variation at a common polymorphism (A1832G) within exon 8 of the Wolfram gene

A number of linkage studies provide evidence consistent with the existence of a bipolar susceptibility gene on chromosome 4p16. The gene for Wolfram syndrome, a rare recessive neurodegenerative disorder, lies in this region and has recently been cloned. Psychiatric disturbances including psychosis, mood disorder, and suicide have been reported at increased frequency in Wolfram patients and in heterozygous carriers of a Wolfram mutation. In the current investigation we have undertaken a case-control association study using a single nucleotide polymorphism (causing an amino acid change) in exon 8 of the Wolfram gene in a UK Caucasian sample of 312 Diagnostic and Statistical Manual of Mental Disorders (fourth edition; DSM IV) bipolar I probands and 301 comparison individuals. We found no evidence that variation at this polymorphism influences susceptibility to bipolar disorder. It remains possible that variation at other sites within or near the Wolfram gene plays important roles in determining susceptibility to affective illness. Am. J. Med. Genet. (Neuropsychiatr. Genet.) 96:154-157, 2000.

Evidence for heterogeneity in recombination in the human pseudoautosomal region: high resolution analysis by sperm typing and radiation-hybrid mapping

Accurate genetic and physical maps for the human pseudoautosomal region were constructed by use of sperm typing and high-resolution radiation-hybrid mapping. PCR analysis of 1,912 sperm was done with a manual, single-sperm isolation method. Data on four donors show highly significant linkage heterogeneity among individuals. The most significant difference was observed in a marker interval located in the middle of the Xp/Yp pseudoautosomal region, where one donor showed a particularly high recombination fraction. Longitudinal models were fitted to the data to test whether linkage heterogeneity among donors was significant for multiple intervals across the region. The results indicated that increased recombination in particular individuals and regions is compensated for by reduced recombination in neighboring intervals. To investigate correspondence between physical and genetic distances within the region, we constructed a high-resolution radiation-hybrid map containing 29 markers. The recombination fraction per unit of physical distance varies between regions ranging from 13- to 70-fold greater than the genome-average rate.

Technical advance: single pollen typing combined with laser-mediated manipulation

We combined single pollen typing with laser-mediated manipulation. After drilling a hole in the wall of a pollen grain from a dioecious plant (Silene latifolia) with a UV-laser microbeam, the single pollen grain was recovered directly in the cap of a PCR tube, using a non-contact method called laser pressure catapulting. The entire genome of the single pollen grain was then amplified with improved primer-extension-preamplification PCR (I-PEP PCR). Nested PCR with sequence tagged site (STS)-specific primers was used to analyze several loci in the haploid genome. The single copy gene MROS1 was detected in most of the single pollen grains analyzed. Bgl10, which is localized on the Y chromosome, was detected in approximately half of the pollen grains. MROS3 is reported to be localized on the X chromosome. Using inverse PCR, we isolated two genomic clones of MROS3: MROS3A and MROS3B. The single pollen analysis using nested PCR showed that MROS3A and MROS3B are derived from different loci that are not located on the X chromosome. Single pollen typing not only reveals sex chromosome-linkage within the haploid genome, but can also discriminate between alleles and different loci. This method should also be useful for measuring recombination frequencies without genetic crossover analysis.

Fine-mapping of a region of variation in recombination rate on BTA23 to the D23S22-D23S23 interval using sperm typing and meiotic breakpoint analysis

Meiotic recombination rate (theta) within chromosome segments of similar physical size is known to vary widely throughout the genome. This variation has a genetic component, occurring between the sexes and among individuals of the same sex. We reported previously the existence of variation in theta between males in the DYA-PRL interval on bovine chromosome 23 (BTA23). This region contains the bovine major histocompatibility complex and has been shown to contain recombination hotspots in humans and mice. The aim of this study was to map more finely the interval(s) on BTA23 where variation in theta occurs using sperm typing and meiotic breakpoint analysis. By adding a marker (DRB3) between DYA and PRL, the DYA-PRL interval was subdivided into two adjacent intervals, thus permitting evaluation and comparison of theta among five bulls. Significant variation in theta was found for both intervals; theta(DYA-DRB3) ranged from 13.2 to 28.1%, and theta(DRB3-PRL) ranged from 2.4 to 13.0%. The variation in theta was individual- and region-specific. A meiotic breakpoint strategy employing PCR amplification products from recombinant sperm was then used to refine the chromosomal location associated with variation in theta within the DYA-DRB3 interval. The subinterval D23S22-D23S23 exhibited the greatest degree of variation among bulls having high and low theta within the DYA-DRB3 interval. To confirm this result, theta(D23S22-D23S23) was directly evaluated in three additional randomly chosen bulls using sperm typing. The region showing variation in theta was narrowed to the D23S22-D23S23 subinterval, ranging from 4.6 to 9.2%. Identification of the molecular basis for variation in theta may be useful for map-dependent applications, such as marker-assisted selection and positional cloning of genes affecting physiologically important traits.