Human recombination hot spots hidden in regions of strong marker association

The fine-scale distribution of meiotic recombination events in the human genome can be inferred from patterns of haplotype diversity in human populations but directly studied only by high-resolution sperm typing. Both approaches indicate that crossovers are heavily clustered into narrow recombination hot spots. But our direct understanding of hot-spot properties and distributions is largely limited to sperm typing in the major histocompatibility complex (MHC). We now describe the analysis of an unremarkable 206-kb region on human chromosome 1, which identified localized regions of linkage disequilibrium breakdown that mark the locations of sperm crossover hot spots. The distribution, intensity and morphology of these hot spots are markedly similar to those in the MHC. But we also accidentally detected additional hot spots in regions of strong association. Coalescent analysis of genotype data detected most of the hot spots but showed significant differences between sperm crossover frequencies and historical recombination rates. This raises the possibility that some hot spots, particularly those in regions of strong association, may have evolved very recently and not left their full imprint on haplotype diversity. These results suggest that hot spots could be very abundant and possibly fluid features of the human genome.     

Intense and highly localized gene conversion activity in human meiotic crossover hot spots

Meiotic gene conversion has an important role in allele diversification and in the homogenization of gene and other repeat DNA sequence families, sometimes with pathological consequences. But little is known about the dynamics of gene conversion in humans and its relationship to meiotic crossover. We therefore developed screening and selection methods to characterize sperm conversions in two meiotic crossover hot spots in the major histocompatibility complex (MHC) and one in the sex chromosomal pseudoautosomal pairing region PAR1 (ref. 9). All three hot spots are active in gene conversion and crossover. Conversion tracts are short and define a steep bidirectional gradient centered at the peak of crossover activity, consistent with crossovers and conversions being produced by the same recombination-initiating events. These initiations seem to be spread over a narrow zone, rather than occurring at a single site, and seem preferentially to yield conversions rather than crossovers. Crossover breakpoints are more broadly diffused than conversion breakpoints, suggesting either differences between conversion and crossover processing after initiation or the existence of a quality control checkpoint at which short interactions between homologous chromosomes are preferentially aborted as conversions.     

Reciprocal crossover asymmetry and meiotic drive in a human recombination hot spot

Human DNA diversity arises ultimately from germline mutation that creates new haplotypes that can be reshuffled by meiotic recombination. Reciprocal crossover generates recombinant haplotypes but should not influence the frequencies of alleles in a population. We demonstrate crossover asymmetry at a recombination hot spot in the major histocompatibility complex, whereby reciprocal exchanges in sperm map to different locations in the hot spot. We identify a single-nucleotide polymorphism at the center of the hot spot and show that, when heterozygous, it seems sufficient to cause this asymmetry, apparently by influencing the efficiency of highly localized crossover initiation. As a consequence, crossovers in heterozygotes are accompanied by biased gene conversion, most likely occurring by gap repair, that can also affect nearby polymorphisms through repair of an extended gap. The result is substantial over-transmission of the recombination-suppressing allele and neighboring markers to crossover products. Computer simulations show that this meiotic drive, although weak at the population level, is sufficient to favor eventual fixation of the recombination-suppressing variant. These findings provide an explanation for the relatively uniform widths of human crossover hot spots and suggest that hot spots may be generally prone to extinction by meiotic drive.     

Crossover clustering and rapid decay of linkage disequilibrium in the Xp/Yp pseudoautosomal gene SHOX

Crossover between the human sex chromosomes during male meiosis is restricted to the terminal pseudoautosomal pairing regions. An obligatory exchange occurs in PAR1, an Xp/Yp pseudoautosomal region of 2.6 Mb, which creates a male-specific recombination 'hot domain' with a recombination rate that is about 20 times higher than the genome average. Low-resolution analysis of PAR1 suggests that crossovers are distributed fairly randomly. By contrast, linkage disequilibrium (LD) and sperm crossover analyses indicate that crossovers in autosomal regions tend to cluster into 'hot spots' of 1-2 kb that lie between islands of disequilibrium of tens to hundreds of kilobases. To determine whether at high resolution this autosomal pattern also applies to PAR1, we have examined linkage disequilibrium over an interval of 43 kb around the gene SHOX. Here we show that in northern European populations, disequilibrium decays rapidly with physical distance, which is consistent with this interval of PAR1 being recombinationally active in male meiosis. Analysis of a subregion of 9.9 kb in sperm shows, however, that crossovers are not distributed randomly, but instead cluster into an intense recombination hot spot that is very similar in morphology to autosomal hot spots. Thus, PAR1 crossover activity may be influenced by male-specific hot spots that are highly suitable for characterization by sperm DNA analysis.     

An initiation site for meiotic crossing-over and gene conversion in the mouse

During meiosis, the reductional segregation of homologous chromosomes at the first meiotic division requires reciprocal exchange (crossing over) between homologs. The number of crossovers is tightly regulated (one to two per homolog in mice), and their distribution in the genome is not random-recombination 'hot' and 'cold' regions can be identified. We developed a molecular assay to study these events directly in mouse germ cells. This analysis was developed with reference to the proteosome subunit beta type 9 (Psmb9, previously called Lmp2) hot-spot region identified through genetic analysis. Here we show that this hot spot is an initiation site of meiotic recombination on the basis of two observations: (i) crossover density is maximal in an interval of 210 bp and decreases on both sides of this region; (ii) a high frequency of gene conversion is found in the region of highest crossover density. We then used this strategy to carry out the first temporal analysis of meiotic recombination in mouse spermatogenesis and demonstrate that crossover events occur during the pachytene stage of meiotic prophase.     

Evidence for substantial fine-scale variation in recombination rates across the human genome

Characterizing fine-scale variation in human recombination rates is important, both to deepen understanding of the recombination process and to aid the design of disease association studies. Current genetic maps show that rates vary on a megabase scale, but studying finer-scale variation using pedigrees is difficult. Sperm-typing experiments have characterized regions where crossovers cluster into 1-2-kb hot spots, but technical difficulties limit the number of studies. An alternative is to use population variation to infer fine-scale characteristics of the recombination process. Several surveys reported 'block-like' patterns of diversity, which may reflect fine-scale recombination rate variation, but limitations of available methods made this impossible to assess. Here, we applied a new statistical method, which overcomes these limitations, to infer patterns of fine-scale recombination rate variation in 74 genes. We found extensive rate variation both within and among genes. In particular, recombination hot spots are a common feature of the human genome: 47% (35 of 74) of genes showed substantive evidence for a hot spot, and many more showed evidence for some rate variation. No primary sequence characteristics are consistently associated with precise hot-spot location, although G+C content and nucleotide diversity are correlated with local recombination rate.     

The roles of MAD1, MAD2 and MAD3 in meiotic progression and the segregation of nonexchange chromosomes

Errors in meiotic chromosome segregation are the leading cause of spontaneous abortions and birth defects. In humans, chromosomes that fail to experience crossovers (or exchanges) are error-prone, more likely than exchange chromosomes to mis-segregate in meiosis. We used a yeast model to investigate the mechanisms that partition nonexchange chromosomes. These studies showed that the spindle checkpoint genes MAD1, MAD2 and MAD3 have different roles. We identified a new meiotic role for MAD3; though dispensable for the segregation of exchange chromosomes, it is essential for the segregation of nonexchange chromosomes. This function of Mad3p could also be carried out by human BubR1. MAD1 and MAD2 act in a surveillance mechanism that mediates a metaphase delay in response to nonexchange chromosomes, whereas MAD3 acts as a crucial meiotic timer, mediating a prophase delay in every meiosis. These findings suggest plausible models for the basis of errant meiotic segregation in humans.     

A high-resolution HLA and SNP haplotype map for disease association studies in the extended human MHC

The proteins encoded by the classical HLA class I and class II genes in the major histocompatibility complex (MHC) are highly polymorphic and are essential in self versus non-self immune recognition. HLA variation is a crucial determinant of transplant rejection and susceptibility to a large number of infectious and autoimmune diseases. Yet identification of causal variants is problematic owing to linkage disequilibrium that extends across multiple HLA and non-HLA genes in the MHC. We therefore set out to characterize the linkage disequilibrium patterns between the highly polymorphic HLA genes and background variation by typing the classical HLA genes and >7,500 common SNPs and deletion-insertion polymorphisms across four population samples. The analysis provides informative tag SNPs that capture much of the common variation in the MHC region and that could be used in disease association studies, and it provides new insight into the evolutionary dynamics and ancestral origins of the HLA loci and their haplotypes.     

Crossover assurance and crossover interference are distinctly regulated by the ZMM proteins during yeast meiosis

Meiotic crossing-over is highly regulated such that each homolog pair typically receives at least one crossover (assurance) and adjacent crossovers are widely spaced (interference). Here we provide evidence that interference and assurance are mechanistically distinct processes that are separated by mutations in a new ZMM (Zip, Msh, Mer) protein from Saccharomyces cerevisiae, Spo16. Like other zmm mutants, spo16 cells have defects in both crossing-over and synaptonemal complex formation. Unlike in previously characterized zmm mutants, the residual crossovers in spo16 cells show interference comparable to that in the wild type. Spo16 interacts with a second ZMM protein, Spo22 (also known as Zip4), and spo22 mutants also show normal interference. Notably, assembly of the MutS homologs Msh4 and Msh5 on chromosomes occurs in both spo16 and spo22, but not in other zmm mutants. We suggest that crossover interference requires the normal assembly of recombination complexes containing Msh4 and Msh5 but does not require Spo16- and Spo22-dependent extension of synaptonemal complexes. In contrast, crossover assurance requires all ZMM proteins and full-length synaptonemal complexes.     

Genetic variation in the 5q31 cytokine gene cluster confers susceptibility to Crohn disease

Linkage disequilibrium (LD) mapping provides a powerful method for fine-structure localization of rare disease genes, but has not yet been widely applied to common disease. We sought to design a systematic approach for LD mapping and apply it to the localization of a gene (IBD5) conferring susceptibility to Crohn disease. The key issues are: (i) to detect a significant LD signal (ii) to rigorously bound the critical region and (iii) to identify the causal genetic variant within this region. We previously mapped the IBD5 locus to a large region spanning 18 cM of chromosome 5q31 (P<10(-4)). Using dense genetic maps of microsatellite markers and single-nucleotide polymorphisms (SNPs) across the entire region, we found strong evidence of LD. We bound the region to a common haplotype spanning 250 kb that shows strong association with the disease (P< 2 x 10(-7)) and contains the cytokine gene cluster. This finding provides overwhelming evidence that a specific common haplotype of the cytokine region in 5q31 confers susceptibility to Crohn disease. However, genetic evidence alone is not sufficient to identify the causal mutation within this region, as strong LD across the region results in multiple SNPs having equivalent genetic evidence-each consistent with the expected properties of the IBD5 locus. These results have important implications for Crohn disease in particular and LD mapping in general.     

The spindle checkpoint rescues the meiotic segregation of chromosomes whose crossovers are far from the centromere

Improper meiotic chromosome segregation causes conditions such as Down's syndrome. Recombination promotes proper chromosome segregation in meiosis I; chromosomes without crossovers near the centromere are more likely to segregate to the same spindle pole (nondisjoin). Here we have used budding yeast to determine whether the spindle checkpoint promotes segregation of such chromosomes. In checkpoint-defective mad2Delta cells, properly segregating chromosomes have more crossovers near the centromere than their wild-type counterparts, and an artificial tether that holds chromosomes together suppresses nondisjunction as long as the tether is near the centromere. The tether partially rescues the segregation of chromosomes that lack crossovers.     

Five amino acids in three HLA proteins explain most of the association between MHC and seropositive rheumatoid arthritis

The genetic association of the major histocompatibility complex (MHC) to rheumatoid arthritis risk has commonly been attributed to alleles in HLA-DRB1. However, debate persists about the identity of the causal variants in HLA-DRB1 and the presence of independent effects elsewhere in the MHC. Using existing genome-wide SNP data in 5,018 individuals with seropositive rheumatoid arthritis (cases) and 14,974 unaffected controls, we imputed and tested classical alleles and amino acid polymorphisms in HLA-A, HLA-B, HLA-C, HLA-DPA1, HLA-DPB1, HLA-DQA1, HLA-DQB1 and HLA-DRB1, as well as 3,117 SNPs across the MHC. Conditional and haplotype analyses identified that three amino acid positions (11, 71 and 74) in HLA-DRβ1 and single-amino-acid polymorphisms in HLA-B (at position 9) and HLA-DPβ1 (at position 9), which are all located in peptide-binding grooves, almost completely explain the MHC association to rheumatoid arthritis risk. This study shows how imputation of functional variation from large reference panels can help fine map association signals in the MHC.     

Juxtaposed regions of extensive and minimal linkage disequilibrium in human Xq25 and Xq28

Linkage disequilibrium (LD), or the non-random association of alleles, is poorly understood in the human genome. Population genetic theory suggests that LD is determined by the age of the markers, population history, recombination rate, selection and genetic drift. Despite the uncertainties in determining the relative contributions of these factors, some groups have argued that LD is a simple function of distance between markers. Disease-gene mapping studies and a simulation study gave differing predictions on the degree of LD in isolated and general populations. In view of the discrepancies between theory and experimental observations, we constructed a high-density SNP map of the Xq25-Xq28 region and analysed the male genotypes and haplotypes across this region for LD in three populations. The populations included an outbred European sample (CEPH males) and isolated population samples from Finland and Sardinia. We found two extended regions of strong LD bracketed by regions with no evidence for LD in all three samples. Haplotype analysis showed a paucity of haplotypes in regions of strong LD. Our results suggest that, in this region of the X chromosome, LD is not a monotonic function of the distance between markers, but is more a property of the particular location in the human genome.     

Islands of linkage disequilibrium

A detailed knowledge of patterns of linkage disequilibrium in human populations is widely seen as a prerequisite for effective population-based disease gene mapping. New data suggest that linkage disequilibrium is highly structured into discrete blocks of sequence separated by hot spots of recombination.     

A predominant role for the HLA class II region in the association of the MHC region with multiple sclerosis

Genetic susceptibility to multiple sclerosis is associated with genes of the major histocompatibility complex (MHC), particularly HLA-DRB1 and HLA-DQB1 (ref. 1). Both locus and allelic heterogeneity have been reported in this genomic region. To clarify whether HLA-DRB1 itself, nearby genes in the region encoding the MHC or combinations of these loci underlie susceptibility to multiple sclerosis, we genotyped 1,185 Canadian and Finnish families with multiple sclerosis (n = 4,203 individuals) with a high-density SNP panel spanning the genes encoding the MHC and flanking genomic regions. Strong associations in Canadian and Finnish samples were observed with blocks in the HLA class II genomic region (P < 4.9 x 10(-13) and P < 2.0 x 10(-16), respectively), but the strongest association was with HLA-DRB1 (P < 4.4 x 10(-17)). Conditioning on either HLA-DRB1 or the most significant HLA class II haplotype block found no additional block or SNP association independent of the HLA class II genomic region. This study therefore indicates that MHC-associated susceptibility to multiple sclerosis is determined by HLA class II alleles, their interactions and closely neighboring variants.     

A humanized model for multiple sclerosis using HLA-DR2 and a human T-cell receptor

Multiple sclerosis (MS) is a complex chronic neurologic disease with a suspected autoimmune pathogenesis. Although there is evidence that the development of MS is determined by both environmental influences and genes, these factors are largely undefined, except for major histocompatibility (MHC) genes. Linkage analyses and association studies have shown that susceptibility to MS is associated with genes in the human histocompatibility leukocyte antigens (HLA) class II region, but the contribution of these genes to MS disease development less compared with their contribution to disorders such as insulin-dependent diabetes mellitus. Due to the strong linkage disequilibrium in the MHC class II region, it has not been possible to determine which gene(s) is responsible for the genetic predisposition. In transgenic mice, we have expressed three human components involved in T-cell recognition of an MS-relevant autoantigen presented by the HLA-DR2 molecule: DRA*0101/DRB1*1501 (HLA-DR2), an MHC class II candidate MS susceptibility genes found in individuals of European descent; a T-cell receptor (TCR) from an MS-patient-derived T-cell clone specific for the HLA-DR2 bound immunodominant myelin basic protein (MBP) 4102 peptide; and the human CD4 coreceptor. The amino acid sequence of the MBP 84-102 peptide is the same in both human and mouse MBP. Following administration of the MBP peptide, together with adjuvant and pertussis toxin, transgenic mice developed focal CNS inflammation and demyelination that led to clinical manifestations and disease courses resembling those seen in MS. Spontaneous disease was observed in 4% of mice. When DR2 and TCR double-transgenic mice were backcrossed twice to Rag2 (for recombination-activating gene 2)-deficient mice, the incidence of spontaneous disease increased, demonstrating that T cells specific for the HLA-DR2 bound MBP peptide are sufficient and necessary for development of disease. Our study provides evidence that HLA-DR2 can mediate both induced and spontaneous disease resembling MS by presenting an MBP self-peptide to T cells.     

Chromosome-wide distribution of haplotype blocks and the role of recombination hot spots

Recent studies of human populations suggest that the genome consists of chromosome segments that are ancestrally conserved ('haplotype blocks'; refs. 1-3) and have discrete boundaries defined by recombination hot spots. Using publicly available genetic markers, we have constructed a first-generation haplotype map of chromosome 19. As expected for this marker density, approximately one-third of the chromosome is encompassed within haplotype blocks. Evolutionary modeling of the data indicates that recombination hot spots are not required to explain most of the observed blocks, providing that marker ascertainment and the observed marker spacing are considered. In contrast, several long blocks are inconsistent with our evolutionary models, and different mechanisms could explain their origins.     

From HLA association to function

A new study refines the association signals for rheumatoid arthritis susceptibility in the major histocompatibility complex (MHC) region to five amino-acid positions encoded in three HLA genes, all within peptide-binding grooves. By adapting statistical methods from genome-wide association studies (GWAS) and using imputation from a large reference panel, they demonstrate the potential for this approach to identify functional variants in associated regions.     

On the subspecific origin of the laboratory mouse

The genome of the laboratory mouse is thought to be a mosaic of regions with distinct subspecific origins. We have developed a high-resolution map of the origin of the laboratory mouse by generating 25,400 phylogenetic trees at 100-kb intervals spanning the genome. On average, 92% of the genome is of Mus musculus domesticus origin, and the distribution of diversity is markedly nonrandom among the chromosomes. There are large regions of extremely low diversity, which represent blind spots for studies of natural variation and complex traits, and hot spots of diversity. In contrast with the mosaic model, we found that most of the genome has intermediate levels of variation of intrasubspecific origin. Finally, mouse strains derived from the wild that are supposed to represent different mouse subspecies show substantial intersubspecific introgression, which has strong implications for evolutionary studies that assume these are pure representatives of a given subspecies.