Conservation of hotspots for recombination in low-copy repeats associated with the NF1 microdeletion

Several large-scale studies of human genetic variation have provided insights into processes such as recombination that have shaped human diversity. However, regions such as low-copy repeats (LCRs) have proven difficult to characterize, hindering efforts to understand the processes operating in these regions. We present a detailed study of genetic variation and underlying recombination processes in two copies of an LCR (NF1REPa and NF1REPc) on chromosome 17 involved in the generation of NF1 microdeletions and in a third copy (REP19) on chromosome 19 from which the others originated over 6.7 million years ago. We find evidence for shared hotspots of recombination among the LCRs. REP19 seems to contain hotspots in the same place as the nonallelic recombination hotspots in NF1REPa and NF1REPc. This apparent conservation of patterns of recombination hotspots in moderately diverged paralogous regions contrasts with recent evidence that these patterns are not conserved in less-diverged orthologous regions of chimpanzees.     

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.     

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.     

Human genome sequence variation and the influence of gene history,mutation and recombination

Variation in the human genome sequence is key to understanding susceptibility to disease in modern populations and the history of ancestral populations. Unlocking this information requires knowledge of the patterns and underlying causes of human sequence diversity. By applying a new population-genetic framework to two genome-wide polymorphism surveys, we find that the human genome contains sizeable regions (stretching over tens of thousands of base pairs) that have intrinsically high and low rates of sequence variation. We show that the primary determinant of these patterns is shared genealogical history. Only a fraction of the variation (at most 25%) is due to the local mutation rate. By measuring the average distance over which genealogical histories are typically preserved, these data provide the first genome-wide estimate of the average extent of correlation among variants (linkage disequilibrium). The results are best explained by extreme variability in the recombination rate at a fine scale, and provide the first empirical evidence that such recombination 'hot spots' are a general feature of the human genome and have a principal role in shaping genetic variation in the human population.     

Recombination and linkage disequilibrium in Arabidopsis thaliana

Linkage disequilibrium (LD) is a major aspect of the organization of genetic variation in natural populations. Here we describe the genome-wide pattern of LD in a sample of 19 Arabidopsis thaliana accessions using 341,602 non-singleton SNPs. LD decays within 10 kb on average, considerably faster than previously estimated. Tag SNP selection algorithms and 'hide-the-SNP' simulations suggest that genome-wide association mapping will require only 40%-50% of the observed SNPs, a reduction similar to estimates in a sample of African Americans. An Affymetrix genotyping array containing 250,000 SNPs has been designed based on these results; we demonstrate that it should have more than adequate coverage for genome-wide association mapping. The extent of LD is highly variable, and we find clear evidence of recombination hotspots, which seem to occur preferentially in intergenic regions. LD also reflects the action of selection, and it is more extensive between nonsynonymous polymorphisms than between synonymous polymorphisms.     

Conserved noncoding sequences are selectively constrained and not mutation cold spots

Noncoding genetic variants are likely to influence human biology and disease, but recognizing functional noncoding variants is difficult. Approximately 3% of noncoding sequence is conserved among distantly related mammals, suggesting that these evolutionarily conserved noncoding regions (CNCs) are selectively constrained and contain functional variation. However, CNCs could also merely represent regions with lower local mutation rates. Here we address this issue and show that CNCs are selectively constrained in humans by analyzing HapMap genotype data. Specifically, new (derived) alleles of SNPs within CNCs are rarer than new alleles in nonconserved regions (P = 3 x 10(-18)), indicating that evolutionary pressure has suppressed CNC-derived allele frequencies. Intronic CNCs and CNCs near genes show greater allele frequency shifts, with magnitudes comparable to those for missense variants. Thus, conserved noncoding variants are more likely to be functional. Allele frequency distributions highlight selectively constrained genomic regions that should be intensively surveyed for functionally important variation.     

Evidence for co-evolution of gene order and recombination rate

There is increasing evidence in eukaryotic genomes that gene order is not random, even allowing for tandem duplication. Notably, in numerous genomes, genes of similar expression tend to be clustered. Are there other reasons for clustering of functionally similar genes? If genes are linked to enable genetic, rather than physical clustering, then we also expect that clusters of certain genes might be associated with blocks of reduced recombination rates. Here we show that, in yeast, essential genes are highly clustered and this clustering is independent of clustering of co-expressed genes and of tandem duplications. Adjacent pairs of essential genes are preferentially conserved through evolution. Notably, we also find that clusters of essential genes are in regions of low recombination and that larger clusters have lower recombination rates. These results suggest that selection acts to modify both the fine-scale intragenomic variation in the recombination rate and the distribution of genes and provide evidence for co-evolution of gene order and recombination rate.     

A high-resolution recombination map of the human genome

Determination of recombination rates across the human genome has been constrained by the limited resolution and accuracy of existing genetic maps and the draft genome sequence. We have genotyped 5,136 microsatellite markers for 146 families, with a total of 1,257 meiotic events, to build a high-resolution genetic map meant to: (i) improve the genetic order of polymorphic markers; (ii) improve the precision of estimates of genetic distances; (iii) correct portions of the sequence assembly and SNP map of the human genome; and (iv) build a map of recombination rates. Recombination rates are significantly correlated with both cytogenetic structures (staining intensity of G bands) and sequence (GC content, CpG motifs and poly(A)/poly(T) stretches). Maternal and paternal chromosomes show many differences in locations of recombination maxima. We detected systematic differences in recombination rates between mothers and between gametes from the same mother, suggesting that there is some underlying component determined by both genetic and environmental factors that affects maternal recombination rates.     

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.     

Microdeletion encompassing MAPT at chromosome 17q21.3 is associated with developmental delay and learning disability

Recently, the application of array-based comparative genomic hybridization (array CGH) has improved rates of detection of chromosomal imbalances in individuals with mental retardation and dysmorphic features. Here, we describe three individuals with learning disability and a heterozygous deletion at chromosome 17q21.3, detected in each case by array CGH. FISH analysis demonstrated that the deletions occurred as de novo events in each individual and were between 500 kb and 650 kb in size. A recently described 900-kb inversion that suppresses recombination between ancestral H1 and H2 haplotypes encompasses the deletion. We show that, in each trio, the parent of origin of the deleted chromosome 17 carries at least one H2 chromosome. This region of 17q21.3 shows complex genomic architecture with well-described low-copy repeats (LCRs). The orientation of LCRs flanking the deleted segment in inversion heterozygotes is likely to facilitate the generation of this microdeletion by means of non-allelic homologous recombination.     

Recombination rates in admixed individuals identified by ancestry-based inference

Studies of recombination and how it varies depend crucially on accurate recombination maps. We propose a new approach for constructing high-resolution maps of relative recombination rates based on the observation of ancestry switch points among admixed individuals. We show the utility of this approach using simulations and by applying it to SNP genotype data from a sample of 2,565 African Americans and 299 African Caribbeans and detecting several hundred thousand recombination events. Comparison of the inferred map with high-resolution maps from non-admixed populations provides evidence of fine-scale differentiation in recombination rates between populations. Overall, the admixed map is well predicted by the average proportion of admixture and the recombination rate estimates from the source populations. The exceptions to this are in areas surrounding known large chromosomal structural variants, specifically inversions. These results suggest that outside of structurally variable regions, admixture does not substantially disrupt the factors controlling recombination rates in humans.     

Evidence for two independent prostate cancer risk-associated loci in the HNF1B gene at 17q12

We carried out a fine-mapping study in the HNF1B gene at 17q12 in two study populations and identified a second locus associated with prostate cancer risk, approximately 26 kb centromeric to the first known locus (rs4430796); these loci are separated by a recombination hot spot. We confirmed the association with a SNP in the second locus (rs11649743) in five additional populations, with P = 1.7 x 10(-9) for an allelic test of the seven studies combined. The association at each SNP remained significant after adjustment for the other SNP.     

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.     

Dichotomy of single-nucleotide polymorphism haplotypes in olfactory receptor genes and pseudogenes

Substantial efforts are focused on identifying single-nucleotide polymorphisms (SNPs) throughout the human genome, particularly in coding regions (cSNPs), for both linkage disequilibrium and association studies. Less attention, however, has been directed to the clarification of evolutionary processes that are responsible for the variability in nucleotide diversity among different regions of the genome. We report here the population sequence diversity of genomic segments within a 450-kb cluster of olfactory receptor (OR) genes on human chromosome 17. We found a dichotomy in the pattern of nucleotide diversity between OR pseudogenes and introns on the one hand and the closely interspersed intact genes on the other. We suggest that weak positive selection is responsible for the observed patterns of genetic variation. This is inferred from a lower ratio of polymorphism to divergence in genes compared with pseudogenes or introns, high non-synonymous substitution rates in OR genes, and a small but significant overall reduction in variability in the entire OR gene cluster compared with other genomic regions. The dichotomy among functionally different segments within a short genomic distance requires high recombination rates within this OR cluster. Our work demonstrates the impact of weak positive selection on human nucleotide diversity, and has implications for the evolution of the olfactory repertoire.     

Lifespan and mitochondrial control of neurodegeneration

We examine the allometric (comparative scaling) relationships between rates of neurodegeneration resulting from equivalent mutations in a diverse group of genes from five mammalian species with different maximum lifespan potentials. In both retina and brain, rates of neurodegeneration vary by as much as two orders of magnitude and are strongly correlated with maximum lifespan potential and rates of formation of mitochondrial reactive oxygen and nitrogen species (RONS). Cell death in these disorders is directly or indirectly regulated by the intrinsic mitochondrial cell death pathway. Mitochondria are the main source of RONS production and integrate cellular stress signals to coordinate the intrinsic pathway. We propose that these two functions are intimately related and that steady-state RONS-mediated signaling or damage to the mitochondrial stress-integration machinery is the principal factor setting the probability of cell death in response to a diverse range of cellular stressors. This provides a new and unifying framework for investigating neurodegenerative disorders.     

Evidence of en bloc duplication in vertebrate genomes

It has been 30 years since it was first proposed that the vertebrate genome evolved through several rounds of genome-wide duplications (polyploidizations). Despite rapid advances in genetics, including sequencing of the complete genomes of several divergent species, this hypothesis has not been tested rigorously and is still a matter of debate. If polyploidizations occurred during chordate evolution, there should be a network of paralogous regions in the present-day jawed vertebrate (Gnathostomata) genomes. Here we present an investigation of the major histocompatibility complex (MHC) paralogous regions, which we accomplished by characterizing the corresponding region in amphioxus by identifying nine anchor genes and sequencing both the anchor genes and the regions that flank them (a total of 400 kb). Phylogenetic analysis of 31 genes (including the anchor genes) in these regions shows that duplications occurred after the divergence of cephalochordates and vertebrates but before the Gnathostomata radiation. The distribution of human and amphioxus orthologs in their respective genomes and the relationship between these distributions support the en bloc duplication events. Our analysis represents the first step towards demonstrating that the human ancestral genome has undergone polyploidization. Moreover, reconstruction of the pre-duplicated region indicates that one of the duplicated regions retains the ancestral organization.     

Genome-wide association study of prostate cancer identifies a second risk locus at 8q24

Recently, common variants on human chromosome 8q24 were found to be associated with prostate cancer risk. While conducting a genome-wide association study in the Cancer Genetic Markers of Susceptibility project with 550,000 SNPs in a nested case-control study (1,172 cases and 1,157 controls of European origin), we identified a new association at 8q24 with an independent effect on prostate cancer susceptibility. The most significant signal is 70 kb centromeric to the previously reported SNP, rs1447295, but shows little evidence of linkage disequilibrium with it. A combined analysis with four additional studies (total: 4,296 cases and 4,299 controls) confirms association with prostate cancer for rs6983267 in the centromeric locus (P = 9.42 x 10(-13); heterozygote odds ratio (OR): 1.26, 95% confidence interval (c.i.): 1.13-1.41; homozygote OR: 1.58, 95% c.i.: 1.40-1.78). Each SNP remained significant in a joint analysis after adjusting for the other (rs1447295 P = 1.41 x 10(-11); rs6983267 P = 6.62 x 10(-10)). These observations, combined with compelling evidence for a recombination hotspot between the two markers, indicate the presence of at least two independent loci within 8q24 that contribute to prostate cancer in men of European ancestry. We estimate that the population attributable risk of the new locus, marked by rs6983267, is higher than the locus marked by rs1447295 (21% versus 9%).