ETS1 variants confer susceptibility to ankylosing spondylitis in Han Chinese
© Shan et al.; licensee BioMed Central Ltd. 2014
Received: 24 January 2014
Accepted: 25 March 2014
Published: 4 April 2014
ETS1 is a negative regulator of the Th17 differentiation gene and plays a central role in the pathogenesis of autoimmune diseases. We aimed to investigate whether polymorphisms in ETS1 confer susceptibility to ankylosing spondylitis (AS) in Han Chinese.
We selected seven single nucleotide polymorphisms (SNPs) within ETS1 based on HapMap data and previous genome-wide association study. Genotyping involved the TaqMan method in 1,015 patients with AS and 1,132 healthy controls from Shandong Province, and 352 AS patients and 400 healthy controls from Ningxia, a northwest region in China. Gene expression was determined by real-time PCR.
The SNP rs1128334 was strongly associated with AS (odds ratio 1.204, 95% confidence interval 1.06-1.37; P = 0.005). This association was confiexrmed in the Ningxia population (P = 0.015). Carriers of the haplotype TAT for rs12574073, rs1128334 and rs4937333 were associated with increased risk of AS and haplotype CGC with reduced risk as compared to controls. In addition, ETS1 expression was lower in AS patients than controls. The risk allele A of rs1128334 and haplotype A-T of rs1128334 and rs4937333 were associated with decreased expression of ETS1.
Common variants in ETS1 may contribute to AS susceptibility in Han Chinese people.
Ankylosing spondylitis (AS) is a common inflammatory arthritis that affects the axial skeleton and peripheral joints ; it is usually accompanied by lower back pain, peripheral arthritis, enthesis and iritis, and even spinal deformity and ankylosis, ultimately limiting the mobility of the spine and other joints . Genetic factors are strongly implicated in the pathogenesis of the disease, and the estimated heritability is up to 97% according to a twin study [3, 4].
The major histocompatibility complex (MHC), mostly from human leukocyte antigen B27 (HLA-B27), accounts for nearly half of the predisposition for AS . Although HLA-B27 is strongly associated with risk of AS, its associated genes account for only 20% to 30% of the overall genetic risk of AS . Thus, other loci may be involved. Besides HLA-B27, non-MHC genes implicated in the risk of AS include ERAP1, KIF21B, IL23R and two intergenic regions at chromosome 2p15 and chromosome 21q22 in a European population [6, 7].
Two recent genome-wide association studies (GWASs) of an Asian population demonstrated a novel region containing ETS1 at chromosome 11q23 strongly associated with systemic lupus erythematosus (SLE) [8, 9]. rs6590330, rs1128334 and rs4937333 were identified as the susceptible variants associated with SLE. Subsequently, many genetic studies reported that rs11221332 of ETS1 is associated with rheumatoid arthritis (RA) and celiac disease, both characterized by excessive activation of the immune system [10, 11]. Also, increasing evidence has suggested that many of the known AS-associated loci overlap with those of other immune-related diseases. Because of clinical and immunological overlap of AS and other autoimmune diseases , ETS1 may be associated with AS as well. In this study, we investigated the association of ETS1 genetic polymorphisms and AS in Han Chinese people.
We recruited 2,899 unrelated subjects in this case-control study from 2008 to 2011, including 1,015 patients (701 male, 69.1%) with AS from Shandong Provincial Hospital and Qilu Hospital in Jinan, Shandong province. All patients, at a mean age of 38 ± 8.3 years, had no history of inflammatory bowel disease and/or psoriasis and met the modified New York criteria for AS. We also recruited 1,132 unrelated, randomly sampled healthy controls (603 male, 53.3%) from a health check-up center at the hospital during the same period; the mean age was 40 ± 7.2 years old. All subjects were ethnic Han living in Shandong, an eastern coastal province of North China. More detailed information about the samples was reported previously . The replication samples of 352 AS patients and 400 controls from Ningxia (a province in northwest China) were provided by Professor Yang Ze, with detailed information . Flow cytometry was used to immunophenotype HLA-B27 in patients and controls, and only HLA-B27-positive patients and HLA-B27-negative controls were recruited.
Genomic DNA from peripheral blood leukocytes was extracted by a standard phenol-chloroform method. The study was approved by the ethics committee for human studies at the Shandong University School of Medicine, and all subjects provided their informed consent.
Sequence analysis of the human ETS1 gene
All nine coding exons and 1,000 bp of the 5′ upstream sequence of ETS1 were sequenced in DNA samples from 100 AS cases to screen for genetic variants. ETS1 was amplified by PCR before sequencing. The PCR reaction was run in a total volume of 50 μL solution containing 5 × PCR buffer, 10 μL; 5 × dNTP (1 mmol/L each of dATP, dTTP, dGTP and dCTP), 10 μL; 50 ng genomic DNA; 0.4 μM of each primer; and 1 U Taq DNA polymerase. The PCR amplification conditions were optimized with an initial denaturation step at 94°C for 5 minutes followed by 35 cycles of 94°C for 40 sec, optimal annealing temperature for 40 sec, extension at 72°C for 40 sec and a final extension step of 72°C for 10 minutes. All sequencing was performed by the Genomic Analysis Facility (Shenzhen Huada, Shenzhen, China) and analyzed with use of Chromas 2.33 software.
Single nucleotide polymorphism (SNP) selection and genotyping
Total peripheral leukocyte RNA was extracted from 50 AS patients and 161 healthy controls by the TRIzol reagent method (Invitrogen, Carlsbad, CA, USA), then underwent reverse transcription by use of a reverse transcript kit (Takara, Otsu, Shiga, Japan). Allelic-specific expression of ETS1 was analyzed only in healthy controls. RNA samples were treated with RNase free DNAase to eliminate genomic DNA contamination before being reverse-transcribed into cDNA. Real-time quantitative RT-PCR to amplify cDNA involved the Roche 480 real-time PCR system with SYBR Green Master Mix (Applied Biosystems, Foster City, CA, USA). Relative expression was analyzed by the comparative threshold cycle (Ct) method and normalized to that of human glyceraldehyde-3-phosphate dehydrogenase (GAPDH), calculated by the 2-△△Ct method and log10 transformed. The primer sequences for ETS1 were sense 5′-TACACAGGCAGTGGACCAATC-3′ and antisense 5′-CCCCGCTGTCTTGTGGATG-3′; and GAPDH, sense 5′-CCAGGTGGTCTCCTCTGACTT-3′ and antisense 5′-GTTGCTGTAGCCAAATTCGTTGT-3′. Melting curve analysis was used to confirm specificity, with four duplicate wells used for each subject .
The SNPs were tested for adherence to Hardy-Weinberg equilibrium (HWE) by chi-square test. Comparison of the allele frequency of the seven SNPs for patients and controls involved Fisher’s exact test with use of Plink 1.07 . The Bonferroni correction was used for multiple test correction. Pairwise linkage disequilibrium (LD) was measured with use of Haploview 4.2. Three-locus haplotypes were calculated by use of a full precise-iteration algorithm with the online SHEsis software . Differences in allelic expression were tested by t- test and analysis of variance (ANOVA) by use of GraphPad Prism v5.01 (GraphPad Software, La Jolla, CA, USA). P <0.05 was considered statistically significant.
To search for new variants in the ETS1 gene, we sequenced all exons including the untranslated region (UTR) and 5′ upstream regulatory region. However, no novel variant was found in any of the nine exons or 1,000 bp of the upstream sequence.
Genotype and allele association analysis of seven single nucleotide polymorphisms (SNPs) in samples of AS in Shandong and rs1128334 in Ningxia
Patients, number (%)
Controls, number (%)
LD evaluation and haplotype analysis
Haplotype frequencies of three single nucleotide polymorphisms in all samples
No. (%) of patients
No. (%) of controls
Odds ratio (95% CI)
T A T
1.24 (1.05, 1.47)
C G C
0.82 (0.70, 0.96)
C G T
1.01 (0.76, 1.35)
mRNA expression of ETS1 in peripheral blood mononuclear cells (PBMCs)
Allelic-specific expression of ETS1
We have demonstrated that ETS1 is associated with AS, thus adding to the list of loci showing overlap between AS, RA and SLE. We identified an SNP, rs1128334, located in the 3′ UTR of ETS1, which was significantly associated with AS in Han Chinese people. We also showed a lower expression of ETS1 in peripheral leukocytes from patients than controls. Furthermore, the risk allele of rs1128334 and the risk haplotype A-T of rs1128334 and rs4937333 were significantly linked to decreased mRNA level of ETS1. To our knowledge, this is the first report demonstrating an association of ETS1 polymorphisms and AS in Han Chinese. Especially, we verified the association of rs1128334 and AS in two populations, Han Chinese from Shandong and Ningxia, which enhances the credibility of our findings.
In previous reports, the common variants rs1128334, rs4937333 and rs6590330 were found to be significantly associated with SLE in east Asian populations. In addition, rs11221332 was suggested to be strongly associated with celiac disease in a European population and RA in a Caucasian population. In agreement with these studies of autoimmune diseases, we observed that ETS1 was also a susceptibility gene for AS at least in northern Han Chinese. Our findings provide support for many autoimmune diseases sharing the same susceptibility loci.
ETS1 encodes a member of the ETS family of transcription factors that activate transcription by binding to cis-regulatory elements in target genes [19, 20]. The ETS1 transcription factor was initially discovered as the proto-oncogene corresponding to v-ets of the avian erythroblastosis virus (E26) [19, 21], which harbors a conserved DNA-binding domain mediating specific DNA binding to the GGAA/T motif . As a crucial transcription factor widely expressed in lymphocytes and vascular endothelial, lacteal glandular epithelium and various invasion tumor cells, ETS1 regulates the development, senescence and death of many immune cells [22, 23]. It also plays a role in both innate and adaptive immune response [23, 24]. Accumulating evidence points to an important role for ETS1 in regulating the differentiation of immune cells such as T-cell differentiation into a helper population, terminal differentiation of B cells, development of natural killer (NK) cells and NK T cells and the expression of cytokine and chemokine genes in a wide variety of different cell lineages [24–27]. Animal experiments showed that autoimmune disease developed in ETS1-knockout mice, as investigated by the production of high titers of autoantibodies, and immune cell infiltration into organs accounted for aberrations in lymphocyte differentiation .
In addition to roles in controlling the production of interferon-gamma in T helper (Th)1 cells  and driving normal Th2 cytokine production, ETS1 is also involved in the Th17 cell lineage as a negative regulator of Th17 cell differentiation . Increasing studies suggest that Th17 cells are involved in the pathogenesis of AS . Recently, Zhang et al. examined serum IL-17 levels from 283 SLE cases and observed a significant synergistic epistatic interaction between two risk variants, rs10893872 and rs1128334, in the ETS1 gene, and the haplotype formed by these variants was significantly associated with serum IL-17 levels in SLE patients , which supports a role for genetic interaction contributing to the complexity of autoimmune disease [31–33]. From our findings, subjects with AS carrying the risk allele A have significantly lower ETS1 expression than those with the G allele. Similar to previous studies, the haplotype formed by rs12574073, rs1128334 and rs4937333 in ETS1 was significantly associated with AS susceptibility. As well, individuals carrying haplotype AT for rs1128334/rs4937333 showed less ETS1 expression than those carrying the GC haplotype. With the HCB data in HapMap, rs10893872 and rs4937333 are in perfect LD (r2 = 1), so our results are consistent with the Zhang et al. study.
Much evidence has confirmed that microRNAs (miRNAs) can regulate gene expression by sequence-specific binding to target mRNAs, but the binding affinity can be affected by SNPs residing in miRNA target sites, which may in turn affect the ability of miRNA to inhibit the mRNA translation into proteins or lead to degradation of the mRNA . Indeed, recent research has demonstrated that SNPs at miRNA binding sites likely affect the expression of the miRNA target genes and thus, may contribute to susceptibility to autoimmune diseases . Therefore, we presumed that variants rs1128334 and rs4937333 located in the 3′ UTR of ETS1 may affect miRNA binding affinity and further influence the expression of the target gene, which may contribute to susceptibility to AS. However, the detailed mechanism of this functional relationship requires further investigation.
In summary, our study identified that the common variants of ETS1 confer susceptibility to AS. This study might provide further evidence to improve our understanding of the exact function of ETS1 in the pathogenesis of autoimmune diseases.
analysis of variance
Endoplasmic reticulum aminopeptidase 1
v-ets avian erythroblastosis virus E26 oncogene homolog 1
genome-wide association study
Han Chinese in Beijing
human leukocyte antigen
interleukin 23 receptor
Kinesin family member 21B
major histocompatibility complex
- NK cell:
natural killer cell
peripheral blood mononuclear cell
reverse transcription PCR
systemic lupus erythematosus
- Th cell:
T helper cell
We thank all the participants in the study. This work was supported by grants from the Natural Science Foundation of China (81072452, 81273281) and the Ph.D. Programs Foundation of the Ministry of Education of China (20100131110035).
- Braun J, Sieper J: Ankylosing spondylitis. Lancet. 2007, 369: 1379-1390. 10.1016/S0140-6736(07)60635-7.View ArticlePubMedGoogle Scholar
- van der Heijde D, Calin A, Dougados M, Khan MA, van der Linden S, Bellamy N: Selection of instruments in the core set for DC-ART, SMARD, physical therapy, and clinical record keeping in ankylosing spondylitis, Progress report of the ASAS Working Group. Assessments in Ankylosing Spondylitis. J Rheumatol. 1999, 26: 951-954.PubMedGoogle Scholar
- Braun J, Bollow M, Remlinger G, Eggens U, Rudwaleit M, Distler A, Sieper J: Prevalence of spondylarthropathies in HLA-B27 positive and negative blood donors. Arthritis Rheum. 1998, 41: 58-67. 10.1002/1529-0131(199801)41:1<58::AID-ART8>3.0.CO;2-G.View ArticlePubMedGoogle Scholar
- Brown MA, Kennedy LG, MacGregor AJ, Darke C, Duncan E, Shatford JL, Taylor A, Calin A, Wordsworth P: Susceptibility to ankylosing spondylitis in twins: the role of genes, HLA, and the environment. Arthritis Rheum. 1997, 40: 1823-1828. 10.1002/art.1780401015.View ArticlePubMedGoogle Scholar
- Reveille JD: The genetic basis of spondyloarthritis. Ann Rheum Dis. 2011, 70: 44-50.View ArticleGoogle Scholar
- Danoy P, Pryce K, Hadler J, Bradbury LA, Farrar C, Pointon J, Ward M, Weisman M, Reveille JD, Wordsworth BP, Stone MA, Maksymowych WP, Rahman P, Gladman D, Inman RD, Brown MA: Association of variants at 1q32 and STAT3 with ankylosing spondylitis suggests genetic overlap with Crohn’s disease. PLoS Genet. 2010, 6: e1001195-10.1371/journal.pgen.1001195.PubMed CentralView ArticlePubMedGoogle Scholar
- Burton PR, Clayton DG, Cardon LR, Craddock N, Deloukas P, Duncanson A, Kwiatkowski DP, McCarthy MI, Ouwehand WH, Samani NJ, Todd JA, Donnelly P, Barrett JC, Davison D, Easton D, Evans DM, Leung HT, Marchini JL, Morris AP, Spencer CC, Tobin MD, Attwood AP, Boorman JP, Cant B, Everson U, Hussey JM, Jolley JD, Knight AS, Koch K, Meech E: Association scan of 14,500 nonsynonymous SNPs in four diseases identifies autoimmunity variants. Nat Genet. 2007, 39: 1329-1337. 10.1038/ng.2007.17.View ArticlePubMedGoogle Scholar
- Yang W, Shen N, Ye DQ, Liu Q, Zhang Y, Qian XX, Hirankarn N, Ying D, Pan HF, Mok CC, Chan TM, Wong RW, Lee KW, Mok MY, Wong SN, Leung AM, Li XP, Avihingsanon Y, Wong CM, Lee TL, Ho MH, Lee PP, Chang YK, Li PH, Li RJ, Zhang L, Wong WH, Ng IO, Lau CS, Sham PC: Genome-wide association study in Asian populations identifies variants in ETS1 and WDFY4 associated with systemic lupus erythematosus. PLoS Genet. 2010, 6: e1000841-10.1371/journal.pgen.1000841.PubMed CentralView ArticlePubMedGoogle Scholar
- Han JW, Zheng HF, Cui Y, Sun LD, Ye DQ, Hu Z, Xu JH, Cai ZM, Huang W, Zhao GP, Xie HF, Fang H, Lu QJ, Li XP, Pan YF, Deng DQ, Zeng FQ, Ye ZZ, Zhang XY, Wang QW, Hao F, Ma L, Zuo XB, Zhou FS, Du WH, Cheng YL, Yang JQ, Shen SK, Li J, Sheng YJ: Genome-wide association study in a Chinese Han population identifies nine new susceptibility loci for systemic lupus erythematosus. Nat Genet. 2009, 41: 1234-1237. 10.1038/ng.472.View ArticlePubMedGoogle Scholar
- Chatzikyriakidou A, Voulgari PV, Georgiou I, Drosos AA: Altered sequence of the ETS1 transcription factor may predispose to rheumatoid arthritis susceptibility. Scand J Rheumatol. 2013, 42: 11-14. 10.3109/03009742.2012.711367.View ArticlePubMedGoogle Scholar
- Dubois PC, Trynka G, Franke L, Hunt KA, Romanos J, Curtotti A, Zhernakova A, Heap GA, Adany R, Aromaa A, Bardella MT, vanden Berg LH, Bockett NA, dela Concha EG, Dema B, Fehrmann RS, Fernandez-Arquero M, Fiatal S, Grandone E, Green PM, Groen HJ, Gwillian R, Houwen RH, Hunt SE, Kaukinen K, Kelleher D, Korponay-Szabo I, Kurppa K, MacMathuna P, Maki M: Multiple common variants for celiac disease influencing immune gene expression. Nat Genet. 2010, 42: 295-302. 10.1038/ng.543.PubMed CentralView ArticlePubMedGoogle Scholar
- Richard-Miceli C, Criswell LA: Emerging patterns of genetic overlap across autoimmune disorders. Genome Med. 2012, 4: 6-10.1186/gm305.PubMed CentralView ArticlePubMedGoogle Scholar
- Liu Y, Zhang H, Li J, Zhao H, Xin Q, Shan S, Dang J, Bian X, Liu Q: Association of common variants in KIF21B and ankylosing spondylitis in a Chinese Han population: a replication study. Immunogenetics. 2013, 65: 835-839. 10.1007/s00251-013-0733-6.View ArticlePubMedGoogle Scholar
- Chen J, Zhou L, Huo ZH, Zhang YH, Yang ZH, Yang BZ, Huang CB, Zhu XQ, Yang Z: Identification of a novel lymphotoxin-alpha (LTA) gene associated with ankylosing spondylitis in Ningxia population. Yi Chuan. 2011, 33: 329-336. 10.3724/SP.J.1005.2011.00329.View ArticlePubMedGoogle Scholar
- The data about the minor allele frequency of all SNPs for Han Chinese in Beijing (HCB) in HapMap. [http://www.hapmap.org]
- Zhao H, Yang W, Qiu R, Li J, Xin Q, Wang X, Feng Y, Shan S, Liu Y, Gong Y, Liu Q: An intronic variant associated with systemic lupus erythematosus changes the binding affinity of Yinyang1 to downregulate WDFY4. Genes Immun. 2012, 13: 536-542. 10.1038/gene.2012.33.View ArticlePubMedGoogle Scholar
- Comparison of the allele frequency of the 7 SNPs for patients and controls involved Fisher’s exact test using Plink 1.07 online. [http://pngu.mgh.harvard.edu/~purcell/plink/]
- Shi YY, He L: SHEsis, a powerful software platform for analyses of linkage disequilibrium, haplotype construction, and genetic association at polymorphism loci. Cell Res. 2005, 15: 97-98. 10.1038/sj.cr.7290272.View ArticlePubMedGoogle Scholar
- Leprince D, Gegonne A, Coll J, de Taisne C, Schneeberger A, Lagrou C, Stehelin D: A putative second cell-derived oncogene of the avian leukaemia retrovirus E26. Nature. 1983, 306: 395-397. 10.1038/306395a0.View ArticlePubMedGoogle Scholar
- Wang CY, Petryniak B, Ho IC, Thompson CB, Leiden JM: Evolutionarily conserved Ets family members display distinct DNA binding specificities. J Exp Med. 1992, 175: 1391-1399. 10.1084/jem.175.5.1391.View ArticlePubMedGoogle Scholar
- Nunn MF, Seeburg PH, Moscovici C, Duesberg PH: Tripartite structure of the avian erythroblastosis virus E26 transforming gene. Nature. 1983, 306: 391-395. 10.1038/306391a0.View ArticlePubMedGoogle Scholar
- Du J, Shen L, Zhao GH, Wang YG, Liao SS, Chen C, Zhou ZF, Luo YY, Jiang H, Xia K, Tang BS: Receptor expression-enhancing protein 1 gene (SPG31) mutations are rare in Chinese Han patients with hereditary spastic paraplegia. Chin Med J (Engl). 2009, 122: 2064-2066.Google Scholar
- Gallant S, Gilkeson G: ETS transcription factors and regulation of immunity. Arch Immunol Ther Exp (Warsz). 2006, 54: 149-163. 10.1007/s00005-006-0017-z.View ArticleGoogle Scholar
- Russell L, Garrett-Sinha LA: Transcription factor Ets-1 in cytokine and chemokine gene regulation. Cytokine. 2010, 51: 217-226. 10.1016/j.cyto.2010.03.006.View ArticlePubMedGoogle Scholar
- Maier H, Colbert J, Fitzsimmons D, Clark DR, Hagman J: Activation of the early B-cell-specific mb-1 (Ig-alpha) gene by Pax-5 is dependent on an unmethylated Ets binding site. Mol Cell Biol. 2003, 23: 1946-1960. 10.1128/MCB.23.6.1946-1960.2003.PubMed CentralView ArticlePubMedGoogle Scholar
- Garvie CW, Hagman J, Wolberger C: Structural studies of Ets-1/Pax5 complex formation on DNA. Mol Cell. 2001, 8: 1267-1276. 10.1016/S1097-2765(01)00410-5.View ArticlePubMedGoogle Scholar
- Pufall MA, Graves BJ: Ets-1 flips for new partner Pax-5. Structure. 2002, 10: 11-14. 10.1016/S0969-2126(01)00701-8.View ArticlePubMedGoogle Scholar
- Nagaleekar VK, Diehl SA, Juncadella I, Charland C, Muthusamy N, Eaton S, Haynes L, Garrett-Sinha LA, Anguita J, Rincon M: IP3 receptor-mediated Ca2+ release in naive CD4 T cells dictates their cytokine program. J Immunol. 2008, 181: 8315-8322.PubMed CentralView ArticlePubMedGoogle Scholar
- Moisan J, Grenningloh R, Bettelli E, Oukka M, Ho IC: Ets-1 is a negative regulator of Th17 differentiation. J Exp Med. 2007, 204: 2825-2835. 10.1084/jem.20070994.PubMed CentralView ArticlePubMedGoogle Scholar
- Shen H, Goodall JC, Hill Gaston JS: Frequency and phenotype of peripheral blood Th17 cells in ankylosing spondylitis and rheumatoid arthritis. Arthritis Rheum. 2009, 60: 1647-1656. 10.1002/art.24568.View ArticlePubMedGoogle Scholar
- Zhang J, Zhang Y, Zhang L, Yang J, Ying D, Zeng S, Lee TL, Lau CS, Chan TM, Leung AM, Mok CC, Wong SN, Lee KW, Ho MH, Lee PP, Chung BH, Chong CY, Wong RW, Mok MY, Wong WH, Lau YL, Yang W: Epistatic interaction between genetic variants in susceptibility gene ETS1 correlates with IL-17 levels in SLE patients. Ann Hum Genet. 2013, 77: 344-350.PubMedGoogle Scholar
- Hughes T, Adler A, Kelly JA, Kaufman KM, Williams AH, Langefeld CD, Brown EE, Alarcon GS, Kimberly RP, Edberg JC, Ramsey-Goldman R, Petri M, Boackle SA, Stevens AM, Reveille JD, Sanchez E, Martin J, Niewold TB, Vila LM, Scofield RH, Gilkeson GS, Gaffney PM, Criswell LA, Moser KL, Merrill JT, Jacob CO, Tsao BP, James JA, Vyse TJ, Alarcon-Riquelme ME: Evidence for gene-gene epistatic interactions among susceptibility loci for systemic lupus erythematosus. Arthritis Rheum. 2011, 64: 485-492.View ArticleGoogle Scholar
- Zhou XJ, Lu XL, Nath SK, Lv JC, Zhu SN, Yang HZ, Qin LX, Zhao MH, Su Y, Shen N, Li ZG, Zhang H: Gene-gene interaction of BLK, TNFSF4, TRAF1, TNFAIP3, and REL in systemic lupus erythematosus. Arthritis Rheum. 2011, 64: 222-231.View ArticleGoogle Scholar
- Wei JC, Yen JH, Juo SH, Chen WC, Wang YS, Chiu YC, Hsieh TJ, Guo YC, Huang CH, Wong RH, Wang HP, Tsai KL, Wu YC, Chang HW, Hsi E, Chang WP, Chang WC: Association of ORAI1 haplotypes with the risk of HLA-B27 positive ankylosing spondylitis. PLoS One. 2011, 6: e20426-10.1371/journal.pone.0020426.PubMed CentralView ArticlePubMedGoogle Scholar
- Hikami K, Kawasaki A, Ito I, Koga M, Ito S, Hayashi T, Matsumoto I, Tsutsumi A, Kusaoi M, Takasaki Y, Hashimoto H, Arinami T, Sumida T, Tsuchiya N: Association of a functional polymorphism in the 3′-untranslated region of SPI1 with systemic lupus erythematosus. Arthritis Rheum. 2011, 63: 755-763. 10.1002/art.30188.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.