The first human disease associations of the major histocompatibility (MHC) genes were reported over half a century ago, providing important insights into the genetic basis of autoimmunity. Since then, numerous association studies have examined the effect of specific MHC variants, also referred to as the human leukocyte antigen (HLA) alleles, on the risk of autoimmune conditions. These initial studies, however, were frequently flawed by low resolution of HLA typing, small cohort sizes, and spurious associations due to the extended linkage disequilibrium across this region. The modern era of genome-wide association studies (GWAS) provided robust validation of the pivotal role of this region in genetic susceptibility to many autoimmune and inflammatory disorders, including SLE, multiple sclerosis, inflammatory bowel disease, type 1 diabetes, rheumatoid arthritis, and many others. In the kidney field, GWAS provided strong evidence for HLA involvement in susceptibility to common forms of immune-mediated glomerular diseases, such as IgA nephropathy (IgAN)1–4 and membranous nephropathy (MN).5 Suggestive HLA signals have also been reported in recent studies of steroid-sensitive nephrotic syndrome6 and lupus nephritis.7 The described effects of HLA are consistently large, with per allele disease risk increasing from 50% to >400% above baseline, highlighting the critical role of HLA in the pathogenesis of these disorders.
The HLA region resides on chromosome 6p21.3 and is among the most gene-dense portions of DNA, with gene products ranging from antigen-binding molecules and receptors to signaling factors. This region is also among the most polymorphic in humans, with 21 core HLA genes coding for proteins critical for the human immune response to infectious pathogens. The region can be subdivided into classic class I, II, and III and extended class I and II. Class I encompasses HLA-A, -B, and -C that function as presenters of peptides to cytotoxic T cells, and class II consists of HLA-DR and HLA-DQ molecules that present epitopes to CD4-positive T cells. The currently known number of unique HLA alleles surpasses 8000 and 2400 for class I and II, respectively, requiring a precise hierarchical nomenclature to discern a specific allele at a given locus: two-digit resolutions differentiate between allele groups, four-digit accuracy between specific alleles, and six- and eight-digit resolutions discern between exonic and intronic HLA variants, respectively. Owning to the complex structure and a high level of linkage disequilibrium in the HLA locus, genetic dissection of specific disease-causing alleles has been challenging. Most frequently, a complex pattern of associations is observed with multiple independent haplotypes conferring risk or protection, such as in type 1 diabetes or IgAN. Moreover, HLA alleles often have pleiotropic effects including opposed effects on the susceptibility to different conditions. For example, the same allele that increases the risk of IgAN may also increase the risk of type 1 diabetes but protect against SLE.4
Idiopathic MN provides another fascinating example of strong genetic effects at the HLA locus. First described as a distinct form of glomerular disease in 1957, MN and its pathogenesis has been the subject of intensive research efforts. Through a series of experiments and animal models, including Heymann nephritis, it has been established that MN represents a kidney-specific antigen-driven autoimmune disease. However, the identity of the human antigen had long been elusive. In 2009, the M-type phospholipase A2 receptor (PLA2R) had finally been identified as the key antigen in the pathogenesis of MN in humans.8 An additional antigen, thrombospondin type-1 domain-containing 7A (THSD7A), was subsequently identified, explaining a smaller fraction of patients with idiopathic MN who were negative for anti-PLA2R antibodies.9
The first report of genetic contributions of the HLA locus to the risk of MN was published in 1979.10 This report was followed by a series of additional HLA association studies, building the case for the involvement of this locus. The first and only GWAS for MN carried out in a relatively small cohorts of French, Dutch, and British patients resulted in an astonishing discovery: two disease association peaks with large effects dominated the Manhattan plot, one at the HLA locus and another at the PLA2R1 locus.5 The PLA2R1 locus may influence the immunogenicity of PLA2R molecule (the primary autoantigen in MN), whereas the production of anti-PLA2R autoantibodies may be further amplified in individuals with permissive MHC haplotypes. Moreover, a strong genetic interaction between these two loci was described. Similar genetic interactions between variants in antigen-presenting and antigen-encoding genes have been reported in only a handful of other autoimmune diseases, mainly in psoriasis11 and type 1 diabetes.12
Although the above GWAS findings have now been replicated in multiple cohorts,13–17 the loci have not yet been fine-mapped, and the exact mechanisms underlying these associations remain unknown. Notably, the original GWAS in Europeans lacked conditional analyses of the region or imputation of classical alleles, thus the precise number of independently associated haplotypes or the identity of specific causal variants had not been defined.18 Two studies in this issue of the Journal of the American Society of Nephrology provide novel insights into the role of specific HLA alleles in idiopathic MN. The first study by Cui et al.19 performed four-digit resolution typing of HLA-DRB1, DQA1, DQB1, and DPB1 genes followed by case-control association analysis in 261 patients and 599 healthy controls, and pointed to two classical alleles, DRB1*1501 and DRB1*0301, with highly significant independent effects on the risk of idiopathic MN among Han Chinese. The second study by Le et al. performed targeted high-throughput sequencing and four-digit resolution HLA analysis in 99 anti-PLA2R–positive MN cases and 100 healthy controls.20 This study confirmed the association of DRB1*1501 with anti-PLA2R–positive MN, and suggested DRB3*0202 as the second independent risk factor, subsequently replicating both of these associations in an independent cohort of 293 cases 285 controls.
In the first study, the DRB1*1501 and DRB1*0301 alleles had large, independent, and genome-wide significant effects with allelic odds ratios of 4.65 and 3.96, respectively. These alleles were also significantly associated with circulating anti-PLA2R antibodies. Similar to single-nucleotide polymorphisms in the original GWAS, both HLA alleles exhibited statistical interaction with the PLA2R1 variant rs4664308. Notably, the homozygosity for risk alleles at the PLA2R1 locus combined with DRB1*1501 or DRB1*0301 positivity conferred up to 30-fold increased risk of kidney disease. This synergistic effect is suggestive of a physical interaction between DRβ1 molecules and the PLA2R epitope(s) during the process of antigen presentation. The analysis of individual amino acid substitutions within MHC proteins further narrowed down the search to the amino acid positions 13 and 71 within the MHC-DRβ1 chain. Both of these positions participate in the formation of the fourth peptide-binding pocket of MHC-DRβ1 chain. Furthermore, in silico modeling of PLA2R peptides presented by this pocket revealed several candidate epitopes that now require experimental validation. Interestingly, modeling sequence variants of PLA2R1 had little effect on T cell epitope prediction, suggesting that coding variants in PLA2R1 are less likely to alter the immunogenicity of its gene product.
In comparison to the first study, Le et al. performed a more comprehensive analysis by targeted high-throughput sequencing of all HLA genes in the region, including class I and class II genes. Despite smaller sample size of the discovery cohort, additional power was gained by stratifying the analysis on the basis of the PLA2R antibody status. Significant associations of two independent alleles, DRB1*1501 and DRB3*0202, were detected in the discovery phase and replicated in additional case-control cohorts. Notably, the effects were significantly larger compared with the first study, with allelic odds ratios of 24.9 and 17.7, respectively. These larger effects are likely explained by differences in the ascertainment of cases: the first study included all cases of idiopathic MN but the second study used “PLA2R-related MN,” defined strictly by anti-PLA2R seropositivity combined with positive glomerular PLA2R antigen staining. Notably, DRB3 alleles were not typed in the first study, but DRB3*0202 resides on the same haplotype as DRB1*0301, thus the results of both studies are technically in full agreement. In fact, the second study suggests that DRB3*0202 is the most likely culprit allele explaining the reported signal at DRB1*0301, considering its stronger statistical significance and the conditional analysis that removes any residual signal at DRB1*0301 after accounting for both DRB1*1501 and DRB3*0202. Interestingly, secondary genotype–phenotype correlation analyses also revealed that DRB1*1501 was strongly associated with an earlier age at kidney disease onset. Among PLA2R-positive MN cases, patients carrying the DRB1*1501 risk allele presented at a median of 35 years of age compared with the median age of 50 for all other patients.
Although these studies present important novel findings, one needs to be mindful of several limitations, including the relatively small sample size of both studies, study cohorts restricted to Han Chinese, and analysis of HLA alleles performed without accounting for genetic rather than self-reported ancestry. Moreover, the study by Cui et al. appears to have missed an important association of DRB3*0202 by not accounting for the full genomic context of the HLA region. Follow-up studies will clearly require a large-scale comprehensive sequencing effort involving the entire locus, also including noncoding segments that may regulate HLA gene expression. Moreover, because the frequencies and effects of individual HLA alleles may vary greatly between populations, such studies would be ideally performed in cohorts of diverse ancestries, including representative cohorts of Asian, European, and African ancestry.
In summary, although the above studies define two critical HLA haplotypes conveying high risk of MN in Han Chinese, it is not yet clear if these haplotypes are equally important in other populations, or if additional variants were missed because of power or coverage limitations. With respect to future studies, an adequately powered multiethnic GWAS combined with sequencing of the HLA region is still greatly needed for MN, not only to refine the HLA signal but also to potentially define additional loci outside of the HLA and PLA2R1 regions. Stratifying cases by PLA2R antibody status clearly reduces heterogeneity and can potentially enhance future discovery efforts. Moreover, it is essential to explore the precise nature of the interaction between PLA2R1 and HLA alleles across different ancestries, and systematically evaluate its role in the context of potential environmental disease triggers.
Disclosures
None.
K.K. is supported, in part, by National Institutes of Health (NIH) grants R01DK105124 and UM1DK100876 from the National Institute of Diabetes and Digestive and Kidney Diseases, and grant U01HG008680 from the National Human Genome Research Institute.
The content is solely the responsibility of the author and does not represent the official views of the NIH.
References
1. Feehally J, Farrall M, Boland A, Gale DP, Gut I, Heath S, Kumar A, Peden JF, Maxwell PH, Morris DL, Padmanabhan S, Vyse TJ, Zawadzka A, Rees AJ, Lathrop M, Ratcliffe PJ: HLA has strongest association with IgA nephropathy in genome-wide analysis. J Am Soc Nephrol 21: 1791–1797, 2010
2. Gharavi AG, Kiryluk K, Choi M, Li Y, Hou P, Xie J, Sanna-Cherchi S, Men CJ, Julian BA, Wyatt RJ, Novak J, He JC, Wang H, Lv J, Zhu L, Wang W, Wang Z, Yasuno K, Gunel M, Mane S, Umlauf S, Tikhonova I, Beerman I, Savoldi S, Magistroni R, Ghiggeri GM, Bodria M, Lugani F, Ravani P, Ponticelli C, Allegri L, Boscutti G, Frasca G, Amore A, Peruzzi L, Coppo R, Izzi C, Viola BF, Prati E, Salvadori M, Mignani R, Gesualdo L, Bertinetto F, Mesiano P, Amoroso A, Scolari F, Chen N, Zhang H, Lifton RP: Genome-wide association study identifies susceptibility loci for IgA nephropathy. Nat Genet 43: 321–327, 2011
3. Kiryluk K, Li Y, Sanna-Cherchi S, Rohanizadegan M, Suzuki H, Eitner F, Snyder HJ, Choi M, Hou P, Scolari F, Izzi C, Gigante M, Gesualdo L, Savoldi S, Amoroso A, Cusi D, Zamboli P, Julian BA, Novak J, Wyatt RJ, Mucha K, Perola M, Kristiansson K, Viktorin A, Magnusson PK, Thorleifsson G, Thorsteinsdottir U, Stefansson K, Boland A, Metzger M, Thibaudin L, Wanner C, Jager KJ, Goto S, Maixnerova D, Karnib HH, Nagy J, Panzer U, Xie J, Chen N, Tesar V, Narita I, Berthoux F, Floege J, Stengel B, Zhang H, Lifton RP, Gharavi AG: Geographic differences in genetic susceptibility to IgA nephropathy: GWAS replication study and geospatial risk analysis. PLoS Genet 8: e1002765, 2012
4. Kiryluk K, Li Y, Scolari F, Sanna-Cherchi S, Choi M, Verbitsky M, Fasel D, Lata S, Prakash S, Shapiro S, Fischman C, Snyder HJ, Appel G, Izzi C, Viola BF, Dallera N, Del Vecchio L, Barlassina C, Salvi E, Bertinetto FE, Amoroso A, Savoldi S, Rocchietti M, Amore A, Peruzzi L, Coppo R, Salvadori M, Ravani P, Magistroni R, Ghiggeri GM, Caridi G, Bodria M, Lugani F, Allegri L, Delsante M, Maiorana M, Magnano A, Frasca G, Boer E, Boscutti G, Ponticelli C, Mignani R, Marcantoni C, Di Landro D, Santoro D, Pani A, Polci R, Feriozzi S, Chicca S, Galliani M, Gigante M, Gesualdo L, Zamboli P, Battaglia GG, Garozzo M, Maixnerová D, Tesar V, Eitner F, Rauen T, Floege J, Kovacs T, Nagy J, Mucha K, Pączek L, Zaniew M, Mizerska-Wasiak M, Roszkowska-Blaim M, Pawlaczyk K, Gale D, Barratt J, Thibaudin L, Berthoux F, Canaud G, Boland A, Metzger M, Panzer U, Suzuki H, Goto S, Narita I, Caliskan Y, Xie J, Hou P, Chen N, Zhang H, Wyatt RJ, Novak J, Julian BA, Feehally J, Stengel B, Cusi D, Lifton RP, Gharavi AG: Discovery of new risk loci for IgA nephropathy implicates genes involved in immunity against intestinal pathogens. Nat Genet 46: 1187–1196, 2014
5. Stanescu HC, Arcos-Burgos M, Medlar A, Bockenhauer D, Kottgen A, Dragomirescu L, Voinescu C, Patel N, Pearce K, Hubank M, Stephens HA, Laundy V, Padmanabhan S, Zawadzka A, Hofstra JM, Coenen MJ, den Heijer M, Kiemeney LA, Bacq-Daian D, Stengel B, Powis SH, Brenchley P, Feehally J, Rees AJ, Debiec H, Wetzels JF, Ronco P, Mathieson PW, Kleta R: Risk HLA-DQA1 and PLA(2)R1 alleles in idiopathic membranous nephropathy. N Engl J Med 364: 616–626, 2011
6. Gbadegesin RA, Adeyemo A, Webb NJ, Greenbaum LA, Abeyagunawardena A, Thalgahagoda S, Kale A, Gipson D, Srivastava T, Lin JJ, Chand D, Hunley TE, Brophy PD, Bagga A, Sinha A, Rheault MN, Ghali J, Nicholls K, Abraham E, Janjua HS, Omoloja A, Barletta GM, Cai Y, Milford DD, O’Brien C, Awan A, Belostotsky V, Smoyer WE, Homstad A, Hall G, Wu G, Nagaraj S, Wigfall D, Foreman J, Winn MP; Mid-West Pediatric Nephrology Consortium: HLA-DQA1 and PLCG2 Are candidate risk loci for childhood-onset steroid-sensitive nephrotic syndrome. J Am Soc Nephrol 26: 1701–1710, 2015
7. Chung SA, Brown EE, Williams AH, Ramos PS, Berthier CC, Bhangale T, Alarcon-Riquelme ME, Behrens TW, Criswell LA, Graham DC, Demirci FY, Edberg JC, Gaffney PM, Harley JB, Jacob CO, Kamboh MI, Kelly JA, Manzi S, Moser-Sivils KL, Russell LP, Petri M, Tsao BP, Vyse TJ, Zidovetzki R, Kretzler M, Kimberly RP, Freedman BI, Graham RR, Langefeld CD; International Consortium for Systemic Lupus Erythematosus Genetics: Lupus nephritis susceptibility loci in women with systemic lupus erythematosus. J Am Soc Nephrol 25: 2859–2870, 2014
8. Beck LH Jr, Bonegio RG, Lambeau G, Beck DM, Powell DW, Cummins TD, Klein JB, Salant DJ: M-type phospholipase A2 receptor as target antigen in idiopathic membranous nephropathy. N Engl J Med 361: 11–21, 2009
9. Tomas NM, Beck LH Jr, Meyer-Schwesinger C, Seitz-Polski B, Ma H, Zahner G, Dolla G, Hoxha E, Helmchen U, Dabert-Gay AS, Debayle D, Merchant M, Klein J, Salant DJ, Stahl RA, Lambeau G: Thrombospondin type-1 domain-containing 7A in idiopathic membranous nephropathy. N Engl J Med 371: 2277–2287, 2014
10. Klouda PT, Manos J, Acheson EJ, Dyer PA, Goldby FS, Harris R, Lawler W, Mallick NP, Williams G: Strong association between idiopathic membranous nephropathy and HLA-DRW3. Lancet 2: 770–771, 1979
11. Strange A, Capon F, Spencer CC, Knight J, Weale ME, Allen MH, Barton A, Band G, Bellenguez C, Bergboer JG, Blackwell JM, Bramon E, Bumpstead SJ, Casas JP, Cork MJ, Corvin A, Deloukas P, Dilthey A, Duncanson A, Edkins S, Estivill X, Fitzgerald O, Freeman C, Giardina E, Gray E, Hofer A, Hüffmeier U, Hunt SE, Irvine AD, Jankowski J, Kirby B, Langford C, Lascorz J, Leman J, Leslie S, Mallbris L, Markus HS, Mathew CG, McLean WH, McManus R, Mössner R, Moutsianas L, Naluai AT, Nestle FO, Novelli G, Onoufriadis A, Palmer CN, Perricone C, Pirinen M, Plomin R, Potter SC, Pujol RM, Rautanen A, Riveira-Munoz E, Ryan AW, Salmhofer W, Samuelsson L, Sawcer SJ, Schalkwijk J, Smith CH, Ståhle M, Su Z, Tazi-Ahnini R, Traupe H, Viswanathan AC, Warren RB, Weger W, Wolk K, Wood N, Worthington J, Young HS, Zeeuwen PL, Hayday A, Burden AD, Griffiths CE, Kere J, Reis A, McVean G, Evans DM, Brown MA, Barker JN, Peltonen L, Donnelly P, Trembath RC; Genetic Analysis of Psoriasis Consortium and the Wellcome Trust Case Control Consortium 2: A genome-wide association study identifies new psoriasis susceptibility loci and an interaction between HLA-C and ERAP1. Nat Genet 42: 985–990, 2010
12. Smyth DJ, Cooper JD, Howson JM, Walker NM, Plagnol V, Stevens H, Clayton DG, Todd JA: PTPN22 Trp620 explains the association of chromosome 1p13 with type 1 diabetes and shows a statistical interaction with HLA class II genotypes. Diabetes 57: 1730–1737, 2008
13. Lv J, Hou W, Zhou X, Liu G, Zhou F, Zhao N, Hou P, Zhao M, Zhang H: Interaction between PLA2R1 and HLA-DQA1 variants associates with anti-PLA2R antibodies and membranous nephropathy. J Am Soc Nephrol 24: 1323–1329, 2013
14. Bullich G, Ballarín J, Oliver A, Ayasreh N, Silva I, Santín S, Díaz-Encarnación MM, Torra R, Ars E: HLA-DQA1 and PLA2R1 polymorphisms and risk of idiopathic membranous nephropathy. Clin J Am Soc Nephrol 9: 335–343, 2014
15. Saeed M, Beggs ML, Walker PD, Larsen CP: PLA2R-associated membranous glomerulopathy is modulated by common variants in PLA2R1 and HLA-DQA1 genes. Genes Immun 15: 556–561, 2014
16. Sekula P, Li Y, Stanescu HC, Wuttke M, Ekici AB, Bockenhauer D, Walz G, Powis SH, Kielstein JT, Brenchley P, Eckardt KU, Kronenberg F, Kleta R, Köttgen A; GCKD Investigators: Genetic risk variants for membranous nephropathy: Extension of and association with other chronic
kidney disease aetiologies. Nephrol Dial Transplant 32: 325–332, 2016
17. Ramachandran R, Kumar V, Kumar A, Yadav AK, Nada R, Kumar H, Kumar V, Rathi M, Kohli HS, Gupta KL, Sakhuja V, Jha V: PLA2R antibodies, glomerular PLA2R deposits and variations in PLA2R1 and HLA-DQA1 genes in primary membranous nephropathy in South Asians. Nephrol Dial Transplant 31: 1486–1493, 2016
18. Kiryluk K: Risk alleles in idiopathic membranous nephropathy. N Engl J Med 364: 2072–2073, 2011
19. Cui Z, Xie LJ, Chen FJ, Pei ZY, Zhang LJ, Qu Z, Huang J, Gu QH, Zhang YM, Wang X, Wang F, Meng LQ, Liu G, Zhou XJ, Zhu L, Lv JC, Liu F, Zhang H, Liao YH, Lai LH, Ronco P, Zhao MH: MHC class II risk alleles and amino acid residues in idiopathic membranous nephropathy. J Am Soc Nephrol 28: 1651–1664, 2017
20. Le WB, Shi JS, Zhang T, Liu L, Qin HZ, Liang S, Zhang YW, Zheng CX, Jiang S, Qin WS, Zhang HT, Liu ZH: HLA-DRB1*15:01 and HLA-DRB3*02:02 in PLA2R-related membranous nephropathy. J Am Soc Nephrol 28: 1642–1650, 2017