A large-scale replication study identifies TNIP1, PRDM1, JAZF1, UHRF1BP1 and IL10 as risk loci for systemic lupus erythematosus.
Journal: 2009/November - Nature Genetics
ISSN: 1546-1718
Abstract:
Genome-wide association studies have recently identified at least 15 susceptibility loci for systemic lupus erythematosus (SLE). To confirm additional risk loci, we selected SNPs from 2,466 regions that showed nominal evidence of association to SLE (P < 0.05) in a genome-wide study and genotyped them in an independent sample of 1,963 cases and 4,329 controls. This replication effort identified five new SLE susceptibility loci (P < 5 x 10(-8)): TNIP1 (odds ratio (OR) = 1.27), PRDM1 (OR = 1.20), JAZF1 (OR = 1.20), UHRF1BP1 (OR = 1.17) and IL10 (OR = 1.19). We identified 21 additional candidate loci with P< or = 1 x 10(-5). A candidate screen of alleles previously associated with other autoimmune diseases suggested five loci (P < 1 x 10(-3)) that may contribute to SLE: IFIH1, CFB, CLEC16A, IL12B and SH2B3. These results expand the number of confirmed and candidate SLE susceptibility loci and implicate several key immunologic pathways in SLE pathogenesis.
Relations:
Content
Citations
(336)
References
(37)
Diseases
(1)
Conditions
(1)
Chemicals
(5)
Genes
(11)
Organisms
(1)
Processes
(2)
Similar articles
Articles by the same authors
Discussion board
Nat Genet 41(11): 1228-1233

A large-scale replication study identifies <em>TNIP1, PRDM1, JAZF1, UHRF1BP1</em> and <em>IL10</em> as risk loci for systemic lupus erythematosus

+27 authors

METHODS

Methods and any associated references are available in the online version of the paper at http://www.nature.com/naturegenetics/.

Supplementary Material

S1

S1

Click here to view.(401K, pdf)

Acknowledgments

We thank the many affected individuals and physicians who contributed DNA samples and clinical data for this study; M.I. Kamboh and P. Davies for the use of Alzheimer’s disease samples as controls in our study; B. Neale for assistance in the percent of genetic variance explained calculation; and S. Sanna and C. Willer for assistance in generating regional association plots. Genotyping of the Swedish samples by the 12K chips was performed using equipment of the SNP technology platform in Uppsala. We thank C. Enström and A.-C. Wiman for assistance with genotyping. Financial support was obtained from the Swedish Research Council for Medicine, the Knut and Alice Wallenberg Foundation the Swedish Rheumatism Association, the King Gustaf V 80th Birthday Foundation, COMBINE, and a Target Identification in Lupus (TIL) grant from the Alliance for Lupus Research, US. This work was supported in part by R01 AR44804, K24 AR02175, the Mary Kirkland Center for Lupus Research, RO1 AR43727 and Institute for Clinical and Translational Research UL1RR025005. These studies were performed in part in the General Clinical Research Center, Moffitt Hospital, University of California, San Francisco, with funds provided by the National Center for Research Resources, 5 M01 RR-00079, US Public Health Service.

Immunology Biomarkers Group, Genentech, South San Francisco, California, USA
Molecular Medicine, Department of Medical Sciences, Uppsala University, Uppsala, Sweden
Rosalind Russell Medical Research Center for Arthritis, Department of Medicine, University of California, San Francisco, California, USA
Rowe Program in Genetics, University of California at Davis, Davis, California, USA
Section of Rheumatology, Department of Medical Sciences, Uppsala University, Uppsala, Sweden
Rheumatology Unit, Department of Medicine, Karolinska Institutet, Karolinska University Hospital, Stockholm, Sweden
Section of Rheumatology, Department of Clinical Sciences, Lund University Hospital, Lund, Sweden
Department of Rheumatology, Umeå University Hospital, Umeå, Sweden
Center for Immunology, University of Minnesota Medical School, Minneapolis, Minnesota, USA
University of Alabama at Birmingham, Birmingham, Alabama, USA
Northwestern University Feinberg School of Medicine, Chicago, Illinois, USA
University of Texas–Houston Health Science Center, Houston, Texas, USA
University of Puerto Rico Medical Science Campus, San Juan, Puerto Rico
University of Pittsburgh Medical Center, Pittsburgh, Pennsylvania, USA
Department of Medicine, Johns Hopkins University School of Medicine, Baltimore, Maryland, USA
The Feinstein Institute for Medical Research, North Shore–Long Island Jewish Health System, Manhasset, New York, USA
Correspondence should be addressed to R.R.G. (moc.eneg@trebor.maharg)

Abstract

Genome-wide association studies have recently identified at least 15 susceptibility loci for systemic lupus erythematosus (SLE). To confirm additional risk loci, we selected SNPs from 2,466 regions that showed nominal evidence of association to SLE (P < 0.05) in a genome-wide study and genotyped them in an independent sample of 1,963 cases and 4,329 controls. This replication effort identified five new SLE susceptibility loci (P < 5 × 10): TNIP1 (odds ratio (OR) = 1.27), PRDM1 (OR = 1.20), JAZF1 (OR = 1.20), UHRF1BP1 (OR = 1.17) and IL10 (OR = 1.19). We identified 21 additional candidate loci with P ≤ 1 × 10. A candidate screen of alleles previously associated with other autoimmune diseases suggested five loci (P < 1 × 10) that may contribute to SLE: IFIH1, CFB, CLEC16A, IL12B and SH2B3. These results expand the number of confirmed and candidate SLE susceptibility loci and implicate several key immunologic pathways in SLE pathogenesis.

Abstract

Systemic lupus erythematosus (SLE) is a chronic inflammatory autoimmune disease characterized by the presence of antibodies to nuclear self-antigens. Many of the lupus autoantibodies recognize nucleic acids and nucleic acid binding proteins, which in turn activate Toll-like receptors, leading to the production of type I interferon1. Despite considerable clinical heterogeneity, SLE ranks among the most heritable of common autoimmune diseases, with a sibling risk ratio of ~30 (ref. 2). Recent genome-wide association (GWA) and candidate gene studies have identified at least 15 common SLE risk alleles that achieve genome-wide significance (P < 5 × 10). These include genes encoding proteins important for adaptive immunity and the production of autoantibodies (HLA class II alleles, BLK, PTPN22 and BANK1) and proteins with roles in innate immunity and interferon signaling (ITGAM, TNFAIP3, STAT4 and IRF5)310. To identify additional risk loci, we performed a targeted replication study of SNPs from 2,466 loci that showed a nominal P value of <0.05 in a recent GWA7 scan of 1,310 individuals with lupus (cases) and 7,859 controls. We also genotyped SNPs from 23 previously reported SLE risk loci, 42 SNPs implicated in other autoimmune diseases and over 7,000 ancestry-informative markers (Fig. 1).

An external file that holds a picture, illustration, etc.
Object name is nihms226372f1.jpg

Overview of experimental design. Variants were selected from loci with P < 0.05 in a genome-wide scan of 1,310 cases and 7,859 controls, previously reported SLE risk loci, confirmed loci from other autoimmune diseases and over 7,000 ancestry-informative markers, and these variants were incorporated into an Illumina custom SNP array. The array was genotyped in independent cases and controls from the United States and Sweden. 823 of the Swedish controls were genotyped using the Illumina 310K SNP array. Variants were analyzed as described in Online Methods.

We designed a custom SNP array (Illumina Infinium II) consisting of over 12,000 variants and genotyped two independent SLE case and control populations from the United States (1,129 SLE cases and 2,991 controls) and Sweden (834 SLE cases and 1,338 controls). Included among the US controls were 2,215 Alzheimer’s disease case-control samples, which were judged to be acceptable as controls because the genetic factors underlying SLE and Alzheimer’s disease are expected to be independent. We next applied data quality filters to remove poorly performing samples and SNPs, population outliers, duplicate samples and related individuals (see Online Methods). Following these quality control measures, we examined a final set of 10,848 SNPs (Fig. 1). Association statistics for 3,735 SNPs were calculated and corrected for population stratification using 7,113 ancestry-informative markers (see Online Methods).

We first examined 25 SNPs (from 23 loci) that were previously reported to be associated with SLE (Table 1 and Supplementary Table 1). We found further evidence of association for 21 of the variants (P < 0.05), including 9 loci that reached genome-wide significance (P < 5 × 10) in the current combined dataset. Among the loci with genome-wide significant results were HLA-DRB1 (HLA*DR3 or DRB1*0301), IRF5, TNFAIP3, BLK, STAT4, ITGAM, PTPN22, PHRF1 (also called KIAA1542) and TNFSF4 (also called OX40L). The analysis also provided additional evidence of association for variants at nine loci for which a single previous study reported genome-wide levels of significance: HLA-DRB1 (HLA*DR2 or DRB1*1501), TNFAIP3 (rs6920220), BANK1, ATG5, PTTG1, PXK, FCGR2A, UBE2L3 and IRAK1-MECP2.

Table 1

Replication results of previously reported SLE risk loci

SNPChrCritical regionP
Gene of interestRisk alleleRisk allele frequencyOR (95% CI)
GWASUSSwedenCombined
Variants with P < 5 × 10 in the current dataset
rs3135394a6p21.3232.027–32.8747.8 × 10−221.8 × 10−268.3 × 10−212.0 × 10−60HLA-DRB1bG0.101.98 (1.84–2.14)
rs7574865a2q32.2191.609–191.6813.0 × 10−196.4 × 10−162.7 × 10−121.4 × 10−41STAT4T0.231.57 (1.49–1.69)
rs2070197a7q32.1128.276–128.476n.a.1.4 × 10−164.1 × 10−95.8 × 10−24IRF5C0.111.88 (1.78–1.95)
rs11860650a16p11.231.195–31.2775.3 × 10−111.8 × 10−59.2 × 10−81.9 × 10−20ITGAMT0.131.43 (1.32–1.54)
rs27363408p23.111.331–11.4885.5 × 10−84.6 × 10−90.00287.9 × 10−17BLKT0.251.35 (1.27–1.43)
rs5029937a6q23.3138.174–138.2841.0 × 10−42.4 × 10−73.1 × 10−55.3 × 10−13TNFAIP3T0.031.71 (1.51–1.95)
rs24766011p13.2113.963–114.2513.3 × 10−54.5 × 10−51.5 × 10−53.4 × 10−12PTPN22A0.101.35 (1.24–1.47)
rs496312811p15.50.485–0.6640.00211.5 × 10−58.7 × 10−44.9 × 10−9PHRF1C0.671.20 (1.13–1.27)
rs22059601q25.1171.454–171.5239.5 × 10−60.0306.7 × 10−46.3 × 10−9TNFSF4T0.231.22 (1.15–1.30)
Variants with a previous report of P < 5 × 10−8
rs9271366a6p21.3232.446–32.6950.00797.4 × 10−48.3 × 10−51.4 × 10−7HLA-DRB1cG0.161.26 (1.18–1.36)
rs6920220a6q23.3138.000–138.0489.9 × 10−45.2 × 10−40.0494.0 × 10−7TNFAIP3A0.211.17 (1.10–1.25)
rs2269368Xq28152.743–152.9432.5 × 10−5n.a.0.00497.5 × 10−7IRAK1-MECP2T0.141.11 (1.01–1.22)
rs24310995q33.3159.813–159.8211.5 × 10−50.160.0471.6 × 10−6PTTG1G0.521.15 (1.09–1.22)
rs575421722q11.220.240–20.3150.00608.4 × 10−40.0182.3 × 10−6UBE2L3T0.191.20 (1.13–1.27)
rs2245214a6q21106.749–106.8760.0324.3 × 10−60.351.2 × 10−5ATG5G0.371.15 (1.09–1.21)
rs105164874q24102.930–103.1340.0970.0910.00158.3 × 10−4BANK1G0.701.11 (1.04–1.18)
rs2176082a3p14.358.214–58.4430.0100.0120.00311.2 × 10−5PXKA0.281.17 (1.10–1.25)
rs18012741q23.3159.724–159.7464.1 × 10−4n.a.n.a.4.1 × 10−4FCGR2AG0.501.16 (1.09–1.20)
Variants with a previous report of P > 5 × 10−8
rs280519a19p13.210.387–10.4307.1 × 10−4n.a.0.0367.4 × 10−5TYK2A0.481.13 (1.06–1.21)
rs101560917p21.38.134–8.1540.0950.00318.7 × 10−46.5 × 10−4ICA1T0.101.16 (1.06–1.27)
rs20220131q25.3181.538–181.6700.262.05 × 10−52.8 × 10−40.0015NMNAT2T0.601.09 (1.03–1.16)
rs78298168q12.156.985–57.0250.490.760.190.17LYNA0.791.05 (0.96–1.17)
rs207172522q13.241.908–41.9700.630.340.290.30SCUBE1G0.861.09 (0.98–1.20)
rs5744168a1q41n.a.n.a.1.000.400.67TLR5G0.941.02 (0.94–1.12)
rs5097491q23.3158.993–159.0670.640.940.930.76LY9G0.961.01 (0.91–1.12)

Critical region here is defined as the minimal region containing variants with r > 0.4 in the HapMap CEU population and is reported in HG18 coordinates (Mb). P values calculated from indicated case-control population (GWAS: 1,310 cases and 7,859 controls; US: 1,129 cases and 2,991 controls; Sweden: 834 cases and 1,338 controls; combined: 3,273 cases and 12,188 controls) and combined P values were calculated as described in Online Methods. Risk allele is reported relative to + reference strand. Risk allele frequency is the frequency in control chromosomes. OR is the combined odds ratio as described in Online Methods.

Indicates markers that were imputed, as described in Online Methods, in the GWAS samples and directly genotyped in the replication samples.
rs3135394 has an r = 0.87 to the HLA*DR3 (DRB1*0301) allele.
rs9271366 has an r = 0.97 to the HLA*DR2 (DRB1*1501) allele. See Supplementary Table 1 for expanded summary statistics. n.a, not available due to failure to pass quality control measures (TYK2, FCGR2A and IRAK1-MECP2), or the specific variant was not present in the genome-wide array (TLR5 and IRF5); however, rs2070197 (IRF5 region) is in strong LD with rs10488631, which had a P = 2 × 10 in the genome scan.

An earlier candidate gene study9 identified MECP2 as a potential risk locus for SLE; however, in the current dataset, SNPs near IRAK1, a gene critical for Toll-like receptor 7 and 9 signaling and located within the identified region of linkage disequilibrium (LD) surrounding MECP2, showed the strongest evidence of association. Similar findings were recently reported11, and further work will be required to determine the causal allele in the IRAK1-MECP2 locus. We found additional evidence of association for three loci (TYK2, ICA1 and NMNAT2) that had previously shown significant but not genome wide–level evidence for association610. For four previously implicated variants (LYN, SCUBE1, TLR5 and LY9), no evidence of association was observed in the combined dataset.

To identify previously unknown SLE risk loci, we examined 3,188 SNPs from 2,446 distinct loci that showed evidence of association to SLE in our genome-wide dataset7, which comprised 502,033 SNPs genotyped in 1,310 SLE cases and an expanded set of 7,859 controls. Using this dataset, we imputed over 2.1 million variants using Phase II HapMap CEU samples as a reference (see Online Methods) and generated a rank-ordered list of association statistics. Variants with P < 0.05 were selected for possible inclusion on the custom replication array. For efficient genotyping, we identified groups of correlated variants (r > 0.2) and then carried out selection of at least two SNPs from each group where all SNPs had P < 0.001. For the remaining groups, the SNP with the lowest P value in the group was included. In the replication samples, we calculated the association statistics (see Online Methods) and observed a significant enrichment of the replication results relative to the expected null distribution (Fig. 2). Excluding previously reported SLE risk alleles, there were 134 loci with P < 0.05 (64 expected; P = 2 × 10) and 12 loci with P < 0.001 (1 expected; P = 1 × 10), suggesting the presence of true positive associations.

An external file that holds a picture, illustration, etc.
Object name is nihms226372f2.jpg

Newly discovered genome-wide significant associations in SLE. (ae) Association results from the GWA scan are plotted on the y axis versus genomic position on the indicated chromosome on the x axis within a 500-kb region surrounding the loci defined by (a) TNIP1, (b) PRDM1, (c) JAZF1, (d) UHRF1BP1 and (e) IL10. The meta-analysis P value for the most strongly associated marker is indicated by a red square. P values from the genome scan are shaded to indicate LD to the genome-wide associated variant: red, r > 0.8; yellow, r > 0.5; green, r > 0.2; gray, r < 0.2. Along the bottom are the recombination rates from the CEU HapMap (light blue line) and the known human genes (blue). A previously reported and independent SLE risk locus at the nearby ATG5 gene is indicated (b; rs2245214). (f) Histogram of P values of 1,256 independent SNPs (r < 0.1 to any other SNP in the array) in the 1,963 case and 4,329 control replication samples. Under a null distribution, the expected density of results is indicated by the dashed line. A significant enrichment of results with P < 0.05 was observed.

The replication study identified five new SLE risk loci with a combined P value that exceeded the genome-wide threshold for significance (P < 5 × 10): TNIP1, PRDM1, JAZF1, UHRF1BP1 and IL10 (Table 2 and Supplementary Table 2). These loci are discussed in more detail below.

Table 2

Newly discovered SLE risk loci in the combined dataset

SNPChr.Critical regionP
Gene of interestRisk alleleRisk allele frequencyOR (95% CI)
GWASUSSwedenCombined
Genome-wide significant loci
rs7708392a5150.419–150.4414.5 × 10−77.7 × 10−41.2 × 10−53.8 × 10−13TNIP1C0.241.27 (1.10–1.35)
rs65684316106.675–106.7056.1 × 10−60.00160.00507.1 × 10−10PRDM1A0.381.20 (1.14–1.27)
rs849142a728.108–28.2234.5 × 10−70.105.4 × 10−41.5 × 10−9JAZF1T0.491.19 (1.13–1.26)
rs11755393a634.658–35.0900.00143.7 × 10−45.1 × 10−42.2 × 10−8UHRF1BP1G0.351.17 (1.10–1.24)
rs30245051205.007–205.0162.6 × 10−60.0621.8 × 10−44.0 × 10−8IL10A0.161.19 (1.11–1.28)
Loci with combined P value < 1 × 10−6
rs10911363a1181.672–181.8162.0 × 10−41.5 × 10−50.529.5 × 10−8NCF2T0.271.19 (1.12–1.26)
rs12444486a1684.548–84.5763.5 × 10−50.0210.0261.9 × 10−7IRF8T0.501.16 (1.10–1.23)
rs11013210a1023.181–23.3371.6 × 10−50.0130.122.0 × 10−7ARMC3T0.211.18 (1.11–1.26)
rs1874791a167.563–67.6873.1 × 10−50.0120.113.4 × 10−7IL12RB2A0.181.18 (1.10–1.26)
rs97829551233.893–234.1076.4 × 10−60.0570.124.6 × 10−7LYSTC0.741.18 (1.11–1.26)
rs7683537a4185.805–185.9141.6 × 10−40.110.00137.6 × 10−7MLF1IPT0.821.23 (1.14–1.33)
rs42807312117.706–117.3151.7 × 10−50.220.00797.7 × 10−7TAOK3T0.691.18 (1.11–1.26)
rs497273a12119.610–119.8915.0 × 10−50.0680.0218.2 × 10−7SPPL3G0.651.14 (1.08–1.21)
Loci with combined P value < 1 × 10−5
rs1861525725.097–25.1838.5 × 10−50.160.00271.9 × 10−6CYCSG0.051.27 (1.12–1.45)
rs921916750.193–50.2054.8 × 10−40.0270.0142.0 × 10−6IKZF1C0.181.15 (1.07–1.23)
rs73336711373.177–73.1982.2 × 10−40.140.00272.2 × 10−6KLF12G0.081.22 (1.11–1.34)
rs12992463222.312–22.4642.1 × 10−50.230.0232.6 × 10−6A0.501.12 (1.06–1.19)
rs126209992237.616–237.7701.6 × 10−50.0400.453.1 × 10−6COPS8C0.191.13 (1.06–1.21)
rs503425a11118.079–118.1980.00123.3 × 10−40.433.3 × 10−6DDX6C0.201.16 (1.08–1.24)
rs10742326a1134.733–34.8091.4 × 10−40.0170.213.6 × 10−6APIPG0.591.14 (1.08–1.21)
rs4766921a12117.835–117.8834.6 × 10−5n.a.0.0364.6 × 10−6KIAA1853G0.671.18 (1.09–1.27)
rs11951576a56.741–6.8662.5 × 10−50.420.0144.6 × 10−6POLS-SRD5AC0.691.14 (1.08–1.22)
rs64387003123.355–123.4547.4 × 10−50.230.0205.5 × 10−6CD86C0.821.18 (1.09–1.27)
rs6486730a12127.830–127.8408.2 × 10−50.160.0496.9 × 10−6SLC15A4G0.411.13 (1.07–1.19)
rs4748857a1023.529–23.6542.2 × 10−40.681.3 × 10−46.9 × 10−6C10orf67C0.731.16 (1.09–1.24)
rs3914167a539.426–39.4541.8 × 10−40.240.00817.6 × 10−6DAB2-C9G0.271.15 (1.09–1.23)

Samples, critical region, P values, risk alleles and ORs are as defined in the Table 1 legend.

Indicates markers that were imputed, as described in Online Methods, from the GWAS samples and directly genotyped in the replication samples. See Supplementary Table 2 for expanded summary statistics.

A variant (rs7708392) on 5q33.1 that resides within an intron of TNIP1 (encoding TNF-α-induced protein 3 (TNFAIP3)-interacting protein 1) was significantly associated with SLE in all three cohorts and had a combined P = 3.8 × 10 (Fig. 2). Variants near TNIP1 were recently found to contribute to risk of psoriasis12; however, the SLE and psoriasis risk variants are separated by 21 kb and appear to have distinct genetic signals (r = 0.001). TNIP1 and TNFAIP3 are interacting proteins13, but the precise role of TNIP1 in regulating TNFAIP3 is unknown. The association of multiple distinct variants near TNFAIP3 with SLE414, rheumatoid arthritis15, psoriasis12 and type 1 diabetes16 suggests that this pathway has an important role in regulating autoimmunity.

A second confirmed risk variant (rs6568431, P = 7.12 × 10) was identified in an intergenic region between PRDM1 (PR domain containing 1, with ZNF domain, also known as BLIMP1) and ATG5 (APG5 autophagy 5-like). The signal at rs6568431 appears to be distinct from the previously reported6 SLE risk allele within ATG5 (rs2245214, Table 1), as rs6568431 has an r < 0.1 with rs2245214, and rs2245214 remains significantly associated with SLE (P < 1 × 10) after conditional logistic regression incorporating rs6568431 (Fig. 2).

The promoter region of JAZF1 (juxtaposed with another zinc finger gene 1) is a third newly confirmed SLE locus (rs849142, P = 1.54 × 10). Of note, this same variant was previously associated with risk of type 2 diabetes17 and with height variation18. A separate prostate cancer risk allele near JAZF1 (rs10486567)19 showed no evidence for association in the current study.

A fourth newly identified SLE risk locus is defined by a nonsynonymous allele (R454Q) of UHRF1BP1 (ICBP90 binding protein 1; rs11755393, P = 2.22 × 10). This allele encodes a nonconservative amino-acid change in a putative binding partner of UHRF1, a transcription and methylation factor linked to multiple pathways20. The UHRF1BP1 risk allele is in a region of extended LD that encompasses multiple genes, including SNRPC (small nuclear ribonucleoprotein polypeptide C), which is part of a RNA processing complex often targeted by SLE autoantibodies.

The fifth newly identified SLE locus is IL10 (interleukin-10; rs3024505, P = 3.95 × 10; Fig. 2). IL10 is an important immunoregulatory cytokine that functions to downregulate immune responses21, and variation in IL10 has inconsistently been reported to be associated with SLE22. The variant associated with SLE is identical to a SNP recently identified as contributing to risk of ulcerative colitis23 and type 1 diabetes24, suggesting the possibility of shared pathophysiology in the IL10 pathway across these disorders.

Using a significance threshold of P < 1 × 10 in the combined replication sample, we identified 21 additional SLE candidate risk loci (Table 2 and Supplementary Table 2). Less than one locus (0.01 loci, specifically) with P < 1 × 10 was expected under a null distribution for the meta-analysis (P = 8 × 10), suggesting that several of these loci are likely to be true positives for association to SLE. Notable candidate genes in this list include: (i) IRF8 (interferon regulatory factor 8), which was implicated in a previous GWA study (GWAS)4 and whose family members IRF5 and IRF7 are within confirmed SLE risk loci; (ii) TAOK3 (TAO kinase 3), a kinase expressed in lymphocytes, and the disease-associated variant is a missense allele (rs428073, N47S); (iii) LYST (lysosomal trafficking regulator), mutations of which cause Chediak-Higashi syndrome in humans, a complex disorder characterized by a lymphoproliferative disorder; and (iv) IL12RB2 (interleukin 12 receptor, beta 2), a locus which includes IL23R and SERPBP1 but appears distinct from the IL23R variants reported in inflammatory bowel disease, psoriasis and ankylosing spondylitis25.

A noteworthy feature of recent GWAS is the large number of overlapping loci found to be shared between different complex diseases26. We tested 42 variants from 35 loci that were previously reported as autoimmune-disease risk alleles for association with SLE (Table 3 and Supplementary Table 3). No single locus had an unadjusted P value < 5 × 10; however, we found an enrichment of previously identified disease-associated alleles. From the 35 loci tested (42 total variants), there were 5 alleles with unadjusted P < 0.0004 (less than 1 expected by chance, P = 4.4 × 10) and with P < 0.05 after a Bonferroni correction for the 35 pre-specified loci. For each of the five variants, the SLE-associated allele matched a previously reported allele and had the same direction of effect (Table 3). We observed a highly significant association to SLE of a missense allele of IFIH1 (rs1990760, P = 3.3 × 10) that has previously been associated with type 1 diabetes and Graves’ disease2728. We also observed an association of SLE with a missense allele (R32Q) of CFB (complement factor B, rs641153) that resides in the HLA class III region and is a validated risk allele for age-related macular degeneration29. This missense allele in CFB is not in significant LD with other HLA region variants associated with SLE (DR2/DR3) and remained significant (P < 0.05) after conditional logistic regression analyses that incorporated DR2 and DR3. The HLA is a complex genetic region, but it is noteworthy that the rs641153 SNP has a protective effect nearly identical to that of the reported age-related macular degeneration (AMD) risk allele29. Further validation of the five candidate disease alleles is required.

Table 3

Candidate autoimmune loci with evidence of association to SLE

SNPGeneChrP
Combined correctedRisk alleleRisk allele frequencyORPhenotypeReferences
GWASUSSwedenCombined
rs1990760IFIH123.2 × 10−50.0150.00393.34 × 10−71.12 × 10−5T0.601.17T1D, Graves’27,28
rs641153aCFB60.0079n.a.0.00111.4 × 10−40.0049G0.911.30AMD29
rs12708716aCLEC16A160.151.3 × 10−40.0621.6 × 10−40.0056A0.641.16T1D, Addison’s, MS3234
rs6887695aIL12B50.0140.040.031.7 × 10−40.0060G0.681.13Psoriasis, IBD12,35
rs17696736SH2B3120.00360.120.194.0 × 10−40.014T0.501.08T1D, Celiac, SLE33,36,37

All alleles in the table either were identical to the reported variants or have r > 0.8 to the reported variant and were the same risk allele with the same direction of effect. Samples, individual and combined P values, risk allele frequency and OR are as described in Table 1 legend. Combined-corrected P value is the Bonferroni-corrected P value for the 35 previously reported risk loci. Other autoimmunity associations: T1D, type 1 diabetes; AMD, age-related macular degeneration; MS, multiple sclerosis; IBD, inflammatory bowel disease. See Supplementary Table 3 for expanded summary statistics and a complete list of variants tested.

Indicates markers that were imputed, as described in Online Methods, from the GWAS samples and directly genotyped in the replication samples.

Using 26 SLE risk alleles (21 previously reported loci in Table 1 plus the 5 newly identified SLE loci), we performed several additional analyses. First, we performed pairwise interaction analysis with the previously confirmed loci, and, consistent with previous literature from SLE6 and other complex diseases30, we observed no evidence for non-additive interactions. Using conditional logistic regression analyses, we found no evidence for multiple independent alleles contributing to risk at any of the individual risk loci. We next estimated the percent of variance explained by each of the confirmed SLE risk alleles using previously described methods30. HLA-DRB1 (HLA*DR3), IRF5 and STAT4 were each estimated to account for over 1% of the genetic variance, whereas the remaining loci each accounted for less than 1% of the variance. Together, the 26 SLE risk loci explain an estimated 8% of the total genetic susceptibility to SLE.

Targeted replication of GWAS results is an efficient study design to confirm additional risk loci31. However, there are few available data as to the probability of replicating results that fall short of accepted P value criteria for genome-wide significance. In the current study, all variants with P < 0.05 from the original GWAS were included for replication. The lower a locus’ P value is in the GWAS, the higher is the probability of that locus reaching candidate or confirmed status in the replication meta-analysis (Fig. 3). Of note, no candidate or confirmed loci were obtained in the current study from the group of variants with a GWAS P value between 0.05 and 0.01, despite accounting for ~50% of all variants tested in the replication. These results may be useful in helping guide future targeted study designs, although clearly the size of the original GWAS population, the replication sample size, the disease architecture and the effect size of the candidate disease-associated variants need to be carefully considered in planning replication efforts.

An external file that holds a picture, illustration, etc.
Object name is nihms226372f3.jpg

Percentage of newly discovered variants reaching candidate (P < 1 × 10) and confirmed (P < 5 × 10) status in the meta-analysis stratified by the P value in the original GWAS.

These data provide further evidence that common variation in genes that function in the adaptive and innate arms of the immune system are important in establishing SLE risk. Although each of the identified alleles accounts for only a fraction of the overall genetic risk, these and other ongoing studies are providing new insight into the pathogenesis of lupus and are suggesting new targets and pathways for drug discovery and development.

Footnotes

Note: Supplementary information is available on the Nature Genetics website.

AUTHOR CONTRIBUTIONS

V.G. and J.K.S. performed the primary statistical analyses and contributed to initial manuscript preparation; J.K.S managed DNA samples and performed genotyping. G.H. contributed to the statistical analyses and experimental design. K.E.T. and S.A.C. performed statistical analyses and contributed to manuscript preparation. X.S., W.O. and R.C.F. managed DNA samples and contributed to experimental design. G.N., I.G., E.S., L.P., G.S., A.J., A.A.B., S.R.-D., E.C.B, E.E.B., G.S.A., J.C.E., R.R.-G., G.M. Jr., J.D.R., L.M.V., R.P.K., S.M. and M.A.P. provided samples and phenotype information. A.L. managed samples and oversaw genotyping efforts. P.K.G. provided samples and contributed to the initial manuscript preparation. M.F.S. and R.K. contributed statistical analyses and contributed to the selection of the ancestry-informative markers. L.R., L.A.C. and A.-C.S. contributed samples, input into experimental design, data interpretation and initial manuscript preparation; A.-C.S. oversaw genotyping efforts. R.R.G. and T.W.B. contributed to experimental design and interpretation, statistical analyses and initial manuscript preparation. All authors contributed to the final paper.

COMPETING INTERESTS STATEMENT

The authors declare competing financial interests: details accompany the full-text HTML version of the paper at http://www.nature.com/naturegenetics/.

Published online at http://www.nature.com/naturegenetics/.

Reprints and permissions information is available online at http://npg.nature.com/reprintsandpermissions/.

Footnotes

References

  • 1. Rönnblom L, Pascual VThe innate immune system in SLE: type I interferons and dendritic cells. Lupus. 2008;17:394–399.[Google Scholar]
  • 2. Vyse TJ, Todd JAGenetic analysis of autoimmune disease. Cell. 1996;85:311–318.[PubMed][Google Scholar]
  • 3. Cunninghame Graham DS, et al Polymorphism at the TNF superfamily gene OX40L confers susceptibility to systemic lupus erythematosus. Nat Genet. 2008;40:83–89.[Google Scholar]
  • 4. Graham RR, et al Genetic variants near TNFAIP3 on 6q23 are associated with systemic lupus erythematosus. Nat Genet. 2008;40:1059–1061.[Google Scholar]
  • 5. Graham RR, Hom G, Ortmann W, Behrens TWReview of recent genome-wide association scans in lupus. J Intern Med. 2009;265:680–688.[PubMed][Google Scholar]
  • 6. Harley JB, et al Genome-wide association scan in women with systemic lupus erythematosus identifies susceptibility variants in ITGAM, PXK, KIAA1542 and other loci. Nat Genet. 2008;40:204–210.[Google Scholar]
  • 7. Hom G, et al Association of systemic lupus erythematosus with C8orf13-BLK and ITGAM-ITGAX. N Engl J Med. 2008;358:900–909.[PubMed][Google Scholar]
  • 8. Kozyrev SV, et al Functional variants in the B-cell gene BANK1 are associated with systemic lupus erythematosus. Nat Genet. 2008;40:211–216. erratum40, 484 (2004) [[PubMed][Google Scholar]
  • 9. Sawalha AH, et al Common variants within MECP2 confer risk of systemic lupus erythematosus. PLoS ONE. 2008;3:e1727.[Google Scholar]
  • 10. Sigurdsson S, et al Polymorphisms in the tyrosine kinase 2 and interferon regulatory factor 5 genes are associated with systemic lupus erythematosus. Am J Hum Genet. 2005;76:528–537.[Google Scholar]
  • 11. Jacob CO, et al Identification of IRAK1 as a risk gene with critical role in the pathogenesis of systemic lupus erythematosus. Proc Natl Acad Sci USA. 2009;106:6256–6261.[Google Scholar]
  • 12. Nair RP, et al Genome-wide scan reveals association of psoriasis with IL-23 and NF-κB pathways. Nat Genet. 2009;41:199–204.[Google Scholar]
  • 13. Heyninck K, Kreike MM, Beyaert RStructure-function analysis of the A20-binding inhibitor of NF-κB activation, ABIN-1. FEBS Lett. 2003;536:135–140.[PubMed][Google Scholar]
  • 14. Musone SL, et al Multiple polymorphisms in the TNFAIP3 region are independently associated with systemic lupus erythematosus. Nat Genet. 2008;40:1062–1064.[Google Scholar]
  • 15. Plenge RM, et al Two independent alleles at 6q23 associated with risk of rheumatoid arthritis. Nat Genet. 2007;39:1477–1482.[Google Scholar]
  • 16. Fung EY, et al Analysis of 17 autoimmune disease-associated variants in type 1 diabetes identifies 6q23/TNFAIP3 as a susceptibility locus. Genes Immun. 2009;10:188–191.[PubMed][Google Scholar]
  • 17. Zeggini E, et al Meta-analysis of genome-wide association data and large-scale replication identifies additional susceptibility loci for type 2 diabetes. Nat Genet. 2008;40:638–645.[Google Scholar]
  • 18. Johansson A, et al Common variants in the JAZF1 gene associated with height identified by linkage and genome-wide association analysis. Hum Mol Genet. 2009;18:373–380.[Google Scholar]
  • 19. Thomas G, et al Multiple loci identified in a genome-wide association study of prostate cancer. Nat Genet. 2008;40:310–315.[PubMed][Google Scholar]
  • 20. Arita K, Ariyoshi M, Tochio H, Nakamura Y, Shirakawa MRecognition of hemi-methylated DNA by the SRA protein UHRF1 by a base-flipping mechanism. Nature. 2008;455:818–821.[PubMed][Google Scholar]
  • 21. Diveu C, McGeachy MJ, Cua DJCytokines that regulate autoimmunity. Curr Opin Immunol. 2008;20:663–668.[PubMed][Google Scholar]
  • 22. Nath SK, Harley JB, Lee YHPolymorphisms of complement receptor 1 and interleukin-10 genes and systemic lupus erythematosus: a meta-analysis. Hum Genet. 2005;118:225–234.[PubMed][Google Scholar]
  • 23. Franke A, et al Sequence variants in IL10, ARPC2 and multiple other loci contribute to ulcerative colitis susceptibility. Nat Genet. 2008;40:1319–1323.[PubMed][Google Scholar]
  • 24. Barrett JC, et al Genome-wide association study and meta-analysis find that over 40 loci affect risk of type 1 diabetes. Nat Genet. 2009;41:703–707.[Google Scholar]
  • 25. Duerr RH, et al A genome-wide association study identifies IL23R as an inflammatory bowel disease gene. Science. 2006;314:1461–1463.[Google Scholar]
  • 26. Zhernakova A, van Diemen CC, Wijmenga CDetecting shared pathogenesis from the shared genetics of immune-related diseases. Nat Rev Genet. 2009;10:43–55.[PubMed][Google Scholar]
  • 27. Smyth DJ, et al A genome-wide association study of nonsynonymous SNPs identifies a type 1 diabetes locus in the interferon-induced helicase (IFIH1) region. Nat Genet. 2006;38:617–619.[PubMed][Google Scholar]
  • 28. Sutherland A, et al Genomic polymorphism at the interferon-induced helicase (IFIH1) locus contributes to Graves’ disease susceptibility. J Clin Endocrinol Metab. 2007;92:3338–3341.[PubMed][Google Scholar]
  • 29. Gold B, et al Variation in factor B (BF) and complement component 2 (C2) genes is associated with age-related macular degeneration. Nat Genet. 2006;38:458–462.[Google Scholar]
  • 30. Barrett JC, et al Genome-wide association defines more than 30 distinct susceptibility loci for Crohn’s disease. Nat Genet. 2008;40:955–962.[Google Scholar]
  • 31. Hirschhorn JN, Daly MJGenome-wide association studies for common diseases and complex traits. Nat Rev Genet. 2005;6:95–108.[PubMed][Google Scholar]
  • 32. Awata T, et al Association of type 1 diabetes with two loci on 12q13 and 16p13 and the influence of coexisting thyroid autoimmunity in Japanese. J Clin Endocrinol Metab. 2009;94:231–235.[PubMed][Google Scholar]
  • 33. Skinningsrud B, et al Polymorphisms in CLEC16A and CIITA at 16p13 are associated with primary adrenal insufficiency. J Clin Endocrinol Metab. 2008;93:3310–3317.[PubMed][Google Scholar]
  • 34. Zoledziewska M, et al Variation within the CLEC16A gene shows consistent disease association with both multiple sclerosis and type 1 diabetes in Sardinia. Genes Immun. 2009;10:15–17.[PubMed][Google Scholar]
  • 35. Fisher SA, et al Genetic determinants of ulcerative colitis include the ECM1 locus and five loci implicated in Crohn’s disease. Nat Genet. 2008;40:710–712.[Google Scholar]
  • 36. Hunt KA, et al Newly identified genetic risk variants for celiac disease related to the immune response. Nat Genet. 2008;40:395–402.[Google Scholar]
  • 37. Smyth DJ, et al Shared and distinct genetic variants in type 1 diabetes and celiac disease. N Engl J Med. 2008;359:2767–2777.[Google Scholar]
Collaboration tool especially designed for Life Science professionals.Drag-and-drop any entity to your messages.