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.     

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.     

Evaluating and improving power in whole-genome association studies using fixed marker sets

Emerging technologies make it possible for the first time to genotype hundreds of thousands of SNPs simultaneously, enabling whole-genome association studies. Using empirical genotype data from the International HapMap Project, we evaluate the extent to which the sets of SNPs contained on three whole-genome genotyping arrays capture common SNPs across the genome, and we find that the majority of common SNPs are well captured by these products either directly or through linkage disequilibrium. We explore analytical strategies that use HapMap data to improve power of association studies conducted with these fixed sets of markers and show that limited inclusion of specific haplotype tests in association analysis can increase the fraction of common variants captured by 25-100%. Finally, we introduce a Bayesian approach to association analysis by weighting the likelihood of each statistical test to reflect the number of putative causal alleles to which it is correlated.     

Population genomic analysis of outcrossing and recombination in yeast

The budding yeast Saccharomyces cerevisiae has been used by humans for millennia to make wine, beer and bread. More recently, it became a key model organism for studies of eukaryotic biology and for genomic analysis. However, relatively little is known about the natural lifestyle and population genetics of yeast. One major question is whether genetically diverse yeast strains mate and recombine in the wild. We developed a method to infer the evolutionary history of a species from genome sequences of multiple individuals and applied it to whole-genome sequence data from three strains of Saccharomyces cerevisiae and the sister species Saccharomyces paradoxus. We observed a pattern of sequence variation among yeast strains in which ancestral recombination events lead to a mosaic of segments with shared genealogy. Based on sequence divergence and the inferred median size of shared segments (approximately 2,000 bp), we estimated that although any two strains have undergone approximately 16 million cell divisions since their last common ancestor, only 314 outcrossing events have occurred during this time (roughly one every 50,000 divisions). Local correlations in polymorphism rates indicate that linkage disequilibrium in yeast should extend over kilobases. Our results provide the initial foundation for population studies of association between genotype and phenotype in S. cerevisiae.     

Genome-wide association study of leaf architecture in the maize nested association mapping population

US maize yield has increased eight-fold in the past 80 years, with half of the gain attributed to selection by breeders. During this time, changes in maize leaf angle and size have altered plant architecture, allowing more efficient light capture as planting density has increased. Through a genome-wide association study (GWAS) of the maize nested association mapping panel, we determined the genetic basis of important leaf architecture traits and identified some of the key genes. Overall, we demonstrate that the genetic architecture of the leaf traits is dominated by small effects, with little epistasis, environmental interaction or pleiotropy. In particular, GWAS results show that variations at the liguleless genes have contributed to more upright leaves. These results demonstrate that the use of GWAS with specially designed mapping populations is effective in uncovering the basis of key agronomic traits.     

Genome-wide association study of quantitative resistance to southern leaf blight in the maize nested association mapping population

Nested association mapping (NAM) offers power to resolve complex, quantitative traits to their causal loci. The maize NAM population, consisting of 5,000 recombinant inbred lines (RILs) from 25 families representing the global diversity of maize, was evaluated for resistance to southern leaf blight (SLB) disease. Joint-linkage analysis identified 32 quantitative trait loci (QTLs) with predominantly small, additive effects on SLB resistance. Genome-wide association tests of maize HapMap SNPs were conducted by imputing founder SNP genotypes onto the NAM RILs. SNPs both within and outside of QTL intervals were associated with variation for SLB resistance. Many of these SNPs were within or near sequences homologous to genes previously shown to be involved in plant disease resistance. Limited linkage disequilibrium was observed around some SNPs associated with SLB resistance, indicating that the maize NAM population enables high-resolution mapping of some genome regions.     

Characterization of single-nucleotide polymorphisms in coding regions of human genes.

A major goal in human genetics is to understand the role of common genetic variants in susceptibility to common diseases. This will require characterizing the nature of gene variation in human populations, assembling an extensive catalogue of single-nucleotide polymorphisms (SNPs) in candidate genes and performing association studies for particular diseases. At present, our knowledge of human gene variation remains rudimentary. Here we describe a systematic survey of SNPs in the coding regions of human genes. We identified SNPs in 106 genes relevant to cardiovascular disease, endocrinology and neuropsychiatry by screening an average of 114 independent alleles using 2 independent screening methods. To ensure high accuracy, all reported SNPs were confirmed by DNA sequencing. We identified 560 SNPs, including 392 coding-region SNPs (cSNPs) divided roughly equally between those causing synonymous and non-synonymous changes. We observed different rates of polymorphism among classes of sites within genes (non-coding, degenerate and non-degenerate) as well as between genes. The cSNPs most likely to influence disease, those that alter the amino acid sequence of the encoded protein, are found at a lower rate and with lower allele frequencies than silent substitutions. This likely reflects selection acting against deleterious alleles during human evolution. The lower allele frequency of missense cSNPs has implications for the compilation of a comprehensive catalogue, as well as for the subsequent application to disease association.     

Intensely punctate meiotic recombination in the class II region of the major histocompatibility complex

There is considerable interest in understanding patterns of linkage disequilibrium (LD) in the human genome, to aid investigations of human evolution and facilitate association studies in complex disease. The relative influences of meiotic crossover distribution and population history on LD remain unclear, however. In particular, it is uncertain to what extent crossovers are clustered into 'hot spots, that might influence LD patterns. As a first step to investigating the relationship between LD and recombination, we have analyzed a 216-kb segment of the class II region of the major histocompatibility complex (MHC) already characterized for familial crossovers. High-resolution LD analysis shows the existence of extended domains of strong association interrupted by patchwork areas of LD breakdown. Sperm typing shows that these areas correspond precisely to meiotic crossover hot spots. All six hot spots defined share a remarkably similar symmetrical morphology but vary considerably in intensity, and are not obviously associated with any primary DNA sequence determinants of hot-spot activity. These hot spots occur in clusters and together account for almost all crossovers in this region of the MHC. These data show that, within the MHC at least, crossovers are far from randomly distributed at the molecular level and that recombination hot spots can profoundly affect LD patterns.     

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.     

Mitotic recombination is suppressed by chromosomal divergence in hybrids of distantly related mouse strains

Mitotic recombination occurs with high frequency in humans and mice. It leads to loss of heterozygosity (LOH) at important gene loci and can cause disease. However, the genetic modulators of mitotic recombination are not well understood. As recombination depends on a high level of nucleotide sequence homology, we postulate that the frequency of somatic variants derived from mitotic recombination should be diminished in progeny from crosses between strains of mice in which nucleotide sequences have diverged. Here we report that mitotic recombination is suppressed, to various degrees in different tissues, in hybrids of distantly related mouse strains. Reintroduction of greater chromosomal homology by backcrossing restores mitotic recombination in offspring. Thus, chromosomal divergence inhibits mitotic recombination and, consequently, may act as a modifier of cancer susceptibility by limiting the rate of LOH. The suppression of mitotic recombination in some F1 hybrids in which meiotic recombination persists indicates that these processes are differentially affected by chromosomal divergence.     

Transferability of tag SNPs in genetic association studies in multiple populations

A general question for linkage disequilibrium-based association studies is how power to detect an association is compromised when tag SNPs are chosen from data in one population sample and then deployed in another sample. Specifically, it is important to know how well tags picked from the HapMap DNA samples capture the variation in other samples. To address this, we collected dense data uniformly across the four HapMap population samples and eleven other population samples. We picked tag SNPs using genotype data we collected in the HapMap samples and then evaluated the effective coverage of these tags in comparison to the entire set of common variants observed in the other samples. We simulated case-control association studies in the non-HapMap samples under a disease model of modest risk, and we observed little loss in power. These results demonstrate that the HapMap DNA samples can be used to select tags for genome-wide association studies in many samples around the world.     

Genome-wide genetic association of complex traits in heterogeneous stock mice

Difficulties in fine-mapping quantitative trait loci (QTLs) are a major impediment to progress in the molecular dissection of complex traits in mice. Here we show that genome-wide high-resolution mapping of multiple phenotypes can be achieved using a stock of genetically heterogeneous mice. We developed a conservative and robust bootstrap analysis to map 843 QTLs with an average 95% confidence interval of 2.8 Mb. The QTLs contribute to variation in 97 traits, including models of human disease (asthma, type 2 diabetes mellitus, obesity and anxiety) as well as immunological, biochemical and hematological phenotypes. The genetic architecture of almost all phenotypes was complex, with many loci each contributing a small proportion to the total variance. Our data set, freely available at http://gscan.well.ox.ac.uk, provides an entry point to the functional characterization of genes involved in many complex traits.     

Genome-wide association studies of 14 agronomic traits in rice landraces

Uncovering the genetic basis of agronomic traits in crop landraces that have adapted to various agro-climatic conditions is important to world food security. Here we have identified ～ 3.6 million SNPs by sequencing 517 rice landraces and constructed a high-density haplotype map of the rice genome using a novel data-imputation method. We performed genome-wide association studies (GWAS) for 14 agronomic traits in the population of Oryza sativa indica subspecies. The loci identified through GWAS explained ～ 36% of the phenotypic variance, on average. The peak signals at six loci were tied closely to previously identified genes. This study provides a fundamental resource for rice genetics research and breeding, and demonstrates that an approach integrating second-generation genome sequencing and GWAS can be used as a powerful complementary strategy to classical biparental cross-mapping for dissecting complex traits in rice.     

A genome-wide association study of metabolic traits in human urine

We present a genome-wide association study of metabolic traits in human urine, designed to investigate the detoxification capacity of the human body. Using NMR spectroscopy, we tested for associations between 59 metabolites in urine from 862 male participants in the population-based SHIP study. We replicated the results using 1,039 additional samples of the same study, including a 5-year follow-up, and 992 samples from the independent KORA study. We report five loci with joint P values of association from 3.2 × 10(-19) to 2.1 × 10(-182). Variants at three of these loci have previously been linked with important clinical outcomes: SLC7A9 is a risk locus for chronic kidney disease, NAT2 for coronary artery disease and genotype-dependent response to drug toxicity, and SLC6A20 for iminoglycinuria. Moreover, we identify rs37369 in AGXT2 as the genetic basis of hyper-β-aminoisobutyric aciduria.     

Efficient mapping of mendelian traits in dogs through genome-wide association

With several hundred genetic diseases and an advantageous genome structure, dogs are ideal for mapping genes that cause disease. Here we report the development of a genotyping array with approximately 27,000 SNPs and show that genome-wide association mapping of mendelian traits in dog breeds can be achieved with only approximately 20 dogs. Specifically, we map two traits with mendelian inheritance: the major white spotting (S) locus and the hair ridge in Rhodesian ridgebacks. For both traits, we map the loci to discrete regions of <1 Mb. Fine-mapping of the S locus in two breeds refines the localization to a region of approximately 100 kb contained within the pigmentation-related gene MITF. Complete sequencing of the white and solid haplotypes identifies candidate regulatory mutations in the melanocyte-specific promoter of MITF. Our results show that genome-wide association mapping within dog breeds, followed by fine-mapping across multiple breeds, will be highly efficient and generally applicable to trait mapping, providing insights into canine and human health.