Trait-stratified genome-wide association study identifies novel and diverse genetic associations with serologic and cytokine phenotypes in systemic lupus erythematosus
© Kariuki et al.; licensee BioMed Central Ltd. 2010
Received: 7 January 2009
Accepted: 26 July 2010
Published: 26 July 2010
Systemic lupus erythematosus (SLE) is a highly heterogeneous disorder, characterized by differences in autoantibody profile, serum cytokines, and clinical manifestations. SLE-associated autoantibodies and high serum interferon alpha (IFN-α) are important heritable phenotypes in SLE which are correlated with each other, and play a role in disease pathogenesis. These two heritable risk factors are shared between ancestral backgrounds. The aim of the study was to detect genetic factors associated with autoantibody profiles and serum IFN-α in SLE.
We undertook a case-case genome-wide association study of SLE patients stratified by ancestry and extremes of phenotype in serology and serum IFN-α. Single nucleotide polymorphisms (SNPs) in seven loci were selected for follow-up in a large independent cohort of 538 SLE patients and 522 controls using a multi-step screening approach based on novel metrics and expert database review. The seven loci were: leucine-rich repeat containing 20 (LRRC20); protein phosphatase 1 H (PPM1H); lysophosphatidic acid receptor 1 (LPAR1); ankyrin repeat and sterile alpha motif domain 1A (ANKS1A); protein tyrosine phosphatase, receptor type M (PTPRM); ephrin A5 (EFNA5); and V-set and immunoglobulin domain containing 2 (VSIG2).
SNPs in the LRRC20, PPM1H, LPAR1, ANKS1A, and VSIG2 loci each demonstrated strong association with a particular serologic profile (all odds ratios > 2.2 and P < 3.5 × 10-4). Each of these serologic profiles was associated with increased serum IFN-α. SNPs in both PTPRM and LRRC20 were associated with increased serum IFN-α independent of serologic profile (P = 2.2 × 10-6 and P = 2.6 × 10-3 respectively). None of the SNPs were strongly associated with SLE in case-control analysis, suggesting that the major impact of these variants will be upon subphenotypes in SLE.
This study demonstrates the power of using serologic and cytokine subphenotypes to elucidate genetic factors involved in complex autoimmune disease. The distinct associations observed emphasize the heterogeneity of molecular pathogenesis in SLE, and the need for stratification by subphenotypes in genetic studies. We hypothesize that these genetic variants play a role in disease manifestations and severity in SLE.
Systemic lupus erythematosus ((SLE) OMIM #152700) is a complex disease characterized by multi-system involvement commonly affecting the skin, renal, musculoskeletal, and hematopoetic systems. SLE is caused by interactions between susceptibility genes and environmental factors resulting in an irreversible loss of immunologic self-tolerance . Incidence is highest in women during the reproductive years ; however, people of all ages, genders, and ancestral backgrounds are susceptible. Disease features range from typically reversible manifestations such as rash or inflammatory arthritis to life-threatening end-organ damage such as glomerulonephritis or thrombosis, and it is difficult to predict which manifestations will affect an individual patient.
Interferon alpha (IFN-α) is a pleiotropic type I interferon with the potential to break immunologic self-tolerance by activating antigen-presenting cells after uptake of self material . Serum IFN-α is elevated in many SLE patients, and elevations often correlate with disease activity [4, 5]. Recombinant human IFN-α administered as a therapy for chronic viral hepatitis and malignancy is thought to cause de novo SLE in some patients . IFN-α-induced SLE typically resolves after the IFN-α is discontinued [7, 8], supporting the idea that IFN-α was causal. We have previously shown that serum IFN-α is abnormally high in 20% of healthy first degree relatives of SLE patients as compared to < 5% of healthy unrelated individuals . Spouses of SLE patients did not have high serum IFN-α. Taken together, these data suggest that high serum IFN-α is a heritable risk factor for SLE. Additionally, serum IFN-α activity is highest during the ages of peak SLE incidence in both patients and their healthy first degree relatives . The high IFN-α trait in SLE families is inherited in a complex fashion, suggesting polygenic inheritance which has not been fully characterized .
Autoantibodies directed at double-stranded DNA (dsDNA) and RNA binding proteins (anti-Ro, anti-La, anti-Sm, and anti-RNP, collectively anti-RBP) are characteristically found in SLE sera, and are the strongest known predictors of high serum IFN-α in SLE patients . The presence of anti-RBP and anti-dsDNA autoantibodies appears to be heritable in SLE families . These autoantibodies are rare in healthy first degree relatives (one to three percent prevalence), and thus heritability of autoantibody traits is not sufficient to explain the heritability of high serum IFN-α in SLE families . Immune complexes formed by these SLE-associated autoantibodies can directly stimulate IFN-α production in vitro, likely via the endosomal Toll-like receptors . In humans in vivo, we find that while anti-RBP antibodies are associated with high serum IFN-α, they are not sufficient . This suggests that although these SLE-associated autoantibodies are correlated and likely mechanistically related to serum IFN-α, other factors such as genetic variation will likely be important to in vivo IFN-α levels.
We have begun to genetically map the high IFN-α trait in SLE, and found that a number of genetic variants which are associated with susceptibility to SLE are also associated with increased serum IFN-α, including variants of interferon regulatory factor 5 (IRF5) , interferon regulatory factor 7 (IRF7) , protein tyrosine phosphatase non receptor type 22 (PTPN22) , and osteopontin (SPP1) . We have also shown that the autoimmune disease associated variant of signal transducer and activator of transcription 4 (STAT4) modulates cellular sensitivity to IFN-α signaling . These studies provide further support for the hypothesis that IFN-α pathway dysregulation is a primary causal factor in human SLE. This work also suggests that the serum IFN-α trait could be a useful tool in identifying additional genetic loci related to SLE pathogenesis.
In this study, we performed a genome-wide association study (GWAS) on 104 SLE patients stratified by extremes of phenotype in serum IFN-α. Given the importance of anti-RBP and anti-dsDNA antibodies on serum IFN-α, patients were also stratified by the presence or absence of these antibodies. The size of the GWAS study did not allow for rigorous exclusion of false positives, so we used a multi-step screening algorithm which included novel metrics derived from the GWAS chip data and expert literature and database review to increase the likelihood of selecting true positive SNPs for validation. SNPs selected were then genotyped in an independent SLE cohort, testing for association with autoantibody profile and serum IFN-α in a case-case study design, as well as overall association with SLE in a case-control study. This approach was chosen to maximize our opportunity for novel gene discovery in a moderate-sized well-characterized cohort with complete molecular phenotype data. The trait-stratified GWAS screening approach followed by selected validation was successful, and we report seven novel loci which are strongly associated with autoantibodies and serum IFN-α, two important heritable subphenotypes in SLE.
Materials and methods
Patients and samples
The initial cohort of SLE patients studied in the GWAS scan was obtained from the Hospital for Special Surgery Lupus Registry, and consisted of 104 SLE patients, including 20 African-American, 36 European-American, 32 Hispanic-American, and 16 Asian-American SLE patients. Patients were selected from the top 33% and bottom 33% of serum IFN-α activity, and were additionally stratified by ancestry and the presence or absence of anti-RBP or anti-dsDNA antibodies. The independent validation cohort of 538 SLE patients was obtained from the University of Chicago Translational Research in the Department of Medicine (TRIDOM) registry and Rush University Medical Center, and consisted of 280 African-American, 173 European-American, and 85 Hispanic-American SLE patients. The validation cohort was not selected for inclusion in the study based upon any particular SLE-associated phenotype. All patients met the revised 1982 American College of Rheumatology criteria for the diagnosis of SLE . Samples from 522 controls were obtained from the TRIDOM registry, including 361 African-American and 161 European-American subjects who were individually screened for the absence of autoimmune disease by medical record review and matched by sex to the cases. The controls were on average older than the cases, although given our ability to screen for the absence of autoimmune disease this should be of benefit as one would assume that older individuals who do not have autoimmune disease are likely a low-risk group for autoimmunity. Clinical characteristics and demographic information for the patients and controls in the validation cohort are summarized in Supplemental table S1 in Additional file 1. Informed consent was obtained from all subjects at each site, and the study was approved by the institutional review board at each institution.
Measurement of serum IFN-α activity and autoantibodies
The reporter cell assay for IFN-α has been described in detail elsewhere [9, 21]. In brief, reporter cells (WISH cells, ATCC #CCL-25) were used to measure the ability of patient sera to cause IFN-induced gene expression. The output of the assay is interpreted as a functional activity of IFN-α in the serum sample. Results from the IFN-α assay were standardized to a healthy multi-ancestral reference population as previously described, and a serum IFN-α activity score was calculated based upon the mean and SD of the reference population . Antibodies to anti-Ro, anti-La, anti-Sm, and anti-RNP were measured in all samples by ELISA methods, and anti-dsDNA antibodies were measured using Crithidia luciliae immunofluorescence, with detectable fluorescence considered positive. Autoantibodies were assayed in the TRIDOM and Rush University samples in the University of Chicago clinical laboratory, and the Hospital for Special Surgery samples used in the GWAS screening were assayed in the Hospital for Special Surgery clinical laboratory. Standard clinical lab cutoff points were used to categorize samples as positive or negative. There were no significant differences in autoantibody prevalence between the two clinical laboratories.
Screening GWAS study
The initial cohort of SLE patients studied in the GWAS scan consisted of 104 SLE patients obtained from the Hospital for Special Surgery Registry. A study design incorporating multiple ancestral backgrounds was chosen, as both autoantibodies and serum IFN-α levels are heritable pathogenic factors which are shared between all ancestral backgrounds [9, 12]. This does not ensure that the same factors will underlie the serum IFN-α trait in all backgrounds, and we considered each ancestral background separately in the GWAS. Patients were selected in an "extremes of phenotype" design from the top 33% and bottom 33% of serum IFN-α activity, and were additionally stratified by ancestry and the presence or absence of anti-RBP or anti-dsDNA antibodies. Genomic DNA from SLE patients in the HSS registry was combined into pools of four individuals matched by ancestry, IFN-α levels, and presence or absence of anti-RBP or anti-dsDNA autoantibodies, respectively. The high IFN-α group was defined as those patients in the top 33% of serum IFN-α who had one or more of the following autoantibody specificities: Ro, La, Sm, RNP, or dsDNA, and the low group was defined as those in the bottom 33% of serum IFN-α who lacked these autoantibodies. These pools containing equal amounts of DNA from four individuals matched by these traits were run on Affymetrix 500K SNP mapping chips (Santa Clara, CA, USA) per manufacturer protocol. Data from the individual probes were then analyzed, and probe intensity was used to represent allelic frequencies present in the pooled DNA. Probe data were summarized and analyzed using Gene Pool by Pearson JV et al. . We considered two metrics to help identify allele frequency differences between the high vs. low IFN-α groups and the present vs. absent anti-RBP or anti-dsDNA groups. The first metric is the Silhouette score proposed by Pearson et al. .
We also developed our own score which is the difference in the weighted mean of the relative allele signals (RAS = intensity of A allele/(intensity of A + intensity of B allele)) for each ethnic group (EA, AA, HI). The weights were constructed as the inverse standard deviation of the RAS score calculated from the 8 (or 12) probes per SNP present on the Affymetrix array. The weights were then normalized so that they sum to one. Our adjusted RAS score is then calculated by summing each group's weighted RAS score. The goal of constructing our own statistic was three-fold: First, it down-weighted SNPs that had highly variable signals even if the mean differences between groups was large; this allows a certain measure of quality control that is not offered as effectively by the Silhouette score; second, it accounts for heterogeneity of allele frequencies among the ancestral groups; third, it creates a statistic that is robust to population stratification which can lead to false positives.
Specifically, the weighted mean RAS score is calculated as follows:
Of the 104 subjects studied in the GWAS, 44 were in the high IFN-α/autoantibody positive group and 60 were in the low IFN-α/autoantibody negative group. It was difficult to match even numbers of subjects between these two groups in each ancestral background due to the stringent double phenotype case selection criteria, and some of the ancestral groups deviated from a 1:1 match (proportions by ancestral background were as follows: Asian-American 8 high and 8 low, African-American 12 high and 8 low, European-American 12 high and 24 low, Hispanic-American 12 high and 20 low). This was not confounding, as the data from each ancestral background were analyzed separately.
Selection of SNPs from GWAS study for validation
We next employed a multiple step screening algorithm to select a subset of SNPs for replication that were highly suggestive based on both their Silhouette and weighted RAS scores. SNPs were first ordered by weighted RAS score, and the top 200 SNPs by RAS score were considered. SNPs on this list with a Silhouette score that was not within the top 200 Silhouette scores were removed (31%), leaving SNPs that ranked in the top 200 by both metrics. These remaining SNPs were again ranked by weighted RAS score, and the location of each of the top 50 SNPs was examined in detail (Supplemental table S2 in Additional file 1). Similar to other autoimmune diseases, SNPs which are associated with SLE in unbiased genome-wide mapping studies are largely within or near genes with immune system function , and we presumed that our study design which stratified patients by two immune system phenotypes would likely yield SNPs in genes related to the immune system. Therefore, we assumed that the prior probability of a true association was increased if 1) the SNP was located within or near a known gene, and 2) if that locus had immune system relevance based upon public databases. Nine SNPs of the top 50 met both of these criteria (Supplementary table S2 in Additional file 1).
SNP genotyping in the validation cohort
Individuals in the validation cohort were genotyped at the following SNPs using ABI Taqman Assays-by-Design primers and probes: rs10762360 (leucine-rich repeat containing 20, LRRC20), rs10784318 (protein phosphatase 1 H, PPM1H), rs10980684 (lysophosphatidic acid receptor 1, LPAR1), rs2820223 (ankyrin repeat and sterile alpha motif domain containing 1A, ANKS1A), rs17494870 (protein tyrosine phosphatase, receptor type M, PTPRM), rs26725 (ephrin A5, EFNA5), and rs112219769 (V-set and immunoglobulin domain containing 2, VSIG2). Scatter plots were all reviewed individually for quality, and none of the SNPs genotyped deviated significantly from the expected Hardy-Weinberg proportions (P > 0.01 for all SNPs in all ancestral backgrounds). Each SNP assay resulted in successful genotype calls in > 98% of the subjects, and none of the cases or controls were excluded due to incomplete or unsuccessful genotyping.
Statistical analysis in validation cohort
Logistic regression models were used to detect associations between the presence of autoantibodies and genotype at candidate SNPs. We performed this analysis because the initial GWAS study was stratified by both serum IFN-α and autoantibodies, and we sought to detect and clarify any potential associations with serology. Screening models incorporated each of the SNPs as predictor variables with a single autoantibody trait as the outcome variable in each ancestral background separately, and backward step-wise regression was performed discarding predictors which did not show evidence for association (P > 0.1). When more than one SNP was significantly associated with a particular antibody, interaction terms were included to detect potential SNP-SNP interactions. SLE-associated autoantibodies can show associations with other autoantibodies, (for example, the presence of anti-Sm is highly correlated with the presence of anti-RNP) and so we undertook additional logistic regression modeling to explore and refine SNP-antibody associations if one allele of a SNP was associated more strongly with a combination of antibodies than with either antibody alone. To test combination models, we generated new outcome variables which represented a combination of the two associated autoantibodies based upon the direction of the associations observed in the original regressions (such as presence of both anti-Ro and anti-dsDNA vs. the absence of one or both of those antibodies). If the antibody combination outcome variable provided significantly improved fit as compared to the individual autoantibody outcome variables, then the combination was reported as the best fit model.
To account for potential differences in admixture or population structure within self-reported ancestral backgrounds in the validation cohort, we performed a principal component analysis which included 13 independent SNPs (including the seven described in this study) that had been genotyped in all of our cases and controls. The first component obtained in this analysis provided a strong separation between subjects of self-reported African-American ancestry and self-reported European-American ancestry. Using the central value for this component to categorically split the population into African-derived vs. European-derived ancestry bins, self-reported African-American subjects were concordant with the African-derived ancestry bin 75.1% of the time, and self-reported European-Americans were concordant with the European-derived ancestry bin 86.3% of the time (P-value for a random distribution of self-reported ancestral backgrounds between these two bins = 1.3 × 10-91). Self-reported Hispanic-American ancestry subjects were largely overlapping with the European-American subjects in this analysis. We retested our autoantibody regression models and case-control analyses incorporating the first principal component as a quantitative covariate designating proportional European- vs. African-derived ancestry, and did not find significant differences in autoantibody associations with or without the inclusion of the first principal component in the regression (data not shown). The lack of a significant impact upon the reported associations provides some reassurance that these associations are not the result of differences between cases and controls in proportional admixture or population structure differences within self-reported ancestral backgrounds.
To analyze serum IFN-levels in the context of SNP genotype, column graphs were constructed to compare quantitative IFN-α levels between genotype categories. Because IFN-α levels were far from normally distributed, we used a non-parametric version of the two sample t-test (non-parametric Mann-Whitney U) to test if IFN-α levels varied significantly between two SNP genotype groups. Case-control allele frequency differences were tested using a chi-square test of independence. Reported P-values are not adjusted for multiple comparisons.
Results of GWAS screening study and multi-step screening algorithm
The top 50 SNPs from the screening GWAS study were determined using a combination of their Silhouette and weighted RAS scores (see Methods, Supplemental table S2 in Additional file 1), and we examined each of these SNPs in detail. When a SNP was located within or near a gene, a comprehensive literature and public database search was performed on each gene. We assumed that the prior probability of a true association was increased if, 1) the SNP was located within or near a known gene, and 2) if that locus had immune system relevance based upon public databases. Nine SNPs located in or near the following eight genes met both criteria: STAT4, LRRC20, PPM1 H, LPAR1, ANKS1A, PTPRM, EFNA5, and VSIG2. Two SNPs in the top 50 were located in the ANKS1A region, and in this case the SNP that was statistically stronger in the GWAS study was chosen for validation. The SNP in STAT4 (rs9967792) was located in the third intron of the gene. This SNP is in the LD block immediately adjacent to the SLE-associated haplotype block of STAT4 [24, 25]. We have previously demonstrated that variation in this third intron of STAT4 was associated with both serum IFN-α and autoantibody profile in multiple ancestral backgrounds in a study which included some of the same individuals , and we presumed that this was a recapitulation of that result, and did not pursue this SNP. Finding STAT4 in the top 50 SNPs demonstrated that true positive SNPs were indeed included among the top 50 in our screening GWAS study. Statistical results from the GWAS in each ancestral background separately for the seven SNPs chosen for validation are shown in Supplemental table S3 in Additional file 1.
Association of GWAS candidates with serologic subsets in SLE
Serologic associations observed in SLE patients in the replication cohort, best fit model from logistic regression
Locus and Allele
Associated Antibody Profile
Odds Ratio (95% CI)
LRRC20 rs10762360 G
(1.51 to 3.81)
1.6 × 10-4
PPM1 H rs10784318 C
EA and HA joint
(1.88 to 7.47)
3.6 × 10-5
LPAR1 rs10980684 T
Anti-Ro with Anti-Sm
AA and HA joint
(1.42 to 3.47)
1.2 × 10-4
ANKS1A rs2820223 C
Anti-Ro with Anti-dsDNA
(1.68 to 5.87)
3.5 × 10-4
PTPRM rs17494870 C
EFNA5 rs26725 C
(1.10 to 2.72)
VSIG2 rs11219769 T
Anti-RNP lacking Anti-Sm
EA and HA joint
(1.81 to 8.75)
2.1 × 10-4
Association of GWAS candidates with serum IFN-α activity in the SLE cohort
Case-control Association Study of the Seven Candidate loci from GWAS
There were no significant differences in allele frequencies between SLE cases and non-autoimmune controls in overall case-control analyses in African- and European-American ancestral backgrounds (Supplemental table S7 in Additional file 1). Differences in allele frequencies for each SNP in each ancestral background were modest (most ORs were in the 0.8 to 1.2 range). The size of our case-control populations do not allow us to exclude the possibility that an association exists, but it appears that the most dramatic genetic effects observed with these SNPs in SLE will be upon subphenotypes.
In this study, we used a detailed subphenotyping strategy to detect genetic factors associated with two important heritable phenotypes in SLE: autoantibody traits and serum IFN-α activity. Both of these traits have been associated with increased disease activity in SLE [5, 26], and both are thought to play a pathogenic role in disease initiation [9, 27]. The number and strength of the genetic associations discovered suggests that this technique could be a powerful tool in other diseases. A version of this study design has been used in rheumatoid arthritis, as patients are often stratified by the presence or absence of anti-cyclic citrullinated peptide antibodies, and large differences in genetic associations between these two serologic subsets have been observed . The seven loci that we found to be strongly associated with subphenotypes in SLE could potentially be discovered in a large case-control study, however these loci would likely be characterized by modest odds ratios for association in a traditional case-control study design based upon our case-control data. A limitation of our study is that we did not have complete clinical data available for all of the cases in the study, and thus we did not attempt to detect associations with clinical features of disease. Another limitation is that we did not have dense ancestry-informative genetic markers available in our cohort, and the pooled design of the GWAS study limited to some degree the amount of ancestry-related genetic information obtained from individual subjects in the GWAS portion of the study. As noted in the Methods, we tested associations within each self-reported ancestral background separately, used statistics that are robust to differences related to ancestry in the GWAS meta-analysis, and also used statistical methods to address potential confounding by differences in proportional ancestry within self-reported ancestral groups in the validation cohort. These analyses did not suggest that our results were confounded by ancestry. Follow-up replication and fine-mapping studies of these seven loci in independent multi-ancestral cohorts with dense ancestry-informative markers would be warranted.
The SNPs examined in this study are located in or near genes which are both novel in their association with SLE and highly interesting as potential immunomodulatory loci. The LRRC20 gene product has not been studied to date, but is known to contain leucine-rich repeats which are similar to the Toll-like receptor family of proteins involved in microbial pattern sensing and stimulation of cytokine pathways, including the IFN-α pathway. PTPRM is a receptor-type protein tyrosine phosphatase , which is of particular interest as a number of protein tyrosine phosphatases have been linked to autoimmune disease, including PTPN22 and PTPN2 [1, 30]. PTPRM mediates cell-cell aggregation via its extracellular domain , and the phosphatase domain allows intra-cellular signaling. PPM1 H is an intra-cellular phosphatase which is involved in regulation of apoptosis, possibly via the p53 pathway . LPAR1 and EFNA5 encode receptors expressed on leukocytes which can mediate apoptosis and chemotaxis , and leukocyte adhesion  respectively. Interestingly, genetic variation 124 kb upstream of LPAR1 was recently associated with circulating monocyte count in European populations . ANKS1A is also expressed in leukocytes and can be phosphorlyated by Lck , suggesting a potential role in T-cell signaling. After completion of our study, a SNP near ANKS1A in the immediately adjacent UHRF1 binding protein 1 (UHRF1BP1) locus was reported to confer modest risk of SLE (OR = 1.17) in a large case-control study in European ancestry . While we observed the major association between ANKS1A and SLE phenotype in African-Americans, it is possible that with a larger European ancestry cohort we may have observed an association between this locus and SLE phenotype, and this genomic region clearly warrants further study in multiple ancestral backgrounds. VSIG2, also known as cortical thymocyte receptor like, is expressed in the thymus and may be involved in antigen presentation . There are no known direct physical interactions between these gene products, although functionally each of these genes converge upon important immune system features involved in SLE pathogenesis, such as cell-cell interaction, antigen responses, and apoptosis.
Of the SLE risk loci which have previously been associated with serum IFN-α levels in SLE patients [15–19, 38], only STAT4 was represented in our top 50 SNP list. It is important to note that our screening GWAS study was not powered to exclude potential associations with serum IFN-α in SLE with any degree of certainty, and absence of a gene from our list does not mean that it cannot be associated with serum IFN-α phenotype in SLE. We were encouraged that one of these previously reported IFN-α-associated loci was present in the top 50 list, supporting the idea that we were able to identify some true positives. We would estimate that many more loci will be associated with the serum IFN-α phenotype in SLE, and this study is by no means exhaustive. It is interesting that the human leukocyte antigen (HLA) locus was not represented in our list of top SNPs, as this locus is known to be one of the strongest associations with SLE and many other autoimmune diseases [23, 39]. Different HLA alleles are known to be associated with particular autoantibody profiles in SLE . The HLA locus is highly complex and shows a large degree of variation by ancestry, and it is possible that the multi-ancestral nature of our screening GWAS cohort did not allow detection of HLA associations.
The success of this study design suggests that application of a similar technique in other diseases could be highly advantageous, allowing for robust gene discovery related to important pathways and subphenotypes using a smaller cohort than is typically studied in case-control GWAS. Additionally, we demonstrate a number of different patterns of association between SNP genotype and molecular subphenotypes (Figure 3). These intricate subphenotype relationships begin to suggest pathogenic models, which can be followed up in mechanistic studies. This is another strength of this study design, as more information regarding the potential function and role in disease pathogenesis is immediately available for the candidate genes as compared to standard case-control designs. The models outlined in Figure 3 can be tested functionally, allowing for integration of these loci into the mechanisms of human SLE pathogenesis.
We compare patients from the extremes of high vs. low interferon and autoantibody levels, and identify a number of novel genetic regions that are associated with high interferon alpha and particular autoantibodies in SLE patients. These genetic associations were very strong in the particular patient groups who demonstrated the specific findings and were not as strong in the overall cohort, suggesting that this approach was necessary to discover these genes. These genotype-phenotype relationships begin to suggest pathogenic models, which can be followed up in mechanistic studies and allow these genes to be integrated into the framework of human SLE pathogenesis. The number and strength of the genetic associations discovered suggests that this technique could be a powerful tool in other diseases. Because both autoantibody profile and serum IFN-α are important heritable pathogenic factors in SLE, we hypothesize that the gene variants we describe will play a role in disease susceptibility and severity in SLE.
- ANKS1A :
ankyrin repeat and sterile alpha motif domain containing 1A
anti-RNA-binding protein antibodies
- EFNA5 :
genome-wide association study
human leukocyte antigen
- IRF5 :
interferon regulatory factor 5
- IRF7 :
interferon regulatory factor 7
- LPAR1 :
lysophosphatidic acid receptor 1
- LRRC20 :
leucine-rich repeat containing 20
- PPM1H :
protein phosphatase 1H
- PTPN22 :
protein tyrosine phosphatase, non-receptor type 22
- PTPRM :
protein tyrosine phosphatase, receptor type, M
relative allele signals
systemic lupus erythematosus
single nucleotide polymorphism
- SPP1 :
- STAT4 :
signal transducer and activator of transcription 4
Translational Research in the Department of Medicine
- UHRF1BP1 :
UHRF1 binding protein 1
- VSIG2 :
V-set and immunoglobulin domain containing 2.
This work was supported by grants from the National Institutes of Health (AI083790, AI071651, and UL1 RR024999 to TBN; AI059893 to MKC; and NIDDK STEP UP Program to JA); the Lupus Research Institute (TBN and MKC); the Lupus Clinical Trials Consortium (TOU); the Alliance for Lupus Research (TBN and MKC), the Mary Kirkland Center for Lupus Research (MKC), and the Arthritis National Research Foundation Eng Tan Scholar Award (TBN). These funding bodies did not have any influence on the study design; collection, analysis, or interpretation of data, or in the writing of the manuscript, or in the decision to submit the manuscript for publication.
- Harley JB, Kelly JA, Kaufman KM: Unraveling the genetics of systemic lupus erythematosus. Springer Semin Immunopathol. 2006, 28: 119-130. 10.1007/s00281-006-0040-5.View ArticlePubMedGoogle Scholar
- Lopez P, Mozo L, Gutierrez C, Suarez A: Epidemiology of systemic lupus erythematosus in a northern Spanish population: gender and age influence on immunological features. Lupus. 2003, 12: 860-865. 10.1191/0961203303lu469xx.View ArticlePubMedGoogle Scholar
- Blanco P, Palucka AK, Gill M, Pascual V, Banchereau J: Induction of dendritic cell differentiation by IFN-alpha in systemic lupus erythematosus. Science. 2001, 294: 1540-1543. 10.1126/science.1064890.View ArticlePubMedGoogle Scholar
- Hooks JJ, Moutsopoulos HM, Geis SA, Stahl NI, Decker JL, Notkins AL: Immune interferon in the circulation of patients with autoimmune disease. N Engl J Med. 1979, 301: 5-8. 10.1056/NEJM197907053010102.View ArticlePubMedGoogle Scholar
- Kirou KA, Lee C, George S, Louca K, Peterson MG, Crow MK: Activation of the interferon-alpha pathway identifies a subgroup of systemic lupus erythematosus patients with distinct serologic features and active disease. Arthritis Rheum. 2005, 52: 1491-1503. 10.1002/art.21031.View ArticlePubMedGoogle Scholar
- Ronnblom LE, Alm GV, Oberg KE: Possible induction of systemic lupus erythematosus by interferon-alpha treatment in a patient with a malignant carcinoid tumour. J Intern Med. 1990, 227: 207-210. 10.1111/j.1365-2796.1990.tb00144.x.View ArticlePubMedGoogle Scholar
- Niewold TB, Swedler WI: Systemic lupus erythematosus arising during interferon-alpha therapy for cryoglobulinemic vasculitis associated with hepatitis C. Clin Rheumatol. 2005, 24: 178-181. 10.1007/s10067-004-1024-2.View ArticlePubMedGoogle Scholar
- Niewold TB: Interferon alpha-induced lupus: proof of principle. J Clin Rheumatol. 2008, 14: 131-132. 10.1097/RHU.0b013e318177627d.PubMed CentralView ArticlePubMedGoogle Scholar
- Niewold TB, Hua J, Lehman TJ, Harley JB, Crow MK: High serum IFN-alpha activity is a heritable risk factor for systemic lupus erythematosus. Genes Immun. 2007, 8: 492-502. 10.1038/sj.gene.6364408.PubMed CentralView ArticlePubMedGoogle Scholar
- Niewold TB, Adler JE, Glenn SB, Lehman TJ, Harley JB, Crow MK: Age- and sex-related patterns of serum interferon-alpha activity in lupus families. Arthritis Rheum. 2008, 58: 2113-2119. 10.1002/art.23619.PubMed CentralView ArticlePubMedGoogle Scholar
- Kariuki SN, Niewold TB: Genetic regulation of serum cytokines in systemic lupus erythematosus. Transl Res. 2010, 155: 109-117. 10.1016/j.trsl.2009.08.012.PubMed CentralView ArticlePubMedGoogle Scholar
- Ramos PS, Kelly JA, Gray-McGuire C, Bruner GR, Leiran AN, Meyer CM, Namjou B, Espe KJ, Ortmann WA, Reichlin M, Langefeld CD, James JA, Gaffney PM, Behrens TW, Harley JB, Moser KL: Familial aggregation and linkage analysis of autoantibody traits in pedigrees multiplex for systemic lupus erythematosus. Genes Immun. 2006, 7: 417-432. 10.1038/sj.gene.6364316.View ArticlePubMedGoogle Scholar
- Lovgren T, Eloranta ML, Bave U, Alm GV, Ronnblom L: Induction of interferon-alpha production in plasmacytoid dendritic cells by immune complexes containing nucleic acid released by necrotic or late apoptotic cells and lupus IgG. Arthritis Rheum. 2004, 50: 1861-1872. 10.1002/art.20254.View ArticlePubMedGoogle Scholar
- Niewold TB, Rivera TL, Buyon JP, Crow MK: Serum type I interferon activity is dependent on maternal diagnosis in anti-SSA/Ro-positive mothers of children with neonatal lupus. Arthritis Rheum. 2008, 58: 541-546. 10.1002/art.23191.PubMed CentralView ArticlePubMedGoogle Scholar
- Niewold TB, Kelly JA, Flesch MH, Espinoza LR, Harley JB, Crow MK: Association of the IRF5 risk haplotype with high serum interferon-alpha activity in systemic lupus erythematosus patients. Arthritis Rheum. 2008, 58: 2481-2487. 10.1002/art.23613.PubMed CentralView ArticlePubMedGoogle Scholar
- Salloum R, Franek BS, Kariuki SN, Rhee L, Mikolaitis RA, Jolly M, Utset TO, Niewold TB: Genetic variation at the IRF7/PHRF1 locus is associated with autoantibody profile and serum interferon-alpha activity in lupus patients. Arthritis Rheum. 2010, 62: 553-561. 10.1002/art.27182.PubMed CentralView ArticlePubMedGoogle Scholar
- Kariuki SN, Crow MK, Niewold TB: The PTPN22 C1858T polymorphism is associated with skewing of cytokine profiles toward high interferon-alpha activity and low tumor necrosis factor alpha levels in patients with lupus. Arthritis Rheum. 2008, 58: 2818-2823. 10.1002/art.23728.PubMed CentralView ArticlePubMedGoogle Scholar
- Kariuki SN, Moore JG, Kirou KA, Crow MK, Utset TO, Niewold TB: Age- and gender-specific modulation of serum osteopontin and interferon-alpha by osteopontin genotype in systemic lupus erythematosus. Genes Immun. 2009, 10: 487-494. 10.1038/gene.2009.15.PubMed CentralView ArticlePubMedGoogle Scholar
- Kariuki SN, Kirou KA, MacDermott EJ, Barillas-Arias L, Crow MK, Niewold TB: Cutting edge: Autoimmune disease risk variant of STAT4 confers increased sensitivity to IFN-alpha in lupus patients in vivo. J Immunol. 2009, 182: 34-38.PubMed CentralView ArticlePubMedGoogle Scholar
- Tan EM, Cohen AS, Fries JF, Masi AT, McShane DJ, Rothfield NF, Schaller JG, Talal N, Winchester RJ: The 1982 revised criteria for the classification of systemic lupus erythematosus. Arthritis Rheum. 1982, 25: 1271-1277. 10.1002/art.1780251101.View ArticlePubMedGoogle Scholar
- Hua J, Kirou K, Lee C, Crow MK: Functional assay of type I interferon in systemic lupus erythematosus plasma and association with anti-RNA binding protein autoantibodies. Arthritis Rheum. 2006, 54: 1906-1916. 10.1002/art.21890.View ArticlePubMedGoogle Scholar
- Pearson JV, Huentelman MJ, Halperin RF, Tembe WD, Melquist S, Homer N, Brun M, Szelinger S, Coon KD, Zismann VL, Webster JA, Beach T, Sando SB, Aasly JO, Heun R, Jessen F, Kolsch H, Tsolaki M, Daniilidou M, Reiman EM, Papassotiropoulos A, Hutton ML, Stephan DA, Craig DW: Identification of the genetic basis for complex disorders by use of pooling-based genomewide single-nucleotide-polymorphism association studies. Am J Hum Genet. 2007, 80: 126-139. 10.1086/510686.PubMed CentralView ArticlePubMedGoogle Scholar
- Harley JB, Alarcon-Riquelme ME, Criswell LA, Jacob CO, Kimberly RP, Moser KL, Tsao BP, Vyse TJ, Langefeld CD: 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. 10.1038/ng.81.PubMed CentralView ArticlePubMedGoogle Scholar
- Sigurdsson S, Nordmark G, Garnier S, Grundberg E, Kwan T, Nilsson O, Eloranta ML, Gunnarsson I, Svenungsson E, Sturfelt G, Bengtsson AA, Jonsen A, Truedsson L, Rantapaa-Dahlqvist S, Eriksson C, Alm G, Goring HH, Pastinen T, Syvanen AC, Ronnblom L: A risk haplotype of STAT4 for systemic lupus erythematosus is over-expressed, correlates with anti-dsDNA and shows additive effects with two risk alleles of IRF5. Hum Mol Genet. 2008, 17: 2868-2876. 10.1093/hmg/ddn184.PubMed CentralView ArticlePubMedGoogle Scholar
- Remmers EF, Plenge RM, Lee AT, Graham RR, Hom G, Behrens TW, de Bakker PI, Le JM, Lee HS, Batliwalla F, Li W, Masters SL, Booty MG, Carulli JP, Padyukov L, Alfredsson L, Klareskog L, Chen WV, Amos CI, Criswell LA, Seldin MF, Kastner DL, Gregersen PK: STAT4 and the risk of rheumatoid arthritis and systemic lupus erythematosus. N Engl J Med. 2007, 357: 977-986. 10.1056/NEJMoa073003.PubMed CentralView ArticlePubMedGoogle Scholar
- Rahman A, Hiepe F: Anti-DNA antibodies--overview of assays and clinical correlations. Lupus. 2002, 11: 770-773. 10.1191/0961203302lu313rp.View ArticlePubMedGoogle Scholar
- Arbuckle MR, McClain MT, Rubertone MV, Scofield RH, Dennis GJ, James JA, Harley JB: Development of autoantibodies before the clinical onset of systemic lupus erythematosus. N Engl J Med. 2003, 349: 1526-1533. 10.1056/NEJMoa021933.View ArticlePubMedGoogle Scholar
- Kallberg H, Padyukov L, Plenge RM, Ronnelid J, Gregersen PK, van der Helm-van Mil AH, Toes RE, Huizinga TW, Klareskog L, Alfredsson L: Gene-gene and gene-environment interactions involving HLA-DRB1, PTPN22, and smoking in two subsets of rheumatoid arthritis. Am J Hum Genet. 2007, 80: 867-875. 10.1086/516736.PubMed CentralView ArticlePubMedGoogle Scholar
- Aricescu AR, Siebold C, Choudhuri K, Chang VT, Lu W, Davis SJ, van der Merwe PA, Jones EY: Structure of a tyrosine phosphatase adhesive interaction reveals a spacer-clamp mechanism. Science. 2007, 317: 1217-1220. 10.1126/science.1144646.View ArticlePubMedGoogle Scholar
- Parkes M, Barrett JC, Prescott NJ, Tremelling M, Anderson CA, Fisher SA, Roberts RG, Nimmo ER, Cummings FR, Soars D, Drummond H, Lees CW, Khawaja SA, Bagnall R, Burke DA, Todhunter CE, Ahmad T, Onnie CM, McArdle W, Strachan D, Bethel G, Bryan C, Lewis CM, Deloukas P, Forbes A, Sanderson J, Jewell DP, Satsangi J, Mansfield JC, Cardon L, et al: Sequence variants in the autophagy gene IRGM and multiple other replicating loci contribute to Crohn's disease susceptibility. Nat Genet. 2007, 39: 830-832. 10.1038/ng2061.PubMed CentralView ArticlePubMedGoogle Scholar
- Sugiura T, Noguchi Y, Sakurai K, Hattori C: Protein phosphatase 1 H, overexpressed in colon adenocarcinoma, is associated with CSE1L. Cancer Biol Ther. 2008, 7: 285-292.View ArticlePubMedGoogle Scholar
- Hu X, Haney N, Kropp D, Kabore AF, Johnston JB, Gibson SB: Lysophosphatidic acid (LPA) protects primary chronic lymphocytic leukemia cells from apoptosis through LPA receptor activation of the anti-apoptotic protein AKT/PKB. J Biol Chem. 2005, 280: 9498-9508. 10.1074/jbc.M410455200.View ArticlePubMedGoogle Scholar
- Wimmer-Kleikamp SH, Nievergall E, Gegenbauer K, Adikari S, Mansour M, Yeadon T, Boyd AW, Patani NR, Lackmann M: Elevated protein tyrosine phosphatase activity provokes Eph/ephrin-facilitated adhesion of pre-B leukemia cells. Blood. 2008, 112: 721-732. 10.1182/blood-2007-11-121681.View ArticlePubMedGoogle Scholar
- Ferreira MA, Hottenga JJ, Warrington NM, Medland SE, Willemsen G, Lawrence RW, Gordon S, de Geus EJ, Henders AK, Smit JH, Campbell MJ, Wallace L, Evans DM, Wright MJ, Nyholt DR, James AL, Beilby JP, Penninx BW, Palmer LJ, Frazer IH, Montgomery GW, Martin NG, Boomsma DI: Sequence variants in three Loci influence monocyte counts and erythrocyte volume. Am J Hum Genet. 2009, 85: 745-749. 10.1016/j.ajhg.2009.10.005.PubMed CentralView ArticlePubMedGoogle Scholar
- Emaduddin M, Edelmann MJ, Kessler BM, Feller SM: Odin (ANKS1A) is a Src family kinase target in colorectal cancer cells. Cell Commun Signal. 2008, 6: 7-10.1186/1478-811X-6-7.PubMed CentralView ArticlePubMedGoogle Scholar
- Gateva V, Sandling JK, Hom G, Taylor KE, Chung SA, Sun X, Ortmann W, Kosoy R, Ferreira RC, Nordmark G, Gunnarsson I, Svenungsson E, Padyukov L, Sturfelt G, Jonsen A, Bengtsson AA, Rantapaa-Dahlqvist S, Baechler EC, Brown EE, Alarcon GS, Edberg JC, Ramsey-Goldman R, McGwin G, Reveille JD, Vila LM, Kimberly RP, Manzi S, Petri MA, Lee A, Gregersen PK, et al: A large-scale replication study identifies TNIP1, PRDM1, JAZF1, UHRF1BP1 and IL10 as risk loci for systemic lupus erythematosus. Nat Genet. 2009, 41: 1228-1233. 10.1038/ng.468.PubMed CentralView ArticlePubMedGoogle Scholar
- Chretien I, Marcuz A, Courtet M, Katevuo K, Vainio O, Heath JK, White SJ, Du Pasquier L: CTX, a Xenopus thymocyte receptor, defines a molecular family conserved throughout vertebrates. Eur J Immunol. 1998, 28: 4094-4104. 10.1002/(SICI)1521-4141(199812)28:12<4094::AID-IMMU4094>3.0.CO;2-2.View ArticlePubMedGoogle Scholar
- Harley IT, Niewold TB, Stormont RM, Kaufman KM, Glenn SB, Franek BS, Kelly JA, Kilpatrick JR, Hutchings D, Divers J, Bruner GR, Edberg JC, McGwin G, Petri MA, Ramsey-Goldman R, Reveille JD, Vila-Pérez LM, Merrill JT, Gilkeson GS, Vyse TJ, Alarcon-Riquelme ME, Cho SK, Jacob CO, Alarcon GS, Moser KL, Gaffney PM, Kimberly RP, Bae SC, Langefeld CD, Harley JB, et al: The role of genetic variation near interferon-kappa in systemic lupus erythematosus. J Biomed Biotechnol. 2010, 2010: Article ID 706825-View ArticleGoogle Scholar
- Plenge RM, Seielstad M, Padyukov L, Lee AT, Remmers EF, Ding B, Liew A, Khalili H, Chandrasekaran A, Davies LR, Li W, Tan AK, Bonnard C, Ong RT, Thalamuthu A, Pettersson S, Liu C, Tian C, Chen WV, Carulli JP, Beckman EM, Altshuler D, Alfredsson L, Criswell LA, Amos CI, Seldin MF, Kastner DL, Klareskog L, Gregersen PK: TRAF1-C5 as a risk locus for rheumatoid arthritis--a genomewide study. N Engl J Med. 2007, 357: 1199-1209. 10.1056/NEJMoa073491.PubMed CentralView ArticlePubMedGoogle Scholar
- Graham RR, Ortmann W, Rodine P, Espe K, Langefeld C, Lange E, Williams A, Beck S, Kyogoku C, Moser K, Gaffney P, Gregersen PK, Criswell LA, Harley JB, Behrens TW: Specific combinations of HLA-DR2 and DR3 class II haplotypes contribute graded risk for disease susceptibility and autoantibodies in human SLE. Eur J Hum Genet. 2007, 15: 823-830. 10.1038/sj.ejhg.5201827.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.