- Research Article
- Open Access
Association of STAT4 with Rheumatoid Arthritis in the Korean Population
Molecular Medicine volume 13, pages455–460(2007)
A recent study in the North American White population has documented the association of a common STAT4 haplotype (tagged by rs7574865) with risk for rheumatoid arthritis (RA) and systemic lupus erythematosus. To replicate this finding in the Korean population, we performed a case-control association study. We genotyped 67 single nucleotide polymorphisms (SNPs) within the STAT1 and STAT4 regions in 1123 Korean patients with RA and 1008 ethnicity-matched controls. The most significant four risk SNPs (rs11889341, rs7574865, rs8179673, and rs10181656 located within the third intron of STAT4) among 67 SNPs are identical with those in the North American study. All four SNPs have modest risk for RA susceptibility (odds ratio 1.21–1.27). A common haplotype defined by these markers (TTCG) carries significant risk for RA in Koreans [34 percent versus 28 percent, P = 0.0027, OR (95 percent CI) = 1.33 (1.10–1.60)]. By logistic regression analysis, this haplotype is an independent risk factor in addition to the classical shared epitope alleles at the HLA-DRB1 locus. There were no significant associations with age of disease onset, radiographic progression, or serologic status using either allelic or haplotypic analysis. Unlike several other risk genes for RA such as PTPN22, PADI4, and FCRL3, a haplotype of the STAT4 gene shows consistent association with RA susceptibility across Whites and Asians, suggesting that this risk haplotype predates the divergence of the major racial groups.
Rheumatoid arthritis (RA) is a chronic autoimmune arthritis characterized by progressive joint destruction and autoantibody formation such as anti-cyclic citrullinated peptide (anti-CCP) and rheumatoid factor (RF). Both genetic and environmental factors, and their interaction, play a role in the development of RA (1–3). A common set of alleles at the HLA-DRB1 locus (the “shared epitope” alleles) has been associated with RA in populations of White European and Asian ancestry (4,5). However, commonality of other risk loci in these population groups has been difficult to demonstrate.
Over the last several years, additional risk genes for RA have been identified in both White and Asian populations. PTPN22 was discovered as a risk factor for RA by genome-wide association scanning of functional SNPs (6) and has been replicated in many White RA cohorts. In addition, PTPN22 confers risk for several other autoimmune diseases such as Type1 diabetes and systemic lupus erythematosus (SLE), providing evidence for common pathways of pathogenesis in these disorders (7).
However, the PTPN22 risk allele (R620W) is extremely rare in the Asian population and, so far, there is no evidence of association of PTPN22 with RA in non-White populations (8,9). In contrast, PADI4, SLC22A4, and FCRL3 have been associated with RA in the Japanese population and replicated in the other Asian groups (10–12), but have given weak or negative results in populations of European ancestry (13,14). These divergent results suggest genetic heterogeneity of RA across the major racial groups (15).
Recently, a study in the North American White population has documented the association of a common STAT4 haplotype with both RA and SLE using a combined positional mapping and candidate-gene approach catalyzed by finding a linkage peak on chromosome 2q (16). In this report we demonstrate that a common STAT4 haplotype confers a similar degree of risk for RA in both Asian and White populations.
Materials and Methods
We included 1128 Korean RA patients who were enrolled consecutively from the outpatient clinic of The Hospital for Rheumatic Diseases, Hanyang University, Seoul, South Korea and 1022 ethnically matched controls. All RA patients met the American College of Rheumatology 1987 classification criteria for RA (17). Written informed consent was obtained from all participants. HLA-DRB1 typing was performed in cases (N = 1120) and controls (N = 1002) by a polymerase chain reaction sequence based typing (PCR-SBT) method using the reference protocol of the Twelfth International Histocompatibility Workshop (18).
Clinical and Laboratory Data of the Studied Subjects
Clinical data such as sex and age at disease onset from medical records and from interviews conducted at the time of enrollment were used in the analysis.
The staging system proposed by Steinbrocker et al was used as a marker of the radiographic severity of RA (19). Subjects were initially classified into stages I–IV and were then dichotomized into two groups: stage I and stages II–IV. Serum levels of rheumatoid factor (RF) were measured nephelometrically. Serum anti-CCP levels were quantitatively measured in duplicate from all RA patients by enzyme-linked immunosorbent assay using the DIASTAT Anti-CCP kit FCCP 200 (Axis-Shield Plc, Scotland, UK) according to the manufacturer’s instructions. The upper limit of the reference range of anti-CCP is five units and anti-CCP levels were not titered out beyond 100 units.
We first confirmed that the most significant risk SNP (rs7574865) of the STAT4 gene in the North American RA population was also polymorphic in Asian population (MAF 0.28 and 0.33 in Japanese and Chinese population, respectively) using HapMap Phase II data (http://www.hapmap.org). We therefore genotyped the same panel of 67 SNPs as in our previous report (16). These were selected from Hapmap Phase II data and included all non-synonymous coding SNPs reported in dbSNP, and covered both the STAT1 and STAT4 genes on chromosome 2q.
The DNAs were genotyped using a multiplexed primer extension method (Sequenom, Inc., San Diego, CA). In brief, multiplex PCR was used to amplify DNA products containing up to 28 SNPs in one reaction from 5 ng genomic DNA. Synthetic oligonucleotides that bind adjacent to the SNP site were then hybridized and extended with nucleotides complementary to the template SNP site using modified nucleotides that terminate the extension reaction at the interrogated SNP, thus generating alternate products of sufficiently different masses to be separated by mass spectrometry. The extended products were separated by MALDI-TOF mass spectrometry and the genotypes determined with SpectroTyper software (Sequenom). Calls were evaluated and edited by cluster analysis performed with the SpectroTyper software.
We excluded all the case or control samples with missing genotypes on more than 50 percent of the 67 studied SNPs. We also excluded SNPs with significant deviation from Hardy Weinberg equilibrium (P values less than 0.005) in the control group, and we did not consider SNPs with minor allele frequency (MAF) less than 0.01. Association tests for each SNP in cases and controls were determined using Haploview version 3.32.
Alleles/genotypes or haplotypes of STAT4 were analyzed for association with clinical or laboratory variables using chi-square test, Mann-Whitney U test, or logistic regression test, using SAS software, version 9.1 (SAS Institute, Cary, NC, USA). P < 0.05 is considered significant.
Association of a STAT4 Haplotype with Rheumatoid Arthritis
We analyzed 2131 samples (1123 cases and 1008 controls) out of an initial 2150 samples (1128 cases and 1022 controls) after excluding 19 samples (five cases and 14 controls) with more than 50 percent missing genotypes. In the remaining subjects, the average call rate in the 67 SNPs was 95.6 percent (84.3–99.7 percent). After filtering all SNPs with minor allele frequency (MAF) less than 0.01 (n = 12) or SNPs with significant deviation from Hardy Weinberg equilibrium, i.e. P values less than 0.005 in the control group (n = 3), we analyzed the distribution of allele frequency of 52 SNPs in the RA cases and controls. We removed 15 SNPs out of initial 67 studied SNPs based on Hardy Weinberg Equilibrium less than P value 0.005 (rs3024839, rs10931481, and rs7596818) or MAF less than 0.01 (rs17749316, rs13005843, rs3024933, rs16833220, rs932169, rs2459611, rs13010752, rs12998748, rs13011805, rs7599504, rs10194402, rs17769459). Call rates for each of the 52 SNPs is shown. (Table 1). The most significant four SNPs (P < 0.005) in the present study are identical with those in the North American study.
The LD pattern of the 52 SNPs within the region of STAT1/STAT4 is shown in Figure 1. The four most significant disease-associated SNPs (rs11889341, rs7574865, rs8179673, and rs10181656) are in strong linkage disequilibrium (D′ > 0.93; r2 > 0.73). A haplotype association analysis of these four risk SNPs revealed that the same haplotype (TTCG) as the one in the North American study showed significant risk for RA [0.34 versus 0.28, P = 0.0027, OR (95 percent CI) = 1.33 (1.10–1.60)] (Table 2).
Stratification Analysis by Clinical or Laboratory Variables
We examined patient subgroups for association with STAT4 using logistic regression. There was no significant genotypic or haplotypic association with sex, radiographic severity, or age of disease onset (data not shown). When we stratified the RA cases into anti-CCP positive and negative groups, we observed significant associations of all four risk SNPs and the risk haplotype of STAT4 with anti-CCP positive RA (Table 3), but not with anti-CCP negative RA compared with controls. However, because we observed the similar trend of distribution between both groups, this lack of significance in anti-CCP negative RA probably resulted from the smaller sample size (anti-CCP negative RA, N = 111). This finding suggests that the STAT4 contribution to disease risk may not be restricted to the anti-CCP + RA subset. Using a logistic regression analysis, the STAT4 risk haplotype is an independent risk factor (P = 0.049) in addition to SE (P < 0.0001) for RA.
To address the relationship of the disease-associated STAT4 variants with anti-CCP and RF titer, we divided the cases with anti-CCP into high (≥ 100) and low titer groups. This was necessary because CCP titers were truncated at the high end. This analysis did not show an effect of the STAT4 risk alleles or haplotype on the frequency of high anti-CCP group (data not shown). There was a suggestion of an effect on RF titer, but this was not significant, as shown in Table 4.
The North American Rheumatoid Arthritis Consortium (NARAC) recently reported on two new RA linkage regions at chromosomes 2q33 and 11p12 with logarithm of odds (LOD) scores of 3.52 and 3.09, respectively (20). Dense SNP mapping of the 2q RA linkage peak led to identification of a new susceptibility gene, STAT4, for RA (16). This association was replicated in several independent White RA case and control populations (16). In the current study, we have now confirmed that STAT4 is associated with RA in a large Korean population dataset, with the same common haplotype, which is more common in Koreans than in North American populations, but nevertheless confers a similar degree of risk.
A significant source of variability in the RA genetics literature has been the inability to replicate genetic findings across the major racial groups, particularly Whites and Asians. An interesting example of this is the association of the intracellular phosphatase, PTPN22, with RA and other autoimmune diseases. These disease associations have been widely replicated in White populations(7), but the PTPN22 risk allele (R620W) is exceedingly rare in Asian populations(15). Furthermore, attempts to identify other risk variants of PTPN22 that might be associated with RA in Asians have been unrevealing (9, H-S Lee unpublished). This has raised the possibility that there is true locus heterogeneity for RA among these major racial groups.
Several other examples have arisen in which associations are observed in Asian populations, but not in Whites. The most robust of these examples is the association of PADI4 with rheumatoid arthritis. PADI4 is a compelling candidate gene, because it encodes one of the enzymes responsible for citrullination of endogenous proteins, and an antibody response to citrullinated peptides is highly specific for RA(21). Numerous studies in Asian populations have demonstrated the association of PADI4 with RA (10,22–25), but these associations are either absent or very weak in populations of European ancestry (13,14). It is possible that this difference reflects an interaction of PADI4 genetic susceptibility with environmental factors, because citrullination may also be related to smoking or other environmental exposures (1,2,26).
The current report is the first clear demonstration of non-MHC related susceptibility gene for RA that confers a similar degree of risk among both White and Asian populations. Furthermore, it appears that the risk haplotype is likely to be identical in the two racial groups, suggesting that the responsible functional variant is ancient in origin. Indeed, the same haploytpe also is found in African populations (TTCG, 0.14 in the Yoruba people of Ibadan, Nigeria, www.hapmap.org) and it will be of great interest to see if the STAT4 associations with RA and lupus are also present in this population group. The associated haplotype is located primarily in the third intron of the STAT4 gene, and the actual functional allele(s) remain to be identified. A full resequencing of the STAT4 gene is in progress, and this will help to direct future studies of splice variation and/or expression differences that may explain the disease associated haplotype.
STAT4 encodes a transcription factor that lies in the signaling pathway of several important cytokines, including IL-12 and type I interferons, as well as IL-23 (27). STAT4 is present in the cytosol and upon cytokine signaling it becomes phosphoryated and translocates to the nucleus. The target genes for STAT4 include γIFN and therefore it plays a key role in the IL12 induced differentiation of T cells into the Th1 pathway. In addition, STAT4 may also be involved in the production of IL17 by Th17 cells, in response to IL23 (28). At the same time, γIFN production tends to inhibit differentiation toward the Th17 pathway(29). Thus, while genetic differences in STAT4 dependent signaling may be involved in regulating the balance of Th1 vs. Th17 responses, the expected effects of increased vs. decreased STAT4 activity are not obvious. Furthermore, compared with the mouse, the production of IL17 is not as clearly restricted to an easily definable Th17 subset in humans (30,31), and recent work suggests that some subsets of CD4 cells can produce both IL17 and γIFN (29). Therefore, it will be important to try to relate the major risk haplotype of STAT4 to phenotypic differences in these various T cell subsets.
Given the current evidence of Th17 involvement in chronic inflammation in RA (29), it would be expected that the STAT4 risk alleles would generally enhance a Th17 response. Several studies in animal models strongly suggest a key regulatory role of STAT4 in experimental arthritis (32). There is evidence that STAT4 may play a role both during initiation of disease as well as in the maintenance of the inflammatory process, leading to the idea that STAT4 may be a useful therapeutic target as has been demonstrated in murine collagen induced arthritis (32). In this context, it may be relevant that STAT4 also is involved in signaling responses to type 1 interferons in activated monocytes, macrophages, and mature dendritic cells. The dependence of type 1 interferon signaling on STAT4 is of particular interest in view of our recent observation that in addition to RA, STAT4 has a strong association with systemic lupus (16), a disease in which dysregulation of interferon pathways is prominent.
The finding of a common risk haplotype for RA among Asian and White populations shows that replication studies across racial boundaries may be useful for confirming at least some risk loci. In view of the large number of potential risk alleles that are coming out of whole genome scans (33), our data suggest that comparison across racial groups can be a reasonable approach to gene identification, and we hope to be able to carry out a genome-wide study of the Korean RA population in the near future.
Klareskog L, Padyukov L, Ronnelid J, Alfredsson L. (2006) Genes, environment and immunity in the development of rheumatoid arthritis. Curr. Opin. Immunol. 18:650–5.
Klareskog L et al. (2006) A new model for an etiology of rheumatoid arthritis: smoking may trigger HLA-DR (shared epitope)-restricted immune reactions to autoantigens modified by citrullination. Arthritis Rheum. 54:38–46.
Seldin MF, Amos CI, Ward R, Gregersen PK. (1999) The genetics revolution and the assault on rheumatoid arthritis. Arthritis Rheum. 42:1071–9.
Hall FC et al. (1996) Influence of the HLA-DRB1 locus on susceptibility and severity in rheumatoid arthritis. Qjm. 89:821–9.
Lee HS, Lee KW, Song GG, Kim HA, Kim SY, Bae SC. (2004) Increased susceptibility to rheumatoid arthritis in Koreans heterozygous for HLA-DRB1*0405 and *0901. Arthritis Rheum. 50: 3468–75.
Begovich AB et al. (2004) A missense single-nucleotide polymorphism in a gene encoding a protein tyrosine phosphatase (PTPN22) is associated with rheumatoid arthritis. Am. J. Hum. Genet. 75:330–7.
Gregersen PK, Lee HS, Batliwalla F, Begovich AB. (2006) PTPN22: setting thresholds for autoimmunity. Semin. Immunol. 18:214–23.
Ikari K et al. (2006) Haplotype analysis revealed no association between the PTPN22 gene and RA in a Japanese population. Rheumatology (Oxford) 45:1345–8.
Kawasaki E et al. (2006) Systematic search for single nucleotide polymorphisms in a lymphoid tyrosine phosphatase gene (PTPN22): association between a promoter polymorphism and type 1 diabetes in Asian populations. Am. J. Med. Genet. A. 140:586–93.
Suzuki A et al. (2003) Functional haplotypes of PADI4, encoding citrullinating enzyme peptidy-larginine deiminase 4, are associated with rheumatoid arthritis. Nat. Genet. 34:395–402.
Tokuhiro S et al. (2003) An intronic SNP in a RUNX1 binding site of SLC22A4, encoding an organic cation transporter, is associated with rheumatoid arthritis. Nat. Genet. 35:341–8.
Kochi Y et al. (2005) A functional variant in FCRL3, encoding Fc receptor-like 3, is associated with rheumatoid arthritis and several autoimmunities. Nat. Genet. 37:478–85.
Plenge RM et al. (2005) Replication of putative candidate-gene associations with rheumatoid arthritis in > 4,000 samples from North America and Sweden: association of susceptibility with PTPN22, CTLA4, and PADI4. Am. J. Hum. Genet. 77:1044–60.
Barton A et al. (2004) A functional haplotype of the PADI4 gene associated with rheumatoid arthritis in a Japanese population is not associated in a United Kingdom population. Arthritis Rheum. 50:1117–21.
Mori M, Yamada R, Kobayashi K, Kawaida R, Yamamoto K. (2005) Ethnic differences in allele frequency of autoimmune-disease-associated SNPs. J. Hum. Genet. 50:264–6.
Remmers EF et al. (2007) STAT4 and risk of rheumatoid arthritis and systemic lupus erythematosus. N. Engl. J. Med. 357:13–22.
Arnett FC et al. (1988) The American Rheumatism Association 1987 revised criteria for the classification of rheumatoid arthritis. Arthritis Rheum. 31:315–24.
Begnon JD, Fernandez-Vina MA. (1997) Protocols of the 12th International Histocompatibility Workshop for typing of HLA class II alleles by DNA amplification by the polymerase chain reaction (PCR) and hybridization with sequence specific oligonucleotide probes (SSOP). In: Charron D (ed.) Genetic Diversity of HLA: Functional and Medical Implication. EDK, Paris, p 584–95.
Steinbrocker O, Traeger CH, Batterman RC. (1949) Therapeutic criteria in rheumatoid arthritis. JAMA 140: 659–62.
Amos CI et al. (2006) High-density SNP analysis of 642 White families with rheumatoid arthritis identifies two new linkage regions on 11p12 and 2q33. Genes Immun. 7:277–86.
Schellekens GA et al. (2000) The diagnostic properties of rheumatoid arthritis antibodies recognizing a cyclic citrullinated peptide. Arthritis Rheum. 43:155–63.
Ikari K et al. (2005) Association between PADI4 and rheumatoid arthritis: a replication study. Arthritis Rheum. 52: 3054–7.
Kang CP, Lee HS, Ju H, Cho H, Kang C, Bae SC. (2006) A functional haplotype of the PADI4 gene associated with increased rheumatoid arthritis susceptibility in Koreans. Arthritis Rheum. 54: 90–6.
Cha S, Choi CB, Han TU, Kang CP, Kang C, Bae SC. (2007) Association of anti-cyclic citrullinated peptide antibody levels with PADI4 haplotypes in early rheumatoid arthritis and with shared epitope alleles in very late rheumatoid arthritis. Arthritis Rheum. 56:1454–63.
Lee YH, Rho YH, Choi SJ, Ji JD, Song GG. (2007) PADI4 polymorphisms and rheumatoid arthritis susceptibility: a meta-analysis. Rheumatol. Int. 27: 827–33.
Lee HS et al. (2007) Interaction between smoking, the shared epitope, and anti-cyclic citrullinated peptide: a mixed picture in three large North American rheumatoid arthritis cohorts. Arthritis Rheum. 56:1745–53.
Murphy KM, Reiner SL. (2002) The lineage decisions of helper T cells. Nat. Rev. Immunol. 2: 933–44.
Mathur AN et al. (2007) Stat3 and Stat4 direct development of IL-17-secreting Th cells. J. Immunol. 178:4901–7.
Miossec P. (2007) Interleukin-17 in fashion, at last: Ten years after its description, its cellular source has been identified. Arthritis Rheum. 56:2111–5.
Acosta-Rodriguez EV et al. (2007) Surface phenotype and antigenic specificity of human interleukin 17-producing T helper memory cells. Nat. Immunol. 8:639–46.
Chen Z, Tato CM, Muul L, Laurence A, O’Shea JJ. (2007) Distinct regulation of IL-17 in human helper T lymphocytes. Arthritis. Rheum. 56:2936–46.
Hildner KM et al. (2007) Targeting of the transcription factor STAT4 by antisense phosphorothioate oligonucleotides suppresses collagen-induced arthritis. J. Immunol. 178:3427–36.
(2007) Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls. Nature. 447:661–8.
This work was made possible by support from the Eileen Ludwig Greenland Center for Rheumatoid Arthritis. Additional support was provided by grants from the Korea Health 21 R&D Project, Ministry of Health and Welfare, Republic of Korea (01-PJ3-PG6-01GN11-0002), the National Institutes of Health, RO1-AR44222 (PKG) and the intramural program of the National Institutes of Arthritis and Musculoskeletal and Skin Diseases.
About this article
Cite this article
Lee, H., Remmers, E.F., Le, J.M. et al. Association of STAT4 with Rheumatoid Arthritis in the Korean Population. Mol Med 13, 455–460 (2007). https://0-doi-org.brum.beds.ac.uk/10.2119/2007-00072.Lee