Genome-wide association (GWA) studies have reported susceptible regions in the human genome for many common diseases and traits, however, these loci only explain a minority of trait heritability. To boost the power of a GWA study, substantial research endeavors have been focused on integrating other available genomic information in the analysis. Advances in high through-put technologies have generated a wealth of genomic data, and made combining SNP and gene expression data become feasible.
In this paper we propose a novel procedure to incorporate gene expression information into GWA analysis. This procedure utilizes weights constructed by gene expression measurements to adjust p values from a GWA analysis. Results from simulation analyses indicate that the proposed procedures may achieve substantial power gains while controlling family-wise type I error rate (FWER) at the nominal level. To demonstrate the implementation of our proposed approach, we apply the weight adjustment procedure to a GWA study for serum interferon-regulated chemokine levels in systemic lupus erythematosus (SLE) patients. The study results can provide valuable insights for the functional interpretation of GWA signals.
The R source code for implementing the proposed weighting procedure is available at http://www.biostat.umn.edu/~yho/research.html
p value weighting; family-wise error rate; statistical power; integrative genomic analysis; SLE
Despite the success of genome-wide association studies (GWAS) in detecting a large number of loci for complex phenotypes such as rheumatoid arthritis (RA) susceptibility, the lack of information on the causal genes leaves important challenges to interpret GWAS results in the context of the disease biology. Here, we genetically fine-map the RA risk locus at 19p13 to define causal variants, and explore the pleiotropic effects of these same variants in other complex traits. First, we combined Immunochip dense genotyping (n = 23,092 case/control samples), Exomechip genotyping (n = 18,409 case/control samples) and targeted exon-sequencing (n = 2,236 case/controls samples) to demonstrate that three protein-coding variants in TYK2 (tyrosine kinase 2) independently protect against RA: P1104A (rs34536443, OR = 0.66, P = 2.3x10-21), A928V (rs35018800, OR = 0.53, P = 1.2x10-9), and I684S (rs12720356, OR = 0.86, P = 4.6x10-7). Second, we show that the same three TYK2 variants protect against systemic lupus erythematosus (SLE, Pomnibus = 6x10-18), and provide suggestive evidence that two of the TYK2 variants (P1104A and A928V) may also protect against inflammatory bowel disease (IBD; Pomnibus = 0.005). Finally, in a phenome-wide association study (PheWAS) assessing >500 phenotypes using electronic medical records (EMR) in >29,000 subjects, we found no convincing evidence for association of P1104A and A928V with complex phenotypes other than autoimmune diseases such as RA, SLE and IBD. Together, our results demonstrate the role of TYK2 in the pathogenesis of RA, SLE and IBD, and provide supporting evidence for TYK2 as a promising drug target for the treatment of autoimmune diseases.
A major challenge in human genetics is to devise a systematic strategy to integrate disease-associated variants with diverse genomic and biological datasets to provide insight into disease pathogenesis and guide drug discovery for complex traits such as rheumatoid arthritis (RA)1. Here, we performed a genome-wide association study (GWAS) meta-analysis in a total of >100,000 subjects of European and Asian ancestries (29,880 RA cases and 73,758 controls), by evaluating ~10 million single nucleotide polymorphisms (SNPs). We discovered 42 novel RA risk loci at a genome-wide level of significance, bringing the total to 1012–4. We devised an in-silico pipeline using established bioinformatics methods based on functional annotation5, cis-acting expression quantitative trait loci (cis-eQTL)6, and pathway analyses7–9 – as well as novel methods based on genetic overlap with human primary immunodeficiency (PID), hematological cancer somatic mutations and knock-out mouse phenotypes – to identify 98 biological candidate genes at these 101 risk loci. We demonstrate that these genes are the targets of approved therapies for RA, and further suggest that drugs approved for other indications may be repurposed for the treatment of RA. Together, this comprehensive genetic study sheds light on fundamental genes, pathways and cell types that contribute to RA pathogenesis, and provides empirical evidence that the genetics of RA can provide important information for drug discovery.
To investigate whether the two subtypes of advanced age-related macular degeneration (AMD), choroidal neovascularization (CNV) and geographic atrophy (GA), segregate separately in families and to identify which genetic variants are associated with these two subtypes.
Sibling correlation study and genome-wide association study (GWAS)
For the sibling correlation study, we included 209 sibling pairs with advanced AMD. For the GWAS, we included 2594 participants with advanced AMD subtypes and 4134 controls. Replication cohorts included 5383 advanced AMD participants and 15,240 controls.
Participants had AMD grade assigned based on fundus photography and/or examination. To determine heritability of advanced AMD subtypes, we performed a sibling correlation study. For the GWAS, we conducted genome-wide genotyping and imputed 6,036,699 single nucleotide polymorphism (SNPs). We then analyzed SNPs with a generalized linear model controlling for genotyping platform and genetic ancestry. The most significant associations were evaluated in independent cohorts.
Main Outcome Measures
Concordance of advanced AMD subtypes in sibling pairs and associations between SNPs with GA and CNV advanced AMD subtypes.
The difference between the observed and expected proportion of siblings concordant for the same subtype of advanced AMD was different to a statistically significant degree (P=4.2 x 10−5) meaning that siblings of probands with CNV or GA are more likely to develop CNV or GA, respectively. In the analysis comparing participants with CNV to those with GA, we observed a statistically significant association at the ARMS2/HTRA1 locus [rs10490924, odds ratio (OR)=1.47, P=4.3 ×10−9] which was confirmed in the replication samples (OR=1.38, P=7.4 x 10−14 for combined discovery and replication analysis).
Whether a patient with AMD develops CNV vs. GA is determined in part by genetic variation. In this large GWAS meta-analysis and replication analysis, the ARMS2/HTRA1 locus confers increased risk for both advanced AMD subtypes but imparts greater risk for CNV than for GA. This locus explains a small proportion of the excess sibling correlation for advanced AMD subtype. Other loci were detected with suggestive associations which differ for advanced AMD subtypes and deserve follow-up in additional studies.
Macrophage stimulating protein (MSP) is a serum growth factor that binds to and activates the receptor tyrosine kinase, Recepteur d'Origine Nantais (RON). A non-synonymous coding variant in MSP (689C) has been associated with genetic susceptibility to both Crohn's disease and ulcerative colitis, two major types of inflammatory bowel disease (IBD) characterized by chronic inflammation of the digestive tract. We investigated the consequences of this polymorphism for MSP-RON pathway activity and IBD pathogenesis.
RON expression patterns were examined on mouse and human cells and tissues under normal and disease conditions to identify cell types regulated by MSP-RON. Recombinant MSP variants were tested for their ability to bind and stimulate RON and undergo proteolytic activation. MSP concentrations were quantified in the serum of individuals carrying the MSP 689R and 689C alleles.
In intestinal tissue, RON was primarily expressed by epithelial cells under normal and disease conditions. The 689C polymorphism had no impact on the ability of MSP to bind to or signal through RON. In a cohort of normal individuals and IBD patients, carriers of the 689C polymorphism had lower concentrations of MSP in their serum.
By reducing the quantities of circulating MSP, the 689C polymorphism, or a variant in linkage disequilibrium with this polymorphism, may impact RON ligand availability and thus receptor activity. Given the known functions of RON in regulating wound healing and our analysis of RON expression patterns in human intestinal tissue, these data suggest that decreased RON activity may impact the efficiency of epithelial repair and thus underlie the increased IBD susceptibility associated with the MSP 689C allele.
African Americans, East Asians, and Hispanics with systemic lupus erythematosus (SLE) are more likely to develop renal disease than SLE patients of European descent. We investigated whether European genetic ancestry protects against the development of lupus nephritis and explored genetic and socioeconomic factors that might explain this effect.
This was a cross-sectional study of 1906 adults with SLE. Participants were genotyped for 126 single nucleotide polymorphisms (SNPs) informative for ancestry. A subset of participants was also genotyped for 80 SNPs in 14 candidate genes for renal disease in SLE. We used logistic regression to test the association between European ancestry and renal disease. Analyses adjusted for continental ancestries, socioeconomic status, and candidate genes.
Participants (n=1906) had on average 62.4% European, 15.8% African, 11.5% East Asian, 6.5% Amerindian, and 3.8% South Asian ancestry. Among participants, 34% (n=656) had renal disease. A 10% increase in European ancestry was associated with a 15% reduction in the odds of having renal disease after adjustment for disease duration and sex (OR 0.85, 95% CI 0.82-0.87, p=1.9 × 10−30). Adjusting for other genetic ancestries, measures of socioeconomic status, or SNPs in genes most associated with renal disease (IRF5 (rs4728142), BLK (rs2736340), STAT4 (rs3024912), ITGAM (rs9937837) and HLA-DRB1*0301 and DRB1*1501, p<0.05) did not substantively alter this relationship.
European ancestry is protective against the development of renal disease in SLE, an effect independent of other genetic ancestries, common risk alleles, and socioeconomic status.
Recent genome-wide association studies (GWASs) conducted in Asian populations have identified novel risk loci for systemic lupus erythematosus (SLE). Here, we genotyped 10 single-nucleotide polymorphisms (SNPs) in eight such loci and investigated their disease associations in three independent Caucasian SLE case–control cohorts recruited from Sweden, Finland and the United States. The disease associations of the SNPs in ETS1, IKZF1, LRRC18-WDFY4, RASGRP3, SLC15A4, TNIP1 and 16p11.2 were replicated, whereas no solid evidence of association was observed for the 7q11.23 locus in the Caucasian cohorts. SLC15A4 was significantly associated with renal involvement in SLE. The association of TNIP1 was more pronounced in SLE patients with renal and immunological disorder, which is corroborated by two previous studies in Asian cohorts. The effects of all the associated SNPs, either conferring risk for or being protective against SLE, were in the same direction in Caucasians and Asians. The magnitudes of the allelic effects for most of the SNPs were also comparable across different ethnic groups. On the contrary, remarkable differences in allele frequencies between Caucasian and Asian populations were observed for all associated SNPs. In conclusion, most of the novel SLE risk loci identified by GWASs in Asian populations were also associated with SLE in Caucasian populations. We observed both similarities and differences with respect to the effect sizes and risk allele frequencies across ethnicities.
systemic lupus erythematosus; genetic-association study; Asian; Caucasian
C-src tyrosine kinase, Csk, physically interacts with the intracellular phosphatase Lyp (PTPN22) and can modify the activation state of downstream Src kinases, such as Lyn, in lymphocytes. We identified an association of Csk with systemic lupus erythematosus (SLE) and refined its location to an intronic polymorphism rs34933034 (OR 1.32, p = 1.04 × 10−9). The risk allele is associated with increased CSK expression and augments inhibitory phosphorylation of Lyn. In carriers of the risk allele, B cell receptor (BCR)-mediated activation of mature B cells, as well as plasma IgM, are increased. Moreover, the fraction of transitional B cells is doubled in the cord blood of carriers of the risk allele compared to non-risk haplotypes due to an expansion of the late transitional cells, a stage targeted by selection mechanisms. This suggests that the Lyp-Csk complex increases susceptibility to lupus at multiple maturation and activation points of B cells.
Systemic lupus erythematosus (SLE) is a multisystem complex autoimmune disease of uncertain etiology (OMIM 152700). Over recent years a genetic component to SLE susceptibility has been established1–3. Recent successes with association studies in SLE have identified genes including IRF5 (refs. 4,5) and FCGR3B6. Two tumor necrosis factor (TNF) superfamily members located within intervals showing genetic linkage with SLE are TNFSF4 (also known as OX40L; 1q25), which is expressed on activated antigen-presenting cells (APCs)7,8 and vascular endothelial cells9, and also its unique receptor, TNFRSF4 (also known as OX40; 1p36), which is primarily expressed on activated CD4+ T cells10. TNFSF4 produces a potent co-stimulatory signal for activated CD4+ T cells after engagement of TNFRSF4 (ref. 11). Using both a family-based and a case-control study design, we show that the upstream region of TNFSF4 contains a single risk haplotype for SLE, which is correlated with increased expression of both cell-surface TNFSF4 and the TNFSF4 transcript. We hypothesize that increased expression of TNFSF4 predisposes to SLE either by quantitatively augmenting T cell–APC interaction or by influencing the functional consequences of T cell activation via TNFRSF4.
The overexpression of interferon (IFN)-inducible genes is a prominent feature of SLE, serves as a marker for active and more severe disease, and is also observed in other autoimmune and inflammatory conditions. The genetic variations responsible for sustained activation of IFN responsive genes are unknown.
We systematically evaluated association of SLE with a total of 1,754 IFN-pathway related genes, including IFN-inducible genes known to be differentially expressed in SLE patients and their direct regulators. We performed a three-stage design where two cohorts (total n=939 SLE cases, 3,398 controls) were analyzed independently and jointly for association with SLE, and the results were adjusted for the number of comparisons.
A total of 16,137 SNPs passed all quality control filters of which 316 demonstrated replicated association with SLE in both cohorts. Nine variants were further genotyped for confirmation in an average of 1,316 independent SLE cases and 3,215 independent controls. Association with SLE was confirmed for several genes, including the transmembrane receptor CD44 (rs507230, P = 3.98×10−12), cytokine pleiotrophin (PTN) (rs919581, P = 5.38×10−04), the heat-shock DNAJA1 (rs10971259, P = 6.31×10−03), and the nuclear import protein karyopherin alpha 1 (KPNA1) (rs6810306, P = 4.91×10−02).
This study expands the number of candidate genes associated with SLE and highlights the potential of pathway-based approaches for gene discovery. Identification of the causal alleles will help elucidate the molecular mechanisms responsible for activation of the IFN system in SLE.
Selective IgA deficiency (IgAD; serum IgA<0.07 g/l) is the most common form of human primary immune deficiency, affecting approximately 1∶600 individuals in populations of Northern European ancestry. The polygenic nature of IgAD is underscored by the recent identification of several new risk genes in a genome-wide association study. Among the characterized susceptibility loci, the association with specific HLA haplotypes represents the major genetic risk factor for IgAD. Despite the robust association, the nature and location of the causal variants in the HLA region remains unknown. To better characterize the association signal in this region, we performed a high-density SNP mapping of the HLA locus and imputed the genotypes of common HLA-B, -DRB1, and -DQB1 alleles in a combined sample of 772 IgAD patients and 1,976 matched controls from 3 independent European populations. We confirmed the complex nature of the association with the HLA locus, which is the result of multiple effects spanning the entire HLA region. The primary association signal mapped to the HLA-DQB1*02 allele in the HLA Class II region (combined P = 7.69×10−57; OR = 2.80) resulting from the combined independent effects of the HLA-B*0801-DRB1*0301-DQB1*02 and -DRB1*0701-DQB1*02 haplotypes, while additional secondary signals were associated with the DRB1*0102 (combined P = 5.86×10−17; OR = 4.28) and the DRB1*1501 (combined P = 2.24×10−35; OR = 0.13) alleles. Despite the strong population-specific frequencies of HLA alleles, we found a remarkable conservation of these effects regardless of the ethnic background, which supports the use of large multi-ethnic populations to characterize shared genetic association signals in the HLA region. We also provide evidence for the location of association signals within the specific extended haplotypes, which will guide future sequencing studies aimed at characterizing the precise functional variants contributing to disease pathogenesis.
The human leukocyte antigen (HLA) locus is robustly associated with many immune-mediated conditions. However, identification of the genetic variants contributing to the disease pathophysiology has been greatly hampered by the extensive chromosomal conservation within this genomic region. To better understand the association of the HLA locus in selective IgA deficiency (IgAD), we used an extensive genotyping database from a recent genome-wide association study (GWAS) to generate a high-density SNP map of this region in a combined sample of >2,700 individuals from 3 independent European populations. In addition, we took advantage of recent methodological advances to impute the more common HLA-B, -DRB1, and -DQB1 alleles in all subjects. We confirmed the strong disease-association of the HLA locus and identified several different signals located in specific conserved HLA haplotypes contributing independent risk or protection for IgAD. Further analysis of the chromosomal sequences associated with the associated HLA alleles allowed us to refine the mapping of the susceptibility variants. These findings represent the most comprehensive high-density SNP mapping of the HLA locus in IgAD to date and provide important new information as to the location of the genetic variants contributing to this common immune deficiency.
Infectious and inflammatory diseases have repeatedly shown strong genetic associations within the major histocompatibility complex (MHC); however, the basis for these associations remains elusive. To define host genetic effects on the outcome of a chronic viral infection, we performed genome-wide association analysis in a multiethnic cohort of HIV-1 controllers and progressors, and we analyzed the effects of individual amino acids within the classical human leukocyte antigen (HLA) proteins. We identified >300 genome-wide significant single-nucleotide polymorphisms (SNPs) within the MHC and none elsewhere. Specific amino acids in the HLA-B peptide binding groove, as well as an independent HLA-C effect, explain the SNP associations and reconcile both protective and risk HLA alleles. These results implicate the nature of the HLA–viral peptide interaction as the major factor modulating durable control of HIV infection.
In spite of the well-known clustering of multiple autoimmune disorders in families, analyses of specific shared genes and polymorphisms between systemic lupus erythematosus (SLE) and other autoimmune diseases (ADs) have been limited. Therefore, we comprehensively tested autoimmune variants for association with SLE, aiming to identify pleiotropic genetic associations between these diseases. We compiled a list of 446 non–Major Histocompatibility Complex (MHC) variants identified in genome-wide association studies (GWAS) of populations of European ancestry across 17 ADs. We then tested these variants in our combined Caucasian SLE cohorts of 1,500 cases and 5,706 controls. We tested a subset of these polymorphisms in an independent Caucasian replication cohort of 2,085 SLE cases and 2,854 controls, allowing the computation of a meta-analysis between all cohorts. We have uncovered novel shared SLE loci that passed multiple comparisons adjustment, including the VTCN1 (rs12046117, P = 2.02×10−06) region. We observed that the loci shared among the most ADs include IL23R, OLIG3/TNFAIP3, and IL2RA. Given the lack of a universal autoimmune risk locus outside of the MHC and variable specificities for different diseases, our data suggests partial pleiotropy among ADs. Hierarchical clustering of ADs suggested that the most genetically related ADs appear to be type 1 diabetes with rheumatoid arthritis and Crohn's disease with ulcerative colitis. These findings support a relatively distinct genetic susceptibility for SLE. For many of the shared GWAS autoimmune loci, we found no evidence for association with SLE, including IL23R. Also, several established SLE loci are apparently not associated with other ADs, including the ITGAM-ITGAX and TNFSF4 regions. This study represents the most comprehensive evaluation of shared autoimmune loci to date, supports a relatively distinct non–MHC genetic susceptibility for SLE, provides further evidence for previously and newly identified shared genes in SLE, and highlights the value of studies of potentially pleiotropic genes in autoimmune diseases.
It is well known that multiple autoimmune disorders cluster in families. However, all of the genetic variants that explain this clustering have not been discovered, and the specific genetic variants shared between systemic lupus erythematosus (SLE) and other autoimmune diseases (ADs) are not known. In order to better understand the genetic factors that explain this predisposition to autoimmunity, we performed a comprehensive evaluation of shared autoimmune genetic variants. First we considered results from 17 ADs and compiled a list with 446 significant genetic variants from these studies. We identified some genetic variants extensively shared between ADs, as well as the ADs that share the most variants. The genetic overlap between SLE and other ADs was modest. Next we tested how important all the 446 genetic variants were in our collection with a minimum of 1,500 SLE patients. Among the most significant variants in SLE, the majority had already been identified in previous studies, but we also discovered variants in two important immune genes. In summary, our data identified diseases with common genetic risk factors and novel SLE effects, and this supports a relatively distinct genetic susceptibility for SLE. This study helps delineate the genetic architecture of ADs.
Systemic lupus erythematosus (SLE) is a complex trait characterised by the production of a range of auto-antibodies and a diverse set of clinical phenotypes. Currently, ∼8% of the genetic contribution to SLE in Europeans is known, following publication of several moderate-sized genome-wide (GW) association studies, which identified loci with a strong effect (OR>1.3). In order to identify additional genes contributing to SLE susceptibility, we conducted a replication study in a UK dataset (870 cases, 5,551 controls) of 23 variants that showed moderate-risk for lupus in previous studies. Association analysis in the UK dataset and subsequent meta-analysis with the published data identified five SLE susceptibility genes reaching genome-wide levels of significance (Pcomb<5×10−8): NCF2 (Pcomb = 2.87×10−11), IKZF1 (Pcomb = 2.33×10−9), IRF8 (Pcomb = 1.24×10−8), IFIH1 (Pcomb = 1.63×10−8), and TYK2 (Pcomb = 3.88×10−8). Each of the five new loci identified here can be mapped into interferon signalling pathways, which are known to play a key role in the pathogenesis of SLE. These results increase the number of established susceptibility genes for lupus to ∼30 and validate the importance of using large datasets to confirm associations of loci which moderately increase the risk for disease.
Genome-wide association studies have revolutionised our ability to identify common susceptibility alleles for systemic lupus erythematosus (SLE). In complex diseases such as SLE, where many different genes make a modest contribution to disease susceptibility, it is necessary to perform large-scale association studies to combine results from several datasets, to have sufficient power to identify highly significant novel loci (P<5×10−8). Using a large SLE collection of 870 UK SLE cases and 5,551 UK unaffected individuals, we firstly replicated ten moderate-risk alleles (P<0.05) from a US–Swedish study of 3,273 SLE cases and 12,188 healthy controls. Combining our results with the US-Swedish data identified five new loci, which crossed the level for genome-wide significance: NCF2 (neutrophil cytosolic factor 2), IKZF1 (Ikaros family zinc-finger 1), IRF8 (interferon regulatory factor 8), IFIH1 (interferon-induced helicase C domain-containing protein 1), and TYK2 (tyrosine kinase 2). Each of these five genes regulates a different aspect of the immune response and contributes to the production of type-I and type-II interferons. Although further studies will be required to identify the causal alleles within these loci, the confirmation of five new susceptibility genes for lupus makes a significant step forward in our understanding of the genetic contribution to SLE.
Genome-wide association (GWA) studies have identified numerous, replicable, genetic associations between common single nucleotide polymorphisms (SNPs) and risk of common autoimmune and inflammatory (immune-mediated) diseases, some of which are shared between two diseases. Along with epidemiological and clinical evidence, this suggests that some genetic risk factors may be shared across diseases—as is the case with alleles in the Major Histocompatibility Locus. In this work we evaluate the extent of this sharing for 107 immune disease-risk SNPs in seven diseases: celiac disease, Crohn's disease, multiple sclerosis, psoriasis, rheumatoid arthritis, systemic lupus erythematosus, and type 1 diabetes. We have developed a novel statistic for Cross Phenotype Meta-Analysis (CPMA) which detects association of a SNP to multiple, but not necessarily all, phenotypes. With it, we find evidence that 47/107 (44%) immune-mediated disease risk SNPs are associated to multiple—but not all—immune-mediated diseases (SNP-wise PCPMA<0.01). We also show that distinct groups of interacting proteins are encoded near SNPs which predispose to the same subsets of diseases; we propose these as the mechanistic basis of shared disease risk. We are thus able to leverage genetic data across diseases to construct biological hypotheses about the underlying mechanism of pathogenesis.
Over the last five years we have found over 100 genetic variants predisposing to common diseases affecting the immune system. In this study we analyze 107 such variants across seven diseases and find that almost half are shared across diseases. We also find that the patterns of sharing across diseases cluster these variants into groups; proteins encoded near variants in the same group tend to interact. This suggests that genetic variation may influence entire pathways to create risk to multiple diseases.
Despite significant progress in the identification of genetic loci for age-related macular degeneration (AMD), not all of the heritability has been explained. To identify variants which contribute to the remaining genetic susceptibility, we performed the largest meta-analysis of genome-wide association studies to date for advanced AMD. We imputed 6 036 699 single-nucleotide polymorphisms with the 1000 Genomes Project reference genotypes on 2594 cases and 4134 controls with follow-up replication of top signals in 5640 cases and 52 174 controls. We identified two new common susceptibility alleles, rs1999930 on 6q21-q22.3 near FRK/COL10A1 [odds ratio (OR) 0.87; P = 1.1 × 10−8] and rs4711751 on 6p12 near VEGFA (OR 1.15; P = 8.7 × 10−9). In addition to the two novel loci, 10 previously reported loci in ARMS2/HTRA1 (rs10490924), CFH (rs1061170, and rs1410996), CFB (rs641153), C3 (rs2230199), C2 (rs9332739), CFI (rs10033900), LIPC (rs10468017), TIMP3 (rs9621532) and CETP (rs3764261) were confirmed with genome-wide significant signals in this large study. Loci in the recently reported genes ABCA1 and COL8A1 were also detected with suggestive evidence of association with advanced AMD. The novel variants identified in this study suggest that angiogenesis (VEGFA) and extracellular collagen matrix (FRK/COL10A1) pathways contribute to the development of advanced AMD.
Systemic Lupus Erythematosus (SLE) is a systemic autoimmune disease in which the type I interferon pathway has a crucial role. We have previously shown that three genes in this pathway, IRF5, TYK2 and STAT4, are strongly associated with risk for SLE. Here, we investigated 78 genes involved in the type I interferon pathway to identify additional SLE susceptibility loci. First, we genotyped 896 single-nucleotide polymorphisms in these 78 genes and 14 other candidate genes in 482 Swedish SLE patients and 536 controls. Genes with P<0.01 in the initial screen were then followed up in 344 additional Swedish patients and 1299 controls. SNPs in the IKBKE, TANK, STAT1, IL8 and TRAF6 genes gave nominal signals of association with SLE in this extended Swedish cohort. To replicate these findings we extracted data from a genomewide association study on SLE performed in a US cohort. Combined analysis of the Swedish and US data, comprising a total of 2136 cases and 9694 controls, implicates IKBKE and IL8 as SLE susceptibility loci (Pmeta=0.00010 and Pmeta=0.00040, respectively). STAT1 was also associated with SLE in this cohort (Pmeta=3.3 × 10−5), but this association signal appears to be dependent of that previously reported for the neighbouring STAT4 gene. Our study suggests additional genes from the type I interferon system in SLE, and highlights genes in this pathway for further functional analysis.
systemic lupus erythematosus; type I interferon system; candidate gene study; single nucleotide polymorphism; IKBKE; IL8
Systemic lupus erythematosus (SLE) is a clinically heterogeneous, systemic autoimmune disease characterized by autoantibody formation. Previously published genome-wide association studies (GWAS) have investigated SLE as a single phenotype. Therefore, we conducted a GWAS to identify genetic factors associated with anti–dsDNA autoantibody production, a SLE–related autoantibody with diagnostic and clinical importance. Using two independent datasets, over 400,000 single nucleotide polymorphisms (SNPs) were studied in a total of 1,717 SLE cases and 4,813 healthy controls. Anti–dsDNA autoantibody positive (anti–dsDNA +, n = 811) and anti–dsDNA autoantibody negative (anti–dsDNA –, n = 906) SLE cases were compared to healthy controls and to each other to identify SNPs associated specifically with these SLE subtypes. SNPs in the previously identified SLE susceptibility loci STAT4, IRF5, ITGAM, and the major histocompatibility complex were strongly associated with anti–dsDNA + SLE. Far fewer and weaker associations were observed for anti–dsDNA – SLE. For example, rs7574865 in STAT4 had an OR for anti–dsDNA + SLE of 1.77 (95% CI 1.57–1.99, p = 2.0E-20) compared to an OR for anti–dsDNA – SLE of 1.26 (95% CI 1.12–1.41, p = 2.4E-04), with pheterogeneity<0.0005. SNPs in the SLE susceptibility loci BANK1, KIAA1542, and UBE2L3 showed evidence of association with anti–dsDNA + SLE and were not associated with anti–dsDNA – SLE. In conclusion, we identified differential genetic associations with SLE based on anti–dsDNA autoantibody production. Many previously identified SLE susceptibility loci may confer disease risk through their role in autoantibody production and be more accurately described as autoantibody propensity loci. Lack of strong SNP associations may suggest that other types of genetic variation or non-genetic factors such as environmental exposures have a greater impact on susceptibility to anti–dsDNA – SLE.
Systemic lupus erythematosus (SLE) is a chronic autoimmune disease that can involve virtually any organ system. SLE patients produce antibodies that bind to their own cells and proteins (autoantibodies) which can cause irreversible organ damage. One particular SLE–related autoantibody directed at double-stranded DNA (anti–dsDNA) is associated with kidney involvement and more severe disease. Previous genome-wide association studies (GWAS) in SLE have studied SLE itself, not particular SLE manifestations. Therefore, we conducted this GWAS of anti–dsDNA autoantibody production to identify genetic associations with this clinically important autoantibody. We found that many previously identified SLE–associated genes are more strongly associated with anti–dsDNA autoantibody production than SLE itself, and they may be more accurately described as autoantibody propensity genes. No strong genetic associations were observed for SLE patients who do not produce anti–dsDNA autoantibodies, suggesting that other factors may have more influence in developing this type of SLE. Further investigation of these autoantibody propensity genes may lead to greater insight into the causes of autoantibody production and organ damage in SLE.
Systemic lupus erythematosus (SLE) is a systemic autoimmune disease characterized by the production of antinuclear autoantibodies and the inflammatory infiltration of many organ systems. SLE is a complex disorder in which multiple genetic variants, together with environmental and hormonal factors, contribute to disease risk. In this article, we summarize our current understanding of the genetic contribution to SLE in light of recent genome-wide association studies, which have brought the total number of confirmed SLE susceptibility loci to 29. In the second section, we explore the functional implications of these risk loci and, in particular, highlight the role that many of these genes play in the Toll-like receptor and type I interferon signaling pathways. Finally, we discuss the genetic overlap between SLE and other autoimmune and inflammatory conditions as several risk loci are shared among multiple disorders, suggesting common underlying pathogenic mechanisms.
genome-wide association studies; pathogenesis; systemic lupus erythematosus; Toll-like receptor signaling; type I interferon signaling
Systemic lupus erythematosus (SLE) is a genetically complex disease with heterogeneous clinical manifestations. Recent studies have greatly expanded the number of established SLE risk alleles, but the distribution of multiple risk alleles in cases versus controls and their relationship to subphenotypes have not been studied. We studied 22 SLE susceptibility polymorphisms with previous genome-wide evidence of association (p<5×10−8) in 1919 SLE cases from 9 independent Caucasian SLE case series and 4813 independent controls. The mean number of risk alleles in cases was 15.1 (SD 3.1) while the mean in controls was 13.1 (SD 2.8), with trend p = 4×10−128. We defined a genetic risk score (GRS) for SLE as the number of risk alleles with each weighted by the SLE risk odds ratio (OR). The OR for high-low GRS tertiles, adjusted for intra-European ancestry, sex, and parent study, was 4.4 (95% CI 3.8–5.1). We studied associations of individual SNPs and the GRS with clinical manifestations for the cases: age at diagnosis, the 11 American College of Rheumatology classification criteria, and double-stranded DNA antibody (anti-dsDNA) production. Six subphenotypes were significantly associated with the GRS, most notably anti-dsDNA (ORhigh-low = 2.36, p = 9e−9), the immunologic criterion (ORhigh-low = 2.23, p = 3e−7), and age at diagnosis (ORhigh-low = 1.45, p = 0.0060). Finally, we developed a subphenotype-specific GRS (sub-GRS) for each phenotype with more power to detect cumulative genetic associations. The sub-GRS was more strongly associated than any single SNP effect for 5 subphenotypes (the above plus hematologic disorder and oral ulcers), while single loci are more significantly associated with renal disease (HLA-DRB1, OR = 1.37, 95% CI 1.14–1.64) and arthritis (ITGAM, OR = 0.72, 95% CI 0.59–0.88). We did not observe significant associations for other subphenotypes, for individual loci or the sub-GRS. Thus our analysis categorizes SLE subphenotypes into three groups: those having cumulative, single, and no known genetic association with respect to the currently established SLE risk loci.
Systemic lupus erythematosus is a chronic disabling autoimmune disease, most commonly striking women in their thirties or forties. It can cause a wide variety of clinical manifestations, including kidney disease, arthritis, and skin disorders. Prognosis varies greatly depending on these clinical features, with kidney disease and related characteristics leading to greater morbidity and mortality. It is also complex genetically; while lupus runs in families, genes increase one's risk for lupus but do not fully determine the outcome. The interactions of multiple genes and/or interactions between genes and environmental factors may cause lupus, but the causes and disease pathways of this very heterogeneous disease are not well understood. By examining relationships between the presence of multiple lupus risk genes, lupus susceptibility, and clinical manifestations, we hope to better understand how lupus is triggered and by what biological pathways it progresses. We show in this work that certain clinical manifestations of lupus are highly associated with cumulative genetic variations, i.e. multiple risk alleles, while others are associated with a single variation or none at all.
Systemic lupus erythematosus (SLE) is an autoimmune disease which behaves as a complex genetic trait. At least 20 SLE risk susceptibility loci have been mapped using both candidate gene and genome-wide association strategies. The gene encoding the pro-inflammatory cytokine, IL18, has been reported as a candidate gene showing an association with SLE. This pleiotropic cytokine is expressed in a range of immune cells and has been shown to induce interferon-γ and tumour necrosis factor-α. Serum interleukin-18 has been reported to be elevated in patients with SLE. Here we aimed to densely map single nucleotide polymorphisms (SNPs) across IL18 to investigate the association across this locus. We genotyped 36 across IL18 by Illumina bead express in 372 UK SLE trios. We also genotyped these SNPs in a further 508 non-trio UK cases and were able to accurately impute a dense marker set across IL18 in WTCCC2 controls with a total of 258 SNPs. To improve the study's power, we also imputed a total of 158 SNPs across the IL18 locus using data from an SLE genome-wide association study and performed association testing. In total, we analysed 1818 cases and 10 770 controls in this study. Our large well-powered study (98% to detect odds ratio = 1.5, with respect to rs360719) showed that no individual SNP or haplotype was associated with SLE in any of the cohorts studied. We conclude that we were unable to replicate the SLE association with rs360719 located upstream of IL18. No evidence for association with any other common variant at IL18 with SLE was found.
Prior evidence demonstrates that the PRL -1149 T (minor) allele decreases prolactin expression and may be associated with autoimmune diseases. Our goal was to determine the role of the PRL -1149 G/T polymorphism (rs1341239) in rheumatoid arthritis (RA) susceptibility.
We examined the association between PRL -1149 G/T and RA risk in four separate collections consisting of 3,405 RA cases and 4,111 controls of self-described “white” European ancestry. Samples were genotyped using one of three genotyping platforms, and strict quality control metrics were applied. We tested for association with a two-tailed Cochran Mantel-Haenszel additive, fixed effects model.
In the individual collections, odds ratios for an association between PRL -1149 T and RA risk ranged from 0.80 to 0.97. In a joint meta-analysis across all four collections, the odds ratio for an association between PRL -1149 T and RA risk was 0.90 (95% CI 0.84-0.96) with p = 0.001.
Our results suggest an association between the PRL -1149 T allele and decreased RA risk. The effect size is small but similar to odds ratios for other genetic polymorphisms associated with complex traits, including RA. Further studies will be necessary to confirm association unequivocally.
TNFAIP3 interacting protein 1, TNIP1 (ABIN-1) is involved in inhibition of nuclear factor-κB (NF-κB) activation by interacting with TNF alpha-induced protein 3, A20 (TNFAIP3), an established susceptibility gene to systemic lupus erythematosus (SLE) and rheumatoid arthritis (RA). Recent genome-wide association studies revealed association of TNIP1 with SLE in the Caucasian and Chinese populations. In this study, we investigated whether the association of TNIP1 with SLE was replicated in a Japanese population. In addition, association of TNIP1 with RA was also examined.
A case-control association study was conducted on the TNIP1 single nucleotide polymorphism (SNP) rs7708392 in 364 Japanese SLE patients, 553 RA patients and 513 healthy controls.
Association of TNIP1 rs7708392C was replicated in Japanese SLE (allele frequency in SLE: 76.5%, control: 69.9%, P = 0.0022, odds ratio [OR] 1.40, 95% confidence interval [CI] 1.13-1.74). Notably, the risk allele frequency in the healthy controls was considerably greater in Japanese (69.9%) than in Caucasians (24.3%). A tendency of stronger association was observed in the SLE patients with renal disorder (P = 0.00065, OR 1.60 [95%CI 1.22-2.10]) than in all SLE patients (P = 0.0022, OR 1.40 [95%CI 1.13-1.74]). Significant association with RA was not observed, regardless of the carriage of human leukocyte antigen DR β1 (HLA-DRB1) shared epitope. Significant gene-gene interaction between TNIP1 and TNFAIP3 was detected neither in SLE nor RA.
Association of TNIP1 with SLE was confirmed in a Japanese population. TNIP1 is a shared SLE susceptibility gene in the Caucasian and Asian populations, but the genetic contribution appeared to be greater in the Japanese and Chinese populations because of the higher risk allele frequency. Taken together with the association of TNFAIP3, these observations underscore the crucial role of NF-κB regulation in the pathogenesis of SLE.
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−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 ≤ 1 × 10−5. A candidate screen of alleles previously associated with other autoimmune diseases suggested five loci (P < 1 × 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.
To determine whether genetic substructure in European-derived populations is associated with specific manifestations of systemic lupus erythematosus (SLE), including mucocutaneous phenotypes, autoantibody production, and renal disease.
SLE patients of European descent (n=1754) from 8 case collections were genotyped for over 1,400 ancestry informative markers that define a north/south gradient of European substructure. Based on these genetic markers, we used the STRUCTURE program to characterize each SLE patient in terms of percent northern (vs. southern) European ancestry. Non-parametric methods, including tests of trend, were used to identify associations between northern European ancestry and specific SLE manifestations.
In multivariate analyses, increasing levels of northern European ancestry were significantly associated with photosensitivity (ptrend=0.0021, OR for highest quartile of northern European ancestry compared to lowest quartile 1.64, 95% CI 1.13–2.35) and discoid rash (ptrend=0.014, ORhigh-low 1.93, 95% CI 0.98–3.83). In contrast, northern European ancestry was protective for anticardiolipin (ptrend=1.6 × 10−4, ORhigh-low 0.46, 95% CI 0.30–0.69) and anti-dsDNA (ptrend=0.017, ORhigh-low 0.67, 95% CI 0.46–0.96) autoantibody production.
This study demonstrates that specific SLE manifestations vary according to northern vs. southern European ancestry. Thus, genetic ancestry may contribute to the clinical heterogeneity and variation in disease outcomes among SLE patients of European descent. Moreover, these results suggest that genetic studies of SLE subphenotypes will need to carefully address issues of population substructure due to genetic ancestry.