Abstract

Polymorphisms in genes involved in toll-like receptor/interferon signalling pathways have been reported previously to be associated with SLE in many populations. This study aimed to investigate the role of seven single nucleotide polymorphisms within TNFAIP3, STAT4, and IRF5, which are involved in upstream and downstream pathways of type I interferon production, in SLE in the South East Asian populations. Genotyping of 360 Malaysian SLE patients and 430 normal healthy individuals revealed that minor alleles of STAT4 rs7574865 and rs10168266 were associated with elevated risk of SLE in the Chinese and Malay patients, respectively (, odds ratio ; , , respectively). Polymorphisms in TNFAIP3 and IRF5 did not show significant associations with SLE in any of the ethnicities. Combined analysis of the Malays, Chinese, and Indians for each SNP indicated that STAT4 rs10168266 was significantly associated with the Malaysian SLE as a whole (; ). The meta-analysis of STAT4 rs10168266, which combined the data of other studies and this study, further confirmed its importance as the risk factor for SLE by having pooled OR of 1.559 and value of <0.001.

1. Introduction

Systemic lupus erythematosus (SLE) is a prototypic autoimmune disease affecting various parts of the body including skin, kidneys, lungs, joints, heart, nervous system, and hematopoietic organs. It is a disease whereby a diverse array of autoantibody production, complement activation, immune complex deposition, and inflammation cause damages in those organs. Although the exact aetiology of SLE still remains unclear, a combination of genetic risk factors and environmental events is believed to contribute to an irreversible break in immunological self-tolerance. With the introduction of genome-wide association studies, a huge breakthrough has been made in the discovery of SLE associated susceptibility genes that in turn advances our understanding of pathogenesis of SLE. Recently, several reviews have categorised the susceptible genes according to their immunological pathways and cell types. Three biological pathways involved in SLE have been forwarded by Harley et al. [1]: (i) innate immune response including toll-like receptor (TLR)/interferon (IFN) signalling pathways; (ii) adaptive immune response including B, T, and antigen-presenting cells immune signal transduction; and (iii) immune complex clearance mechanism.

Defects in TLR/IFN signalling pathways cause immune complexes containing self-nucleic acids to interact with TLR7 and TLR9 inside plasmacytoid dendritic cells and B cells endosomes, resulting in the secretion of type I IFN and interleukin (IL)-6. The combined triggering of both B cell receptors and TLR leads to autoreactive B-cell proliferation. Their further differentiation into plasmablasts and autoantibody-secreting plasma cells is induced by type I IFN and IL-6, respectively [2]. TNFAIP3, STAT4, IRF5, TREX, and IRAK1 are the genes involved in upstream and downstream pathways of type I IFN production that have been recently identified. The STAT4 gene consists of 24 exons that spread over a 120 kb region on chromosome 2q32.3. It encodes a transcription factor that mediates signals induced by IL-12, IL-23, and type I IFN and activates the production of IFN-γ and IL-17. It also directs the differentiation of helper T cells toward the proinflammatory T-helper type 1 and T-helper type 17 lineages that have been shown to play a critical role in the pathogenesis of SLE. The STAT4 null allele in lupus-prone mouse model confers reduced autoantibody production and glomerulonephritis, indicating that STAT4 may be involved in multiple SLE-associated phenotypes [3]. There are a few studies involving STAT4-deficient lupus-prone mice which demonstrate the role of STAT4 in autoantibody production only [4, 5]. Polymorphisms in the STAT4 gene have been found to be strongly associated with SLE susceptibility, in particular rs7574865 [6, 7]. The simultaneous association of the risk allele T of STAT4 rs7574865 with both lower serum IFN-α activity and increased IFN-α-induced gene expression has been reported, confirming that this polymorphism was associated with increased IFN-α sensitivity [8, 9].

TNFAIP3, or tumour necrosis factor alpha-induced protein 3 gene, encodes the A20 protein which is a negative regulator of the NF-κB signalling pathway, an essential pathway in the pathogenesis of SLE. A20 is an ubiquitin-editing enzyme required for effective termination of NF-κB-mediated proinflammatory responses induced by TLRs, TNF receptor, IL-1 receptor, and NOD2 [10]. A meta-analysis and imputation study identified a 109 kb risk haplotype spanning TNFAIP3 region with lupus nephritis and hematologic manifestation [11]. A nonsynonymous mutation (), rs223092, in TNFAIP3 gene which causes phenylalanine-to-cysteine change at position 127 of A20 protein has been consistently linked with SLE various ethnic groups.

The final candidate gene, IRF5, which is IFN regulatory factor 5, is a transcription factor that mediates inflammatory and immune responses [12]. This factor stimulates the production of the proinflammatory cytokines TNF-α, IL-12, and IL-6 following TLR signalling as well as transactivation of type I IFN and IFN-induced genes [13, 14]. Polymorphisms in IRF5 cause functional changes in messenger RNA, which in turn alter IFR5-mediated transcription resulting in elevated SLE risk [15]. It was also suggested that SLE patients who carry IRF5 risk haplotype and are positive for either anti-RBP or anti-dsDNA potentially have higher serum IFN-α activity [16]. In this study, we aimed to investigate the association between seven single nucleotide polymorphisms (SNPs) in STAT4, TNFAIP3, IRF5 genes, and SLE in the South East Asian scenario, particularly in the Malaysian participants. We also attempted to compare and pool the ORs of SNPs which were significant in the Malaysian SLE with the other studies through meta-analysis.

2. Materials and Methods

2.1. Sample Collection and DNA Extraction

A total of 790 Malaysians were included in this study, which is comprised of 360 SLE patients and 430 healthy controls. Blood samples were collected from patients diagnosed with SLE according to 4 out of ACR criteria and healthy volunteers recruited at the University of Malaya Medical Centre (UMMC), Kuala Lumpur, in compliance with requirements as stipulated by the UMMC Medical Ethics Committee (UMMC Ethics Approval Code: 733.19). The distribution of samples from Malays, Chinese, and Indians, as well as the ratio of females to males, is shown in Table 1. Genomic DNA was isolated from the peripheral blood samples by using the standard DNA extraction method as described previously [17]. The concentration and purity of the extracted DNA were further quantified by measuring the absorbance values at 260 nm and 280 nm via a spectrophotometer.

2.2. Genotyping with Tetraprimer ARMS-PCR

SNPs that were included in this study are listed in Table 2. Tetraprimer ARMS-PCR was performed in the genotyping of rs10168266 and rs7601754 in STAT4 region, rs2230926 and rs3757173 in TNFAIP3 region, and rs4728142 in IRF5 region. Primers were designed using computer software accessible through the Internet at http://cedar.genetics.soton.ac.uk/public_html/primer1.html, developed by Ye and team [18]. In silico PCR as described previously was further carried out to ensure the self-designed primers were targeted to the gene regions of interest [1921]. Each PCR reaction was carried out in a total of 10 μL, containing 50 ng of template DNA, appropriate concentration of inner and outer primers and MgCl2 (Table 3), 200 μM dNTP, 20 mM Tris-HCl pH8.4, 50 mM KCl, and 0.15 U Taq polymerase (Fermentas, Vilnius, Lithuania.). The PCR mixture was then subjected to touchdown PCR, whereby it was incubated for 5 min at 95°C, followed by 30 cycles of 45 s denaturation at 95°C, 45 s of annealing (started at temperature 10°C higher than annealing temperature, decreasing by 1°C per cycle, maintained at annealing temperature for the remaining 20 cycles) and 45 s of extension at 72°C, and a final extension at 72°C for 10 min at the end of the cycles. The annealing temperatures for different PCRs are stated in Table 3. Five microlitres of PCR amplicons was electrophoresed on a 2% (w/v) agarose gel. The agarose gel was viewed under UV illumination and image was recorded using a gel documentation system. The results obtained were further verified by sequencing.

2.3. Genotyping with Real-Time PCR

Predesigned TaqMan SNP genotyping assays were used to genotype SNPs where tetraprimers could not be designed for ARMS-PCR (probe ID: rs7574865 in STAT4 region, C_29882391_10; rs729302 in IRF5 region, C_2691216_10; Applied Biosystems, NY, USA). Fifty nanograms of template DNA was mixed with 2X Taqman GTXpress master mix (Applied Biosystems) and 20X Taqman genotyping assay (Applied Biosystems) to make up to a total volume of 10 μL. Real-time PCR reaction was initiated with pre-PCR read step at 60°C for 1 min, followed by DNA polymerase activation at 95°C for 20 s, 40 cycles of denaturation (95°C for 3 s) and annealing/extension (60°C for 30 s), and ended with a final extension step at 60°C for 1 min. Fluorescence was detected using an Applied Biosystems 7500 Fast Real-Time PCR System. The results were verified by sequencing.

2.4. Association Test

Allele and genotype frequencies were calculated, followed by performing a goodness-of-fit test to evaluate whether or not the observed genotype frequencies of each polymorphisms were departures from Hardy-Weinberg equilibrium (HWE) in control subjects (  values > 0.05). An open access HWE calculator developed by Rodriguez et al. [22] was used. Fisher’s exact test was conducted on 2 × 2 contingency table using SPSS software to assess the association of each SNP with SLE susceptibility in Malays, Chinese, and Indians. values were adjusted according to Bonferroni correction and was regarded as significant. Odds ratios (ORs) with 95% confidence intervals (CIs) were calculated. Adjusted ORs were computed using logistic regression, whereby major allele and major homozygous genotype of each SNP were set as reference group and their ORs were adjusted to 1.

The ORs of the three ethnicities were combined using Mantel-Haenszel test to evaluate the overall association of each SNP in the Malaysian population. The Mantel-Haenszel ORs were calculated using Comprehensive Meta-Analysis Version 2.0 software (Biostat, NJ, USA). The between-subgroup heterogeneity was tested using Cochran’s Q statistic. Random effect model was used when heterogeneity was significant ; otherwise, the fixed effect model was used.

2.5. Meta-Analysis

Meta-analysis was conducted using Comprehensive Meta-Analysis Version 2.0 software (Biostat) for the SNP(s) which was/were significantly associated with SLE in the Malaysian population by including data from other studies as well as the current study. Due to insufficient number of Indian subjects, analysis was not carried out for this ethnic group.

We examined the association between STAT4 polymorphisms and SLE fully and rigorously, with the use of the key words “STAT4,” “polymorphisms,” “systemic lupus erythematosus,” and “SLE.” Electronic databases including Pubmed, Embase, and Web of Science were thoroughly searched until December 2013. Only fully published articles were included and the eligible studies were identified based on the following criteria: (a) the study was original, (b) the patients were sporadic cases, (c) having available allele and genotype frequency data, and (d) having sufficient published data to determine OR with 95% CI.

Data extraction was performed by collecting the following information from each study: the first author’s name, year of publication, ethnicity, the number of cases and controls, and the frequency of minor allele (MAF) of each polymorphism in both cases and controls. For studies including several independent case-control populations, each case-control population was extracted separately. Malay and Chinese populations from this study were also included in the meta-analysis.

The heterogeneity across studies was evaluated by using Cochran’s Q statistic. Random effect model was used for meta-analysis when heterogeneity was significant (); otherwise, fixed effect model was used. By inputting the study name, total number of cases and controls, and MAF of cases and controls, Comprehensive Meta-Analysis Version 2.0 software calculated the ORs, 95% CI, and values for each study, as well as the pooled OR, 95% CI, and values for the meta-analysis.

3. Results

3.1. Polymorphisms and SLE Risk

The goodness-of-fit test demonstrated that all polymorphisms investigated in this study fulfilled HWE in the control group. The Malay SLE patients were significantly associated with minor alleles of STAT4 rs7574865, rs10168266, and TNFAIP3 rs2230926, and heterozygous genotype TG of TNFAIP3 rs2230926 (Table 4). Increased SLE susceptibility in Chinese population was significantly conferred by minor alleles and minor homozygous genotypes TT of STAT4 rs7574865 and rs10168266. However, after Bonferroni adjustment, significant associations were only observed between minor allele T of STAT4 rs7574865 and Chinese SLE patients (, , 95% CI: 1.12–1.82) and between minor allele T of STAT4 rs10168266 and Malay SLE patients (, , 95% CI: 1.20–2.71). STAT4 rs7601754, TNFAIP3 rs3757173, and IRF5 rs4728142 and rs729302 did not show significant association with SLE in any of the ethnicities.

Combined analysis of the three ethnicities was carried out to represent the association of each SNP with SLE in the Malaysian population as a whole. The analysis revealed that only minor allele T of STAT4 rs10168266 was significantly associated with the Malaysian SLE (, , 95% CI: 1.143–1.802) (Table 5).

3.2. Meta-Analysis of STAT4 rs10168166

Since STAT4 rs10168166 showed significant association with the Malaysian SLE, meta-analysis was carried out to combine the data from other studies with the current one [6, 2325]. Four relevant articles were identified eligible and a total of 5 subgroups were included for comparison. Four subgroups were from Asian population, while one was from European population. Data extracted from these articles is shown in Table 6.

In the overall analysis, significant association of STAT4 rs10168166 with SLE was observed. The fixed effect model was used as the heterogeneity test did not appear significant (). The pooled OR for the minor allele T was 1.559, with 95% CI of 1.459–1.665 and value of <0.001 (Figure 1).

4. Discussion

Understanding the full molecular pathology of SLE remains a great challenge, although many insights have been revealed. Recognition of self-nucleic acids by toll-like receptors TLR7 and TLR9 on plasmacytoid dendritic cells and B cells is believed to be an important step in the pathogenesis of this disease [26]. Increased antinuclear antibodies and production of type I IFN are both correlated with the severity of disease. STAT4, TNFAIP3, and IRF5 are genes involved in regulating TLR/IFN signalling pathways. SNPs investigated in this study have consistently shown associations with SLE susceptibility in many populations, especially in Asians [7, 24, 2729]. When the various ethnic groups were considered, rs7574865 and rs10168266 of STAT4 gene were significant in Chinese and Malays, respectively. However, only rs10168266 of STAT4 was observed to have correlations with SLE in the Malaysians generally. None of the SNPs seemed to influence SLE in Indians. Due to population demographics and lower SLE risk predisposition, the fewer Indians recruited in this study may have impacted the results. The SNPs of IRF5 were not significant. This suggests that the IRF5 genetic variants tested for in this study are not linked to SLE in our cohort and that there may be other variants that are more important in the Indian ethnic group.

Therefore, it may be concluded from the present study that STAT4 gene polymorphisms feature more prominently as the genetic risk factors in the Malaysian SLE rather than those polymorphisms in TNFAIP3 and IRF5. Rs10168266 which is located in intron 5 of STAT4 gene has been frequently related to SLE susceptibility in the Asian population, particularly in Korean population, and also in the European population [6, 2325, 30]. This was also reflected in the findings of this study. Nevertheless, not many studies were done on this SNP and thus only four studies were included in the meta-analysis of this study. After the analysis, the pooled OR and value once again showed that this SNP was overall an important risk factor for SLE and more attention should be taken.

Rs7574865, which was significantly associated with SLE in Chinese in this study, is located in the third intron of the STAT4 gene. The minor/risk allele T has reported associations with other immune-mediated diseases such as rheumatoid arthritis, primary Sjögren’s syndrome, type-1 diabetes, Crohn’s disease, and ulcerative colitis [3134]. The association of this particular SNP with SLE susceptibility was observed in many populations, including both European and Asian populations [7]. SNP haplotype in the third intron of STAT4 marked by rs7574865 was found to be associated with SLE susceptibility and it could be responsible for splice variation or regulatory effects of STAT4 [31, 35].

The next SNP that may be important in the Malaysian SLE was TNFAIP3 rs2230936. This coding SNP is a nonsynonymous variant causing a phenylalanine-to-cysteine change at residue 127 of the A20 protein. It has been already proven that minor Cys127 is relatively stable compared to the Phe127 protein, causing it to be less effective at inhibiting TNF-induced NF-κB activity [36]. This reduced autoinflammatory activity of A20 could result in excessive cellular response to TNF. Interestingly, as opposed to other findings suggesting that minor allele G was the risk factor of SLE, the results of this study demonstrated that it conferred protection against SLE in our cohort [11, 24, 27, 3640]. This study speculates that apart from rs223093, other factors such as adjacent SNPs may possibly alter the structure of A20 protein. This may play a role in SLE susceptibility in the Malays. We suggest that multiple (at least four) genes may collectively play critical roles in the development of this disease [4143]. This speculation has yet to be validated.

Finally, both SNPs in IRF5 gene investigated in our study were not significant in the Malaysian patients although both are fairly established SLE risk factors for Europeans and some Asians [4447]. Presumably, other SNPs of this gene would have to be considered.

5. Conclusion

The present study was relatively small in contrast to larger studies of SLE by other researchers. Nevertheless, we present evidence to suggest that the genes involved in TLR/IFN signalling pathways, especially STAT4 rs10168266 polymorphisms, contribute to the development of SLE in Malays and Chinese.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors would like to thank all the participants and nurses at University of Malaya Medical Centre (UMMC) who assisted in sample collection, especially Madam Julita. They are also grateful to physicians at the Renal Unit of UMMC for providing full cooperation during patient recruitment and sample collection and colleagues for helpful discussions. The study was funded by Monash University (Sunway Campus) Major Grant (no.: 5140060) awarded to Maude Elvira Phipps, and University of Malaya Research Grant (UMRG) (no.: RG018-09HTM) awarded to Kek Heng Chua.