Search tips
Search criteria 


Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
Nat Genet. Author manuscript; available in PMC 2012 June 2.
Published in final edited form as:
PMCID: PMC3366038

Meta-analysis identifies six new susceptibility loci for atrial fibrillation

Patrick T Ellinor,1,2,3,4,74 Kathryn L Lunetta,5,6,74 Christine M Albert,4,7,8,74 Nicole L Glazer,9,74 Marylyn D Ritchie,10,11,74 Albert V Smith,12,74 Dan E Arking,13,74 Martina Müller-Nurasyid,14,15,16,74 Bouwe P Krijthe,17,74 Steven A Lubitz,1,7,74 Joshua C Bis,9,74 Mina K Chung,18,19,74 Marcus Dörr,20,74 Kouichi Ozaki,21,74 Jason D Roberts,22 J Gustav Smith,23,24 Arne Pfeufer,25,26 Moritz F Sinner,1,6,14 Kurt Lohman,27 Jingzhong Ding,28 Nicholas L Smith,29,30,31 Jonathan D Smith,18,32 Michiel Rienstra,1,33 Kenneth M Rice,34 David R Van Wagoner,18,19 Jared W Magnani,6,35 Reza Wakili,14 Sebastian Clauss,14 Jerome I Rotter,36 Gerhard Steinbeck,14 Lenore J Launer,37 Robert W Davies,38 Matthew Borkovich,22 Tamara B Harris,37 Honghuang Lin,35 Uwe Völker,39 Henry Völzke,40 David J Milan,1,2 Albert Hofman,17,41 Eric Boerwinkle,42 Lin Y Chen,43 Elsayed Z Soliman,44 Benjamin F Voight,24 Guo Li,9 Aravinda Chakravarti,13 Michiaki Kubo,45 Usha B Tedrow,4,7,8 Lynda M Rose,7 Paul M Ridker,4,7,8 David Conen,46 Tatsuhiko Tsunoda,47 Tetsushi Furukawa,48 Nona Sotoodehnia,9,49 Siyan Xu,5,50 Naoyuki Kamatani,51 Daniel Levy,6 Yusuke Nakamura,52 Babar Parvez,53 Saagar Mahida,1 Karen L Furie,4,54 Jonathan Rosand,3,4,54 Raafia Muhammad,53 Bruce M Psaty,9,30,31,55 Thomas Meitinger,25,26 Siegfried Perz,56 H-Erich Wichmann,57,58,59 Jacqueline C M Witteman,17,60 W H Linda Kao,60,61 Sekar Kathiresan,24,62 Dan M Roden,53,63 Andre G Uitterlinden,16,41,64 Fernando Rivadeneira,17,41,64 Barbara McKnight,34 Marketa Sjögren,65 Anne B Newman,66 Yongmei Liu,67 Michael H Gollob,22 Olle Melander,65 Toshihiro Tanaka,21,75 Bruno H Ch Stricker,17,41,64,68,69,75 Stephan B Felix,20,75 Alvaro Alonso,70,75 Dawood Darbar,53,75 John Barnard,71,75 Daniel I Chasman,4,7,75 Susan R Heckbert,30,31,75 Emelia J Benjamin,6,35,50,72,75 Vilmundur Gudnason,12,75 and Stefan Kääb14,73,75


Atrial fibrillation is a highly prevalent arrhythmia and a major risk factor for stroke, heart failure and death1. We conducted a genome-wide association study (GWAS) in individuals of European ancestry, including 6,707 with and 52,426 without atrial fibrillation. Six new atrial fibrillation susceptibility loci were identified and replicated in an additional sample of individuals of European ancestry, including 5,381 subjects with and 1 0,030 subjects without atrial fibrillation (P < 5 × 10−8). Four of the loci identified in Europeans were further replicated in silico in a GWAS of Japanese individuals, including 843 individuals with and 3,350 individuals without atrial fibrillation. The identified loci implicate candidate genes that encode transcription factors related to cardiopulmonary development, cardiac-expressed ion channels and cell signaling molecules.

Genome-wide association studies in individuals of European descent have identified three genomic regions associated with atrial fibrillation on chromosomes 4q25 (PITX2)25, 16q22 (ZFHX3)3,5 and 1q21 (KCNN3)2. Studies of electrocardiographic traits have also identified a number of loci associated with atrial fibrillation6,7. However, despite these findings, much of the heritability of atrial fibrillation remains unexplained, justifying the search for additional genetic variants underlying atrial fibrillation risk8. Large-scale meta-analysis of GWAS results is a powerful method to identify additional genetic variation underlying traits and conditions. We therefore conducted a meta-analysis of multiple well-phenotyped GWAS samples of European ancestry to identify additional atrial fibrillation susceptibility loci.

Six prospective cohort and 10 prevalent study samples contributed to the discovery analysis, which was adjusted for age and sex (Table 1, Online Methods and Supplementary Note). Atrial fibrillation status was systematically ascertained in each sample (Online Methods and Supplementary Note). After application of quality control SNP exclusion criteria in each study (Supplementary Table 1), meta-analysis was performed, applying genomic control to each study. The genomic control inflation factor for the meta-analysis was 1.042 for the full set of SNPs and 1.040 after omitting all SNPs within 500 kb of the association signals that reached genome-wide significance. The quantile-quantile plot of the expected versus observed P-value distributions for association of the 2,609,549 SNPs analyzed is shown (Supplementary Fig. 1). We identified ten loci that exceeded our preset threshold for genome-wide significance (P < 5 × 10−8) (Fig. 1). The three loci most significantly associated with atrial fibrillation were at previously identified atrial fibrillation susceptibility loci on chromosomes 4q25 in PITX2 (rs6817105; P = 1.8 × 10−74)4, 16q22 in ZFHX3 (rs2106261; P = 3.2 × 10−16)3,5 and 1q21 in KCNN3 (rs6666258; P = 2.0 × 10−14)2 (Table 2).

Figure 1
Manhattan plot of meta-analysis results for genome-wide association with atrial fibrillation. The −log10 (P value) is plotted against the physical position of each SNP on each chromosome. The threshold for genome-wide significance, P < ...
Table 1
Subject characteristics
Table 2
Summary of GWAS meta-analysis results with P < 5 × 10−8

Seven new genomic loci were associated with atrial fibrillation with P < 5 × 10−8 in the discovery stage (Table 2). The most significantly associated SNP in each of the seven new loci was genotyped and tested for association with atrial fibrillation in an additional 3,132 to 5,289 independent individuals with atrial fibrillation and 8,159 to 11,148 referent individuals derived from six studies of individuals of European ancestry (Supplementary Table 2). Six of the loci associated with atrial fibrillation in the discovery stage met our criteria for independent replication. Study-specific replication results are detailed (Supplementary Table 3). The results from meta-analysis of the discovery and replication results are shown (Table 2), as are regional plots (Fig. 2). Recognizing that the genes in closest physical proximity to the associated SNPs are not always the causative genes, we report below the genetic associations in order of statistical significance along with the nearest gene.

Figure 2
Regional plots for seven new atrial fibrillation loci in the discovery sample with P < 1 × 10−8. SNPs are plotted by meta-analysis P value and genomic position (NCBI Build 36). The SNP of interest is labeled. The strength of LD ...

The most significant new association in the discovery stage was on chromosome 1q24 (rs3903239; overall P = 8.4 × 10−14) in PRRX1, which encodes a homeodomain transcription factor highly expressed in the developing heart, particularly in connective tissue. Biological interaction between PRRX1 and a related homeobox transcription factor gene, PRRX2, results in abnormalities of great vessel development in a mouse knockout model9. In a separate PRRX1 knockout, fetal pulmonary vasculature development was impaired10.

A second locus was identified on chromosome 7q31 (rs3807989; overall P = 3.6 × 10−12) in CAV1, which encodes caveolin-1, a cellular membrane protein involved in signal transduction. CAV1 is selectively expressed in the atria11, and its knockout has been associated with dilated cardiomyopathy12. The CAV1 protein colocalizes with and negatively regulates the activity of KCNH2 (ref. 13), a potassium channel involved in cardiac repolarization; the corresponding KCNH2 gene was found to be associated with atrial fibrillation in a candidate gene association study, although not in our present analysis14. The top SNP at the CAV1 locus identified in the current study, rs3807989, was previously identified in a GWAS of the PR and QRS intervals and related to atrial fibrillation6,7. The relationship between other previously reported PR-associated loci and atrial fibrillation are reported (Supplementary Table 4). Of note, significant associations with atrial fibrillation were observed for SNPs related to the PR interval in SOX5, TBX5, SCN5A and SCN10A6,7.

The third locus on chromosome 14q23 (rs1152591; overall P = 5.8 × 10−13) was located in an intron of SYNE2, which encodes numerous nesprin-2 isoforms, some of which are highly expressed in the heart and skeletal muscle. Nesprin-2 localizes throughout the sarcomere and is involved in maintaining nuclear structural integrity by anchoring the nucleus to the cytoskeleton. In a candidate gene approach, mutations in SYNE2 were found to segregate in some families with Emery-Dreifuss muscular dystrophy15, which is characterized by skeletal muscle atrophy, cardiomyopathy and cardiac conduction defects.

The fourth locus on chromosome 9q22 (rs10821415; overall P = 4.2 × 1011) was located in an ORF on chromosome 9. Genes at this locus include FBP1 and FBP2, which are important for gluconeogenesis. Autosomal recessive FBP1 deficiency has been described, but cardiovascular features did not seem to be prominent16. Variants at 9q22 have been implicated in the regulation of height, pulmonary function, angiogenesis and attention deficit hyperactivity disorder, although rs10821415 is not in substantial linkage disequilibrium (LD) with any of these SNPs (r2 < 0.30).

A fifth signal was located at 15q24 (rs7164883; overall P = 2.8 × 10−17) in the first intron of HCN4. The HCN4 protein is the predominant cardiac hyperpolarization-activated cyclic nucleotide–gated channel and is highly expressed in the sinoatrial node. HCN4 activity underlies the funny current (If) that governs cardiac pacemaking, and mutations in HCN4 have been associated with various forms of sinus nodal dysfunction17,18.

A sixth locus on chromosome 10q22 (rs10824026; overall P = 4.0 × 10−9) was located 5 kb upstream of SYNPO2L and 20 kb upstream of MYOZ1. The proteins encoded by SYNPO2L and MYOZ1 are both expressed in skeletal and cardiac muscle, localize to the Z-disc and interact with numerous other proteins. However, the precise role of either gene in cardiovascular physiology is unknown19,20. A mouse knockout of MYOZ1 showed increased calcineurin activity and cardiac hypertrophy in response to pressure overload. However, candidate gene approaches have not supported a prominent role for MYOZ1 mutations in causing familial dilated cardiomyopathy21. Of note, the SYNPO2L locus is located within a previously reported atrial fibrillation susceptibility locus identified in a family with autosomal dominant atrial fibrillation22.

One other locus was identified in the meta-analysis of atrial fibrillation in the WNT8A gene (rs2040862; P = 3.2 × 10−8); however, this association failed to replicate in additional independent cohorts with atrial fibrillation (replication P = 0.36; combined P = 2.5 × 10−7).

There was evidence of significant heterogeneity in the discovery meta-analysis at the previously published atrial fibrillation susceptibility signals at 4q25 in PITX2 and at 16q22 in ZFHX3 (Table 2). Effect heterogeneity at the PITX2 locus has already been observed4,23,24.

We then sought to determine whether the top SNPs or their proxies at each locus were associated with alterations in gene expression in an expression quantitative trait locus (eQTL) database, the Genotype-Tissue Expression eQTL browser. The top SNP at the SYNPO2L locus, rs10824026, is in strong LD with a SNP, rs12570126 (R2 = 0.932), that was found to correlate with the expression of MYOZ1 and SYNPO2L (P = 1.5 × 10−6; data derived from lymphoblastoid cell lines in 270 individuals from the HapMap Consortium)25. Furthermore, the top SNP at the SYNPO2L locus is in LD with a nonsynonymous SNP in SYNPO2L, rs3812629 (R2 = 0.8; encoding the SYNPO2L P707L variant with respect to transcript Q9H987-1), that is predicted to be damaging by both the PolyPhen-2 and the SIFT algorithms. None of the other identified atrial fibrillation risk SNPs were associated with variations in gene expression in the Genotype-Tissue Expression eQTL browser.

We next examined the generalizability of our findings by examining our results in a separate GWAS of individuals of Japanese ancestry, including 843 subjects with and 3,350 subjects without atrial fibrillation in the Japan BioBank study (Supplementary Fig. 2a,b). Within the Japanese GWAS, only the chromosome 4q25 (PITX2) locus exceeded the preset threshold for genome-wide significance (rs2634073; odds ratio (OR) = 1.84, 95% confidence interval (CI) = 1.59–2.13; 3.7 × 10−17; Supplementary Figs. 2a and 3 and Supplementary Table 5). At the previously published locus at 16q22 (ZFHX3)3,5, rs12932445 was associated with atrial fibrillation in participants of Japanese ancestry (P = 6.8 × 10−4). The relationship between atrial fibrillation and variants at the KCNN3-PMVK locus on chromosome 1q21 failed to replicate (P = 0.17); however, a regional plot of this locus revealed a distinct signal at the rs7514452 SNP—approximately 375 kb away from rs6666258, the top SNP in the European ancestry sample—that was modestly associated with atrial fibrillation in the Japanese sample (P = 4.9 × 10−5; R2 = 0.002). At PRRX1 and CAV1, the top SNPs in the European samples were also associated with atrial fibrillation in Japanese individuals. For the loci near PRRX1 and C9orf3, alternate SNPs in the Japanese cohort were more significantly associated with atrial fibrillation than the top SNP at each locus in Europeans (Supplementary Fig. 3 and Supplementary Table 5).

Our study was subject to a number of limitations. To maximize both the power and the generalizability of our study, we included all available individuals with atrial fibrillation; thus, some individuals had comorbidities, such as systolic dysfunction and hypertension. However, none of the identified risk variants for atrial fibrillation were strongly associated with systolic dysfunction in the EchoGen Consortium26, a meta-analysis of echocardiographic data from 5 community-based cohorts consisting of over 12,000 individuals of European descent (P value < 1 × 10−5). Further, when our replication results were adjusted for hypertension status, the identified variants remained significantly associated with atrial fibrillation (Supplementary Table 3). Ultimately, the development of a comprehensive risk score incorporating clinical, biochemical and genetic marker data will be necessary to clarify the incremental benefit of our findings in clinical care. Our eQTL analyses were limited to data available within the Genotype-Tissue Expression eQTL browser; future eQTL analyses in cardiac tissue may be helpful in identifying a relationship between the SNPs associated with atrial fibrillation risk and gene expression. Finally, we acknowledge that the identified variants may not be causal but may represent causal elements in the same or different molecular pathways; future statistical, bioinformatic and biological analyses investigating potential genetic interactions are warranted. Fine mapping and deep resequencing will be necessary to uncover the genetic architecture accompanying the identified common atrial fibrillation susceptibility signals.

In summary, our GWAS meta-analysis for atrial fibrillation has identified six new susceptibility loci in or near plausible candidate genes involved in pacemaking activity, signal transduction and cardiopulmonary development. Our results show that atrial fibrillation has multiple genetic associations and identifies new targets for biological investigation.


Methods and any associated references are available in the online version of the paper.

Supplementary Material

Supplementary material



Acknowledgments are contained in the Supplementary Note.


Note: Supplementary information is available in the online version of the paper.


Study concept and design was determined by P.T.E., K.L.L., C.M.A., B.P.K., M.K.C., J.D., H.V., A.H., Y.N., D.M.R., A.G.U., A.B.N., Y.L., T. Tanaka, B.H.C.S., S.B.F., D.D., S.R.H., E.J.B., V.G. and S.K. Acquisition of data was performed by P.T.E., C.M.A., N.L.G., J.C.B., M.K.C., M.D., J.D.R., A.P., M.F.S., J.D., N.L.S., J.D.S., R.W., J.I.R., L.J.L., T.B.H., U.V., A.H., E.Z.S., M.K., U.B.T., D.C., T.F., R.M., B.M.P., T.M., S.P., H.-E.W., J.C.M.W., A.G.U., F.R., M.S., A.B.N., Y.L., M.H.G., B.H.C.S., A.A., S.R.H., E.J.B., V.G. and S.K. Analysis and interpretation of data were performed by P.T.E., K.L.L., N.L.G., M.D.R., A.V.S., D.E.A., M.M.-N., B.P.K., S.A.L., J.C.B., M.D., K.O., J.D.R., J.G.S., M.F.S., K.L., J.D., N.L.S., K.M.R., J.I.R., R.D.W., U.V., G.L., T. Tsunoda, N.S., N.K., B.M.P., S.K., D.M.R., B.M., Y.L., M.H.G., D.D., S.R.H. and E.J.B. Drafting of the manuscript was carried out by P.T.E., K.L.L., B.P.K., S.A.L. and S.K. Critical revision of the manuscript for important intellectual content was carried out by C.M.A., N.L.G., A.V.S., D.E.A., B.P.K., J.C.B., M.K.C., M.D., J.G.S., A.P., M.F.S., J.D., N.L.S., J.D.S., M.R., K.M.R., D.R.V.W., J.W.M., R.W., J.I.R., L.J.L., T.B.H., U.V., H.V., D.J.M., A.H., L.Y.C., E.Z.S., G.L., U.B.T., P.M.R., D.C., N.S., S.M., K.L.F., B.M.P., T.M., S.P., J.C.M.W., D.M.R., A.G.U., F.R., B.M., M.S., A.B.N., Y.L., M.H.G., O.M., B.H.C.S., S.B.F., A.A., D.D., D.I.C., S.R.H., E.J.B. and V.G. Statistical analysis was performed by K.L.L., N.L.G., M.D.R., A.V.S., D.E.A., M.M.-N., B.P.K., J.C.B., J.G.S., M.F.S., K.L., J.D., K.M.R., R.W.D., U.V., G.L., L.M.R., Y.L. and D.I.C. Funding was obtained by P.T.E., C.M.A., M.K.C., J.G.S., N.L.S., J.D.S., J.I.R., L.J.L., T.B.H., H.V., A.H., E.B., U.B.T., P.M.R., B.M.P., T.M., H.-E.W., J.C.M.W., D.M.R., A.G.U., Y.L., O.M., B.H.C.S., S.B.F., A.A., D.D., S.R.H., E.J.B., V.G. and S.K. Study supervision was performed by P.T.E., K.L.L., C.M.A., M.K.C., N.L.S., J.I.R., A.H., P.M.R., B.M.P., J.C.M.W., A.G.U., Y.L., B.H.C.S., S.B.F., S.R.H., E.J.B., V.G. and S.K. P.T.E., K.L.L., C.M.A., N.L.G., M.D.R., A.V.S., D.E.A., M.M.-N., B.P.K., N.A.L., J.C.B., M.K.C., M.D., K.O., T. Tanaka, B.H.C.S., S.B.F., A.A., D.D., J.B., D.I.C., S.R.H., V.G., E.J.B. and S. Kääb had full access to all of the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis.


The authors declare no competing financial interests.


1. Fuster V, et al. ACC/AHA/ESC 2006 Guidelines for the Management of Patients with Atrial Fibrillation: a report of the American College of Cardiology/American Heart Association Task Force on Practice Guidelines and the European Society of Cardiology Committee for Practice Guidelines (Writing Committee to Revise the 2001 Guidelines for the Management of Patients With Atrial Fibrillation): developed in collaboration with the European Heart Rhythm Association and the Heart Rhythm Society. Circulation. 2006;114:e257–e354. [PubMed]
2. Ellinor PT, et al. Common variants in KCNN3 are associated with lone atrial fibrillation. Nat. Genet. 2010;42:240–244. [PMC free article] [PubMed]
3. Benjamin EJ, et al. Variants in ZFHX3 are associated with atrial fibrillation in individuals of European ancestry. Nat. Genet. 2009;41:879–881. [PMC free article] [PubMed]
4. Gudbjartsson DF, et al. Variants conferring risk of atrial fibrillation on chromosome 4q25. Nature. 2007;448:353–357. [PubMed]
5. Gudbjartsson DF, et al. A sequence variant in ZFHX3 on 16q22 associates with atrial fibrillation and ischemic stroke. Nat. Genet. 2009;41:876–878. [PMC free article] [PubMed]
6. Pfeufer A, et al. Genome-wide association study of PR interval. Nat. Genet. 2010;42:153–159. [PMC free article] [PubMed]
7. Holm H, et al. Several common variants modulate heart rate, PR interval and QRS duration. Nat. Genet. 2010;42:117–122. [PubMed]
8. Lubitz SA, et al. Association between familial atrial fibrillation and risk of new-onset atrial fibrillation. J. Am. Med. Assoc. 2010;304:2263–2269. [PMC free article] [PubMed]
9. Bergwerff M, et al. Loss of function of the Prx1 and Prx2 homeobox genes alters architecture of the great elastic arteries and ductus arteriosus. Virchows Arch. 2000;436:12–19. [PubMed]
10. Ihida-Stansbury K, et al. Paired-related homeobox gene Prx1 is required for pulmonary vascular development. Circ. Res. 2004;94:1507–1514. [PubMed]
11. Volonte D, McTiernan CF, Drab M, Kasper M, Galbiati F. Caveolin-3 form heterooligomeric complexes in atrial cardiac myocytes that are required for doxorubicin-induced apoptosis. Am. J. Physiol. Heart Circ. Physiol. 2008;294:H392–H401. [PubMed]
12. Zhao YY, et al. Defects in caveolin-1 cause dilated cardiomyopathy and pulmonary hypertension in knockout mice. Proc. Natl. Acad. Sci. USA. 2002;99:11375–11380. [PubMed]
13. Lin J, et al. The regulation of the cardiac potassium channel (HERG) by caveolin-1. Biochem. Cell Biol. 2008;86:405–415. [PubMed]
14. Sinner MF, et al. The non-synonymous coding IKr-channel variant KCNH2-K897T is associated with atrial fibrillation: results from a systematic candidate gene-based analysis of KCNH2 (HERG) Eur. Heart J. 2008;29:907–914. [PubMed]
15. Zhang Q, et al. Nesprin-1 and -2 are involved in the pathogenesis of Emery Dreifuss muscular dystrophy and are critical for nuclear envelope integrity. Hum. Mol. Genet. 2007;16:2816–2833. [PubMed]
16. Matsuura T, et al. Two newly identified genomic mutations in a Japanese female patient with fructose-1,6-bisphosphatase (FBPase) deficiency. Mol. Genet. Metab. 2002;76:207–210. [PubMed]
17. Schulze-Bahr E, et al. Pacemaker channel dysfunction in a patient with sinus node disease. J. Clin. Invest. 2003;111:1537–1545. [PMC free article] [PubMed]
18. Milanesi R, Baruscotti M, Gnecchi-Ruscone T, DiFrancesco D. Familial sinus bradycardia associated with a mutation in the cardiac pacemaker channel. N. Engl. J. Med. 2006;354:151–157. [PubMed]
19. Frey N, Olson EN. Calsarcin-3, a novel skeletal muscle–specific member of the calsarcin family, interacts with multiple Z-disc proteins. J. Biol. Chem. 2002;277:13998–14004. [PubMed]
20. Beqqali A, et al. CHAP is a newly identified Z-disc protein essential for heart and skeletal muscle function. J. Cell Sci. 2010;123:1141–1150. [PubMed]
21. Arola AM, et al. Mutations in PDLIM3 and MYOZ1 encoding myocyte Z line proteins are infrequently found in idiopathic dilated cardiomyopathy. Mol. Genet. Metab. 2007;90:435–440. [PubMed]
22. Brugada R, et al. Identification of a genetic locus for familial atrial fibrillation. N. Engl. J. Med. 1997;336:905–911. [PubMed]
23. Lubitz SA, et al. Independent susceptibility markers for atrial fibrillation on chromosome 4q25. Circulation. 2010;122:976–984. [PMC free article] [PubMed]
24. Kääb S, et al. Large scale replication and meta-analysis of variants on chromosome 4q25 associated with atrial fibrillation. Eur. Heart J. 2009;30:813–819. [PMC free article] [PubMed]
25. Stranger BE, et al. Population genomics of human gene expression. Nat. Genet. 2007;39:1217–1224. [PMC free article] [PubMed]
26. Vasan RS, et al. Genetic variants associated with cardiac structure and function: a meta-analysis and replication of genome-wide association data. J. Am. Med. Assoc. 2009;302:168–178. [PMC free article] [PubMed]
27. Johnson AD, et al. SNAP: a web-based tool for identification and annotation of proxy SNPs using HapMap. Bioinformatics. 2008;24:2938–2939. [PMC free article] [PubMed]