Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Gender-Specific Effects of Genetic Variants within Th1 and Th17 Cell-Mediated Immune Response Genes on the Risk of Developing Rheumatoid Arthritis

  • Rafael Cáliz ,

    Contributed equally to this work with: Rafael Cáliz, Luz María Canet

    Affiliation Rheumatology Department, Virgen de las Nieves University Hospital, Granada, Spain

  • Luz María Canet ,

    Contributed equally to this work with: Rafael Cáliz, Luz María Canet

    Affiliation Genomic Oncology Area, GENYO, Centre for Genomics and Oncological Research: Pfizer/University of Granada/Andalusian Regional Government, Granada, Spain

  • Carmen Belén Lupiañez,

    Affiliation Genomic Oncology Area, GENYO, Centre for Genomics and Oncological Research: Pfizer/University of Granada/Andalusian Regional Government, Granada, Spain

  • Helena Canhão,

    Affiliations Rheumatology Research Unit, Instituto de Medicina Molecular, Faculdade de Medicina da Universidade de Lisboa, Lisbon, Portugal, Rheumatology Department, Santa Maria Hospital–CHLN, Lisbon, Portugal

  • Alejandro Escudero,

    Affiliation Rheumatology Department, Reina Sofía Hospital, Córdoba, Spain

  • Ileana Filipescu,

    Affiliation Rheumatology Department, Reina Sofía Hospital, Córdoba, Spain

  • Juana Segura-Catena,

    Affiliation Genomic Oncology Area, GENYO, Centre for Genomics and Oncological Research: Pfizer/University of Granada/Andalusian Regional Government, Granada, Spain

  • María José Soto-Pino,

    Affiliation Rheumatology Department, Virgen de las Nieves University Hospital, Granada, Spain

  • Manuela Expósito-Ruiz,

    Affiliation FIBAO Fundation, Virgen de las Nieves University Hospital, Granada, Spain

  • Miguel Ángel Ferrer,

    Affiliation Rheumatology Department, Virgen de las Nieves University Hospital, Granada, Spain

  • Antonio García,

    Affiliation Rheumatology Department, Virgen de las Nieves University Hospital, Granada, Spain

  • Lurdes Romani,

    Affiliation Rheumatology Department, Virgen de las Nieves University Hospital, Granada, Spain

  • Alfonso González-Utrilla,

    Affiliation Rheumatology Department, Virgen de las Nieves University Hospital, Granada, Spain

  • Teresa Vallejo,

    Affiliation Rheumatology Department, Virgen de las Nieves University Hospital, Granada, Spain

  • Eva Pérez-Pampin,

    Affiliation Rheumtology Unit, Instituto de Investigacion Sanitaria-Hospital Clinico Universitario de Santiago, Santiago de Compostela, Spain

  • Kari Hemminki,

    Affiliations Division of Molecular Genetic Epidemiology, German Cancer Research Center, Heidelberg, Germany, Center for Primary Health Care Research, Clinical Research Center, Malmö, Sweden

  • Asta Försti,

    Affiliations Division of Molecular Genetic Epidemiology, German Cancer Research Center, Heidelberg, Germany, Center for Primary Health Care Research, Clinical Research Center, Malmö, Sweden

  • Eduardo Collantes,

    Affiliation Rheumatology Department, Reina Sofía Hospital, Córdoba, Spain

  • João Eurico Fonseca,

    Affiliations Rheumatology Research Unit, Instituto de Medicina Molecular, Faculdade de Medicina da Universidade de Lisboa, Lisbon, Portugal, Rheumatology Department, Santa Maria Hospital–CHLN, Lisbon, Portugal

  •  [ ... ],
  • Juan Sainz

    juan.sainz@genyo.es

    Affiliations Genomic Oncology Area, GENYO, Centre for Genomics and Oncological Research: Pfizer/University of Granada/Andalusian Regional Government, Granada, Spain, Division of Molecular Genetic Epidemiology, German Cancer Research Center, Heidelberg, Germany

  • [ view all ]
  • [ view less ]

Abstract

The present study was conducted to explore whether single nucleotide polymorphisms (SNPs) in Th1 and Th17 cell-mediated immune response genes differentially influence the risk of rheumatoid arthritis (RA) in women and men. In phase one, 27 functional/tagging polymorphisms in C-type lectins and MCP-1/CCR2 axis were genotyped in 458 RA patients and 512 controls. Carriers of Dectin-2rs4264222T allele had an increased risk of RA (OR = 1.47, 95%CI 1.10–1.96) whereas patients harboring the DC-SIGNrs4804803G, MCP-1rs1024611G, MCP-1rs13900T and MCP-1rs4586C alleles had a decreased risk of developing the disease (OR = 0.66, 95%CI 0.49–0.88; OR = 0.66, 95%CI 0.50–0.89; OR = 0.73, 95%CI 0.55–0.97 and OR = 0.68, 95%CI 0.51–0.91). Interestingly, significant gender-specific differences were observed for Dectin-2rs4264222 and Dectin-2rs7134303: women carrying the Dectin-2rs4264222T and Dectin-2rs7134303G alleles had an increased risk of RA (OR = 1.93, 95%CI 1.34–2.79 and OR = 1.90, 95%CI 1.29–2.80). Also five other SNPs showed significant associations only with one gender: women carrying the MCP-1rs1024611G, MCP-1rs13900T and MCP-1rs4586C alleles had a decreased risk of RA (OR = 0.61, 95%CI 0.43–0.87; OR = 0.67, 95%CI 0.47–0.95 and OR = 0.60, 95%CI 0.42–0.86). In men, carriers of the DC-SIGNrs2287886A allele had an increased risk of RA (OR = 1.70, 95%CI 1.03–2.78), whereas carriers of the DC-SIGNrs4804803G had a decreased risk of developing the disease (OR = 0.53, 95%CI 0.32–0.89). In phase 2, we genotyped these SNPs in 754 RA patients and 519 controls, leading to consistent gender-specific associations for Dectin-2rs4264222, MCP-1rs1024611, MCP-1rs13900 and DC-SIGNrs4804803 polymorphisms in the pooled sample (OR = 1.38, 95%CI 1.08–1.77; OR = 0.74, 95%CI 0.58–0.94; OR = 0.76, 95%CI 0.59–0.97 and OR = 0.56, 95%CI 0.34–0.93). SNP-SNP interaction analysis of significant SNPs also showed a significant two-locus interaction model in women that was not seen in men. This model consisted of Dectin-2rs4264222 and Dectin-2rs7134303 SNPs and suggested a synergistic effect between the variants. These findings suggest that Dectin-2, MCP-1 and DC-SIGN polymorphisms may, at least in part, account for gender-associated differences in susceptibility to RA.

Introduction

Rheumatoid arthritis (RA) is an autoimmune disease characterized by chronic inflammatory activity in the synovial joints often leading to progressive cartilage and bone destruction [1]. Although etiology of RA is still unknown, it has been suggested that activated macrophages and dendritic cells (DCs), rather than T-cells, may be important for the onset of the disease. In the RA synovium, macrophages mature into tissue macrophages and differentiate into DCs leading to massive release of a wide number of pro-inflammatory mediators [2][4] and to the induction of adhesion (C-type lectins) and co-stimulatory molecules (CD40, CD80, CD86) that participate in the activated macrophage- and DC-induced T-cell proliferation [5], [6]. Among all these immune mediators, MCP-1/CCL2, a chemokine highly expressed in the synovial fluid of RA patients [7], binds to CCR2 and promotes the recruitment of antigen-presenting cells and T cells [8] whereas C-type lectins such as Dectin-1, Dectin-2 and DC-SIGN regulate the monocyte-induced T-cell activation in the RA synovium [5], [6]. C-type lectins also modulate adaptive immune responses, of which especially Th17 responses are implicated in the stimulation of osteoclastogenesis and bone destruction during RA [9].

RA is three times more frequent in women than men [10] and emerging data also suggest that women are more likely to present a worse course of the disease and to become severely disabled [11], [12]. Although different hypotheses concerning sex-related differences in RA incidence and severity have been generated, the hypothesis suggesting a sexual dimorphism in the intensity of immune responses remains as one of the most probable mechanisms in promoting and establishing a different synovial membrane inflammation and, subsequently, different levels of cartilage and bone erosion [13], [14]. Women have higher immunoglobulin levels than men, they show stronger Th1-type cell-mediated immune responses and they have higher absolute numbers of CD4+ lymphocytes and a more proactive cytokine profile [15], [16], which likely contribute to their increased autoimmune responses [17]. Experimental studies in rodents have also shown different immune responses between female and male animals and an equivalent sexual dimorphism in the incidence of RA [18]. The underlying reasons for this gender bias are still unknown, but studies in monozygotic twins have suggested that genetic factors may, at least in part, account for sex-related differences in the immune responses [19] and, consequently, in the susceptibility to autoimmune diseases. Genetic factors implicated in RA have been widely studied using both candidate genes [20] and whole-genome screens [21] but, so far, only a few studies have investigated the link between SNPs and the gender-associated differences in susceptibility to RA [22], [23]. Considering these facts, the present study was designed to evaluate the influence of 27 tagging and potentially functional SNPs within the MCP-1, CCR2, DC-SIGN, Dectin-1, and Dectin-2 genes in the risk to develop RA in women and men, separately.

Materials and Methods

Study Population

In phase 1, all participants enrolled were Caucasian and recruited at the department of Rheumatology of the Virgen de las Nieves (Granada, Spain) and Reina Sofia (Córdoba, Spain) hospitals. All participants gave their written informed consent to participate in the study, which was approved by the ethical review committee of participant institutions (Virgen de las Nieves University Hospital, Granada, Spain; Reina Sofia Hospital, Córdoba, Spain). The study was performed according to the principles of the Declaration of Helsinki. The population consisted of 970 participants, 458 RA patients (360 women and 98 men) and 512 healthy controls (217 women and 295 men). Rheumatoid patients were treated at their respective departments of Rheumatology from January 2004 to January 2010. The diagnosis of RA was assigned by physician investigators and fulfilled the 1987 American College of Rheumatology (ACR) criteria. We chose DAS28 as a measure of disease activity as it is a validated score for established RA. Moderate to high activity disease was defined as DAS28≥3.2 while low disease activity was defined as DAS28<3.2. Controls were blood donor subjects randomly recruited at the Regional Blood Transfusion and Tissue Centre (Granada-Almería, Spain).

In phase 2, in order to increase the statistical power of the study and confirm both overall and gender-specific associations, we extended the study by recruiting additional RA cases (n = 831) and controls (n = 550) from our own institution as well as from other collaborating institutions. Seven hundred and seventy-three RA patients and 201 controls were recruited from the Santa Maria Hospital–CHLN (Lisbon, Portugal). Fifty-eight additional RA patients were recruited from the University Clinical Hospital of Santiago de Compostela (Santiago de Compostela, Spain). We also recruited 349 healthy controls from our own institution (Virgen de las Nieves University hospital, Granada, Spain; n = 260) as well as from the Reina Sofía hospital (Córdoba, Spain; n = 89). All participants gave their written informed consent to participate in the study, which was approved by the respective ethical review committee of participant institutions. Thirty-four non-Caucasian subjects were excluded from the statistical analyses (n = 34) and some additional patients and controls were also removed for technical reasons (low DNA quality, unknown age or gender, etc.). Seven hundred and fifty-four RA patients (626 women and 128 men) and 519 healthy controls (348 women and 171 men) were finally available for genotyping.

SNP Selection and Genotyping

Twenty-seven tagging/functional SNPs within DC-SIGN, Dectin-1, Dectin-2, MCP-1 and CCR2 were selected to genotype the entire panel of individuals (Table 1; Supplementary material). The aim of the SNP tagging was to identify a set of SNPs that efficiently tags all the known SNPs while the functional approach was used to determine the net impact of potentially functional variants within DC-SIGN, Dectin-1, Dectin-2, MCP-1 and CCR2 genes on RA risk. Tagging SNPs were selected using Haploview Tagger program (http://www.broad.mit.edu/mpg/haploview/; http://www.broad.mit.edu/mpg/tagger/). SNPs with a MAF>0.05 were included in the selection of tag SNPs using a pairwise tagging with a minimum r2 of 0.8. In this selection we forced the inclusion of the DC-SIGNrs4804803, MCP-1rs1024610 and MCP-1rs1024611 polymorphisms as their functionality has been suggested [24][26]. One extra SNP within DC-SIGN (rs11465384) and in strong linkage disequilibrium (LD) with at least 4 SNPs was selected as redundant SNP in case of genotyping failure. Genomic DNA was extracted from peripheral blood mononuclear cells (PBMCs) using Qiagen Mini Kit (Qiagen, Valencia, CA, USA). Genotyping of DC-SIGN, Dectin-1, Dectin-2, MCP-1 and CCR2 polymorphisms was carried out using allele-specific KASPar® assays (LGC Genomics KBioscience, London, United Kingdom) in a 384-well plate format (Applied Biosystems, Foster City, CA, USA) where population samples (RA and healthy controls) were randomly distributed. KASPar reactions were performed using KASPar assay mix (containing probes) and KASPar kit containing 2X Reaction Mix and MgCl2 (50 mM). Touch-down PCR conditions were: denaturation at 94°C for 15 min, 10 cycles of denaturation at 94°C for 20 sec, annealing at 61°C for 60 sec (dropping –0.6°C per cycle) and 26 cycles of denaturation at 94°C for 10 sec, annealing at 55°C for 60 sec. Recycling conditions were 94°C for 10 sec, annealing and elongation at 60°C for 60 sec. PCR products were analyzed with the ABI Prism 7900HT detection system using the SDS 2.4 software (Applied Biosystems). For internal quality control, 5% of samples were randomly selected and included as duplicate. Concordance between the original and the duplicate samples for the 27 SNPs analyzed was ≥99.0%. Call rates for all SNPs were ≥95.0% with the exception of the Dectin-1rs11053599 SNP with a call rate <90.0%. This latter SNP was excluded from further analysis.

thumbnail
Table 1. Selected SNPs within DC-SIGN, Dectin-1, Dectin-2, MCP-1 and CCR2 genes.

https://doi.org/10.1371/journal.pone.0072732.t001

Statistical Analysis

The Hardy-Weinberg Equilibrium (HWE) tests were performed in the control group by a standard observed-expected chi-square (χ2) test at each polymorphic site (http://ihg2.helmholtz-muenchen.de/cgi-bin/hw/hwa1.pl). Binary logistic regression was used to assess the main effects of the genetic polymorphisms on RA risk using co-dominant and dominant inheritance models. For each SNP, the more common allele in the control group was assigned as the reference category. Binary logistic regression analyses were adjusted for age and gender whereas gender-stratified analyses were adjusted for age. In the pooled analysis, overall binary logistic regression analyses were adjusted by age, gender and centre whereas gender-stratified analyses were adjusted for age and centre. Statistical power was estimated using Quanto software (http://hydra.usc.edu/gxe/). All analyses were conducted using the statistical software SSPS (version 14.0, SPSS Inc., Chicago, USA). All tests were considered to be statistically significant when p<0.05. In order to correct for multiple testing, we calculated an adjusted significance level using Meff correction. We estimated an individual Meff value for each gene (i.e. the effective number of independent variables) and the study-wise Meff value was estimated by adding up the independent gene Meff values by using the SNP Spectral Decomposition Lite program (http://gump.qimr.edu.au/general/daleN/SNPSpDsuperlite/) [27]. The resulting number of independent marker loci was applied to correct for multiple testing.

Functional Prediction of Associated SNPs

We used a web-based tool FastSNP (available at http://fastsnp.ibms.sinica.edu.tw) for predicting the functional significance of the SNPs associated with RA. A detailed explanation on the FastSNP tool has been described elsewhere [28]. FastSNP utilizes information from different web resources (Hapmap, dbSNP, NCBI Blast, Polyphen, SNPs3D, TRANSFAC, ESEfinder, Rescue-ESE and FAS-ESS) to determine whether SNPs are located at exonic splicing regulatory sites or whether they alter the transcription factor-binding site of a gene (for instance, acting as intronic enhancer) or even whether SNPs are affecting Micro-RNA binding sites. In addition, these tools are able to identify either nonsense or non-synonymous SNPs that lead to premature termination of translation or amino acid changes and that are, therefore, very likely to affect protein function.

SNP-SNP Interaction Analysis

We also analyzed high-order interactions between significant SNPs using the multifactor dimensionality reduction (MDR) constructive induction algorithm. A detailed explanation on the MDR method has been described elsewhere [29]. Cross-validation and permutation testing were used to identify the best models. All possible two-way SNP interactions were tested using 10-fold cross-validation and the exhaustive search. The model with the highest testing balanced accuracy (TA) and cross validation consistency (CVC) was selected as “best model”. Statistical significance was evaluated by the Sign test and confirmed using a 1.000-fold permutation test to compare observed testing balanced accuracies with those expected under the null hypothesis of no association (using the MDR permutation testing module 0.4.9 alpha) [29]. MDR results were considered statistically significant at the level of 0.05. MDR software and MDR permutation testing module are open-source and freely available from http://www.epistasis.org. Logistic regression analyses were also performed to confirm significant interaction results from MDR analyses.

Results

Demographic and clinical characteristics of the RA patients analysed in phase 1 are described in Table 2. Controls were slightly younger than RA patients (53.46±9.69 vs. 58.51±13.13; p<0.001). Seventy-eight percent of RA patients were female, 70.7% were positive for anti-CCP antibodies and 71.2% were positive for RF (rheumatoid factor). According to the DAS28, male patients showed low active disease (DAS28Male = 2.89) whereas women had moderately active disease (DAS28Female = 3.42). Women tended to have more often a DAS28≥3.2 compared to men (56.18% vs. 43.48%, p = 0.078). The overall average disease duration was 12.39 years.

thumbnail
Table 2. Demographic and clinical characteristics of RA patients.

https://doi.org/10.1371/journal.pone.0072732.t002

All analyzed SNPs were in HWE in the control group with the exception of Dectin-1rs16910631 (p>0.01). Five SNPs showed overall association with RA (Table 3 and Table S1). Carriers of Dectin-2rs4264222T allele had an increased risk of RA (OR = 1.47, 95%CI 1.10–1.96, P = 0.009) whereas patients harboring the DC-SIGNrs4804803G, MCP-1rs1024611G, MCP-1rs13900T and MCP-1rs4586C alleles had a decreased risk of developing the disease (OR = 0.66, 95%CI 0.49–0.88, P = 0.004; OR = 0.66, 95%CI 0.50–0.89, P = 0.006; OR = 0.73, 95%CI 0.55–0.97, P = 0.03 and OR = 0.68, 95%CI 0.51–0.91, P = 0.009; Table 3). In addition, patients bearing the Dectin-2rs7134303G allele showed a trend to have an increased risk for RA (OR = 1.35, 95%CI 1.00–1.83). Most importantly, binary logistic regression analyses revealed gender-specific associations with RA for six SNPs (Table 3). There was a significant effect modification by gender for Dectin-2rs4264222 and Dectin-2rs7134303 SNPs (Pinteraction = 0.041 and 0.011, respectively; Table 3). Women carrying either the Dectin-2rs4264222T or Dectin-2rs7134303G alleles had an increased risk of RA (OR = 1.93, 95%CI 1.34–2.79 and OR = 1.90, 95%CI 1.29–2.80). Additionally, women carrying the MCP-1rs1024611G, MCP-1rs13900T and MCP-1rs4586C alleles had a decreased risk of RA compared with those carrying the wild-type genotype (OR = 0.61, 95%CI 0.43–0.87; OR = 0.67, 95%CI 0.47–0.95 and OR = 0.60, 95%CI 0.42–0.86). As also men had decreased, but statistically not significant ORs, no effect modification by gender was observed. A borderline gender-specific effect was observed for the DC-SIGNrs2287886 SNP (Pinteraction = 0.071; Table 3). Here, male carriers of the DC-SIGNrs2287886A allele had an increased risk of RA (OR = 1.70, 95%CI 1.03–2.78). Additionally, the association of the DC-SIGNrs4804803G with a decreased risk was stronger in men than in women (OR = 0.53, 95%CI 0.32–0.89 vs. OR = 0.73, 95%CI 0.51–1.04).

thumbnail
Table 3. Dectin-2, DC-SIGN and MCP-1 polymorphisms associated with rheumatoid arthritis.

https://doi.org/10.1371/journal.pone.0072732.t003

After correction for multiple testing using the SNP Spectral Decomposition Lite program (number of independent marker loci, 21; p = 0.05/21 = 0.002), only the Dectin-2rs4264222 and Dectin-2rs7134303 associations in women retained significance (P = 3.00 E-04 and P = 9.00 E-04), whereas the associations with DC-SIGNrs4804803, MCP-1rs1024611 and MCP-1rs4586 showed borderline significance in the whole population (P = 0.004, P = 0.006 and P = 0.009; Table 3).

In phase 2, in order to increase the statistical power of this study and to confirm the significant associations, we further genotyped the significant SNPs in 754 Caucasian RA patients (626 women and 128 men) and 519 healthy controls (348 women and 171 men). Clinical characteristics of the RA patient population are shown in Table S2 and characteristics of the pooled population are shown in Table S3. In the phase 2, RA patients had a similar age and showed slightly higher percentages of anti-CCP (77.72 vs. 70.70) and RF positivity (75.98 vs. 71.20). As before, male patients showed low active disease when compared to women (DAS28Male = 3.21 vs. DAS28Female = 3.78). Thus, we end up with a total population of 2.243 individuals including 1.212 RA patients (986 women and 226 men) and 1.031 controls (565 women and 466 men). After this recruitment the study had over 80% power (co-dominant model) to detect an odds ratio of 1.19 at alpha = 0.008 (multiple testing threshold) for a polymorphism with a minor allele frequency of 0.25. Although the gender-stratified analysis reduced the statistical power of the study, we still had power to detect reasonably small risks (ORWOMEN = 1.23 and ORMEN = 1.36).

Binary logistic regression analysis adjusted for age, gender and center in the pooled population showed that MCP-1rs1024611, MCP-1rs13900, MCP-1rs4586 and DC-SIGNrs4804803 polymorphisms were associated with a decreased risk of developing RA (OR = 0.76, 95%CI 0.61–0.95, P = 0.01; OR = 0.77, 95%CI 0.62–0.96, P = 0.02; OR = 0.80, 95%CI 0.64–0.99, P = 0.04; and OR = 0.77, 95%CI 0.62–0.96, P = 0.02; Table 4). Although there was no significant effect modification by gender, Dectin-2 and MCP-1 polymorphisms seemed to have an stronger effect in women compared to men, whereas variants within DC-SIGN showed more evident effects in men. Thus, women carrying the Dectin-2rs4264222T allele had an increased risk of RA (OR = 1.38, 95%CI 1.08–1.77, P = 0.01) whereas women harbouring the MCP-1rs1024611G and MCP-1rs13900T alleles had a decreased risk of RA in comparison with those carrying the wild-type genotype (OR = 0.74, 95%CI 0.58–0.94, P = 0.02 and OR = 0.76, 95%CI 0.59–0.97, P = 0.03; Table 4). Women bearing the Dectin-2rs7134303G showed also a trend towards an increased risk of the disease (OR = 1.30, 95%CI 1.00–1.69, P = 0.05) whereas those women carrying the MCP-1rs4586C allele had a decreased risk for RA (OR = 0.78, 95%CI 0.61–1.00, P = 0.05; Table 4). We could confirm that none of these effects was observed in men. We also confirmed that men carrying the DC-SIGNrs4804803G allele had a decreased risk to develop RA (OR = 0.56, 95%CI 0.34–0.93, P = 0.02; Table 4) whereas men harboring DC-SIGNrs2287886A allele had trend to have an increased risk of RA (OR = 1.54, 95%CI 0.96–2.46, P = 0.07). None of these associations resisted multiple testing adjustments in the pooled sample (number of independent marker loci was estimated including only those SNPs genotyped in the pooled sample; p = 0.05/6 = 0.008) and all require further replication in independent populations.

thumbnail
Table 4. Overall analysis of Dectin-2, DC-SIGN and MCP-1 polymorphisms with rheumatoid arthritis.

https://doi.org/10.1371/journal.pone.0072732.t004

For predicting the effect of polymorphisms found associated with RA we used FastSNP, which predicts the possible effect of genetic variants on the protein function and/or structure. Among the SNPs associated either with increased or decreased risk of RA, 2 SNPs were predicted to have functional impact. The predictive functional analysis suggested an intronic enhancer function for Dectin-2rs7134303 due to its location in a transcription factor-binding site (risk score 1–2) and a function as gene expression regulator for the MCP-1rs1024611 (risk score 1–3). The presence of the Dectin-2rs7134303G allele creates a nuclear protein-binding site for the transcription factor GATA-1 (GATA-binding factor 1) whereas the presence of the MCP-1rs1024611T allele, which confers protection against RA, disrupts a binding site for both GATA-1 and GATA-2 in the promoter region of the MCP-1 gene. Whereas GATA-1 is a transcription factor expressed in circulating inflammatory monocytes and indispensable for effective DC maturation and survival, GATA-2 expression inhibits monocyte differentiation. These data suggest a central role of the Dectin-2rs7134303 and MCP-1rs1024611 polymorphisms in the gender-associated susceptibility to RA and point out a possible participation of GATA-1 and GATA-2 in the modulation of Dectin-2 and MCP-1 expression and, consequently, in the activation of both monocytes and DCs.

We also investigated the epistatic effect of the significant SNPs in order to find possible gender-specific high order interactions. A summary of the results for the models that had maximum testing accuracy and maximum cross validation consistency in women and men is presented in Table 5. In women, the overall best model consisted of Dectin-2rs4264222 and Dectin-2rs7134303 SNPs that interacted in a synergistic or non-additive manner. Logistic regression analysis applied to this two -locus model confirmed our results (P = 0.007). In men, no SNP-SNP interactions were shown to be significant (Table 5).

thumbnail
Table 5. MDR analysis to detect two-locus disease models.

https://doi.org/10.1371/journal.pone.0072732.t005

Finally, we explored the effect of selected SNPs on disease activity and severity. No differences were found in relation to DAS28, the positivity for anti-CCP antibodies or RF when RA patients were grouped by gender and genotype (data not shown).

Discussion

The idea that genetic factors have an impact on the predisposition to RA is well supported by our understanding of the disease biology and by data from a wide range of genetic epidemiologic studies [30][32]. This fact, along with studies showing sex-specific differences in incidence of RA and in immune response [10], [33], suggests that genetic factors modulating the immune response may contribute to the sex-specific incidence rates. Considering this hypothesis, the present study was carried out to assess whether genetic variants within immune-related genes may differentially play a role in determining the risk of RA in women and men. Although to date there have been too few studies addressing the effect of gender on the risk of RA, some previous studies have reported female-specific associations between genetic variants in immune genes (TNF, TNFR2, IL4R and CD4) and risk of RA [23], [34], [35].

The present population-based case-control study included 2.243 individuals (1.212 RA patients and 1.031 controls) and confirmed the role of 4 SNPs within Dectin-2, MCP-1 and DC-SIGN genes in determining the risk of RA. These results were in agreement with a previous study performed by Platinga et al. (2010), who reported no overall association of Dectin-1rs16910526 SNP with RA [36] but were in contrast with those results described by Dieguez-González et al. (2009) reporting no association of the DC-SIGNrs4804803 and RA [37]. These controversial results might be attributed to differences in population size and environmental factors but also confounding factors such as gender.

We also identified, for the first time, genetic variants in Dectin-2, DC-SIGN and MCP-1 as contributing to the gender-specific risk for RA. Females with the Dectin-2rs4264222T allele had an increased risk to develop RA whereas females carrying any of the MCP-1rs1024611G and MCP-1rs13900T alleles had a decreased risk to develop RA compared with carriers of the wild-type genotype. We also found that women bearing the Dectin-2rs7134303G tended to have an increased risk of the disease whereas those women carrying the MCP-1rs4586C allele showed a borderline significant decreased risk for RA. In males, none of these effects were found but carriers of the DC-SIGNrs4804803G had a significantly decreased risk of developing the disease whereas those harboring the DC-SIGNrs2287886A allele had trend to have an increased risk of RA. Although there was not significant effect modification by gender, Dectin-2, MCP-1 and DC-SIGN polymorphisms showed consistent associations in women that were not seen in men and vice versa. Of note is that, after adjusting the significance level to account for multiple comparisons (study-wise significant P-threshold = 0.008), Dectin-2rs4264222, MCP-1rs1024611, MCP-1rs13900 and MCP-1rs4586 SNPs showed both overall and gender-specific borderline significant associations. These results underlie the potential importance of analyzing RA data both with and without gender as a stratifying factor.

Gender-specific SNP-SNP interactions among the selected variants indicated that the combined effect of the polymorphisms with or without significant main effects conferred risk for RA. The MDR approach used in this study identified a two-locus significant model associated with high risk of RA in women while no significant models were found in men. In women, the strongest model for predicting RA risk was a model including Dectin-2rs4264222 and Dectin-2rs7134303 SNPs that interact in synergistic manner to increase the risk of RA. Although interesting, these results should be interpreted with caution given that the statistical modeling of SNP-SNP interactions may not be assumed as a true biological interaction. However, based on in silico tools, hypotheses concerning the molecular mechanisms resulting in differential activity of genes can be created. In this regard, these two SNPs (Dectin-2rs7134303 and Dectin-2rs4264222) were found to be important in this study, both in the single SNP or the epistatic analysis. MDR analysis was further confirmed by logistic regression analysis. Dectin-2rs7134303 creates an allele-specific change in transcriptional binding site for GATA-1 and could contribute, acting as intronic enhancer, to the transcriptional regulation of Dectin-2 gene and, consequently, to affect its intracellular signaling and downstream gene expression.

According to these results, it is conceivable to suggest that Dectin-2 variants, independently and/or through interactions may, at least partially, account for gender differences in susceptibility to RA. Dectin-2, a type II transmembrane receptor, is mainly expressed in DCs, macrophages and B-cells and its expression may be induced by inflammatory stimuli [38]. It is well known that Dectin-2 recognizes carbohydrate motifs of self and non-self antigens. However, it has also been suggested that Dectin-2 binds to endogenous ligands in CD4+CD25+ [39] and CD8+ T-cells [40] thus triggering inflammatory responses. Most recently, Dectin-2 has been involved in promoting Syk- and CARD9-dependent NFkB activation [41], inducing the expression of a wide variety of pro-inflammatory mediators including cytokines, chemokines and co-stimulatory molecules and, thereby enhancing T-cell mediated immune responses. In line with this, Sato et al. also reported that, contrary to Dectin-1 that binds Syk directly, Dectin-2 couples to Syk through its association with Fc receptor gamma (FcRγ) chain leading to the activation of NFkB pathway and production of cytokines such as TNF and IL1RA [42]. Interestingly, it has also been described that Dectin-2 is involved in the modulation of Th17 immune responses and that the production of IL1β and IL23 upon Dectin-2 requires c-Rel-dependent activation of NFkB [43]. These findings, along with those suggesting the role of CD4+ T cells secreting IL17 (Th17 cells) [44] and NFkB activation in the pathogenesis and/or progression of RA [45], may indicate a key role for Dectin-2 in the exacerbation of the immune response that define the RA. Nonetheless, the molecular mechanisms by which the Dectin-2rs7134303 and Dectin-2rs4264222 polymorphisms might increase the risk for developing RA remain to be clarified. Based on the findings reported by Gutierrez et al. (2007) that demonstrate the expression of GATA-1 in circulating inflammatory monocytes and its role in the maturation and survival of DCs [46] and given our in silico predictions that the Dectin-2rs7134303_G allele might create a binding site for this transcription factor, we might hypothesize that this SNP could regulate Dectin-2 expression and, consequently, to promote Dectin-2-induced NFkB activation and production of cytokines such as TNF and IL1RA, the key molecules in the ethiopathogenesis of RA. However, we cannot rule out the possibility that Dectin-2rs7134303 and Dectin-2rs4264222 polymorphisms may be in linkage disequilibrium with unidentified susceptible variants that are responsible for the significant associations observed.

MCP-1/CCL2 is a modulator of monocyte/macrophage recruitment to the site of inflammation and subsequent T-cell activation [8] that acts through the CC chemokine receptor 2 (CCR2) [47]. MCP-1 is involved in the promotion of leukocyte infiltration during the inflammation process and it has been found highly expressed in the synovial fluid of RA patients [7]. Several polymorphisms have been identified in the gene encoding MCP-1 and some of them have previously been studied in relation to RA susceptibility [48],[49]. So far, most of these studies have been focused only on a potentially functional SNP, the MCP-1rs1024611 (-2518A/G) [24], [25], but no association with RA was found in Caucasian [48] or non-Caucasian populations [49], [50]. Conversely to these studies, we found an overall protective effect of the MCP-1rs1024611, MCP-1rs13900 and MCP-1rs4586 polymorphisms on the development of RA. However, given that the effect of these SNPs was restricted to women and given that these previous studies were based on a small population size, these studies might not have had sufficient statistical power to detect this association. Our predictive functional analysis also revealed that MCP-1rs1024611 disrupts the binding site for GATA-1 and GATA-2, which support our hypothesis suggesting a central role of these transcription factors mediating monocyte and DC maturation by controlling MCP-1 and Dectin-2 expression.

Similarly to other case-control studies, this study has several potential limitations. First, the statistical power was limited when a gender-stratified analysis was performed. To overcome this problem, we reported dominant rather than only co-dominant model. This limitation has also to be considered when analyzing SNP-SNP interactions as the combined analysis reduced some genotype groups to a few individuals. Finally, another limitation of this study is that no in vivo mechanisms were clarified in relation to Dectin-2, DC-SIGN and MCP-1 variants.

Conclusions

In conclusion, this study represents a preliminary step to account for reported gender differences in RA incidence by demonstrating that genetic polymorphisms within immune-related genes may have different effects in women and men, hence determining the way their immune systems respond to autoimmune stimuli. Nonetheless, further studies using independent populations are warranted to validate our findings.

Supporting Information

Table S1.

Genotype frequencies and risk estimates of polymorphic loci in genes related to the macrophage/dendritic cell-induced immune response. 1Models adjusted for age and gender. 2Models adjusted for age. 3p value for testing of effect modification by gender was calculated utilizing an interaction term of gender and genetic polymorphism assuming a co-dominant model of inheritance. Results in bold show p<0.05. Abbreviations: OR, odds ratio; CI, confidence interval. All analyzed SNPs were in HWE in the control group with the exception of Dectin-1rs16910631 (p>0.01). This SNP was excluded from the analysis.

https://doi.org/10.1371/journal.pone.0072732.s001

(DOCX)

Table S2.

Demographic and clinical characteristics of the RA population (Phase 2). Data are means ± standard deviation. Abbreviations: RF, rheumatoid factor; Anti-CCP: anti-cyclic citrullinated peptide antibodies; DAS28, disease activity score; DMARDs, disease-modifying antirheumatic drugs. * Anti-CCP value was available only in 314 patients (254 women and 60 men).

https://doi.org/10.1371/journal.pone.0072732.s002

(DOCX)

Table S3.

Demographic and clinical characteristics of the pooled population. Data are means ± standard deviation. Abbreviations: RF, rheumatoid factor; Anti-CCP: anti-cyclic citrullinated peptide antibodies; DAS28, disease activity score; DMARDs, disease-modifying antirheumatic drugs. † Rheumatoid factor and anti-CCP values were available for 1.100 (907 women and 193 men) and 709 patients (582 women and 127 men), respectively. *Of those 1.212 RA patients (986 women and 226 men) were genotyped.

https://doi.org/10.1371/journal.pone.0072732.s003

(DOCX)

Acknowledgments

We thank all participants who have agreed to participate in this study. Authors also thank Dr. Marta Alarcón Riquelme for her suggestions that substantially improved the manuscript and María Dolores Casares, Ángeles Molina and Carmen Olóriz for sample collection and technical support.

Author Contributions

Conceived and designed the study: RC AF KH JS. Responsible for genotyping: JS LC CBL. Did the statistical data analysis: JS MER. Coordinated the sample collection: HC AE IF JSC MAF AG LR AGU TV AG EPP EC JEF JS. Were involved in the records review and data acquisition: MJSP HC. Drafted the manuscript: RC AF KH JS.

References

  1. 1. Scrivo R, Di Franco M, Spadaro A, Valesini G (2007) The immunology of rheumatoid arthritis. Ann N Y Acad Sci 1108: 312–322.
  2. 2. Miossec P (2004) An update on the cytokine network in rheumatoid arthritis. Curr Opin Rheumatol 16: 218–222.
  3. 3. Burmester GR, Stuhlmuller B, Keyszer G, Kinne RW (1997) Mononuclear phagocytes and rheumatoid synovitis. Mastermind or workhorse in arthritis? Arthritis Rheum 40: 5–18.
  4. 4. Kinne RW, Brauer R, Stuhlmuller B, Palombo-Kinne E, Burmester GR (2000) Macrophages in rheumatoid arthritis. Arthritis Res 2: 189–202.
  5. 5. Gringhuis SI, den Dunnen J, Litjens M, van der Vlist M, Wevers B, et al. (2009) Dectin-1 directs T helper cell differentiation by controlling noncanonical NF-kappaB activation through Raf-1 and Syk. Nat Immunol 10: 203–213.
  6. 6. van Lent PL, Figdor CG, Barrera P, van Ginkel K, Sloetjes A, et al. (2003) Expression of the dendritic cell-associated C-type lectin DC-SIGN by inflammatory matrix metalloproteinase-producing macrophages in rheumatoid arthritis synovium and interaction with intercellular adhesion molecule 3-positive T cells. Arthritis Rheum 48: 360–369.
  7. 7. Akahoshi T, Wada C, Endo H, Hirota K, Hosaka S, et al. (1993) Expression of monocyte chemotactic and activating factor in rheumatoid arthritis. Regulation of its production in synovial cells by interleukin-1 and tumor necrosis factor. Arthritis Rheum 36: 762–771.
  8. 8. Daly C, Rollins BJ (2003) Monocyte chemoattractant protein-1 (CCL2) in inflammatory disease and adaptive immunity: therapeutic opportunities and controversies. Microcirculation 10: 247–257.
  9. 9. LeibundGut-Landmann S, Gross O, Robinson MJ, Osorio F, Slack EC, et al. (2007) Syk- and CARD9-dependent coupling of innate immunity to the induction of T helper cells that produce interleukin 17. Nat Immunol 8: 630–638.
  10. 10. Myasoedova E, Crowson CS, Kremers HM, Therneau TM, Gabriel SE (2010) Is the incidence of rheumatoid arthritis rising?: results from Olmsted County, Minnesota, 1955–2007. Arthritis Rheum 62: 1576–1582.
  11. 11. Iikuni N, Sato E, Hoshi M, Inoue E, Taniguchi A, et al. (2009) The influence of sex on patients with rheumatoid arthritis in a large observational cohort. J Rheumatol 36: 508–511.
  12. 12. Sokka T, Toloza S, Cutolo M, Kautiainen H, Makinen H, et al. (2009) Women, men, and rheumatoid arthritis: analyses of disease activity, disease characteristics, and treatments in the QUEST-RA study. Arthritis Res Ther 11: R7.
  13. 13. Muller-Ladner U, Pap T, Gay RE, Neidhart M, Gay S (2005) Mechanisms of disease: the molecular and cellular basis of joint destruction in rheumatoid arthritis. Nat Clin Pract Rheumatol 1: 102–110.
  14. 14. Karouzakis E, Neidhart M, Gay RE, Gay S (2006) Molecular and cellular basis of rheumatoid joint destruction. Immunol Lett 106: 8–13.
  15. 15. Ackerman LS (2006) Sex hormones and the genesis of autoimmunity. Arch Dermatol 142: 371–376.
  16. 16. Libert C, Dejager L, Pinheiro I (2010) The X chromosome in immune functions: when a chromosome makes the difference. Nat Rev Immunol 10: 594–604.
  17. 17. Bouman A, Schipper M, Heineman MJ, Faas MM (2004) Gender difference in the non-specific and specific immune response in humans. Am J Reprod Immunol 52: 19–26.
  18. 18. Talal N (1992) Sjogren's syndrome: historical overview and clinical spectrum of disease. Rheum Dis Clin North Am 18: 507–515.
  19. 19. Rose NR (1997) Autoimmune diseases: tracing the shared threads. Hosp Pract (Minneap) 32: 147–154.
  20. 20. Gutierrez-Roelens I, Lauwerys BR (2008) Genetic susceptibility to autoimmune disorders: clues from gene association and gene expression studies. Curr Mol Med 8: 551–561.
  21. 21. Julia A, Ballina J, Canete JD, Balsa A, Tornero-Molina J, et al. (2008) Genome-wide association study of rheumatoid arthritis in the Spanish population: KLF12 as a risk locus for rheumatoid arthritis susceptibility. Arthritis Rheum 58: 2275–2286.
  22. 22. Padyukov L, Hytonen AM, Smolnikova M, Hahn-Zoric M, Nilsson N, et al. (2004) Polymorphism in promoter region of IL10 gene is associated with rheumatoid arthritis in women. J Rheumatol 31: 422–425.
  23. 23. Hussein YM, Mohamed RH, Pasha HF, El-Shahawy EE, Alzahrani SS (2011) Association of tumor necrosis factor alpha and its receptor polymorphisms with rheumatoid arthritis in female patients. Cell Immunol 271: 192–196.
  24. 24. Rovin BH, Lu L, Saxena R (1999) A novel polymorphism in the MCP-1 gene regulatory region that influences MCP-1 expression. Biochem Biophys Res Commun 259: 344–348.
  25. 25. McDermott DH, Yang Q, Kathiresan S, Cupples LA, Massaro JM, et al. (2005) CCL2 polymorphisms are associated with serum monocyte chemoattractant protein-1 levels and myocardial infarction in the Framingham Heart Study. Circulation 112: 1113–1120.
  26. 26. Sakuntabhai A, Turbpaiboon C, Casademont I, Chuansumrit A, Lowhnoo T, et al. (2005) A variant in the CD209 promoter is associated with severity of dengue disease. Nat Genet 37: 507–513.
  27. 27. Stephens M, Scheet P (2005) Accounting for decay of linkage disequilibrium in haplotype inference and missing-data imputation. Am J Hum Genet 76: 449–462.
  28. 28. Yuan HY, Chiou JJ, Tseng WH, Liu CH, Liu CK, et al. (2006) FASTSNP: an always up-to-date and extendable service for SNP function analysis and prioritization. Nucleic Acids Res 34: W635–641.
  29. 29. Moore JH (2004) Computational analysis of gene-gene interactions using multifactor dimensionality reduction. Expert Rev Mol Diagn 4: 795–803.
  30. 30. Bowes J, Barton A (2008) Recent advances in the genetics of RA susceptibility. Rheumatology (Oxford) 47: 399–402.
  31. 31. Worthington J (2005) Investigating the genetic basis of susceptibility to rheumatoid arthritis. J Autoimmun 25 Suppl: 16–20
  32. 32. Raychaudhuri S, Thomson BP, Remmers EF, Eyre S, Hinks A, et al. (2009) Genetic variants at CD28, PRDM1 and CD2/CD58 are associated with rheumatoid arthritis risk. Nat Genet 41: 1313–1318.
  33. 33. Gabriel SE (2001) The epidemiology of rheumatoid arthritis. Rheum Dis Clin North Am 27: 269–281.
  34. 34. Hussein YM, El Tarhouny SA, Mohamed RH, El-Shal AS, Abul-Saoud AM, et al.. (2011) Association of CD4 enhancer gene polymorphisms with rheumatoid arthritis in Egyptian female patients. Rheumatol Int.
  35. 35. Hussein Y, El-Tarhouny S, Mohamed R, Pasha H, Abul-Saoud A (2011) Association of interleukin-4 receptor gene polymorphisms with rheumatoid arthritis in Egyptian female patients. Joint Bone Spine.
  36. 36. Plantinga TS, Fransen J, Takahashi N, Stienstra R, van Riel PL, et al. (2010) Functional consequences of DECTIN-1 early stop codon polymorphism Y238X in rheumatoid arthritis. Arthritis Res Ther 12: R26.
  37. 37. Dieguez-Gonzalez R, Akar S, Calaza M, Gonzalez-Alvaro I, Fernandez-Gutierrez B, et al. (2009) Lack of association with rheumatoid arthritis of selected polymorphisms in 4 candidate genes: CFH, CD209, eotaxin-3, and MHC2TA. J Rheumatol 36: 1590–1595.
  38. 38. Taylor PR, Reid DM, Heinsbroek SE, Brown GD, Gordon S, et al. (2005) Dectin-2 is predominantly myeloid restricted and exhibits unique activation-dependent expression on maturing inflammatory monocytes elicited in vivo. Eur J Immunol 35: 2163–2174.
  39. 39. Aragane Y, Maeda A, Schwarz A, Tezuka T, Ariizumi K, et al. (2003) Involvement of dectin-2 in ultraviolet radiation-induced tolerance. J Immunol 171: 3801–3807.
  40. 40. Carter RW, Thompson C, Reid DM, Wong SY, Tough DF (2006) Induction of CD8+ T cell responses through targeting of antigen to Dectin-2. Cell Immunol 239: 87–91.
  41. 41. Robinson MJ, Osorio F, Rosas M, Freitas RP, Schweighoffer E, et al. (2009) Dectin-2 is a Syk-coupled pattern recognition receptor crucial for Th17 responses to fungal infection. J Exp Med 206: 2037–2051.
  42. 42. Sato K, Yang XL, Yudate T, Chung JS, Wu J, et al. (2006) Dectin-2 is a pattern recognition receptor for fungi that couples with the Fc receptor gamma chain to induce innate immune responses. J Biol Chem 281: 38854–38866.
  43. 43. Gringhuis SI, Wevers BA, Kaptein TM, van Capel TM, Theelen B, et al. (2011) Selective C-Rel activation via Malt1 controls anti-fungal T(H)-17 immunity by dectin-1 and dectin-2. PLoS Pathog 7: e1001259.
  44. 44. Miossec P (2007) Interleukin-17 in fashion, at last: ten years after its description, its cellular source has been identified. Arthritis Rheum 56: 2111–2115.
  45. 45. Makarov SS (2001) NF-kappa B in rheumatoid arthritis: a pivotal regulator of inflammation, hyperplasia, and tissue destruction. Arthritis Res 3: 200–206.
  46. 46. Gutierrez L, Nikolic T, van Dijk TB, Hammad H, Vos N, et al. (2007) Gata1 regulates dendritic-cell development and survival. Blood 110: 1933–1941.
  47. 47. Rollins BJ (1997) Chemokines. Blood 90: 909–928.
  48. 48. Gonzalez-Escribano MF, Torres B, Aguilar F, Rodriguez R, Garcia A, et al. (2003) MCP-1 promoter polymorphism in Spanish patients with rheumatoid arthritis. Hum Immunol 64: 741–744.
  49. 49. Lee YH, Kim HJ, Rho YH, Choi SJ, Ji JD, et al. (2003) Functional polymorphisms in matrix metalloproteinase-1 and monocyte chemoattractant protein-1 and rheumatoid arthritis. Scand J Rheumatol 32: 235–239.
  50. 50. Hwang SY, Cho ML, Park B, Kim JY, Kim YH, et al. (2002) Allelic frequency of the MCP-1 promoter -2518 polymorphism in the Korean population and in Korean patients with rheumatoid arthritis, systemic lupus erythematosus and adult-onset Still's disease. Eur J Immunogenet 29: 413–416.