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

Dectin-1 and DC-SIGN Polymorphisms Associated with Invasive Pulmonary Aspergillosis Infection

  • Juan Sainz ,

    juan.sainz@genyo.es

    Affiliations Genomic Oncology Area, Genyo (Pfizer-University of Granada-Andalusian Government Centre for Genomics and Oncological Research), Granada, Spain, Molecular Genetic Epidemiology Department, German Cancer Research Center, Heidelberg, Germany, Haematology Department, University Hospital Virgen de las Nieves, Granada, Spain

  • Carmen Belén Lupiáñez,

    Affiliations Genomic Oncology Area, Genyo (Pfizer-University of Granada-Andalusian Government Centre for Genomics and Oncological Research), Granada, Spain, Haematology Department, University Hospital Virgen de las Nieves, Granada, Spain

  • Juana Segura-Catena,

    Affiliations Genomic Oncology Area, Genyo (Pfizer-University of Granada-Andalusian Government Centre for Genomics and Oncological Research), Granada, Spain, Haematology Department, University Hospital Virgen de las Nieves, Granada, Spain

  • Lourdes Vazquez,

    Affiliation Hematology Department, University Hospital of Salamanca, Salamanca, Spain

  • Rafael Ríos,

    Affiliations Genomic Oncology Area, Genyo (Pfizer-University of Granada-Andalusian Government Centre for Genomics and Oncological Research), Granada, Spain, Haematology Department, University Hospital Virgen de las Nieves, Granada, Spain

  • Salvador Oyonarte,

    Affiliation Blood Transfusion Regional Centre and Sectorial Tissue Bank, Granada, Spain

  • Kari Hemminki,

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

  • Asta Försti,

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

  • Manuel Jurado

    Affiliations Genomic Oncology Area, Genyo (Pfizer-University of Granada-Andalusian Government Centre for Genomics and Oncological Research), Granada, Spain, Haematology Department, University Hospital Virgen de las Nieves, Granada, Spain

Abstract

The recognition of pathogen-derived structures by C-type lectins and the chemotactic activity mediated by the CCL2/CCR2 axis are critical steps in determining the host immune response to fungi. The present study was designed to investigate whether the presence of single nucleotide polymorphisms (SNPs) within DC-SIGN, Dectin-1, Dectin-2, CCL2 and CCR2 genes influence the risk of developing Invasive Pulmonary Aspergillosis (IPA). Twenty-seven SNPs were selected using a hybrid functional/tagging approach and genotyped in 182 haematological patients, fifty-seven of them diagnosed with proven or probable IPA according to the 2008 EORTC/MSG criteria. Association analysis revealed that carriers of the Dectin-1rs3901533 T/T and Dectin-1rs7309123 G/G genotypes and DC-SIGNrs4804800 G, DC-SIGNrs11465384 T, DC-SIGN7248637 A and DC-SIGN7252229 C alleles had a significantly increased risk of IPA infection (OR = 5.59 95%CI 1.37–22.77; OR = 4.91 95%CI 1.52–15.89; OR = 2.75 95%CI 1.27–5.95; OR = 2.70 95%CI 1.24–5.90; OR = 2.39 95%CI 1.09–5.22 and OR = 2.05 95%CI 1.00–4.22, respectively). There was also a significantly increased frequency of galactomannan positivity among patients carrying the Dectin-1rs3901533_T allele and Dectin-1rs7309123_G/G genotype. In addition, healthy individuals with this latter genotype showed a significantly decreased level of Dectin-1 mRNA expression compared to C-allele carriers, suggesting a role of the Dectin-1rs7309123 polymorphism in determining the levels of Dectin-1 and, consequently, the level of susceptibility to IPA infection. SNP-SNP interaction (epistasis) analysis revealed significant interactions models including SNPs in Dectin-1, Dectin-2, CCL2 and CCR2 genes, with synergistic genetic effects. Although these results need to be further validated in larger cohorts, they suggest that Dectin-1, DC-SIGN, Dectin-2, CCL2 and CCR2 genetic variants influence the risk of IPA infection and might be useful in developing a risk-adapted prophylaxis.

Introduction

Invasive Pulmonary Aspergillosis (IPA) is a life-threatening infection caused mainly by Aspergillus fumigatus, an opportunistic fungal pathogen that frequently colonizes respiratory tracts and rapidly spreads to blood vessels and tissues [1]. The incidence of IPA infection has increased in the last years due to the use of immunosuppressive and immunomodulatory drugs and is still causing significant morbidity and mortality worldwide, especially in immunosuppressed and hematopoietic stem cell transplantation (HSCT) patients [2].

Early recognition of A. fumigatus by myeloid leukocytes is a crucial for down-stream immune response and conidia clearance and, therefore, in the development and progression of IPA infection [3]. Myeloid leukocytes are activated by pattern-recognition receptors (PRRs), which directly recognize fungal cell wall structures and pathogen-associated molecular patterns (PAMPs) [4]. Among the PRRs expressed in neutrophils, pulmonary macrophages and dendritic cells, C-type lectin family members such as DC-SIGN (Dendritic Cell-Specific Intercellular adhesion molecule-3-Grabbing Non-integrin), Dectin-1 (Dendritic cell-associated C-type lectin-1) and Dectin-2 (Dendritic cell-associated C-type lectin-2) seem to be important in the effectiveness of innate immune response against A. fumigatus [5], [6]. Recent studies have reported that the C-type lectins mediate A. Fumigatus recognition, capture and internalization [7] and that their binding to fungal cell wall carbohydrates is highly specific and selective [8]. DC-SIGN and Dectin-1 recognize β-(1-3)-glucans and galactomannans (GM) in the cell wall of A. fumigatus [7], [9] while Dectin-2 interacts with α-mannans [10]. However, numerous studies have shown that C-type lectins are not only involved in the recognition of fungal pathogens but also in the induction of anti-fungal Th1 and Th17 immune responses [5], [11].

Although the mechanisms underlying DC-SIGN-mediated immune response are still highly speculative, it has been suggested that DC-SIGN promotes dendritic cell (DC) migration and T-cell activation through the ICAM-3 binding [12], [13] and modulates TLR signalling by targeting Raf-1, which regulates NFκB p65 activation and, consequently the production of pro-inflammatory cytokines [14]. Likewise, Dectin-1 and Dectin-2 induce Syk-dependent canonical and noncanonical NFκB pathways [15], [16] promoting the production of some pro-inflammatory cytokines (IL1α, IL1β, IL12 and IL8), chemokines (CCL2/MCP-1, CCL3/MIP-1, CXCL2/MIP-2 and CXCL10) [17] and the IL17-mediated neutrophil recruitment [18]. It has also been described that Dectin-1 and Dectin-2 are able to work in collaboration with TLRs (mainly TLR2, TLR4 and TLR6) modulating immune responses through the production of IL6, IL12p70 and TNF [3], [16]. These findings support the hypothesis that C-type lectins, cytokines or chemokines are important mediators during infection with A. fumigatus.

Several polymorphisms in human PRRs [19][24] as well as in some cytokines [25][27], chemokines [28] and their receptors [29], [30] have hitherto been associated with an increased risk of invasive aspergillosis in susceptible hosts. Based on these observations, the objective of the present study was to investigate the role of tagging and potentially functional single-nucleotide polymorphisms (SNPs) located within the DC-SIGN, Dectin-1, Dectin-2, MCP-1/CCL2 and CCR2 genes on IPA susceptibility.

Materials and Methods

Study population and clinical diagnosis of IPA infection

All participants enrolled were Caucasian and recruited in the University Hospital Virgen de las Nieves (Granada, Spain) or in the University Hospital of Salamanca (Salamanca, Spain). All determinations and genetic analyses in hematological patients were performed with fully informed written consent, and anonymity of the data was guaranteed. The study protocol was approved by the Ethical Review Committee of Virgen de las Nieves University Hospital, Granada, Spain. The population included 182 hematological patients recruited between January 2004 and January 2011. All hematological patients in this study received a prolonged chemotherapy treatment or underwent HSCT and were therefore considered susceptible to develop IPA infection. Demographic information and clinical data were obtained by detailed review of hospital records. Data were gathered on: site of infection; host factor criteria (severe neutropenia for >10 days, persistent fever for >96 h refractory to appropriate broad-spectrum antibacterial treatment, signs and symptoms indicating graft-versus-host disease [GVHD], corticoid therapy [>0.3 mg/kg per day], and invasive fungal infection during a previous episode of neutropenia), microbiological criteria (positive result for Aspergillus antigen in 2 consecutive blood samples when considering an index of 0.5 or in only 1 sample when the index was higher than 0.8), and clinical criteria of lower respiratory tract infection [major criteria: any of the following new infiltrates on computed tomography (CT) imaging: halo sign, air-crescent sign, or cavity within area of consolidation; minor criteria: cough, thoracic pain, hemoptysis, pathologic pulmonary sound, and radiological evidence suggestive of invasive infection]. Laboratory data were also recorded. Proven and probable IPA was diagnosed based on the updated criteria (2008) reported by the European Organization for Research and Treatment of Cancer/Invasive Fungal Infections Cooperative Group (EORTC/IFICG) [31].

Quantification of serum galactomannan antigen

Serum GM detection has been shown to be a useful test for the early diagnosis and follow-up of IPA and is now included in IPA diagnosis criteria [31]. In the present study, serum GM antigen was determined twice weekly during the hospital stay and at each outpatient visit until the end of their immunosuppressant or chemotherapeutic treatment. Serum GM concentrations were determined by Platelia Aspergillus ELISA kit (Bio-Rad, Marnes-la-Coquette, France) according to the manufacturer's instructions. This commercial kit has proven to offer good sensitivity to detect GM [32], and GM concentration was found to correlate with the fungal tissue burden [33], [34]. A test sample was classified as positive when the optical density ratio was ≥0.5 in two consecutive positive samples or >0.8 in a unique serum sample. A careful review of concomitant treatments (piperacillin-tazobactam or amoxicillin-clavulonic acid) in each patient was necessary to detect possible false-positive GM determinations. Likewise, tests were performed on the same day to avoid sample contamination and to ensure accuracy of results.

SNP selection and genotyping

Twenty-seven tagging/functional SNPs within DC-SIGN, Dectin-1, Dectin-2, CCL2 and CCR2 were selected to genotype the entire panel of individuals (Table 1). 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, CCL2 and CCR2 genes on IPA risk. Tagging SNPs were selected using Haploview Tagger program (http://www.broad.mit.edu/mpg/haploview/; http://www.broad.mit.edu/mpg/tagger/) and a pairwise tagging with a minimum r2 of 0.8. In this selection we forced the inclusion of the DC-SIGNrs4804803, CCL2rs1024610 and CCL2rs1024611 polymorphisms as their functionality has been suggested [35][37]. Genomic DNA was extracted from peripheral blood mononuclear cells (PBMCs) Qiagen Mini Kit (Qiagen, Valencia, CA, USA). Genotyping of DC-SIGN, Dectin-1, Dectin-2, CCL2 and CCR2 polymorphisms was carried out using KASPar assays (KBiosciences, Hoddesdon, Hertfordshire, UK) in a 384-well plate format (Applied Biosystems, Foster City, CA, USA) where hematological patient samples were randomly distributed. KASPar reactions were performed using KASPar assay mix (containing probes) and KASPar kit containing 2× Reaction Mix and MgCl2 (50 mM). PCR conditions were: denaturation at 94°C for 15 min, 10 cycles of denaturation at 94°C for 10 sec, annealing at 57°C for 5 sec and elongation at 72°C for 10 sec and 20 cycles of denaturation at 94°C for 10 sec, annealing at 57°C for 20 sec and elongation at 72°C for 40 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, about 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.5%. Call rates for all SNPs were ≥97.8% with the exception of the Dectin-1rs11053599 SNP with a call rate of 94.5%.

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

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

Statistical analysis

The Hardy-Weinberg Equilibrium (HWE) tests were performed in the control group (non-IPA patients) by a standard observed-expected chi-square (χ2) test at each polymorphic site (http://ihg2.helmholtz-muenchen.de/cgi-bin/hw/hwa1.pl). Unconditional logistic regression was used to assess the main effects of the genetic polymorphisms on IPA risk using co-dominant, dominant and recessive inheritance models. For each SNP, the more common allele in the control group was assigned as the reference category. All analyses were adjusted for age, gender, hematological malignancy and established risk factors for IPA infection (HSCT, neutropenia, GVHD and corticoid therapy use) and were conducted using the statistical software SSPS (version 14.0, SPSS Inc., Chicago, USA). All tests were considered to be statistically significant with a p value of <0.05.

Adjustment for multiple testing was carried out following a conservative threshold for statistical significance, based on a revised version of the Bonferroni method: we calculated for each gene the “number of effective independent variables” (Meff) by use of the SNP Spectral Decomposition approach (http://gump.qimr.edu.au/general/daleN/SNPSpDsuperlite/) [38]. We obtained a study-wise Meff value by adding up the gene Meff's.

SNPtool (http://www.dkfz.de/de/molgen_epidemiology/tools/SNPtool.html) [39] and the Haploview v4.2 software were used for LD blocks reconstruction and haplotype association statistics. Block structures were determined according to the method of Gabriel et al. [40].

Functional prediction of associated SNPs

We used the Web-based tool FastSNP [41] available at http://fastsnp.ibms.sinica.edu.tw for predicting the functional significance of the SNPs associated with IPA infection. FastSNP utilizes information from another online program PolyPhen (http://www.bork.embl-heidelberg.de/PolyPhen/) and from four different web resources (TFSearch, ESEfinder, Rescue-ESE and FAS-ESS) to determine whether SNPs are located at exonic splicing regulatory sites, cause a non-conservative amino acid change or whether they alter the transcription factor-binding site of a gene (for instance, acting as intronic enhancer). The score was given by this tool on the basis of levels of risk with a ranking of 0 (no effect), 1 (very low), 2 (low), 3 (medium), 4 (high), or 5 (very high).

RNA extraction, reverse transcription, and qPCR

Whole blood samples from 21 healthy donors were collected into PAXGENE RNA tubes and stored at −80°C until use. Total RNAs were extracted using a PAXGENE Blood RNA Isolation Kit (PreAnalytiX) and reverse transcribed to cDNA using QuantiTect Reverse Transcription Kit (catalog number: 205311; Qiagen). Real-time quantitative PCR was carried out using an ABI PRISM® 7500 HT Sequence Detection System (Applied Biosystems) according to manufacturer's instructions. Briefly, 2 µl of the cDNA were loaded in a PCR reaction containing 12.5 µl of 2× QuantiTect SYBR Green PCR Master Mix with an appropriate concentration of MgCl2, 2.5 µl of primers (Hs_CLEC7A_1_SG QuantiTect Primer Assay, catalog number: QT00024059; Geneglobe, Qiagen) and 8 µl of RNase-free water. PCR cycling conditions were as follows: 95°C for 15 minutes, followed by 40 cycles of denaturation at 95°C for 15 seconds combined with annealing at 60°C for 30 seconds and extension at 72°C for 30 seconds. All samples were run in duplicate. Relative quantification of Dectin-1 mRNA expression was calculated with the 2−ΔΔCt method. We obtained the fold changes in gene expression normalized to an internal control gene (GAPDH, Hs_GAPDH_2_SG QuantiTect Primer Assay, Catalog number: QT01192646; Qiagen) and relative to one calibrator (FirstChoice®Human Brain Reference RNA, catalog number: 6050; Applied Biosystems).

SNP-SNP interaction analysis

We also analyzed high-order interactions between SNPs using the multifactor dimensionality reduction (MDR) constructive induction algorithm. A detailed explanation on the MDR method has been described elsewhere [42], [43]. SNPs in LD (cut-off of r2∼0.8) were excluded from the MDR analysis. Cross-validation and permutation testing were used to identify the best models. All possible two-, three- and four-way SNP interactions were tested using 100-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 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). MDR results were considered statistically significant at the 0.05 level. Finally, interactions were visualized by constructing an interaction dendrogram according to the method described by Moore et al [43]. MDR software and MDR permutation testing module are open-source and freely available from http://www.epistasis.org.

Results

Population characteristics are described in Table 2. Compared to non-IPA, IPA cases were more likely to have cough and pathologic pulmonary sound (p<0.001 and 0.011, respectively) and presented more often pathological chest radiographies and CT scans (p<0.001). Established risk factors for IPA infection such as HSCT, neutropenia, GVHD and corticoid therapy use were homogenously distributed between IPA and non-IPA patient groups.

thumbnail
Table 2. Demographic and clinical data of IPA and non-IPA patients.

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

Fifty-seven patients were diagnosed with proven or probable IPA infection according to the revised EORTC/MSG criteria (2008). The association of risk of IPA infection with the individual 27 SNPs in the C-type lectin and chemokine genes is shown in the Table S1. All analyzed polymorphisms fulfilled Hardy-Weinberg expectations for the control group (non-IPA patients).

Several polymorphisms were found to be associated with IPA infection (Table 3). The Dectin-1rs3901533, Dectin-1rs7309123 and DC-SIGNrs4804800 polymorphisms showed the strongest association with risk of IPA infection (ORT/T = 5.59, 95%CI1.37–22.77; ORG/G = 4.91, 95%CI1.52–15.89 and ORA/G+G/G = 2.75, 95%CI 1.27–5.95, respectively). Additionally, carriers of the DC-SIGNrs11465384_T, DC-SIGN7248637_A and DC-SIGN7252229_C alleles had a significantly increased risk of IPA infection (ORC/T+T/T = 2.70, 95%CI 1.24–5.90; ORG/A+A/A = 2.39, 95%CI 1.09–5.22; ORG/C+C/C = 2.05, 95%CI 1.00–4.22, respectively) whereas patients carrying the DC-SIGNrs2287886_A allele showed a tendency to be associated with a decreased risk of IPA infection (per-allele OR = 0.60, 95%CI 0.34–1.05) (Table 3 and Table S1). After correction for multiple testing using SNPSpD (number of independent marker loci, 21; p = 0.05/21 = 0.002), none of the SNPs retained significance although Dectin-1rs7309123 showed significance when the carrier status was analyzed (CC+CG vs. GG, P = 0.001; Table 3).

thumbnail
Table 3. DC-SIGN and Dectin-1 polymorphisms associated with susceptibility to Invasive Pulmonary Aspergillosis.

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

In order to assess the degree to which the selected SNPs and the positivity of GM correlated, we also performed lineal regression analysis. In the whole population, 3784 assays were performed in duplicate and 531 were considered as positive (14.03%). We found a significantly higher percentage of positive GM among patients carrying the Dectin-1rs3901533_T allele and among those patients bearing the Dectin-1rs7309123_G/G genotype suggesting a role of these polymorphisms in determining a defective recognition and clearance of Aspergillus conidia (Figure 1). No correlation was observed for the polymorphisms within DC-SIGN that were associated with the infection.

thumbnail
Figure 1. Distribution of positive galactomannan percentage by DC-SIGN and Dectin-1 genotypes.

https://doi.org/10.1371/journal.pone.0032273.g001

We subsequently assessed whether SNPs associated with IPA infection showed capacity to change putative transcription factor binding sites using FastSNP. The predictive functional analysis suggested an intronic enhancer function for the Dectin-1rs7309123 SNP due to its location in transcription factor Cdxa (caudal type homeo-box transcription factor 1) binding site (GAAAGAC; score 1–2). These data suggest a central role of the rs7309123 in the susceptibility to IPA infection. None of the remaining SNPs associated with IPA infection was predicted to affect transcription factor binding sites or splicing or to introduce a damaging amino acid change.

To explore the potential consequences of Dectin-1rs7309123 SNP, Dectin-1 mRNA expression was measured in 21 healthy donors and was correlated with genotypes. Our results showed that subjects harbouring GG genotype showed a significantly decreased level of Dectin-1 mRNA expression compared to C-allele carriers (p<0.001; Figure 2). These results further supported our hypothesis suggesting that Dectin-1rs7309123_G allele may disrupt binding sites for potential transcription factors.

thumbnail
Figure 2. Correlation between Dectin-1rs7903123 genotype and expression of Dectin-1 as measured by quantitative PCR in total RNA from healthy donors (n = 21).

Dectin-1 mRNA levels were normalized for GAPDH mRNA levels. FirstChoice®Human Brain Reference RNA was used as Calibrator (Ambion; Catalog number AM6050).

https://doi.org/10.1371/journal.pone.0032273.g002

The MDR analysis of the SNPs tested revealed statistically significant interactions. Two SNPs (rs1024611 and rs4264222) were excluded from the MDR analysis as they were in LD with other SNPs using a cut-off value of r2∼0.80 in our sample set. The best interaction models selected by the TuRF filter algorithm along with its testing accuracy and cross-validation consistency are shown in Table 4. The overall best model with the highest cross-validation consistency (CVC) consisted of a model that included the CCL2rs4586, Dectin-1rs3901533, CCR2rs3918358 and Dectin-2rs7134303 SNPs. This model had a significant testing accuracy of 0.7735 (permutation p<0.001) and a cross-validation consistency of 100/100. Of note is that two SNPs showing genetic interaction in this model were not significantly associated with an increased risk of IPA infection in the univariate analysis (CCR2rs3918358 and Dectin-2rs7134303). The best two-locus model consisted of the Dectin-1rs3901533 and DC-SIGNrs4804800 SNPs, two variants that showed also association in the single-locus analysis. This model had a testing accuracy of 0.6409 and a cross-validation consistency of 76/100. This model was not any more significant after 1.000-fold permutation testing. However, the entropy based information gain calculated for this pair of SNPs indicated strong synergy, which may be interpreted as the two SNPs acting together to increase the risk of IPA infection. The best three-locus model included the CCL2rs4586, Dectin-1rs3901533 and Dectin-2rs7134303 SNPs. In this model, testing accuracy raised to 0.7085 (permutation p = 0.025) whereas the cross-validation consistency was of 68/100. Figure 3 illustrates an interaction dendogram that summarizes the estimates of interactions.

thumbnail
Figure 3. Interaction dendrogram generated by the MDR software.

The interaction dendrogram was used to confirm, visualize, and interpret the interaction model. The MDR analysis was performed by using the open-source MDR software package. The colors used depict the degree of synergy, ranging from red (highest information gain) to blue (highest information redundancy). Note that the interaction between Dectin-1 (rs3901533) and DC-SIGN (rs4804800) SNPs showed the highest degree of synergy (gain of information).

https://doi.org/10.1371/journal.pone.0032273.g003

thumbnail
Table 4. Multifactor dimensionality reduction analysis summary.

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

Discussion

The marked differences in susceptibility to IPA infection among hematological patients (with or without allo-HSCT) suggest that the effective immune response against Aspergillus is determined by both environmental and host genetic factors [44], [45]. Studies on genetic polymorphisms in genes coding for components of the innate immunity have supported this hypothesis [19], [23][30], [46].

In this report, we studied the influence of the tagging and potentially functional polymorphisms of DC-SIGN, Dectin-1, Dectin-2, CCL2/MCP-1 and CCR2 in susceptibility to IPA infection in a Spanish population. Polymorphisms in these genes have been reported to influence a number of infectious diseases including HIV-1 [47][49], HTLV-1 [50], CMV [51], Tuberculosis [52][55], HCV [56], [57], HBV [58], Dengue [35] and SARS [59] among others, revealing their potential role in host defense against pathogens.

We found that polymorphisms in Dectin-1 (rs3901533 and rs7309123) and DC-SIGN (rs4804800, rs11465384, rs7248637 and rs7252229) were associated with an increased risk to develop IPA infection, which points towards their critical involvement in the pathogenesis of this invasive fungal infection. The highest risk of IPA infection was found for carriers of the Dectin-1rs3901533_T/T and Dectin-1rs7309123_G/G genotypes and the DC-SIGNrs4804800_G allele. Patients carrying these genotypes/alleles had from 2 to 6 times increased risk of IPA infection. Additionally, patients carrying the DC-SIGNrs11465384_T, DC-SIGN7248637_A and DC-SIGNrs7252229_C alleles showed a 2-fold increased risk in comparison with patients harboring the wild-type allele. Although it was not statistically significant, we also found that the DC-SIGNrs2287886 SNP may be associated with a reduced risk of IPA infection. Interestingly, this latter result was in agreement with our previously reported findings using the former EORTC/MSG classification criteria, 2005 [60]. The apparent effect of these SNPs on IPA susceptibility persisted even after adjustment for age, gender, hematological malignancy and several known risk factors (HSCT, neutropenia, GVHD and corticoid therapy use), indicating that Dectin-1 and DC-SIGN variants contribute independently to the risk of infection.

Another interesting finding of this study was the significantly greater positive GM percentage of patients carrying the Dectin-1rs3901533_T allele than those with the wild-type allele. Additionally, patients harboring the G/G genotype for the Dectin-1rs7903123 SNP showed an increased percentage of positive GM compared to those carrying the C allele (C/C+C/G). No differences were found when positive GM determinations were correlated with DC-SIGN polymorphisms. These data along with the remarkable degree of association of Dectin-1 and DC-SIGN variants with risk of IPA infection provides a compelling evidence for a critical role for these PRRs in immune response to IPA infection. In this regard, numerous studies have shown that Dectin-1 and DC-SIGN are not only involved in the recognition of fungal pathogens but also in the induction of anti-fungal Th1 and Th17 immune responses [5], [11]. Mezger et al. also demonstrated that Dectin-1 is involved in the induction of several pro-inflamatory cytokines, chemokines and immune receptors [18] while Werner et al. showed that Dectin-1 is also regulating Th17-mediated immune response in the lungs [61]. Furthermore, Dennehy and Brown suggested a role of Dectin-1 mediating its own signaling, as well as synergizing with TLRs to trigger NFκB-mediated immune response against fungal pathogens [62].

Although it is now well recognized that SNPs in genes modulating immune response are likely to be determinants of host susceptibility to fungal infections, so far, little is known regarding the biological significance of these variants. In order to shed light into the potential functionality of Dectin-1 (rs3901533 and rs7309123) and DC-SIGN (rs4804800, rs11465384, rs7248637, rs7252229 and rs2287886) variants, we investigated whether they were involved in disruption of a binding site for critical transcription factors that might influence transcription level of these genes. Our predictive analysis showed that the carriage of the C allele for the Dectin-1rs7309123 SNP creates a putative binding site for Cdxa, a relatively unknown transcription factor, which might be involved in the control of Dectin-1 gene expression. To assess whether the Dectin-1rs7309123 polymorphism might be associated with a decreased expression of Dectin-1, we correlated Dectin-1 mRNA expression levels with Dectin-1rs7309123 genotypes. Interestingly, we observed that individuals harbouring the GG genotype showed a relatively lower expression than those carrying the C allele (CC+GC). These results further supported our hypothesis suggesting that Dectin-1rs7309123 SNP may have an effect on the Dectin-1-mediated recognition of Aspergillus conidia and subsequent immune responses. However, this predicted change in transcription activity is only suggestive at this stage and will need further validation using in vitro functional assays.

Several lines of evidence point to the relevance of epistatic effects in the etiology of complex diseases but, up to now, no studies have been carried out to analyze the presence of SNP-SNP interactions in IPA infection. For this reason, we decide to assess interactions among genetic polymorphisms within DC-SIGN, Dectin-1, Dectin-2, CCL2 and CCR2, genes and the risk of IPA infection. The MDR approach used in this study was able to determine two multilocus combinations associated with high risk to develop IPA infection. Of the interactions identified, MDR indicated that the type of interaction in the two significant models was synergistic. These results support the hypothesis that multiple SNP–SNP interactions may play a role in determining the risk of IPA infection. This hypothesis is biologically plausible since the immune system would warrant prevention of fungal infection even when some genetic variants were present.

Recent population-based studies have even led to the identification of several SNPs involved in the early recognition of Aspergillus and associated them with an increased risk for invasive fungal infection. It has previously been suggested that SNPs within C-type lectin genes are associated with fungal infections. Platinga et al. (2009) and Cunha et al. (2010) suggested that patients carrying the Y238X (rs16910526) polymorphism in the Dectin-1 gene were more likely to be colonized with Aspergillus and Candida species [19], [63], compared with those harboring the wild-type allele. However, these results were not replicated in a very recent study [24]. In our study, such findings were neither evidenced even when HSCT and non-HSCT patients were analyzed separately (data not shown). This puzzling finding suggests that, although intronic polymorphisms within Dectin-1 and DC-SIGN are indeed themselves a strong indication that these genes play an important role in the susceptibility to IPA infection, we cannot rule out the possibility that these SNPs are part of a bigger haplotype containing important other genetic variants in the neighboring genes. In any case, because all these population-based studies have been conducted using relatively small cohorts, additional studies in larger set of patients are needed to definitely establish the role of these variants in the susceptibility to invasive fungal infection.

In summary, this study provides evidence of association between Dectin-1 and DC-SIGN polymorphisms and the risk of IPA infection. By the inclusion of functional prediction analyses, the correlation of genotypes with positive GM determinations and Dectin-1 mRNA expression levels, the study strongly supported the role of Dectin-1 gene variants in determining susceptibility to IPA infection. Epistatic analyses also suggested the presence of a gene-gene interaction involving Dectin-1 with CCL2 and CCR2 variants to determine IPA infection. Despite all these evidences, additional studies using larger cohorts will be necessary to confirm the effect of these polymorphisms on the susceptibility to IPA infection.

Supporting Information

Table S1.

Associations of polymorphisms involved in the phagocyte-immune related response against Aspergillus. 1Models adjusted for age, gender, hematological malignancy, HSCT, neutropenia (defined as absolute neutrophil count <500 cells/mm3 for a period of more than 10 days), GVHD and corticoid therapy use (>0.3 mg/Kg/day). ‡ Assuming a recessive model of inheritance. Abbreviations: OR, odds ratio; CI, confidence interval. Differences in samples numbers are due to failures in genotyping.

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

(DOC)

Acknowledgments

We thank all participants in this study and Rosa Cano (laboratory technician) for sample collection and technical support.

Author Contributions

Conceived and designed the experiments: JS MJ. Performed the experiments: JS. Analyzed the data: JS CBL. Contributed reagents/materials/analysis tools: KH AF. Wrote the paper: JS AF MJ. Sample and data collection: JSC CBL LV RR SO.

References

  1. 1. Soubani AO, Chandrasekar PH (2002) The clinical spectrum of pulmonary aspergillosis. Chest 121: 1988–1999.
  2. 2. Marr KA, Carter RA, Boeckh M, Martin P, Corey L (2002) Invasive aspergillosis in allogeneic stem cell transplant recipients: changes in epidemiology and risk factors. Blood 100: 4358–4366.
  3. 3. Luther K, Torosantucci A, Brakhage AA, Heesemann J, Ebel F (2007) Phagocytosis of Aspergillus fumigatus conidia by murine macrophages involves recognition by the dectin-1 beta-glucan receptor and Toll-like receptor 2. Cell Microbiol 9: 368–381.
  4. 4. Aderem A, Underhill DM (1999) Mechanisms of phagocytosis in macrophages. Annu Rev Immunol 17: 593–623.
  5. 5. Geijtenbeek TB, Gringhuis SI (2009) Signalling through C-type lectin receptors: shaping immune responses. Nat Rev Immunol 9: 465–479.
  6. 6. Steele C, Rapaka RR, Metz A, Pop SM, Williams DL, et al. (2005) The beta-glucan receptor dectin-1 recognizes specific morphologies of Aspergillus fumigatus. PLoS Pathog 1: e42.
  7. 7. Serrano-Gomez D, Dominguez-Soto A, Ancochea J, Jimenez-Heffernan JA, Leal JA, et al. (2004) Dendritic cell-specific intercellular adhesion molecule 3-grabbing nonintegrin mediates binding and internalization of Aspergillus fumigatus conidia by dendritic cells and macrophages. J Immunol 173: 5635–5643.
  8. 8. Adams EL, Rice PJ, Graves B, Ensley HE, Yu H, et al. (2008) Differential high-affinity interaction of dectin-1 with natural or synthetic glucans is dependent upon primary structure and is influenced by polymer chain length and side-chain branching. J Pharmacol Exp Ther 325: 115–123.
  9. 9. Feinberg H, Mitchell DA, Drickamer K, Weis WI (2001) Structural basis for selective recognition of oligosaccharides by DC-SIGN and DC-SIGNR. Science 294: 2163–2166.
  10. 10. Saijo S, Ikeda S, Yamabe K, Kakuta S, Ishigame H, et al. (2010) Dectin-2 recognition of alpha-mannans and induction of Th17 cell differentiation is essential for host defense against Candida albicans. Immunity 32: 681–691.
  11. 11. 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.
  12. 12. Geijtenbeek TB, Torensma R, van Vliet SJ, van Duijnhoven GC, Adema GJ, et al. (2000) Identification of DC-SIGN, a novel dendritic cell-specific ICAM-3 receptor that supports primary immune responses. Cell 100: 575–585.
  13. 13. Bozza S, Gaziano R, Spreca A, Bacci A, Montagnoli C, et al. (2002) Dendritic cells transport conidia and hyphae of Aspergillus fumigatus from the airways to the draining lymph nodes and initiate disparate Th responses to the fungus. J Immunol 168: 1362–1371.
  14. 14. den Dunnen J, Gringhuis SI, Geijtenbeek TB (2009) Innate signaling by the C-type lectin DC-SIGN dictates immune responses. Cancer Immunol Immunother 58: 1149–1157.
  15. 15. 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.
  16. 16. Saijo S, Iwakura Y (2011) Dectin-1 and Dectin-2 in innate immunity against fungi. Int Immunol.
  17. 17. Reid DM, Gow NA, Brown GD (2009) Pattern recognition: recent insights from Dectin-1. Curr Opin Immunol 21: 30–37.
  18. 18. Mezger M, Kneitz S, Wozniok I, Kurzai O, Einsele H, et al. (2008) Proinflammatory response of immature human dendritic cells is mediated by dectin-1 after exposure to Aspergillus fumigatus germ tubes. J Infect Dis 197: 924–931.
  19. 19. Cunha C, Di Ianni M, Bozza S, Giovannini G, Zagarella S, et al. (2010) Dectin-1 Y238X polymorphism associates with susceptibility to invasive aspergillosis in hematopoietic transplantation through impairment of both recipient- and donor-dependent mechanisms of antifungal immunity. Blood 116: 5394–5402.
  20. 20. Bochud PY, Chien JW, Marr KA, Leisenring WM, Upton A, et al. (2008) Toll-like receptor 4 polymorphisms and aspergillosis in stem-cell transplantation. N Engl J Med 359: 1766–1777.
  21. 21. Carvalho A, Pasqualotto AC, Pitzurra L, Romani L, Denning DW, et al. (2008) Polymorphisms in toll-like receptor genes and susceptibility to pulmonary aspergillosis. J Infect Dis 197: 618–621.
  22. 22. Kesh S, Mensah NY, Peterlongo P, Jaffe D, Hsu K, et al. (2005) TLR1 and TLR6 polymorphisms are associated with susceptibility to invasive aspergillosis after allogeneic stem cell transplantation. Ann N Y Acad Sci 1062: 95–103.
  23. 23. de Boer MG, Jolink H, Halkes CJ, van der Heiden PL, Kremer D, et al. (2011) Influence of polymorphisms in innate immunity genes on susceptibility to invasive aspergillosis after stem cell transplantation. PLoS One 6: e18403.
  24. 24. Chai LY, de Boer MG, van der Velden WJ, Plantinga TS, van Spriel AB, et al. (2011) The Y238X stop codon polymorphism in the human beta-glucan receptor dectin-1 and susceptibility to invasive aspergillosis. J Infect Dis 203: 736–743.
  25. 25. Sainz J, Perez E, Gomez-Lopera S, Jurado M (2008) IL1 gene cluster polymorphisms and its haplotypes may predict the risk to develop invasive pulmonary aspergillosis and modulate C-reactive protein level. J Clin Immunol 28: 473–485.
  26. 26. Sainz J, Hassan L, Perez E, Romero A, Moratalla A, et al. (2007) Interleukin-10 promoter polymorphism as risk factor to develop invasive pulmonary aspergillosis. Immunol Lett 109: 76–82.
  27. 27. Seo KW, Kim DH, Sohn SK, Lee NY, Chang HH, et al. (2005) Protective role of interleukin-10 promoter gene polymorphism in the pathogenesis of invasive pulmonary aspergillosis after allogeneic stem cell transplantation. Bone Marrow Transplant 36: 1089–1095.
  28. 28. Mezger M, Steffens M, Beyer M, Manger C, Eberle J, et al. (2008) Polymorphisms in the chemokine (C-X-C motif) ligand 10 are associated with invasive aspergillosis after allogeneic stem-cell transplantation and influence CXCL10 expression in monocyte-derived dendritic cells. Blood 111: 534–536.
  29. 29. Sainz J, Salas-Alvarado I, Lopez-Fernandez E, Olmedo C, Comino A, et al. (2010) TNFR1 mRNA expression level and TNFR1 gene polymorphisms are predictive markers for susceptibility to develop invasive pulmonary aspergillosis. Int J Immunopathol Pharmacol 23: 423–436.
  30. 30. Sainz J, Perez E, Hassan L, Moratalla A, Romero A, et al. (2007) Variable number of tandem repeats of TNF receptor type 2 promoter as genetic biomarker of susceptibility to develop invasive pulmonary aspergillosis. Hum Immunol 68: 41–50.
  31. 31. De Pauw B, Walsh TJ, Donnelly JP, Stevens DA, Edwards JE, et al. (2008) Revised definitions of invasive fungal disease from the European Organization for Research and Treatment of Cancer/Invasive Fungal Infections Cooperative Group and the National Institute of Allergy and Infectious Diseases Mycoses Study Group (EORTC/MSG) Consensus Group. Clin Infect Dis 46: 1813–1821.
  32. 32. Marr KA, Balajee SA, McLaughlin L, Tabouret M, Bentsen C, et al. (2004) Detection of galactomannan antigenemia by enzyme immunoassay for the diagnosis of invasive aspergillosis: variables that affect performance. J Infect Dis 190: 641–649.
  33. 33. Maertens J, Theunissen K, Verbeken E, Lagrou K, Verhaegen J, et al. (2004) Prospective clinical evaluation of lower cut-offs for galactomannan detection in adult neutropenic cancer patients and haematological stem cell transplant recipients. Br J Haematol 126: 852–860.
  34. 34. Maertens J, Van Eldere J, Verhaegen J, Verbeken E, Verschakelen J, et al. (2002) Use of circulating galactomannan screening for early diagnosis of invasive aspergillosis in allogeneic stem cell transplant recipients. J Infect Dis 186: 1297–1306.
  35. 35. 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.
  36. 36. 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.
  37. 37. 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.
  38. 38. 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.
  39. 39. Chen B, Wilkening S, Drechsel M, Hemminki K (2009) SNP_tools: A compact tool package for analysis and conversion of genotype data for MS-Excel. BMC Res Notes 2: 214.
  40. 40. Gabriel SB, Schaffner SF, Nguyen H, Moore JM, Roy J, et al. (2002) The structure of haplotype blocks in the human genome. Science 296: 2225–2229.
  41. 41. 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.
  42. 42. Ritchie MD, Hahn LW, Moore JH (2003) Power of multifactor dimensionality reduction for detecting gene-gene interactions in the presence of genotyping error, missing data, phenocopy, and genetic heterogeneity. Genet Epidemiol 24: 150–157.
  43. 43. Moore JH (2004) Computational analysis of gene-gene interactions using multifactor dimensionality reduction. Expert Rev Mol Diagn 4: 795–803.
  44. 44. Ok M, Einsele H, Loeffler J (2011) Genetic susceptibility to Aspergillus fumigatus infections. Int J Med Microbiol 301: 445–452.
  45. 45. Asakura Y, Komatsu T (2009) Toll-like receptor 4 polymorphisms and aspergillosis. N Engl J Med 360–635; author reply 635–636.
  46. 46. Zaas AK, Liao G, Chien JW, Weinberg C, Shore D, et al. (2008) Plasminogen alleles influence susceptibility to invasive aspergillosis. PLoS Genet 4: e1000101.
  47. 47. Martin MP, Lederman MM, Hutcheson HB, Goedert JJ, Nelson GW, et al. (2004) Association of DC-SIGN promoter polymorphism with increased risk for parenteral, but not mucosal, acquisition of human immunodeficiency virus type 1 infection. J Virol 78: 14053–14056.
  48. 48. Smith MW, Dean M, Carrington M, Winkler C, Huttley GA, et al. (1997) Contrasting genetic influence of CCR2 and CCR5 variants on HIV-1 infection and disease progression. Hemophilia Growth and Development Study (HGDS), Multicenter AIDS Cohort Study (MACS), Multicenter Hemophilia Cohort Study (MHCS), San Francisco City Cohort (SFCC), ALIVE Study. Science 277: 959–965.
  49. 49. Koizumi Y, Kageyama S, Fujiyama Y, Miyashita M, Lwembe R, et al. (2007) RANTES -28G delays and DC-SIGN - 139C enhances AIDS progression in HIV type 1-infected Japanese hemophiliacs. AIDS Res Hum Retroviruses 23: 713–719.
  50. 50. Kashima S, Rodrigues ES, Azevedo R, da Cruz Castelli E, Mendes-Junior CT, et al. (2009) DC-SIGN (CD209) gene promoter polymorphisms in a Brazilian population and their association with human T-cell lymphotropic virus type 1 infection. J Gen Virol 90: 927–934.
  51. 51. Mezger M, Steffens M, Semmler C, Arlt EM, Zimmer M, et al. (2008) Investigation of promoter variations in dendritic cell-specific ICAM3-grabbing non-integrin (DC-SIGN) (CD209) and their relevance for human cytomegalovirus reactivation and disease after allogeneic stem-cell transplantation. Clin Microbiol Infect 14: 228–234.
  52. 52. Vannberg FO, Chapman SJ, Khor CC, Tosh K, Floyd S, et al. (2008) CD209 genetic polymorphism and tuberculosis disease. PLoS One 3: e1388.
  53. 53. Ganachari M, Ruiz-Morales JA, Gomez de la Torre Pretell JC, Dinh J, Granados J, et al. (2010) Joint effect of MCP-1 genotype GG and MMP-1 genotype 2G/2G increases the likelihood of developing pulmonary tuberculosis in BCG-vaccinated individuals. PLoS One 5: e8881.
  54. 54. Flores-Villanueva PO, Ruiz-Morales JA, Song CH, Flores LM, Jo EK, et al. (2005) A functional promoter polymorphism in monocyte chemoattractant protein-1 is associated with increased susceptibility to pulmonary tuberculosis. J Exp Med 202: 1649–1658.
  55. 55. Barreiro LB, Neyrolles O, Babb CL, Tailleux L, Quach H, et al. (2006) Promoter variation in the DC-SIGN-encoding gene CD209 is associated with tuberculosis. PLoS Med 3: e20.
  56. 56. Ryan EJ, Dring M, Ryan CM, McNulty C, Stevenson NJ, et al. (2010) Variant in CD209 promoter is associated with severity of liver disease in chronic hepatitis C virus infection. Hum Immunol 71: 829–832.
  57. 57. Muhlbauer M, Bosserhoff AK, Hartmann A, Thasler WE, Weiss TS, et al. (2003) A novel MCP-1 gene polymorphism is associated with hepatic MCP-1 expression and severity of HCV-related liver disease. Gastroenterology 125: 1085–1093.
  58. 58. Park BL, Kim YJ, Cheong HS, Kim LH, Choi YH, et al. (2006) Association of common promoter polymorphisms of MCP1 with hepatitis B virus clearance. Exp Mol Med 38: 694–702.
  59. 59. Chan KY, Xu MS, Ching JC, So TM, Lai ST, et al. (2010) CD209 (DC-SIGN) -336A>G promoter polymorphism and severe acute respiratory syndrome in Hong Kong Chinese. Hum Immunol 71: 702–707.
  60. 60. Sainz J, Segura-Catena J, Jurado M (2010) [Association between genetic polymorphism in the promotor region of CD209 and propensity to develop invasive pulmonary aspergillosis]. Methods Find Exp Clin Pharmacol 32: Suppl A9–13.
  61. 61. Werner JL, Gessner MA, Lilly LM, Nelson MP, Metz AE, et al. (2011) Neutrophils produce IL-17A in a Dectin-1 and IL-23 dependent manner during invasive fungal infection. Infect Immun.
  62. 62. Dennehy KM, Brown GD (2007) The role of the beta-glucan receptor Dectin-1 in control of fungal infection. J Leukoc Biol 82: 253–258.
  63. 63. Plantinga TS, van der Velden WJ, Ferwerda B, van Spriel AB, Adema G, et al. (2009) Early stop polymorphism in human DECTIN-1 is associated with increased candida colonization in hematopoietic stem cell transplant recipients. Clin Infect Dis 49: 724–732.