Interferon regulatory factor 5 genetic variants are associated with cardiovascular disease in patients with rheumatoid arthritis
- Mercedes García-Bermúdez†1,
- Raquel López-Mejías†2,
- Fernanda Genre†2,
- Santos Castañeda3,
- Javier Llorca4,
- Carlos González-Juanatey5,
- Alfonso Corrales2,
- Begoña Ubilla2,
- José A Miranda-Filloy6,
- Trinitario Pina2,
- Carmen Gómez-Vaquero7,
- Luis Rodríguez-Rodríguez8,
- Benjamín Fernández-Gutiérrez8,
- Alejandro Balsa9,
- Dora Pascual-Salcedo9,
- Francisco J López-Longo10,
- Patricia Carreira11,
- Ricardo Blanco2,
- Javier Martín1 and
- Miguel A González-Gay2, 12Email author
© García-Bermúdez et al.; licensee BioMed Central Ltd. 2014
Received: 11 March 2014
Accepted: 17 June 2014
Published: 10 July 2014
Rheumatoid arthritis (RA) is a complex polygenic inflammatory disease associated with accelerated atherosclerosis and increased cardiovascular (CV) disease risk. Interferon regulatory factor 5 (IRF5) is a regulator of type I interferon induction. Recently, researchers have described an association between multiple single-nucleotide polymorphisms of the IRF5 gene and some rheumatic disorders. In this study, we aimed to evaluate whether three different haplotype blocks within the IRF5 locus which have been shown to alter the protein function are involved in the risk of CV events occurring in Spanish RA patients.
Three IRF5 polymorphisms (rs2004640, rs2070197 and rs10954213) representative of each haplotype group were genotyped by performing TaqMan assays using a 7900HT Fast Real-Time PCR System with tissue from a total of 2,137 Spanish patients diagnosed with RA. Among them, 390 (18.2%) had experienced CV events. The relationship of IRF5 genotypes and haplotypes to CV events was tested using Cox regression.
Male sex, age at RA diagnosis and most traditional risk factors (hypertension, dyslipidemia and smoking habit) were associated with increased risk for CV events in the RA population. Interestingly, a protective effect of both IRF5 rs2004640 GG and IRF5 rs10954213 GG genotypes against the risk for CV events after adjusting the results for sex, age at RA diagnosis and traditional CV disease risk factors was observed (hazard ratio (HR) = 0.6, 95% confidence interval (CI) = 0.38 to 0.92, P = 0.02; and HR = 0.58, 95% CI = 0.36 to 0.95, P = 0.03, respectively). Moreover, we detected a protective effect of the GTG haplotype against the risk for CV events after adjusting the results for potential confounding factors (HR = 0.72, 95% CI = 0.56 to 0.93, P = 0.012).
Our results reveal that IRF5 gene variants are associated with risk of CV events in patients with RA.
Rheumatoid arthritis (RA) is a chronic inflammatory rheumatic disease associated with an increased risk for cardiovascular (CV) events and CV disease–related deaths compared to the general population . Because of that, adequate stratification of CV disease risk has special relevance in RA patients. Researchers in several studies have demonstrated that RA is an independent risk factor for premature heart disease . This process is the result of a combined effect of traditional CV disease risk factors , the magnitude and severity of a chronic inflammatory response  and genetic factors located inside  and outside the human leukocyte antigen (HLA) region [5–8].
Type I interferons (IFNs) are signaling molecules involved in both innate and adaptive immunity. Misregulated expression of type I IFN genes has been observed in peripheral white blood cells of patients with several autoimmune diseases . In this regard, interferon regulatory factor 5 (IRF5) has been involved in the regulation of type I IFN gene transcription [10, 11]. This protein is critical for the production of proinflammatory cytokines , such as TNF-α, interleukin 12 (IL-12) and IL-6, following Toll-like receptor signaling. IRF5 also acts as a molecular switch that controls inflammatory mechanisms mediated by macrophage cells .
IRF5 protein is encoded by the IRF5 gene, which is located in chromosome 7q32.1. Several genetic studies have described the relevance of multiple single-nucleotide polymorphisms (SNPs) of the IRF5 gene in different rheumatic disorders, such as RA and lupus erythematosus [13, 14]. In this context, three groups of correlated IRF5 variants, designated as groups 1, 2 and 3, have been found to be independently associated with these inflammatory disorders and with different functional roles. Group 1 includes SNPs tagging a 30-bp in-frame insertion/deletion variant of exon 6 that alters protein stability. The association of group 2 SNPs seems to be explained by the T allele of IRF5 rs2004640 polymorphism that allows the expression of an alternative isoform and is associated with significantly higher levels of IRF5 expression. The association of group 3 SNPs is probably due to the rs10954213 A allele that creates an early polyadenylation site which leads to higher IRF5 expression .
Taking into account all of these considerations together in the present study, we analyzed the potential role of one tagging SNP of each group of independent signals (rs2004640, rs2070197 and rs10954213) in the risk for CV disease in a large and well-characterized cohort of patients with RA.
Patients and study protocol
A set of 2,137 Spanish patients with RA were included in the present study. Blood samples were obtained from patients recruited from University Hospital Lucus Augusti (Lugo), Marqués de Valdecilla University Hospital (Santander), Hospital de Bellvitge (Barcelona) Hospital Clínico San Carlos, La Paz University Hospital, Hospital de La Princesa, University Hospital Gregorio Marañón and 12 de Octubre University Hospital (Madrid). Informed written consent was obtained from all participants. The study was approved by the ethics committees of Galicia for University Hospital Lucus Augusti; of Cantabria for Marqués de Valdecilla University Hospital; of Cataluña for Hospital de Bellvitge; and of Madrid for Hospital Clínico San Carlos, La Paz University Hospital, Hospital de La Princesa, University Hospital Gregorio Marañón and 12 de Octubre University Hospital. All the patients fulfilled the RA classification criteria published by the American College of Rheumatology in 1987 and 2010 [15, 16]. In all the cases, patients were assessed for the IRF5 rs2004640, rs2070197 and rs10954213 polymorphisms.
Demographic and clinical characteristics a
% ( n/ N) or mean ± SD
Age at disease onset (yr)
53 ± 14.8
Follow-up duration (yr)
12.5 ± 8.25
Cardiovascular risk factors
Patients with cardiovascular events (total)
Ischemic heart disease
Patient DNA was obtained from peripheral blood using standard methods. The IRF5 rs2004640, rs2070197 and rs10954213 polymorphisms were genotyped with predesigned TaqMan SNP genotyping assays in a 7900HT Fast Real-Time PCR System according to the conditions recommended by the manufacturer (Applied Biosystems, Foster City, CA, USA). Negative controls and duplicate samples were included to check the accuracy of genotyping.
The genotype data were checked for deviation from Hardy-Weinberg equilibrium (HWE) using a previously described method . The relationship of genotypes, alleles and haplotypes to CV events that occurred in the follow-up was tested using Cox regression adjusted for sex, age at RA diagnosis and traditional CV disease risk factors. For that purpose, we used the most frequent genotype, allele and haplotype as reference variables. The end of follow-up was the first date of occurrence of the following: end of the study period, date of death or date of CV event. Follow-up time was estimated as the difference between the RA diagnosis date and the end of follow-up. Patients without CV events during the follow-up time and those who died by any non–CV-event–related cause were censored. The results are expressed as hazard ratios (HRs) with 95% confidence intervals (CIs). In order to have a reference population, 10,000 replications were generated by assigning CV events in randomized fashion to the actual studied population. The proportion of the risk for CV disease due to age at disease diagnosis, smoking history, hypertension and genetic variants at the IRF5 locus among the RA patients was estimated as the Nelson-Aalen cumulative HR. Comparative values between the different genetic models of inheritance for IRF5 polymorphisms were estimated using the Akaike Information Criterion (AIC). Statistical significance was defined as P ≤ 0.05, and all analyses were performed using Stata SE/12 statistical software (Stata Corp, College Station, TX, USA).
The IRF5 rs2004640, rs2070197 and rs10954213 polymorphism genotype distributions were in Hardy-Weinberg equilibrium. The genotyping success rate was greater than 98% in all cases. Genotype and allele frequencies of the IRF5 rs2004640, rs2070197 and rs10954213 polymorphisms were in agreement with the data in the HapMap database .
Factors associated with cardiovascular events in a series of 2,137 rheumatoid patients a
HR (95% CI)
Males (reference: females)
1.42 (1.07 to 1.87)
Age at RA diagnosis (by each year)
1.07 (1.06 to 1.09)
1.40 (1.029 to 1.92)
1.18 (0.83 to 1.68)
0.93 (0.67 to 1.31)
1.42 (1.05 to 1.93)
1.58 (1.26 to 1.99)
1.07 (0.79 to 1.45)
0.59 (0.38 to 0.92)
1.16 (0.69 to 1.95)
1.32 (0.32 to 5.52)
0.95 (0.71 to 1.26)
0.58 (0.36 to 0.95)
Results of haplotype analysis in rheumatoid arthritis patients with vs. without cardiovascular disease a
HR (95% CI)b
0.72 (0.56 to 0.93)
0.85 (0.61 to 1.18)
1.14 (0.75 to 1.74)
Additional file 1: Table S1 describes the proportion of the CV disease risk disease for RA patients due to age at disease diagnosis, smoking history, hypertension and the genetic variants at the IRF5 locus. Additional file 2: Table S2 displays the genetic model of inheritance of our analysis. Interestingly, as shown in Additional file 1: Table S1, the presence of GG genotypes in the IRF5 rs2004640 and rs1095213 polymorphisms decreased CV disease risk in 0.7% and 1% of patients, respectively (Additional file 1: Table S1). In addition, smoking, age at disease diagnosis and presence of hypertension increased the risk of CV disease (Additional file 1: Table S1). Further, as shown in Additional file 2: Table S2, the AIC estimation for the recessive model showed the lower level in the IRF5 rs2004640 and rs1095213 polymorphisms, which made it the preferable model for our study (Additional file 2: Table S2).
CV disease is the most common cause of premature mortality in patients with RA [1, 2]. The augmented CV disease mortality observed in this pathology is the result of a compound effect mediated by traditional CV disease risk factors and chronic inflammation [3, 4]. Because genes have also been associated with an increased risk of CV disease in RA, research in the past several years has been focused on the search for genetic markers that may improve the identification of RA patients at risk of experiencing CV events [4–8].
Outside the HLA region, IFN pathway genes, which encode cytokines with critical modulatory effects on innate and adaptive immunity, have been shown to represent a key component of the genetic network leading to autoimmune processes. In this context, several polymorphisms in the IRF5 gene (a regulator of type I IFN induction) [10, 11] are associated with an increased risk of immune-mediated diseases [13, 14]. Because of that, in this study we analyzed three genetic variants (IRF5 rs2004640, rs2070197 and rs10954213) as representatives of three different IRF5 haplotype blocks . To the best of our knowledge, our study constitutes the largest series of RA patients assessed for the potential influence of IRF5 polymorphisms in the risk of CV disease. Interestingly, when we studied each of the polymorphisms separately, our results revealed a protective effect of the IRF5 rs2004640 GG and IRF5 rs10954213 GG genotype variants against the risk of CV events in RA. Moreover, when we analyzed all the genetic variants together to create haplotypes, our results revealed a protective effect of the GTG haplotype (the one that carries both G alleles of IRF5 rs2004640 and rs10954213 polymorphisms) against the risk of CV disease. Because the protective effect of the mutant allele of the IRF5 rs2004640 and rs1095213 polymorphisms is seen only in homozygosis, we can conclude that our study is adjusted to a recessive model of inheritance. Additionally, the results derived from our study are in accordance with those obtained in other pathologies [20, 21]. In this context, the results of a recent study of individuals diagnosed with RA revealed an association between the IRF5 rs2004640 polymorphism and subclinical atherosclerosis .
The results obtained in this study provide additional evidence on the potential role that genetic factors may play in the development of CV disease in RA. The search for genetic markers associated with CV disease in RA may be important to establishing a better characterization of RA patients at risk of CV disease. Improved understanding of these underlying genetic may be useful to establish future therapeutic targets to decrease the risk of CV disease in RA patients.
Our results reveal that IRF5 gene variants are associated with CV disease in RA patients.
JM and MAGG share senior authorship of this article.
American College of Rheumatology
Anti–cyclic citrullinated protein/peptide antibody
Human leukocyte antigen
Interferon regulatory factor
Systemic lupus erythematosus
Tumor necrosis factor.
We wish to thank all the patients with RA who participated and made this study possible. We want to especially thank Sofía Vargas, Rodrigo Ochoa, M Luisa López, M Jesús Ibañez and Sara Olavarria for their technical assistance. This study was supported by European Union FEDER funds (European Fund for Regional Development) and the Health Research Fund (Fondo de Investigación Sanitaria grants PI06/0024, PS09/00748 and PI12/00060) from Instituto de Salud Carlos III (ISCIII, Health Ministry, Spain). It was also partially supported by Cooperative Health Research Thematic Network (RETICS) Program RD12/0009 (RIER) from Instituto de Salud Carlos III (ISCIII, Health Ministry, Spain) and by grants from the European Innovative Medicines Initiative Be the Cure (IMI BTCure) program. RLM is a recipient of a Sara Borrell postdoctoral fellowship from the Instituto Carlos III de Salud at the Spanish Ministry of Health (CD12/00425). FG and BU are supported by funds from the RETICS Program (RIER) (RD12/0009/0013).
- Avina-Zubieta JA, Thomas J, Sadatsafavi M, Lehman AJ, Lacaille D: Risk of incident cardiovascular events in patients with rheumatoid arthritis: a meta-analysis of observational studies. Ann Rheum Dis. 2012, 71: 1524-1529.View ArticlePubMedGoogle Scholar
- Solomon DH, Goodson NJ, Katz JN, Weinblatt ME, Avorn J, Setoguchi S, Canning C, Schneeweiss S: Patterns of cardiovascular risk in rheumatoid arthritis. Ann Rheum Dis. 2006, 65: 1608-1612.PubMed CentralView ArticlePubMedGoogle Scholar
- Dessein PH, Norton GR, Woodiwiss AJ, Joffe BI, Wolfe F: Influence of nonclassical cardiovascular risk factors on the accuracy of predicting subclinical atherosclerosis in rheumatoid arthritis. J Rheumatol. 2007, 34: 943-951.PubMedGoogle Scholar
- Gonzalez-Gay MA, Gonzalez-Juanatey C, Lopez-Diaz MJ, Pineiro A, Garcia-Porrua C, Miranda-Filloy JA, Ollier WE, Martin J, Llorca J: HLA-DRB1 and persistent chronic inflammation contribute to cardiovascular events and cardiovascular mortality in patients with rheumatoid arthritis. Arthritis Rheum. 2007, 57: 125-132.View ArticlePubMedGoogle Scholar
- García-Bermúdez M, González-Juanatey C, López-Mejías R, Teruel M, Corrales A, Miranda-Filloy JA, Castañeda S, Balsa A, Fernández-Gutierrez B, González-Álvaro I, Gómez-Vaquero C, Blanco R, Llorca J, Martín J, González-Gay MA: Study of association of CD40-CD154 gene polymorphisms with disease susceptibility and cardiovascular risk in Spanish rheumatoid arthritis patients. PLoS One. 2012, 7: e49214-PubMed CentralView ArticlePubMedGoogle Scholar
- García-Bermúdez M, López-Mejías R, Genre F, Castañeda S, González-Juanatey C, Llorca J, Corrales A, Miranda-Filloy JA, Rueda-Gotor J, Gómez-Vaquero C, Rodríguez-Rodríguez L, Fernández-Gutiérrez B, Pascual-Salcedo D, Balsa A, López-Longo FJ, Carreira P, Blanco R, González-Álvaro I, Martín J, González-Gay MA: SMAD3 rs17228212 gene polymorphism is associated with reduced risk to cerebrovascular accidents and subclinical atherosclerosis in anti-CCP negative Spanish rheumatoid arthritis patients. PLoS One. 2013, 8: e77695-PubMed CentralView ArticlePubMedGoogle Scholar
- López-Mejías R, García-Bermúdez M, González-Juanatey C, Castañeda S, Miranda-Filloy JA, Gómez-Vaquero C, Fernández-Gutiérrez B, Balsa A, Pascual-Salcedo D, Blanco R, González-Álvaro I, Llorca J, Martín J, González-Gay MA: NFKB1-94ATTG ins/del polymorphism (rs28362491) is associated with cardiovascular disease in patients with rheumatoid arthritis. Atherosclerosis. 2012, 224: 426-429.View ArticlePubMedGoogle Scholar
- López-Mejías R, Genre F, García-Bermúdez M, Corrales A, González-Juanatey C, Llorca J, Miranda-Filloy JA, Rueda-Gotor J, Blanco R, Castañeda S, Martín J, González-Gay MA: The ZC3HC1 rs11556924 polymorphism is associated with increased carotid intima-media thickness in patients with rheumatoid arthritis. Arthritis Res Ther. 2013, 15: R152-PubMed CentralView ArticlePubMedGoogle Scholar
- Higgs BW, Liu Z, White B, Zhu W, White WI, Morehouse C, Brohawn P, Kiener PA, Richman L, Fiorentino D, Greenberg SA, Jallal B, Yao Y: Patients with systemic lupus erythematosus, myositis, rheumatoid arthritis and scleroderma share activation of a common type I interferon pathway. Ann Rheum Dis. 2011, 70: 2029-2036.View ArticlePubMedGoogle Scholar
- Schoenemeyer A, Barnes BJ, Mancl ME, Latz E, Goutagny N, Pitha PM, Fitzgerald KA, Golenbock DT: The interferon regulatory factor, IRF5, is a central mediator of Toll-like receptor 7 signaling. J Biol Chem. 2005, 280: 17005-17012.View ArticlePubMedGoogle Scholar
- Barnes BJ, Moore PA, Pitha PM: Virus-specific activation of a novel interferon regulatory factor, IRF-5, results in the induction of distinct interferon α genes. J Biol Chem. 2001, 276: 23382-23390.View ArticlePubMedGoogle Scholar
- Krausgruber T, Blazek K, Smallie T, Alzabin S, Lockstone H, Sahgal N, Hussell T, Feldmann M, Udalova IA: IRF5 promotes inflammatory macrophage polarization and TH1-TH17 responses. Nat Immunol. 2011, 12: 231-238.View ArticlePubMedGoogle Scholar
- Eyre S, Bowes J, Diogo D, Lee A, Barton A, Martin P, Zhernakova A, Stahl E, Viatte S, McAllister K, Amos CI, Padyukov L, Toes RE, Huizinga TW, Wijmenga C, Trynka G, Franke L, Westra HJ, Alfredsson L, Hu X, Sandor C, de Bakker PI, Davila S, Khor CC, Heng KK, Andrews R, Edkins S, Hunt SE, Langford C, Symmons D: High-density genetic mapping identifies new susceptibility loci for rheumatoid arthritis. Nat Genet. 2012, 44: 1336-1340.PubMed CentralView ArticlePubMedGoogle Scholar
- Graham RR, Kozyrev SV, Baechler EC, Reddy MV, Plenge RM, Bauer JW, Ortmann WA, Koeuth T, González Escribano MF, Pons-Estel B, Petri M, Daly M, Gregersen PK, Martín J, Altshuler D, Behrens TW, Alarcón-Riquelme ME, Argentine and Spanish Collaborative Groups: A common haplotype of interferon regulatory factor 5 (IRF5) regulates splicing and expression and is associated with increased risk of systemic lupus erythematosus. Nat Genet. 2006, 38: 550-555.View ArticlePubMedGoogle Scholar
- Arnett FC, Edworthy SM, Bloch DA, McShane DJ, Fries JF, Cooper NS, Healey LA, Kaplan SR, Liang MH, Luthra HS, Medsger TA, Mitchell DM, Neustadt DH, Pinals RS, Schaller JG, Sharp JT, Wilder RL, Hunder GG: The American Rheumatism Association 1987 revised criteria for the classification of rheumatoid arthritis. Arthritis Rheum. 1988, 31: 315-324.View ArticlePubMedGoogle Scholar
- Aletaha D, Neogi T, Silman AJ, Funovits J, Felson DT, Bingham CO, Birnbaum NS, Burmester GR, Bykerk VP, Cohen MD, Combe B, Costenbader KH, Dougados M, Emery P, Ferraccioli G, Hazes JM, Hobbs K, Huizinga TW, Kavanaugh A, Kay J, Kvien TK, Laing T, Mease P, Ménard HA, Moreland LW, Naden RL, Pincus T, Smolen JS, Stanislawska-Biernat E, Symmons D: 2010 rheumatoid arthritis classification criteria: an American College of Rheumatology/European League Against Rheumatism collaborative initiative. Arthritis Rheum. 2010, 62: 2569-2581.View ArticlePubMedGoogle Scholar
- Gonzalez-Juanatey C, Llorca J, Martin J, Gonzalez-Gay MA: Carotid intima-media thickness predicts the development of cardiovascular events in patients with rheumatoid arthritis. Semin Arthritis Rheum. 2009, 38: 366-371.View ArticlePubMedGoogle Scholar
- Institut für Humangenetik: Case–Control Studies: Tests for Deviation from Hardy-Weinberg Equilibrium and Tests for Association. [http://ihg.gsf.de/cgi-bin/hw/hwa1.pl]
- International HapMap Project. [http://hapmap.ncbi.nlm.nih.gov/]
- Carmona FD, Martin JE, Beretta L, Simeón CP, Carreira PE, Callejas JL, Fernández-Castro M, Sáez-Comet L, Beltrán E, Camps MT, Egurbide MV, Airó P, Scorza R, Lunardi C, Hunzelmann N, Riemekasten G, Witte T, Kreuter A, Distler JHW, Madhok R, Shiels P, van Laar JM, Fonseca C, Denton C, Herrick A, Worthington J, Schuerwegh AJ, Vonk MC, Voskuyl AE, the Spanish Scleroderma Group: The systemic lupus erythematosus IRF5 risk haplotype is associated with systemic sclerosis. PLoS One. 2013, 8: e54419-PubMed CentralView ArticlePubMedGoogle Scholar
- Fan JH, Gao LB, Pan XM, Li C, Liang WB, Liu J, Li Y, Zhang L: Association between IRF-5 polymorphisms and risk of acute coronary syndrome. DNA Cell Biol. 2010, 29: 19-23.View ArticlePubMedGoogle Scholar
- Vosslamber S, van Sijl A, Bos CL, Lubbers J, de Ridder S, Voskuyl AE, Nurmohamed MT, Verweij CL: Interferon regulatory factor 5 (IRF5) gene variant rs2004640 is associated with carotid intima media thickness in rheumatoid arthritis patients [abstract A1.77]. Ann Rheum Dis. 2014, 73: A34-View ArticleGoogle 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/4.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.