Recombination hotspots rather than population history dominate linkage disequilibrium in the MHC class II region

Recombination, demographic history, drift and selection influence the extent of linkage disequilibrium (LD) in the human genome, but their relative contributions remain unclear. To investigate the effect of meiotic recombination versus population history on LD, three populations with different demographic histories (UK north Europeans, Saami and Zimbabweans) were genotyped for high-frequency single-nucleotide polymorphisms (SNPs) across a 75 kb DNA segment of the MHC class II region. This region spans three well-characterized recombination hotspots and a 60 kb long LD block. Despite a high level of underlying haplotype diversity and considerable divergence in haplotype composition between populations, all three populations showed very similar patterns of LD. Surprisingly, the entire 60 kb LD block was present even in Africans, although it was relatively difficult to detect owing to a systematic deficiency of high frequency SNPs. In contrast, DNA within recombination hotspots did not show this low nucleotide diversity in Africans. Thus, while population history has some influence on LD, our findings suggest that recombination hotspots play a major global role in shaping LD patterns as well as helping to maintain localized SNP diversity in this region of the MHC.     

Polymorphism in the activity of human crossover hotspots independent of local DNA sequence variation

Meiotic crossovers in the human genome cluster into highly localized hotspots identifiable indirectly from patterns of DNA diversity and directly by high-resolution sperm typing. Little is known about factors that control hotspot activity and the apparently rapid turnover of hotspots during recent evolution. Clues can, however, be gained by characterizing variation in sperm crossover activity between men. Previous studies have identified single nucleotide polymorphisms within hotspots that appear to suppress crossover activity and which may be involved in hotspot attenuation/extinction. We now analyse a closely spaced pair of hotspots (MSTM1a, MSTM1b) on chromosome 1q42.3, the former being a candidate for a young hotspot that has failed to leave a significant mark on haplotype diversity. Extensive surveys of different men revealed substantial polymorphism in sperm crossover frequencies at both hotspots, but with very different patterns of variation. Hotspot MSTM1b was active in all men tested but with widely differing crossover frequencies. In contrast, MSTM1a was active in only a few men and appeared to be recombinationally inert in the remainder, providing the first example of presence/absence polymorphism of a human hotspot. Haplotype analysis around both hotspots identified active and suppressed men sharing identical haplotypes, establishing that these major variations in the presence/absence of a hotspot and in quantitative activity are not caused by local DNA sequence variation. These findings suggest a role for distal regulators or epigenetic factors in hotspot activity and provide the first direct evidence for the rapid evolution of recombination hotspots in humans.     

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.     

Mapping recombination sites in the HLA class II region.

Studies of linkage disequilibrium across the HLA class II region have been useful in predicting where recombination is most likely to occur. Strong associations exist within the 100kb region from DQB1 to DRB1, suggesting a low frequency of recombination in this region. Conversely, a lack of association between alleles of TAP1 and TAP2 (about 15kb) has been observed, suggesting that recombination occurs here with relative frequency. Analysis of familial-derived recombinant chromosomes is likely to enrich the current knowledge of recombination across the class II region which has been derived from disequilibrium studies of known class II haplotypes. Families containing recombinant chromosomes within the 700kb interval between the DPB1 and DRB1 genes have been identified for such purposes. Using single-strand conformation polymorphism analysis, 122 novel polymorphic markers distributed throughout the class II region were identified and used to delineate the site of crossover events in 26 class II recombinant chromosomes. The three regions where recombination occurred most frequently are as follows: the 45kb interval between DN-A and RING3 in 10 cases, the 40kb interval between DQB3 and DQB1 in 6 cases, and an 8.8kb segment of the TAP2 gene in 3 cases. The final seven recombinant chromosomes require further resolution, but all of these potentially fall within one of the three intervals mentioned. Two of the three TAP2 recombinant chromosomes have been mapped to overlapping 139bp and 850bp segments of intron 2. Associations between polymorphic markers immediately flanking each of the 3 defined regions using data from unrelated individuals will be compared with the familial recombination data observed in these regions.     

A recombination hotspot leads to sequence variability within a novel gene (AK005651) and contributes to type 1 diabetes susceptibility

More than 25 loci have been linked to type 1 diabetes (T1D) in the nonobese diabetic (NOD) mouse, but identification of the underlying genes remains challenging. We describe here the positional cloning of a T1D susceptibility locus, Idd11, located on mouse chromosome 4. Sequence analysis of a series of congenic NOD mouse strains over a critical 6.9-kb interval in these mice and in 25 inbred strains identified several haplotypes, including a unique NOD haplotype, associated with varying levels of T1D susceptibility. Haplotype diversity within this interval between congenic NOD mouse strains was due to a recombination hotspot that generated four crossover breakpoints, including one with a complex conversion tract. The Idd11 haplotype and recombination hotspot are located within a predicted gene of unknown function, which exhibits decreased expression in relevant tissues of NOD mice. Notably, it was the recombination hotspot that aided our mapping of Idd11 and confirms that recombination hotspots can create genetic variation affecting a common polygenic disease. This finding has implications for human genetic association studies, which may be affected by the approximately 33,000 estimated hotspots in the genome.     

Recombination within the human MHC

Summary: Recombination (crossing over) in the human MHC is thought to have played a role in generation of novel alleles at various HLA loci. It is also responsible for the diversity observed at the haplotype level, although the functional consequences of this activity are not clear. Historic and family studies of recombination have provided estimations of recombination fractions across the MHC and identified potential hotspots for recombination in the class II region. Other characteristics of recombination in the human MHC such as haplotype specificity in recombination frequency and localized sequence motifs involved in recombination have been considered, but have been difficult to address given the constraints of human population studies, Single-sperm typing holds promise in overcoming some of the limitations inherent in the study of recombination in human populations. Both family-based and sperm typing analyses of recombination, along with our knowledge of linkage disequilibrium patterns in the MHC, may provide novel information regarding the evolution of HLA haplotypes chat will be difficult to obtain by other means.     

Recombination hotspots and block structure of linkage disequilibrium in the human genome exemplified by detailed analysis of PGM1 on 1p31

The distribution of linkage disequilibrium (LD) in the human genome has important consequences for the design of experiments that infer susceptibility genes for complex disease using association studies. Recent studies have shown a non-random distribution of human meiotic recombination associated with intervening tracts of LD. Little is known about the processes, patterns and frequency of reciprocal meiotic recombination in humans. However, this phenomenon can be better understood by the fine structure analysis of several genomic regions by mapping hotspots and characterizing regions with variable LD. Here, we report clustered hotspot activity with intervening blocks of LD within the human PGM1 gene (1p31) using data derived from meiotic and population studies. Earlier work has suggested a high recombination rate in two regions within the PGM1 gene, site A (exons 4-8) and site B (exons 1A-4). Sequencing of eight individuals across 6 kb of targeted regions in site B identified 18 informative SNPs. Individuals from three distinct populations, Caucasian (n=264), Chinese (n=222) and Vietnamese (n=187), were genotyped, and haplotypes were determined using estimate of haplotypes, ldmax and Arlequin. Allelic association and haplotype analysis in these samples revealed variable recombination rates across PGM1, demonstrating the presence of: (i) three hotspots and (ii) three haplotype blocks. The spatial arrangement of haplotype blocks was identical in all populations studied. The pattern of association within PGM1 represents a region decomposed into small blocks of LD, where increased recombination activity has disrupted the ancestral chromosome. Additionally, crossovers in phased data mapped preferentially to regions where LD collapses, which also overlap with sequence motifs.     

Non-hotspot-related breakpoints of common deletions in Sotos syndrome are located within destabilised DNA regions

Background: Sotos syndrome (SoS) is a disorder characterised by excessive growth, typical craniofacial features, and developmental retardation. It is caused by haploinsuffiency of NSD1 at 5q35. There is a 3.0 kb recombination hotspot in which the breakpoints of around 80% of SoS patients with a common deletion can be mapped. Objective: To identify deletion breakpoints located outside the SoS recombination hotspot. Methods: A screening system for the directly orientated segments of the SoS LCRs was developed for 10 SoS patients with a common deletion who were negative for the SoS hotspot. Deletion-junction fragments were analysed for DNA duplex stability and their relation to scaffold/matrix attachment regions (S/MARs). These features were compared with the SoS hotspot and recombination hotspots of other genomic disorders. Results: The breakpoint was mapped in four SoS patients, two with a deletion in the maternally derived chromosome. These breakpoint regions were located ∼2.5 kb, ∼9.6 kb, ∼27.2, and ∼27.7 kb telomeric to the SoS hotspot and were confined to 164 bp, 46 bp, 256 bp, and 124 bp, respectively. Two of the regions were mapped within Alu elements. All crossover events were found to have occurred within or adjacent to a highly destabilised DNA duplex with a high S/MAR probability. In contrast, the SoS hotspot and other genomic disorders'' recombination hotspots were mapped to stabilised DNA helix regions, flanked by destabilised regions with high probability of containing S/MAR elements. Conclusions: The data suggest that a specific chromatin structure may increase susceptibility for recurrent crossover events and thus predispose to recombination hotspots in genomic disorders.     

Sequence features in regions of weak and strong linkage disequilibrium

We use genotype data generated by the International HapMap Project to dissect the relationship between sequence features and the degree of linkage disequilibrium in the genome. We show that variation in linkage disequilibrium is broadly similar across populations and examine sequence landscape in regions of strong and weak disequilibrium. Linkage disequilibrium is generally low within approximately 15 Mb of the telomeres of each chromosome and noticeably elevated in large, duplicated regions of the genome as well as within approximately 5 Mb of centromeres and other heterochromatic regions. At a broad scale (100-1000 kb resolution), our results show that regions of strong linkage disequilibrium are typically GC poor and have reduced polymorphism. In addition, these regions are enriched for LINE repeats, but have fewer SINE, DNA, and simple repeats than the rest of the genome. At a fine scale, we examine the sequence composition of "hotspots" for the rapid breakdown of linkage disequilibrium and show that they are enriched in SINEs, in simple repeats, and in sequences that are conserved between species. Regions of high and low linkage disequilibrium (the top and bottom quartiles of the genome) have a higher density of genes and coding bases than the rest of the genome. Closer examination of the data shows that whereas some types of genes (including genes involved in immune response and sensory perception) are typically located in regions of low linkage disequilibrium, other genes (including those involved in DNA and RNA metabolism, response to DNA damage, and the cell cycle) are preferentially located in regions of strong linkage disequilibrium. Our results provide a detailed analysis of the relationship between sequence features and linkage disequilibrium and suggest an evolutionary justification for the heterogeneity in linkage disequilibrium in the genome.     

Death of PRDM9 coincides with stabilization of the recombination landscape in the dog genome

Analysis of diverse eukaryotes has revealed that recombination events cluster in discrete genomic locations known as hotspots. In humans, a zinc-finger protein, PRDM9, is believed to initiate recombination in >40% of hotspots by binding to a specific DNA sequence motif. However, the PRDM9 coding sequence is disrupted in the dog genome assembly, raising questions regarding the nature and control of recombination in dogs. By analyzing the sequences of PRDM9 orthologs in a number of dog breeds and several carnivores, we show here that this gene was inactivated early in canid evolution. We next use patterns of linkage disequilibrium using more than 170,000 SNP markers typed in almost 500 dogs to estimate the recombination rates in the dog genome using a coalescent-based approach. Broad-scale recombination rates show good correspondence with an existing linkage-based map. Significant variation in recombination rate is observed on the fine scale, and we are able to detect over 4000 recombination hotspots with high confidence. In contrast to human hotspots, 40% of canine hotspots are characterized by a distinct peak in GC content. A comparative genomic analysis indicates that these peaks are present also as weaker peaks in the panda, suggesting that the hotspots have been continually reinforced by accelerated and strongly GC biased nucleotide substitutions, consistent with the long-term action of biased gene conversion on the dog lineage. These results are consistent with the loss of PRDM9 in canids, resulting in a greater evolutionary stability of recombination hotspots. The genetic determinants of recombination hotspots in the dog genome may thus reflect a fundamental process of relevance to diverse animal species.     

Allelic Association and Recombination Hotspots in the Mucin Gene (MUC) Complex on Chromosome 11p15.5

A family of four highly polymorphic genes encoding secreted gel forming mucins is located in the middle of a recombination rich region of the short arm of chromosome 11 (11p15.5; tel MUC6 - MUC2 - MUC5AC - MUC5B cen; approx. 400 kb). These genes are of interest as risk factors for inflammatory diseases of the epithelia, and we have for example reported association of a VNTR polymorphism of the major mucin domain of MUC2 with asthma, despite the fact that MUC2 is not a major respiratory mucin. To understand the significance of this and other mucin gene associations it is important to describe the patterns of linkage disequilibrium (LD) across this chromosomal region, which is still incomplete on HapMap and the UCSC Golden Path sequence. Our previous studies on the 40 core CEPH families provided direct evidence for several recombination events within the immediate region of the gene complex. This study examines these recombination events in more detail, and also the patterns of LD across the gene complex. We refine the location of the breakpoints, and the combined data suggest two probable recombination hotspots. Three breakpoints are located between MUC6 and MUC2 : there is no association between MUC6 and MUC2 , and the data collected here, combined with that publicly available, maps a hotspot to a region of 4 kb. The other recombinants map between MUC2 and intron 8 of MUC5B. Relatively strong LD is detected between MUC2 and MUC5AC , and although 10/70 of the chromosomes tested shared a common haplotype, which extends from MUC2 to MUC5B, statistically significant association was not detected between MUC2 and the markers tested in MUC5B . We discuss the possibility that the previously reported association between MUC2 and asthma is most likely attributable to association with functional variation in MUC5AC, which encodes one of the two major mucins expressed in both healthy and diseased airways.     

Factors influencing recombination frequency and distribution in a human meiotic crossover hotspot

Little is known about the factors that influence the frequency and distribution of meiotic recombination events within human crossover hotspots. We now describe the detailed analysis of sperm recombination in the NID1 hotspot. Like the neighbouring MS32 hotspot, the NID1 hotspot is associated with a minisatellite, suggesting that hotspots predispose DNA to tandem repetition. Unlike MS32, crossover resolution breakpoints in NID1 avoid the minisatellite, producing a cold spot within the hotspot. This avoidance may be related to the palindromic nature of the minisatellite interfering with the generation and/or processing of recombination intermediates. The NID1 hotspot also contains a single nucleotide polymorphism (SNP) close to the centre, which appears to directly influence the frequency of crossover initiation. Quantitative gene conversion assays show that this SNP affects the frequency of gene conversion and crossover to a very similar extent, providing evidence that conversions and crossovers are triggered by the same recombination initiating events. The recombination-suppressing allele is over-transmitted to recombinant progeny, and provides the most dramatic example to date of recombination-mediated meiotic drive, of a magnitude sufficient to virtually guarantee that the recombination suppressor will eventually replace the more active allele in human populations.     

Consideration of the haplotype diversity at nonallelic homologous recombination hotspots improves the precision of rearrangement breakpoint identification

Precise characterization of nonallelic homologous recombination (NAHR) breakpoints is key to identifying those features that influence NAHR frequency. Until now, analysis of NAHR‐mediated rearrangements has generally been performed by comparison of the breakpoint‐spanning sequences with the human genome reference sequence. We show here that the haplotype diversity of NAHR hotspots may interfere with breakpoint‐mapping. We studied the transmitting parents of individuals with germline type‐1 NF1 deletions mediated by NAHR within the paralogous recombination site 1 (PRS1) or paralogous recombination site 2 (PRS2) hotspots. Several parental wild‐type PRS1 and PRS2 haplotypes were identified that exhibited considerable sequence differences with respect to the reference sequence, which also affected the number of predicted PRDM9‐binding sites. Sequence comparisons between the parental wild‐type PRS1 or PRS2 haplotypes and the deletion breakpoint‐spanning sequences from the patients (method #2) turned out to be an accurate means to assign NF1 deletion breakpoints and proved superior to crude reference sequence comparisons that neglect to consider haplotype diversity (method #1). The mean length of the deletion breakpoint regions assigned by method #2 was 269‐bp in contrast to 502‐bp by method #1. Our findings imply that paralog‐specific haplotype diversity of NAHR hotspots (such as PRS2) and population‐specific haplotype diversity must be taken into account in order to accurately ascertain NAHR‐mediated rearrangement breakpoints.     

Crossover breakpoint mapping identifies a subtelomeric hotspot for male meiotic recombination

Segregation analysis of CEPH and other pedigrees yielded six paternal crossover breakpoints in the approximately 85 kb interval between the minisatellite loci D16S309 (MS205) and D16S83 (EKMDA2) in 16p13.3. Three crossovers were mapped to within the same small (<3 kb) interval, which does not co-localize with any tandem repeat array or expressed sequence. Haplotyping of loci harbouring single nucleotide polymorphism (SNP) markers in this interval confirmed the exchange of flanking markers in the three recombinant individuals. Sequence analysis revealed the presence of recombination-associated motifs and binding sites for the protein translin. Haplotyping of 108 individuals from three European populations at four loci harbouring SNPs showed substantial linkage equilibrium across this interval. Hence molecular and population genetic data are consistent with the presence of an intense male-specific recombination hotspot at this locus.     

Strong linkage disequilibrium of TAP1*0301 and TAP2D alleles with the HLA A1-B35-DRB1*1104-DQA1*0103-DQB1*0603 extended haplotype in Ashkenazi Jews.

On the basis of data collected during the 12th International Histocompatibility Workshop, we postulated a possible linkage disequilibrium between the TAP1C allele and the DRB1*1104-DQA1*0103-DQB1*0603 haplotype characteristic of Ashkenazi Jews. In order to confirm and extend this preliminary observation, a group of 34 subjects carrying this haplotype was analysed for TAP1 and TAP2 polymorphisms and compared with two control groups sharing either the DRB1 or the DQA1 and DQB1 alleles with the test group. The TAP1*0301 and TAP2D alleles were found to be strongly associated with the entire haplotype, but not with the DRB1*1104 or the DQA1*0103-DQB1*0603 alleles when carried separately. These data show a strong linkage disequilibrium of TAP1*0301 and TAP2D alleles with the DRB1*1104-DQA1*0103-DQB1*0603 haplotype of Ashkenazi Jews, extended on the centromeric and telomeric side to the DPB1*0201 and A1-B35 alleles, respectively.     

Compound haplotypes at Xp11.23 and human population growth in Eurasia

To investigate patterns of diversity and the evolutionary history of Eurasians, we have sequenced a 2.8 kb region at Xp11.23 in a sample of African and Eurasian chromosomes. This region is in a long intron of CLCN5 and is immediately flanked by a highly variable minisatellite, DXS255, and a human‐specific Ta0 LINE. Compared to Africans, Eurasians showed a marked reduction in sequence diversity. The main Euro‐Asiatic haplotype seems to be the ancestral haplotype for the whole sample. Coalescent simulations, including recombination and exponential growth, indicate a median length of strong linkage disequilibrium, up to ～9kb for this area. The K_a/K_s ratio between the coding sequence of human CLCN5 and its mouse orthologue is much less than 1. This implies that the region sequenced is unlikely to be under the strong influence of positive selective processes on CLCN5, mutations in which have been associated with disorders such as Dent's disease. In contrast, a scenario based on a population bottleneck and exponential growth seems a more likely explanation for the reduced diversity observed in Eurasians. Coalescent analysis and linked minisatellite diversity (which reaches a gene diversity value greater than 98% in Eurasians) suggest an estimated age of origin of the Euro‐Asiatic diversity compatible with a recent out‐of‐Africa model for colonization of Eurasia by modern Homo sapiens.     

Strong linkage disequilibrium of TAP1*0301 and TAP2D alleles with the HLA A1-B35-DRB1*1104-DQA1*0103-DQB1*0603 extended haplotype in Ashkenazi Jews

On the basis of data collected during the 12th International Histocompatibility Workshop, we postulated a possible linkage disequilibrium between the TAP1C allele and the DRB1*1104-DQA1*0103-DQB1*0603 haplotype characteristic of Ashkenazi Jews. In order to confirm and extend this preliminary observation, a group of 34 subjects carrying this haplotype was analysed for TAP1 and TAP2 polymorphisms and compared with two control groups sharing either the DRB1 or the DQA1 and DQB1 alleles with the test group. The TAP1*0301 and TAP2D alleles were found to be strongly associated with the entire haplotype, but not with the DRB1*1104 or the DQA1*0103-DQB1*0603 alleles when carried separately. These data show a strong linkage disequilibrium of TAP1*0301 and TAP2D alleles with the DRB1*1104-DQA1*0103-DQB1*0603 haplotype of Ashkenazi Jews, extended on the centromeric and telomeric side to the DPB1*0201 and A1-B35 alleles, respectively.     

Haplotypic relationship between SNP and microsatellite markers at the NOS2A locus in two populations

The density of genetic markers required for successful association mapping of complex diseases depends on linkage disequilibrium (LD) between non-functional markers and functional variants. The haplotypic relationship between stable markers and potentially unstable but highly informative markers (e.g. microsatellites) indicates that LD might be maintained over considerable genetic distance in non-African populations, supporting the use of such 'mixed marker haplotypes' in LD-based mapping, and allowing inferences to be drawn about human origins. We investigated sequence variation in the proximal 2.6 kb of the inducible nitric oxide synthase (NOS2A) promoter and the relationship between SNP haplotypes and a pentanucleotide microsatellite (the 'NOS2A(-2.6) microsatellite') in Gambians and UK Caucasians. UK Caucasians exhibited a subset of sequence diversity observed in Gambians, sharing four of 11 SNPs and a similar haplotypic structure. Five SNPs were found in the sequence of interspersed repetitive DNA elements. In both populations, there was dramatic loss of LD between SNP haplotypes and microsatellite alleles across a very short physical distance, suggesting a high intrinsic mutation rate of the NOS2A(-2.6) microsatellite, the SNP haplotypes are relatively ancient, or that this was a region of frequent recombination. Understanding locus- and population-specific LD is essential when designing and interpreting genetic association studies.     

A chloroplast DNA mutational hotspot and gene conversion in a noncoding region near rbcL in the grass family (Poaceae)

The noncoding DNA region of the chloroplast genome, flanked by the genes rbcL and psaI (ORF36), has been sequenced for seven species of the grass family (Poaceae). This region had previously been observed as a hotspot area for length mutations. Sequence comparison reveals that short duplications, probably resulting from slipped-strand mispairing, account for many small length differences between sequences but that major mutational hotspots are localized in three small areas, two of which show potential secondary structure. Mutation in one of these hotspots appears to be a result of more complex recombination events. All seven species contain a pseudogene for rpl23 and evidence is presented that this pseudogene is being maintained by gene conversion with the functional gene. Different transition/transversion biases and AT contents between the pseudogene and the surrounding noncoding sequences are noted. In the subfamily Panicoideae there is a deletion in which almost 1 kb of ancestral sequence, including the 3 end of the rpl23 pseudogene, has been replaced by a non-homologous 60-base sequence of unknown origin. Two other deletions of almost the same region have occurred in the grass family. The deleted noncoding region has mutational and compositional properties similar to the rbcL coding sequence and the rpl23 pseudogene. The three independent deletions, as well as the pattern of mutation in the localized hotspots, indicate that such noncoding DNA may be misleading for studies of phylogenetic inference.