DNA microarray analysis of rheumatoid arthritis susceptibility genes identified by genome-wide association studies
Arthritis Research & Therapy volume 12, Article number: 401 (2010)
In a recent interesting review, Alex Clarke and Timothy Vyse described the genetics of rheumatic disease . In the past several years, genome-wide association studies (GWAS) have led to the identification of six high-risk rheumatoid arthritis (RA) susceptibility genes - namely, CD244, PADI4, SLC22A2, PTPN22, CTLA4, and STAT4 (summarized in ). In vitro studies using mutant alleles and cultured cells have revealed the individual upregulation of CD244, PADI4, SLC22A2, and PTPN22 [2–6]; however, studies on the expression of RA susceptibility genes in RA patients are rare. We therefore investigated the expression of the above-mentioned six RA susceptibility genes in 112 RA patients using DNA microarray analysis. This study aims to clarify whether DNA microarray analysis and GWAS produce comparable results with respect to RA susceptibility genes.
Total RNA extracted from total peripheral blood cells obtained from 112 RA patients and 45 healthy individuals was used to prepare aminoallyl RNA. As a reference, mixed RNA from 45 healthy individuals was used. The aminoallyl RNA of each individual and the reference was subjected to Cy3 and Cy5 labeling, respectively, and was hybridized with an oligonucleotide-based DNA microarray. The data obtained were analyzed by nonparametric statistical group comparison. The intensities of the noprobe spots were used as the background. The median and standard deviation of the background intensity were calculated. The genes with an intensity value that was less than the median plus 2 standard deviation of the background intensity were identified as null. The Cy3/Cy5 ratios of all spots on the DNA microarray were normalized using the global ratio median. Only gene expression data that were collected from at least 80% of samples from each group were selected for further analysis. The unpaired Mann-Whitney test was used to determine statistically significant differences in the mRNA expression levels between the RA and healthy groups. Statistical significance was set at P < 0.05.
The results of our DNA microarray analysis showed that the expressions of four out of the six RA susceptibility genes were significantly higher in RA patients than in healthy individuals (1.0 × 10-16 to 2.32 × 10-5) (Table 1). As described above, the upregulation of these four genes (CD244, PADI4, SLC22A2, and PTPN22) has been previously confirmed in in vitro studies. We found, however, that CTLA4 expression levels were similar between the RA and control groups, whereas STAT4 expression was significantly downregulated in the RA group (1.38 × 10-8). We investigated the expression of other RA susceptibility genes - namely, TRF1/C5 , CD40 , and CCL21  - and found that their expressions were similar in both groups. The genetic risk factors for RA were recently reported to differ between Caucasian and Asian (Korean) populations . The samples used in our microarray analysis were derived from the same Asian (Japanese) cohort. The expression profiles for these three genes may therefore not be consistent with the profiles determined by GWAS.
In this study, we revealed the correlation between five out of the six high-risk RA susceptibility genes using DNA microarray analysis. Prostate cancer susceptibility genes identified by GWAS were recently reported to be consistent with those identified by microarray analysis . We therefore concluded that the combination of microarray analysis and GWAS would be a more effective approach for gene identification than the analysis of individual datasets. Moreover, the simultaneous use of both methods would allow for more accurate identification of RA candidate genes.
genome-wide association studies
Clarke A, Vyse TJ: Genetics of rheumatic disease. Arthritis Res Ther. 2009, 11: 248-10.1186/ar2781.
Suzuki A, Yamada R, Kochi Y, Sawada T, Okada Y, Matsuda K, Kamatani Y, Mori M, Shimane K, Hirabayashi Y, Takahashi A, Tsunoda T, Miyatake A, Kubo M, Kamatani N, Nakamura M, Yamamoto K: Functional SNPs in CD244 increase the risk of rheumatoid arthritis in Japanese population. Nat Genet. 2008, 40: 1224-1229. 10.1038/ng.205.
Suzuki A, Yamada R, Chang X, Tokuhiro S, Sawada T, Suzuki M, Nagasaki M, Nakayama-Hamada M, Kawaida R, Ono M, Ohtsuki M, Furukawa H, Yoshino S, Yukioka M, Tohma S, Matsubara T, Wakitani S, Teshima R, Nishioka Y, Sekine A, Iida A, Takahashi A, Tsunoda T, Nakamura Y, Yamamoto K: Functional haplotypes of PADI4, encoding citrullinating enzyme peptidylarginine deiminase 4, are associated with rheumatoid arthritis. Nat Genet. 2003, 34: 395-402. 10.1038/ng1206.
Chang X, Zhao Y, Sun S, Zhang Y, Zhu Y: The expression of PADI4 in synovium of rheumatoid arthritis. Rheumatol Int. 2009, 29: 1411-1416. 10.1007/s00296-009-0870-2.
Tokuhiro S, Yamada R, Chang X, Suzuki A, Kochi Y, Sawada T, Suzuki M, Nagasaki M, Ohtsuki M, Ono M, Furukawa H, Nagashima M, Yoshino S, Mabuchi A, Sekine A, Saito S, Takahashi A, Tsunoda T, Nakamura Y, Yamamoto K: An intronic SNP in a RUNX1 binding site of SLC22A4, encoding an organic cation transporter, is associated with rheumatoid arthritis. Nat Genet. 2003, 35: 341-348. 10.1038/ng1267.
Begovich AB, Carlton VE, Honigberg LA, Schrodi SJ, Chokkalingam AP, Alexander HC, Ardlie KG, Huang Q, Smith AM, Spoerke JM, Conn MT, Chang M, Chang SY, Saiki RK, Catanese JJ, Leong DU, Garcia VE, McAllister LB, Jeffery DA, Lee AT, Batliwalla F, Remmers E, Criswell LA, Seldin MF, Kastner DL, Amos CI, Sninsky JJ, Gregersen PK: A missense single-nucleotide polymorphism in a gene encoding a protein tyrosine phosphatase (PTPN22) is associated with rheumatoid arthritis. Am J Hum Genet. 2004, 75: 330-337. 10.1086/422827.
Plenge RM, Seielstad M, Padyukov L, Lee AT, Remmers EF, Ding B, Liew A, Khalili H, Chandrasekaran A, Davies LR, Li W, Tan AK, Bonnard C, Ong RT, Thalamuthu A, Pettersson S, Liu C, Tian C, Chen WV, Carulli JP, Beckman EM, Altshuler D, Alfredsson L, Criswell LA, Amos CI, Seldin MF, Kastner DL, Klareskog L, Gregersen PK: TRAF1-C5 as a risk locus for rheumatoid arthritis-A genomewide study. N Engl J Med. 2007, 357: 1199-1209. 10.1056/NEJMoa073491.
Raychaudhuri S, Remmers EF, Lee AT, Hackett R, Guiducci C, Burtt NP, Gianniny L, Korman BD, Padyukov L, Kurreeman FA, Chang M, Catanese JJ, Ding B, Wong S, Helm-van Mil van der AH, Neale BM, Coblyn J, Cui J, Tak PP, Wolbink GJ, Crusius JB, Horst-Bruinsma van der IE, Criswell LA, Amos CI, Seldin MF, Kastner DL, Ardlie KG, Alfredsson L, Costenbader KH, Altshuler D, Huizinga TW, Shadick NA, Weinblatt ME, de Vries N, Worthington J, Seielstad M, Toes RE, Karlson EW, Begovich AB, Klareskog L, Gregersen PK, Daly MJ, Plenge RM: Common variants at CD40 and other loci confer risk of rheumatoid arthritis. Nat Genet. 2008, 40: 1216-1223. 10.1038/ng.233.
Lee HS, Korman BD, Le JM, Kastner DL, Remmers EF, Gregersen PK, Bae SC: Genetic risk factors for rheumatoid arthritis differ in Caucasian and Korean population. Arthritis Rheum. 2009, 60: 364-371. 10.1002/art.24245.
Gorlov IP, Gallick GE, Gorlova OY, Amos C, Logothetis CJ: GWAS meets microarray: are the results of genome wide association studies and gene-expression profiling consistent? Prostate cancer as an example. Plos One. 2009, 4: e6511-10.1371/journal.pone.0006511.
This study was financially supported partly by the grant from the Ministry of Health, Labor and Welfare of Japan.
An erratum to this article is available at http://dx.doi.org/10.1186/ar3023.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Sugino, H., Lee, HM. & Nishimoto, N. DNA microarray analysis of rheumatoid arthritis susceptibility genes identified by genome-wide association studies. Arthritis Res Ther 12, 401 (2010). https://doi.org/10.1186/ar2937
- Rheumatoid Arthritis
- Rheumatoid Arthritis Patient
- Susceptibility Gene
- Rheumatoid Arthritis Group
- CTLA4 Expression