400 likes | 548 Views
Signals of natural selection in the HapMap project data The International HapMap Consortium. Gil McVean Department of Statistics, Oxford University. The International HapMap Project. To facilitate the design and analysis of association studies
E N D
Signals of natural selection in the HapMap project dataThe International HapMap Consortium Gil McVean Department of Statistics, Oxford University
The International HapMap Project • To facilitate the design and analysis of association studies • A genome-wide map of genetic variation across 270 individuals from four populations • CEPH families from Utah • Yoruba from Nigeria • Han Chinese from Beijing • Japanese from the Tokyo region • Phase I collected data on approximately 1.2 million SNPs • Phase II increases SNP density to more than one per kb • All data publicly available at www.hapmap.org
Looking for selection • A genome-wide map of variation can also be used to hunt for regions of the genome where natural selection has acted • Selective sweeps • Balancing selection • Local adaptation • Why? • Interest • Functional polymorphism • The signal of selection we observe tells us about the genetic architecture of traits
Methods for mapping selection • Model-based • Compare genetic variation to ‘neutral’ model • Purely empirical • Consider the ‘most extreme’ genomic regions • ‘Calibrated’ • Compare to examples of (very few) proven selective importance
In what way are selected regions unusual? (in the HapMap data)
HLA and resistance to infectious disease HLA The HLA region shows extremely high levels of polymorphism
17q21 inversion and reproductive success The inversion has multiple (66) SNPs in perfect association (r2 = 1)
LCT and lactase persistence The LCT gene shows an extended haplotype structure in European populations
The Duffy locus and resistance to Plasmodium vivax The FY gene shows extreme population differentiation
Different selective histories leave different footprints in genetic variation
How much of the genome looks as ‘unusual’ as these selected loci?
Sets of perfect proxies as extreme as the 17q21 inversion Inversion
EHH as extreme as LCT Lactase
Differentiation as extreme as the Duffy locus (NB not FY*O) Duffy
For ¾ cases, the selected locus is at the very extreme of the genome-wide distribution
Heterozygosity across the genome Top 1% Top 5% Top 10% Bottom 10% Bottom 5% Bottom 1%
Elevated heterozygosity on 8p Chromosome 6 MHC Chromosome 8 8p23 inversion
Distribution of long runs of perfect proxies ≥ 50 SNPs 20 – 50 SNPs 10-20 SNPs 17q21 Inversion
Distribution of EHH Top 0.1% Top 1% Top 10%
Distribution of differentiation Top 0.1% Top 1% Top 10%
SLC24A5 Lamason et al (Science 2005)
Unusual regions of the genome suggest interesting biology BUT The hypothesis of historical selection is fundamentally untestable
What hypothesis can we test? Signals of selection should tend to occur near regions of known functional importance i.e. genes
Are genes over-represented in regions of high heterozygosity?
Are genes over-represented in regions of high differentiation?
Only differentiation shows a tendency for an increased density of ‘selection’ near genes
Selection on standing variation • Why should we see an excess of one type of signal of adaptive evolution near genes, but not another? • Perhaps the signals are sensitive to assumptions about selection occurs? • EHH methods will be most powerful for identifying selection on a single, novel mutation • Differentiation will pick cases where an already polymorphic mutation, present on multiple haplotype backgrounds, becomes favoured in one geographic region • Perhaps most selection has been on standing variation?
Acknowledgements • The International HapMap Consortium • Oxford Statistics • Peter Donnelly, Simon Myers, Chris Spencer, Raphaelle Chaix • Funding agencies • NIH, TSC, The Wellcome Trust, BBSRC, the Fyssen Foundation
Distribution of Fay and Wu’s H statistic Bottom 0.1% Bottom 1% Bottom 10%
Distribution of Tajima D statistic Top 1% Top 5% Top 10% Bottom 10% Bottom 5% Bottom 1%
Fay and Wu H (negative) Tajima D (negative)