Worldwide patterns of haplotype diversity at 9p21.3, a locus associated with type 2 diabetes and coronary heart disease
Genome Medicine volume 1, Article number: 51 (2009)
A 100 kb region on 9p21.3 harbors two major disease susceptibility loci: one for type 2 diabetes (T2D) and one for coronary heart disease (CHD). The single nucleotide polymorphisms (SNPs) associated with these two diseases in Europeans reside on two adjacent haplotype blocks with independent effects on disease. To help delimit the regions that likely harbor the disease-causing variants in populations of non-European origin, we studied the haplotype diversity and allelic history of the 9p21.3 region using 938 unrelated individuals from 51 populations (Human Genome Diversity Panel). We used SNP data from Illumina's 650Y SNP arrays supplemented with five additional SNPs within the region of interest. Haplotype frequencies were analyzed with the EM algorithm implemented in PLINK. For the T2D locus, the TT risk haplotype of SNPs rs10811661 and rs10757283 was present at similar frequencies in all global populations, while a shared 6-SNP haplotype that carries the protective C allele of rs10811661 was found at a frequency of 2.9% in Africans and 41.3% in East Asians and was associated with low haplotype diversity. For the CHD locus, all populations shared a core risk haplotype spanning >17.5 kb, which shows dramatic increase in frequency between African (11.5%) and Middle Eastern (63.7%) populations. Interestingly, two SNPs (rs2891168 and rs10757278) tagging this CHD risk haplotype are most strongly associated with CHD disease status according to independent clinical fine-mapping studies. The large variation in linkage disequilibrium patterns identified between the populations demonstrates the importance of allelic background data when selecting SNPs for replication in global populations. Intriguingly, the protective allele for T2D and the risk allele for CHD show an increase in frequency in non-Africans compared to Africans, implying different population histories for these two adjacent disease loci.
A 100 kb region on chromosome 9p21.3 has been recently identified as harboring susceptibility variants to both coronary heart disease (CHD)/myocardial infarction [1–5] and to type 2 diabetes (T2D) [6, 7] in study populations of European origin. The associated variants reside on two adjacent haplotype blocks, as defined using HapMap data from 30 trios of Northern and Western European ancestry (CEU) collected by the Centre d'Etude du Polymorphisme Humain (CEPH) . The effects of the variants are independent: the T2D risk variant does not confer increased risk of CHD, and vice versa [9, 10]. The variants associated with CHD also contribute to the risk of other disease phenotypes, such as abdominal aortic and intracranial aneurysms  and ischemic stroke . The extensive linkage disequilibrium (LD) seen in this genomic region in the CEU HapMap population is diminished in the other HapMap populations of African, Chinese and Japanese ancestry, both for pair-wise LD levels as well as size of the haplotype blocks. In addition, the disease-associated variants are located in a genomic region of unknown function. These two issues contribute to the challenge of identifying the possible causative variants, and studying their effects in populations of non-European ancestry.
To help delimit the regions that likely harbor the disease-causing variants in populations of non-European origin, we studied the haplotype diversity and allelic history of the 9p21.3 region using existing genotype data from 938 unrelated individuals from 51 populations, from Sub-Saharan Africa, North Africa, Europe, the Middle East, South/Central Asia, East Asia, Oceania and the Americas (the Human Genome Diversity Panel (HGDP-CEPH) ). Descriptions of genome-wide single nucleotide polymorphism (SNP) variation across these populations have recently been published by two independent groups, using Illumina's HumanHap550 BeadChip in 29 populations  and HumanHap650Y BeadChip in 51 populations . Here we present an analysis of the 650Y SNP data , supplemented with five additional SNPs: rs11790231, rs10965227, rs7045889 and rs10811661 typed using Sequenom iPLEX chemistry (Sequenom, San Diego, CA, USA) and rs1333049 typed using KASPar chemistry (Kbioscience, Hoddesdon, Herts, UK) . The additional SNPs were selected to include disease-associated SNPs and haplotype-tagging SNPs that were not present on the 650Y chip. Genotype quality controls included eight duplicates of a CEPH sample and eight water controls in every 384-well plate. Genotype clusters were manually reviewed and genotyping success rate for each SNP was >98.9%. The genotype data included the variants associated with T2D (rs10811661 and rs10757283), as well as eight variants associated with CHD or found in high LD (r2 > 0.9) with associated variants across a 44 kb region (rs10116277, rs1537370, rs10738607, rs4977574, rs944797, rs2383207, rs1537375 and rs1333049) in the HapMap CEU population. Haplotype frequencies were analyzed with the EM algorithm implemented in PLINK , which estimates the frequencies of probabilistically inferred sets of haplotypes within a population-based sample set. Haplotype structure was visualized with Haploview . The analysis was performed for each geographic region separately, including in the analysis all individuals from the various populations in that geographic region. Due to the small number of unrelated individuals studied here and the uncertainty of phasing, we omit region-specific haplotypes with frequencies <5%.
Our analyses show that the T2D and CHD loci have different allelic histories, which is in agreement with their independent effects on disease. The haplotype structure of the critical region containing the CHD- and T2D-associated SNPs in the HGDP European sample is shown in Figure 1, and for comparison the same region is shown in the HGDP African sample (Additional data file 1). For the T2D locus, the T allele of rs10811661 was found to be associated with disease risk , while a larger meta-analysis identified haplotype TT of rs10811661 and rs10757283 as most strongly associated with disease . The TT risk haplotype is present in similar frequencies in all global populations, while a shared 6-SNP haplotype that carries the protective C allele of rs10811661 is found at a frequency of 2.9% in Africans and 41.3% in East Asians and is associated with low haplotype diversity (Table 1). This frequency difference between populations and lack of haplotype diversity of the protective allele is reminiscent of the TCF7L2 T2D locus, in which the protective allele is found at a frequency of 10 to 31% in Africans but at 95% in East Asians [15, 18]. Such large allele frequency differences and lack of haplotype diversity are indications of the past action of positive natural selection . However, the degree of population differentiation for rs10811661 is not unusual compared to random SNPs in the genome (Fst = 0.126 (P = 0.224) across the 51 HGDP populations) , suggesting a neutrally selected region, while the protective allele of rs7901695 at the TCF7L2 locus was likely driven to high frequency in East Asians (global Fst = 0.213 (P = 0.08) across the 51 HGDP populations) by positive selection .
The risk allele frequencies of four of the CHD-associated SNPs are shown in Figure 2. Although these SNPs show highly similar allele frequencies and are in almost perfect LD in European populations (r2 > 0.9), they show dramatic differences in allele frequencies across other populations, most notably in African populations. In order to decipher which of these risk alleles might be the true causative variant (or in high LD with it) and thus may be suitable for testing in non-European populations, we studied the haplotype diversity across the different geographic regions for eight highly correlated CHD-associated SNPs (r2 > 0.9 in CEU HapMap population; Table 2). All populations appear to share a core risk haplotype as a part of the longer risk haplotype identified in Europeans. This risk haplotype (GGGC, for SNPs rs4977574, rs944797, rs2383207, rs1537375) spans >17.5 kb, and is tagged by the risk allele G of SNP rs4977574. The G allele of rs4977574 is also the best tag SNP for the longer risk haplotype (>44.1 kb) that is most common in all populations. All the other CHD-associated risk alleles were also found on other haplotypes in non-European populations. The risk allele of rs4977574 shows a dramatic change in frequency between African and Middle Eastern populations (Figure 2), and tags the only 8-SNP haplotype of African origin that becomes common in European populations (Table 2). Interestingly, two comprehensive fine-mapping studies of this region in case-control samples have identified SNPs in the same haplotype block as rs4977574 (rs2891168 and rs10757278, shown in Figure 1) as most strongly associated with disease [1, 9]. These three SNPs (rs4977574, rs2891168, and rs10757278) are highly correlated with each other in all four HapMap populations and are the most appropriate for further analyses in non-Europeans. The 44 kb LD region harbors the ANRIL (antisense noncoding RNA in the INK4 locus) gene, which codes for a large antisense non-coding RNA, and was found to be expressed in tissues involved in atherosclerosis . The three CHD risk haplotype tagging SNPs are located in regions of regulatory potential, as defined from alignments of several mammalian sequences [20, 21], and thus may be representing the actual functional domains associated with disease risk.
A handful of replication studies in non-European populations for CHD-related phenotypes have been published to date. Most of these studies make use of populations of East Asian ancestry [22–25] in which patterns of LD are similar to LD patterns in Europeans. Not surprisingly, these studies confirm previously described associations with disease phenotypes discovered in populations of European ancestry. A replication study in a multi-ethnic sample  that included relatively small numbers of cases and controls per ethnic origin confirmed association in Hispanics, but found no association in African Americans, possibly due to the small sample size and the low frequency of the alleles studied. Studies from populations of diverse ancestry are generally lacking. Our results demonstrate the importance of ancestry-specific allelic background when selecting SNPs for replication in global populations, and demonstrate that this approach can complement fine-mapping studies to possibly identify novel putative causative variant/s. Intriguingly, our data imply very different population histories for these two adjacent disease loci, with an increase in the prevalence of the T2D protective allele, most notably in East Asian populations, versus an increase in the prevalence of the CHD risk allele already in Middle Eastern populations. The HGDP SNP data we used here are publicly available and represent a valuable resource for studies of other complex diseases.
Additional data files
The following additional data are available with the online version of this paper. Additional data file 1 is a Powerpoint file showing the pattern of linkage disequilibrium in the African HGDP samples on chromosome 9:22071397-22124172, a region approximately 53 kb long.
Centre d'Etude du Polymorphisme Humain
HapMap data from 30 trios of Northern and Western European ancestry
coronary heart disease
Human Genome Diversity Panel
single nucleotide polymorphism
type 2 diabetes.
Helgadottir A, Thorleifsson G, Manolescu A, Gretarsdottir S, Blondal T, Jonasdottir A, Jonasdottir A, Sigurdsson A, Baker A, Palsson A, Masson G, Gudbjartsson DF, Magnusson KP, Andersen K, Levey AI, Backman VM, Matthiasdottir S, Jonsdottir T, Palsson S, Einarsdottir H, Gunnarsdottir S, Gylfason A, Vaccarino V, Hooper WC, Reilly MP, Granger CB, Austin H, Rader DJ, Shah SH, Quyyumi AA, et al.: A common variant on chromosome 9p21 affects the risk of myocardial infarction. Science. 2007, 316: 1491-1493. 10.1126/science.1142842
Larson MG, Atwood LD, Benjamin EJ, Cupples LA, D'Agostino RB, Fox CS, Govindaraju DR, Guo CY, Heard-Costa NL, Hwang SJ, Murabito JM, Newton-Cheh C, O'Donnell CJ, Seshadri S, Vasan RS, Wang TJ, Wolf PA, Levy D: Framingham Heart Study 100K project: genome-wide associations for cardiovascular disease outcomes. BMC Med Genet. 2007, 8 (Suppl 1): S5- 10.1186/1471-2350-8-S1-S5
McPherson R, Pertsemlidis A, Kavaslar N, Stewart A, Roberts R, Cox DR, Hinds DA, Pennacchio LA, Tybjaerg-Hansen A, Folsom AR, Boerwinkle E, Hobbs HH, Cohen JC: A common allele on chromosome 9 associated with coronary heart disease. Science. 2007, 316: 1488-1491. 10.1126/science.1142447
Samani NJ, Erdmann J, Hall AS, Hengstenberg C, Mangino M, Mayer B, Dixon RJ, Meitinger T, Braund P, Wichmann HE, Barrett JH, Konig IR, Stevens SE, Szymczak S, Tregouet DA, Iles MM, Pahlke F, Pollard H, Lieb W, Cambien F, Fischer M, Ouwehand W, Blankenberg S, Balmforth AJ, Baessler A, Ball SG, Strom TM, Braenne I, Gieger C, Deloukas P, et al.: Genomewide association analysis of coronary artery disease. N Engl J Med. 2007, 357: 443-453. 10.1056/NEJMoa072366
Wellcome Trust Case Control Consortium: Genome-wide association study of 14, 000 cases of seven common diseases and 3, 000 shared controls. Nature. 2007, 447: 661-678. 10.1038/nature05911
Scott LJ, Mohlke KL, Bonnycastle LL, Willer CJ, Li Y, Duren WL, Erdos MR, Stringham HM, Chines PS, Jackson AU, Prokunina-Olsson L, Ding CJ, Swift AJ, Narisu N, Hu T, Pruim R, Xiao R, Li XY, Conneely KN, Riebow NL, Sprau AG, Tong M, White PP, Hetrick KN, Barnhart MW, Bark CW, Goldstein JL, Watkins L, Xiang F, Saramies J, et al.: A genome-wide association study of type 2 diabetes in Finns detects multiple susceptibility variants. Science. 2007, 316: 1341-1345. 10.1126/science.1142382
Zeggini E, Weedon MN, Lindgren CM, Frayling TM, Elliott KS, Lango H, Timpson NJ, Perry JR, Rayner NW, Freathy RM, Barrett JC, Shields B, Morris AP, Ellard S, Groves CJ, Harries LW, Marchini JL, Owen KR, Knight B, Cardon LR, Walker M, Hitman GA, Morris AD, Doney AS, McCarthy MI, Hattersley AT: Replication of genome-wide association signals in UK samples reveals risk loci for type 2 diabetes. Science. 2007, 316: 1336-1341. 10.1126/science.1142364
Frazer KA, Ballinger DG, Cox DR, Hinds DA, Stuve LL, Gibbs RA, Belmont JW, Boudreau A, Hardenbol P, Leal SM, Pasternak S, Wheeler DA, Willis TD, Yu F, Yang H, Zeng C, Gao Y, Hu H, Hu W, Li C, Lin W, Liu S, Pan H, Tang X, Wang J, Wang W, Yu J, Zhang B, Zhang Q, Zhao H, et al.: A second generation human haplotype map of over 3.1 million SNPs. Nature. 2007, 449: 851-861. 10.1038/nature06258
Broadbent HM, Peden JF, Lorkowski S, Goel A, Ongen H, Green F, Clarke R, Collins R, Franzosi MG, Tognoni G, Seedorf U, Rust S, Eriksson P, Hamsten A, Farrall M, Watkins H: Susceptibility to coronary artery disease and diabetes is encoded by distinct, tightly linked SNPs in the ANRIL locus on chromosome 9p. Hum Mol Genet. 2008, 17: 806-814. 10.1093/hmg/ddm352
Helgadottir A, Thorleifsson G, Magnusson KP, Gretarsdottir S, Steinthorsdottir V, Manolescu A, Jones GT, Rinkel GJ, Blankensteijn JD, Ronkainen A, Jaaskelainen JE, Kyo Y, Lenk GM, Sakalihasan N, Kostulas K, Gottsater A, Flex A, Stefansson H, Hansen T, Andersen G, Weinsheimer S, Borch-Johnsen K, Jorgensen T, Shah SH, Quyyumi AA, Granger CB, Reilly MP, Austin H, Levey AI, Vaccarino V, et al.: The same sequence variant on 9p21 associates with myocardial infarction, abdominal aortic aneurysm and intracranial aneurysm. Nat Genet. 2008, 40: 217-224. 10.1038/ng.72
Matarin M, Brown WM, Singleton A, Hardy JA, Meschia JF: Whole genome analyses suggest ischemic stroke and heart disease share an association with polymorphisms on chromosome 9p21. Stroke. 2008, 39: 1586-1589. 10.1161/STROKEAHA.107.502963
Cann HM, de Toma C, Cazes L, Legrand MF, Morel V, Piouffre L, Bodmer J, Bodmer WF, Bonne-Tamir B, Cambon-Thomsen A, Chen Z, Chu J, Carcassi C, Contu L, Du R, Excoffier L, Ferrara GB, Friedlaender JS, Groot H, Gurwitz D, Jenkins T, Herrera RJ, Huang X, Kidd J, Kidd KK, Langaney A, Lin AA, Mehdi SQ, Parham P, Piazza A, et al.: A human genome diversity cell line panel. Science. 2002, 296: 261-262. 10.1126/science.296.5566.261b
Jakobsson M, Scholz SW, Scheet P, Gibbs JR, VanLiere JM, Fung HC, Szpiech ZA, Degnan JH, Wang K, Guerreiro R, Bras JM, Schymick JC, Hernandez DG, Traynor BJ, Simon-Sanchez J, Matarin M, Britton A, Leemput van de J, Rafferty I, Bucan M, Cann HM, Hardy JA, Rosenberg NA, Singleton AB: Genotype, haplotype and copy-number variation in worldwide human populations. Nature. 2008, 451: 998-1003. 10.1038/nature06742
Li JZ, Absher DM, Tang H, Southwick AM, Casto AM, Ramachandran S, Cann HM, Barsh GS, Feldman M, Cavalli-Sforza LL, Myers RM: Worldwide human relationships inferred from genome-wide patterns of variation. Science. 2008, 319: 1100-1104. 10.1126/science.1153717
Myles S, Davison D, Barrett J, Stoneking M, Timpson N: Worldwide population differentiation at disease-associated SNPs. BMC Med Genomics. 2008, 1: 22- 10.1186/1755-8794-1-22
Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, Maller J, Sklar P, de Bakker PI, Daly MJ, Sham PC: PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007, 81: 559-575. 10.1086/519795
Barrett JC, Fry B, Maller J, Daly MJ: Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005, 21: 263-265. 10.1093/bioinformatics/bth457
Helgason A, Palsson S, Thorleifsson G, Grant SF, Emilsson V, Gunnarsdottir S, Adeyemo A, Chen Y, Chen G, Reynisdottir I, Benediktsson R, Hinney A, Hansen T, Andersen G, Borch-Johnsen K, Jorgensen T, Schafer H, Faruque M, Doumatey A, Zhou J, Wilensky RL, Reilly MP, Rader DJ, Bagger Y, Christiansen C, Sigurdsson G, Hebebrand J, Pedersen O, Thorsteinsdottir U, Gulcher JR, et al.: Refining the impact of TCF7L2 gene variants on type 2 diabetes and adaptive evolution. Nat Genet. 2007, 39: 218-225. 10.1038/ng1960
Sabeti PC, Schaffner SF, Fry B, Lohmueller J, Varilly P, Shamovsky O, Palma A, Mikkelsen TS, Altshuler D, Lander ES: Positive natural selection in the human lineage. Science. 2006, 312: 1614-1620. 10.1126/science.1124309
King DC, Taylor J, Elnitski L, Chiaromonte F, Miller W, Hardison RC: Evaluation of regulatory potential and conservation scores for detecting cis-regulatory modules in aligned mammalian genome sequences. Genome Res. 2005, 15: 1051-1060. 10.1101/gr.3642605
Kolbe D, Taylor J, Elnitski L, Eswara P, Li J, Miller W, Hardison R, Chiaromonte F: Regulatory potential scores from genome-wide three-way alignments of human, mouse, and rat. Genome Res. 2004, 14: 700-707. 10.1101/gr.1976004
Chen Z, Qian Q, Ma G, Wang J, Zhang X, Feng Y, Shen C, Yao Y: A common variant on chromosome 9p21 affects the risk of early-onset coronary artery disease. Mol Biol Rep. 2008, 36: 889-893. 10.1007/s11033-008-9259-7
Hinohara K, Nakajima T, Takahashi M, Hohda S, Sasaoka T, Nakahara K, Chida K, Sawabe M, Arimura T, Sato A, Lee BS, Ban JM, Yasunami M, Park JE, Izumi T, Kimura A: Replication of the association between a chromosome 9p21 polymorphism and coronary artery disease in Japanese and Korean populations. J Hum Genet. 2008, 53: 357-359. 10.1007/s10038-008-0248-4
Hiura Y, Fukushima Y, Yuno M, Sawamura H, Kokubo Y, Okamura T, Tomoike H, Goto Y, Nonogi H, Takahashi R, Iwai N: Validation of the association of genetic variants on chromosome 9p21 and 1q41 with myocardial infarction in a Japanese population. Circ J. 2008, 72: 1213-1217. 10.1253/circj.72.1213
Shen GQ, Li L, Rao S, Abdullah KG, Ban JM, Lee BS, Park JE, Wang QK: Four SNPs on chromosome 9p21 in a South Korean population implicate a genetic locus that confers high cross-race risk for development of coronary artery disease. Arterioscler Thromb Vasc Biol. 2008, 28: 360-365. 10.1161/ATVBAHA.107.157248
Assimes TL, Knowles JW, Basu A, Iribarren C, Southwick A, Tang H, Absher D, Li J, Fair JM, Rubin GD, Sidney S, Fortmann SP, Go AS, Hlatky MA, Myers RM, Risch N, Quertermous T: Susceptibility locus for clinical and subclinical coronary artery disease at chromosome 9p21 in the multi-ethnic ADVANCE study. Hum Mol Genet. 2008, 17: 2320-2328. 10.1093/hmg/ddn132
We are grateful for the skilled laboratory work of Anne Vikman. LP and KS have been supported by the Academy of Finland Centre of Excellence in Complex Disease Genetics, the Biocentrum Helsinki Foundation, Helsinki, Finland and the Finnish Foundation for Cardiovascular Research. NT has been supported by MRC Centre grant, number #G0600705.
The authors declare that they have no competing interests.
LP and KS conceived of the study, and were in charge of study design and coordination. KS was in charge of the additional genotyping of markers on the Sequenom system, performed the bulk of the statistical analyses, and drafted the manuscript. SM and NT provided genotype data for SNP rs1333049, and SM helped in creating quality figures. HT and LCS provided the Illumina 650Y HGDP cleaned data and provided the HGDP DNA samples. HT and EJ were involved in the statistical analyses of the data. All authors participated in discussing study design and results interpretation, and read, commented and approved the final manuscript.
Electronic supplementary material
Additional data file 1: R2 values between each SNP pair are shown in shades of grey (black R2 = 1, white R2 = 0) and within each box. The SNPs best tagging the disease-associating haplotypes (rs4977574 and rs10811661) are in bold. The positions of two SNPs that have been identified as most strongly associated with CHD in two separate fine-mapping studies of Europeans, rs2891168 and rs10757278 (see main text), are shown above the genomic sequence line. The position of the ANRIL gene is shown in the upper panel, while the CDKN2B gene is located 72 kb upstream of the first SNP shown, rs10116277. (PPT 177 KB)
About this article
Cite this article
Silander, K., Tang, H., Myles, S. et al. Worldwide patterns of haplotype diversity at 9p21.3, a locus associated with type 2 diabetes and coronary heart disease. Genome Med 1, 51 (2009). https://doi.org/10.1186/gm51