Figures
Abstract
Background
Systemic lupus erythematosus (SLE) is a systemic multisystem autoimmune disorder influenced by genetic background and environmental factors. Our aim here was to replicate findings of associations between 7 of the implicated single nucleotide polymorphisms (SNPs) in IRF5, BLK, STAT4, TNFAIP3, SPP1, TNIP1 and ETS1 genes with susceptibility to childhood-onset SLE in the Japanese population. In particular, we focused on gender differences in allelic frequencies.
Methodology/Principal Findings
The 7 SNPs were genotyped using TaqMan assays in 75 patients with childhood-onset SLE and in 190 healthy controls. The relationship between the cumulative number of risk alleles and SLE manifestations was explored in childhood-onset SLE. Logistic regression was used to test the effect of each polymorphism on susceptibility to SLE, and Wilcoxon rank sum testing was used for comparison of total risk alleles. Data on rs7574865 in the STAT4 gene and rs9138 in SPP1 were replicated for associations with SLE when comparing cases and controls (corrected P values ranging from 0.0043 to 0.027). The rs2230926 allele of TNFAIP3 was associated with susceptibility to SLE in males, but after Bonferroni correction there were no significant associations with any of the other four SNPs in IRF5, BLK, TNIP1 and ETS1 genes. The cumulative number of risk alleles was significantly increased in childhood-onset SLE relative to healthy controls (P = 0.0000041). Male SLE patients had a slightly but significantly higher frequency of the TNFAIP3 (rs2230926G) risk allele than female patients (odds ratio [OR] = 4.05, 95% confidence interval [95%CI] = 1.46–11.2 P<0.05).
Conclusions
Associations of polymorphisms in STAT4 and SPP1 with childhood-onset SLE were confirmed in a Japanese population. Although these are preliminary results for a limited number of cases, TNFAIP3 rs2230926G may be an important predictor of disease onset in males. We also replicated findings that the cumulative number of risk alleles was significantly increased in childhood-onset SLE.
Citation: Kadota K, Mori M, Yanagimachi M, Miyamae T, Hara T, Kanetaka T, et al. (2013) Analysis of Gender Differences in Genetic Risk: Association of TNFAIP3 Polymorphism with Male Childhood-Onset Systemic Lupus Erythematosus in the Japanese Population. PLoS ONE 8(8): e72551. https://doi.org/10.1371/journal.pone.0072551
Editor: David L. Boone, University of Chicago, United States of America
Received: March 24, 2013; Accepted: July 11, 2013; Published: August 30, 2013
Copyright: © 2013 Kadota et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This work was supported by a Grant-in-Aid for Scientific Research from the Japan Society for the Promotion of Science (No. 16790583). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Systemic lupus erythematosus (SLE) is an autoimmune connective tissue disease which presents with various clinical symptoms and disorders in many internal organs. Disease onset in about 60% of patients occurs between the ages of 16 and 55 years, and approximately 15–20% are diagnosed in adolescence [1]. SLE occurs predominantly in women, with a female-to-male ratio of 10∶1. However, for childhood-onset SLE, this ratio is lower, at about 3∶1.
Childhood-onset SLE represents 10–20% of all SLE cases, and is associated with greater disease severity than adult-onset SLE, including more rapid renal damage. The majority of childhood-onset SLE patients will have developed renal damage within 5–10 years of disease onset [2]. Other studies of clinical characteristics indicated that men younger than 33 at the time of diagnosis were more likely to develop lupus nephritis, and that patients of Asian descent tend to have more severe disease, including the renal manifestations [3]. Recent large-scale association studies revealed a relationship between genetic risk factors and patient age at disease onset [4]. Childhood-onset SLE is associated with more frequent occurrence of the various SLE symptoms than seen in adult-onset disease, and has a more serious disease course [1]. Relative to adult-onset SLE, there is a higher genetic risk in childhood-onset SLE; this influences age of disease onset, which is an important predictor of disease severity [4]. Therefore, we hypothesize that gender-specific genetic effects of childhood-onset SLE are involved in its etiology. To test this hypothesis, here we explore the distribution of known susceptibility loci in a cohort of Japanese childhood-onset SLE patients.
Recent studies have identified the following genes as affecting susceptibility to SLE in Japanese: human leukocyte antigen DRB1 (HLA-DRB1) *15∶01, interferon regulatory factor 5 (IRF5) rs41298401C, signal transducer and activator of transcription-4 (STAT4) rs7574865T, B lymphoid tyrosine kinase (BLK) rs13277113A, tumor necrosis factor alpha-induced protein 3 (TNFAIP3) rs2230926G and TNFAIP3-interacting protein 1 (TNIP1) rs7708392C [5]–[10]. The same risk alleles are associated with SLE in the Caucasian and Chinese populations [11]–[14]. We previously evaluated associations of polymorphisms in IRF5, BLK, STAT4, TNFAIP3 and TNIP1 genes with susceptibility to SLE in Japanese children [7], [9].
Osteopontin (SPP1: secreted phosphoprotein 1) located on 4q22, also called early T lymphocyte activation 1, is a bone matrix mediator with important roles in inflammation and immunity [15]. An association between SPP1 genetic polymorphisms and increased osteopontin protein levels has been reported [16]. Two SPP1 single nucleotide polymorphisms (SNPs) (rs1126772 and rs9138) in the 3′ untranslated region were found to be associated with autoimmune/lymphoproliferative syndrome (ALPS). Osteopontin serum levels are associated both with these polymorphisms and an elevated risk of developing ALPS [17]. The same investigators reported that rs7687316 and rs9138 were significantly associated with SLE in an Italian population [18]. Another report implicated rs1126616 and rs9138 in conferring a significantly higher risk of SLE on men than women [19]. Therefore, it was suggested that rs9138 in the SPP1 gene is a specific autosomal gene in male SLE.
E26 transformation-specific 1 (ETS1) is a member of the E26 transformation-specific family of transcription factors. ETS1 plays a role in the negative regulation of Th17 cells and in B cell differentiation [20]. Recently, SLE susceptibility genes have been identified in two genome-wide association studies (GWAS) in Asian populations. It was found that genetic variants in ETS1 were associated with susceptibility to SLE; reduced expression of ETS1 may increase differentiation and activity of both plasma cells and Th17 cells [21], [22]. Two SNPs in ETS1 (rs6590330 and rs4937333) were found to be significantly related to SLE, and reported to be associated with age of onset [11].
Recently, an association of the cumulative number of risk alleles and SLE has been reported in the United States and Japan [4], [23]. Koga et al. analyzed the cumulative number of risk alleles at eight established susceptibility loci, namely HLA-DRB1, IRF5, STAT4, BLK, TNFAIP3, TNIP1, FCGR2B and TNFSF13, in a Japanese population. Patients were found to carry more risk alleles than healthy controls. The odds ratio (OR) for SLE was significantly increased in individuals carrying 10 or more risk alleles. In contrast, subjects with 4 or less risk alleles had significantly decreased SLE occurrence.
The present study was conducted to determine whether the 7 candidate genes IRF5, BLK, STAT4, TNFAIP3, SPP1, TNIP1 and ETS1 were associated with childhood-onset SLE in the Japanese population. We investigated sex-specific differences in allelic frequencies and the cumulative number of risk alleles. We further examined the relationship between genetic risk for SLE and many clinical manifestations of the disease.
Materials and Methods
Ethics Statement
This study was performed in accordance with the Declaration of Helsinki and approved by the Ethics Committee of Yokohama City University School of Medicine. Written informed consent was given by parents/guardians for all children that participated in this study.
Subjects
The aim of this study was to explore associations between SNPs in the seven genes IRF5, BLK, STAT4, TNFAIP3, SPP1, TNIP1 and ETS1 and susceptibility to childhood-onset SLE in a Japanese population. The distribution of SNPs in these genes confirmed that they acted as SLE susceptibility loci (Table S1) [6]–[11], [19]. We analyzed the genotype data at the seven risk SNPs and investigated associations between cumulative number of risk alleles and susceptibility to SLE. In addition, we examined gender differences in allelic frequencies and cumulative risk for SLE. Moreover, we investigated associations between these SNPs and the different clinical characteristics of the patients.
Study Population and Clinical Data
A total of 75 children, 15 males (20.0%) and 60 females (80.0%), enrolled in this study were followed-up at Yokohama City University Hospital between December 2007 and December 2010. The mean age of patients was about 12 years at onset of SLE (Table 1). All patients fulfilled the American College of Rheumatology (ACR) Criteria score for SLE of 4 or higher [24], [25]. We explored gender differences in the clinical manifestations (malar rash, arthritis, serositis, thrombocytopenia, proteinuria and anti-dsDNA antibodies, according to the ACR classification criteria).
Genotyping
SNPs in the BLK (rs13277113), STAT4 (rs7574865), TNFAIP3 (rs2230926), SPP1 (rs9138), TNIP1 (rs7708392) and ETS1 (rs4937333) genes were selected based on previous research [7]–[10], [21], [26]. Seventy five patients with SLE and 190 healthy controls were genotyped. Genomic DNA was isolated from peripheral blood using the QIAamp DNA Mini kit (Qiagen K.K, Tokyo, Japan). Genotyping was performed using TaqMan SNP Genotyping Assays (Applied Biosystems, Foster, CA, USA assay ID: C___1886916_10 for rs13277113, C__29882391_10 for rs7574865, C___7701116_10 for rs2230926, C___8826997_10 for rs9138, C__29349759_10 for rs7708392, C____278298_10 for rs4937333). We developed a custom-made assay for rs41298401 in IRF5 (forward primer, 5′-GCTGCGCCTGGAAAGC-3′; reverse primer, 5′-GGGAGGCGCTTTGGAAGT-3′; probes, VIC-CTGCTGTAGGCACCC and FAM-CTGCTGTAGCCACCC; Custom TaqMan Applied Biosystems, Foster, CA, USA). To confirm the accuracy of genotyping, we conducted bidirectional sequencing on 6 randomly selected individuals for each genotype at rs41298401 and undertook direct sequencing of the SNP. No genotyping errors were found. These SNPs were analyzed by real-time PCR using the AB7500 Real Time PCR system (Applied Biosystems, Foster, CA, USA) under the conditions recommended by the manufacturer.
Statistical Analysis
Hardy-Weinberg Equilibrium (HWE) was measured by chi-squared testing for the control groups of each study (significance set at P<0.05). All analyses were performed using IBM SPSS software (SPSS for Windows, Version 19.0, Chicago, IL, USA). The multiple comparison P-values were corrected by the Bonferroni method, in which P values were multiplied by the number of tested SNPs, and P<0.05 considered significant (2-sided test).
We compared the clinical manifestations of SLE cases between male and female patients using Fisher’s exact probability test. Logistic regression was used to test the effect of each of the 7 SNPs on differences between SLE cases and healthy controls. A genetic additive model (with values: AA = 2, Aa = 1 and aa = 0, respectively) was considered. We analyzed the correlation between the risk for each SNP using a dominant or recessive model in logistic regression analysis. Dominant (risk allele homozygotes = 1, heterozygotes = 1 and nonrisk allele homozygotes = 0) and recessive (risk allele homozygotes = 1, heterozygotes = 0 and nonrisk allele homozygotes = 0) models were tested. Results are expressed as odds ratios (OR) with 95% confidence intervals (CI). Comparison of male and female patients used logistic regression [27], which was also used to test the effect of each polymorphism separately. Gender-specific associations were assessed by comparing the additive model between cases and controls for each gender separately.
The cumulative number of risk alleles was determined in each patient. Wilcoxon rank sum testing was used for comparisons of groups.
Results
Study Population and Clinical Data
We studied 75 cases and 190 controls. The median age at SLE diagnosis was 11.9 years (Table 1). There were 15 male and 60 female patients. All patients and controls (male n = 91, female n = 99) were Japanese. Female gender was a risk factor for SLE susceptibility (P<0.001, Fisher’s exact test). However, proteinuria was higher in males than females with SLE (P = 0.026, Fisher’s exact test) (Table 2), while malar rash, arthritis, serositis, thrombocytopenia and the presence of anti-dsDNA antibodies were not significantly difference between the sexes.
Associations between the 7 Tested SNPs and Susceptibility to SLE
We conducted an analysis of Hardy-Weinberg equilibrium and genotypic and allelic associations of 7 SNPs with susceptibility to SLE. The genotype frequencies for the SNPs in both patients with SLE and healthy controls were in Hardy-Weinberg equilibrium (P>0.05).
A significant difference in allele distribution of STAT4 rs7574865T between SLE patients and controls (OR = 1.77, 95%CT = 1.20–2.62, P = 0.027) was observed (Table 3). The genotype frequency of rs9138 in SPP1 was also significantly different between SLE patients and controls (OR = 3.31, 95%CI = 1.67–6.56, P = 0.0043) in the recessive model (AA vs AC/CC) (Table 3). The genotype AA was significantly higher in patients than in controls.
Clinical Subsets
Next, we explored whether these risk alleles were associated with the development of particular clinical manifestations. Frequencies of the rs9138 A-allele in SPP1 (OR = 0.48, 95%CI = 0.25–0.92, P = 0.028) and the rs4937333 T-allele in ETS1 (OR = 2.28, 95%CI = 1.03–5.04, P = 0.042) were significantly increased in SLE patients with proteinuria as compared to those without. There were no significant associations between TNFAIP3 rs2230926G and proteinuria in SLE patients. Malar rash was associated with STAT4 rs7574865T (OR = 3.08, 95%CI = 1.21–7.85, P = 0.019), whereas no significant associations were observed in the other subsets of patients with arthritis, serositis, thrombocytopenia and anti-dsDNA antibodies.
The Cumulative Number of Risk Alleles at the 7 SNPs and Gender Differences
We had selected these 7 established SLE susceptibility alleles due to their reported significance in Asian populations or because they conferred a higher risk of SLE to males than females (Table S1). The cumulative risk was calculated in each patient by summing the risk scores. It was found that the cumulative number of the risk alleles was significantly increased in patients with childhood-onset SLE relative to the healthy controls (P = 0.0000041, Wilcoxon rank sum test). In contrast, there were no significant associations with gender itself (P = 0.085, Wilcoxon rank sum test).
Secondly, we evaluated the effect of these SNPs separately on the gender difference in SLE susceptibility. A significant gender-gene interaction was observed for rs2230926, a polymorphism located in TNFAIP3. The frequency of the risk allele in TNFAIP3 (rs2230926G) was slightly but significantly higher in men than women with SLE (OR = 4.05, 95%CI = 1.46–11.23, P<0.05) (Table 4), whereas there were no differences between men (n = 91) and women (n = 99) in the controls. To further assess the effect of TNFAIP3 (rs2230926), we calculated associations with SLE in gender-specific case-control analysis. Rs2230926 was associated with SLE in males (case n = 15, control n-91) but not in females (case n = 60, control n = 99) (OR male = 6.17, 95%CI = 2.10–18.09, P = 0.0064 corrected; OR female = 1.25, 95%CI = 0.52–3.02, P = 0.62 uncorrected) (Table 5).
Discussion
SLE is a systemic multisystem autoimmune disorder, susceptibility to which is influenced by a complex genetic background and potentially by environmental exposures. Because chronologically, childhood has less time to accumulate environmental exposures, the impact of genetic factors may be easier to disentangle than in adult-onset SLE. Recent studies identified many susceptibility genes relevant to SLE. We have replicated results for 3 of these (STAT4, TNFAIP3 and SPP1) in the present study, but we were unable to replicate the findings of susceptibility for four others (IRF5, BLK, ETS1 and TNIP1). The frequencies of the risk alleles IRF5, BLK, ETS1 and TNIP1 in the present study are similar to those published in other recent Asian studies [6], [8], [10], [11]. We may not have seen a significant difference because of the small size of our study.
The cumulative number of risk alleles is reported to be significantly increased in SLE relative to healthy controls in the Japanese population [23]. We did replicate this association in Japanese childhood-onset SLE. A recent study of cumulative genetic risk scores for SLE in men and women concluded that there is a higher genetic risk factor in men [28] and that male patients had a higher frequency of clinical manifestations than women [29]. In the present study, we found no significant differences in the number of SLE risk alleles in male and female patients. Female predominance in SLE is clearly evident in general, but the male: female ratio is lower in children with prepubertal onset and onset recently after puberty than at reproductive ages [30]. We suggest that genetic susceptibility is more clearly associated with childhood-onset SLE than adult-onset SLE. We did not find a significant difference in the cumulative number of risk alleles between male and female SLE patients. Therefore, the increased frequency of adult-onset SLE among women may be attributed to differences in environmental factors and the metabolism of sex hormones in particular.
Our findings suggest that susceptibility risks conferred by polymorphisms in the three genes STAT4, TNFAIP3, SPP1 participate in the pathogenesis of childhood-onset SLE, but identified the TNFAIP3 rs2230926G allele as the only one significantly different in male and female patients. This finding suggests that TNFAIP3 rs2230926G may be an important predictor of disease in males. However, this should still be considered a tentative conclusion because there were only 15 male patients in our study, and previous large-scale studies of TNFAIP3 have not identified a similar gender effect. Therefore, a multicenter study of the same ethnic group is needed to validate the results from this study. A20 (also known as TNFAIP3) encodes a ubiquitin-editing protein that acts in the negative regulation of nuclear factor-kappaB (NF-κB) signaling via several pathways including those involving tumor necrosis factor (TNF) and Toll-like receptors (TLR) [31]. Resent genome-wide association studies (GWAS) have revealed associations between TNFAIP3 and SLE [32]. The SNPs rs2230926 in TNFAIP3 introduces the amino acid substitution Phe127 to Cys127. Previous studies demonstrated that TNFAIP3 risk allele rs2230926G leads to reduced NF-κB inhibition. Levels of TNFAIP3 mRNA are lower in carriers of the rs2230926 G-allele (GT subjects) than in non-carriers (TT), resulting in less protein expression and less A20 inhibitory activity [9], [33]. Further research is needed to determine the influence of gene polymorphisms in the TNFAIP3 gene on the potential role of A20 protein in the pathogenesis of SLE.
A recent review reported that a high a percentage of childhood-onset SLE patients have renal involvement [34]. In the present study, we observed that SPP1 rs9138C was associated with proteinuria in SLE patients (OR = 2.10, 95%CI = 1.09–4.08, P = 0.028). On the other hand, the frequency of SPP1 rs9138C tended to be decreased in the SLE patients studied here (OR = 0.62, 95%CI = 0.42–0.91, P = 0.15). Comparing the data from our Japanese population with those previously published from a Caucasian population, the frequency of the C allele at rs9138 was found to be significantly increased [26]. SPP1 rs9138C was reported to contribute to SLE susceptibility in a population including both African-American and European-American males, but not females [19]. It can be suggested that SPP1 rs9138C is only associated with increased risk of renal disorders in childhood- but not adult-onset SLE patients. Kariuki et al reported age- and gender- dependent effects of SPP1 (rs9138C) on cytokine profiles in SLE patients. The rs9138C allele is associated with higher serum osteopontin and interferon-alpha in male and female patients ≤23 years old. [35] These results are consistent with our findings.
Our study is the first to report an association between risk alleles and clinical manifestations in Japanese children with SLE. In particular, both geography and race affect the prevalence of SLE and the frequency and severity of clinical manifestations. Although our study is of a small single-center cohort, this minimizes any discrepancies in management and assessment, and it reflects the homogeneity of the ethnic and socio-demographic background. In addition, because it focuses exclusively on patients with childhood-onset SLE, it is more appropriate for defining detailed genetic susceptibility in this group than cohorts investigating all SLE patients. Our hospital is a pediatric rheumatology center treating many patients with severe SLE. However, we acknowledge that the current study suffers from some shortcomings. The sample size for male SLE patients is small (15 males). In addition, further research is needed to study relationships between adult-onset SLE and childhood-onset patients. Therefore, it is important that the association of TNFAIP3 and SPP1 polymorphisms with childhood-onset SLE susceptibility is confirmed by other groups.
Conclusions
Associations of STAT4 and SPP1 with childhood-onset SLE were confirmed in a Japanese population. For TNFAIP3, this genetic effect may be confined to males. We replicated previous findings that the cumulative number of risk alleles of seven known SNPs was significantly increased in childhood-onset SLE relative to controls, while there was no association with the observed gender difference.
Supporting Information
Table S1.
Single nucleotide polymorphisms (SNPs) representing the 7 systemic lupus erythematosus (SLE) susceptibility loci genotyped in our study.
https://doi.org/10.1371/journal.pone.0072551.s001
(XLS)
Author Contributions
Conceived and designed the experiments: KK MM MY SY MK. Performed the experiments: KK MM MY TM TH T. Kanetaka TN RH TI. Analyzed the data: KK MM MY T. Kaneko. Contributed reagents/materials/analysis tools: KK MM MY. Wrote the paper: KK MM MY T. Kaneko SY.
References
- 1. Brunner HI, Gladman DD, Ibanez D, Urowitz MD, Silverman ED (2008) Difference in disease features between childhood-onset and adult-onset systemic lupus erythematosus. Arthritis Rheum 58: 556–562.
- 2. Kamphuis S, Silverman ED (2010) Prevalence and burden of pediatric-onset systemic lupus erythematosus. Nat Rev Rheumatol 6: 538–546.
- 3. Seligman VA, Lum RF, Olson JL, Li H, Criswell LA (2002) Demographic differences in the development of lupus nephritis: a retrospective analysis. Am J Med 112: 726–729.
- 4. Webb R, Kelly JA, Somers EC, Hughes T, Kaufman KM, et al. (2011) Early disease onset is predicted by a higher genetic risk for lupus and is associated with a more severe phenotype in lupus patients. Ann Rheum Dis 70: 151–156.
- 5. Tsuchiya N, Ohashi J, Tokunaga K (2002) Variations in immune response genes and their associations with multifactorial immune disorders. Immunol Rev 190: 169–181.
- 6. Kawasaki A, Kyogoku C, Ohashi J, Miyashita R, Hikami K, et al. (2008) Association of IRF5 polymorphisms with systemic lupus erythematosus in a Japanese population: support for a crucial role of intron 1 polymorphisms. Arthritis Rheum 58: 826–834.
- 7. Kawasaki A, Ito I, Hikami K, Ohashi J, Hayashi T, et al. (2008) Role of STAT4 polymorphisms in systemic lupus erythematosus in a Japanese population: a case-control association study of the STAT1–STAT4 region. Arthritis Res Ther 10: R113.
- 8. Ito I, Kawasaki A, Ito S, Hayashi T, Goto D, et al. (2009) Replication of the association between the C8orf13-BLK region and systemic lupus erythematosus in a Japanese population. Arthritis Rheum 60: 553–558.
- 9. Kawasaki A, Ito I, Ito S, Hayashi T, Goto D, et al. (2010) Association of TNFAIP3 polymorphism with susceptibility to systemic lupus erythematosus in a Japanese population. J Biomed Biotechnol 2010: 207578.
- 10. Kawasaki A, Ito S, Furukawa H, Hayashi T, Goto D, et al. (2010) Association of TNFAIP3 interacting protein 1, TNIP1 with systemic lupus erythematosus in a Japanese population: a case-control association study. Arthritis Res Ther 12: R174.
- 11. Zhong H, Li XL, Li M, Hao LX, Chen RW, et al. (2011) Replicated associations of TNFAIP3, TNIP1 and ETS1 with systemic lupus erythematosus in a southwestern Chinese population. Arthritis Res Ther 13: R186.
- 12. Taylor KE, Remmers EF, Lee AT, Ortmann WA, Plenge RM, et al. (2008) Specificity of the STAT4 genetic association for severe disease manifestations of systemic lupus erythematosus. PLoS Genet 4: e1000084.
- 13. Hom G, Graham RR, Modrek B, Taylor KE, Ortmann W, et al. (2008) Association of systemic lupus erythematosus with C8orf13-BLK and ITGAM-ITGAX. N Engl J Med 358: 900–909.
- 14. Musone SL, Taylor KE, Lu TT, Nititham J, Ferreira RC, et al. (2008) Multiple polymorphisms in the TNFAIP3 region are independently associated with systemic lupus erythematosus. Nat Genet 40: 1062–1064.
- 15. Ashkar S, Weber GF, Panoutsakopoulou V, Sanchirico ME, Jansson M, et al. (2000) Eta-1 (osteopontin): an early component of type-1 (cell-mediated) immunity. Science 287: 860–864.
- 16. Wong CK, Lit LC, Tam LS, Li EK, Lam CW (2005) Elevation of plasma osteopontin concentration is correlated with disease activity in patients with systemic lupus erythematosus. Rheumatology (Oxford) 44: 602–606.
- 17. Chiocchetti A, Indelicato M, Bensi T, Mesturini R, Giordano M, et al. (2004) High levels of osteopontin associated with polymorphisms in its gene are a risk factor for development of autoimmunity/lymphoproliferation. Blood 103: 1376–1382.
- 18. D’Alfonso S, Barizzone N, Giordano M, Chiocchetti A, Magnani C, et al. (2005) Two single-nucleotide polymorphisms in the 5′ and 3′ ends of the osteopontin gene contribute to susceptibility to systemic lupus erythematosus. Arthritis Rheum 52: 539–547.
- 19. Han S, Guthridge JM, Harley IT, Sestak AL, Kim-Howard X, et al. (2008) Osteopontin and systemic lupus erythematosus association: a probable gene-gender interaction. PLoS One 3: e0001757.
- 20. Leng RX, Pan HF, Chen GM, Feng CC, Fan YG, et al. (2011) The dual nature of Ets-1: focus to the pathogenesis of systemic lupus erythematosus. Autoimmun Rev 10: 439–443.
- 21. Yang W, Shen N, Ye DQ, Liu Q, Zhang Y, et al. (2010) Genome-wide association study in Asian populations identifies variants in ETS1 and WDFY4 associated with systemic lupus erythematosus. PLoS Genet 6: e1000841.
- 22. Han JW, Zheng HF, Cui Y, Sun LD, Ye DQ, et al. (2009) Genome-wide association study in a Chinese Han population identifies nine new susceptibility loci for systemic lupus erythematosus. Nat Genet 41: 1234–1237.
- 23. Koga M, Kawasaki A, Ito I, Furuya T, Ohashi J, et al. (2011) Cumulative association of eight susceptibility genes with systemic lupus erythematosus in a Japanese female population. J Hum Genet 56: 503–507.
- 24. Tan EM, Cohen AS, Fries JF, Masi AT, McShane DJ, et al. (1982) The 1982 revised criteria for the classification of systemic lupus erythematosus. Arthritis Rheum 25: 1271–1277.
- 25. Hochberg MC (1997) Updating the American College of Rheumatology revised criteria for the classification of systemic lupus erythematosus. Arthritis Rheum 40: 1725.
- 26. Trivedi T, Franek BS, Green SL, Kariuki SN, Kumabe M, et al. (2011) Osteopontin alleles are associated with clinical characteristics in systemic lupus erythematosus. J Biomed Biotechnol 2011: 802581.
- 27. Gonzalez JR, Armengol L, Sole X, Guino E, Mercader JM, et al. (2007) SNPassoc: an R package to perform whole genome association studies. Bioinformatics 23: 644–645.
- 28. Hughes T, Adler A, Merrill JT, Kelly JA, Kaufman KM, et al. (2012) Analysis of autosomal genes reveals gene-sex interactions and higher total genetic risk in men with systemic lupus erythematosus. Ann Rheum Dis 71: 694–699.
- 29. Stefanidou S, Benos A, Galanopoulou V, Chatziyannis I, Kanakoudi F, et al. (2011) Clinical expression and morbidity of systemic lupus erythematosus during a post-diagnostic 5-year follow-up: a male:female comparison. Lupus 20: 1090–1094.
- 30. Yacoub Wasef SZ (2004) Gender differences in systemic lupus erythematosus. Gend Med 1: 12–17.
- 31. Coornaert B, Carpentier I, Beyaert R (2009) A20: central gatekeeper in inflammation and immunity. J Biol Chem 284: 8217–8221.
- 32. Graham RR, Cotsapas C, Davies L, Hackett R, Lessard CJ, et al. (2008) Genetic variants near TNFAIP3 on 6q23 are associated with systemic lupus erythematosus. Nat Genet 40: 1059–1061.
- 33. Vereecke L, Beyaert R, van Loo G (2009) The ubiquitin-editing enzyme A20 (TNFAIP3) is a central regulator of immunopathology. Trends Immunol 30: 383–391.
- 34. Lu LJ, Wallace DJ, Ishimori ML, Scofield RH, Weisman MH (2010) Review: Male systemic lupus erythematosus: a review of sex disparities in this disease. Lupus 19: 119–129.
- 35. Kariuki SN, Moore JG, Kirou KA, Crow MK, Utset TO, et al. (2009) Age- and gender-specific modulation of serum osteopontin and interferon-alpha by osteopontin genotype in systemic lupus erythematosus. Genes Immun 10: 487–494.