A Multilocus Genetic Study in a Cohort of Italian SLE Patients Confirms the Association with STAT4 Gene and Describes a New Association with HCP5 Gene

Abstract

Background

Systemic lupus erythematosus (SLE) is an autoimmune disease with complex pathogenesis in which genes and environmental factors are involved. We aimed at analyzing previously identified loci associated with SLE or with other autoimmune and/or inflammatory disorders (STAT4, IL10, IL23R, IRAK1, PSORS1C1, HCP5, MIR146a, PTPN2, ERAP1, ATG16L1, IRGM) in a sample of Italian SLE patients in order to verify or confirm their possible involvement and relative contribution in the disease.

Materials and methods

Two hundred thirty-nine consecutive SLE patients and 278 matched healthy controls were enrolled. Study protocol included complete physical examination, and clinical and laboratory data collection. Nineteen polymorphisms were genotyped by allelic discrimination assays. A case-control association study and a genotype-phenotype correlation were performed.

Results

STAT4 was the most associated gene [P = 3×10⁻⁷, OR = 2.13 (95% CI: 1.59–2.85)]. IL10 confirmed its association with SLE [rs3024505: P = 0.02, OR = 1.52 (95% CI: 1.07–2.16)]. We describe a novel significant association between HCP5 locus and SLE susceptibility [rs3099844: P = 0.01, OR = 2.06 (95% CI: 1.18–3.6)]. The genotype/phenotype correlation analysis showed several associations including a higher risk to develop pericarditis with STAT4, and an association between HCP5 rs3099844 and anti-Ro/SSA antibodies.

Conclusions

STAT4 and IL10 confirm their association with SLE. We found that some SNPs in PSORS1C1, ATG16L1, IL23R, PTPN2 and MIR146a genes can determine particular disease phenotypes. HCP5 rs3099844 is associated with SLE and with anti-Ro/SSA. This polymorphism has been previously found associated with cardiac manifestations of SLE, a condition related with anti-Ro/SSA antibodies. Thus, our results may provide new insights into SLE pathogenesis.

About the Authors:

Cinzia Ciccacci

Affiliation: Department of Biomedicine and Prevention, Section of Genetics, School of Medicine, University of Rome “Tor Vergata”, Rome, Italy

Carlo Perricone

Affiliation: Lupus Clinic, Reumatologia, Dipartimento di Medicina Interna e Specialità Mediche, Sapienza Università di Roma, Rome, Italy

Fulvia Ceccarelli

Affiliation: Lupus Clinic, Reumatologia, Dipartimento di Medicina Interna e Specialità Mediche, Sapienza Università di Roma, Rome, Italy

Sara Rufini

Affiliation: Department of Biomedicine and Prevention, Section of Genetics, School of Medicine, University of Rome “Tor Vergata”, Rome, Italy

Davide Di Fusco

Affiliation: Department of Biomedicine and Prevention, Section of Genetics, School of Medicine, University of Rome “Tor Vergata”, Rome, Italy

Cristiano Alessandri

Affiliation: Lupus Clinic, Reumatologia, Dipartimento di Medicina Interna e Specialità Mediche, Sapienza Università di Roma, Rome, Italy

Francesca Romana Spinelli

Affiliation: Lupus Clinic, Reumatologia, Dipartimento di Medicina Interna e Specialità Mediche, Sapienza Università di Roma, Rome, Italy

Enrica Cipriano

Affiliation: Lupus Clinic, Reumatologia, Dipartimento di Medicina Interna e Specialità Mediche, Sapienza Università di Roma, Rome, Italy

Giuseppe Novelli

Affiliation: Department of Biomedicine and Prevention, Section of Genetics, School of Medicine, University of Rome “Tor Vergata”, Rome, Italy

Guido Valesini

Affiliation: Lupus Clinic, Reumatologia, Dipartimento di Medicina Interna e Specialità Mediche, Sapienza Università di Roma, Rome, Italy

Paola Borgiani

Contributed equally to this work with: Paola Borgiani, Fabrizio Conti

* E-mail: [email protected]

Affiliation: Department of Biomedicine and Prevention, Section of Genetics, School of Medicine, University of Rome “Tor Vergata”, Rome, Italy

Fabrizio Conti

Contributed equally to this work with: Paola Borgiani, Fabrizio Conti

Affiliation: Lupus Clinic, Reumatologia, Dipartimento di Medicina Interna e Specialità Mediche, Sapienza Università di Roma, Rome, Italy

Introduction

Systemic lupus erythematosus (SLE) is a chronic autoimmune disease characterized by the production of autoantibodies against intracellular, nucleic acid and cell surface antigens. The disease is relapsing and affects multiple organ systems including the skin, kidneys, lungs, heart and the central nervous system. The pathogenesis of SLE is complex, with several susceptibility genes and environmental factors involved in its development and clinical manifestations. Familial aggregation and a higher rate of concordance for SLE in monozygotic than in dizygotic twins [1], [2] provide support for a strong role of genetic factors in the pathogenesis of this disorder.

During last decade, many genes were described for their potential role in predisposing to SLE [3]. The list of involved genes has increased with the advent of Genome Wide Association (GWA) studies. These studies confirmed genetic associations of over 40 genes [3]. The number of genes emphasizes the complexity of the inheritability of SLE. Most of the identified genes encode proteins that participate in key pathogenic pathways, including Toll-like receptor and type I interferon signaling pathways, immune regulation pathways and those that control the clearance of immune complexes [4].

Interestingly, SLE and other autoimmune diseases share many of the risk loci identified, suggesting the involvement of common pathways in autoimmune diseases [5]. Moreover, the existence of common genetic susceptibility was already suggested by the clustering of multiple autoimmune diseases within families, including families with SLE [2], [6]. This genetic overlap is exemplified by the well-known associations of certain HLA loci with multiple autoimmune disease as well as non-HLA risk loci in different pathways such as IL2RA, STAT4, PTPN22 and IFIH1 [7], [8]. In this context, in a previous paper, we have described for the first time, an association between SLE and the TRAF3IP2 gene [9], a gene involved in psoriasis and psoriatic arthritis [10], but also described in cutaneous complications in inflammatory bowel diseases [11].

Here, our aim was to analyze a number of previously identified loci associated with SLE or with other autoimmune and/or inflammatory diseases (STAT4, IL10, IL23R, IRAK1, PSORS1C1, HCP5, MIR146a, PTPN2, ERAP1, ATG16L1, IRGM) in a sample of Italian SLE patients in order to confirm or verify their possible involvement and their relative contribution in SLE susceptibility. The genes were chosen due their role in different pathways involved in SLE pathogenesis. We therefore performed a case/control study and a genotype-phenotype correlation analysis.

Materials and Methods

Sample collection

Two hundred and thirty-nine SLE consecutive Caucasian patients ≥18 years of age affected with SLE were retrospectively enrolled at the Lupus Clinic of the Rheumatology Unit, Sapienza University of Rome (Sapienza Lupus Cohort). SLE diagnosis was performed according to the revised 1997 American College of Rheumatology (ACR) criteria. This cohort was already described in Perricone et al. 2013 [9]. Briefly, 21 (8.8%) were males and 218 (91.2%) females, mean age ± SD was 42.7±11.4 years, mean age at diagnosis ±SD was 30.9±11.6 years. Musculo-skeletal was the most frequent involvement (71.0% of patients), 96.6% had positive ANA, and 73.3% had positive anti-dsDNA. Study protocol included complete physical examination and blood drawing. The clinical and laboratory data were collected in a standardized, computerized, electronically-filled form including demographics, past medical history with date of diagnosis, co-morbidities, and previous and concomitant treatments.

The evaluation of clinical and laboratory parameters was performed at the Lupus Clinic of Sapienza University of Rome. Clinical and laboratory features were assessed with a dichotomous score (present = 1; absent = 0).

Written informed consent was obtained from each patient and the ethic committee of Sapienza University of Rome approved the study design.

Two hundred seventy-eight age- and ethnicity-matched healthy subjects, enrolled at University of Rome Tor Vergata, served as controls.

Clinical laboratory analysis

Among the laboratory tests, the following parameters were evaluated: complete blood count [CBC], including erythrocytes, leucocytes and platelets count, erythrocyte sedimentation rate [ESR], C-reactive protein [CRP] levels, serum levels of creatinine and complement C3 and C4, anti-nuclear (ANA), anti-double-stranded DNA (anti-dsDNA), ENA, anti-Sm, anti-RNP, anti-Ro/SSA, anti-La/SSB, anti-CL IgG and IgM, anti-β₂GPI IgG and IgM, and LA. To be positive, each of these tests had to be positive in at least 2 occasions 12 weeks apart all along patient clinical history.

ANA were determined by means of indirect immunofluorescence (IIF) on Hep-2 (titer ≥1:160 or ++ on a scale from + to ++++); anti-dsDNA with ELISA assays (considering levels higher than the cut-off of the reference laboratory) or IIF on Crithidia Luciliae, ENA (including anti-Ro/SSA, anti-La/SSB, anti-Sm, anti-RNP) by ELISA assay considering titers above the cut-off of the reference laboratory, anti-CL of IgG or IgM isotype, by ELISA, in serum or plasma at medium or high titers [e.g.,>40 standardized international units (GPL/MPL units) or above the 99^th percentile], anti-β₂GPI of IgG or IgM isotype, by ELISA, in serum or plasma (above the 99^th percentile), and finally LAC was assessed according to the guidelines of the International Society on Thrombosis and Hemostasis (Scientific Subcommittee on lupus anticoagulant/phospholipid-dependent antibodies) [12].

All subjects underwent blood drawing (5 ml supplemented with 0.5% EDTA) in order to perform genetic analysis.

DNA extraction and genotyping

Genomic DNA was isolated from peripheral blood mononuclear cells using a Qiagen blood DNA mini kit. A panel of SNPs was selected on the basis of literature data in 11 genes involved in immune response, autophagy and inflammation, already described as associated with any kind of autoimmune disease. The selected SNPs were the following: rs7574865 (STAT4), rs1800872 and rs3024505 (IL10), rs9469003 and rs3099844 (HCP5), rs11209026, rs11803505 and rs7517847 (IL23R), rs2233945 (PSORS1C1), rs13361189 and rs4958847 (IRGM), rs2241880 (ATG16L1), rs2542151 and rs7234029 (PTPN2), rs30187 and rs27524 (ERAP1), rs2910164 and rs2431697 (MIR146A) and rs3027898 (IRAK1). Genotyping was performed by allelic discrimination assay by TaqMan assays (Applied Biosystems, Foster City, CA, USA) and ABI PRISM 7000. Each assay was run with positive (samples previously confirmed by direct sequencing as heterozygous and/or variant homozygous) and negative controls.

Statistical analysis

The Hardy–Weinberg equilibrium was verified for all SNPs by the Pearson χ² test.

Differences in alleles/genotypes/haplotypes frequencies between cases and control were evaluated by Pearson χ² test. Haplotypes were inferred using Arlequin, version 3.11 [13]. Odds ratios (ORs) with 95% CI were calculated. A genotype-phenotype correlation analysis has been performed considering the heterozygotes and variant homozygotes together (1degree of freedom [df]). A binary logistic regression analysis (stepwise) was carried out considering the presence/absence of SLE as dependent variable and all genetics associated factors as independent variables. All statistical analyses were performed by the SPSS 13.0 program (SSPS Inc., IL, USA).

Results

Case/control association analysis

We studied 19 polymorphisms in 11 genes. Deviations from Hardy–Weinberg equilibrium were not observed. In Table 1, we reported the comparisons of the allele frequencies between cases and controls for all studied polymorphisms. Among the studied genes, STAT4 resulted, as expected, the most associated with a P = 3×10⁻⁷ and OR = 2.13 (95% CI 1.59–2.85). Moreover, the Mantel–Haenszel test for linear association highlights that there was an additive effect of STAT4 risk allele (P = 5.3×10⁻⁷): subjects with one risk allele presented an OR = 1.57, while subjects with two risk alleles reached an OR = 8.4 (95% CI 3.42–20.72, Table 2). Furthermore, IL10 confirmed the association with SLE in our cohort with the rs3024505 SNP (P = 0.02 and OR = 1.52 [95% CI 1.07–2.16]).

[Figure omitted. See PDF.]

Table 1. Results of allelic association analysis.

https://doi.org/10.1371/journal.pone.0111991.t001

[Figure omitted. See PDF.]

Table 2. Genotypes distribution for rs7574865 STAT4 SNP revealed an allele dose additive effect.

https://doi.org/10.1371/journal.pone.0111991.t002

Moreover, for the first time, we described an association between HCP5 locus and SLE susceptibility: the variant allele of rs3099844 resulted a risk allele with P = 0.01 and OR = 2.06 (95% CI 1.18–3.6). No significant associations were observed for all the other studied SNPs.

A binary logistic regression analysis was carried out considering the three associated SNPs as independent variables and the presence/absence of SLE as dependent variable (Table 3). This analysis confirmed the involvement of STAT4 and HCP5. This model explains the 5.9% of genetic susceptibility (Cox and Snell R²).

[Figure omitted. See PDF.]

Table 3. Final model of binary logistic regression analysis (stepwise method) considering only associated genetic factors.

https://doi.org/10.1371/journal.pone.0111991.t003

Haplotypes analysis

We inferred the haplotypes distribution for IL10, IL23R, HCP5, IRGM, PTPN and ERAP1 genes. However, we did not find any significant association comparing the distribution in cases and controls (data not shown). Only for IL10 we observed a significant difference comparing the TC (composed by the variant allele of rs3024505 and wildtype allele of rs1800872) vs CA (composed by the wild-type allele of rs3024505 and the variant allele of rs1800872) haplotype (P = 0.012 and OR = 1.7 [95% CI 1.12–2.56]), but this evidence does not improve consistently the statistical significance obtained considering only the rs3024505 SNP.

Number of risk alleles

We counted the number of risk alleles, considering as risk alleles the variant alleles of rs7574865 (STAT4), rs3024505 (IL10), and rs3099844 (HCP5) SNPs, and we evaluated the different risk allele number distribution between cases and controls (Figure 1, Table 4). As expected the class with no risk allele was significantly more present in controls than in patients (P = 0.0008, OR = 0.45). On the contrary, the class with two risk alleles was significantly more present in cases than in controls (P = 0.0003, OR = 3.78). Finally, genetic profiles with three or more risk alleles were present only in SLE patients (P<0.0001).

[Figure omitted. See PDF.]

Figure 1. Counting of risk alleles in cases and controls.

The risk alleles considered are the variant alleles of rs7574865, rs3099844 and rs3024505 SNPs.

https://doi.org/10.1371/journal.pone.0111991.g001

[Figure omitted. See PDF.]

Table 4. Counting of risk alleles in cases and controls.

https://doi.org/10.1371/journal.pone.0111991.t004

Genotype-phenotype correlation analysis

A genotype/phenotype correlation analysis between the genotypes and clinical characteristics of SLE patients was performed. Results of this analysis were reported in Table 5. STAT4, that showed the strongest association with the disease, was also associated with a higher risk to develop pericarditis, while HCP5 SNPs resulted associated with anti-Ro/SSA, anti-dsDNA and aCL.

[Figure omitted. See PDF.]

Table 5. Association between disease phenotypes and studied polymorphisms.

https://doi.org/10.1371/journal.pone.0111991.t005

Interestingly, some SNPs that did not associate with the disease, showed association with particular disease phenotypes. For example, PTPN2 SNP resulted associated with a higher risk to develop pericarditis and serositis, while ATG16L1 SNP resulted protective towards pericarditis. IL23R SNPs resulted associated with a higher risk to develop aphthous/ulcer and with a minor risk to develop malar rash. MIR146A SNP resulted associated with a higher risk to develop arthritis and with C3 consumption. PSORS1C1 resulted associated with the presence of LAC and decreased C4.

Moreover, PSORS1C1 resulted also associated with age at diagnosis (table 6): indeed the average age at onset is significantly lower in patients with rs2233945 SNP variant (28.6 vs 32.2 years, P = 0.042).

[Figure omitted. See PDF.]

Table 6. PSORS1C1 SNP correlates with age at onset.

https://doi.org/10.1371/journal.pone.0111991.t006

Discussion

In the last years, the advent of genome wide association studies and the technological possibility to detect polymorphisms in large samples have revolutionized the field of human genetics by identification and confirmation of common gene variants that confer susceptibility to a plethora of diseases. A number of studies support the concept that common risk genes underlie multiple autoimmune disorders and suggest the involvement of common pathways of pathogenesis among different diseases [14]. In this study, we have examined specific genes, involved in autoimmunity, in inflammation and in autophagy to evaluate their involvement in SLE susceptibility. We aimed at confirming the association between STAT4, IL10 and IRAK1 [15]–[19] in a monocentric Italian cohort with SLE. Other genes, such as IL23R, PSORS1C1, HCP5, were found to be associated with psoriasis and rheumatoid arthritis [20]–[23]. Regarding genes involved in autophagy, we have chosen IRGM and ATG16L1 as candidate genes, since they were associated with Crohn's disease [24]. PTPN2 was associated with Crohn's disease [25] and type I diabetes [26]; ERAP1 was associated with psoriasis [27] and multiple sclerosis [28].

In our study, only three genes resulted associated with the disease development: STAT4, IL10 and HCP5. Indeed, we confirmed the STAT4 rs7574865 variant allele as a risk factor for SLE susceptibility, with an OR even higher than previously described (OR = 2.13 [95% CI 1.59–2.85]). Moreover, we observed an allele-dose risk effect, whereby individuals with two STAT4 alleles have 8.4-fold risk to develop the disease compared with wild-type subjects, while heterozygous subjects show 1.57-fold risk. Considering the fact that, usually, each risk factor gives a moderate contribution to SLE susceptibility, in our cohort, this SNP has a considerable weight. STAT4 is a member of the STAT family of transcription factors, and is involved in mediating responses to IL12 in lymphocytes, and regulating the differentiation of T helper cells. The STAT4 rs7574865 polymorphism, located within the third intron of the gene, shows the strongest association with several autoimmune conditions [29], [30], such as SLE, RA, type 1 and 2 diabetes, systemic sclerosis, inflammatory bowel diseases, primary Sjogren's syndrome (SS), juvenile idiopathic arthritis, primary antiphospholipid syndrome, autoimmune thyroid diseases, multiple sclerosis, psoriasis, granulomatosis with polyangiitis, giant cell arteritis [31].

The STAT4 gene can be activated and phosphorylated upon ligation of the type I IFN receptor by IFN-α, and subsequently induce downstream transcription of IFN-α-induced genes [32]. It has been recently showed that a great number of genes target for STAT4 are enriched to functional pathways in the type I interferon system, and have key roles in the inflammatory response. Thus, by participating in transcription complex with other co-factors, STAT4 harbors the potential of regulating a large number of target genes, which may contribute to their strong association with SLE [33].

The other interesting result is the association between a gene variant in HCP5 gene and SLE susceptibility that was never described before. Indeed, this is the first study, by our knowledge, that evaluated the involvement of this gene in SLE susceptibility. Moreover, HCP5 SNPs resulted associated with anti-Ro/SSA and anti-dsDNA antibodies. The rs3099844 polymorphism has been previously described as involved in Steven Johnson Syndrome and Toxic Epidermal Necrolysis susceptibility [34]–[36], but also in another autoimmune condition, primary sclerosing cholangitis [37]. In this condition, anti-ENA including anti-SSA can be found in up to 11.5% of patients and anti-dsDNA at a lesser extent (1–3% of patients). However, the most interesting study is a genome-wide association study showing an association between the rs3099844 and the cardiac manifestations of systemic lupus erythematosus [38]. As it is well known, this clinical condition is associated with the presence of anti-Ro/SSA antibodies [39]–[43]. Altogether, these data suggest a strong evidence of an association of such polymorphism with the anti-Ro/SSA antibodies, and further studies are awaited to clarify the specific underlying mechanisms. Another SNP in the HCP5 gene (rs2395029) was previously associated with psoriasis and psoriatic arthritis [23]. HCP5 gene (major histocompatibility complex P5), located in HLA region, is expressed primarily in cells of the immune system such as spleen, blood and thymus consistent with a potential role in autoimmunity [23]. This allele was associated with lower viral load set point following HIV-1 infection. This is of interest, since anti-Ro/SSA antibodies have been described in HIV infected subjects [44], [45]. Hence, it is possible that HCP5 carriers mount a strong immune reaction to viral infection that, in genetically susceptible individuals, could lead to excessive inflammation. On the other hand, sera from SLE and SS patients were also found to present reactivity with HIV-1 p24 antigen [46], suggesting a specific pattern of nonrandom cross-reactivity between virus p24 and autoimmune sera. Moreover, both a direct interaction between the SSB and HIV-1 trans-activation response element as well as the involvement of the SSB autoantigen in HIV-1 gene expression have been reported previously [47], [48]. Thus, there exist a relationship between the anti-Ro/SSA, anti-La/SSB antibodies and HIV [49] of which HCP5 may represent the bridge.

Finally, besides the confirmation of STAT4 and the newly described association with HCP5, we confirmed the risk effect of IL10 rs3024505 variant allele on SLE development. Polymorphisms in this gene were described as associated with several autoimmune disease, including SLE, rheumatoid arthritis, psoriasis and Crohn's disease, highlighting that IL10 is one of the common risk genes. IL10, primarily produced by monocytes and lymphocytes, is a multifunctional cytokine in immune-regulation and inflammation. A number of studies have been conducted to investigate the association of the IL10 gene polymorphisms with SLE susceptibility, although with conflicting results [18]. Considering the three associated SNPs, we also observed a cumulative genetic effect that distinguish cases from controls (Figure 1). A similar analysis in a previous study, described a different genetic risk scores in males and females [50].

Regarding the other genes, although we did not find any association with the development of the disease, some of the examined polymorphisms resulted associated with particular phenotypes, suggesting that these genes probably do not affect the disease susceptibility but modulate the phenotype. Of course, this lack of association does not mean that these genes are not involved with the disease; indeed, we cannot exclude that other SNPs in these same genes, not investigated in our study, could be associated. Interestingly, PSORS1C1 gene, which maps in the HLA region, resulted associated with LAC, decreased C4, and with an early-onset of the disease. PSORS1C1 gene is one of psoriasis susceptibility loci [51], but it is also associated with systemic sclerosis [52] and has been suggested to play an important role in the development of RA [53]. We also described an association between MiR-146a and the development of arthritis complications. It is known that MiR-146a is a negative regulator of the IFN pathway, and under-expression of MiR-146a conduces to alterations in the type I IFN pathway in lupus patients by targeting the key signaling proteins [54]. Autophagy has been described as a critical homeostatic mechanism in T lymphocytes, influencing proliferation and differentiation [55], and there are studies in the literature which highlighted the association between SLE and autophagy genes [56], [57]. However, in our study, we did not find any association of the disease with ATG16L1 and IRGM genetic variants. Of course, this lack of association could be due to the SNPs selection and other SNPs in the same gene could associate with the disease. Nonetheless, it is interesting to note that rs2241880 in ATG16L1 resulted associated with a minor risk to develop pericarditis as a complication. Admittedly, a limitation of our study is the relatively small number of samples; however, it should be underlined that this a monocentric cohort of patients of the same ethnicity, strictly followed and clinically well characterized.

In conclusion, our study confirms the role of STAT4 and IL10 in SLE susceptibility and describes for the first time the association between a HCP5 genetic variant and SLE susceptibility. These results confirm the importance of the adaptive and innate arms of the immune system in establishing SLE risk. However, each of the identified alleles accounts for only a fraction of the overall genetic risk, and functional studies are necessary to clearly explain their involvement with SLE pathogenesis.

Author Contributions

Conceived and designed the experiments: CC CP GN GV PB F. Conti. Performed the experiments: CC CP F. Ceccarelli SR DDF CA FRS EC. Analyzed the data: CC CP PB F. Ceccarelli. Contributed reagents/materials/analysis tools: CC SR DDF PB. Wrote the paper: CC CP PB F. Ceccarelli.

Citation: Ciccacci C, Perricone C, Ceccarelli F, Rufini S, Di Fusco D, Alessandri C, et al. (2014) A Multilocus Genetic Study in a Cohort of Italian SLE Patients Confirms the Association with STAT4 Gene and Describes a New Association with HCP5 Gene. PLoS ONE 9(11): e111991. https://doi.org/10.1371/journal.pone.0111991

References

1. Deapen D, Escalante A, Weinrib L, Horwitz D, Bachman B, et al. (1992) A revised estimate of twin concordance in systemic lupus-erythematosus. Arthritis Rheum 35: 311–318.

2. Alarcón-Segovia D, Alarcón-Riquelme ME, Cardiel MH, Caeiro F, Massardo L, et al. (2005) Grupo Latinoamericano de Estudio del Lupus Eritematoso (GLADEL) (2005) Familial aggregation of systemic lupus erythematosus, rheumatoid arthritis, and other autoimmune diseases in 1177 lupus patients from the GLADEL cohort. Arthritis Rheum 52: 1138–1147.

3. Cui Y, Sheng Y, Zhang X (2013) Genetic susceptibility to SLE: recent progress from GWAS. J Autoimmun 41: 25–33.

4. Deng Y1, Tsao BP (2010) Genetic susceptibility to systemic lupus erythematosus in the genomic era. Nat Rev Rheumatol 6: 683–92.

5. Delgado-Vega A, Sánchez E, Löfgren S, Castillejo-López C, Alarcón-Riquelme ME (2010) Recent findings on genetics of systemic autoimmune diseases. Curr Opin Immunol 22: 698–705

6. Ramos PS, Criswell LA, Moser KL, Comeau ME, Williams AH, et al. (2011) International Consortium on the Genetics of Systemic Erythematosus A comprehensive analysis of shared loci between systemic lupus erythematosus (SLE) and sixteen autoimmune diseases reveals limited genetic overlap. PLoS Genet 7: e1002406.

7. Gregersen PK, Olsson LM (2009) Recent advances in the genetics of autoimmune disease. Annu Rev Immunol 27: 363–391.

8. Perricone R, Perricone C, Shoenfeld Y (2011) Autoimmunity: when the immune system becomes the self-ish giant. Autoimmun Rev 10: 575–6.

9. Perricone C, Ciccacci C, Ceccarelli F, Di Fusco D, Spinelli FR, et al. (2013) TRAF3IP2 gene and systemic lupus erythematosus: association with disease susceptibility and pericarditis development. Immunogenetics 65: 703–9.

10. Hüffmeier U, Uebe S, Ekici AB, Bowes J, Giardina E, et al. (2010) Common variants at TRAF3IP2 are associated with susceptibility to psoriatic arthritis and psoriasis. Nat Genet 42: 996–9.

11. Ciccacci C, Biancone L, Di Fusco D, Ranieri M, Condino G, et al. (2013) TRAF3IP2 gene is associated with cutaneous extraintestinal manifestations in inflammatory bowel disease. J Crohns Colitis 7: 44–52.

12. Miyakis S, Lockshin MD, Atsumi T, Branch DW, Brey RL, et al. (2006) International consensus statement on an update of the classification criteria for definite antiphospholipid syndrome (APS). J Thromb Haemost 4: 295–306.

13. Excoffier L, Lischer HE (2010) Arlequin suite ver. 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour 10: 564–7.

14. Gregersen PK, Behrens TW (2006) Genetics of autoimmune disease — disorders of immune homeostasis. Nat Rev Genet 7: 917–28.

15. Zheng J, Yin J, Huang R, Petersen F, Yu X (2013) Meta-analysis reveals an association of STAT4 polymorphisms with systemic autoimmune disorders and anti-dsDNA antibody. Hum Immunol 74: 986–92.

16. Wang C, Sandling JK, Hagberg N, Berggren O, Sigurdsson S, et al. (2013) Genome-wide profiling of target genes for the systemic lupus erythematosus-associated transcription factors IRF5 and STAT4. Ann Rheum Dis 72: 96–103.

17. Schotte H, Willeke P, Becker H, Poggemeyer J, Gaubitz M, et al. (2014) Association of extended interleukin-10 promoter haplotypes with disease susceptibility and manifestations in German patients with systemic lupus erythematosus. Lupus 23: 378–85.

18. Liu P, Song J, Su H, Li L, Lu N, et al. (2013) IL-10 gene polymorphisms and susceptibility to systemic lupus erythematosus: a meta-analysis. PLoS One 8: e69547.

19. Kaufman KM, Zhao J, Kelly JA, Hughes T, Adler A, et al. (2013) Argentine Collaborative Group; BIOLUPUS network, Fine mapping of Xq28: both MECP2 and IRAK1 contribute to risk for systemic lupus erythematosus in multiple ancestral groups. Ann Rheum Dis 72: 437–44.

20. Song GG, Bae SC, Choi SJ, Ji JD, Lee YH (2012) Associations between interleukin-23 receptor polymorphisms and susceptibility to rheumatoid arthritis: a meta-analysis. Mol Biol Rep 39: 10655–63.

21. Zhai Y, Xu K, Huang F, Peng H, Feng CC, et al. (2012) Association of interleukin 23 receptor gene polymorphisms (rs10489629, rs7517847) with rheumatoid arthritis in European population: a meta-analysis. Mol Biol Rep 39: 8987–94.

22. Sun H, Xia Y, Wang L, Wang Y, Chang X (2013) PSORS1C1 may be involved in rheumatoid arthritis. Immunol Lett 153: 9–14

23. Liu Y, Helms C, Liao W, Zaba LC, Duan S, et al. (2008) A genome-wide association study of psoriasis and psoriatic arthritis identifies new disease loci. PLoS Genet 4: e1000041.

24. Latiano A, Palmieri O, Valvano MR, D'Incà R, Cucchiara S, et al. (2008) Replication of interleukin 23 receptor and autophagy-related 16-like 1 association in adult- and pediatric-onset inflammatory bowel disease in Italy. World J Gastroenterol 14: 4643–51.

25. Glas J, Wagner J, Seiderer J, Olszak T, Wetzke M, et al. (2012) PTPN2 gene variants are associated with susceptibility to both Crohn's disease and ulcerative colitis supporting a common genetic disease background. PLoS One 7: e33682.

26. Espino-Paisan L, de la Calle H, Fernández-Arquero M, Figueredo MA, de la Concha EG, et al. (2011) A polymorphism in PTPN2 gene is associated with an earlier onset of type 1 diabetes. Immunogenetics 63: 255–8.

27. Genetic Analysis of Psoriasis Consortium & the Wellcome Trust Case Control Consortium 2 (2010) Strange A, Capon F, Spencer CC, Knight J, et al. (2010) A genome-wide association study identifies new psoriasis susceptibility loci and an interaction between HLA-C and ERAP1. Nat Genet 42: 985–90.

28. Guerini FR, Cagliani R, Forni D, Agliardi C, Caputo D, et al. (2012) A functional variant in ERAP1 predisposes to multiple sclerosis. PLoS One 7: e29931.

29. Korman BD, Kastner DL, Gregersen PK, Remmers EF (2008) STAT4: genetics, mechanisms, and implications for autoimmunity. Curr Allergy Asthma Rep 8: 398–403.

30. Perricone C, Ceccarelli F, Valesini G (2011) An overview on the genetic of rheumatoid arthritis: a never-ending story. Autoimmun Rev 10: 599–608.

31. Liang YL, Wu H, Shen X, Li PQ, Yang XQ, et al. (2012) Association of STAT4 rs7574865 polymorphism with autoimmune diseases: a meta-analysis. Mol Biol Rep 39: 8873–82.

32. Kariuki SN, Niewold TB (2010) Genetic regulation of serum cytokines in systemic lupus erythematosus. Transl Res 155: 109–17.

33. Wang C, Sandling JK, Hagberg N, Berggren O, Sigurdsson S, et al. (2013) Genome-wide profiling of target genes for the systemic lupus erythematosus-associated transcription factors IRF5 and STAT4. Ann Rheum Dis 72: 96–103

34. Borgiani P, Di Fusco D, Erba F, Marazzi MC, Mancinelli S, et al. (2014) HCP5 genetic variant (rs3099844) contributes to Nevirapine-induced Stevens Johnsons Syndrome/Toxic Epidermal Necrolysis susceptibility in a population from Mozambique. Eur J Clin Pharmacol 70: 275–8

35. Tohkin M, Kaniwa N, Saito Y, Sugiyama E, Kurose K, et al. (2013) Japan Pharmacogenomics Data Science Consortium. A whole-genome association study of major determinants for allopurinol-related Stevens-Johnson syndrome and toxic epidermal necrolysis in Japanese patients. Pharmacogenomics J 13: 60–9.

36. Génin E, Schumacher M, Roujeau JC, Naldi L, Liss Y, et al. (2011) Genome-wide association study of Stevens-Johnson Syndrome and Toxic Epidermal Necrolysis in Europe. Orphanet J Rare Dis 6: 52.

37. Karlsen TH, Franke A, Melum E, Kaser A, Hov JR, et al. (2010) Genome-wide association analysis in primary sclerosing cholangitis. Gastroenterology 138: 1102–11.

38. Clancy RM, Marion MC, Kaufman KM, Ramos PS, Adler A (2010) International Consortium on Systemic Lupus Erythematosus Genetics, Identification of candidate loci at 6p21 and 21q22 in a genome-wide association study of cardiac manifestations of neonatal lupus. Arthritis Rheum 62: 3415–24.

39. Julkunen H, Kurki P, Kaaja R, Heikkila R, Immonen I, et al. (1993) Isolated congenital heart block. Long-term outcome of mothers and characterization of the immune response to SS-A/Ro and to SS-B/La. Arthritis Rheum 36: 1588–1598.

40. Salomonsson S, Dorner T, Theander E, Bremme K, Larsson P, et al. (2002) A serologic marker for fetal risk of congenital heart block. Arthritis Rheum 46: 1233–1241.

41. Fritsch C, Hoebeke J, Dali H, Ricchiuti V, Isenberg DA, et al. (2006) 52-kDa Ro/SSA epitopes preferentially recognized by antibodies from mothers of children with neonatal lupus and congenital heart block. Arthritis Res Ther 8: R4.

42. Eronen M, Miettinen A, Walle TK, Chan EK, Julkunen H (2004) Relationship of maternal autoimmune response to clinical manifestations in children with congenital complete heart block. Acta Paediatr 93: 803–809.

43. Salomonsson S, Dzikaite V, Zeffer E, Eliasson H, Ambrosi A, et al. (2011) A population-based investigation of the autoantibody profile in mothers of children with atrioventricular block. Scand J Immunol 74: 511–517.

44. Muller S, Richalet P, Laurent-Crawford A, Barakat S, Rivière Y, et al. (1992) Autoantibodies typical of non-organ-specific autoimmune diseases in HIV- seropositive patients. AIDS 6: 933–942.

45. Argov S, Schattner A, Burstein R, Handzel ZT, Shoenfeld Y, et al. (1991) autoantibodies in male homosexuals and HIV infection. Immunol Lett 30: 31–35.

46. Deas JE, Liu LG, Thompson JJ, Sander DM, Soble SS, et al. (1998) Reactivity of sera from systemic lupus erythematosus and Sjogren syndrome patients with peptides derived from human immunodeficiency virus p24 capsid antigen. Clin Diagn Lab Immunol 5: 181–185.

47. Svitkin YV, Pause A, Sonenberg N (1994) La autoantigen alleviates translational repression by the 5′ leader sequence of the human immunodeficiency virus type 1 mRNA. J Virol 68: 7001–7007.

48. Chang YN, Kenan DJ, Keene JD, Gatignol A, Jeang KT (1995) Direct interactions between autoantigen La and human immunodeficiency virus leader RNA. J Virol 69: 618–619.

49. Billaut-Mulot O, Cocude C, Kolesnitchenko V, Truong MJ, Chan EK, et al. (2001) SS-56, a novel cellular target of autoantibody responses in Sjögren syndrome and systemic lupus erythematosus. J Clin Invest 108: 861–869.

50. 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–9.

51. Chang YT, Liu HN, Shiao YM, Lin MW, Lee DD, et al. (2005) A study of PSORS1C1 gene polymorphisms in Chinese patients with psoriasis. Br J Dermatol 153: 90–6.

52. Allanore Y, Saad M, Dieudé P, Avouac J, Distler JH, et al. (2011) Genome-wide scan identifies TNIP1, PSORS1C1, and RHOB as novel risk loci for systemic sclerosis. PLoS Genet 7: e1002091.

53. Sun H, Xia Y, Wang L, Wang Y, Chang X (2013) PSORS1C1 may be involved in rheumatoid arthritis. Immunol Lett 153: 9–14

54. Tang Y, Luo X, Cui H, Ni X, Yuan M, et al. (2009) MicroRNA-146A contributes to abnormal activation of the type I interferon pathway in human lupus by targeting the key signaling proteins. Arthritis Rheum 60: 1065–75

55. Pierdominici M, Vomero M, Barbati C, Colasanti T, Maselli A, et al. (2012) Role of autophagy in immunity and autoimmunity, with a special focus on systemic lupus erythematosus. FASEB J 26: 1400–12

56. Zhou XJ, Lu XL, Lu JC, Yang HZ, Qin LX, et al. (2011) Genetic association of PRDM1-ATG5 intergenic region and autophagy with systemic lupus erythematosus in a Chinese population. Ann Rheum Dis 70: 1330–7.

57. Järvinen TM, Hellquist A, Zucchelli M, Koskenmies S, Panelius J, et al. (2012) Replication of GWAS-identified systemic lupus erythematosus susceptibility genes affirms B-cell receptor pathway signalling and strengthens the role of IRF5 in disease susceptibility in a Northern European population. Rheumatology (Oxford) 51: 87–92.