Abstract
We genotyped 2,861 cases of primary biliary cirrhosis (PBC) from the UK PBC Consortium and 8,514 UK population controls across 196,524 variants within 186 known autoimmune risk loci. We identified 3 loci newly associated with PBC (at P < 5 × 10−8), increasing the number of known susceptibility loci to 25. The most associated variant at 19p12 is a low-frequency nonsynonymous SNP in TYK2, further implicating JAK-STAT and cytokine signaling in disease pathogenesis. An additional five loci contained nonsynonymous variants in high linkage disequilibrium (LD; r2 > 0.8) with the most associated variant at the locus. We found multiple independent common, low-frequency and rare variant association signals at five loci. Of the 26 independent non–human leukocyte antigen (HLA) signals tagged on the Immunochip, 15 have SNPs in B-lymphoblastoid open chromatin regions in high LD (r2 > 0.8) with the most associated variant. This study shows how data from dense fine-mapping arrays coupled with functional genomic data can be used to identify candidate causal variants for functional follow-up.
This is a preview of subscription content, access via your institution
Access options
Subscribe to this journal
Receive 12 print issues and online access
206,07 € per year
only 17,17 € per issue
Buy this article
- Purchase on SpringerLink
- Instant access to full article PDF
Prices may be subject to local taxes which are calculated during checkout


Similar content being viewed by others
References
Kaplan, M.M. & Gershwin, M.E. Primary biliary cirrhosis. N. Engl. J. Med. 353, 1261–1273 (2005).
James, O.F. et al. Primary biliary cirrhosis once rare, now common in the United Kingdom? Hepatology 30, 390–394 (1999).
Jones, D.E., Watt, F.E., Metcalf, J.V., Bassendine, M.F. & James, O.F. Familial primary biliary cirrhosis reassessed: a geographically-based population study. J. Hepatol. 30, 402–407 (1999).
Hirschfield, G.M. et al. Primary biliary cirrhosis associated with HLA, IL12A, and IL12RB2 variants. N. Engl. J. Med. 360, 2544–2555 (2009).
Liu, X. et al. Genome-wide meta-analyses identify three loci associated with primary biliary cirrhosis. Nat. Genet. 42, 658–660 (2010).
Mells, G.F. et al. Genome-wide association study identifies 12 new susceptibility loci for primary biliary cirrhosis. Nat. Genet. 43, 329–332 (2011).
Zhernakova, A., van Diemen, C.C. & Wijmenga, C. Detecting shared pathogenesis from the shared genetics of immune-related diseases. Nat. Rev. Genet. 10, 43–55 (2009).
Cortes, A. & Brown, M. Promise and pitfalls of the Immunochip. Arthritis Res. Ther. 13, 101 (2011).
Trynka, G. et al. Dense genotyping identifies and localizes multiple common and rare variant association signals in celiac disease. Nat. Genet. 43, 1193–1201 (2011).
Cordell, H.J. & Clayton, D.G. A unified stepwise regression procedure for evaluating the relative effects of polymorphisms within a gene using case/control or family data: application to HLA in type 1 diabetes. Am. J. Hum. Genet. 70, 124–141 (2002).
Ban, M. et al. Replication analysis identifies TYK2 as a multiple sclerosis susceptibility factor. Eur. J. Hum. Genet. 17, 1309–1313 (2009).
Wallace, C. et al. The imprinted DLK1-MEG3 gene region on chromosome 14q32.2 alters susceptibility to type 1 diabetes. Nat. Genet. 42, 68–71 (2010).
Strange, A. et al. A genome-wide association study identifies new psoriasis susceptibility loci and an interaction between HLA-C and ERAP1. Nat. Genet. 42, 985–990 (2010).
Franke, A. et al. Genome-wide meta-analysis increases to 71 the number of confirmed Crohn's disease susceptibility loci. Nat. Genet. 42, 1118–1125 (2010).
Hunt, K.A. et al. Newly identified genetic risk variants for celiac disease related to the immune response. Nat. Genet. 40, 395–402 (2008).
Stahl, E.A. et al. Genome-wide association study meta-analysis identifies seven new rheumatoid arthritis risk loci. Nat. Genet. 42, 508–514 (2010).
Barrett, J.C. et al. Genome-wide association study and meta-analysis find that over 40 loci affect risk of type 1 diabetes. Nat. Genet. 41, 703–707 (2009).
Wang, D. et al. CD4+CD25+ but not CD4+Foxp3+ T cells as a regulatory subset in primary biliary cirrhosis. Cell. Mol. Immunol. 7, 485–490 (2010).
Rong, G. et al. Imbalance between T helper type 17 and T regulatory cells in patients with primary biliary cirrhosis: the serum cytokine profile and peripheral cell population. Clin. Exp. Immunol. 156, 217–225 (2009).
Couturier, N. et al. Tyrosine kinase 2 variant influences T lymphocyte polarization and multiple sclerosis susceptibility. Brain 134, 693–703 (2011).
Zhernakova, A. et al. Evolutionary and functional analysis of celiac risk loci reveals SH2B3 as a protective factor against bacterial infection. Am. J. Hum. Genet. 86, 970–977 (2010).
Donaldson, P.T. et al. HLA class II alleles, genotypes, haplotypes, and amino acids in primary biliary cirrhosis: a large-scale study. Hepatology 44, 667–674 (2006).
Mullarkey, M.E. et al. Human leukocyte antigen class II alleles in Caucasian women with primary biliary cirrhosis. Tissue Antigens 65, 199–205 (2005).
Wassmuth, R. et al. HLA class II markers and clinical heterogeneity in Swedish patients with primary biliary cirrhosis. Tissue Antigens 59, 381–387 (2002).
Invernizzi, P. et al. Human leukocyte antigen polymorphisms in Italian primary biliary cirrhosis: a multicenter study of 664 patients and 1992 healthy controls. Hepatology 48, 1906–1912 (2008).
Leslie, S., Donnelly, P. & McVean, G. A statistical method for predicting classical HLA alleles from SNP data. Am. J. Hum. Genet. 82, 48–56 (2008).
Dilthey, A.T., Moutsianas, L., Leslie, S. & McVean, G. HLA*IMP—an integrated framework for imputing classical HLA alleles from SNP genotypes. Bioinformatics 27, 968–972 (2011).
Gaffney, D.J. et al. Dissecting the regulatory architecture of gene expression QTLs. Genome Biol. 13, R7 (2012).
Song, L. et al. Open chromatin defined by DNaseI and FAIRE identifies regulatory elements that shape cell-type identity. Genome Res. 21, 1757–1767 (2011).
ENCODE Project Consortium. A user's guide to the encyclopedia of DNA elements (ENCODE). PLoS Biol. 9 e1001046 (2011).
Shah, T.S. et al. optiCall: a robust genotype-calling algorithm for rare, low frequency and common variants. Bioinformatics 28, 1598–1603 (2012).
Teo, Y.Y. et al. A genotype calling algorithm for the Illumina BeadArray platform. Bioinformatics 23, 2741–2746 (2007).
Giannoulatou, E., Yau, C., Colella, S., Ragoussis, J. & Holmes, C.C. GenoSNP: a variational Bayes within-sample SNP genotyping algorithm that does not require a reference population. Bioinformatics 24, 2209–2214 (2008).
Price, A.L. et al. Principal components analysis corrects for stratification in genome-wide association studies. Nat. Genet. 38, 904–909 (2006).
Devlin, B. & Roeder, K. Genomic control for association studies. Biometrics 55, 997–1004 (1999).
Howie, B.N., Donnelly, P. & Marchini, J. A flexible and accurate genotype imputation method for the next generation of genome-wide association studies. PLoS Genet. 5, e1000529 (2009).
Purcell, S. et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 81, 559–575 (2007).
Morris, J.A., Randall, J.C., Maller, J.B. & Barrett, J.C. Evoker: a visualization tool for genotype intensity data. Bioinformatics 26, 1786–1787 (2010).
So, H.C., Gui, A.H., Cherny, S.S. & Sham, P.C. Evaluating the heritability explained by known susceptibility variants: a survey of ten complex diseases. Genet. Epidemiol. 35, 310–317 (2011).
Falconer, D. & Mackay, T. Introduction to Quantitative Genetics 464 (Longman, London, 1996).
Acknowledgements
We are grateful to the PBC Foundation for helping us to establish the PBC Genetics Study, for endorsing it and for encouraging members of the Foundation to contribute samples. We thank all of the research nurses who assisted with participant recruitment in collaborating centers. We thank the staff in the UK National Institute for Health Research Clinical Research Network (NIHR CRN) and Clinical Research Collaboration (CRC) Cymru for providing invaluable support. We are grateful to K. Chittock and his colleagues at Source Bioscience for performing DNA extraction. We thank O. Burren for designing the participant database and for providing information technology support. We are grateful to A. Dilthey for providing support regarding HLA*IMP. We thank J. Stone for coordinating Immunochip design and production at Illumina. We also thank the members of each disease consortium who initiated and sustained the cross-disease Immunochip project and shared control genotypes. Finally, we thank the individuals who contributed the DNA samples used in this study.
The PBC sample collection was funded by the Isaac Newton Trust, the PBC Foundation, The Addenborooke's Charitable Trust and the Wellcome Trust (085925/Z/08/Z). The PBC Genetics Study is a portfolio study of the UK NIHR CRN (portfolio reference 5630). The project is also supported by the Wellcome Trust (WT090355/A/09/Z, WT090355/B/09/Z and 098051). Genotyping of samples at the University of Virginia used resources provided by the Type 1 Diabetes Genetics Consortium, a collaborative clinical study sponsored by the US National Institute of Diabetes and Digestive and Kidney Diseases, the National Institute of Allergy and Infectious Diseases, the National Human Genome Research Institute, the National Institute of Child Health and Human Development and the Juvenile Diabetes Research Foundation International and is supported by U01-DK-062418 from the US National Institutes of Health. We would like to thank the UK Medical Research Council and Wellcome Trust for funding the collection of DNA for the British 1958 Birth Cohort (MRC grant G0000934, Wellcome Trust grant 068545/Z/02). We acknowledge use of DNA from The UK Blood Services collection of Common Controls (UKBS collection), funded by Wellcome Trust grant 076113/C/04/Z, by Wellcome Trust/Juvenile Diabetes Research Foundation grant 061858 and by the National Institute of Health Research of England. The collection was established as part of the Wellcome Trust Case Control Consortium. G.F.M. is a Clinical Research Training Fellow of the Medical Research Council (G0800460). G.F.M. is also supported by a Raymond and Beverly Sackler Studentship.
Author information
Authors and Affiliations
Consortia
Contributions
Study concept and design were performed by D.J.G., G.F.M., L.J., H.J.C., M.A.H., J.M.N., P.T.D., D.E.J., G.J.A., A.J.B., A.B., M.H.D., J.C.B., the WTCCC3 Management Committee (Supplementary Note), R.N.S. and C.A.A. Case ascertainment and phenotyping were performed by S.J.D., D.B.D. and The UK PBC Consortium (Supplementary Note). Control samples were ascertained by The UK Blood Service Controls group (Supplementary Note) and The 1958 Birth Cohort Controls group (Supplementary Note). Genotyping was performed by The WTCCC3 DNA, Genotyping and Informatics group (Supplementary Note). J.Z.L., M.A.A., D.J.G., L.J., The WTCCC3 Data Analysis group (Supplementary Note) and C.A.A. performed statistical analysis. The manuscript was prepared by J.Z.L., M.A.A., D.J.G. and C.A.A. All authors reviewed the final manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Additional information
A full list of members and affiliations is provided in the Supplementary Note.
A full list of members and affiliations is provided in the Supplementary Note.
Supplementary information
Supplementary Text and Figures
Supplementary Tables 1–9, Supplementary Figures 1–12 and Supplementary Note (PDF 6160 kb)
Rights and permissions
About this article
Cite this article
Liu, J., Almarri, M., Gaffney, D. et al. Dense fine-mapping study identifies new susceptibility loci for primary biliary cirrhosis. Nat Genet 44, 1137–1141 (2012). https://doi.org/10.1038/ng.2395
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/ng.2395
This article is cited by
-
Causal relationship between systemic lupus erythematosus and primary liver cirrhosis based on two-sample bidirectional Mendelian randomization and transcriptome overlap analysis
Arthritis Research & Therapy (2024)
-
How genetic risk contributes to autoimmune liver disease
Seminars in Immunopathology (2022)
-
Environmental factors, medical and family history, and comorbidities associated with primary biliary cholangitis in Japan: a multicenter case–control study
Journal of Gastroenterology (2022)
-
Regional heritability mapping identifies several novel loci (STAT4, ULK4, and KCNH5) for primary biliary cholangitis in the Japanese population
European Journal of Human Genetics (2021)
-
rs1944919 on chromosome 11q23.1 and its effector genes COLCA1/COLCA2 confer susceptibility to primary biliary cholangitis
Scientific Reports (2021)